Polyhedra, Moduli,
and Quantum Geometry

Web working edition

On this page
  1. 23.1 Broken geodesics and the effect of reconnection
  2. 23.2 The finite coefficient distributions
  3. 23.3 Length triangularity and the first atom
  4. 23.4 An exponential coefficient bound from bracelet positivity
  5. 23.5 Convergent heat and Laplace transforms
  6. 23.6 All iterates of a simple curve
  7. 23.7 A figure-eight family and an exact cancellation
  8. 23.8 Metric variation and the twist derivative
  9. 23.9 The logarithmic product and primitive extraction
  10. 23.10 Scope and further questions

23 Simple multicurve coefficients of a twisted trace formula

This section develops a spectral version of the crossing-resolution idea. We work throughout with a marked, closed, oriented hyperbolic surface XX of genus g≥2g\geq2 and curvature −1-1. All closed curves are unoriented, unless conjugacy classes are explicitly being counted with orientation. A multicurve may have parallel components; its length and number of components count those multiplicities. The empty multicurve is included. Cusps and boundary require additional trace terms and are not included in the assertions below.

The inputs are the twisted Selberg formula, the simple-multicurve basis of the character ring, and Dylan Thurston’s sign-coherence theorem for the bracelet basis. We prove the intervening length and coefficient estimates explicitly. The resulting construction is a deduction from these established inputs, with no novelty claim.

23.1 Broken geodesics and the effect of reconnection

