Polyhedral billiards, eigenfunction concentration and almost periodic control

TL;DR

This paper investigates the dynamics of billiard flows in convex polyhedra, establishes finiteness and length estimates for periodic tubes, and applies these results to eigenfunction concentration, introducing new control estimates for almost-periodic boundary conditions.

Contribution

It extends known results on billiard dynamics to higher dimensions, provides quantitative estimates for periodic tubes, and introduces a novel control estimate for almost-periodic boundary conditions.

Findings

Finitely many immersed periodic tubes exist outside pockets.

Quantitative estimates for lengths of periodic tubes.

Eigenfunction mass concentrates near pockets in convex polyhedral billiards.

Abstract

We study dynamical properties of the billiard flow on convex polyhedra away from a neighbourhood of the non-smooth part of the boundary, called ``pockets''. We prove there are only finitely many immersed periodic tubes missing the pockets and moreover establish a new quantitative estimate for the lengths of such tubes. This extends well-known results in dimension $2$. We then apply these dynamical results to prove a quantitative Laplace eigenfunction mass concentration near the pockets of convex polyhedral billiards. As a technical tool for proving our concentration results on irrational polyhedra, we establish a control-theoretic estimate on a product space with an almost-periodic boundary condition. This extends previously known control estimates for periodic boundary conditions, and seems to be of independent interest.

Equations

$${\int }_U|u{|}^2\ge c{\int }_P|u{|}^2.$$

$${\int }_U|u{|}^2\ge c.$$

$$U$$

$${\int }_U|u{|}^2\ge c.$$

$$\underset{t\in \mathbb{R}}{\sup }\Vert f(t+\tau )-f(t)\Vert \le \epsilon .$$

$$\mathcal{M}\{f\}:=\underset{T\to \infty }{\lim }\frac{1}{2T}{\int }_{-T}^{T}f(t)dt.$$

$$a(\lambda ;f):=\underset{T\to \infty }{\lim }\frac{1}{2T}{\int }_{-T}^{T}f(t)e^{-i\lambda t}dt=\mathcal{M}\{f(t)e^{-i\lambda t}\}.$$

$$f(t)\sim \sum _{k=1}^{\infty }{a}_k{e}^{i{\lambda }_{k}t}.$$

$$\mathcal{M}\{\Vert f(t){\Vert }^2\}=\sum _{k=1}^{\infty }\Vert a_k{\Vert }^2<\infty .$$

$$\mathcal{M}\{\Vert f(t)-\sum _{k=1}^n{a}_k{e}^{i{\lambda }_{k}t}{\Vert }^2\}=\sum _{k=n+1}^{\infty }\Vert a_k{\Vert }^2\to 0.$$

$$u(x,t+L)=u(\phi (x),t),(x,t)\in M_x\times \mathbb{R}.$$

$$u(x,L)=u(\phi (x),0),x\in M_x.$$

$$S(\phi ,\epsilon )=\{k\in \mathbb{Z}\mid \mathrm{dist}({\phi }^k,\mathrm{id})<\epsilon \}.$$

$$g:\mathbb{R}\ni t\mapsto u(\cdot ,t)\in H^s(M_x)$$

$$\mathrm{dist}(\psi ,\mathrm{id})<\epsilon ⟹\Vert u(\psi (x),t)-u(x,t){\Vert }_{H^s(M_x)}<\delta$$

$$\underset{t\in \mathbb{R}}{\sup }\Vert g(t+Lk)-g(t){\Vert }_{H^s(M_x)}=\underset{t\in \mathbb{R}}{\sup }\Vert u({\phi }^k(x),t)-u(x,t){\Vert }_{H^s(M_x)}<\delta$$

$$\phi =R_{{\alpha }_1,\dots ,{\alpha }_m}:(x_1,\dots ,x_m)\mapsto (x_1+{\alpha }_1,\dots ,x_m+{\alpha }_m)\mathrm{mod}{\mathbb{Z}}^m$$

$$\mathcal{T}=\{k({\alpha }_1,\dots ,{\alpha }_m)\mathrm{mod}{\mathbb{Z}}^m:k\in \mathbb{Z}\}.$$

$$t=k({\alpha }_1,\dots ,{\alpha }_m)\mathrm{mod}{\mathbb{Z}}^m.$$

$$l:=2\max \{|k_1|,\dots ,|k_N|\}+1$$

$$S_{C^r}(\phi ,\epsilon )=\{k\in \mathbb{Z}\mid {\mathrm{dist}}_{C^r}({\phi }^k,\mathrm{id})<\epsilon \}$$

$${\mathcal{O}}_x:=\overline{\{f^{m}x\mid m\in \mathbb{Z}\}}.$$

$$\epsilon >d(f^{m}x,f^{n_i}x)=d(f^{m-n_i}x,x).$$

$$(f,f,\dots ,f):(M_x{)}^N\to (M_x{)}^N,$$

$$S(m):=\{a\in C^{\infty }({\mathbb{R}}^{2n}):|{\partial }^{\beta }a(x,\xi )|\le C_{\beta }⟨\xi {⟩}^m,\mathrm{for}\mathrm{each}\mathrm{multiindex}\beta \}.$$

$${\mathcal{F}}_{h}u(\xi )={\int }_{{\mathbb{R}}^n}{e}^{-i\frac{y\cdot \xi }{h}}u(y)dy$$

$$a\big{(}x_{1},\dotso,x_{n-1},x_{n}+L,\xi_{1},\dotso,\xi_{n-1},\xi_{n}\big{)}=a\big{(}\mathcal{R}(x_{1},\dotso,x_{n-1}),x_{n},\mathcal{R}^{-1}(\xi_{1},\dotso,\xi_{n-1}),\xi_{n}\big{)}.$$

$${\mathcal{F}}_h(u\circ \phi )(\xi )$$

$$=e^{iL\frac{{\xi }_n}{h}}{\int }_{{\mathbb{R}}^n}{e}^{-i\frac{y^{\prime }}{h}\cdot {\mathcal{R}}^{-1}(\xi )}u(y^{\prime })d{y}^{\prime }=e^{iL\frac{{\xi }_n}{h}}{\mathcal{F}}_h(u)({\mathcal{R}}^{-1}(\xi )),$$

$$\begin{gathered} (2\pi h{)}^n\Phi ({\phi }^{\ast }u)(x)={\int }_{{\mathbb{R}}^n}{e}^{i\frac{x\cdot \xi }{h}}{\mathcal{F}}_h(u\circ \phi )(\xi )a(\xi )d\xi \\ ={\int }_{{\mathbb{R}}^n}{e}^{iL\frac{{\xi }_n^{\prime }}{h}}{e}^{i\frac{x}{h}\cdot \mathcal{R}({\xi }^{\prime })}{\mathcal{F}}_h(u)({\xi }^{\prime })a({\xi }^{\prime })d{\xi }^{\prime }={\int }_{{\mathbb{R}}^n}{e}^{i\phi (x)\cdot \frac{{\xi }^{\prime }}{h}}{\mathcal{F}}_h(u)({\xi }^{\prime })a({\xi }^{\prime })d{\xi }^{\prime }, \end{gathered}$$

$$a(hD)u(x)=\sum _{k\in \mathbb{Z}}{A}_{k}u(x),x\in {\mathbb{R}}^{n-1}\times [0,L),$$

Full text

1. Introduction

In the present paper our interests are twofold. The first line of studies is related to several dynamical properties of the billiard flow on convex polyhedra in $\mathbb{R}^{n}$. The second line of questions addresses an allied spectral problem motivated from quantum physics: namely, the problem of high energy eigenfunction concentration.

Concerning the dynamical perspective, there has been a lot of interest in the billiard flow on polygons – for example, we refer to [14, 10, 6, 9, 7, 19] etc. For a more comprehensive overview, one might also wish to consult the extensive work [20]. As is well-known, this area of research contains questions that might seem quite simple looking at first, but turn out to be rather deep on closer investigation: an illustration is given by Galperin’s example of non-periodic and not everywhere dense billiard trajectories in convex two-dimensional polygons (see [13]).

Moreover, when one increases the number of dimensions and studies the billiard flow on polyhedra (either convex or non-convex) further challenges arise. To get an insight into some of the issues involved in generalising reasonably well-known two-dimensional results to three dimensions, we refer the reader to [14]; see also [3], which gives a three-dimensional analogue of a result by Katok (see [19]) regarding the entropy of the symbolic dynamics of the billiard flow on polygons to higher dimensions.

In this work we consider the billiard flow on polyhedra, particularly the interaction of the flow with the pockets of the billiard (i.e. an open neighbourhood around the singular set of the polytope - cf. Section 3.1). Our main interest is in the finer properties of billiard dynamics away from the pockets. In this direction our results say that for convex polyhedra, the “non-singular part” of the phase space away from the pockets essentially decomposes into a collection of finitely many immersed “periodic tubes”. Moreover, we show that the lengths of the immersed “periodic tubes” satisfy a certain quantitative estimate.

We now turn to the second line of research devoted to the study of the underlying spectral problem. Let $-\Delta$ denote the Dirichlet (or Neumann, as the case might be) Laplacian on a convex polyhedron $P\subset\mathbb{R}^{n}$. In [16] concentration of eigenfunctions on certain flat polygonal surfaces has been investigated. The result states that given a polygon $P$ whose vertex set is denoted by $V$ and an open neighbourhood $U$ of $V$, there exists a positive constant $c(U)$ such that for any Dirichlet eigenfunction $u$ (satisfying $-\Delta u=\lambda u$) the following concentration estimate holds:

[TABLE]

In other words, a certain amount of eigenfunction mass collects near the “singular points” of the boundary $\partial P$, or, $U$ analytically controls $P$ in a certain sense (see [26]). Clearly, this can be thought of as an ergodicity phenomenon. We study this phenomenon in further settings, namely higher dimensional polytopes. Since the Laplace eigenfunctions are eigenfunctions of the wave propagator $U(t):=e^{it\sqrt{-\Delta}}$, it is expected (at least heuristically) that the high-frequency limits $\lambda\to\infty$ should reflect the dynamics of the classical geodesic flow. In particular, when the geodesic flow is ergodic, one should expect the eigenfunctions to diffuse in phase space, which gives an intuitive feeling as to why estimates like (1.1) should be expected to hold (see also [28]). However, note that a typical neighbourhood of the singular points of $\partial P$ is not large enough to geometrically control $P$ (in the sense of [1, 7, 8]), and special properties of the billiard flow on polyhedra need to be used. For other studies of eigenfunctions on polygons see [24, 23]; see also [5, 27] for general overviews of a somewhat different flavour.

For our main result, consider a convex polyhedron $P\subset\mathbb{R}^{n}$ and write $\mathcal{S}$ for the singular set of the boundary $\partial P$. By definition $\mathcal{S}\subset\partial P$ is the union of faces of dimension $\leq n-2$, or in other words, the $(n-2)$-skeleton of $\partial P$. This brings us to our main result

Theorem 1.1**.**

Given a neighbourhood $U$ of $\mathcal{S}$ inside the polyhedron $P$, there exists a positive constant $c=c(U)$ such that for any $L^{2}$-normalised Dirichlet (or Neumann) eigenfunction $u$, we have that

[TABLE]

Again, we observe that the claim of Theorem 1.1 holds for all eigenfunctions, as opposed to “almost all” eigenfunctions (i.e. a density $1$ subsequence) as usually established in quantum ergodicity statements.

Before getting into more details, we outline the strategy presented in [16] that resolves the two-dimensional case $n=2$. The idea is based upon the following two special properties of the billiard flow on polygonal domains:

[TABLE]

[TABLE]

Now, let us assume to the contrary that a sequence of eigenfunctions concentrates away from $U$. The two dynamical properties mentioned above, together with propagation results for eigenfunctions, imply that the corresponding semiclassical measure on the sphere bundle must concentrate along the families of periodic geodesics that sweep out such cylinders. This is ruled out using an argument from [24], which in turn is based on an estimate by Burq and Zworski (see [7, 8]). This concludes the discussion and proves the concentration estimate over $U$.

1.1. A brief outline of our results

We first discuss the dynamical side of our results. Our main innovations in this direction will be to discover suitable replacements and generalisations of Properties (P1) and (P2) to higher dimensional polyhedra.

It is known (see [14]) that a billiard trajectory which avoids a neighbourhood of the singular points need not be periodic itself, but it is contained in an immersed tube that “closes up”, with a suitable cross-section. Such a tube has the property that (a certain iterate of) the Poincaré return map on the cross-section of the tube is a rotation. This can be used as a substitute for Property (P1) above, and we will call such tubes “periodic tubes” with slight abuse of language (see Definition 3.4 below and the explanations that follow).

As we will demonstrate below, Property (P2) will also generalise in our setting, in the sense that the number of such immersed periodic tubes missing a neighbourhood of $\mathcal{S}$ is also finite, which we call the “finite tube condition”. Moreover, we further strengthen this claim, as we obtain a quantitative estimate for the lengths of such tubes. A novelty here is that one needs to quantify the formation of singular points on the boundary of these tubes, which involves a certain dynamical invariant of the associated rotation.

However, we will show that a much weaker property than the finite tube condition is sufficient for the application to eigenfunctions, and we call this Property (P2’) (see Section 3.4 below). It roughly says that for each periodic tube, the flow-out of a neighbourhood of the singular set contains a conical vicinity of the tube direction, near the boundary of the tube. We emphasise that despite the fact that the finite tube condition could be used in our proofs, we regard the Property (P2’) as a less strict and more natural condition, which may lead to further applications of our analytic argument.

From the analytical standpoint, we prove that the Burq-Zworski type control result used by [16] extends easily and applies to the case when $P$ is rational. However, in the case of an irrational polyhedron, further complications arise as our periodic tubes consist of billiard trajectories that are no longer individually periodic, but are “almost periodic”, except for one “central” trajectory. This motivates us to introduce the new notion of an almost periodic boundary condition on a cylinder. Our main contribution in this analytical part is proving a version of a control result due to Burq and Zworski (see [8, Proposition 6.1]) that holds in higher dimensions for almost periodic boundary conditions. This is enough to address the case of irrational polyhedra and could potentially have other applications. We mention in passing that our methods also yield a generalisation of [24, Theorem 2], see Theorem 6.2 below. This basically says that given any immersed (almost) periodic tube $T$ in $P$, no eigenfunction can concentrate in $T$ and away from $\partial T$.

Finally, we mention that it is possible to state a more general version of Theorem 1.1 in a somewhat more abstract setting, as follows. Take a convex domain $\Omega\subset\mathbb{R}^{n}$ with piecewise-smooth boundary, and which is such that $\partial\Omega$ consists of a finite disjoint union of convex $(n-1)$-dimensional polyhedra $\cup_{i=1}^{k}\mathcal{F_{i}}$, which are called the “faces” of $\Omega$, and the set $\mathcal{S}:=\partial\Omega\setminus\cup_{i=1}^{k}\mathcal{\operatorname{int}F_{i}}$. A variant of the standard billiard flow can be introduced on $\Omega$ which is the usual billiard flow with the stipulation that the billiard particle is stopped when it hits $\mathcal{S}$. Then it can be checked that the proof of Theorem 1.1 goes through verbatim to give us the following corollary:

Corollary 1.2**.**

Let $U$ be any neighbourhood of $\mathcal{S}$ inside $\Omega$. Then, there exists a positive constant $c=c(U)$ such that for any $L^{2}$-normalised Dirichlet (or Neumann) eigenfunction $u$, we have that

[TABLE]

As an application of this corollary, one can say that if one has a cube with smoothened edges (eg., a die), then a certain amount of $L^{2}$-mass collects at any neighbourhood of the smoothened edges. As another application, consider the $3$-dimensional equivalent of the Bunimovich stadium, that is, a rectangular parallelepiped, with two topological hemispheres fitted smoothly at the ends. The above corollary will dictate that the $L^{2}$-mass of eigenfunctions cannot totally concentrate away from the “wings” (hemispheres) of the stadium.

The main upshot is that polyhedra are in no way privileged objects, and the structure of the “singularity” of the $(n-2)$-skeleton does not really play a special role. However, for aesthetic reasons, we prefer to state our main result in a more intuitive and geometrically appealing setting, as in Theorem 1.1.

1.2. Open questions and further work

