Polyhedra, Moduli,
and Quantum Geometry

Web working edition

On this page
  1. 29.1 The norm-ball obstruction and the correct error scale
  2. 29.2 The Fourier and spectral meaning of the obstruction
  3. 29.3 A coefficient observable from primitive pole-to-pole arcs
  4. 29.4 Uniform rotation norms, including a sharp hyperbolic cutoff
  5. 29.5 Control of the terminal derivative
  6. 29.6 What the bridge retains
  7. 29.7 Affine displacement sums: an exact isometric model
  8. 29.8 Two marks: the rank-one divergence and its removal
  9. 29.9 An exact Gaussian two-mark correlation
  10. 29.10 Sharp boundaries and the rotation norm

29 From simple-curve coefficients to controlled observables

The preceding section supplies an explicit analytic norm for its mixing estimate. We now connect that norm to a selected family of simple curves and their coefficient distributions. The genus-one norm-ball obstruction first fixes the attainable scale of sharp remainders. Primitive pole-to-pole arcs then give geometric coefficient observables with uniform rotation and terminal-derivative bounds. Exact affine Gaussian calculations provide independent tests of the analytic estimates. Throughout, a fixed-surface sharp count, a moduli-averaged covariance, and the conditioning of coefficient extraction remain separate problems.

29.1 The norm-ball obstruction and the correct error scale

The genus-one length norm gives a necessary benchmark for any passage from smooth return estimates to sharp simple-curve counting. Let XX be a fixed complete hyperbolic once-punctured torus, and let BXB_X be the unit ball of the McShane–Rivin length norm on H1(X;ℝ)≃ℝ2H_1(X;\mathbb R)\simeq\mathbb R^2. Put AX(L)=#(LBX∩ℤ2),PX(L)=#(LBX∩ℤprim2).A_X(L)=\#(LB_X\cap\mathbb Z^2),\qquad P_X(L)=\#(LB_X\cap\mathbb Z^2_{\mathrm{prim}}). The first count includes the origin and nonprimitive classes. The second counts oriented primitive classes; the count of unoriented simple closed geodesics is NX(L)=12PX(L)N_X(L)=\tfrac12P_X(L). Elementary planar lattice counting and Möbius inversion give AX(L)=Area⁡(BX)L2+OX(L),(29.1)PX(L)=Area⁡(BX)ζ(2)L2+OX(Llog⁡L).(29.2)\begin{aligned} A_X(L)&=\operatorname{Area}(B_X)L^2+O_X(L),\\ P_X(L)&=\frac{\operatorname{Area}(B_X)}{\zeta(2)}L^2 +O_X(L\log L). \end{aligned} Indeed, PX(L)=∑dμ(d)(AX(L/d)−1)P_X(L)=\sum_d\mu(d)(A_X(L/d)-1), with d≪XLd\ll_X L; the boundary errors sum as L∑d≪XLd−1L\sum_{d\ll_X L}d^{-1}. These are the upper estimates behind [16, Theorem 4.1], with the orientation factor made explicit here.

The next theorem in that paper concerns a dense set of hyperbolic structures. Its mechanism is an infinitely flat boundary point whose supporting normal is rational with respect to the homology lattice. The subsequent Conjecture 4.3 proposes an O(L1/2+ε)O(L^{1/2+\varepsilon}) error for the modular torus; thus the obstruction is not asserted there for every individual surface. The phrase “at least O(L)O(L)” in Theorem 4.2 should be separated from the literal modern assertion limsup⁡|EX(L)|/L>0\limsup|E_X(L)|/L>0. The elementary argument below proves the precise consequence needed here: no fixed improvement of the exponent below one is possible under the stated flat-point hypothesis. We do not use it to claim a literal Ω(L)\Omega(L) bound.

Proposition 29.1 (A primitive-row obstruction). Let B⊂ℝ2B\subset\mathbb R^2 be a strictly convex compact body with the origin in its interior. Suppose that p∈∂Bp\in\partial B has irrational radial direction, a rational supporting normal, and infinite graph-flatness. Then, for EB(L)=#(LB∩ℤprim2)−Area⁡(B)ζ(2)L2,E_B(L)=\#(LB\cap\mathbb Z^2_{\mathrm{prim}}) -\frac{\operatorname{Area}(B)}{\zeta(2)}L^2, and every ε>0\varepsilon>0, there are cε>0c_\varepsilon>0 and arbitrarily large LL such that |EB(L)|≥cεL1−ε.|E_B(L)|\geq c_\varepsilon L^{1-\varepsilon}.(29.3) In particular EB(L)=O(L1−η)E_B(L)=O(L^{1-\eta}) fails for every fixed η>0\eta>0.

Proof. An integral unimodular change of coordinates makes the supporting normal vertical. Write p=(α,h)p=(\alpha,h), with h>0h>0 and α/h\alpha/h irrational. Near pp, the upper boundary has the form y=h−f(x−α),f(0)=0,y=h-f(x-\alpha),\qquad f(0)=0, where f≥0f\geq0 is convex. Infinite graph-flatness means that for every integer N>2N>2 there are CN,ρN>0C_N,\rho_N>0 such that 0≤f(s)≤CN|s|N(|s|≤ρN).0\leq f(s)\leq C_N|s|^N\qquad (|s|\leq\rho_N). Fix NN, let nn run through primes, and set Ln=nh,δn=Ln−1,Ln+=Ln+δn.L_n=\frac nh,\qquad \delta_n=L_n^{-1},\qquad L_n^+=L_n+\delta_n. The line y=ny=n supports LnBL_nB at the single point (Lnα,n)(L_n\alpha,n). Its first coordinate is irrational, so this row contains no lattice point of LnBL_nB.

