1 Introduction
A billiard trajectory in a Euclidean triangle moves at unit speed and reflects specularly at an open side: its tangential velocity is preserved and its normal velocity changes sign. The natural probability measure is normalized area in the table times uniform angular measure. Ergodicity asks whether every measurable set invariant under this flow has probability zero or one. Straight sides supply no dispersing curvature, while trajectories meeting a vertex have no prescribed continuation. The effect of repeated reflections therefore carries the essential dynamical difficulty.
Straightening reflected trajectories by unfolding connects polygonal billiards to flat geometry; the phase-space construction and null vertex-hit exception appear explicitly in Zemlyakov–Katok [12]. For a rational polygon, whose angles are rational multiples of \pi, the finite reflection group separates the full phase space into invariant directional surfaces. Kerckhoff, Masur, and Smillie proved unique ergodicity in almost every direction [6]; see also the account of Masur–Tabachnikov [7, Sections 3.2 and 3.5]. The same work yields full-flow ergodicity for a dense G_\delta family of tables. Its rational-approximation and Baire-category argument does not assert ergodicity for each irrational triangle.
Vorobets gave a quantitative sufficient approximation condition and explicit ergodic examples [8, Definition 1.3 and Theorem 1.1]. The condition requires exceptionally fast simultaneous rational approximation of the angles; his examples include irrational right triangles. More recently, Chaika and Forni proved weak mixing of the full billiard flow for a dense G_\delta set of n-gons for each fixed n\ge3 [2]. These results distinguish generic or specially approximable tables from a conclusion for every triangle having an irrational angle. Numerical studies of irrational right triangles and fully irrational isosceles triangles have raised doubts about ergodicity [9, 10]. Their finite-time evidence does not establish a positive-measure invariant set and leaves room for very slow convergence of averages.
Theorem 1. Let Q be the interior of a nondegenerate Euclidean triangle with angles \alpha,\beta,\gamma. Suppose that at least one of
is irrational. The unit-speed billiard flow \Phi_t, with specular reflection at the open sides, is ergodic on Q\times S^1 for the probability measure
Here trajectories hitting a vertex in either time direction are discarded. In particular, if a measurable set A satisfies \mu(\Phi_t^{-1}A\mathbin\triangle A)=0 for every t\in\R, then \mu(A)\in\{0,1\}.
In particular, the theorem applies to every triangle for which all three angle ratios are irrational. It also includes triangles with one rational angle and two irrational angles. No genericity or Diophantine condition is imposed. The conclusion concerns the full unit-speed phase space; it is an ergodicity statement for its specified probability measure, not unique ergodicity, mixing, or a claim about each individual orbit.
For each fixed triangle in Theorem 1, the classical conclusion also gives interior and boundary quantum ergodicity for every Dirichlet or Neumann eigenbasis. For every real such eigenbasis, the number of nodal domains tends to infinity along a density-one subsequence. We state the precise spectral conclusions and verify their hypotheses in Section 7.
The main issue is low regularity. Forni and Moll provide an analytic framework for flat surfaces with cone points [3]. Their angular Fourier decomposition and Cauchy–Riemann identities relate transport invariance to a recurrence among coefficient derivatives. Under non-rational holonomy, their Theorem 4.1 proves constancy of an L^2 function invariant under both horizontal directions. Their Corollary 4.2 proves constancy of a geodesic-invariant function in the horizontal Sobolev space H^{1,0}, defined by completion of a cone-adapted smooth domain with square-integrable derivatives in both horizontal directions [3, Section 3.4]. A bounded invariant indicator has no such regularity a priori. We work directly with bounded invariant functions, including invariant indicators. The localized energy estimate in Section 4 first gives an L^2 gradient for each angular coefficient. This is a bound on individual coefficients, not a square-summable family of gradients. The remaining step is to prove that each such gradient vanishes.
Proof strategy. Glue two copies of the triangle along corresponding sides to form its oriented flat double, whose singular points lie at the three vertices. On its unit tangent bundle, let X denote differentiation along a geodesic and Y differentiation in the perpendicular parallel direction. The equation Xf=0 relates neighboring angular Fourier coefficients of an invariant function f. For each integer m, the sum of angular modes m,m+2,m+4,\ldots has an X-derivative consisting of just one endpoint coefficient derivative. We must show that this endpoint is zero.
We smooth only in position, away from the three cone points. Averaging over long segments that stay clear of the cones gives bounded test functions close to f, with a small transport derivative. Integrating the smoothed endpoint equation along these segments yields a pairing that tends to zero. The energy estimate controls the two coefficient gradients in that pairing while the smoothing scale tends to zero. The resulting weak limit satisfies the exact local transport equation. Applying the bounded-coefficient estimate again supplies a new bound independent of the remaining approximation parameter. Only then do we let the weak limits approach f. This second passage proves that the endpoint derivative, and hence every coefficient gradient, vanishes. No horizontal Sobolev hypothesis is imposed on f.
Parallel coefficients must respect the rotation acquired around a cone point. A single irrational angle eliminates all nonzero angular modes; connectedness makes the remaining mode constant. This isolates the angle hypothesis from the analytic argument.
All the analytic steps needed for this passage are proved below. Constants may depend on the fixed triangle, but not on Fourier indices or small smoothing parameters.
Section 2 sets up the double and coefficient bundles. Section 3 supplies the supported smoothing identities, and Section 4 proves the individual-coefficient estimates. Section 5 carries out the two successive limits. Section 6 completes the holonomy argument, and Section 7 derives the spectral consequences.
2 The double and angular coefficients
We pass to a surface on which reflection becomes straight continuation. Throughout the analytic argument, the triangle is fixed and nondegenerate; no irrationality assumption is needed until Section 6.
Glue two closed copies of the triangle along corresponding sides, give the copies opposite orientations, and remove the three vertices. The result is a connected oriented flat surface M. Its completion is compact and has cone angles 2\alpha,2\beta,2\gamma at the finite set V of removed points. Open edges are regular: continuation across an edge is described by reflecting a Euclidean chart on one copy. Transitions between oriented flat charts are rotations and translations.
Let SM be its unit tangent bundle, \nu its normalized Liouville measure, and \Psi_t its geodesic flow. Folding to the original triangle sends this flow to the billiard flow. Away from edges, there are two equally weighted sheets, so the folding map sends \nu to \mu. Each regular billiard segment lifts from either initial sheet, and its time map permutes the two sheet choices. The set of trajectories hitting a vertex has measure zero. Indeed, for each fixed finite sequence of crossed sides, unfolding makes a vertex hit the condition that a straight trajectory pass through a fixed point. This is a null condition in position–direction space. The countable union over side sequences, vertices, and both time directions remains null.
There is no additional obstruction from infinitely many reflections in finite time. Such an accumulation could occur only at a vertex: an open side is regular upstairs, and away from the vertices compactness gives regular neighborhoods. Near a vertex, unfold successive reflections in the wedge of its fixed positive angle. A straight line avoiding the tip sweeps at most \pi in polar angle, so it crosses only finitely many successive wedge sectors. A line ending at the tip likewise crosses only finitely many before reaching it. Thus the excluded set accounts for all finite-time failures of the flow. On its complement the regular geodesic flow is defined for every time. In flat charts it is translation at fixed angle and preserves area times angular measure. Its local preservation, followed along regular trajectories, proves preservation of \nu; folding gives the same assertion for \mu. Edge crossings introduce no singularity upstairs.
In a flat chart write v=(\cos\theta,\sin\theta), and set
These are the forward and perpendicular horizontal vector fields on SM; they commute locally. All distributional equations below are on the regular bundle. In particular, a bounded function invariant almost everywhere under each \Psi_t satisfies Xf=0, by testing against compactly supported smooth functions. We use these equalities in L^2 separately for each time, and under time integrals by Fubini. We do not require a representative invariant pointwise outside one common exceptional set for all real times.
We use the angular Fourier framework for flat geodesic flows developed by Forni and Moll [3, Sections 3.2–3.5], deriving the coordinate rules needed here directly. Fiberwise Fourier projection is globally defined. Locally write
For a change of chart z'=e^{i\phi}z+c, one has \theta'=\theta+\phi and u'_j=e^{-ij\phi}u_j. Thus the coefficients are sections of flat unitary line bundles. Their local derivatives transform by
The flat gradient of a coefficient means its weak gradient in these trivializations, tested against smooth sections compactly supported in the regular surface. An L^2 gradient will always mean a finite global norm on that surface. It does not specify a trace at a cone tip or membership in a completed operator domain. We use the normalized area measure for coefficient norms and \nu for bundle norms; Parseval then has no extra factor. Flat gradients are denoted by \nabla. For fields having horizontal derivatives,
Inner products are complex and linear in the first entry. The mode identities, also valid distributionally, are
Both terms in each sum have the same unitary transition factor.
3 Smoothing away from vertices
The next construction allows differentiation and integration by parts without assigning values at a cone point. Its support margin will be essential when a smoothed field is extended by zero.
Put r(x)=d(x,V) in the completed double. There is K>0 such that, for all sufficiently small \eps>0, every center with r(x)>K\eps has an embedded flat disk of radius 4\eps. Near a fixed cone point a disk of radius c r(x) is embedded for some c>0; away from the ends, compactness supplies the assertion. Also, for each fixed L,
Choose a nonnegative smooth radial kernel of mass one in the unit disk, denoted by \rho, and put \rho_\eps(y)=\eps^{-2}\rho(y/\eps). At a center with r>K\eps, define
in an embedded flat disk, keeping the direction parallel. The formula uses ordinary Euclidean area inside that chart. Radiality makes it independent of the oriented flat chart. Call \{r>K\eps\} the safe region. When a norm is taken over a base set, its inverse image in SM is understood.
Lemma 2. The operators just defined have the following properties.
- 1.
On the safe region, S_\eps and \eps\nabla S_\eps are uniformly bounded from L^2(SM) to L^2. They commute locally with X,Y, including on distributions, and \norm{S_\eps u}_\infty\le\norm{u}_\infty.
- 2.
For every fixed u\in L^2(SM),
\eps\norm{Y S_\eps u}_{L^2(r>K\eps)}\longrightarrow0. \tag{4}On compact regular sets, S_\eps u\to u locally in L^2.
- 3.
Fix L>K+4. If g\in L^2(SM) vanishes on \{r<L\eps\}, then S_\eps g is supported in \{r\ge(L-1)\eps\} and extends by zero to a field with compact spatial support in M. For u\in L^2(SM),
\ip{Y S_\eps u}{g}=-\ip{u}{Y S_\eps g}. \tag{5}The left pairing uses only safe centers. The analogous identity without Y is symmetric. If Xg\in L^2 distributionally, then
X(S_\eps g)=S_\eps(Xg),\qquad \norm{X S_\eps g}_2\le C\norm{Xg}_2, \tag{6}with the same zero-extension convention.
Proof. In a flat chart containing the convolution disks, the operation is ordinary spatial convolution at fixed angle. Its kernel, and its first spatial derivative multiplied by \eps, have uniformly bounded absolute integrals in either spatial variable. Parallel transport preserves angular measure. The usual kernel bound on L^2 therefore proves the asserted norms. Local commutation, positivity, and the L^\infty estimate follow from the same chart description.
For a smooth compactly supported input, the expression in (4) tends to zero directly: for small scales its convolution remains away from the tips and its first derivatives are uniformly bounded in L^2. Such inputs are dense in L^2(SM); the uniform scaled-derivative bound proves the claim for any fixed input. The ordinary approximate-identity argument gives local convergence.
For the third assertion, a contributing kernel pair has base points at distance at most \eps. If one point is in the support of g, both lie strictly inside the safe region. Reversal of the short segment reverses its parallel transport and leaves the radial kernel unchanged. This gives symmetry. In a common flat chart the kernel satisfies Y_x K_\eps(x-y)=-Y_y K_\eps(x-y), with the angles identified by parallel transport. Kernel differentiation and Fubini give (5). No derivative of g is required.
The support estimate follows from the triangle inequality for r. In particular, the convolution vanishes in a neighborhood of the boundary of the safe region, so its extension by zero creates no boundary derivative. Distributional local commutation proves (6); Xg is supported in the same closed spatial support as g, and the L^2 kernel bound gives the final inequality.
4 A localized coefficient estimate
The estimate below uses the full gradient only where transport invariance fails. This will confine the cost of cutting off an invariant function to small cone annuli. The square summability used in its proof belongs to the regularized test field; it will not be assumed for the original invariant function. The equal Cauchy–Riemann norms and parity recurrence also underlie Forni–Moll’s Sobolev arguments [3, Lemma 3.3, Corollary 4.2, and Lemma 5.1]. We prove the support-localized version needed here on compact subsets of the regular surface.
Lemma 3. Let w\in L^2(SM) have compact spatial support in M, with Xw,Yw\in L^2(SM) in the sense of distributions. If Xw vanishes over M\setminus E, where E\subset M is measurable, then for every j\in\Z,
Proof. For each coefficient, integration by parts gives
Indeed, the difference of the first two squared norms is proportional to the integral of dw_j\wedge d\overline{w_j}. The primitive w_j\,d\overline{w_j} is global because its constant unitary transition factors cancel, and its support is compact in the regular surface. Stokes gives zero. Spatially smooth fields suffice initially; ordinary local Sobolev approximation gives the stated generality.
Put D_j=\norm{\pa w_j}_2^2. By (2) and (8),
The pairing is over E, since its first factor vanishes elsewhere. Moreover D_k\to0 as k\to+\infty, because the horizontal derivatives of w are square summable in their Fourier modes. Telescoping in steps of two and applying Cauchy–Schwarz gives
Pointwise in the base variable, the sum inside is bounded by
Lemma 4. On the regular double M of a fixed nondegenerate triangle, suppose h\in L^\infty(SM), \norm{h}_\infty\le H, and Xh=0 locally in distributions. Then every coefficient has a global square-integrable flat gradient, and
Proof. Take a smooth spatial cutoff \chi_\eps equal to zero for r\le(K+1)\eps, equal to one for r\ge(K+2)\eps, with |\nabla\chi_\eps|\le C/\eps. For small \eps, its transition regions are disjoint cone annuli. Set w=\chi_\eps S_\eps h, extended by zero near the vertices. This field has compact spatial support and square-integrable horizontal derivatives. Local commutation gives
Thus Xw is supported in the transition annuli E_\eps, whose total area is O(\eps^2). On these annuli the kernel derivative estimate and boundedness of h give
Both norms on the right of (7) are at most CH. Hence \norm{\nabla w_j}_2\le CH, independently of \eps,j.
For each fixed j, local L^2 convergence of w_j to h_j and weak compactness of its globally bounded gradients identify the weak gradient limit with \nabla h_j, by testing in regular charts. Lower semicontinuity gives (9). This argument neither assigns traces at the cone tips nor asserts square summability over j.
The hypothesis of Lemma 4 is the local distributional equation, rather than invariance under the globally defined measurable flow. This distinction will let us apply it to a weak limit in the next section without first establishing a new global invariance statement.
5 From flow invariance to parallel coefficients
We now show that the individual gradients just obtained are zero. The proof uses the regularized coefficient estimate twice: first to pass a fixed-mode pairing through the smoothing limit, and then to pass it through a second approximation limit. The two limits have different uniform estimates.
Proposition 5. If f\in L^\infty(SM) is invariant almost everywhere under \Psi_t for each t\in\R, then
as distributions on the regular surface.
Proof. Write H=\norm{f}_\infty. The case H=0 is immediate. We know Xf=0, so Lemma 4 applies. Fix m\in\Z throughout the following argument.
The endpoint equation. Define
The latter expression is a global field of mode m-1. The recurrence (2) implies
distributionally. One precise verification uses finite tails U_N:
The terminal fields R_N are uniformly bounded in L^2 by Lemma 4, and their angular modes tend to infinity. For any fixed \varphi\in L^2(SM),
Thus they converge weakly to zero; their norms need not tend to zero. Since U_N\to U in L^2, (11) follows. On safe disks, smoothing gives
Our aim is to prove \ip{B}{P_{m-1}Yf}=0. The recurrence (2) gives P_{m-1}Yf=-2iB, so that identity forces B=0; the construction below obtains it without assuming Yf\in L^2.
Long segments with clearance. Fix s>0, put T=s/\eps, and fix L>K+4. This time scale keeps the discarded mass of order s, while the endpoint error will vanish as \eps\to0 at fixed s. Let G_\eps be the set of regular states whose trajectories satisfy r(\Psi_t z)>L\eps throughout 0\le t\le T. The distance function is 1-Lipschitz along a unit-speed path. A time mesh of spacing at most \eps, with at most T/\eps+2 points, detects every failure of this condition inside \{r\le(L+1)\eps\}. Measure preservation and (3) therefore give
Null sets of singular states may be ignored.
Set a_\eps=f\ind_{G_\eps} and form the negative-time average
It has absolute value at most H. If an integrand is nonzero at z, then z lies at time t on a good segment starting at \Psi_{-t}z. Consequently g_\eps=0 on \{r<L\eps\}. Invariance of f, Minkowski’s inequality, and (13) show that
For the second assertion, finite-time averaging gives the weak identity
It follows by testing (14) against a compactly supported smooth regular field, changing variables by the flow, and integrating its derivative in t. It does not differentiate the indicator of G_\eps.
Integrating the smoothed equation. We justify use of the fundamental theorem in (12) along good segments, despite the absence of angular smoothing. For each fixed \eps, choose a countable cover by flat charts enlarged to contain the required convolution disks. For almost every angular slice in each chart, the smoothed fields have smooth spatial representatives and their differential identity holds everywhere in position. Overlap identities have the same property; include their exceptional angular sets in this countable collection.
During a visit to a flat chart a regular orbit has constant chart angle. If it met one of the exceptional slices, it would remain in that slice for an open time interval, which contains a rational time. For each chart and rational time the set of such initial states is null, by measure preservation. Removing their countable union makes the calculus valid along every chart visit of almost every good segment. Smooth spatial representatives agree on overlaps, so their identities can be integrated along the whole segment. This argument also handles its endpoints by regular continuation there.
Changing variables in the pairing with (14), and using the fundamental theorem, now gives
Both endpoints are safe. By Cauchy–Schwarz and measure preservation, the absolute value of the right side is at most
by (4). Here U is fixed.
Put w_\eps=S_\eps g_\eps, with the supported extension of Lemma 2. Its adjoint identity converts (16) into
The same lemma and (15) give
It has compact spatial support in M. Applying Lemma 3 with E=M yields
The first limit: fixed s. Take a sequence \eps\to0 and a subsequence along which w_\eps converges weakly in L^2 to g^{(s)}. The closed convex set \{|u|\le H\}\subset L^2 is weakly closed, so \norm{g^{(s)}}_\infty\le H. For a fixed smooth compactly supported test field \varphi, symmetry of convolution gives, for sufficiently small \eps,
The uniform L^2 bound and density of these tests show that g_\eps has the same weak L^2 limit. Therefore
The first equality follows either from (15) or the bound for Xw_\eps; the second follows by weak lower semicontinuity.
Only coefficient derivatives at modes m and m-2 contribute to (17), because B has mode m-1. For fixed s, (18) allows a further subsequence on which both these gradients converge weakly in their global L^2 spaces. Testing in regular charts identifies the limits as the gradients of the corresponding coefficients of g^{(s)}. Since B is fixed in L^2, passing to the limit gives
Here this projected derivative means the combination of the two coefficient derivatives in (2); no full L^2 bound on Yg^{(s)} is asserted.
The second limit: s\to0. We do not use the deteriorating bound (18) for this limit. Instead, g^{(s)} is bounded by H and solves the exact local equation in (19), so Lemma 4 gives
uniformly in s. Choose s_n\downarrow0 and the limits just constructed for each s_n. By (19), they converge strongly to f in L^2. Their relevant gradients are bounded in L^2, and every weak subsequential limit is identified locally with the corresponding gradient of f. Thus (20) implies
The equation Xf=0 in mode m-1 says \pa f_{m-2}=-\pb f_m. By (2),
Hence \pb f_m=0. Since m was arbitrary, the recurrence also gives \pa f_j=0 for every j, proving the proposition.
6 Holonomy and the conclusion
Proof of Theorem 1. Relabel the vertices so that \alpha/\pi is irrational. Let A be a measurable invariant set as in the theorem, and lift its indicator to a bounded invariant function f on SM. By Proposition 5, each coefficient has zero flat gradient. In every regular chart it therefore has a constant representative, and these constants match under the unitary transitions. Each coefficient is a parallel section.
Continuation around the vertex of angle \alpha has holonomy rotation by -2\alpha modulo 2\pi, as follows by developing its punctured cone of angle 2\alpha. A parallel coefficient of mode j must consequently be fixed by the phase e^{2ij\alpha}. If j\ne0, irrationality of \alpha/\pi makes that phase different from one, so the coefficient vanishes. The zeroth coefficient is constant on the connected surface M. Fourier completeness shows that f is constant almost everywhere. As an indicator with the same mass as A, it has integral zero or one. Thus \mu(A)\in\{0,1\}.
Remark 6. The precise angle condition is visible in the last step. All the analytic arguments apply to the double of a fixed nondegenerate triangle. A single irrational angle is enough to eliminate every nonzero parallel mode, which is why Theorem 1 includes triangles with one rational angle as well as triangles with three irrational angles. Since the three angles sum to \pi, exactly one irrational angle and two rational angles cannot occur. The constants in the analytic estimates are allowed to depend on the triangle.
7 Quantum ergodicity and nodal domains
The full-flow theorem has spectral consequences for each fixed triangle in its scope. The following statement keeps the boundary-condition dependence of the Cauchy data and their limiting measure explicit.
Corollary 7. Let Q satisfy the hypotheses of Theorem 1, and fix either Dirichlet (B=D) or Neumann (B=N) boundary conditions. Let (u_j)_{j\ge1} be any complete orthonormal eigenbasis of L^2(Q) for the corresponding Laplacian, with
where frequencies are listed with multiplicity. A set S\subset\mathbb N has density one if \#(S\cap\{1,\ldots,N\})/N\to1. Then the following conclusions hold.
- 1.
There is a density-one set S_{\mathrm{int}} such that, for every classical order-zero pseudodifferential operator A whose Schwartz kernel is compactly supported in Q\times Q,
\lim_{\substack{j\to\infty\\j\in S_{\mathrm{int}}}} \langle Au_j,u_j\rangle_{L^2(Q)} =\int_{Q\times S^1}\sigma_0(A)\dd\mu.Here the Euclidean metric identifies unit covectors with Q\times S^1, and \sigma_0(A) is restricted to them.
- 2.
Let \Gamma be the regular boundary, the union of the three open sides. For \lambda_j>0, let h_j=\lambda_j^{-1} and define the Cauchy data on the open sides by
u_j^b=\begin{cases} h_j\,\partial_\nu u_j|_\Gamma,&B=D,\\ u_j|_\Gamma,&B=N, \end{cases}where \nu is the outward unit normal. Write B^*\Gamma=\{(s,\eta)\in T^*\Gamma:|\eta|<1\}, using arc length s and signed tangential covector \eta on each side. There is a density-one set S_\partial\subset\{j:\lambda_j>0\} such that, for every a\in C_c^\infty(T^*\Gamma) and a properly supported order-zero semiclassical quantization A_h=\operatorname{Op}_h(a),
\lim_{\substack{j\to\infty\\j\in S_\partial}} \langle A_{h_j}u_j^b,u_j^b\rangle_{L^2(\Gamma)} =\frac{2}{\pi\operatorname{Area}(Q)} \int_{B^*\Gamma}a(s,\eta)w_B(\eta)\dd s\,d\eta,where
w_D(\eta)=\sqrt{1-\eta^2},\qquad w_N(\eta)=(1-\eta^2)^{-1/2}.The Cauchy data are not further normalized to have unit boundary norm.
- 3.
If the eigenbasis is real-valued, there is a density-one set S_{\mathrm{nod}} such that the number of nodal domains of u_j, the connected components of Q\setminus\{u_j=0\}, tends to infinity as \lambda_j\to\infty with j\in S_{\mathrm{nod}}.
Proof. Place \overline Q in a sufficiently large flat torus. Its boundary is Lipschitz, with three smooth open sides and a finite codimension-two singular set; incident sides meet transversely. Thus the triangle meets the piecewise-smooth hypotheses in [1, Section 2], [4, Definition 3.1], and [5, Definition 1.16], including the last reference’s metric-extension condition. Theorem 1 supplies the required billiard-flow ergodicity. The first conclusion is the interior quantum-ergodicity theorem of Zelditch and Zworski [11], in the Dirichlet and Neumann formulation recalled in [1, Theorem 2 and Remark 2.1].
For the boundary conclusion, discard the null sets of vertex hits and tangential collisions. The collision section is B^*\Gamma, with invariant measure ds\,d\eta. Boundary-flight coordinates identify the flow, up to null sets, with the suspension of its collision map under the positive free-flight time \tau, and Liouville measure is proportional to ds\,d\eta\,dt for 0\le t<\tau. The roof \tau is positive almost everywhere and bounded by the diameter of Q. Hence an invariant collision set of positive measure with positive-measure complement would suspend to a nontrivial invariant flow set, contradicting Theorem 1. The collision map is therefore ergodic. Hassell–Zelditch’s boundary theorem [4, Theorem 1.1, equations (1.3)–(1.5), and Section 7] now gives the second conclusion: its prefactor is 4/(\operatorname{vol}(S^1)\operatorname{Area}(Q)) =2/(\pi\operatorname{Area}(Q)). Their inward normal instead of our outward normal changes the Dirichlet data by a sign and leaves the quadratic pairing unchanged.
Finally, Hezari’s nodal-domain theorem for two-dimensional ergodic billiards with piecewise-smooth boundary [5, Theorem 1.6] gives the third conclusion for real eigenfunctions. The finitely many zero-frequency terms have no effect on density. The three density-one sets may be replaced by their finite intersection when a simultaneous subsequence is desired.
The subsequences may depend on the fixed triangle, boundary condition, and eigenbasis. Zero-density exceptional subsequences remain possible: these conclusions assert neither quantum unique ergodicity nor an equidistribution or nodal-count rate, mixing, or a spectral-statistics law.
References
- [1]Nicolas Burq. Quantum ergodicity of boundary values of eigenfunctions: A control theory approach. 2003. arXiv:math/0301349v1, January 30, 2003. arXiv:math/0301349v1, https://arxiv.org/abs/math/0301349v1.
- [2]Jon Chaika and Giovanni Forni. Weakly mixing polygonal billiards. Annals of Mathematics, 203(3):695–735, 2026. doi:10.4007/annals.2026.203.3.1. doi:10.4007/annals.2026.203.3.1, arXiv:2003.00890.
- [3]Giovanni Forni and Nelson Moll. Cohomological equation for geodesic flows on flat surfaces. 2025. arXiv:2510.18128v1, October 20, 2025. arXiv:2510.18128v1, https://arxiv.org/abs/2510.18128v1.
- [4]Andrew Hassell and Steve Zelditch. Quantum ergodicity of boundary values of eigenfunctions. Communications in Mathematical Physics, 248(1):119–168, 2004. doi:10.1007/s00220-004-1070-2. Cited numbering follows arXiv:math/0211140v4, November 24, 2003. doi:10.1007/s00220-004-1070-2, arXiv:math/0211140v4.
- [5]Hamid Hezari. Quantum ergodicity and L^p norms of restrictions of eigenfunctions. 2016. arXiv:1606.08066v1, June 26, 2016. arXiv:1606.08066v1, https://arxiv.org/abs/1606.08066v1.
- [6]Steven Kerckhoff, Howard Masur, and John Smillie. Ergodicity of billiard flows and quadratic differentials. Annals of Mathematics, 124(2):293–311, 1986. doi:10.2307/1971280. doi:10.2307/1971280.
- [7]Howard Masur and Serge Tabachnikov. Rational billiards and flat structures. In Boris Hasselblatt and Anatole Katok, editors, Handbook of Dynamical Systems, pages 1015–1089, North-Holland, 2002. Author version dated December 16, 1999, https://math.uchicago.edu/~masur/handbook1.pdf. doi:10.1016/S1874-575X(02)80015-7. doi:10.1016/S1874-575X(02)80015-7.
- [8]Ya. B. Vorobets. Ergodicity of billiards in polygons. Sbornik: Mathematics, 188(3):389–434, 1997. doi:10.1070/SM1997v188n03ABEH000211. doi:10.1070/SM1997v188n03ABEH000211.
- [9]Jiao Wang, Giulio Casati, and Tomaž Prosen. Nonergodicity and localization of invariant measure for two colliding masses. Physical Review E, 89(4):042918, 2014. doi:10.1103/PhysRevE.89.042918. doi:10.1103/PhysRevE.89.042918, arXiv:1309.7617.
- [10]Katerina Zahradova, Julia Slipantschuk, Oscar F. Bandtlow, and Wolfram Just. Impact of symmetry on ergodic properties of triangular billiards. Physical Review E, 105(1):L012201, 2022. doi:10.1103/PhysRevE.105.L012201. doi:10.1103/PhysRevE.105.L012201.
- [11]Steven Zelditch and Maciej Zworski. Ergodicity of eigenfunctions for ergodic billiards. Communications in Mathematical Physics, 175(3):673–682, 1996. doi:10.1007/BF02099513. doi:10.1007/BF02099513.
- [12]A. N. Zemlyakov and A. B. Katok. Topological transitivity of billiards in polygons. Mathematical Notes, 18(2):760–764, 1975. doi:10.1007/BF01818045. doi:10.1007/BF01818045.