There are quite a few interesting questions left to consider. It would be interesting to speculate if any of the dynamical results stated here have analogues for polyhedra which are not necessarily convex. As a starting point, one can check whether statements like Theorem 3.7 work for not so badly non-convex objects, like a non-convex polyhedron which is formed from two convex polyhedra attached at a common face. From the analytical point of view, an exciting question is to determine the dependence of the constant $c=c(\varepsilon)$ in Theorem 1.1 in the case when $U$ is an $\varepsilon$-neighbourhood of the singular set, i.e. to quantify the estimates in [16] and this paper, at least asymptotically as $\varepsilon\searrow 0$.

1.3. Structure of the paper

In Section 2 we introduce and study the new notion of functions satisfying an almost periodic boundary condition on a cylinder. We also recall the theory of almost periodic functions and study pseudo-differential operators on a mapping torus.

As mentioned before, our paper splits naturally into a dynamical and an analytical part. Section 3 addresses some of the main dynamical components of our paper: we prove the finite tube condition (see Theorem 3.12) and the Property (P2’) (see Lemma 3.14). In Section 4, we prove the other main dynamical result: a qualitative bound on lengths of periodic tubes missing a neighbourhood of the singular set (see Theorem 4.5). Then, in Section 5, we discuss the analytical results required to prove Theorem 1.1, concerning control-theoretic estimates with an almost periodic boundary condition (see Theorem 5.4). In Section 6 we prove Theorem 1.1 first in the (easier) case of rational polyhedra and then go on to prove it in full generality.

In the appendices, for the convenience of the reader and to make this work self-contained, we give important details about billiard dynamics and control theory. In Appendix A we discuss theorems about billiard dynamics proved in [14] and present proofs that work in any dimension. In Appendix B we present a proof of a result of N. Burq in control theory on a product space with a periodic boundary condition.

2. Preliminaries

In the first part of this section, we discuss and recall the theory of almost periodic functions with values in a Banach space. We then introduce the notion of an admissible isometry and explain how admissible isometries give rise to almost periodic boundary conditions. In fact, we show that every isometry is admissible. In the last part, we discuss the theory of pseudo-differential operators on a mapping torus.

2.1. Almost periodic functions

We introduce the theory of almost periodic functions, as developed by H. Bohr in 1920s and later generalised by others. We will follow mostly the first two chapters of [22]. For this purpose, let $X$ be a Banach space with norm $\lVert{\cdot}\rVert$. We will say a number $\tau\in\mathbb{R}$ is an $\varepsilon$-almost period of $f:\mathbb{R}\to X$ if

[TABLE]

We also say that a subset $E\subset\mathbb{R}$ is relatively dense if there is an $l>0$ such that for any $\alpha\in\mathbb{R}$, the interval $(\alpha,\alpha+l)\subset\mathbb{R}$ of length $l$ contains an element of $E$.

We start with a basic definition:

Definition 2.1**.**

A continuous function $f:\mathbb{R}\to X$ is called almost periodic if for every $\varepsilon>0$, there is an $l=l(\varepsilon)>0$ such that for each $\alpha\in\mathbb{R}$, the interval $(\alpha,\alpha+l)\subset\mathbb{R}$ contains a number $\tau=\tau(\varepsilon)$ such that (2.1) holds.

An immediate observation is that if $f$ is periodic, then it is almost periodic. A simple example of a non-periodic, but almost periodic function is given by $f(t)=\sin t+\sin(\sqrt{2}t)$. Also, any almost periodic function is uniformly continuous [22, Chapter 1]. Another equivalent definition is due to S. Bochner and says that (cf. [22, p. 4])

Definition 2.2**.**

Let $f:\mathbb{R}\to X$ be continuous. For $h\in\mathbb{R}$, we define $f^{h}(t):=f(t+h)$. Then $f$ is almost periodic if and only if the family of functions $\{f^{h}\mid h\in\mathbb{R}\}$ is compact in the topology of uniform convergence on $\mathbb{R}$.

Next, we discuss expansion into trigonometric polynomials, i.e. harmonic analysis, similarly to the case of periodic functions. The fundamental theorem in this area is the Approximation Theorem [22, p. 17] which says that every almost periodic $f$ is a uniform limit of sums of trigonometric polynomials. In other words, for every $\varepsilon>0$, there is a sum $\sum_{k=1}^{n_{\varepsilon}}e^{i\lambda_{k}t}a_{k}(\varepsilon)$ that is $\varepsilon$-close to $f(t)$ in the uniform norm, where $a_{k}(\varepsilon)\in X$ and $\{\lambda_{k}\}_{k=1}^{\infty}\subset\mathbb{R}$ is a countable set of exponents called the spectrum of $f$. Clearly for periodic functions on $[0,1]$, we may take $\lambda_{0}=0$ and $\lambda_{2k-1}=2k\pi,\lambda_{2k}=-2k\pi$ for positive integer values of $k$. Several main properties of almost periodic functions can be deduced from the Approximation Theorem.

We will denote the mean value of an almost periodic function $f$ (it may be shown that it exists, see [22, p. 22])

[TABLE]

We now define the Bohr transformation $a(\lambda;f)$ of $f$ for $\lambda\in\mathbb{R}$ as a “mean value Fourier transform”

[TABLE]

The values $\lambda=\lambda_{k}$ in the spectrum of $f$ are then exactly the values for which $a(\lambda;f)\neq 0$. We will write formally, with no convergence implied, for $a_{k}:=a(\lambda_{k};f)$

[TABLE]

However, one can show that certain Bochner-Fejér sums converge uniformly to $f$, obtained by taking partial sums of the right hand side of (2.3) and applying suitable multipliers (see [22, Section 2.4]). One may also show, using the Approximation Theorem, that the Fourier coefficients $a_{k}(\lambda_{k};f)$ are uniquely associated to $f$. In other words, if $f$ and $g$ are two almost periodic functions, then $a(\lambda;f)=a(\lambda;g)$ for every $\lambda\in\mathbb{R}$ implies $f\equiv g$ (see [22, p. 24]).

Next, we assume that $X$ is a Hilbert space and state a Parseval-type identity, which says that if (2.3) holds, then

[TABLE]

For a proof, see [22, p. 31]. From (2.4) it also follows that, as $n\to\infty$

[TABLE]

2.2. Admissible isometries and almost periodic boundary conditions

We consider a compact Riemannian manifold $(M_{x},g_{x})$ with Lipschitz boundary $\partial M_{x}$, where we have introduced the lower index notation to indicate explicitly that points on the manifold will be denoted by $x$. In the rest of the paper, we will deal with functions $u:M_{x}\times\mathbb{R}\to\mathbb{C}$, satisfying some invariance properties

[TABLE]

Here $L>0$ is a positive number (length) and $\varphi:M_{x}\to M_{x}$ is an isometry. It is clear that the invariance (2.5) can be interpreted as a boundary condition for $u$ on $M_{x}\times[0,L]$

[TABLE]

In order to address the desired control estimates (cf. Theorem 5.4 below), we first need to impose a condition on the isometry $\varphi$ that generalises the periodic case ($\varphi=\operatorname{id}$).

We call an isometry $\varphi:M_{x}\to M_{x}$ with a corresponding induced isometry of the boundary $\varphi|_{\partial M_{x}}:\partial M_{x}\to\partial M_{x}$ admissible, if for every $\varepsilon>0$, the set

[TABLE]

is relatively dense in $\mathbb{Z}$. In this definition, we include the possibility that the boundary $\partial M_{x}$ is empty. A set $A\subset\mathbb{Z}$ ($A\subset\mathbb{N}$) is relatively dense if there exists an $N\in\mathbb{N}$ such that every consecutive $N$ integers (positive integers) contain an element of $A$, i.e. every set of the form $\{k,k+1,\dotso,k+N-1\}$ for $k\in\mathbb{Z}$ ($k\in\mathbb{N}$) contains an element of $A$. Here $\operatorname{dist}(\cdot,\cdot)$ denotes distance between mappings in $C^{\infty}(M_{x},M_{x})$. If the isometry $\varphi$ is admissible, we call the boundary condition in (2.6) almost periodic.

Denote by $H^{s}(M_{x})$ the Sobolev space of index $s\in\mathbb{R}$. The relation between almost periodic functions and admissible isometries is given by

Lemma 2.3**.**

Let $\varphi:M_{x}\to M_{x}$ be an admissible isometry and $u:M_{x}\times\mathbb{R}\to\mathbb{C}$ be such that $u\in C(\mathbb{R},H^{s}(M_{x}))$ for some $s\in\mathbb{R}$, satisfying that $u(x,t+L)=u(\varphi(x),t)$ for all $(x,t)$ and some $L>0$ fixed. Then the map

[TABLE]

is almost periodic.

Proof.

Let $\varepsilon>0$. By assumption, $S(\varphi,\varepsilon)\subset\mathbb{N}$ is relatively dense, so we may pick $N_{\varepsilon}\in\mathbb{N}$ such that $S(\varphi,\varepsilon)$ has non-empty intersection with any $N_{\varepsilon}$ consecutive positive integers.

Fix now $\delta>0$ and choose $\varepsilon>0$ small enough such that, for $\psi:M_{x}\to M_{x}$ smooth

[TABLE]

for all $t\in[0,L]$. By the invariance property of $u$ and since $\varphi^{*}$ is an isometric isomorphism on $H^{s}(M_{x})$, this also holds for $t\in\mathbb{R}$.

Define $\mathcal{I}(\delta):=S(\varphi,\varepsilon)L\subset\mathbb{R}$. By the defining properties of $S(\varphi,\varepsilon)$, any interval of length $LN_{\varepsilon}$ contains an element of $\mathcal{I}(\delta)$, i.e. $\mathcal{I}(\delta)$ is relatively dense in $\mathbb{R}$. Also, for any $k\in S(\varphi,\varepsilon)$

[TABLE]

by (2.9). Thus the set of $\delta$-almost periods is relatively dense for any $\delta>0$ and so $g$ is almost periodic. ∎

In a moment, we are going to demonstrate that all isometries are in fact admissible. Before proving the general case, we first discuss why any arbitrary rotation $R\in SO(n)$ acting on the closed unit ball $D^{n}\subset\mathbb{R}^{n}$ is admissible.

We first show that for any $m\in\mathbb{N}$ the higher dimensional rotations on the $m$-torus \mathbb{T}^{m}:=\big{(}\frac{\mathbb{R}}{\mathbb{Z}}\big{)}^{m}, generated by an $m$-tuple $(\alpha_{1},\dotso,\alpha_{m})\in\mathbb{R}^{m}$ and defined by

[TABLE]

are admissible. By ergodic theory of $\mathbb{T}^{m}$ (using Fourier expansions), we know that $\varphi$ is ergodic iff $1,\alpha_{1},\dotso,\alpha_{m}$ are linearly independent (l.i.) over $\mathbb{Q}$, i.e. iff $\alpha_{1},\dotso,\alpha_{m}$ are l.i. over $\mathbb{Z}$ taken modulo $\mathbb{Z}$. We denote $\alpha:=(\alpha_{1},\dotso,\alpha_{m})\mod\mathbb{Z}^{m}$ and consider the orbit of $\alpha$

[TABLE]

By Kronecker’s theorem, we know that the closure $\overline{\mathcal{T}}\subset\mathbb{T}^{m}$ is a torus.111By using the ergodicity of rotations and Kronecker’s theorem one should be able to deduce the claim about $\varphi$ directly.

We give an elementary proof that $\varphi$ is admissible, not relying on Kronecker’s theorem or ergodicity of rotations. First observe that $\overline{\mathcal{T}}$ is a group under addition. Fix an $\varepsilon>0$. By compactness, we may take elements $t_{1}=k_{1}\alpha\mod\mathbb{Z}^{m},\dotso,t_{N}=k_{N}\alpha\mod\mathbb{Z}^{m}\in\mathcal{T}$, where $k_{i}\in\mathbb{Z}$, such that they are $\varepsilon$-dense, i.e. for any $t\in\overline{\mathcal{T}}$, there is an $i$ such that $t_{i}$ is $\varepsilon$-close to $t$. Consider now an arbitrary $k\in\mathbb{Z}$ and an element

[TABLE]

Then $-t\mod\mathbb{Z}^{m}\in\mathcal{T}$ and there is an index $i$ such that $t_{i}$ is $\varepsilon$-close to $-t$. Equivalently, we may say that $(k_{i}+k)\alpha\mod\mathbb{Z}^{m}$ is $\varepsilon$-close to zero. Therefore, if we put

[TABLE]

we see that in any $l$ consecutive integers we may find one, say $r$, such that $r\alpha\mod\mathbb{Z}^{m}$ is $\varepsilon$-close to zero. Thus $\varphi$ is admissible.

Now, for a rotation $R\in SO(n)$, one may take a unitary matrix $P$ so that $P^{-1}RP=Q$ is diagonal and has eigenvalues coming in pairs $(e^{i\alpha},e^{-i\alpha})$ for some generalised angles $\alpha\in\mathbb{R}$. By the discussion above and since the action of $Q$ is conjugate to a rotation on a torus as in (2.10), we get that for any $\varepsilon>0$, the set of indices $k$ such that $Q^{k}$ is $\varepsilon$-close to $\operatorname{id}$ is relatively dense. This implies that $R:D^{n}\to D^{n}$ is admissible.

We now generalise the main ideas in the preceding paragraphs to give a proof which will not require a precise representation for the orbits as done above. For a more abstract point of view of topological dynamics and also for the converse claim, see [15, Chapter 4].

Proposition 2.4**.**

Let $\varphi:M_{x}\to M_{x}$ be an isometry. Then $\varphi$ is admissible.

Proof.

We will inductively show that $\varphi$ is $C^{r}$-admissible for any $r\geq 0$, by which we mean that for each $\varepsilon>0$ the set

[TABLE]

is relatively dense, where the metric $\operatorname{dist}_{C^{r}}(\cdot,\cdot)$ is now the $C^{r}$-metric (cf. (2.7)). We start with the case $r=0$.

The claim hinges on proving that in compact, isometric systems, points are uniformly recurrent. For $x\in M$, consider the closure of both the forward and backward orbits

[TABLE]

By compactness, there is a set of indices $n_{1},\dotso,n_{k}$ such that the points $f^{n_{1}}x,\dotso,f^{n_{k}}x$ are $\varepsilon$-dense in $\mathcal{O}_{x}$. Observe first that the set of indices $m\in\mathbb{Z}$ for which $f^{m}x$ is $\varepsilon$-close to $x$ is relatively dense. For, by the hypothesis there is an $i$ with

[TABLE]

Therefore, $\tilde{N}=\max(n_{1},\dotso,n_{k})+1$ would work in the definition of relatively dense. Consider now the $\varepsilon$-ball $U=B(x,\varepsilon)$ around $x$ and its iterates and take any $m\in\mathbb{Z}$. Then by (2.12), there is an $i$ such that $f^{m-n_{i}}x\in f^{m-n_{i}}U\cap U\neq\emptyset$.

We come back to the main proof and let $U_{i}=B(x_{i},\varepsilon)$, $i=1,\dotso,N$ be a cover of $M$ by $\varepsilon$-balls. Apply now the previous paragraph to the isometry

[TABLE]

where on the left hand side $f$ appears $N$ times, and the point $(x_{1},x_{2},\dotso,x_{N})$. Thus, there is a relatively dense set of integers $\mathcal{T}\subset\mathbb{Z}$ with the property that $U_{i}\cap f^{k}(U_{i})\neq\emptyset$ for every $k\in\mathcal{T}$ and $i=1,\dotso,N$. We claim that $\mathcal{T}\subset S_{C^{0}}(\varphi,4\varepsilon)$. Let $k\in\mathcal{T}$, choose $y_{i}\in U_{i}\cap f^{k}(U_{i})$ and take any $x\in M$. By construction there is an $i$ with $x\in U_{i}$, so $d(x,y_{i})<2\varepsilon$. Since $f$ is an isometry and $f^{k}x\in f^{k}U_{i}$, we have $d(f^{k}x,y_{i})<2\varepsilon$. By triangle inequality $d(x,f^{k}x)<4\varepsilon$, which implies $\operatorname{dist}_{C^{0}}(\varphi^{k},\operatorname{id})\leq 4\varepsilon$ and proves the claim.

Next, we prove the inductive hypothesis. Observe that $f$ induces a map $df:SM_{x}\to SM_{x}$ between unit sphere bundles. Moreover, this map is an isometry for the restriction of the Sasaki metric to $SM_{x}$ (see e.g. [25, Chapter 1]). Thus we may apply the same argument as in the previous paragraph to prove the claim for $r=1$. Similarly, by taking unit sphere bundles inductively, we may prove the claim for any $r$. Therefore, $\varphi$ is admissible. ∎