In Ln+BL_n^+B, the same row lies a normalized depth hδn/Ln+h\delta_n/L_n^+ below the supporting line. Therefore it contains the horizontal interval centered at Ln+αL_n^+\alpha with half-width rn=min⁡{(h2CN)1/N(Ln+)1−1/Nδn1/N,n4}.r_n=\min\left\{ \left(\frac h{2C_N}\right)^{1/N} (L_n^+)^{1-1/N}\delta_n^{1/N},\ \frac n4 \right\}. For sufficiently large nn the local graph description applies, the first quantity realizes the minimum, and rn≍B,NLn1−2/N.r_n\asymp_{B,N}L_n^{1-2/N}. The interval has length at most n/2n/2. Among its integer points (k,n)(k,n), at most one has n∣kn\mid k; every other point is primitive because nn is prime. Thus PB(Ln+)−PB(Ln)≥2rn−O(1).P_B(L_n^+)-P_B(L_n)\geq 2r_n-O(1). The change in the quadratic main term is bounded, since (Ln+)2−Ln2=2+Ln−2.(L_n^+)^2-L_n^2=2+L_n^{-2}. It follows that max⁡{|EB(Ln)|,|EB(Ln+)|}≥cB,NLn1−2/N\max\{|E_B(L_n)|,|E_B(L_n^+)|\} \geq c_{B,N}L_n^{1-2/N} for all sufficiently large primes nn. Choose N>2/εN>2/\varepsilon to obtain the assertion. To rule out a proposed O(L1−η)O(L^{1-\eta}) estimate, apply the assertion with 0<ε<η0<\varepsilon<\eta. ◻

The proof works directly with primitive classes, so no cancellation in a Möbius sum needs to be controlled. The short interval [Ln,Ln+Ln−1][L_n,L_n+L_n^{-1}] acquires many new primitive points along a nearly flat supporting row, while its predicted main-term increment stays bounded. Infinite flatness alone does not force the row length to be a positive fraction of LnL_n; that is why the proof gives an exponent-one barrier, not a literal linear lower bound. The distinction between a fixed-surface estimate, an estimate uniform over moduli, and a smoothed estimate must be retained.

29.2 The Fourier and spectral meaning of the obstruction

The flatness obstruction has a direct oscillatory-integral interpretation. For illustration, suppose a localized smooth boundary arc has graph y=h−f(x),f(x)=axk+O(xk+1),a>0,y=h-f(x),\qquad f(x)=a x^k+O(x^{k+1}),\quad a>0, where k≥2k\geq2 is even. Apply smooth cutoffs supported near this arc and equal to one at the supporting point. Integrating the localized indicator in the normal coordinate first gives, up to rapidly decreasing terms, −e−2πirh2πir∫χ(x)e2πirf(x)⁡dx.-\frac{e^{-2\pi i r h}}{2\pi i r} \int \chi(x)e^{2\pi i r f(x)}\,\mathop{}\!\mathrm dx. The substitution x=r−1/kux=r^{-1/k}u in the stationary-phase calculation gives χ0𝟏B̂(rν)=ca,ke−2πirhr−1−1/k+o(r−1−1/k),ca,k≠0.\widehat{\chi_0\mathbf1_B}(r\nu) =c_{a,k}\,e^{-2\pi i r h}r^{-1-1/k} +o(r^{-1-1/k}),\qquad c_{a,k}\ne0.(29.4) where ν\nu is the chosen normal direction and χ0\chi_0 is the localized cutoff. The familiar nondegenerate exponent is r−3/2r^{-3/2} when k=2k=2; larger contact order weakens it toward r−1r^{-1}. Localization matters: the contributions of different supporting points can interfere in the Fourier transform of the whole ball. No global lower bound follows from discarding those phases.

Under dilation, 𝟏LB̂(m)=L2𝟏B̂(Lm).\widehat{\mathbf1_{LB}}(m)=L^2\widehat{\mathbf1_B}(Lm). Thus a finite-type flat point with a lattice-rational normal contributes on the scale L1−1/kL^{1-1/k} to the corresponding localized Poisson terms. For an infinitely flat point there is no fixed finite type to insert in this formula. Proposition 29.1 supplies a direct counting argument in that case, without assuming smooth boundary regularity at all nearby points or invoking a formal stationary-phase expansion there.

There is also an exact elementary spectral model. On an auxiliary flat torus ℝ2/ℤ2\mathbb R^2/\mathbb Z^2, define the positive self-adjoint Fourier multiplier 𝖯Xe2πi⟨m,x⟩=∥m∥Xe2πi⟨m,x⟩.\mathsf P_X e^{2\pi i\langle m,x\rangle} =\|m\|_X e^{2\pi i\langle m,x\rangle}. Norm equivalence gives a discrete spectrum tending to infinity. Its counting function is AX(L)A_X(L), and restriction to the closed span of primitive Fourier modes has counting function PX(L)P_X(L). Restricting further to the even part identifies mm and −m-m and gives NX(L)N_X(L). This multiplier is defined spectrally; the nonsmooth length norm is not being asserted to be a classical smooth pseudodifferential symbol. At a differentiable point pp of its unit ball, a supporting covector ν\nu satisfies ∇∥⋅∥X(p)=ν⟨ν,p⟩.\nabla\|\cdot\|_X(p)=\frac{\nu}{\langle\nu,p\rangle}. Thus a lattice-rational normal is exactly a rational direction for the translation velocity of the associated Hamiltonian. This explains the periodic-phase interpretation, without assuming that a global smooth trace formula applies to the length norm.