Colin de Verdière’s construction approximates loop space by broken geodesics and applies stationary phase to a complex heat parametrix [19]. For an equal time subdivision the finite-dimensional energy is EN(x0,…,xN−1)=N∑i=0N−1d(xi,xi+1)2,xN=x0,E_N(x_0,\ldots,x_{N-1}) =N\sum_{i=0}^{N-1}d(x_i,x_{i+1})^2,\qquad x_N=x_0, with consecutive points within the injectivity radius. The critical equations match incoming and outgoing velocities; the critical value for a constant-speed closed geodesic is ℓγ2\ell_\gamma^2. On a hyperbolic surface every closed geodesic, including a self-intersecting one, has positive transverse length Hessian. For a unit-speed parametrization and a periodic normal variation f(s)N(s)f(s)N(s) this is δ2Lγ(f,f)=∫0ℓγ((f′)2+f2)⁡ds.\delta^2L_\gamma(f,f) =\int_0^{\ell_\gamma}\bigl((f')^2+f^2\bigr)\,\mathop{}\!\mathrm ds. Consequently simplicity cannot be selected by retaining only positive Hessians in the ordinary loop space.

Reconnection enlarges that space. In a small embedded disk of radius rr about a transverse crossing of angle θ∈(0,π)\theta\in(0,\pi), the two original diameters have total length 4r4r. The two noncrossing pairings have lengths L0=4arsinh⁡(sinh⁡rsin⁡(θ/2)),L∞=4arsinh⁡(sinh⁡rcos⁡(θ/2)).(23.1)\begin{aligned} L_0&=4\mathop{\mathrm{arsinh}}\bigl(\sinh r\sin(\theta/2)\bigr),\\ L_\infty&=4\mathop{\mathrm{arsinh}}\bigl(\sinh r\cos(\theta/2)\bigr). \end{aligned} Both are strictly less than 4r4r. In particular, 4r−L0=4r(1−sin⁡(θ/2))+O(r3).4r-L_0=4r\bigl(1-\sin(\theta/2)\bigr)+O(r^3). The descent is first order in the surgery radius. Tightening the resolved components to geodesics, and discarding contractible components, cannot increase total length. This is a change of connectivity, not a negative direction of the original Hessian. An unconstrained global minimization would select the empty cycle; we instead keep the different multicurves as algebraic observables.

23.2 The finite coefficient distributions

Let 𝒞(S)\mathscr C(S) be the set of essential simple multicurves on the underlying marked surface. For ρ:π1(S)→SU(2)\rho:\pi_1(S)\to\mathrm{SU}(2) set TC(ρ)=∏α⊂Ctr⁡ρ(α),T⌀=1.T_C(\rho)=\prod_{\alpha\subset C}\mathop{\mathrm{tr}}\rho(\alpha), \qquad T_\varnothing=1.(23.2) The phrase “unnegated trace basis” refers to this sign convention, not to positivity of the function on SU(2)\mathrm{SU}(2). The functions TCT_C are linearly independent on the unitary character variety and form a basis of the complex character ring after complexification; see Charles–Marché [20, Theorem 1.1]. The usual Kauffman convention uses −tr-\mathop{\mathrm{tr}} per component, so differs by a basis sign. The relation tr⁡Atr⁡B=tr⁡(AB)+tr⁡(AB−1)\mathop{\mathrm{tr}}A\,\mathop{\mathrm{tr}}B=\mathop{\mathrm{tr}}(AB)+\mathop{\mathrm{tr}}(AB^{-1}) resolves crossings. For every nontrivial closed curve γ\gamma there is a unique finite expansion tr⁡ρ(γ)=∑C∈𝒞(S)aC(γ)TC(ρ),aC(γ)∈ℤ.\mathop{\mathrm{tr}}\rho(\gamma)=\sum_{C\in\mathscr C(S)}a_C(\gamma)T_C(\rho), \qquad a_C(\gamma)\in\mathbb Z.(23.3) These integers depend only on topology, not on XX.

Let Δρ\Delta_\rho be the rank-two flat unitary bundle Laplacian, and write λj(ρ)=1/4+rj(ρ)2\lambda_j(\rho)=1/4+r_j(\rho)^2. For an even φ∈Cc∞(ℝ)\varphi\in C_c^\infty(\mathbb R) put h(r)=∫ℝφ(L)eirL⁡dL,h(r)=\int_\mathbb R\varphi(L)e^{irL}\,\mathop{}\!\mathrm dL, and define 𝒯φX(ρ)=∑jh(rj(ρ))−2Area⁡(X)4π∫ℝrh(r)tanh⁡(πr)⁡dr.\begin{aligned} \mathscr T^X_\varphi(\rho) &=\sum_j h(r_j(\rho)) -\frac{2\operatorname{Area}(X)}{4\pi} \int_\mathbb Rr h(r)\tanh(\pi r)\,\mathop{}\!\mathrm dr. \end{aligned} With the two orientations grouped, the twisted Selberg formula is 𝒯φX(ρ)=∑δ∈𝒫un∑k≥1ℓδsinh⁡(kℓδ/2)φ(kℓδ)tr⁡ρ(δk).\mathscr T^X_\varphi(\rho) =\sum_{\delta\in\mathcal P_{\mathrm{un}}} \sum_{k\geq1} \frac{\ell_\delta}{\sinh(k\ell_\delta/2)} \varphi(k\ell_\delta)\mathop{\mathrm{tr}}\rho(\delta^k).(23.4) Here 𝒫un\mathcal P_{\mathrm{un}} consists of unoriented primitive closed geodesics. Naud’s statement [21, Theorem 2.2] sums over oriented primitive conjugacy classes with denominator 2sinh2\sinh; grouping inverses gives (23.4). The compact length support makes the geometric sum finite.

Definition 23.1 (Coefficient distribution). For C∈𝒞(S)C\in\mathscr C(S) define the locally finite signed measure on (0,∞)(0,\infty) σCX=∑δ∈𝒫un∑k≥1ℓδaC(δk)sinh⁡(kℓδ/2)δkℓδ.\boxed{ \sigma_C^X =\sum_{\delta\in\mathcal P_{\mathrm{un}}}\sum_{k\geq1} \frac{\ell_\delta a_C(\delta^k)} {\sinh(k\ell_\delta/2)}\,\delta_{k\ell_\delta}. }(23.5) Thus [TC]𝒯φX=⟨σCX,φ⟩.[T_C]\mathscr T^X_\varphi=\langle\sigma_C^X,\varphi\rangle. The notation δL\delta_L denotes the unit point mass at LL.

Local finiteness follows from the finiteness of the closed-geodesic length set below any bound. The identity term in the trace formula is a constant function of ρ\rho, hence affects only the empty multicurve coefficient. The input is the labelled family ρ↦Spec⁡Δρ\rho\mapsto\operatorname{Spec}\Delta_\rho on the marked character variety. A scalar spectrum alone does not retain these coefficients.

23.3 Length triangularity and the first atom

Lemma 23.2 (Length ordering of resolutions). If δ\delta is primitive and k≥1k\geq1, then aC(δk)≠0⇒ℓX(C)≤kℓX(δ).a_C(\delta^k)\neq0\quad\Longrightarrow\quad \ell_X(C)\leq k\ell_X(\delta). Equality can occur only when δ\delta is simple and C=kδC=k\delta, the multicurve consisting of kk parallel copies. Its coefficient in that case is 11. If δ\delta is nonsimple, the inequality is strict.

Proof. For a primitive nonsimple geodesic, completely resolve its crossings. Each state has a representative obtained by reconnecting the original arcs in small disks. Formula (23.1) gives a strict total saving; tightening and removing contractible loops cannot increase length. At a multiple intersection, resolve the endpoints in a small convex disk by a noncrossing pairing. The original pairing uses intersecting diameters, so at least one new chord is shorter than a diameter. The same strict conclusion holds. Coefficient cancellation cannot introduce a state not produced by a resolution.

The product version is also useful: every multicurve in the expansion of TC1⋯TCjT_{C_1}\cdots T_{C_j} has length at most ∑iℓX(Ci)\sum_i\ell_X(C_i). Superimpose representatives, separate parallel components, resolve crossings, and tighten.

For powers use the Chebyshev polynomials P0(x)=2,P1(x)=x,Pk+1(x)=xPk(x)−Pk−1(x),tr⁡(Ak)=Pk(tr⁡A).P_0(x)=2,\quad P_1(x)=x,\quad P_{k+1}(x)=xP_k(x)-P_{k-1}(x), \qquad \mathop{\mathrm{tr}}(A^k)=P_k(\mathop{\mathrm{tr}}A). If δ\delta is simple, Pk(Tδ)P_k(T_\delta) contains only parallel copies mδm\delta with m≤km\leq k, and its leading term is TδkT_\delta^k. If δ\delta is nonsimple, every basis term in tr⁡ρ(δ)\mathop{\mathrm{tr}}\rho(\delta) has length strictly below ℓδ\ell_\delta. Apply the product inequality to each monomial in PkP_k: every nonconstant output has length strictly below kℓδk\ell_\delta; the constant term has length zero. ◻

Proposition 23.3 (First atom for a simple support). For a simple primitive curve α\alpha and m≥1m\geq1, σmαX=ℓαsinh⁡(mℓα/2)δmℓα+a measure supported strictly above mℓα.\boxed{ \sigma_{m\alpha}^X =\frac{\ell_\alpha}{\sinh(m\ell_\alpha/2)} \delta_{m\ell_\alpha} +\text{a measure supported strictly above }m\ell_\alpha. }(23.6) For each fixed CC the gap is locally uniform in XX, and uniform on compact subsets of marked Teichmüller space. If CC has two or more nonparallel support components, every contributing orbit has length strictly greater than ℓX(C)\ell_X(C); this statement does not assert that the resulting signed measure is nonzero.

Proof. Lemma 23.2 gives the lower support bound. At mℓαm\ell_\alpha the only contributing pair is (δ,k)=(α,m)(\delta,k)=(\alpha,m), with coefficient 11. This gives the stated mass, whose numerator is the primitive length. Every other contributing pair has strictly greater length. Local finiteness gives a gap at each XX.

For uniformity on a compact set of marked metrics, use uniform bilipschitz comparison with one reference metric. The target lengths are bounded, and only finitely many orbit labels can have lengths within one unit above them somewhere on the compact set. The strictly positive differences for those finitely many labels have a positive minimum; all remaining labels have a gap at least one. The argument also applies to nonparallel CC, without requiring a nonzero first mass after cancellation. ◻

For m=1m=1 the first atom determines both the simple length and its normalization. Choose a fixed smooth cutoff χ\chi, equal to one near ℓα\ell_\alpha and supported away from every other atom in a small neighborhood of XX. Then ℓX(α)=⟨σαX,Lχ(L)⟩⟨σαX,χ(L)⟩,⟨σαX,χ⟩=ℓαsinh⁡(ℓα/2).\boxed{ \ell_X(\alpha) =\frac{\langle\sigma_\alpha^X,L\chi(L)\rangle} {\langle\sigma_\alpha^X,\chi(L)\rangle},\qquad \langle\sigma_\alpha^X,\chi\rangle =\frac{\ell_\alpha}{\sinh(\ell_\alpha/2)}. }(23.7) This is an exact localized identity, with no limiting interchange. Figure 7 relates the local reconnections to this separate first-atom statement for a simple support.

Crossing resolution and the first coefficient atom
Figure 7. Reconnection and coefficient extraction are different operations. The upper pictures show the two noncrossing pairings in a hyperbolic metric disk; both shorten the original pair of diameters. Reconnection changes connectivity and is not a negative direction of the length Hessian within the original loop space. For a simple primitive support α\alpha, the lower schematic records the exact first atom of σαX\sigma_\alpha^X at ℓα\ell_\alpha, with the stated mass; all remaining atoms lie strictly above that length. Their depicted positions and weights are illustrative, not numerical data.

23.4 An exponential coefficient bound from bracelet positivity

The passage from compactly supported tests to heat kernels needs a growth estimate. A direct crossing count gives an exponential in ℓγ2\ell_\gamma^2 and only a small-time conclusion. Bracelet positivity gives a stronger estimate.

At skein parameter +1+1, Dylan Thurston proves that every diagram has a sign-coherent expansion in the bracelet basis [22, Theorem 2]. A bracelet supported on a simple curve β\beta with winding nn is the Chebyshev element Pn(zβ)P_n(z_\beta) [22, Proposition 4.4]. Distinct bracelet support curves are disjoint; their windings are included in the length ℓX(D)=∑iniℓX(βi)\ell_X(D)=\sum_i n_i\ell_X(\beta_i). The canonical twisted Fuchsian evaluation associated to XX is positive on bracelet basis elements: BD(X)=∏i2cosh⁡(niℓX(βi)/2)>0.B_D(X)=\prod_i2\cosh\bigl(n_i\ell_X(\beta_i)/2\bigr)>0. For a one-component diagram representing γ\gamma, its evaluation has absolute value 2cosh⁡(ℓX(γ)/2)2\cosh(\ell_X(\gamma)/2). Sign coherence therefore gives ∑D|bD(γ)|BD(X)=2cosh⁡(ℓX(γ)/2).\sum_D |b_D(\gamma)| B_D(X) =2\cosh\bigl(\ell_X(\gamma)/2\bigr).(23.8) A spin structure identifies the +1+1 and −1-1 skein algebras by signs on diagrams. Changing to (23.2) also changes only basis signs. On a simple generator this replaces zβz_\beta by ϵβTβ\epsilon_\beta T_\beta, with ϵβ∈{1,−1}\epsilon_\beta\in\{1,-1\}; the identity Pn(ϵx)=ϵnPn(x)P_n(\epsilon x)=\epsilon^nP_n(x) shows that the bracelet-to-multicurve coefficient magnitudes are unchanged. All the absolute estimates below are consequently valid in our unnegated trace basis. The positive Fuchsian evaluation is an auxiliary algebraic estimate; the Laplacians defining σCX\sigma_C^X still use unitary twists.

Proposition 23.4 (Fixed-coefficient bound). Write C=∑j=1qrjβjC=\sum_{j=1}^q r_j\beta_j with distinct support curves. For r≥1r\geq1 and ℓ>0\ell>0 put Mr(ℓ)=supn≥1|[xr]Pn(x)|2cosh⁡(nℓ/2),KC,X=∏j=1qMrj(ℓX(βj)).M_r(\ell)=\sup_{n\geq1} \frac{|[x^r]P_n(x)|}{2\cosh(n\ell/2)},\qquad K_{C,X}=\prod_{j=1}^q M_{r_j}(\ell_X(\beta_j)). Set K⌀,X=1K_{\varnothing,X}=1. These constants are finite, and |aC(γ)|≤KC,X2cosh⁡(ℓX(γ)/2)\boxed{ |a_C(\gamma)|\leq K_{C,X}\,2\cosh\bigl(\ell_X(\gamma)/2\bigr) }(23.9) for every nontrivial closed curve, including iterates.

Proof. When n≥r≥1n\geq r\geq1 have the same parity, |[xr]Pn(x)|=2nn+r((n+r)/2r)≤2nrr!;|[x^r]P_n(x)| =\frac{2n}{n+r}\binom{(n+r)/2}{r} \leq\frac{2n^r}{r!}; otherwise the coefficient is zero. Consequently Mr(ℓ)≤2r!(2reℓ)r<∞.M_r(\ell)\leq \frac{2}{r!}\left(\frac{2r}{e\ell}\right)^r<\infty.(23.10) To contribute to CC, a bracelet element must contain every βj\beta_j in its support. Additional support curves must contribute constant terms, which are either zero or ±2\pm2. Their absolute constant coefficient divided by 2cosh⁡(nℓ/2)2\cosh(n\ell/2) is at most one. It follows that |[TC]BD|≤KC,XBD(X).|[T_C]B_D|\leq K_{C,X}B_D(X). Apply this bound to (23.8) and use the triangle inequality. Spin and basis signs do not affect the estimate. ◻

For a single simple α\alpha one can use Kα,X=supn≥1n oddn2cosh⁡(nℓα/2)≤2eℓα.K_{\alpha,X} =\sup_{\substack{n\geq1\\n\text{ odd}}} \frac{n}{2\cosh(n\ell_\alpha/2)} \leq\frac{2}{e\ell_\alpha}. The estimates are uniform for a fixed target CC on compact subsets of Teichmüller space. Related positivity and Chebyshev estimates for analytic lamination expansions on punctured planar surfaces appear in Basok–Chelkak [23, Section 4.4, Theorems 4.7 and 4.9]. Their full analytic expansion theorem is not needed here.

Corollary 23.5 (Growth of the signed measure). For each fixed X,CX,C there is a constant BC,XB_{C,X} such that |σCX|([0,R])≤BC,X(1+R)eR.|\sigma_C^X|([0,R])\leq B_{C,X}(1+R)e^R.(23.11)

Proof. Put ε=sys⁡(X)>0\varepsilon=\operatorname{sys}(X)>0 and L=kℓδL=k\ell_\delta. Each orbital summand has absolute mass at most ℓδsinh⁡(L/2)KC,X2cosh⁡(L/2)≤2KC,Xcoth⁡(ε/2)L.\frac{\ell_\delta}{\sinh(L/2)} K_{C,X}2\cosh(L/2) \leq2K_{C,X}\coth(\varepsilon/2)L. The number of primitive/iterate classes of length at most RR is OX(eR)O_X(e^R). One elementary proof conjugates each class so its axis passes a bounded distance from a fixed point in ℍ2\mathbb H^2; the chosen element then moves that point by at most R+OX(1)R+O_X(1). Orbit points are uniformly separated, and hyperbolic disk area is O(eR)O(e^R). Summing the displayed mass bound proves the claim. No prime-geodesic asymptotic is required. ◻

23.5 Convergent heat and Laplace transforms

Define GCX(t)=∫0∞e−L2/(4t)⁡dσCX(L),HCX(t)=e−t/44πtGCX(t),(23.12)DCX(s)=∫0∞e−sL⁡dσCX(L).(23.13)\begin{aligned} G_C^X(t)&=\int_0^\infty e^{-L^2/(4t)}\,\mathop{}\!\mathrm d\sigma_C^X(L), &H_C^X(t)&=\frac{e^{-t/4}}{\sqrt{4\pi t}}G_C^X(t), \\ D_C^X(s)&=\int_0^\infty e^{-sL}\,\mathop{}\!\mathrm d\sigma_C^X(L). && \end{aligned} By (23.11), GCG_C and HCH_C are absolutely convergent for every t>0t>0, and holomorphic for ℜt>0\Re t>0 with the principal square root. The Laplace transform is absolutely convergent and holomorphic for ℜs>1\Re s>1. Differentiation in these variables is permitted locally uniformly in the stated domains.

There is also a convergent expansion of the whole heat trace. For a finite character polynomial set ∥f∥2=∑C|[TC]f|2#C,cX=log⁡(1+2)sys⁡(X).\|f\|_2=\sum_C|[T_C]f|\,2^{\#C},\qquad c_X=\frac{\log(1+\sqrt2)}{\operatorname{sys}(X)}. The subscript 22 denotes the trace bound |tr⁡ρ(α)|≤2|\mathop{\mathrm{tr}}\rho(\alpha)|\leq2; this is a weighted coefficient norm, not an L2L^2 norm. The Chebyshev recurrence gives ∑r|[xr]Pn(x)|2r≤2(1+2)n.\sum_r|[x^r]P_n(x)|2^r\leq2(1+\sqrt2)^n. A bracelet has at most 3g−33g-3 distinct support curves, so ∥BD∥2≤23g−3ecXℓX(D)\|B_D\|_2\leq2^{3g-3}e^{c_X\ell_X(D)}. In (23.8) the nonzero coefficients are integers, and BD(X)≥eℓX(D)/2B_D(X)\geq e^{\ell_X(D)/2}. Hence ℓX(D)≤ℓX(γ)+2log⁡2\ell_X(D)\leq\ell_X(\gamma)+2\log2 for every occurring DD, and ∑D|bD|≤2eℓX(γ)/2\sum_D|b_D|\leq2e^{\ell_X(\gamma)/2}. Therefore ∥tr⁡γ∥2≤CXe(1/2+cX)ℓX(γ),CX=23g−2e2cXlog⁡2.\|\mathop{\mathrm{tr}}_\gamma\|_2\leq C_X e^{(1/2+c_X)\ell_X(\gamma)},\qquad C_X=2^{3g-2}e^{2c_X\log2}.(23.14) Gaussian decay dominates this exponential bound and the orbit count. It follows that Tr⁡e−tΔρ−2IX(t)=∑C∈𝒞(S)HCX(t)TC(ρ),t>0,\boxed{ \operatorname{Tr}e^{-t\Delta_\rho}-2I_X(t) =\sum_{C\in\mathscr C(S)}H_C^X(t)T_C(\rho),\qquad t>0, }(23.15) where IX(t)=Area⁡(X)4π∫ℝrtanh⁡(πr)e−t(r2+1/4)⁡dr,∑C|HCX(t)|2#C<∞.I_X(t)=\frac{\operatorname{Area}(X)}{4\pi} \int_\mathbb Rr\tanh(\pi r)e^{-t(r^2+1/4)}\,\mathop{}\!\mathrm dr, \qquad \sum_C|H_C^X(t)|2^{\#C}<\infty. The series is uniform on the unitary character variety and locally uniform for ℜt>0\Re t>0. It is obtained by taking the limit of the finite length-cutoff expansions; the preceding absolute estimate justifies the rearrangement. For a compact set of such tt, ℜ(1/t)\Re(1/t) has a positive lower bound, so the same Gaussian majorant proves local uniformity and termwise holomorphy. On the spectral side, smooth cutoffs of the Gaussian converge in the Schwartz topology on the real spectral axis and at the finitely many imaginary spectral parameters; the compact-surface Weyl bound justifies passing to the heat trace for each unitary twist.

This construction does not assert uniqueness among arbitrary infinite expansions of a continuous function. For example, ∥tr⁡ρ(α2n+1)∥∞≤2\|\mathop{\mathrm{tr}}\rho(\alpha^{2n+1})\|_\infty\leq2, while its TαT_\alpha coefficient is (−1)n(2n+1)(-1)^n(2n+1). Coefficient extraction is therefore not bounded for the uniform norm on unitary characters. Our coefficients are fixed by the finite spectral tests first; the estimates then justify their heat transform.

Proposition 23.6 (Simple-curve semiclassical term). Let α\alpha be simple and ℓ=ℓX(α)\ell=\ell_X(\alpha). Choose b>ℓb>\ell so that every other orbital contribution to σαX\sigma_\alpha^X has length at least bb. Then GαX(t)=ℓsinh⁡(ℓ/2)e−ℓ2/(4t)+O(e−b2/(4t)),HαX(t)=e−t/44πtℓsinh⁡(ℓ/2)e−ℓ2/(4t)[1+O(e−(b2−ℓ2)/(4t))],(23.16)DαX(s)=ℓsinh⁡(ℓ/2)e−sℓ[1+O(e−s(b−ℓ))](s→+∞).(23.17)\begin{aligned} G_\alpha^X(t) &=\frac{\ell}{\sinh(\ell/2)}e^{-\ell^2/(4t)} +O\bigl(e^{-b^2/(4t)}\bigr),\\ H_\alpha^X(t) &=\frac{e^{-t/4}}{\sqrt{4\pi t}} \frac{\ell}{\sinh(\ell/2)}e^{-\ell^2/(4t)} \left[1+O\!\left(e^{-(b^2-\ell^2)/(4t)}\right)\right], \\ D_\alpha^X(s) &=\frac{\ell}{\sinh(\ell/2)}e^{-s\ell} \left[1+O\bigl(e^{-s(b-\ell)}\bigr)\right] \quad(s\to+\infty). \end{aligned} In particular, ℓ2=−limt↓04tlog⁡GαX(t),ℓ=−lims→+∞s−1log⁡DαX(s).\boxed{ \ell^2=-\lim_{t\downarrow0}4t\log G_\alpha^X(t),\qquad \ell=-\lim_{s\to+\infty}s^{-1}\log D_\alpha^X(s). }(23.18) Both logarithms are real for sufficiently small tt, respectively sufficiently large positive ss.

Proof. Separate the first atom in (23.6). For its remaining measure τ\tau, supported on [b,∞)[b,\infty), fix t0>0t_0>0. For 0<t≤t00<t\leq t_0, ∫e−L2/(4t)⁡d|τ|(L)≤e−b2/(4t)+b2/(4t0)∫e−L2/(4t0)⁡d|τ|(L).\int e^{-L^2/(4t)}\,\mathop{}\!\mathrm d|\tau|(L) \leq e^{-b^2/(4t)+b^2/(4t_0)} \int e^{-L^2/(4t_0)}\,\mathop{}\!\mathrm d|\tau|(L). The last integral is finite. The same argument with a fixed real s0>1s_0>1 gives the Laplace remainder. The leading masses are positive, so the leading terms eventually dominate the signed tails. ◻

After removing the explicit e−t/4e^{-t/4} factor, the correction to the single-geodesic term is exponentially small, rather than a new power series in tt. Thus the constant-curvature heat coefficient has an exact isolated-orbit prefactor, followed by contributions at larger length actions. For mαm\alpha replace the action by (mℓα)2/(4t)(m\ell_\alpha)^2/(4t) and the mass by ℓα/sinh⁡(mℓα/2)\ell_\alpha/\sinh(m\ell_\alpha/2).

23.6 All iterates of a simple curve

Write σCX,[δ]\sigma_C^{X,[\delta]} for the contribution of one primitive unoriented class δ\delta and all of its iterates. The polynomial identity Pk(x)=∑j=0⌊k/2⌋(−1)jk(k−j−1)!j!(k−2j)!xk−2j,k≥1,P_k(x)=\sum_{j=0}^{\lfloor k/2\rfloor} (-1)^j\frac{k(k-j-1)!}{j!(k-2j)!}x^{k-2j},\qquad k\geq1,(23.19) follows from the recurrence. If α\alpha is simple, the coefficient of mαm\alpha vanishes unless k=m+2jk=m+2j, and for m≥1m\geq1 is bm+2j,m=(−1)j(m+2j)(m+j−1)!j!m!.b_{m+2j,m}=(-1)^j \frac{(m+2j)(m+j-1)!}{j!\,m!}.(23.20) In particular, with ℓ=ℓα\ell=\ell_\alpha, σαX,[α]=ℓ∑j≥0(−1)j(2j+1)sinh⁡((2j+1)ℓ/2)δ(2j+1)ℓ.\boxed{ \sigma_\alpha^{X,[\alpha]} =\ell\sum_{j\geq0} \frac{(-1)^j(2j+1)}{\sinh((2j+1)\ell/2)} \delta_{(2j+1)\ell}. }(23.21) The first masses have multipliers 1,−3,5,−7,…1,-3,5,-7,\ldots at ℓ,3ℓ,5ℓ,7ℓ,…\ell,3\ell,5\ell,7\ell,\ldots. For 2α2\alpha the multipliers are 1,−4,9,−16,…1,-4,9,-16,\ldots at 2ℓ,4ℓ,6ℓ,8ℓ,…2\ell,4\ell,6\ell,8\ell,\ldots. These are explicit later atoms, not changes to the first-atom mass.

For m≥1m\geq1 the coefficient generating function is Bm(q):=∑k≥1bk,mqk=qm(1−q2)(1+q2)m+1.B_m(q):=\sum_{k\geq1}b_{k,m}q^k =\frac{q^m(1-q^2)}{(1+q^2)^{m+1}}.(23.22) Indeed, ∑k≥1Pk(x)qk=(xq−2q2)/(1−xq+q2)\sum_{k\geq1}P_k(x)q^k=(xq-2q^2)/(1-xq+q^2), and extracting xmx^m gives this formula. Expanding the hyperbolic denominator yields, for ℜs>−1/2\Re s>-1/2, DmαX,[α](s)=2ℓ∑n≥0Bm(qn),qn=e−(s+n+1/2)ℓ,=21−mℓ∑n≥0sech⁡m(un)tanh⁡(un),un=(s+n+1/2)ℓ.(23.23)\begin{aligned} D_{m\alpha}^{X,[\alpha]}(s) &=2\ell\sum_{n\geq0}B_m(q_n), &q_n&=e^{-(s+n+1/2)\ell},\\ &=2^{1-m}\ell\sum_{n\geq0} \mathop{\mathrm{sech}}^m(u_n)\tanh(u_n), &u_n&=(s+n+1/2)\ell. \end{aligned} This larger convergence domain concerns the one fixed primitive family; the proved half-plane for the full DCD_C is ℜs>1\Re s>1.

The heat contribution of the same family is HαX,[α](t)=e−t/44πtℓe−ℓ2/(4t)sinh⁡(ℓ/2)×[1−3sinh⁡(ℓ/2)sinh⁡(3ℓ/2)e−2ℓ2/t+5sinh⁡(ℓ/2)sinh⁡(5ℓ/2)e−6ℓ2/t−⋯].(23.24)\begin{aligned} H_\alpha^{X,[\alpha]}(t) ={}&\frac{e^{-t/4}}{\sqrt{4\pi t}} \frac{\ell e^{-\ell^2/(4t)}}{\sinh(\ell/2)}\\ &\mathrel{}\times\left[ 1-3\frac{\sinh(\ell/2)}{\sinh(3\ell/2)}e^{-2\ell^2/t} +5\frac{\sinh(\ell/2)}{\sinh(5\ell/2)}e^{-6\ell^2/t} -\cdots\right]. \end{aligned} Other primitive classes may supply corrections before these particular iterate actions. Once ℓα\ell_\alpha is recovered, subtracting the j≥1j\geq1 terms of (23.21) removes exactly the nonprimitive powers of the simple curve α\alpha; it leaves powers of nonsimple curves that resolve with an α\alpha coefficient.

23.7 A figure-eight family and an exact cancellation

Take an embedded geodesic pair of pants with three distinct essential cuffs a,b,ca,b,c, and c=(ab)−1c=(ab)^{-1}. Set x=Ta,y=Tb,z=Tc,γ=a−1b,Tγ=xy−z.x=T_a,\qquad y=T_b,\qquad z=T_c,\qquad \gamma=a^{-1}b, \qquad T_\gamma=xy-z. These are functions of arbitrary unitary holonomy. The minus sign therefore differs from the positive hyperbolic length convention in (18.1). Powers are given exactly by Tγk=Pk(xy−z)T_{\gamma^k}=P_k(xy-z). For example, Tγ=xy−z,Tγ2=x2y2−2xyz+z2−2,Tγ3=x3y3−3x2y2z+3xyz2−z3−3xy+3z.\begin{aligned} T_\gamma&=xy-z,\\ T_{\gamma^2}&=x^2y^2-2xyz+z^2-2,\\ T_{\gamma^3}&=x^3y^3-3x^2y^2z+3xyz^2-z^3-3xy+3z. \end{aligned} Only multicurves C=ra+rb+mcC=ra+rb+mc occur. Put d=r+m≥1d=r+m\geq1. Their coefficients vanish unless k=d+2jk=d+2j, and in that case ara+rb+mc(γk)=(−1)j+mk(k−j−1)!j!r!m!=(−1)m(r+mr)bk,r+m.\boxed{ a_{ra+rb+mc}(\gamma^k) =(-1)^{j+m}\frac{k(k-j-1)!}{j!\,r!\,m!} =(-1)^m\binom{r+m}{r}b_{k,r+m}. }(23.25) This follows by substituting xy−zxy-z into (23.19) and expanding the relevant monomial. Consequently ∑k≥1ara+rb+mc(γk)qk=(−1)m(r+mr)qr+m(1−q2)(1+q2)r+m+1.\sum_{k\geq1}a_{ra+rb+mc}(\gamma^k)q^k =(-1)^m\binom{r+m}{r} \frac{q^{r+m}(1-q^2)}{(1+q^2)^{r+m+1}}.

For r+m=1r+m=1 the two target coefficients are opposite for every iterate. If L=ℓγL=\ell_\gamma, then σa+bX,[γ]=−σcX,[γ]=L∑j≥0(−1)j(2j+1)sinh⁡((2j+1)L/2)δ(2j+1)L.\boxed{ \sigma_{a+b}^{X,[\gamma]} =-\sigma_c^{X,[\gamma]} =L\sum_{j\geq0} \frac{(-1)^j(2j+1)}{\sinh((2j+1)L/2)} \delta_{(2j+1)L}. }(23.26) Thus σcX+σa+bX\sigma_c^X+\sigma_{a+b}^X cancels the entire chosen figure-eight family, including all its powers. This does not assert that the full sum vanishes: other primitive classes contribute to both coefficients. It gives an exact family-by-family cancellation which is preserved by either transform in (23.12)–(23.13).

The distinct-cuff assumption makes the monomials above distinct multicurve basis elements. If cuff images coincide, as in some one-holed-torus configurations of Section 18, the corresponding variables must first be identified and equal basis terms collected; the displayed coefficient indexing cannot be carried over unchanged.

23.8 Metric variation and the twist derivative

For a fixed topological CC, the map X↦σCXX\mapsto\sigma_C^X is smooth as a distribution-valued map. After applying any compactly supported test function, it is locally a finite sum of smooth length functions. It is not asserted to be real analytic for arbitrary smooth test functions.

For an orbital label (δ,k)(\delta,k) write L=kℓδ,w=ℓδcsch⁡(L/2).L=k\ell_\delta,\qquad w=\ell_\delta\operatorname{csch}(L/2). For a tangent vector vv to Teichmüller space one has DvσCX=∑δ,kaC(δk)[(Dvw)δL−w(DvL)δL′],D_v\sigma_C^X =\sum_{\delta,k}a_C(\delta^k) \left[(D_vw)\delta_L-w(D_vL)\delta'_L\right],(23.27) where ⟨δL′,φ⟩=−φ′(L)\langle\delta'_L,\varphi\rangle=-\varphi'(L) and DvL=kDvℓδ,Dvw=csch⁡(L/2)(1−L2cothL2)Dvℓδ.D_vL=kD_v\ell_\delta,\qquad D_vw=\operatorname{csch}(L/2) \left(1-\frac L2\coth\frac L2\right)D_v\ell_\delta. The coefficients aCa_C are constant under metric variation. This formula describes movement of the atoms as well as variation of their masses; omitting the δL′\delta'_L term would lose the length variation itself.

For a twist of displacement τβ\tau_\beta about a simple geodesic β\beta, Wolpert’s cosine formula gives ∂ℓδ∂τβ=∑p∈δ∩βcos⁡θp,\frac{\partial\ell_\delta}{\partial\tau_\beta} =\sum_{p\in\delta\cap\beta}\cos\theta_p,(23.28) with the consistent right-twist angle convention and intersections counted with multiplicity [24]. Substitution in (23.27) gives an exact variation formula for every coefficient distribution. For a simple α\alpha, one can equivalently differentiate the localized ratio (23.7); it recovers (23.28) from two spectral coefficient moments, without an infinite limit. The first-variation argument also applies when δ\delta is immersed: cut it at its crossings with β\beta. Interior variations of the geodesic arcs cancel by stationarity; the relative tangential displacement at each paired endpoint contributes the projection cos⁡θp\cos\theta_p. Summing over the passages through the cut gives the formula with precisely the stated multiplicities. If δ\delta is a power of β\beta, its twist derivative is zero and the sum is interpreted as empty.

The compact-set gap in Proposition 23.3 also permits differentiating the leading heat asymptotic locally. The leading logarithmic response is Dvlog⁡HαX(t)=[−ℓα2t+1ℓα−12cothℓα2]Dvℓα+O(t−1e−c/t)D_v\log H_\alpha^X(t) =\left[-\frac{\ell_\alpha}{2t} +\frac1{\ell_\alpha} -\frac12\coth\frac{\ell_\alpha}{2}\right] D_v\ell_\alpha +O\bigl(t^{-1}e^{-c/t}\bigr)(23.29) for some local c>0c>0. To justify the differentiated tail, metric length derivatives on a compact parameter neighborhood are bounded by constants times the length, while (23.9) is locally uniform. Gaussian decay absorbs the additional polynomial factors.

23.9 The logarithmic product and primitive extraction

A useful connection to the Selberg product keeps the orientation normalization explicit. Define, for ℜv>1\Re v>1, 𝒵ρX(v)=∏δ∈𝒫un∏n≥0det⁡(I−e−(v+n)ℓδρ(δ)).\mathcal Z_\rho^X(v) =\prod_{\delta\in\mathcal P_{\mathrm{un}}}\prod_{n\geq0} \det\left(I-e^{-(v+n)\ell_\delta}\rho(\delta)\right).(23.30) This definition uses unoriented primitives; a product over both oriented conjugacy classes squares it. Expanding the logarithm gives 2∂∂vlog𝒵ρX(v)|v=s+1/2=∑δ,kℓδe−skℓδsinh⁡(kℓδ/2)tr⁡ρ(δk).2\left.\frac{\partial}{\partial v} \log\mathcal Z_\rho^X(v)\right|_{v=s+1/2} =\sum_{\delta,k} \frac{\ell_\delta e^{-sk\ell_\delta}} {\sinh(k\ell_\delta/2)}\mathop{\mathrm{tr}}\rho(\delta^k).(23.31) With coefficients defined by finite length cutoffs, the TCT_C coefficient on the right is DCX(s)D_C^X(s) for ℜs>1\Re s>1. The coefficient estimate proves this assertion; it is not an appeal to continuity of coefficient extraction in the uniform norm. No meromorphic continuation of the individual DCD_C is established by this calculation.

Removing all iterates by a global Chebyshev substitution is not legitimate in the rank-two character ring. If x=tr⁡Ax=\mathop{\mathrm{tr}}A, y=tr⁡By=\mathop{\mathrm{tr}}B, u=tr⁡(AB)u=\mathop{\mathrm{tr}}(AB), and v=tr⁡(AB−1)v=\mathop{\mathrm{tr}}(AB^{-1}), then u+v=xyu+v=xy, but P2(u)+P2(v)−P2(x)P2(y)=2(tr⁡[A,B]−2),P_2(u)+P_2(v)-P_2(x)P_2(y) =2\bigl(\mathop{\mathrm{tr}}[A,B]-2\bigr),(23.32) which is generally nonzero. Thus sending every trace to its Chebyshev polynomial does not define a multiplicative operation on this character ring.

There is an exact repair if extra twisting data are supplied. For the original SU(2)\mathrm{SU}(2) representation ρ\rho, work in the virtual representation ring generated by its symmetric powers. Let ψkρ\psi^k\rho denote its Adams operation, characterized by χψkρ(γ)=χρ(γk)\chi_{\psi^k\rho}(\gamma)=\chi_\rho(\gamma^k). For the standard rank-two representation, ψkρ=Sym⁡kρ−Sym⁡k−2ρ(k≥2),ψ1ρ=ρ.\psi^k\rho=\operatorname{Sym}^k\rho- \operatorname{Sym}^{k-2}\rho\quad(k\geq2), \qquad \psi^1\rho=\rho. Trace formulas extend additively to virtual representations. All twists used here have inverse-invariant characters, so the unoriented normalization is preserved. Let Σρ\Sigma_\rho denote the character-valued geometric measure on the right of (23.4), and set Qρ=sinh⁡(L/2)Σρ,Pρ=∑δ∈𝒫unℓδχρ(δ)δℓδ.Q_\rho=\sinh(L/2)\Sigma_\rho,\qquad P_\rho=\sum_{\delta\in\mathcal P_{\mathrm{un}}} \ell_\delta\chi_\rho(\delta)\delta_{\ell_\delta}. For DkδL=δkLD_k\delta_L=\delta_{kL}, Qρ=∑k≥1DkPψkρ,Pρ=∑k≥1μ(k)DkQψkρ.Q_\rho=\sum_{k\geq1}D_kP_{\psi^k\rho},\qquad \boxed{P_\rho=\sum_{k\geq1}\mu(k)D_kQ_{\psi^k\rho}.} These sums are locally finite, and the second identity follows from ψjψk=ψjk\psi^j\psi^k=\psi^{jk} and divisor Möbius inversion. It requires spectra for higher symmetric-power twists, which are additional data beyond the original rank-two family. The explicit subtraction in (23.21) needs no such enlargement.

23.10 Scope and further questions

The construction now supplies signed coefficient distributions, an isolated first simple-curve atom, globally convergent heat transforms, a Laplace half-plane, explicit iterate families, a figure-eight cancellation, and metric derivatives. It has not produced a positive measure with the weights of simple-geodesic counting, nor a pointwise constant McShane series.

The scalar limitation is genuine: length-isospectral closed hyperbolic surfaces can have different simple length spectra [25]. The labelled holonomy family retains the information used by (23.7).

For a counting comparison, let d=6g−6d=6g-6 and normalize Thurston measure by the integral multicurve lattice. If B(X)=μTh{λ:ℓX(λ)≤1},B(X)=\mu_{\mathrm{Th}}\{\lambda:\ell_X(\lambda)\leq1\}, ordinary lattice counting gives #{C∈ℳℒℤ:ℓX(C)≤R}∼B(X)Rd\#\{C\in\mathcal{ML}_{\mathbb Z}:\ell_X(C)\leq R\}\sim B(X)R^d [77, Proposition 3.1]. Its Laplace consequence is ∑C∈ℳℒℤe−tℓX(C)∼Γ(d+1)B(X)t−d(t↓0).\sum_{C\in\mathcal{ML}_{\mathbb Z}}e^{-t\ell_X(C)} \sim\Gamma(d+1)B(X)t^{-d}\qquad(t\downarrow0). This is a target for a measure comparison, not a consequence of the coefficient heat expansion. It counts all weighted multicurves; connected or fixed-orbit counts require additional information. The large-length limit here is also different from the minimum-selection limits in (23.18).

Two concrete further problems are to continue DC(s)D_C(s) beyond its proved half-plane, with control of the coefficient operation, and to find combinations of coefficients whose cancellations organize globally over pants. Formula (23.26) provides an exact local model for the second problem. An identity constant in XX would require an additional normalization or geometric partition; the existence of these transforms alone does not supply that constancy. Section 24 strengthens the lattice benchmark and gives explicit orbital gaps, uniform recovery estimates, and a conditional counting criterion.