2.3. PDOs on a mapping torus

We will denote by $\mathcal{R}\in SO(n)$ a rotation in $\mathbb{R}^{n}$ that leaves the $x_{n}$-axis fixed. The rotation $\mathcal{R}$ is identified with the rotation it induces in $(x_{1},\dotso,x_{n-1})$ coordinates on $\mathbb{R}^{n-1}\subset\mathbb{R}^{n}$. Also, define $\varphi(x_{1},x_{2},\dotso,x_{n}):=(\mathcal{R}^{-1}(x_{1},\dotso,x_{n-1}),x_{n}+L)$ for some positive $L$.

Consider a bounded open set $\Omega\subset\mathbb{R}^{n-1}$ invariant under the rotation $\mathcal{R}$ and define the mapping torus $\mathcal{C}_{\varphi}:=(\Omega\times[0,L])/\left((x,L)\sim(\mathcal{R}(x),0)\right)$. The study of PDOs here is similar in spirit to the study of PDOs on the $n$-torus [30, Chapter 5.3].

We will consider symbol classes, for $m\in\mathbb{R}$

[TABLE]

Here $C_{\beta}>0$ is a positive constant. We denote $S(0)$ simply by $S$. The symbols are quantised by the formula $a(x,hD)u(x)=\mathcal{F}_{h}^{-1}a(x,\xi)\mathcal{F}_{h}u$, where

[TABLE]

is the semiclassical Fourier transform. Now a symbol $a\in S(m)$ defines a map $a(x,hD):\mathscr{S}(\mathbb{R}^{n})\to\mathscr{S}(\mathbb{R}^{n})$ and by duality $a(x,hD):\mathscr{S}^{\prime}(\mathbb{R}^{n})\to\mathscr{S}^{\prime}(\mathbb{R}^{n})$, where $\mathscr{S}(\mathbb{R}^{n})$ are Schwartz functions. For $a\in S$, by standard theory we have $a(x,hD):L^{2}(\mathbb{R}^{n})\to L^{2}(\mathbb{R}^{n})$ uniformly in $h$. Our symbols will satisfy an additional invariance relation under $\varphi$, for $(x,\xi)\in\mathbb{R}^{2n}$

[TABLE]

From now on for simplicity we assume $a=a(\xi)\in S$ and satisfying (2.13), so $a=a\circ\mathcal{R}$. Then we have

Proposition 2.5**.**

The following properties hold for $\Phi$ a semiclassical PDO in $\mathbb{R}^{n}$ with symbol $a=a(\xi)\in S$, satisfying $a=a\circ\mathcal{R}$:

$\varphi^{*}\Phi=\Phi\varphi^{*}$.

• 2.

$P\Phi=\Phi P$* for $P$ a constant coefficient differential operator.*

Proof.

For the first item above, we have by definition and the change of coordinates $y^{\prime}=\varphi(y)$

[TABLE]

where $\mathcal{F}_{h}$ denotes the semiclassical Fourier transform. This further implies, after a change of coordinate $\xi^{\prime}=\mathcal{R}^{-1}(\xi)$ and using $a\circ\mathcal{R}=a$,

[TABLE]

which is interpreted as $(2\pi h)^{n}\varphi^{*}(\Phi u)(x)$.

For the second point, simply recall that $\mathcal{F}_{h}(D^{\alpha}u)=\frac{\xi^{\alpha}}{h^{|\alpha|}}\mathcal{F}_{h}(u)$, where $D=-i\partial$ and $\alpha$ is any multiindex. The proof then follows from a straightforward computation. ∎

The first conclusion of Proposition 2.5 says that $a(hD)u$ is $\varphi$-invariant if $u$ is so, if we assume $a\in S$ satisfies the invariance property (2.13). In this sense, we may study the mapping properties of $a(hD)$ on $L^{2}(\mathcal{C}_{\varphi})$:

Proposition 2.6**.**

The symbol $a=a(\xi)\in S$ satisfying (2.13) induces a map $a(hD):L^{2}(\mathcal{C}_{\varphi})\to L^{2}(\mathcal{C}_{\varphi})$.

Proof.

We follow the method of [30, Theorem 5.5]. Assume w.l.o.g. that $\Omega=\mathbb{R}^{n-1}$ and let $u\in L^{2}(\mathcal{C}_{\varphi})$. Then by a computation similar to the Proposition above, we obtain

[TABLE]

where we write

[TABLE]

We use the notation $\mathbbm{1}_{S}$ for the characteristic function of a set $S$. Now we claim that for $|k|\geq 2$

[TABLE]

as $h\to 0$, with a constant uniform in $k$. To prove this, notice that for any $N\in\mathbb{N}$ we have

[TABLE]

Using this formula, we may write $A_{k}=\mathbbm{1}_{\mathbb{R}^{n-1}\times[0,L)}(\varphi^{-k})^{*}\widetilde{A}_{k}\mathbbm{1}_{\mathbb{R}^{n-1}\times[0,L)}$, where

[TABLE]

Here we introduced

[TABLE]

Also, we write $\chi\in C^{\infty}$ for the cut-off such that $\chi=1$ near $\mathbb{R}^{n-1}\times[0,L]$ and zero outside $\mathbb{R}^{n-1}\times[-L,2L]$. Now by [30, Theorem 4.20] we may write $\widetilde{A}_{k}=b_{k}(x,hD)$ for a decaying symbol $b_{k}$ and then the boundedness of $b_{k}(x,hD)$ on $L^{2}(\mathbb{R}^{n})$ gives the claim. The main result then follows from the expansion (2.14). ∎

Now by using Proposition 2.6 and using the standard theory on $\mathbb{R}^{n}$, we may obtain the usual properties of semiclassical measures on the mapping torus $\mathcal{C}_{\varphi}$: existence under an $L^{2}$ bound, properties of the support and invariance under flow if a suitable equation is satisfied.

3. Billiard dynamics on polyhedra

In this section we discuss dynamical properties of the billiard flow on polyhedra. In the first two parts, we give basic definitions of the objects under study and revise the known results. In the third part, we prove a decomposition of non-singular directions on arbitrary polyhedra into finitely many tubes, i.e. the “finite tube condition”. Finally, we prove a property of maximal periodic tubes in a polyhedron that generalises the Property (P2) from the introduction.

3.1. Billiard dynamics preliminaries

Let $P\subset\mathbb{R}^{n}$ be a polyhedron and let us define the singular set $\mathcal{S}$ as the union of all $(n-2)$-dimensional “edges” of the polyhedron - in other words, the $(n-2)$-skeleton of $\partial P$. This represents the higher dimensional analogues of what are “corners” in polygonal billiards. Let $\mathcal{S}\subset U\subset P$ be an open neighbourhood of the singular set. The billiard flow on $P$ is the usual one, a particle (or a point mass) travels to $P\setminus\mathcal{S}$ with unit velocity, and then on striking one of the faces of the boundary, instantaneously changes direction according to the law of light reflection and continues along the reflected line. Trajectories which strike a singular point stop right there - such trajectories are also called singular. If the above does not happen, then the motion is determined for all time. Now, we introduce some notations and definitions, where we largely follow the exposition in [14].

Let us denote $\Gamma:=\partial P$, and let $T\Gamma$ be the set of all unit tangent vectors with base points in $\Gamma$ and which are directed inside $P$. Now, define the regular (non-singular) part of $T\Gamma$ as $T\Gamma_{1}:=\{x\in T\Gamma:\text{ the forward orbit of }x\text{ never hits }\mathcal{S}\}$. Denote by $f$ the first return (Poincaré) map of the billiard flow to the set $T\Gamma$. Then $f$ and its iterates are defined and smooth everywhere except for the vectors whose billiard orbits hit $\mathcal{S}$, which is a set of measure zero. Let the polyhedron $P$ have $l$ faces denoted by $\mathcal{F}_{1},\dotso,\mathcal{F}_{l}$. Define $\Sigma^{+}_{l}:=\{1,2,...,l\}^{\mathbb{N}}$ to represent the set of all forward strings for the symbolic dynamics of the billiard flow.

Also, given a trajectory starting from $x\in T\Gamma_{1}$, the symbolic string for the forward orbit is given by $w(x)$, defined by $w(x)_{i}=j$ iff the basepoint of $f^{i}(x)$ lies in $\mathcal{F}_{j}$. This gives us the symbolic encoding $\Sigma^{+}_{P}:=\{w\in\Sigma^{+}_{l}:\exists x\in T\Gamma_{1}\text{ such that }w=w(x)\}$, i.e., the set of all possible observable infinite strings. For each such possible string $w\in\Sigma^{+}_{P}$, we define $X(w):=\{x\in T\Gamma_{1}:w(x)=w\}$. In other words, $X(w)$ represents all tangent vectors whose billiard trajectories have the same symbolic representation $w$. An arbitrary element of $X(w)$ is denoted by $x(w)$.

We also take the space to make the following important distinction between rational and irrational polyhedra. As previously remarked, most of the literature focuses on dimension $n=2$, where the definition of rational (irrational) polyhedra can simply be given in terms of the rationality (irrationality, respectively) of the angles at the vertices, namely, the polygon is called rational if all its angles are rational, otherwise it is called irrational. For reasons which will become clear in the course of this paper, many statements which we will concern ourselves with are easier to prove for rational polygons than irrational ones. In dimensions $n\geq 3$, there does not seem to be a standard definition of rational/irrational polyhedra. We use the following definition (see [3, 4]):

Definition 3.1**.**

Let $P\subset\mathbb{R}^{n}$ be a polyhedron and let $\rho_{i}$ represent the linear reflection determined by the $i^{\text{th}}$-face of $P$. Then $P$ is called rational if the group $G$ generated by the $\rho_{i}$ is finite, otherwise the polyhedron is called irrational.

A crucial geometric difference exists between rational/irrational polyhedra, which we explain below. Let us first recall the well-known tool or method of unfolding a trajectory. Let $\gamma$ be a billiard trajectory of $x\in T\Gamma_{1}$. Starting from an initial point, we follow $\gamma$ until it strikes a face, say, $\mathcal{F}_{1}$ of $P$. Then we reflect $P$ about $\mathcal{F}_{1}$, and keep following $\gamma$ inside the reflected polyhedron until it strikes another face, whence we reflect the polyhedron again. Continuing this process indefinitely gives a sequence $P^{\infty}=P^{0}:=P,P^{1},P^{2},....,P^{m},...$ of polyhedra which are skewered on the forward ray determined by $\gamma$. We call this object an infinite corridor or an unfolding along the ray $\gamma$. Obviously, on refolding the corridor the line folds back (immerses) to $\gamma$.

Remark 3.2**.**

We bring to readers’ attention that in Definition 3.1, the reflections in faces are taken to be linear maps, whereas in the definition of a corridor, we take them as affine maps. In the remainder of the paper, this distinction will be addressed in the notation.

We record the following property related to unfolding: a symbolic encoding uniquely determines a direction.

Lemma 3.3**.**

In any polyhedron $P\subset\mathbb{R}^{n}$, if $w(x)=w(y)$, then $x$ and $y$ are parallel vectors.

Proof.

Unfold the polyhedron $P$ along trajectories determined by $x$ and $y$. Unless they are parallel, the distance between their trajectories grows linearly, and cannot be contained in the same corridor. The first time they move into different corridors, their symbolic representations $w(x)$ and $w(y)$ must also differ. ∎

We introduce the notation $D:=(P\sqcup\sigma P)/\sim$ to denote the double of $P$, where $\sigma$ denotes a reflection of $\mathbb{R}^{n}$; $\sigma P$ is the copy of $P$ under this reflection and $\sim$ glues the boundaries of $P$ and $\sigma P$ by a pointwise identification of the corresponding faces in $\partial P$ and $\partial(\sigma P)$.

This space plays the role of Euclidean surfaces with conical singularities (ESCS) as outlined in [16]. The singular set $\mathcal{S}\subset D$ is now of codimension at most $2$ and the space $D_{0}:=D\setminus\mathcal{S}$ can be given a structure of an open Euclidean $n$-manifold. Observe that in certain cases, taking finite covers can further reduce the singular set. For example, in the case $n=3$ a polyhedron $P$ might contain a vertex $v$ or an edge $e$ which are such that forming an $m$-fold cover $\tilde{P}$ around them converts them into removable singularities on $\tilde{P}$. This is the three-dimensional analogue of polygons having an angle of $\frac{2\pi}{m}$ at a vertex. In such cases, one can work with $\tilde{P}$ (and its double) instead of $P$ (and $D$), and the claim of Theorem 1.1 can be sharpened by replacing $\mathcal{S}$ by $\mathcal{S}\setminus\{v\}$ (or $\mathcal{S}\setminus\{e\}$, as the case may be). As an example, if $P$ is a rectangular parallelepiped, one can sharpen the statement of Theorem 1.1 to say that any neighbourhood of one of its vertices and the three edges emanating from it contain a certain fraction of the mass.222Note that we may write down eigenfunctions explicity as products of sines and cosines in this case, but there is no equidistribution in the high energy limit.

Finally, we introduce some more notation for the properties of the period, different types of tubes, etc. Given a set $U\subset\mathbb{R}^{n-1}$ and a local isometry $F:U\times\mathbb{R}\to D_{0}$ , we call $F(U\times\mathbb{R})$ an immersed tube or just a tube. We will often identify $F$ with its image $T:=F(U\times\mathbb{R})$. We also call $U$ the cross-section of the tube $T$. A lifted tube is the image of $T$ in the unit sphere bundle $SD_{0}$, determined by the unit vector in the positive direction of the tube.

Clearly, an immersed tube may be specified by a subset $Q$ of $T\Gamma_{1}$, consisting of parallel vectors whose base points form a convex set on one of the faces of the polyhedron. In dimension $n=3$ in particular, the tube is polygonal (elliptical) if there is an open polygon (ellipse) $V$ such that $V\subset Q\subset\overline{V}$. Given a point $x\in Q$, the image in an unfolding of the tube $T=T(x)$ that is generated by $Q$ is denoted by $T^{\infty}(x)\subset P^{\infty}\subset\mathbb{R}^{n}$. Depending on the use which will be clear from the context, sometimes we simply write $T^{\infty}$ instead of $T^{\infty}(x)$. Unless otherwise stated, we will always assume $Q=X(w)$ for some $w\in\Sigma_{P}^{+}$ and we will refer to such tubes $T(x)$ or $T^{\infty}(x)$ as maximal tubes, since they cannot be enlarged (see Proposition 3.5 below). We now formally define periodic tubes in our context.

Definition 3.4**.**

Let $U\subset\mathbb{R}^{n-1}$ be and let $F:U\times\mathbb{R}\to D_{0}$ be a local isometry. Then an immersed tube $F(U\times\mathbb{R})$ is called periodic if there is a positive number $L$ and a rotation $\mathcal{R}$ in $\mathbb{R}^{n}$ fixing the $\mathbb{R}$ direction, such that $F(x,t+L)=F(\mathcal{R}x,t)$ for all $x\in U$ and $t\in\mathbb{R}$.

We refer to $L$ as the length of the tube, which is also a period of the closed geodesic given by $F(\{x_{0}\}\times\mathbb{R})$, where $x_{0}\in U$ is the centre of mass of $U$. We call $F(\{x_{0}\}\times\mathbb{R})$ the central geodesic of $T$ and $\mathcal{R}$ the rotation associated to the tube $T$. Note that periodicity of an immersed tube $T$, as defined above, is not the same as saying that all the parallel billiard trajectories contained in the tube are individually periodic.

We gather a few basic properties of maximal periodic tubes in a proposition

Proposition 3.5**.**

The cross-section of a maximal tube $T$ is convex. Every trajectory on the boundary of $T$ comes arbitrarily close to the singular set.

Proof.

Let $\Omega$ be the cross-section of $T$. Consider two points $x,y\in\Omega$ and the unfolding $T^{\infty}$. By Lemma 3.3, trajectories determined by $x$ and $y$ hit the same face each time, so by the convexity of the faces of $P$, the trajectories determined by the segment $[x,y]$ always hit the interior of the faces.