The norm-ball discrepancy is therefore a spectral-projector remainder for this particular length multiplier. Its flat normal directions produce degenerate phases in a regularized Poisson or trace calculation. This is distinct from complementary-series spectrum for the SL2(ℝ)SL_2(\mathbb R) action on the marked-torus moduli space. In particular, it does not contradict the temperedness statement or the smooth correlation estimate in Section 28. The rotation-Sobolev issue depends on what is held fixed. A sharp indicator of a nonradial region in fixed physical vector coordinates generally has rotational jumps, so smoothing it introduces growing derivative norms. By contrast, hyperbolic length is unchanged when the flat charts are rotated while the underlying conformal surface stays fixed. A sharp cutoff in that invariant length can therefore lie in the rotation generator’s domain, and a radial flat cutoff can be rotation invariant as well. Sharpness alone is not an obstruction to the stated mixing bound. For those invariant observables, the remaining boundary and regularity issues enter when taking flow derivatives or passing from an averaged correlation to a fixed-surface counting remainder.

29.3 A coefficient observable from primitive pole-to-pole arcs

The marked-torus construction supplies an explicit family of simple curves for which the coefficient data become observables with controlled norms. This family carries a flat-length cutoff as well as a hyperbolic-length weight. The two length parameters remain distinct.

Work on the finite labeled model Y2Y_2 of Section 28.5. Write E=ℝ2/gℤ2E=\mathbb R^2/g\mathbb Z^2, with its origin one of the four ramification points of the pillowcase map. Fix a nonzero parity vector δ∈(ℤ/2ℤ)2\delta\in(\mathbb Z/2\mathbb Z)^2. For a nonzero primitive integer vector m≡δ(mod⁡2)m\equiv\delta\pmod2, project the straight segment [0,gm/2]⊂ℝ2[0,gm/2]\subset\mathbb R^2 first to EE and then to E/{±1}E/\{\pm1\}. It joins the pole at the origin to the pole corresponding to δ/2\delta/2.

Lemma 29.2 (Primitive directions give simple curves). Outside a Haar-null invariant subset of Y2Y_2, the preceding segment projects to an embedded arc avoiding the two regular marked points. Its lift to the genus-two double cover is a nonseparating simple closed regular flat geodesic αm\alpha_m, with ℓqX(αm)=|gm|.\ell_{q_X}(\alpha_m)=|gm|.(29.5) Opposite vectors give the same curve; distinct primitive directions modulo sign give distinct free homotopy classes.

Proof. If two parameters 0<s,t<10<s,t<1 on the segment identify on the pillowcase, then (s−t)m∈2ℤ2(s-t)m\in2\mathbb Z^2 or (s+t)m∈2ℤ2(s+t)m\in2\mathbb Z^2. Primitivity forces the corresponding scalar to be an even integer. The first case gives s=ts=t, and the second is impossible. Likewise the interior contains no two-torsion point. The segment hits one of the additional regular marks only if one of its lifts p,qp,q lies on the closed rational circle ℝgm/gℤ2\mathbb Rgm/g\mathbb Z^2. For each mm this is a null condition in the point fiber. The union over the countable set of primitive mm is still null and is invariant under the entire SL2(ℝ)SL_2(\mathbb R) action.

The preimage of an embedded arc joining two branch values of a double cover is a circle. Removing it leaves a connected cover, since the four other branch values still have nontrivial monodromy; thus the circle is nonseparating. At its endpoints the sphere flat angle π\pi doubles to 2π2\pi, so the two lifted segments join to a regular flat geodesic. Its length is twice |gm|/2|gm|/2. Different directions give intersecting regular geodesics with distinct tangents at the Weierstrass point over the origin. Freely homotopic regular geodesics in a nonpositively curved flat metric lie in a flat cylinder and cannot intersect transversely. This proves the last assertion. ◻

Write ℓm(qX)=ℓX(αm)\ell_m(q_X)=\ell_X(\alpha_m) for the hyperbolic length of the same topological curve. Proposition 23.3 identifies it as the position of the first atom of σαmX\sigma^X_{\alpha_m}. In particular, if a local cutoff χm\chi_m is one near that atom and vanishes on the remaining support, then for a bounded length weight HH, H(ℓm)=⟨σαmX,H(s)χm(s)⟩⟨σαmX,χm(s)⟩.H(\ell_m)= \frac{\langle\sigma^X_{\alpha_m},H(s)\chi_m(s)\rangle} {\langle\sigma^X_{\alpha_m},\chi_m(s)\rangle}.(29.6) For arbitrary bounded Borel HH, the same statement is interpreted using the already recovered first-atom location. No assertion that first-atom extraction is a bounded linear functional is being made.

Choose an even f∈Cc2(ℝ2)f\in C_c^2(\mathbb R^2) supported in the unit disk, a bounded Borel function h:(0,∞)→ℂh:(0,\infty)\to\mathbb C, and parameters R≥1R\geq1, L>0L>0. Define ℱδ;f,hR,L(qX)=12R2∑m∈ℤ2primitivem≡δ(2)f(gm/R)h(ℓm(qX)/L).\mathcal F_{\delta;f,h}^{R,L}(q_X) =\frac1{2R^2} \sum_{\substack{m\in\mathbb Z^2\ \mathrm{primitive}\\m\equiv\delta\ (2)}} f(gm/R)\,h\bigl(\ell_m(q_X)/L\bigr).(29.7) The factor 1/21/2 identifies opposite directions. The function is well defined on Y2Y_2: changes of lattice basis preserve the parity class and relabel the sum. More precisely, define the geometric formula on the invariant conull admissible set of Lemma 29.2 and set it to zero on its complement; this fixes a measurable representative of the same L2L^2 class. Averaging over the finite labeling choices gives observables on the unlabelled stratum, with the same norm bounds. Such an average is an average of selected families; it need not count their union without multiplicities.

29.4 Uniform rotation norms, including a sharp hyperbolic cutoff

Put J=(0−110),A=(100−1),DBf(x)=(Bx)⋅∇f(x).J=\begin{pmatrix}0&-1\\1&0\end{pmatrix},\qquad A=\begin{pmatrix}1&0\\0&-1\end{pmatrix},\qquad D_Bf(x)=(Bx)\cdot\nabla f(x). The notation DBD_B here is a linear vector field on ℝ2\mathbb R^2.

Lemma 29.3 (Uniform primitive-vector bound). For every unimodular lattice gℤ2g\mathbb Z^2 and every R>0R>0, #{m∈ℤ2:mprimitive,|gm|≤R}≤2+2πR2.\#\{m\in\mathbb Z^2:\ m\ \mathrm{primitive},\ |gm|\leq R\} \leq 2+2\pi R^2.(29.8)

Proof. Two nonparallel lattice vectors u,vu,v have |det⁡(u,v)|≥1|\det(u,v)|\geq1. If both have length at most RR, their unoriented angular separation therefore has sine at least R−2R^{-2}. For R<1R<1 there is at most one primitive direction. For R≥1R\geq1, placing the directions on the circle of length π\pi shows that their number is at most π/arcsin⁡(R−2)≤πR2\pi/\arcsin(R^{-2})\leq\pi R^2. Each direction contributes only the two opposite primitive vectors. The displayed bound covers both cases without rounding conventions. ◻

Rotation of the flat charts fixes the underlying conformal surface and the intrinsic pole-to-pole arcs. Consequently Θℓm=0,\Theta\ell_m=0, and differentiation of (29.7) under rotation only replaces ff by DJfD_Jf. Set C*=1+π,Nf,h=C*∥h∥∞(∥f∥∞2+∥DJf∥∞2)1/2.C_*=1+\pi,\qquad N_{f,h}=C_*\|h\|_\infty \bigl(\|f\|_\infty^2+\|D_Jf\|_\infty^2\bigr)^{1/2}. The primitive-vector bound proves, uniformly in R≥1R\geq1 and L>0L>0, 𝒮1(ℱδ;f,hR,L)≤Nf,h.\mathcal S_1(\mathcal F_{\delta;f,h}^{R,L})\leq N_{f,h}.(29.9) This is global, including the modular cusp. Although the curve labels can jump as a regular mark crosses a selected arc, the exceptional walls are SL2SL_2-invariant. Along almost every group orbit the labels are fixed under the natural marking, and the preceding derivatives are legitimate. Uniform bounds on the difference quotients give the corresponding strong L2L^2 generator statements.

In particular, Theorem 28.1 now applies to two observables of the form (29.7), with |Corr⁡t(ℱ,𝒢)|≤4coth⁡π(1+t)e−tNf,hNf̃,h̃,|\operatorname{Corr}_t(\mathcal F,\mathcal G)| \leq4\coth\pi\,(1+t)e^{-t}\,N_{f,h}N_{\widetilde f,\widetilde h},(29.10) uniformly in their two flat cutoffs and two hyperbolic-length scales. The hyperbolic weight may be the sharp cutoff h(s)=𝟏s≤1h(s)=\mathbf1_{s\leq1}: its discontinuity costs no rotation derivative. If the flat weight is radial as well, the whole observable is KK-invariant. In that case a sharp disk cutoff is also allowed directly, and the single-KK-type matrix-coefficient bound gives the sharper estimate |Corr⁡t(F,G)|≤Ξ(at)∥F∘∥2∥G∘∥2≤4π(1+t)e−t∥F∘∥2∥G∘∥2.|\operatorname{Corr}_t(F,G)| \leq\Xi(a_t)\|F^\circ\|_2\|G^\circ\|_2 \leq\frac4\pi(1+t)e^{-t}\|F^\circ\|_2\|G^\circ\|_2.(29.11) No rotational smoothing is needed for these radial tests.

29.5 Control of the terminal derivative

For the derivative counting converse, require h∈C1h\in C^1 and H0=∥h∥∞<∞,H1=sups>0s|h′(s)|<∞.H_0=\|h\|_\infty<\infty,\qquad H_1=\sup_{s>0}s|h'(s)|<\infty. Wolpert’s quasiconformal length inequality gives |Vℓm|≤2ℓm,|ΘVℓm|≤4ℓm.|V\ell_m|\leq2\ell_m, \qquad |\Theta V\ell_m|\leq4\ell_m.(29.12) For completeness, the first inequality follows by lifting a KK-quasiconformal map to the cyclic annular covers associated with the curve. Their conformal moduli are π/ℓX\pi/\ell_X and π/ℓY\pi/\ell_Y; modulus distortion gives K−1≤ℓY/ℓX≤KK^{-1}\leq\ell_Y/\ell_X\leq K. The map associated with ata_t has dilatation e2|t|e^{2|t|}. For the second inequality use Θℓm=0\Theta\ell_m=0 and [Θ,V]=±2W[\Theta,V]=\pm2W, where WW is the generator of a rotated diagonal subgroup. The same first inequality applies to WW. See also [110, Equation (1.1)] for the quasiconformal inequality and its attribution to Wolpert.