For the second claim, assume the trajectory determined by $x\in\partial\Omega$ is at a distance $\varepsilon>0$ away from the singular set. Thus there is a neighbourhood $V$ of $x$, whose trajectories do not hit the singular set. Then the tube $T^{\prime}$ with the cross-section $\Omega^{\prime}:=\Omega\cup V$ contradicts the maximality of $T$. ∎

If $n=3$, we refer to the immersed periodic tube $T$ as a rational tube if $\mathcal{R}$ is a rotation in a rational multiple of $\pi$ and as an irrational tube otherwise. Of course, if the rotation $\mathcal{R}$ is by a rational multiple of $\pi$, then there exists some $\tilde{L}$ such that $F(x,t+\tilde{L})=F(x,t)$. In other words, each parallel billiard trajectory contained in the tube is now individually periodic. For irrational tubes, this is clearly not the case, and only the central geodesic is periodic.

Given a periodic trajectory $\gamma:\mathbb{R}\to P$, we will distinguish between a period and the minimal period, the latter being the minimal $T_{0}>0$ such that $\gamma(t+T_{0})=\gamma(t)$ for all $t$ and the former being any such $T_{0}$.

Remark 3.6**.**

Let us relate periodic tubes with almost periodic boundary conditions (2.6) defined in the previous section. Let $f\in C^{\infty}(D_{0})$ and take a periodic tube $F:\Omega\times\mathbb{R}\to D_{0}$ of length $L$ and rotation $\mathcal{R}$. Consider the pullback $F^{*}u$ to $\Omega\times\mathbb{R}$. By definition, we have

[TABLE]

which by Proposition 2.4 means that $F^{*}u$ satisfies the required condition.

3.2. A revision of known results

Here we collect the relevant preliminary results on polyhedral dynamics that will be of use later. For completeness, we have provided proofs in Appendix A. The proofs are quite instructive and we encourage the reader to go through them.

Consider a polyhedron $P$. As mentioned before, in higher dimensions, depending on the type of the polyhedron, an orbit generated by $x\in T\Gamma_{1}$ with a periodic symbol $w(x)$ might or might not be periodic. This is in contrast to the polygon case (see Theorem A.1). Now, observe that if the forward closure of the orbit of $x$, called $\gamma$ henceforth, does not intersect $\mathcal{S}$, then $\gamma$ can be “thickened” to form a tubular neighbourhood $T$ around $\gamma$ such that each trajectory in $T$ parallel to $\gamma$ also has the symbolic representation $w(x)$. It turns out that such tubes $T$ can themselves be periodic (see below).

With that in place, we have the following result proved in [14, Theorem 5]

Theorem 3.7**.**

Let $P\subset\mathbb{R}^{n}$ be an arbitrary convex polyhedron and $w\in\Sigma^{+}_{P}$ is a periodic sequence with minimal period $k$. The following hold:

1. (1)

There exists $x(w)$ so that $x(w)$ is periodic with minimal period $k$. 2. (2)

In addition, one of the following two cases holds:

1. (a)

There exists $q\geq 1$ such that all $y(w)\in X(w)\setminus x(w)$ are periodic with period $qk$ and the cross-section of the tube generated by $X(w)$ is an open polyhedron. 2. (b)

The set $X(w)$ generates a periodic tube with a convex cross-section $\Omega\subset\mathbb{R}^{n-1}$ and an associated isometry $\mathcal{R}_{0}\in O(n-1)$ keeping $\Omega$ invariant.333 By slightly abusing the notation, we extend here Definition 3.4 to include Poincaré maps in $O(n-1)$. Note that at the cost of doubling the length of the periodic tube, we may always assume this map is in $SO(n-1)$. 3. (3)

If $n=2$ or $3$ and $k$ is odd, then only the first case (2).(a) above can happen and $q=2$. 4. (4)

If $P$ is rational then only the first case (2).(a) above can happen.

Remark 3.8**.**

We remark that in [14] a precise statement of the theorem above was given only in the case $n=3$. In the case $n=3$, one additionally has that the periodic tubes in $(2)$.(b) are either polygonal or elliptical, or in other words $\Omega\subset\mathbb{R}^{2}$ is either a convex polygon or a disc. They respectively correspond to rational and irrational periodic tubes.

For a general value of $n$, the part $(2)$.(b) of the theorem needs to be modified and we have many more options for the cross-section $\Omega$ of a periodic tube. Moreover, $(3)$ now holds only for $n=2,3$, since in higher dimensions it is false that the product of three linear reflections in a general position is a linear reflection.

We include a proof of Theorem 3.7 in Appendix A. We now discuss a couple of explicit examples in the case $n=3$. As an example of when $(2)$.(a) above might occur, consider a right prism whose horizontal cross-section is an equilateral triangle. Consider a trajectory which lies on a plane perpendicular to the height of the prism, and strikes an equilateral triangular cross-section exactly at the mid-point of the three sides. This is periodic with minimal period $3$, as well as any vertical translate of such a trajectory, whereas any other $y(w)\in X(w)$ is periodic with minimal period $6$.

For an example when $(2)$.(b) might occur, it is enough to consider a regular tetrahedron and the closed orbit corresponding to the word $w=(abcd)$, where $a,b,c$ and $d$ encode the faces of the tetrahedron. Then it is possible to show that there is a unique closed orbit $x(w)$, that $X(w)$ generates an elliptical, irrational, periodic tube and the nearby parallel trajectories “come back” rotated by an irrational angle. See [2, Section 8] for explicit computations. Observe that although $X(w)$ generates an immersed solid torus $T$ with disc cross-section, the faces of $P$ intersect $T$ obliquely, and hence the footprint of $T$ on a face of $P$ looks like an ellipse.

Theorem 3.7 shows that the result for the two-dimensional polygonal case in Theorem A.1 does not generalise to higher dimensions in a completely straightforward way. In other words, there exist trajectories whose closure does not contain any singular point, but which are themselves not periodic - this might occur in the case of irrational polyhedra. Moreover, such trajectories are contained in periodic tubes.

However, the following result says that Theorem 3.7 contains all such possible obstructions:

Theorem 3.9**.**

Let $P\subset\mathbb{R}^{n}$ be a convex polyhedron and $w\in\Sigma_{P}^{+}$ non-periodic. Let $x\in X(w)$. Then the closure of the trajectory generated by $x$ intersects the singular set $\mathcal{S}$.

When $n=3$, the above result basically says that $X(w)$ is at best a codimension $1$ “strip”, which includes the case that $X(w)$ consists of a single point. We have included a proof of Theorem 3.9 in Appendix A. As a consequence we obtain the following dichotomy:

Corollary 3.10**.**

For any billiard trajectory $\gamma$ in a convex polyhedron $P$, either $\gamma$ is contained in an immersed periodic tube, or the closure of $\gamma$, $C_{\gamma}$, meets $\mathcal{S}$.

3.3. The finite tube condition

We begin this section by introducing a condition on the finiteness of periodic tubes missing $U_{\varepsilon}$, that will be relevant for our further discussion. Following [16], we have

Definition 3.11**.**

Let $D$ be the double of $P$ as defined above. A region $U\subset D$ is said to satisfy the finite tube condition if there exists a finite collection of periodic tubes $T_{i}$ for $i=1,\dotso,N$ for some $N$, such that any orbit that avoids $U$ belongs to some $T_{i}$.

We will sometimes simply say that $D$ satisfies the finite tube condition if every neighbourhood of $\mathcal{S}$ does so. In this section we prove the finite tube condition in full generality.

Theorem 3.12**.**

Let $P$ be a convex polyhedron and let $D$ be its double. Then, for any $\varepsilon>0$, the $\varepsilon$-neighbourhood $U_{\varepsilon}$ of the singular set $\mathcal{S}$ of $D$ satisfies the finite tube condition.

Proof.

Assume the finite tube condition is false, i.e. assume there are points $x_{1},x_{2},\dotso\in T\Gamma_{1}$ generating trajectories that stay in the future in $P\setminus U_{\varepsilon}$. By Theorem 3.7, for each $i\in\mathbb{N}$, the trajectories generated by $x_{i}$ belong to maximal immersed periodic tubes $T_{i}$, and we assume that the $T_{i}$ are distinct tubes.

By compactness and after relabelling we may assume $x_{i}\to x$ as $i\to\infty$. By continuity, the trajectory generated by $x$ is neither tangent to a face nor has basepoint in the singular set, and stays in $P\setminus U_{\varepsilon}$ in the future. Therefore $x$ belongs to a maximal periodic tube $T$ by Theorem 3.7, where $T$ is the image of a local isometry $F:\Omega\times\mathbb{R}\to D_{0}$. Note that for $i$ large enough, each $x_{i}$ points to a direction not parallel to $x$, by the fact that the tubes $T_{i}$ are distinct.

By Lemma 3.13, there are points $p_{1},\dotso,p_{N}\in\partial\Omega$ that are $\varepsilon/100$-dense, such that for some $t_{1},\dotso,t_{N}$ we have $F(p_{i},t_{i})\in\mathcal{S}$. By the same lemma, such points are uniformly recurrent: there is a $C>0$, such that any slice $F^{-1}(\mathcal{S})\cap\partial\Omega\times[x,x+C]$ of height $C$, contains $N$ points whose projections onto $\partial\Omega$ are $\varepsilon/100$-dense. In other words, $\partial T^{\infty}$ contains singular points whose projections onto the cross-section are $\varepsilon/100$-dense uniformly often.

Consider now $T(x_{i})$, i.e. the maximal tube generated by $x_{i}$. Note that such tubes contain immersed tubes with disc cross-section of radius $\varepsilon$. Taking $i\to\infty$, we obtain that $T^{\infty}(x_{i})$ intersects $\partial T^{\infty}$ in tubes of increasingly large height. By the construction, such tubes must contain singular points eventually. This contradicts the fact that $T^{\infty}(x_{i})$ do not contain singular points in the interior of their unfolding. ∎

3.4. Property 2’.

As is discussed in Section 1, the key to our proof are appropriate generalisations of Properties (P1) and (P2), together with our main estimate Theorem 5.4 (see Section 5 below). The generalised Property (P1) is contained in Theorems 3.7 and 3.9, so we discuss a suitable substitute of Property (P2) that we will call Property 2’. In order to state it, let $F:\Omega\times\mathbb{R}\to D_{0}$ be a maximal periodic tube and call $\Omega_{\varepsilon}$ the complement of the $\varepsilon$-neighbourhood of $\partial\Omega$.

[TABLE]

We will prove the Property (P2’) in this section. We begin by a lemma studying the formation of singularities on the boundary of a periodic tube; to this end assume $L$ is the length of $T$ and $\mathcal{R}$ the associated rotation.

Lemma 3.13**.**

Denote by $\pi_{1}:\mathbb{R}^{n-1}\times\mathbb{R}\to\mathbb{R}$ the projection to the first coordinate. Then the projections of singular points on $\partial T$ under $\pi_{1}$ are dense

[TABLE]

Moreover, singular points are uniformly recurrent: for every $\varepsilon>0$, there is an $L(\varepsilon)>0$ depending on $T$, such that for every $t\in\mathbb{R}$

[TABLE]

is $\varepsilon$-dense in $\partial\Omega$.

Proof.

By Proposition 3.5, for every $p\in\partial\Omega$ and for every $\gamma>0$ there is a $t\in\mathbb{R}$ such that $F(p,t)$ is $\gamma$-close to $\mathcal{S}$ (here we extended $F$ to $\overline{\Omega}\times\mathbb{R}$ as a continuous map to $D$). In fact, more is true and we can study the formation of singular points on the boundary of $T$. The idea of the proof of (3.1) is to use the periodicity of the tube and of singular points. Recall that the local isometry $F$ is invariant under the map $(x,t)\mapsto(\mathcal{R}x,t+L)$. Since this is a rigid motion, the nearby singular points in an unfolding can only approach $\partial T^{\infty}$ arbitrarily close, if there appears a singular point on $\partial T^{\infty}$.

More precisely, consider an unfolding $P^{\infty}=P^{0},P^{1},P^{2},\dotso$ along $T$ and denote by $(\partial P^{i})^{(n-2)}$ the $(n-2)$-skeleton of the boundary of $\partial P^{i}$. Define

[TABLE]

for a sufficiently large $k$ to be specified later. Let us define the forward orbit of $\mathcal{S}_{1}$ under the rotation $\mathcal{R}$ as

[TABLE]

By the observation above and by taking $k$ larger than the period of the symbol associated to $T$, we know that $\partial\Omega\subset\overline{\pi_{1}(\mathcal{S}_{\infty})}$. Also, observe that for any point $x\in\mathbb{R}^{n}$, since $\mathcal{R}$ is an isometry fixing $\partial T^{\infty}$

[TABLE]

Consider now $\pi_{1}(\mathcal{R}^{k_{i}}x_{i})\to y\in\partial\Omega$, for some $k_{i}\geq 0$ and $x_{i}\in\mathcal{S}_{1}$. By compactness of $\mathcal{S}_{\infty}$, we may assume without loss of generality $\mathcal{R}^{k_{i}}x_{i}\to x\in\partial T^{\infty}$, so $\pi_{1}(x)=y$. By compactness of $\mathcal{S}_{1}$ and since by (3.2) we have $d(x_{i},\partial T^{\infty})\to 0$, after re-labelling we may assume $x_{i}\to x^{\prime}\in\partial T^{\infty}\cap\mathcal{S}_{1}$. But then since $\mathcal{R}$ an isometry and by triangle inequality we have $\mathcal{R}^{k_{i}}x^{\prime}\to x$.

The previous argument shows that

[TABLE]

which proves the first claim. The second claim follows from the first claim, the periodicity of $T$ and the fact that $\mathcal{R}$ is an isometry of $\partial\Omega$.∎

Now we may prove this property in full

Lemma 3.14**.**

Property (P2’) holds for any $n\geq 2$.

Proof.

Define $\eta=\delta=\varepsilon/6$. Consider a point $p\in\partial\Omega$ closest to $p_{0}$, at a distance at most $2\varepsilon/3$. By Lemma 3.13, we can find a $z$ with $|z-z_{0}|\leq L(\varepsilon/6)$, such that there is a point $q\in\partial\Omega$ that is $\varepsilon/6$-close to $p$ and $F(q,z)\in\mathcal{S}$. Therefore, $(p_{0},z)$ is by triangle inequality $5\varepsilon/6$-close to $(q,z)$, and so its $\varepsilon/6$-neighbourhood is contained in $F^{-1}(U_{\varepsilon})$. ∎

Remark 3.15**.**

In other words, Property (P2’) is saying that the flow-out of $U_{\varepsilon}$ under the billiard flow contains, for each maximal periodic tube $T$, a conical neighbourhood of the tube direction in a neighbourhood of $\partial T$.

4. A quantitative estimate on periodic tube lengths

In this section we prove our main quantitative estimate on the lengths of periodic tubes. For this, we first establish an estimate on the angle of intersecting tubes. A new feature compared to the polygonal case, is that these estimates depend on dynamical invariants of associated rotations, that we are about to define.

We start with the case $n=2$. Given an irrational rotation $\mathcal{R}\in SO(2)$, $\varepsilon>0$ and a radius $r>0$, we write $N=N(\mathcal{R},\varepsilon,r)\in\mathbb{N}$ for the smallest positive integer such that for every $x$ lying on the circle of radius $r$, denoted by $S_{r}^{1}$, we have

[TABLE]

is $\varepsilon$-dense in $S^{1}_{r}$. If the rotation $\mathcal{R}$ is rational of minimal order $o$, we set $N(\mathcal{R},\varepsilon,r)=o$ for $\varepsilon$ sufficiently small. We specially set $N(\mathcal{R},\varepsilon,0)=1$. By the scaling property of the circle we obtain $N(\mathcal{R},\varepsilon,r)=N(\mathcal{R},\frac{\varepsilon}{r},1)$.

More generally, assume $\mathcal{R}\in SO(n)$. Given an $r\geq 0$, we will denote by $B_{r}\subset\mathbb{R}^{n}$ the ball of radius $r$; we define $B_{0}=\{0\}$. Clearly $\mathcal{R}$ acts on $\overline{B_{r}}$. Write $\{\mathcal{O}_{i}(r):i\in I(r)\}$ for the closed minimal444Minimal in the sense that $\mathcal{O}_{i}(r)$ does not contain any non-empty proper closed $\mathcal{R}$-invariant set. orbits of the action of $\mathcal{R}$ on $\overline{B_{r}}$; the $\mathcal{O}_{i}(r)$ partition $\overline{B_{r}}$ into disjoint pieces.