Let F=ℱδ;f,hR,LF=\mathcal F_{\delta;f,h}^{R,L}, and abbreviate f0=∥f∥∞,fK=∥DJf∥∞,fA=∥DAf∥∞,fKA=∥DJDAf∥∞.f_0=\|f\|_\infty,\quad f_K=\|D_Jf\|_\infty,\quad f_A=\|D_Af\|_\infty,\quad f_{KA}=\|D_JD_Af\|_\infty. Direct differentiation and (29.12) give ∥VF∥∞≤C*(fAH0+2f0H1),∥ΘVF∥∞≤C*(fKAH0+2fKH1+4f0H1).(29.13)\begin{aligned} \|VF\|_\infty &\leq C_*\bigl(f_AH_0+2f_0H_1\bigr),\\ \|\Theta VF\|_\infty &\leq C_*\bigl(f_{KA}H_0+2f_KH_1+4f_0H_1\bigr). \end{aligned} These constants are again independent of R≥1R\geq1 and L>0L>0. No second derivative of hh occurs, because rotation fixes ℓm\ell_m. Thus 𝒮1(VF)\mathcal S_1(VF), and hence 𝒮1((4+V)F)\mathcal S_1((4+V)F), have explicit uniform bounds. In particular the terminal identity extends to these observables: ∫1Te4tCorr⁡t(F,(4+V)G)⁡dt=e4TCorr⁡T(F,G)−e4Corr⁡1(F,G).\int_1^T e^{4t}\operatorname{Corr}_t(F,(4+V)G)\,\mathop{}\!\mathrm dt =e^{4T}\operatorname{Corr}_T(F,G) -e^4\operatorname{Corr}_1(F,G).(29.14) Here GG has a smooth length weight as above; invariance of probability measure gives μ(VG)=0\mu(VG)=0. This proves applicability of the integrated counting converse to these endpoints. It does not identify an individual curve sum with the finite affine return sum of the preceding section. Sharp cutoffs require a separate smoothing or boundary argument when the terminal derivative is used.

29.6 What the bridge retains

For h≡1h\equiv1, the observable reduces to a primitive parity Siegel transform on the modular base and has no nonzero marked-point Fourier mode. For a disk cutoff and fixed gg, elementary lattice counting and Möbius inversion give 12#{mprimitive:m≡δ(2),|gm|≤R}=π6ζ(2)R2+Og(Rlog⁡R)=1πR2+Og(Rlog⁡R).\frac12\#\{m\ \mathrm{primitive}:m\equiv\delta\ (2),\ |gm|\leq R\} =\frac{\pi}{6\zeta(2)}R^2+O_g(R\log R) =\frac1\pi R^2+O_g(R\log R). Indeed, only odd common divisors occur in a nonzero parity class, and ∑doddμ(d)/d2=4/(3ζ(2))\sum_{d\ \mathrm{odd}}\mu(d)/d^2=4/(3\zeta(2)). Adding the hyperbolic weight h(ℓm/L)h(\ell_m/L) can retain dependence on the two regular marks and uses the simple-curve first-atom data. For h≡1h\equiv1 the exact mean is μ(ℱδ;f,1R,L)=16ζ(2)∫ℝ2f(x)⁡dx,\mu(\mathcal F_{\delta;f,1}^{R,L}) =\frac1{6\zeta(2)}\int_{\mathbb R^2}f(x)\,\mathop{}\!\mathrm dx, by the primitive Siegel formula and the equal densities of the three nonzero parity classes.

There is a useful fixed-surface sparsity bound. On a fixed compact (X,qX)(X,q_X) put C(X,qX)=supX|qX|1/2⁡dshyp<∞.C(X,q_X)=\sup_X\frac{|q_X|^{1/2}}{\mathop{}\!\mathrm ds_{\mathrm{hyp}}}<\infty. The selected flat geodesic minimizes flat length in its free homotopy class, so comparison with the hyperbolic geodesic gives |gm|=ℓqX(αm)≤C(X,qX)ℓX(αm).|gm|=\ell_{q_X}(\alpha_m) \leq C(X,q_X)\ell_X(\alpha_m). Consequently the whole selected family with hyperbolic length at most LL lies in the flat window R=C(X,qX)LR=C(X,q_X)L and has size OX,qX(L2)O_{X,q_X}(L^2). It is a sparse family relative to the full genus-two L6L^6 simple-curve count. The bound is for a fixed surface and differential; no uniform comparison constant over moduli is asserted.

Equations (29.6)–(29.13) are a norm-controlled bridge for this selected family. They do not identify the full connected-simple-curve count with an affine return sum. Nor do they prove uniform stability of extracting first atoms from perturbed spectral or coefficient data: the shrinking-gap issue of Section 24.4 remains. The output observables have controlled geometric norms, while the quantitative conditioning of the nonlinear coefficient extraction is a different question.

29.7 Affine displacement sums: an exact isometric model