Definition 4.1**.**

Given an $\varepsilon>0$, we introduce

[TABLE]

Note that the orbits $\mathcal{O}_{i}(r)$ are in fact tori, as can be seen by conjugating $\mathcal{R}$ to a canonical form, which is a direct sum of rotations in $SO(2)$ with possibly some further $\pm 1$’s at the diagonal, and looking at the conjugated action on a suitable torus (cf. discussion before Proposition 2.4). Note also that $N=N(\mathcal{R},\varepsilon,r)$ in (4.1) is well-defined: the minimum always exists by the minimality and compactness of orbits and the maximum is well-defined since there are only finitely many orbit types (up to scaling and translation), by the diagonalisation argument above. In fact, for sufficiently small $\varepsilon>0$, the maximum is achieved by a “typical orbit”, i.e. the one with non-zero components corresponding to each rotation in $SO(2)$ or $\pm 1$ on the diagonal. Note that we also have the scaling property $N(\mathcal{R},\varepsilon,r)=N(\mathcal{R},\frac{\varepsilon}{r},1)$. Given a convex polyhedron $P$, we define $N(\mathcal{R},\varepsilon):=N(\mathcal{R},\varepsilon,\operatorname{diam}P)$.

Intuitively, Definition 4.1 and (4.1) give us a smallest positive integer $N$ such that the orbit of any point $x\in\overline{B_{r}}$ is $\varepsilon$-dense in the corresponding closed minimal orbit $\mathcal{O}_{i}(r)$. It is clear that this more general definition agrees with the one given for $n=2$. We will need these notions when quantifying the formation of the singular set on the boundary of a periodic tube. To be more precise, consider a maximal tube $T=F(\Omega\times\mathbb{R})$ of length $L$ and associated rotation $\mathcal{R}$. We will use the number $N=N(\mathcal{R},\varepsilon,\operatorname{diam}P)$ in Definition 4.1 above to say that the intersection of a slice of $T$ of height $NL$ with the singular set

[TABLE]

is $\varepsilon$-dense when projected orthogonally onto $\partial\Omega$, for any $x$. This will be crucial to exclude the intersection of long maximal tubes in the angle estimate Lemma 4.3 below. In other words, $NL$ quantifies the second part of Proposition 3.5 for periodic tubes.

Example 4.2**.**

Here we consider a few examples shedding some light on Definition 4.1. Firstly, since $SO(1)$ is trivial, we have $N=1$ in the case $n=1$ always. Then, if $\mathcal{R}\in SO(n)$ is rational, i.e. it is of finite order $o$, then the orbits $\mathcal{O}_{i}(r)$ consist of up to $o$ points and one easily sees that $N(\mathcal{R},\varepsilon,r)=o$ for $\varepsilon$ small enough and $r\neq 0$. Finally, if we take $\mathcal{R}\in SO(2)$ to be a rotation in an irrational angle $\alpha$, the speed at which the orbit of a point fills the circle depends on the arithmetic properties of $\alpha$ – see [21, Chapter 4.3.] for precise results. We also remark that we could have taken in (4.1) the $\mathcal{R}$ to be a reflection – then for $\varepsilon$ small enough, we would have $N(\mathcal{R},\varepsilon,r)=2$ since the orbits are of size $1$ or $2$.

Let now $D$ be the doubled polyhedron as before and let $U_{\varepsilon}$ be the $\varepsilon$-neighbourhood of $\mathcal{S}$. By Theorems 3.7 and 3.9, any trajectory not hitting $U_{\varepsilon}$ is contained in an immersed periodic tube with a certain convex cross-section. Consider any immersed tube $T$ of length $L$ with cross-section being a disc of radius $\frac{\varepsilon}{10}$. Recall that these tubes are given by immersing solid cylinders $F:B_{\frac{\varepsilon}{10}}\times\mathbb{R}\to D_{0}$. By definition, $F$ extends to an immersion $F:B_{\varepsilon}\times\mathbb{R}\to D_{0}$. Then we have the following lower bound on the angle of intersections of such periodic tubes in $D$.

Lemma 4.3**.**

Let $T_{1}$ and $T_{2}$ be two intersecting, periodic tubes of radii $\varepsilon/10$, lengths $L_{1},L_{2}$ and associated rotations $\mathcal{R}_{1},\mathcal{R}_{2}$, with geodesics $\gamma_{1}$ and $\gamma_{2}$ in their centre respectively, such that $\gamma_{1}$ and $\gamma_{2}$ do not intersect $U_{\varepsilon}$. More precisely, assume that a geodesic $\gamma_{1}^{\prime}$ lying in $T_{1}$ parallel to $\gamma_{1}$ intersects a geodesic $\gamma_{2}^{\prime}$ lying in $T_{2}$ parallel to $\gamma_{2}$ at a point $m\in D_{0}$, at an angle $\alpha$. Assume also $\gamma_{1}^{\prime}$ is not parallel to $\gamma_{2}^{\prime}$. Then the following bounds hold

We have

[TABLE]

• 2.

In particular, if $n=2$ then

[TABLE]

• 3.

In particular, if $P$ is rational and $o_{1}$ and $o_{2}$ are the finite orders of $\mathcal{R}_{1}$ and $\mathcal{R}_{2}$

[TABLE]

Proof.

Enlarge $T_{1}$ to an immersed tube $T_{1}^{\prime}$ of radius $\frac{3\varepsilon}{5}$, not intersecting $U_{\varepsilon/5}$, such that $T_{1}^{\prime}$ intersects $\gamma_{2}$. This is possible since the distance between $\gamma_{1}$ and $\gamma_{2}$ is at most $\frac{\varepsilon}{5}$. Similarly, extend $T_{2}$ to an immersed tube $T_{2}^{\prime}$ intersecting $\gamma_{1}$. By a slight abuse of notation, we will sometimes identify the tubes $T_{j}^{\prime}$ with their corresponding unfoldings $(T_{j}^{\prime})^{\infty}$ in $\mathbb{R}^{n}$ for $j=1,2$. Following this rule, let $p\in\gamma_{2}$ the unique, closest point to $\gamma_{1}$. It is unique since $\gamma_{1}$ and $\gamma_{2}$ are not parallel by the hypothesis.

Consider the strip $S_{2}$ contained in $(T_{2}^{\prime})^{\infty}$, determined by $\gamma_{1}$ and the direction of the tube $T_{2}^{\prime}$ and denote its width by $h_{2}$. Consider the points $p_{1}$ and $p_{2}$ defined as intersections of the strip $S_{2}$ with the boundary of the cross-section of $(T_{2}^{\prime})^{\infty}$ going through $p$ and consider the midpoint $H$ of the segment $p_{1}p_{2}$ (see Figure 1). By the triangle inequality applied to the triangle $p_{1}pH$, we have

[TABLE]

We introduce the parameter $l_{1}:=\frac{h_{2}}{\sin\alpha}$, i.e. the length of the portion of $\gamma_{1}$ lying in the strip $S_{2}$ (note that $\alpha\neq 0$). Introduce the notation $T_{1}^{\prime\prime}$ and $T_{2}^{\prime\prime}$ for the maximal periodic tubes corresponding to $T_{1}$ and $T_{2}$.

We treat the case $n=3$ first. Assume that $l_{1}>N(\mathcal{R}_{1},\varepsilon/5)L_{1}$ for the sake of contradiction. We may choose a segment $\gamma_{1}^{\prime\prime}$ of length $N(\mathcal{R}_{1},\varepsilon/5)L_{1}$ parallel to $\gamma_{1}$, lying on the boundary of $(T_{1}^{\prime\prime})^{\infty}$, such that $\gamma_{1}^{\prime\prime}$ is at a distance of (say) at least $\varepsilon/5$ to $\partial(T_{2}^{\prime\prime})^{\infty}$. This can be done e.g. by looking at $S_{2}\cap\partial(T_{1}^{\prime\prime})^{\infty}$, i.e. translating $\gamma_{1}$ in the strip $S_{2}$ until we hit the boundary of $(T_{1}^{\prime\prime})^{\infty}$. But by Lemma 3.13 below (or the proof of Theorem 3.7), there exists a singular point lying on the boundary of each slice of the periodic tube $(T_{1}^{\prime\prime})^{\infty}$ of height $L_{1}$. By the invariance of the singular set under the map $(x,t)\mapsto(\mathcal{R}_{1}x,t+L_{1})$ and by the definition of $N(\mathcal{R}_{1},\varepsilon/5)$, such a singular point leaves a $\varepsilon/5$-dense trace when projected to the boundary of the cross-section $\Omega_{1}^{\prime\prime}$ of $T_{1}^{\prime\prime}$, the projections taken over the union of $N(\mathcal{R}_{1},\varepsilon/5)$ copies of periodic slices in $(T_{1}^{\prime\prime})^{\infty}$ of height $L_{1}$. Therefore, the $\varepsilon/5$-neighbourhood of the segment $\gamma_{1}^{\prime\prime}$ must contain singular points, which also lie in the interior of $(T_{2}^{\prime\prime})^{\infty}$, giving a contradiction.

Therefore, we must have $l_{1}\leq N(\mathcal{R},\varepsilon/5)L_{1}$ and by an analogous argument, we obtain $l_{2}:=\frac{h_{1}}{\sin\alpha}\leq N(\mathcal{R},\varepsilon/5)L_{2}$. We combine these two inequalities together with (4.5) to get the estimate 1. for $n=3$.

For an arbitrary value of $n$, the proof is the same up to the point where we consider the intersection of the singular set with $\partial(T_{1}^{\prime\prime})^{\infty}$. Denote by $\Omega_{1}^{\prime\prime}$ the cross-section of $T_{1}^{\prime\prime}$. Since $\Omega_{1}^{\prime\prime}$ is invariant under $\mathcal{R}_{1}$, we know that $\mathcal{R}_{1}$ acts on $\partial\Omega_{1}^{\prime\prime}$ and we may write

[TABLE]

where $\mathcal{O}_{i}$ are minimal closed orbits of the action of $\mathcal{R}_{1}$ restricted to $\partial\Omega_{1}^{\prime\prime}$. By the proof of Lemma 3.13, if we denote by $\mathcal{T}_{1}$ the projection to $\partial\Omega_{1}^{\prime\prime}$ of the singular points in a slice of height $L_{1}$ in $\partial(T_{1}^{\prime\prime})^{\infty}$, the projections of all singular points on $\partial(T_{1}^{\prime\prime})^{\infty}$ are given by the forward orbit of $\mathcal{T}_{1}$ under $\mathcal{R}_{1}$. By the result of Lemma 3.13, this means that for each point $x\in\partial\Omega_{1}^{\prime\prime}$ and a $\delta>0$, there exists a $y\in\mathcal{T}_{1}$ whose forward orbit under $\mathcal{R}_{1}$ is $\delta$-close to $x$. In other words, the union of orbits $\mathcal{O}_{i}$ containing singular points is dense in $\partial\Omega_{1}^{\prime\prime}$. Observe also that $\Omega_{1}^{\prime\prime}\subset B_{\operatorname{diam}P}$, so the $(\varepsilon/5+\delta)$-neighbourhood of the $\gamma_{1}^{\prime\prime}$ as constructed in the previous paragraph, for any $\delta>0$, will contain singular points under the assumption that $l_{1}>N(\mathcal{R}_{1},\varepsilon/5)L_{1}$. But for a small enough $\delta>0$, the $(\varepsilon/5+\delta)$-neighbourhood of $\gamma_{1}^{\prime\prime}$ is contained in $T_{2}^{\prime\prime}$, giving a contradiction. This completes the missing step in the proof above, and proves the estimate 1. in full.

If $n=2$, we have $N(\mathcal{R},\varepsilon)=1$ always since $SO(1)$ is trivial, so (4.3) is a consequence of (4.2). If $P$ is rational, for $\varepsilon>0$ sufficiently small we have $N(\mathcal{R},\varepsilon)=o$ where $o$ is the order of the element $\mathcal{R}$ (cf. Example 4.2), so (4.4) is a consequence of (4.2). ∎

Remark 4.4**.**

By the argument in the second to last paragraph in the proof of Lemma 4.3, we deduce that one may take $L(\varepsilon)=N(\mathcal{R},\varepsilon)L$ in Lemma 3.13.

We are in a position to prove the main asymptotic estimate on the lengths of closed orbits missing an $\varepsilon$-neighbourhood $U_{\varepsilon}$ of the singular set $\mathcal{S}$ in $D$. Note that by Theorem 3.12, the number of periodic tubes whose central geodesic misses $U_{\varepsilon}$ is finite – denote this number by $M(P,\varepsilon)=M(\varepsilon)\in\mathbb{N}_{0}$.

Theorem 4.5**.**

Let $\varepsilon>0$. Enumerate the immersed periodic tubes whose central geodesics do not hit $U_{\varepsilon}$ by $\{T_{i}:i=1,\dotso,M(\varepsilon)\}$. Denote the length of $T_{i}$ by $L_{i}$ and the associated rotation by $\mathcal{R}_{i}$. Then, there exists a universal constant $C=C(n)>0$, such that for $\varepsilon$ small enough

[TABLE]

In particular, if $P$ is rational or $\dim P=2$ we have as $\varepsilon\to 0$

[TABLE]

Proof.

We first consider the case $n=3$. The idea of the proof is to partition the phase space into small distinct volumes around periodic tubes and to estimate these suitably using Lemma 4.3. For this purpose, we introduce the sets

[TABLE]

where $\dot{\gamma}_{i}(x)$ is the unit speed of the geodesic at $x$ in the direction of the immersed tube $T_{i}$. The distance $|\cdot|$ is taken to be the spherical distance on the sphere $S_{x}D_{0}$ over $x$. We may take $T_{i}$ to have a disc cross-section of radius $\varepsilon/10$. Observe that the tubes are disjoint and have no self-intersections by Lemma 4.3. Then $V_{i}$ have volume equal to the one of a tube of radius $\varepsilon/10$ and length $L_{i}$ times a spherical cap of spherical radius $\frac{4\varepsilon/5}{2N(\mathcal{R}_{i},\varepsilon/5)L_{i}}$. Thus we may estimate the volume of $V_{i}$ in phase space

[TABLE]

Here $C^{\prime}>0$ is a constant changing from line to line. The last inequality in (4.9) follows by observing that the lengths of closed orbits in $P$ are bounded from below by some $c(P)>0$, so for $\varepsilon$ small enough we may apply $\sin t\geq t/2$ valid for $t\in[0,\frac{\pi}{3}]$.

Now taking the disjoint sets $V_{i}$ as in (4.8), for $\varepsilon$ small enough and $i=1,\dotso,M(\varepsilon)$, summing the volumes and using the estimate (4.9), we obtain the main asymptotic estimate (4.6). In the general case, the proof is completely analogous to the proof above and we will omit it. If $P$ is rational or $\dim P=2$, the constants $N(\mathcal{R},\varepsilon)$ are explicitly known (see Example 4.2), so one obtains (4.7). ∎

Remark 4.6**.**

In the case of polygons, i.e. $n=2$, this lemma shows the cylinder (finite tube) condition of [16] directly. More precisely, for the sum on the left of (4.7) to be finite, we would need it to contain finitely many elements, since for the case $n=2$ every summand is exactly equal to one. Note that for $n=3$, we do not immediately get the analogous claim since we have the $\frac{1}{(N(\mathcal{R}_{i},c\varepsilon))^{2}L_{i}}$ “weights” in (4.6). However, one may use a slightly different argument (see below).

Using the angle estimate in Lemma 4.3, we may re-prove the finite tube condition obtained in Theorem 3.12.

Proof of Theorem 3.12 using the angle estimate.

To prove the claim, let us assume the contrary. That is, assume there are infinitely many periodic tubes $T_{i}$ of length $L_{i}$ with associated rotations $\mathcal{R}_{i}$, for $i\in\mathbb{N}$, that contain trajectories missing $U_{\varepsilon}$. Consider such trajectories $\gamma_{i}$ lying in $T_{i}$ and missing $U_{\varepsilon}$, generated by some $(x_{i},\theta_{i})\in SD_{0}$. Since there are infinitely many distinct tubes, due to the compactness of the phase space we have that a subsequence, still called $(x_{i},\theta_{i})$ after relabelling, converges to $(x,\theta)\in SD_{0}$. By a straightforward continuity argument, the future trajectory $\gamma$ determined by $(x,\theta)$ also stays in $D_{0}\setminus U_{\varepsilon}$. By Theorems 3.7 and 3.9 such a geodesic $\gamma$ is contained in a periodic, limiting tube $T$ of length $L$ and with the associated rotation $\mathcal{R}$.

Consider the immersed tubes $T_{i}^{\prime},T^{\prime}$ of radius $\varepsilon/10$ centred at $x_{i},x$ and with directions given by $\theta_{i},\theta$ (directions of $\gamma_{i},\gamma$), respectively, for any $i\in\mathbb{N}$. If $\alpha_{i}=\angle(T_{i}^{\prime},T^{\prime})=\angle(\theta_{i},\theta)$ is the angle of these intersecting tubes for each $i\in\mathbb{N}$, by applying Lemma 4.2 we obtain

[TABLE]

Here $N(\mathcal{R},\varepsilon/5)$ is the quantity defined in (4.1). But since $(x_{i},\theta_{i})\to(x,\theta)$, we must have $\alpha_{i}\to 0$, which is a contradiction to (4.10).∎

Remark 4.7**.**

In the statements above, for simplicity we mostly ignored the situation where the associated map $\mathcal{R}$ to a periodic tube is in $O(n)\setminus SO(n)$ and instead dealt with the tube of doubled length whose associated map is orientation preserving. However, the analogous statement of Lemma 4.2 holds for such tubes, and this may be applied to subsequent proofs of Lemma 4.2 and Theorem 3.12 to obtain analogous results.

5. A control estimate with an almost periodic boundary condition

In this section we prove the main control estimate, on a product space with an almost periodic boundary condition. Since we were unable to locate an appropriate reference, we start by establishing some control-theoretic preliminaries. As a general reference for the semiclassical analysis used in this section, the reader is referred to [30] or [11, Appendix E].

To this end, let $(M_{x},g_{x})$ be a compact Riemannian manifold with Lipschitz boundary. Denote by $-\Delta_{g}$ the positive-definite Laplace-Beltrami operator. We say a subset $A\subset M_{x}$ satisfies the geometric control condition or just (GCC) if every geodesic $\gamma$ in $M_{x}$ hits $A$ in finite time.

Lemma 5.1**.**

Let $\omega\subset M_{x}$ be an open neighbourhood of the boundary, satisfying (GCC). There exists a $C=C(g_{x},\omega)>0$, such that for any $s\in\mathbb{R}$ and any $v$ satisfying

[TABLE]

with $v\in H_{0}^{1}(M_{x})$ and $g\in H^{-1}(M_{x})$, we have the apriori estimate

[TABLE]

Proof.

Assume without loss of generality $g,v$ are real-valued. We split the proof to cases according to the value of $s\in\mathbb{R}$. We will use semiclassical defect measures in the case of large $s$.

Case 0: fixed $s$. Assume inequality (5.2) does not hold for this fixed $s$, so there is a sequence $v_{k}\in H_{0}^{1}(M_{x})$ with $\lVert{v_{k}}\rVert_{L^{2}(M_{x})}=1$ and $g_{k}\in H^{-1}(M_{x})$ with

[TABLE]

for any $k\in\mathbb{N}$. Then clearly $\lVert{g_{k}}\rVert_{H^{-1}(M_{x})}\to 0$ and $\lVert{v_{k}|_{\omega}}\rVert_{L^{2}(\omega)}\to 0$ as $k\to\infty$. Elliptic estimates give us that $\lVert{v_{k}}\rVert_{H^{1}(M_{x})}$ is bounded uniformly in $k$,555By elliptic estimates, we have $\lVert{v}\rVert_{H^{1}(M_{x})}\leq C^{\prime}(\lVert{g}\rVert_{H^{-1}(M_{x})}+\lVert{v}\rVert_{L^{2}(M_{x})})$ for some $C^{\prime}>0$, for $v$ and $g$ as in the statement of the Lemma. so by Rellich’s theorem we may assume $v_{k}\to v$ in $L^{2}(M_{x})$ after a possible re-labelling. By the assumptions, we have in the sense of distributions

[TABLE]

Elliptic regularity gives $v$ is $C^{\infty}$ in the interior of $M_{x}$ and therefore by the unique continuation principle $v\equiv 0$, which contradicts $\lVert{v}\rVert_{L^{2}(M_{x})}=1$.

Case 1: bounded $s$. Here we show that the constant $C=C(s)$ in (5.2) is locally bounded. The proof is by contradiction and is very similar to the previous case, so we omit it.

Case 2: $s<-\varepsilon<0$. In this case $1>\varepsilon>0$ is fixed. Integrating by parts, we see

[TABLE]

Estimating right hand side using the boundedness of the $H^{-1}\times H_{0}^{1}$ pairing and using $s<-\varepsilon$, we obtain $\lVert{v}\rVert_{L^{2}(M_{x})}\leq\frac{1}{\varepsilon}\lVert{g}\rVert_{H^{-1}(M_{x})}$, proving (5.2) for $s<-\varepsilon$.

Case 3: $s\to\infty$. Now it suffices to argue by contradiction and assume there is a sequence $s_{k}\to\infty$ with $C(s_{k})\to\infty$, as $k\to\infty$. Then, there is a sequence $v_{k}\in H_{0}^{1}(M_{x})$, with $\lVert{v_{k}}\rVert_{L^{2}(M_{x})}=1$ such that

[TABLE]

with $\lVert{g_{k}}\rVert_{H^{-1}(M_{x})}\to 0$ and $\lVert{v_{k}}\rVert_{L^{2}(\omega)}\to 0$ as $k\to\infty$. We introduce a (small) semiclassical parameter $h_{k}>0$ by $h_{k}^{2}:=\frac{1}{s_{k}}$. Then we have666Recall that in $\mathbb{R}^{n}$, $\lVert{u}\rVert_{H^{-1}_{h}}=\lVert{\langle{h\xi}\rangle^{-1}\widehat{u}(\xi)}\rVert_{L^{2}}$ and $\lVert{u}\rVert_{H^{-1}}=\lVert{\langle{\xi}\rangle^{-1}\widehat{u}(\xi)}\rVert_{L^{2}}$, where $\langle{\xi}\rangle=(1+|\xi|^{2})^{\frac{1}{2}}$. From definitions and since $g_{k}=o_{H^{-1}}(1)$, we have locally $\lVert{h_{k}^{2}g_{k}}\rVert_{H^{-1}_{h_{k}}}\leq h_{k}\lVert{g_{k}}\rVert_{L^{2}}=o(h_{k})$ for $h_{k}<1$ and as $k\to\infty$.

[TABLE]

So by [11, Theorems E.42, E.43, E.44], we know there is a semiclassical (Radon) measure $\mu$ on $S^{*}M_{x}$ associated to a subsequence $v_{k}$ (after re-labelling), such that it is invariant under the geodesic flow777In this proof we only work with the invariance of the semiclassical measure under the flow in the interior. For the invariance under the billiard (broken geodesic) flow, the reader is referred to [29]. of $g_{x}$ and

[TABLE]

Here $a\in C_{0}^{\infty}(T^{*}M_{x})$, compactly supported in the interior of $M_{x}$ and $a(x,hD)$ is a $\Psi$DO obtained by the semiclassical quantisation procedure on (the interior of) $M_{x}$.

Now, since $\lVert{v_{k}}\rVert_{L^{2}(\omega)}\to 0$ as $k\to\infty$, we have by (5.7) that $\mu=0$ on $\pi^{-1}(\omega)$, where   $\pi:S^{*}M_{x}\to M_{x}$ is the projection. Since $\mu$ is invariant by the geodesic flow and $\omega$ satisfies (GCC), we thus have $\mu\equiv 0$. But since we assumed no mass escapes to the boundary, we have that $\mu$ is a probability measure on $S^{*}M_{x}$ (cf. Proposition 6.1), which is a contradiction. ∎

Remark 5.2**.**

The case of Lemma 5.1 for $M_{x}=[0,a]$ was considered in [7], by using elementary means to prove the inequality directly. Note also that if $M_{x}\subset\mathbb{R}^{n}$, then any open neighbourhood of the boundary automatically satisfies the (GCC).

Remark 5.3**.**

In the theorem above, the condition that $\omega$ is a neighbourhood of the boundary is used to prevent the mass to escape from the interior (see [17] for a very similar condition, but in a different context). An example where Dirichlet eigenfunctions concentrate at the boundary is given by the upper hemisphere of $S^{2}\subset\mathbb{R}^{3}$ and eigenfunctions

[TABLE]

for $l\in\mathbb{N}$, of eigenvalue $\lambda=l(l+1)$. Here $c_{l}\sim l^{3/4}$ is the $L^{2}$-normalisation constant and $(\theta,\varphi)$ are the spherical coordinates. Then one may show that $u_{l}$ concentrates on the equator $z=0$. However, it seems that using the techniques as in [8, Theorem 8], one could deal with this and prove a more general version of Lemma 5.1, for $\omega$ any open set satisfying (GCC).

We are now in shape to prove the main theorem of this section, which generalises a result of N. Burq given in the Appendix B. As before, let $(M_{x},g_{x})$ be a compact Riemannian manifold with Lipschitz boundary.

Theorem 5.4**.**

Let $\varphi:M_{x}\to M_{x}$ be an isometry. Assume that $u\in H^{1}_{\operatorname{loc}}(M_{x}\times\mathbb{R})\cap C(\mathbb{R},H^{1}(M_{x}))$, such that $u(x,t+L)=u(\varphi(x),t)$ for all $(x,t)\in M_{x}\times\mathbb{R}$, for some $L>0$. Define $\mathcal{C}_{\varphi}:=M_{x}\times[0,L]/(x,L)\sim(\varphi(x),0)$ to be the mapping torus determined by $\varphi$, with the inherited Riemannian metric from $M_{x}\times\mathbb{R}$. Assume $u$ satisfies, for some $s\in\mathbb{R}$

[TABLE]

where $f\in H_{\operatorname{loc}}^{-1}(M_{x}\times\mathbb{R})\cap C(\mathbb{R},H^{-1}(M_{x}))$. Let $\omega\subset M_{x}$ be an open neighbourhood of the boundary satisfying (GCC) and assume $\omega$ invariant under $\varphi$. Denote the mapping torus over $\omega$ by $\omega_{\varphi}$. Then there exists a constant $C=C(M_{x},g_{x},\omega)>0$, such that the following observability estimate holds:

[TABLE]

Proof.

We use the theory of almost periodic functions outlined in Section 2.1. By Proposition 2.4 we know $\varphi$ is admissible, so by Lemma 2.3 we have

[TABLE]

is almost periodic. We will use the shorthand notation $u(t):=u(\cdot,t)\in L^{2}(M_{x})$ and $f(t):=f(\cdot,t)$. Therefore, there exists a countable set $\{\lambda_{n}\}_{n=1}^{\infty}\subset\mathbb{R}$ such that

[TABLE]

Recall the notation $\sim$ denotes a formal association of the series on the right hand side to the almost periodic function $u(t)$. The Fourier-Bohr transformation is given by

[TABLE]

where $\mathcal{M}\{\cdot\}$ denotes taking the mean value in $\mathbb{R}$. Therefore, we have $a(\lambda_{n};u)=:u_{n}(x)\in H_{0}^{1}(M_{x})$.

Now, since $\varphi$ an isometry so the pullback $\varphi^{*}$ commutes with $\Delta_{g_{x}}$, we have $f(x,t+L)=f(\varphi(x),t)$ for all $(x,t)$, so $f(t)$ is also almost periodic with range in $H^{-1}(M_{x})$. One can check that, using mapping properties of $-\Delta_{g_{x}}-s:H_{0}^{1}(M_{x})\to H^{-1}(M_{x})$ and integration by parts for the $\partial_{t}^{2}$ factor, together with the convergence of the Bohr transform (5.11)

[TABLE]

Similarly, one has $f(t)\sim\sum_{n}f_{n}(x)\cdot e^{i\lambda_{n}t}$ where $f_{n}(x)$ are the Fourier-Bohr coefficients of $f(t)$, and by uniqueness of Fourier-Bohr expansions of almost periodic functions, we have for all $n\in\mathbb{N}$

[TABLE]

Intuitively, the above takes the place of a Fourier series expansion for the periodic case. Now for every $n\in\mathbb{N}$ and some $C>0$, the estimate in Lemma 5.1 yields

[TABLE]

Now we sum the inequalities (5.12) and use Parseval’s identity for almost periodic functions (2.4)

[TABLE]

More precisely, we applied Parseval’s identity in the first line to $u$, in the third line to $f$ and $u|_{\omega\times\mathbb{R}}$. Here we also used the identity

[TABLE]

which follows upon restricting (5.10) to $\omega\times\mathbb{R}$ and recalling that $\varphi:\omega\to\omega$ is an isometry. We used (5.12) in the second line. In the first and last lines we also used the fact that $\lVert{u(t)}\rVert^{2}_{H^{s}(M_{x})}$ and $\lVert{f(t)}\rVert^{2}_{H^{s}(M_{x})}$ are periodic with period $L$ for any allowed $s$ (as $\varphi^{*}$ is an isometric isomorphism on $H^{s}(M_{x})$), so applying $\mathcal{M}\{\cdot\}$ gives their mean over the interval $[0,L]$. This finishes the proof. ∎

Remark 5.5**.**

It might be possible to prove the statement above by using a more direct approach of approximating $u$ by periodic functions on $M_{x}\times[0,kL]$ for $k$ large, applying Theorem B.2 and then taking a limiting procedure.

Remark 5.6**.**

It is interesting to compute the spectrum of $\mathcal{C}_{\varphi}$ in terms of the spectrum of $M_{x}$ and the map $\varphi$. In fact, as the pullback $\varphi^{*}$ and $-\Delta_{g_{x}}$ commute, and since $\varphi^{*}$ is orthogonal as a map on $L^{2}$, we may restrict $\varphi^{*}$ to eigenspaces $E_{\mu_{n}}$ of $-\Delta_{g_{x}}$ and thus choose an orthonormal Dirichlet eigenbasis of $M_{x}$, call it $e_{n}$, such that $-\Delta_{g_{x}}e_{n}=\mu_{n}e_{n}$ and $\varphi^{*}e_{n}=e^{i\nu_{n}L}e_{n}$. The reals $\nu_{n}$ are defined modulo $\frac{2\pi\mathbb{Z}}{L}$. We set for $n\in\mathbb{N}$ and $k\in\mathbb{Z}$

[TABLE]

and observe that $e_{n,k}$ descend to $\mathcal{C}_{\varphi}$. Also

[TABLE]

and one may show that $e_{n,k}$ form an orthogonal eigenbasis of $L^{2}(\mathcal{C}_{\varphi})$, with $\lVert{e_{n,k}}\rVert_{L^{2}}^{2}=L$. Observe that a function $u$ in (5.10), since it descends to $\mathcal{C}_{\varphi}$, satisfies $u_{n}(\varphi(x))=u_{n}(x)e^{i\lambda_{n}L}$. In fact, density of $e_{n}$ may be used to show $e^{i\lambda_{n}L}=e^{i\nu_{k_{n}}L}$ for some $k_{n}\in\mathbb{Z}$ and so $\lambda_{n}\in\nu_{k_{n}}+\frac{2\pi\mathbb{Z}}{L}$. This determines the spectrum of $\mathcal{C}_{\varphi}$ as advertised.

Given $u$ and $f$ as in Theorem 5.4, the argument above gives us a way to write $u=\sum_{n}u_{n}(x)e^{i\lambda_{n}t}$ and $f=\sum_{n}f_{n}(x)e^{i\lambda_{n}t}$. However, the $e^{i\lambda_{n}t}$ are in general not orthogonal with respect to $L^{2}(0,L)$, which would take us eventually to the operation $\mathcal{M}\{\cdot\}$. To prove (5.12), we would need a Parseval’s identity (cf. (2.4)). Therefore, in order to avoid the use of almost periodic theory, we would have to re-prove some essential parts of it.

Having the above estimate with almost periodic boundary conditions at hand, we can recover the control theoretic result with periodic boundary conditions proved and used in [16], which stems from the work of Burq-Zworski [7, 8]. More precisely, we have

Corollary 5.7**.**

If we choose $\varphi=\operatorname{id}$ in Theorem 5.4, we recover [16, Proposition 14] and [8, Proposition 6.1]. At the same time, we recover Theorem B.1 for the case $M_{y}=[0,L]$.