The homogeneous calculation can be tested on explicit geometric sums. These are the relative-period Siegel–Veech transforms discussed in [108, Section 6.1]. On the marked-torus model write a point as (g,v1,v2)(g,v_1,v_2), with the physical marked points understood modulo gℤ2g\mathbb Z^2. For f∈𝒮(ℝ2)f\in\mathscr S(\mathbb R^2) define Sf,j(g,v1,v2)=∑m∈ℤ2f(gm+vj),Ff,j=Sf,j−∫ℝ2f(x)⁡dx.S_{f,j}(g,v_1,v_2)=\sum_{m\in\mathbb Z^2}f(gm+v_j),\qquad F_{f,j}=S_{f,j}-\int_{\mathbb R^2}f(x)\,\mathop{}\!\mathrm dx. The summands test the displacement vectors of straight segments from the torus origin to its jjth marked point. These are flat geometric observables; a hyperbolic simple-curve interpretation requires a separate choice of a lift and of curve weights.

Use the Fourier transform f̂(ξ)=∫f(x)e−2πix⋅ξ⁡dx\widehat f(\xi)=\int f(x)e^{-2\pi i x\cdot\xi}\mathop{}\!\mathrm dx. With vj=gujv_j=gu_j, Poisson summation gives Ff,j(g,u1,u2)=∑n∈ℤ2\{0}f̂(g−Tn)e2πin⋅uj.F_{f,j}(g,u_1,u_2)= \sum_{n\in\mathbb Z^2\setminus\{0\}} \widehat f(g^{-T}n)e^{2\pi i n\cdot u_j}. Fiber Parseval followed by the usual Siegel mean formula for all nonzero lattice vectors gives the exact identity ∥Ff,j∥L2(Y2)2=∥f∥L2(ℝ2)2.\|F_{f,j}\|_{L^2(Y_2)}^2=\|f\|_{L^2(\mathbb R^2)}^2.(29.15) This is also the one-relative-mark case of [108, Proposition 6.11]. The level-two probability cover changes neither side: the integrand descends to the full modular quotient. By polarization and equivariance, ∫Y2Ff,j(x)¯Fh,j(atx)⁡dμ(x)=∫ℝ2f(x)¯h(atx)⁡dx.\boxed{ \int_{Y_2}\overline{F_{f,j}(x)}F_{h,j}(a_tx)\,\mathop{}\!\mathrm d\mu(x) =\int_{\mathbb R^2}\overline{f(x)}h(a_tx)\,\mathop{}\!\mathrm dx. }(29.16) Indeed the Fourier-side integral is ∫f̂(ξ)¯ĥ(a−tξ)⁡dξ\int\overline{\widehat f(\xi)}\widehat h(a_{-t}\xi)\mathop{}\!\mathrm d\xi. These identities extend by density to L2L^2 tests, using the isometric definition of the centered transform when the original sum is unavailable. For indicator tests the variance identity is also [109, Lemma 4]; the mean is Lemma 3 there. For different relative marks, the analogous centered cross-correlation vanishes by integration in the two independent torus fibers.

For the Gaussian f(x)=e−π|x|2f(x)=e^{-\pi|x|^2}, the formula becomes ∫Y2Ff,j¯Ff,j∘at⁡dμ=1(1+e2t)(1+e−2t)=12cosh⁡t.\int_{Y_2}\overline{F_{f,j}}F_{f,j}\circ a_t\,\mathop{}\!\mathrm d\mu =\frac1{\sqrt{(1+e^{2t})(1+e^{-2t})}} =\frac1{2\cosh t}.(29.17) Thus an e−te^{-t} correlation rate is attained by a natural smooth observable. More generally, substitution x1=e−tux_1=e^{-t}u in (29.16) gives, for Schwartz tests, et∫f¯(x)h(atx)⁡dx→(∫ℝf(0,y)¯dy)(∫ℝh(u,0)du).e^t\int\overline f(x)h(a_tx)\,\mathop{}\!\mathrm dx \longrightarrow \left(\int_\mathbb R\overline{f(0,y)}\,\mathop{}\!\mathrm dy\right) \left(\int_\mathbb Rh(u,0)\,\mathop{}\!\mathrm du\right). Schwartz domination justifies the limit.

29.8 Two marks: the rank-one divergence and its removal

The raw product Ff1,1Ff2,2F_{f_1,1}F_{f_2,2} need not be in L2L^2, even for smooth rapidly decreasing tests. Its conditional squared norm is A1(g)A2(g),Aj(g)=∑n≠0|f̂j(g−Tn)|2.A_1(g)A_2(g),\qquad A_j(g)=\sum_{n\ne0}|\widehat f_j(g^{-T}n)|^2. For Gaussians f1=f2=e−π|x|2f_1=f_2=e^{-\pi|x|^2} the divergence can be seen exactly. Write the collinear frequency pairs as (an,bn)(an,bn), where nn is primitive modulo sign and a,b∈ℤ\{0}a,b\in\mathbb Z\setminus\{0\}. Primitive Siegel unfolding gives the rank-one contribution 12ζ(2)∑a,b≠0∫ℝ2e−2π(a2+b2)|x|2⁡dx=14ζ(2)∑a,b≠01a2+b2=∞.\frac1{2\zeta(2)}\sum_{a,b\ne0} \int_{\mathbb R^2}e^{-2\pi(a^2+b^2)|x|^2}\,\mathop{}\!\mathrm dx =\frac1{4\zeta(2)}\sum_{a,b\ne0}\frac1{a^2+b^2} =\infty. The divergence is logarithmic. More generally, if both ∫fj\int f_j are nonzero, Fourier continuity gives Aj(g)≫ε−1A_j(g)\gg\varepsilon^{-1} when the dual lattice has a primitive vector of length ε≪1\varepsilon\ll1; the cusp measure is comparable to ε⁡dε\varepsilon\mathop{}\!\mathrm d\varepsilon. Centering has removed zero frequencies but leaves precisely this collinear-frequency obstruction.