6. Proof of the Main Theorem

In this section, we give the proof of the Main Theorem 1.1. The proof is based on a careful analysis of the closed orbits missing a neighbourhood of the singular set, semiclassical measures and the control estimate from the previous section.

Throughout the rest of the paper, we will slightly abuse notation, and treat $U$ both as a neighbourhood of the singular set $\mathcal{S}$ in the double $D$ and in $P$. It will be clear from the proofs that the main eigenfunction concentration result will also be established for the double $D$, but it will readily imply Theorem 1.1. In fact, an eigenfunction $u_{n}$ on $P$ induces a $C^{\infty}$ eigenfunction on $D_{0}$ (denoted by the same letter), that is equal to $u_{n}$ on $P$ and to $-u_{n}\circ\sigma$ in $\sigma P$ for the Dirichlet Laplacian and $u_{n}\circ\sigma$ for the Neumann Laplacian. A lower bound on the mass of $u_{n}$ near $\mathcal{S}$ in $D$ translates to a similar bound in $P$.

Also, to clarify, one can define the Laplacian on $D_{0}:=D\setminus\mathcal{S}$, by taking the Friedrichs extension of the operator with domain $C_{c}^{\infty}(D_{0})$. It is known that this is self-adjoint with compact resolvent, so has discrete spectrum $\lambda_{j}$ going to infinity, and a complete orthonormal $L^{2}$-basis $u_{j}$ of Laplace eigenfunctions.

To begin the proof of our Main Theorem 1.1, let us assume to the contrary that there is no concentration in the neighbourhood $U$, that is, there exists a subsequence $u_{n}$ satisfying

[TABLE]

Let $\mu$ be an arbitrary semiclassical measure associated to the sequence $u_{n}$ in the standard way. That is, for any $a\in C^{\infty}_{0}(S^{*}D_{0})$ (unit cotangent bundle on $D_{0}=D\setminus\mathcal{S}$) we have

[TABLE]

Here $h_{n}^{2}=\frac{1}{\lambda_{n}}$ and $a(x,hD)$ is a choice of quantisation on $D_{0}$. We will look at the interaction of this semiclassical measure with the geodesic flow on $D$. To start with, we have the following result. Let $U_{0}:=U\setminus\mathcal{S}$ and denote by $\pi:S^{*}D_{0}\to D_{0}$ the footpoint projection.

Proposition 6.1**.**

The support of $\mu$ is disjoint from $\pi^{-1}(U_{0})$ and $\mu$ is a probability measure which is invariant under the geodesic flow.

Proof.

Suppose to the contrary that there is a point $q\in\operatorname{supp}\mu$ with $\pi(q)\in U_{0}$. Hence by definition, there exists a neighbourhood $V$ with $q\in V\subseteq\pi^{-1}(U_{0})$ so that $\mu(V)>0$. Now let $\phi\in C_{0}^{\infty}(S^{*}D_{0})$ be a function taking values in $[0,1]$ with $\phi=1$ on $V$ and $\operatorname{supp}\phi\subseteq\pi^{-1}(U_{0})$. From the nonnegativity of $\mu$ and the assumption on $V$, it follows that $\langle\mu,\phi\rangle>0$. Hence, using the definition of a semiclassical measure, we obtain the existence of a fixed positive constant $\gamma$ such that

[TABLE]

Now, we can decompose

[TABLE]

To estimate the right hand side further from above one can use a bound on the operator norm (see, e.g., [30, Theorem 5.1]):

[TABLE]

By the pseudolocality of PDOs (see, e.g., the Remark on p. 81, [30]) and eventually slightly enlarging $V$, one can assume that the second term on the right in (6.4) can be made negligible, whereas the first term on the right in (6.4) is bounded above as:

[TABLE]

where we have used the assumption (6.1). This contradicts (6.3).

Further, the invariance of $\mu$ under the geodesic flow is a well-known fact for semiclassical defect measures. For further discussion on semiclassical correspondence and Egorov’s theorem, we refer to [30, Chapter 11].

To prove that $\mu$ is a probability measure, one follows a similar line of reasoning as above, i.e. by taking an appropriate cut-off functions supported in $D\setminus U_{0}$ and using (6.2) and the fact that $u_{n}$ are $L^{2}$-normalised. Technically, one needs to argue in a similar fashion as above (using pseudolocality) to estimate the different terms. We emphasise that it is crucial for the arguments that we work away from the singular points. ∎

Now we begin with the proof of Theorem 1.1. Before tackling the general case, we first give an easier proof that works for rational polyhedra, and contains some preliminary ideas. This proof is essentially a higher dimensional analogue of [16].

Proof of Theorem 1.1 for Rational Polyhedra.

Let $\mu$ be a limit semiclassical measure associated to the eigenfunction sequence $u_{n}$ as above, and let $(x_{0},\zeta)\in{S^{*}D_{0}}$ be in the support of $\mu$. According to Proposition 6.1 and Corollary 3.10, the finite tube condition (see Theorem 3.12) implies that $x_{0}$ belongs to a tube periodic in the direction of $\zeta$. Thus, the support of $\mu$ is included in the union of the lifted maximal periodic tubes $T_{i}$, for $i=1,2,\dotso,N$. Furthermore, assume $U$ is the $\varepsilon$-neighbourhood $U_{\varepsilon}$ of $\mathcal{S}$ for an $\varepsilon>0$.

Choose one such tube $T:=T_{i}$ and choose its length $L$ such that every geodesic in $T$ is periodic with period $L$. By definition, $T$ is the image under a local isometry $F:\Omega\times\mathbb{R}\to D_{0}$, where $\Omega$ is a convex polygon. Denote the coordinates on $\Omega$ by $(x,y)$ and on $\mathbb{R}$ by $z$. Using this local isometry $F$, we pull back the eigenfunction $u_{n}$ to $T$. We now apply the black box concentration techniques of [7, 8, 16] to this setting. We choose an “inner shell” of thickness $\varepsilon$ (that is, the complement of a tubular neighbourhood of $\partial\Omega$ of thickness $\varepsilon$), denoted by $\Omega_{\varepsilon}\subset\Omega$, so that $\Omega_{\varepsilon}\times\mathbb{R}$ does not intersect $F^{-1}(U_{\varepsilon})$. Now consider a smooth cut-off function $\chi$ such that $\chi\equiv 1$ inside $\Omega_{\varepsilon}$ and $\chi=0$ outside $\Omega_{\varepsilon/2}$. Then, $v_{n}:=\chi u_{n}$ vanishes near the walls of the tube $T$. The $v_{n}$ are bounded in $L^{2}_{\operatorname{loc}}$, so there exists at least one associated semiclassical measure $\nu$ on $T$. Since $\mu$ is supported on a finite number of lifted tubes, there are only a finite number of directions in the support of $\nu$. So we can find a (constant-coefficient) semiclassical pseudodifferential operator $\Phi=a(hD)$ on $T$ that is microlocally $1$ in a neighbourhood of $\zeta=dz$ (the direction of $T$), but vanishes microlocally in a neighbourhood of every other direction in the support of $\nu$. Let $\Phi_{n}:=a(h_{n}D)$.

Now, consider the sequence of functions $w_{n}:=\Phi_{n}v_{n}$ on $T$. The semiclassical measure $\nu^{\prime}$ associated to this sequence is related to the one for the sequence $v_{n}$ by $\nu^{\prime}=|\sigma(\Phi)|^{2}\nu$, where $a=\sigma(\Phi)$ is the principal symbol of the operator $\Phi$. Thus, the support of $\nu^{\prime}$ is restricted to directions parallel to $dz$ and to the geodesics parametrised by $x,y$ such that $\chi(x,y)=1$.

Then, since $\Phi$ commutes with constant coefficient differential operators, we have on $\Omega\times[0,L]$

[TABLE]

We are now in a position to apply Theorem B.2 (or Theorem 5.4 for $\varphi=\operatorname{id}$), for $u=w_{n}$, $s=\lambda_{n}$, $f=f_{n}:=-\Phi_{n}((\Delta\chi)u_{n})-2\Phi_{n}(\partial_{x}\chi\partial_{x}u_{n})-2\Phi_{n}(\partial_{y}\chi\partial_{y}u_{n})$ and $\omega$ contained in the set $\{\chi=0\}$. By our choice of $\omega$, we have $\|w_{n}\|^{2}_{L^{2}(\omega\times[0,L])}\to 0$. Also, it is clear from the choice of our cut-off that since $\operatorname{supp}(\nabla\chi)$ is disjoint from $\operatorname{supp}\nu^{\prime}$, $\|f_{n}\|_{H^{-1}_{x,y}L^{2}_{z}(\Omega\times[0,L])}\to 0$ as $n\to\infty$. This implies that $\nu^{\prime}$, and hence $\nu$ does not have any mass in the direction of the tube $T$.

We may apply the above argument to each tube $T_{i}$ to obtain $\mu=0$, which contradicts the latter conclusion of Proposition 6.1 and finally proves the claim. ∎

Now we give the proof in the general case, that is, the case of irrational polyhedra. Observe that in this setting there are several issues that render the above proof invalid. Firstly, the boundary condition, as used in Theorem B.2, is not periodic, but almost periodic in the irrational case. We recall and emphasise that in a periodic tube, trajectories need not be individually periodic. Moreover, we improve the proof by relying on the Property (P2’), that is significantly weaker than the finite tube condition.

Proof of Theorem 1.1.

As before, let $\mu$ be the semiclassical measure associated to the eigenfunction sequence $u_{n}$, and let $(x_{0},\zeta)\in S^{*}D_{0}$ be in the suppport of $\mu$. According to Proposition 6.1 and Corollary 3.10, this generates a trajectory contained in a maximal periodic tube $T$ with the direction of $\zeta$, of length $L$ and associated rotation $\mathcal{R}$. By definition, there is a local isometry $F:\Omega\times\mathbb{R}\to T$, where $\Omega$ is a convex set, invariant under $\mathcal{R}$. Using this local isometry, by a slight abuse of notation we identify $\Omega\times\mathbb{R}$ with $T$ and $F^{*}u_{n}$ with $u_{n}$; we also pull back the measure $\mu$ to $T$. Write $x$ for the coordinate on $\Omega$ and $t$ for the $\mathbb{R}$ coordinate. Note that $F\circ\varphi=F$, where $\varphi(x,t)=(\mathcal{R}x,t+L)$; thus $\mu$ is invariant under $\varphi$. Recall $\Omega_{\varepsilon}\subset\Omega$ denotes the inner shell of thickness $\varepsilon$ inside $\Omega$, i.e. the complement of the $\varepsilon$-neighbourhood of $\partial\Omega$.

Let $\chi\in C_{0}^{\infty}(\Omega)$ be invariant by $\mathcal{R}$ and such that

[TABLE]

Denote $v_{n}:=\chi u_{n}$ and by $\mathcal{C}_{\varphi}$ the mapping torus determined by $\varphi$. Since $T$ does finitely many reflections in its one period, we have $\lVert{v_{n}}\rVert_{L^{2}(\mathcal{C}_{\varphi})}$ bounded uniformly in $n$, and the associated semiclassical measure is $\nu=\chi^{2}\mu$.

Take now $\Phi$ a constant coefficient semiclassical PDO in $\mathbb{R}^{n}$, microlocally cutting off near $dt$, constructed as follows. Using the Property (P2’), we may take its symbol $a(\xi)$ supported in a small cone $\Gamma\subset\mathbb{R}^{n}\setminus 0$ around $dt$, such that all lines in the direction of $\Gamma$ with basepoint $x\in\Omega_{\varepsilon/2-\eta}\setminus\Omega_{\varepsilon/2+\eta}$ hit the set $F^{-1}(U_{\varepsilon})$ in finite time (see Figure 2). Moreover, we may take $a(\xi)$ equal to $1$ near $dt$ and invariant under rotations around $dt$, i.e. so that we have $a\circ\mathcal{R}=a$.

Let $w_{n}:=\Phi_{n}(v_{n})$, where $\Phi_{n}=a(h_{n}D)$ and observe that by Proposition 2.5, $w_{n}$ descends to $\mathcal{C}_{\varphi}$. The semiclassical measure $\nu^{\prime}$ associated to $w_{n}$ satisfies $\nu^{\prime}=|a|^{2}\chi^{2}\mu$. By our choices of $\Phi$, $\chi$ and the invariance of $\mu$ under the geodesic flow, we have $\nu^{\prime}=0$ on \Big{(}\Omega\setminus\Omega_{\varepsilon/2+\eta}\Big{)}\times\mathbb{R}. By a computation and since $\Phi$ commutes with constant coefficient differential operators, we have on $\Omega\times[0,L]$ the following equation with an almost periodic boundary condition

[TABLE]

Choose $\omega\subset\Omega$ to be a small enough neighbourhood of $\partial\Omega$, for example $\omega=\Omega\setminus\Omega_{\varepsilon/2-\eta}$ would work. We are now in a position to apply Theorem 5.4 to $u=w_{n}$, $s=\lambda_{n}$, $f=f_{n}:=-\Phi_{n}((\Delta\chi)u_{n})-2\Phi_{n}(\nabla_{x}\chi\cdot\nabla_{x}u_{n})$ and $\omega$ as above to get888Note that a convex domain always has Lipschitz boundary, so we may apply the analytical results of previous section.

[TABLE]

By our choice of $\omega$, i.e. since $\chi=0$ on a neighbourhood of $\overline{\omega}$ and so $\nu^{\prime}=0$ near $\overline{\omega}$, we have $\|w_{n}\|_{L^{2}(\omega_{\alpha})}\to 0$. The $H^{-1}_{x}$ norm of the term for a fixed $t$ in the second bracket on the right hand side of (6.10) is bounded by a constant times the following expression, since $\Phi_{n}$ commutes with constant coefficient differential operators by Proposition 2.5

[TABLE]

This is further bounded by, for some $C>0$

[TABLE]

The $L^{2}_{t}$ norm of the expression in (6.11) over $[0,L]$ goes to zero as $n\to\infty$, since the semiclassical measure $\nu^{\prime}$ vanishes on \Big{(}\Omega_{\varepsilon/2}\setminus\Omega_{\varepsilon/2+\eta}\Big{)}\times\mathbb{R} and $\operatorname{supp}\nabla\chi\subset\Omega_{\varepsilon/2}\setminus\Omega_{\varepsilon/2+\eta}$. This shows that the right hand side of (6.10) goes to zero as $n\to\infty$, which implies $\nu^{\prime}\equiv 0$, and further implies $\mu\equiv 0$ since the choice of $(x_{0},\zeta)$ was arbitrary. But this is a contradiction with Proposition 6.1, which gives that $\mu$ is a probability measure. ∎

6.1. Another application

In proving Theorem 1.1, our main technical result was a version of a control theory estimate in [8], which also works for periodic boundary conditions, as observed in [16]. We have proved an analogous version in arbitrary dimensions, but more importantly, for almost periodic boundary conditions (see Theorem 5.4). Among other consequences of such a control result is a stronger and more general version of [24, Theorem 2], which has been stated in the particular context of partially rectangular billiards.

Theorem 6.2**.**

For any periodic tube $T$ immersed in a convex $\Omega\subset\mathbb{R}^{n}$, and a neighbourhood $U$ of $\partial T$ inside $\Omega$, there exists a constant $C>0$ such that

[TABLE]

that is, no eigenfunction can concentrate in $T$ and away from $\partial T$.

Here by immersing the periodic tube $T$ into $\Omega$, we assume that the parts of the boundary $\partial\Omega$ where $T$ hits reflects are flat, i.e. consist of a hypersurface. The proof of Theorem 6.2 is absolutely similar in spirit to the proof in [24], except that in higher dimensions, we need to use our version of the control result, given by Theorem 5.4. We also note that the $L^{2}(T)$ norm is equivalent to the norm of the pullback over the mapping cylinder corresponding to $T$. We will skip the details.

7. Acknowledgements

The work was started when the third author was at the Technion, Israel supported by the Israeli Higher Council and the research of the third author leading to these results is part of a project that has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 637851). The work was continued when the third author was at the Courant Institute of Mathematical Sciences funded by Fanghua Lin and the Max Planck Institute for Mathematics, Bonn, and finished at IIT Bombay. He wishes to deeply thank them all. The first two authors gratefully acknowledge the Max Planck Institute for Mathematics, Bonn for providing ideal working conditions and for financial support. During the last stages of this project the first author was in Orsay, where he received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 725967). The second author would like to thank the Fraunhofer Institute, IAIS, for the perfect working conditions and environment. The authors also thank Steve Zelditch and Luc Hillairet for very helpful conversations. Lastly, the authors are deeply grateful to both the Referees for significant improvements in the presentation. In particular, the idea of Proposition 2.4, was suggested by one of the Referees, whose comments also led us to add Remark 5.6 and a new proof of Theorem 3.9 in Appendix A.

Appendix A Some dynamical statements

In this appendix for reader’s convenience we prove the statements from Section 3.2 and give original references. We emphasise that our proof of Theorem 3.9 contains our own original ideas.

We start by recording the fact that any trajectory in a polygon whose forward closure does not contain any singular point must be periodic. In other words, no non-periodic trajectory can continue indefinitely inside the polygon without coming arbitrarily close to the pockets. The proof appears in [14, Theorem 2 and Corollary 2] for polygons; the statement follows from the more general Theorem 3.9.

Theorem A.1**.**

For an arbitrary polygon $P$, and for any $x\in T\Gamma_{1}$, if the closure of the set $\pi\{f^{i}x:i\geq 0\}$ does not contain any vertex of the polygon, then the orbit of $x$ is periodic, where $\pi:TP\to P$ is the usual projection map.

We proceed with the proof of Theorem 3.7. Here we adapt the proof of [14, Theorem 5], written down in [14] for $n=3$, to higher dimensions. One distinction to the $n=3$ and $n\geq 4$ cases in the theorem is that we cannot simply classify all the possible cross-sections of periodic tubes (cf. Remark 3.8).

Proof of Theorem 3.7.

Let us enumerate the faces of $P$ by $\mathcal{F}_{1},\dotso,\mathcal{F}_{l}$. Given an unfolding $P=P^{0},P^{1},P^{2},\dotso$, we write $R_{i}^{(j)}$ for the (affine) reflection in the $i$-th face in $P^{j}$ and denote the corresponding (linear) reflection through the origin by $R^{(j)}_{i,0}$. By Proposition 3.5, we have that $X(w)$ is convex and consists of parallel vectors $v\in\mathbb{R}^{n}$ with basepoints on one particular face; by Lemma 3.3, we have $f^{k}:X(w)\to X(w)$. Here $f$ is the first return map to a face. We assume first $k$ is even. Let $w=(i_{0},\dotso,i_{k-1},i_{0},\dotso)$ where $(i_{0},\dotso,i_{k-1})$ repeats itself and define the rotation

[TABLE]

Since $f^{k}:X(w)\to X(w)$, we must have $\mathcal{R}_{0}v=v$. We will also write

[TABLE]

Here $\tau\in\mathbb{R}^{n}$ is a translation vector.

We may extend $f^{k}$ by continuity to $f^{k}:\overline{X(w)}\to\overline{X(w)}$. By Brouwer’s fixed point theorem, we must have at least one fixed point $x(w)$, i.e. $f^{k}(x(w))=x(w)$. We may write $x(w)=(x_{0},v)$ where $x_{0}\in\mathcal{F}_{i_{0}}$. In fact, by unfolding the trajectories we see, for some $L>0$

[TABLE]

We may choose $x_{0}$ as the origin, to obtain $\tau=Lv$. The argument above also shows $f^{k}(X(w))=X(w)$, since the image of $f^{k}$ is obtained by a rigid motion.

Let the projection to the plane through origin normal to $v$ be denoted by $\pi$. Let $\mathcal{O}\subset\mathcal{F}_{i_{0}}$ denote the basepoints of $X(w)$. The cross-section $\Omega$ is obtained by projecting $\mathcal{O}$ to the direction normal to $v$. Since in an unfolding, the parallel trajectories hit $\mathcal{R}_{0}(\mathcal{O})$, we must have

[TABLE]

The last equality holds since $\pi$ and $\mathcal{R}_{0}$ commute. This shows $\Omega$ is invariant under the rotation $\mathcal{R}_{0}$.999By a slight abuse of notation, here we identify $\mathcal{R}_{0}$ with the rotation in $\pi(\mathbb{R}^{n})$. This is possible as $\mathcal{R}_{0}$ fixes $v$. If we parametrise the tube $T$ in an unfolding by $\Omega\times\mathbb{R}$, we also notice that a point $(x,t)\in\Omega\times\mathbb{R}$ corresponds to the same point in $P$ as does the point $(\mathcal{R}_{0}x,t+L)$. If we think of doubles, this means that $T$ is given by the local isometry $F:\Omega\times\mathbb{R}\to D_{0}$ such that $F(\mathcal{R}_{0}x,t+L)=F(x,t)$ for any $(x,t)$.

In the case when $k$ is odd, we may repeat the previous argument for $\mathcal{R}_{0}\in O(n-1)$ that is orientation reversing. If $n=2,3$, this implies $\mathcal{R}_{0}$ is a reflection in a plane containing $v$ and note that upon applying $\mathcal{R}_{0}^{2}=\operatorname{id}$, we get back to the original position, so every individual trajectory in the tube is periodic.

In the case where $P$ is rational, we notice that the rotation $\mathcal{R}_{0}$ in the above argument is of finite order. Thus we may assume that each trajectory is periodic at the cost of multiplying $k$ by an integer $q$. We have $\Omega$ polyhedral, since it is an intersection of a finite number of faces projected to the plane perpendicular to $v$. ∎

Remark A.2**.**

Observe that the above proof will not work when $P$ is not convex, as one can easily come up with examples (for example, consider a cube with a “needle” poking inward from the ceiling) of non-convex polyhedra where the set $X(w)$ is not simply-connected, which disallows the application of the Brouwer fixed point theorem.

Now we give a proof of Theorem 3.9, proved in [14]. Similarly to Theorem 3.7 above, the proof was written down for $n=3$ in [14] and we give a slightly different argument that avoids the use of “the $N_{\delta}$ sets” and is fully self-contained. The main novelty is in Lemma A.3, that asserts that under a bi-Lipschitz type condition, a map on a compact set is uniformly recurrent, or in other words $C^{0}$-admissible in the sense of Proposition 2.4.

Proof of Theorem 3.9.

Assume the contrary, i.e. assume the trajectory $\gamma$ generated by $x$ is a positive distance $\varepsilon$ away from $\mathcal{S}$. Consider the maximal tube $T=T(x)$ corresponding to $X(w)$.

Consider now an arbitrary limit point $f^{m_{i}}x\to x^{*}\in T\Gamma_{1}$ with $m_{i}\to\infty$. Since the trajectory of $x$ is at $\varepsilon$ distance to $\mathcal{S}$, so is the trajectory of $x^{*}$, and so we may consider the maximal tube $T(x^{*})$ associated to $x^{*}$. Since $w$ is not periodic, eventually $f^{m_{i}}x$ is not parallel to $x^{*}$, as otherwise $f^{m_{i}}x$ would eventually generate the same tube as $T(x^{*})=T(x)$ and thus $w$ is periodic, contradiction. We will show this to be absurd, thus the closure of the trajectory of $x$ must hit $\mathcal{S}$. Write $w^{*}$ for the word associated to $x^{*}$.

Fix a $\delta>0$ and define the (closed) set $X_{\delta}(w^{*})$ to consist of points at a distance at least $\delta$ from $\partial X(w^{*})$ inside $X(\omega^{*})$. We introduce the forward orbit of $X_{\delta}(w^{*})$ as follows

[TABLE]

We denote the closure of this forward orbit by $Y_{\delta}:=\overline{Z_{\delta}}$. We claim that $Y_{\delta}\subset T\Gamma_{1}$ is invariant by $f$.

To prove the claim, consider a sequence $f^{n_{i}}x_{i}\to y\in Y_{\delta}$. Note that $y$ is not tangent to a face, as otherwise this would contradict the assumption that trajectories starting in $\operatorname{int}X(w^{*})$ are a positive distance away from $\mathcal{S}$. Since $f^{n_{i}}x_{i}$ are at a distance at least $\delta$ from $\mathcal{S}$, we have $f(y)$ satisfies the same property and the continuity of $f$ on $T\Gamma_{1}$ proves the claim as $f^{n_{i+1}}x_{i}\to f(y)$.

Next, we claim that there are $c,C>0$, such that for any $m\geq 0$ we have for all $y_{1},y_{2}\in Y_{\delta}$

[TABLE]

Here $d(\cdot,\cdot)$ is the natural distance on $T\Gamma_{1}$. The existence of such constants (independent of $m$), for $y_{1},y_{2}\in Z_{\delta}$ follows since the angle between the unfolded tube $T^{\infty}$ and the faces it hits is bounded from below, as otherwise a trajectory starting in $\operatorname{int}X(w^{*})$ would come arbitrarily close to $\mathcal{S}$, contradicting the definition of a tube. The full claim follows by density and $f$ satisfies this “strong Lipschitz property”.

Now we prove a statement similar to Proposition 2.4.101010Compared to [12, Theorem 1.16] and to the proof found in [14], the following lemma proves that all points in $Y_{\delta}$ are uniformly recurrent at the same time. The key to this is that $f$ satisfies the strong property (A.1) and the idea from the proof of Proposition 2.4.

Lemma A.3**.**

For each $\eta>0$, the set of $k\in\mathbb{N}_{0}$ such that $\operatorname{dist}_{C^{0}}(f^{k}|_{Y_{\delta}},\operatorname{id})<\eta$ is relatively dense.

Proof.

The proof is similar to the proof of Proposition 2.4 and we just sketch the main differences. The main idea is that instead of having $f$ an isometry, we will use the property (A.1).

Take $y\in Y_{\delta}$ and consider the orbit $\mathcal{O}_{y}=\{f^{i}y:i\geq 0\}$. By compactness, there are $n_{1},\dotso,n_{N}\in\mathbb{N}_{0}$ such that $f^{n_{1}}y,\dotso,f^{n_{N}}y$ are $\eta$-dense in $\overline{\mathcal{O}_{y}}$. Then the set of $m\in\mathbb{N}_{0}$ such that $f^{m}y$ is $\eta$-close to $y$ is relatively dense, as for each $m$ there is an $i$ such that, using (A.1)

[TABLE]

If we set $U=B(y,\eta)$, by above we see $f^{m-n_{i}}y\in f^{m-n_{i}}U\cap B(y,\frac{\eta}{c})$.

Consider a cover of $Y_{\delta}$ by balls $U_{i}=B(x_{i},\eta)$ for $i=1,\dotso,N$ and apply the previous claim to the map $(f,\dotso,f):(Y_{\delta})^{N}\to(Y_{\delta})^{N}$ (it satisfies the property (A.1) with the same constants) and the point $(x_{1},\dotso,x_{N})$. Thus there is a relatively dense set $\mathcal{T}\subset\mathbb{N}_{0}$ such that $B(x_{i},\frac{\eta}{c})\cap f^{k}U_{i}\neq\emptyset$ for every $k\in\mathcal{T}$ and $i=1,\dotso,N$.

Take $y\in Y_{\delta}$, $k\in\mathcal{T}$ and choose $y_{i}\in B(x_{i},\frac{\eta}{c})\cap f^{k}U_{i}$. There is an $i$ such that $y\in U_{i}$ and so

[TABLE]

Since $f^{k}y\in f^{k}U_{i}$ and using (A.1), we have $d(f^{k}y,y_{i})<2C\eta$. By triangle inequality

[TABLE]

This finishes the proof. ∎

By the proof of Proposition 3.5, $X(w^{*})$ is convex and every trajectory generated by a point on $\partial X(w^{*})$ comes arbitrarily close to the singular set. Take $\delta=\varepsilon/4$ and pick $\varepsilon/4$-dense points $x_{1},\dotso,x_{N}\in\partial X_{\delta}(w^{*})$ such their trajectories come $\varepsilon/2$-close to the singular set in time (length) $L_{0}$. Denote by $T_{\delta}^{\infty}(f^{k}x^{*})$ the unfolded tube generated by $X_{\delta}(w(f^{k}x^{*}))$, i.e. the complement of the $\delta$-neighbourhood of $\partial T^{\infty}(f^{k}x^{*})$. Take $\eta$ in Lemma A.3 small enough, such that $\operatorname{dist}_{C^{0}}(f^{k}|_{Y_{\delta}},\operatorname{id})<\eta$ for $k\in\mathcal{T}$ relatively dense, so that the following holds. For $k\in\mathcal{T}$ and each $i=1,\dotso,N$, the points $f^{k}x_{i}$ are close enough to $x_{i}$, so that the trajectories generated by $f^{k}x_{i}$ come $\varepsilon/2$-close to the singular set in time $L_{0}$. Since $\mathcal{T}$ is relatively dense, this means that there is an $L>0$ such that for each $i$, the trajectories generated by $x_{i}$ come $\varepsilon/2$-close to the singular set in every segment of length $L$.

Finally, consider $T^{\infty}(f^{m_{i}}x)$, i.e. the maximal unfolded tube generated by $f^{m_{i}}x$; it has a minimal diameter $\varepsilon$ by the assumptions. Taking $i\to\infty$ and by the construction above, a segment of length $L$ of one the trajectories generated by $x_{j}$, together with its $\varepsilon/2$-neighbourhood, will eventually be contained in the interior of $T^{\infty}(f^{m_{i}}x)$. Thus the interior of $T^{\infty}(f^{m_{i}}x)$ will contain a singular point, contradicting the definition of maximal tubes. ∎

Appendix B On a control estimate

In this Appendix we prove a result on control and observability for an elliptic PDE on a product space, which has been generalised in Theorem 5.4 above. The result is proved in unpublished lecture notes of N. Burq (see the footnote of [8, pp. 17]).

To this end, let $(M_{x},g_{x})$ and $(M_{y},g_{y})$ be two compact Riemannian manifolds with Lipschitz boundary. We will consider the product manifold $(M_{x}\times M_{y},g_{x}\oplus g_{y})=:(M,g)$. Denote by $-\Delta_{g}$ the positive-definite Laplace-Beltrami operator.

Recall that we say a subset $A\subset M_{x}$ satisfies the geometric control condition or just (GCC) if every geodesic $\gamma$ in $M_{x}$ hits $A$ in finite time.

Theorem B.1**.**

Let $s\in\mathbb{R}$ and assume $u\in H_{0}^{1}(M)$ satisfies

[TABLE]

Now let $\omega\subset M_{x}$ be an open, non-empty set, satisfying (GCC). Then there exists a constant $C=C(M,g)>0$ such that the following observability estimate holds:

[TABLE]

Proof.

Consider the basis of orthonormal Dirichlet eigenfunctions $\{e_{k}\}_{k=1}^{\infty}\subset L^{2}(M_{y})$, such that $-\Delta_{g_{y}}e_{k}=\lambda_{k}e_{k}$, where $0<\lambda_{1}\leq\lambda_{2}\leq\dotso\nearrow\infty$ are eigenvalues. We consider $u$ and $f$ satisfying (B.1), and consider the Fourier expansion into $e_{k}$:

[TABLE]

The first sum converges in $H^{1}(M)$ as $u\in H_{0}^{1}(M)$ and thus the second sum converges in $H^{-1}(M)$. Moreover, we have that $u_{k}\in H_{0}^{1}(M_{x})$ and $f_{k}\in H^{-1}(M_{x})$ as $u\in H_{0}^{1}(M)$, and we may show the following formulas for all $k\in\mathbb{N}$ hold

[TABLE]

Here $d\operatorname{vol}_{y}$ denotes the volume density on $M_{y}$. Then we have, by using (B.1) and the eigenfunction condition:

[TABLE]

Therefore, we conclude for $k\in\mathbb{N}$

[TABLE]

The remaining ingredient is the observability estimate on $(M_{x},g_{x})$ given in Lemma 5.1. The claim follows by applying the lemma to equation (B.3) for each $k\geq 1$

[TABLE]

∎

We can actually handle periodic conditions similarly, so we have

Theorem B.2**.**

Assume the same setup as in the previous theorem, just with $M_{y}=[0,1]$. This time we take periodic boundary conditions in (B.1), i.e.

[TABLE]

Then the estimate (5.9) holds.

Proof.

The proof is completely analogous to the proof of the previous theorem, by just noting that we may look at $u$ as a function on $M_{x}\times S^{1}$. Then we expand $u$ into Fourier basis on $S^{1}$ and apply Lemma 5.1. ∎