There is an explicit L2L^2 replacement. Retain only the full-rank frequency matrices: Hf1,f2(2)(g,u1,u2)=∑n1,n2∈ℤ2det⁡(n1,n2)≠0f̂1(g−Tn1)f̂2(g−Tn2)e2πi(n1⋅u1+n2⋅u2).H_{f_1,f_2}^{(2)}(g,u_1,u_2) =\sum_{\substack{n_1,n_2\in\mathbb Z^2\\\det(n_1,n_2)\ne0}} \widehat f_1(g^{-T}n_1)\widehat f_2(g^{-T}n_2) e^{2\pi i(n_1\cdot u_1+n_2\cdot u_2)}.(29.18) The series is absolutely convergent on each fixed torus fiber. The rank restriction is invariant under the modular transition and under the SL2SL_2 action. It is a fiber Fourier projection, not a positivity-preserving deletion of physical segments. It is defined fiberwise before global integrability has been proved; no orthogonal projection on the global L2L^2 space is being applied to the raw product.

Here is a norm computation which proves global integrability. Put aj=|f̂j|2a_j=|\widehat f_j|^2 and let σ1(k)=∑d∣kd\sigma_1(k)=\sum_{d\mid k}d. Hermite normal form gives σ1(|k|)\sigma_1(|k|) full-modular orbits of integer matrices of determinant k≠0k\ne0, each with trivial stabilizer. Thus Parseval and unfolding give ∥Hf1,f2(2)∥22=∑k≠0σ1(|k|)Ik(a1,a2),\|H_{f_1,f_2}^{(2)}\|_2^2 =\sum_{k\ne0}\sigma_1(|k|)I_k(a_1,a_2),(29.19) where, with probability-covolume Haar normalization, Ik(a1,a2)=1|k|ζ(2)∫ℝ2\{0}∫ℝa1(x)a2(kJx|x|2+sx)⁡ds⁡dx,J(x1,x2)=(−x2,x1).I_k(a_1,a_2) =\frac1{|k|\zeta(2)} \int_{\mathbb R^2\setminus\{0\}}\int_\mathbb R a_1(x)a_2\left(\frac{kJx}{|x|^2}+sx\right) \,\mathop{}\!\mathrm ds\,\mathop{}\!\mathrm dx, \qquad J(x_1,x_2)=(-x_2,x_1). For completeness the Haar-coordinate formula behind this expression is g=[x,Jx|x|2+sx],⁡dg=1ζ(2)⁡dx⁡ds.g=\left[x,\frac{Jx}{|x|^2}+sx\right],\qquad \mathop{}\!\mathrm dg=\frac1{\zeta(2)}\,\mathop{}\!\mathrm dx\,\mathop{}\!\mathrm ds. The normalization follows by unfolding the primitive first-column Siegel transform: the stabilizer shifts ss by integers, and its mean is ζ(2)−1∫f\zeta(2)^{-1}\int f. At level two each full-rank modular orbit splits into six orbits while Haar measure is divided by six, leaving (29.19) unchanged.

An explicit Schwartz bound is convenient. Define A(f)=supξ∈ℝ2(1+|ξ|2)3/2|f̂(ξ)|.A(f)=\sup_{\xi\in\mathbb R^2} (1+|\xi|^2)^{3/2}|\widehat f(\xi)|. For k>0k>0 the orbit integral is bounded by I±k≤A(f1)2A(f2)2∫G(1+ktr⁡(gTg))−3⁡dg=π2A(f1)2A(f2)22ζ(2)k(1+2k)2.I_{\pm k}\leq A(f_1)^2A(f_2)^2 \int_G(1+k\operatorname{tr}(g^Tg))^{-3}\,\mathop{}\!\mathrm dg =\frac{\pi^2 A(f_1)^2A(f_2)^2} {2\zeta(2)k(1+2k)^2}. One way to verify the last integral is the Gaussian identity ∫Ge−αktr⁡(gTg)⁡dg=π2αkζ(2)e−2αk,α,k>0,\int_G e^{-\alpha k\operatorname{tr}(g^Tg)}\,\mathop{}\!\mathrm dg =\frac{\pi^2}{\alpha k\zeta(2)}e^{-2\alpha k}, \qquad \alpha,k>0, followed by the Laplace representation of (1+u)−3(1+u)^{-3}. Since ∑k≥1σ1(k)k−3=ζ(3)ζ(2)\sum_{k\geq1}\sigma_1(k)k^{-3}=\zeta(3)\zeta(2), ∥Hf1,f2(2)∥2≤π2ζ(3)A(f1)A(f2).\boxed{ \|H_{f_1,f_2}^{(2)}\|_2 \leq\frac\pi2\sqrt{\zeta(3)}\,A(f_1)A(f_2). }(29.20)

Let Df(x)=(Jx)⋅∇f(x)Df(x)=(Jx)\cdot\nabla f(x) and set B(f)=A(f)+A(Df)B(f)=A(f)+A(Df). Rotation commutes with the rank projection, so ΘHf1,f2(2)=HDf1,f2(2)+Hf1,Df2(2).\Theta H_{f_1,f_2}^{(2)} =H_{Df_1,f_2}^{(2)}+H_{f_1,Df_2}^{(2)}. Therefore 𝒮1(Hf1,f2(2))≤π2ζ(3)B(f1)B(f2).\mathcal S_1(H_{f_1,f_2}^{(2)}) \leq\frac\pi2\sqrt{\zeta(3)}\,B(f_1)B(f_2).(29.21) The preceding mixing theorem now applies, with constants controlled explicitly by the four test functions. In particular, |⟨Hf1,f2(2),Hh1,h2(2)∘at⟩|≤π2ζ(3)coth⁡π(1+t)e−tB(f1)B(f2)B(h1)B(h2).\left|\left\langle H_{f_1,f_2}^{(2)}, H_{h_1,h_2}^{(2)}\circ a_t\right\rangle\right| \leq\pi^2\zeta(3)\coth\pi\,(1+t)e^{-t} B(f_1)B(f_2)B(h_1)B(h_2). The mean is zero because every retained fiber frequency is nonzero.

29.9 An exact Gaussian two-mark correlation

Take f1=f2=f=e−π|x|2f_1=f_2=f=e^{-\pi|x|^2} and write H(2)=Hf,f(2)H^{(2)}=H_{f,f}^{(2)}. The determinant-orbit calculation can be done at every time: ⟨H(2),H(2)∘at⟩=πζ(2)cosh⁡t∑k=1∞σ1(k)ke−4πkcosh⁡t.\boxed{ \left\langle H^{(2)},H^{(2)}\circ a_t\right\rangle =\frac{\pi}{\zeta(2)\cosh t} \sum_{k=1}^\infty\frac{\sigma_1(k)}k e^{-4\pi k\cosh t}. }(29.22) To see this, the product of the Fourier coefficients is exp⁡[−πtr⁡(ZTCtZ)],Ct=I+a−2t,det⁡Z=k.\exp\bigl[-\pi\operatorname{tr}(Z^TC_tZ)\bigr], \qquad C_t=I+a_{-2t},\qquad \det Z=k. Since det⁡Ct=2cosh⁡t\sqrt{\det C_t}=2\cosh t, an SL2SL_2 change of variables reduces this to the radial Gaussian integral with α=2πcosh⁡t\alpha=2\pi\cosh t. Summing the positive and negative determinant orbits gives (29.22). At t=0t=0 this gives the exact squared norm. The identity ∑k≥1σ1(k)kqk=−∑n≥1log⁡(1−qn),0<q<1,\sum_{k\geq1}\frac{\sigma_1(k)}k q^k =-\sum_{n\geq1}\log(1-q^n),\qquad 0<q<1, also expresses the answer by the ordinary partition Euler product. The resulting correlation decays faster than exponentially in tt. This strong decay is a feature of these particular Gaussian tests; it is consistent with the e−te^{-t} one-mark example and the uniform tempered bound. For a finite evaluation, put q=e−4πcosh⁡tq=e^{-4\pi\cosh t}. The elementary bound σ1(k)/k≤k\sigma_1(k)/k\leq k gives the explicit truncation estimate 0≤∑k>Nσ1(k)kqk≤qN+1((N+1)−Nq)(1−q)2.0\leq\sum_{k>N}\frac{\sigma_1(k)}kq^k \leq\frac{q^{N+1}((N+1)-Nq)}{(1-q)^2}. This controls truncation of the analytic series. Floating-point evaluation of this bound is not an interval enclosure of the total computational error, and checking finitely many times is not a proof of the general mixing theorem.

29.10 Sharp boundaries and the rotation norm

The exact isometry also displays a regularity cost of sharp cutoffs. For a bounded convex body B⊂ℝ2B\subset\mathbb R^2, ∥F𝟏B,j∘kθ−F𝟏B,j∥22=area⁡(B△kθB).\|F_{\mathbf1_B,j}\circ k_\theta-F_{\mathbf1_B,j}\|_2^2 =\operatorname{area}(B\mathbin\triangle k_\theta B). The first variation of area under rotation, interpreted using the almost-everywhere outer normal of a convex boundary, gives area⁡(B△kθB)=|θ|∫∂B|Jx⋅n(x)|⁡ds+o(|θ|).\operatorname{area}(B\mathbin\triangle k_\theta B) =|\theta|\int_{\partial B}|Jx\cdot n(x)|\,\mathop{}\!\mathrm ds+o(|\theta|). Unless BB is rotationally invariant, this rules out a finite L2L^2 rotation derivative. Equivalently, the sharp transform has angular Sobolev regularity of order s<1/2s<1/2, but not s≥1/2s\geq1/2. The direct summation of tempered KK-type estimates requires an order greater than 1/21/2; the earlier theorem uses order one.

This is a generic boundary-regularity phenomenon. It is not by itself the more specific arithmetic obstruction from the unusually flat McShane–Rivin norm ball. It explains why a smooth correlation bound cannot be applied to a sharp norm-ball count while leaving the observable norm unchanged. Here the body BB is held fixed in physical Euclidean coordinates. A cutoff depending only on the hyperbolic length ℓX(α)\ell_X(\alpha) is unchanged when the quadratic differential is rotated over the same conformal surface XX, and does not incur this particular rotation loss. These two kinds of cutoffs must be distinguished. The fixed-surface primitive-lattice count, its Möbius inversion, and its unsmoothing error remain distinct from the affine mean-square identities proved here.

Finite computational companion.

The companion experiments/counting_mixing/ checks exact finite algebra and selected one- and two-mark Gaussian formulas at 6060 and 9090 decimal working digits, including rank-two series/orbit comparisons. Analytic truncation bounds do not certify quadrature or roundoff. These fixtures prove neither general mixing nor the nonsquare genus-two model hypotheses; finite collinear growth does not prove divergence. Its README links results, review and generation records.