==== TABLE OF CONTENTS ==== Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//asymptotics.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//chords.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//common-twist.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//counting.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//cylinder.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//higher-jets.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//introduction.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//local.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//main.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//physical-amplitude.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//projection.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//references.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//repeated-moments.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//residue.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//seam.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//site-moments.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//slit.tex Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//vacuum.tex =========================== ==== BEGIN Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//asymptotics.tex ==== \section{Asymptotics of the cap coefficient}\label{sec:asymptotics} The physical coefficient in \eqref{eq:physical-coefficient} is a normalized ratio \(W/J\). We compute its logarithmic growth. The product supplies the twist dependence; the leading logarithmic contribution of the determinant will be the same for every twist. At the homogeneous consecutive roots \(t_j=iy_j\), set \[ L_a(y)=\log(\sinh^2y+\sin^2(al)),\qquad g=L_0+L_2-L_1, \qquad a_i=N\rho(y_i). \] We use \(x,T,E,r,\rho\) from Section~\ref{sec:counting}. Up to additive constants independent of \(\delta\), the product and Jacobian formulas give \[ \begin{aligned} \log|W|&=\sum_{i\eta N\) has total mass \(O_\eta(1)\), because \(w_i\le C/(\eta N)\) there. On the left remaining range \(M\le i\le\eta N\), the boxes give \(x_i=i+O(1)\) uniformly in \(\delta\). Given \(\varepsilon>0\), choose \(\eta\) small enough that the left-tail ratio in \eqref{eq:quantile-density} is within \(\varepsilon\) of one whenever \(u/N\le2\eta\), and choose \(M\) large enough that \(x_i/i\) is within \(\varepsilon\) of one for \(i\ge M\). For large \(N\), throughout this range, \[ w_i=\frac{1+O(\varepsilon)}{\alpha_{\rm r}i} \] uniformly in \(\delta\). The right remaining range has the same estimate with \(N+1-i\). Since \(\sum_{M\le i\le\eta N}i^{-1}=\log N+O_{M,\eta}(1)\), first \(N\to\infty\) and then \(\varepsilon\to0\) prove \eqref{eq:harmonic-mass} uniformly. We will use the following consequence. Suppose arrays \(z_{N,i}(\delta)\) and a function \(z_*(\delta)\) are uniformly bounded and satisfy \begin{equation}\label{eq:uniform-bulk-condition} \lim_{M\to\infty}\limsup_{N\to\infty} \sup_{\substack{|\delta|\le1/4\\ \kappa_i\ge M}} |z_{N,i}(\delta)-z_*(\delta)|=0. \end{equation} Then \begin{equation}\label{eq:harmonic-averaging} \sup_{|\delta|\le1/4} \left|\frac1{\log N}\sum_iw_i z_{N,i}(\delta) -\frac2{\alpha_{\rm r}}z_*(\delta)\right|\longrightarrow0. \end{equation} To prove it, choose \(M\) so the error in \eqref{eq:uniform-bulk-condition} is small for large \(N\). Its weighted sum outside the fixed edge layers is small times \(\sum_iw_i\); inside those layers it is \(O_M(1)\). Divide by \(\log N\) and use \eqref{eq:harmonic-mass}. This proves the assertion with the supremum in \(\delta\). A scalar bulk limit proved along every sequence with \(N\to\infty\) and \(\kappa_i\to\infty\), allowing arbitrary twists in the compact interval, implies \eqref{eq:uniform-bulk-condition} by contradiction and subsequence extraction. Apply this rule to \(z_{N,i}=h_iB_i\). It is uniformly bounded by the inverse estimate and \eqref{eq:uniform-force-bound}; by \eqref{eq:bulk-counting-limits} and \eqref{eq:bulk-force}, its bulk limit is \(-4s\,ds/d\delta\). Therefore \eqref{eq:product-force} gives \begin{equation}\label{eq:uniform-product-derivative} \frac{\Phi_N'(\delta)}{\log N} =-\frac8{\alpha_{\rm r}}s\,\frac{ds}{d\delta}+o(1) \end{equation} uniformly for \(|\delta|\le1/4\). \subsection{The determinant contribution} It remains to show that the leading logarithmic contribution of \(\det H\) is independent of \(\delta\). Factor \[ H=\operatorname{diag}(d_i)(I-A),\qquad d_i=\frac{f_i}{a_i},\qquad A_{ij}=\frac1{d_i}\frac{r(y_i-y_j)}{a_j}. \] The strict margin \eqref{eq:counting-dominance} and the uniform upper bound on \(d_i\) imply \[ d_i\ge\epsilon,\qquad \sum_j|A_{ij}|\le q_0<1 \] with absolute \(\epsilon,q_0\), uniformly throughout the rectangles. The diagonal logarithms are negligible at the scale \(\log N\). In fact \(f_i-a_i=(r'*E)(y_i)\) is uniformly bounded, and the logarithm is Lipschitz on the fixed interval containing all \(d_i\). Hence \[ |a_i\log d_i|\le C|(r'*E)(y_i)|. \] The right side is uniformly bounded and tends to zero in every bulk limit by \eqref{eq:bulk-smooth-correction}. The averaging rule \eqref{eq:harmonic-averaging}, applied to \(z_{N,i}=a_i\log d_i\), gives \begin{equation}\label{eq:diagonal-log-negligible} \frac1{\log N}\sum_i\log d_i=o(1) \end{equation} uniformly in \(\delta\). The path \(I-tA\), \(0\le t\le1\), is invertible because \(\|tA\|_\infty<1\). Its real determinant remains positive, and the convergent trace series is \[ \log\det(I-A)=-\sum_{k\ge1}\frac{\operatorname{tr}A^k}{k}. \] For all \(N,i,\delta\) and \(k\ge1\), its diagonal terms have the summable bound \begin{equation}\label{eq:trace-term-bound} a_i|(A^k)_{ii}|\le Cq_0^{\,k-1}. \end{equation} Indeed \(a_i|A_{ji}|=|r(y_j-y_i)|/d_j\le C\), while \(\sum_j|(A^{k-1})_{ij}|\le q_0^{k-1}\); multiply these two bounds in the final factor of \(A^k\). For fixed \(k\), the bulk limit of the same quantity is \begin{equation}\label{eq:trace-bulk-limit} a_i(A^k)_{ii}\longrightarrow (r^{*k})(0). \end{equation} For \(k=1\) this is \(r(0)/d_i\to r(0)\). For larger \(k\), expand the product as \[ a_i(A^k)_{ii} =\sum_{j_1,\ldots,j_{k-1}} \frac{r(y_i-y_{j_1})r(y_{j_1}-y_{j_2})\cdots r(y_{j_{k-1}}-y_i)} {d_i d_{j_1}\cdots d_{j_{k-1}}\, a_{j_1}\cdots a_{j_{k-1}}}. \] First restrict each rapidity step to a fixed bounded interval. All intermediate atoms then lie in a fixed bulk window, where \(d_j\to1\) uniformly and the weighted quantile grids converge to Lebesgue measure. Iterated quadrature gives the corresponding truncated convolution. The restriction can be removed uniformly: every row of \(|A|\) has negligible mass outside a large rapidity step by \eqref{eq:counting-kernel-tails}, and the terminal factor \(a_i|A_{ji}|=|r(y_j-y_i)|/d_j\) is bounded and decays exponentially with its step. For fixed \(k\), these bounds control the union of the finitely many omitted step ranges. This proves \eqref{eq:trace-bulk-limit} along every bulk sequence. Apply \eqref{eq:harmonic-averaging} for each fixed \(k\) to \(z_{N,i}=a_i(A^k)_{ii}\). Equations \eqref{eq:trace-term-bound}--\eqref{eq:trace-bulk-limit} give \[ \frac{\operatorname{tr}A^k}{\log N} =\frac2{\alpha_{\rm r}}(r^{*k})(0)+o(1) \] uniformly in \(\delta\). Passing this limit through the trace series is justified by the same bound: \[ \frac{|\operatorname{tr}A^k|}{\log N} \le Cq_0^{k-1}\frac{\sum_iw_i}{\log N} \le C' q_0^{k-1} \] for all large \(N\), uniformly in \(\delta\). After division by \(k\), the right side is summable. Thus \[ \frac{\log\det H}{\log N} =-\frac2{\alpha_{\rm r}}\sum_{k\ge1}\frac{(r^{*k})(0)}k+o(1) \] uniformly in \(\delta\), where \eqref{eq:diagonal-log-negligible} was used. The series limit is independent of the twist. Integrate \eqref{eq:uniform-product-derivative} from the vacuum value \(\delta=1/4\) to \(\delta=\theta/(2\pi)\), and subtract the determinant contributions at those two values. Uniformity of both estimates gives \[ \frac{\log\left|(W/J)_\theta/(W/J)_{\pi/2}\right|}{\log N} =-\frac4{\alpha_{\rm r}} \bigl(s_\theta^2-s_{\pi/2}^2\bigr)+o(1), \qquad 0\le\theta\le\frac\pi2. \] By \eqref{eq:physical-coefficient} this is the logarithmic asymptotic of \(c_{\rm phys}(\theta)\). The bulk row has spin twist \(K^{-1}=e^{-i\theta}\) and winding fugacity \(K+K^{-1}=2\cos\theta\), while the two vacuum caps remain fixed. Since \(s_\theta=2\theta/(3\pi)\) and \(\alpha_{\rm r}=8/3\), we obtain \begin{equation}\label{eq:amplitude-exponent} \frac{\log|c_{\rm phys}(\theta)|}{\log N} =\frac16-\frac{2\theta^2}{3\pi^2}+o(1), \qquad 0\le\theta\le\frac\pi2, \end{equation} with the error uniform on this interval. At zero winding fugacity, \(\theta=\pi/2\), the coefficient is exactly one. At physical fugacity \(2\), \(\theta=0\), it is \(N^{1/6+o(1)}\). The next section uses the geometric interpretation and height estimates to turn this coefficient into a planar nesting law. ==== END Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//asymptotics.tex ==== ==== BEGIN Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//chords.tex ==== \section{From polygon nests to triangle chords}\label{sec:chords} The planar nesting estimate is now proved. We finish by using positive flux identities to normalize the macroscopic triangle chords and calculate their total marked mass. We first record a separate consequence of the positive estimates: the weight $x$ is the critical step fugacity. \subsection{An independent criticality consequence} Let $c_n$ be the number of ordinary centre-to-centre self-avoiding walks of $n$ edges from a fixed centre. We show that the series $\sum_n c_n z^n$ has radius of convergence $x$. The decomposition is the decreasing-span unfolding method associated with Hammersley and Welsh~\cite{HammersleyWelsh1962}, in the present port convention. Along the normal to a horizontal row, centre heights lie on a uniform discrete mesh, and no dual bond is horizontal. Split a walk at its global minimum height and reverse the first part. Each of the two resulting walks is nonnegative relative to its starting height. In either part, take the segment ending at the last vertex of maximum displacement from its start, and repeat on the remainder with the sign reversed. Every segment has positive span and has its two ends at its two extreme heights. The spans strictly decrease after the first comparison, which may be an equality. Their sum is $O(n)$, since the total variation of height along an $n$-edge walk is $O(n)$. There are consequently $O(\sqrt n)$ segments. For a prescribed start, sign, and span, the total weight of such a segment is bounded by an absolute constant. To see this, add boundary legs at its two extremes to make it a crossing between parallel horizontal triangle-side lines. If the outward bond at an extreme is vertical, use its midpoint port. Otherwise first use a fixed slanted outward bond to reach a new triangle with a vertical outward bond. The additional centre lies strictly beyond the old extreme, so it is unvisited. The two boundaries and the source port are determined by the start, sign, and span up to finitely many lattice types, and deleting the added legs has bounded multiplicity and bounded weight cost. Exhausting the strip by finite convex parallelograms, the positive real flux of Section~\ref{sec:cylinder} bounds the total crossing mass by an absolute constant. The conversion between $x$ per edge and $x$ per visited centre costs one further constant for each segment. At a fixed starting centre, the minimum centre has only $O(n^2)$ possible positions. There are at most $n^{O(\sqrt n)}$ possible span lists, and the product of the crossing bounds costs $C^{O(\sqrt n)}$. Dropping endpoint and total-length constraints within the segment sums therefore gives \[ c_nx^n\le \exp\bigl(O(\sqrt n\log(n+2))\bigr). \] The radius of convergence is at least $x$. For the opposite inequality, half-plane saturation and the definition of $D_N$ give the positive omitted mass \[ D_N=2\cos(3l)\left( \sum_{g\ge N}K_\infty(g) +\sum_{g=1}^{N-1}\bigl(K_\infty(g)-K_N(g)\bigr)\right) \asymp N^{-1/4}. \] Each omitted arc uses at least $a_{\rm len}N$ centres for an absolute $a_{\rm len}>0$. An arc in the first sum has endpoints separated by at least $N$; one in the second meets a nonzero horizontal translate and hence has two vertices separated by at least $N$. Deleting the port legs maps these arcs to centre walks with bounded multiplicity, so for large $N$, \[ D_N\le C\sum_{n\ge a_{\rm len}N/2}c_nx^n. \] If the counting series had radius larger than $x$, choosing a fugacity strictly between them would make this tail exponentially small in $N$. This contradicts the polynomial lower bound. Thus the radius is exactly $x$, at either centre type by lattice symmetry. The denominator proof that follows uses the positive deficit and depth estimates directly. \subsection{The macroscopic chord denominator} Put $\chi=\cos(3l)>0$ and write \[ B_\infty=2\sum_{g>0}K_\infty(g)=\chi^{-1} \] for the total half-plane arch mass from one boundary port, including both endpoint directions. We will compare this saturated mass with the arches that fit in a strip or a triangle. \paragraph{Strip crossings.} Let $B_d$ be the same-side arch mass in the infinite strip of width $d$ rows from a fixed source port, and let $S(d)$ be its mass of crossings to the opposite parallel boundary. In a finite convex parallelogram $\Pi$ exhausting this strip, the same-side exits have cosine $\chi$, the opposite-side exits have cosine one, and the real flux through the two end sides is nonnegative. With the analogous finite masses this is \[ 1=\chi B_\Pi+S_\Pi+E_\Pi^{\rm end}, \qquad E_\Pi^{\rm end}\ge0. \] Monotone exhaustion of the two indicated path classes gives the strip missing-mass inequality \begin{equation}\label{eq:strip-missing-mass} S(d)\le 1-\chi B_d=\chi(B_\infty-B_d). \end{equation} The strip contains every quotient-self-avoiding arc of positive span $g0$ so that $(d-C_{\rm lat})/N-C_*>1$ for all large $d$. The depth sum is then $O(N^{-1/4})$ by \eqref{eq:arc-depth-moments}, already with $k=1$. Proposition~\ref{prop:cylinder} bounds $D_N$ by the same order. Increasing the constant for bounded $d$ proves \begin{equation}\label{eq:strip-crossing-bound} S(d)\le C d^{-1/4}. \end{equation} The bound holds for each of the finitely many source types and their lattice rotations. \paragraph{Triangle flux and a parallel cut.} For a source port $a$ on one side of an equilateral lattice triangle $\Delta$, let $B_\Delta(a)$ be the unsigned mass of paths returning to that side, and let $E_\Delta(a)$ be the unsigned mass of paths exiting on either of the other sides. A same-side chord has total turn $\pm\pi$, while a chord to another side has total turn $\pm\pi/3$. The finite real flux identity is therefore \begin{equation}\label{eq:triangle-flux} \chi B_\Delta(a)+\cos(l)E_\Delta(a)=1,\qquad \cos(l)E_\Delta(a)=\chi\bigl(B_\infty-B_\Delta(a)\bigr). \end{equation} This identifies the other-side exits with the half-plane arch mass lost by restricting the paths to the triangle. Clip $\Delta$ by a parallel lattice cut at inward depth $d$ rows from the source side, and retain the component $\Delta_{\le d}$ containing that side. Let $\mathcal L_d(a)$ be the paths from $a$ to the original boundary that use a removed triangle, and let $\mathcal Q_d(a)$ be the paths in $\Delta_{\le d}$ from $a$ to the new cut. Subtract the real flux identities before and after clipping. Every path that remains in the clipped domain and exits on an original side cancels. The new cut is parallel to the source side, so its exits have total turn zero. The exact remaining identity and its consequence are \begin{equation}\label{eq:triangle-cut-flux} \begin{gathered} \sum_{\gamma\in\mathcal L_d(a)}w(\gamma)\cos(cW(\gamma)) =\sum_{\gamma\in\mathcal Q_d(a)}w(\gamma)\le S(d),\\ \sum_{\gamma\in\mathcal L_d(a)}w(\gamma) \le \chi^{-1}S(d)\le C d^{-1/4}. \end{gathered} \end{equation} The last inequality uses $\cos(cW)\ge\chi$ on every convex triangle chord. In particular it bounds all exits whose paths reach beyond the parallel cut. \paragraph{The half-plane diameter tail.} Let $\mathcal H(r)$ be the total mass of half-plane arches from a fixed port, in both directions, with diameter at least $r$. For large $r$, inscribe an equilateral lattice triangle of side $\lfloor r/2\rfloor$ in the half-plane with its base on the boundary. From any port $a$ in the central third of the base, every such arch is missing from the triangle, whose diameter is less than $r$. Equation \eqref{eq:triangle-flux} gives \[ \chi\mathcal H(r)\le \cos(l)E_\Delta(a). \] Sum over these order-$r$ base ports and reverse the paths on the right. From each source on a slanted side, the reversed paths reach a central base port, whose distance from that side is at least a fixed positive multiple of $r$. They therefore reach beyond a parallel cut of depth $c_2r$ rows, with an absolute $c_2>0$ and harmless integer rounding. Equation \eqref{eq:triangle-cut-flux}, summed over the $O(r)$ slanted ports, yields \[ c_3r\,\mathcal H(r)\le C r S(c_2r)\le C r^{3/4}. \] Translation invariance makes $\mathcal H(r)$ the same for every central base source. We have proved \begin{equation}\label{eq:halfplane-diameter-tail} \mathcal H(r)\le C r^{-1/4}. \end{equation} \paragraph{Integrating both endpoints.} For a fixed $\delta>0$, let $\mathcal C_R(\delta)$ be the ordered boundary-to-boundary paths in $\Delta_R$ of diameter at least $\delta R$, and put $Z_R(\delta)=\sum_{\gamma\in\mathcal C_R(\delta)}w(\gamma)$. We claim \begin{equation}\label{eq:chord-cutoff-denominator} Z_R(\delta)\asymp_\delta R^{3/4} \qquad\text{for every fixed }0<\delta\le1/50. \end{equation} The comparison constants may depend on the fixed cutoff $\delta$. For the upper bound, same-side paths from each port have mass at most $\mathcal H(\delta R)=O_\delta(R^{-1/4})$. For paths to another side, let $r$ be the source port's distance along its side to the nearer corner. Every half-plane arch missing from the triangle has reached a slanted side, at distance at least $c_4r-C_{\rm lat}$ from the source. It thus has diameter of that order. By \eqref{eq:triangle-flux} and \eqref{eq:halfplane-diameter-tail}, all other-side exits from that port have mass at most $Cr^{-1/4}$ when $r$ exceeds a fixed constant. The bounded number of ports closer to corners have bounded mass by finite flux. Summing over the $O(R)$ sources gives \[ Z_R(\delta)\le C_\delta R^{3/4} +C\sum_{r=1}^{R}r^{-1/4} \le C_\delta R^{3/4}. \] For the lower bound take a source in the central third of the base and put $N=3R$. Every triangle base-to-base arch projects self-avoidingly to this quotient: the triangle has horizontal width $R0$, equality $B=M$ occurs only for a unit singleton or a unit adjacent pair. Here is a proof. If some count is zero, the variation estimate in the residue proof gives $\Delta=\frac12\sum(m_k-m_{k+1})^2\ge M$. Equality for $B=M$ also forces all distance-two products to vanish. The contiguous, unit-or-flat rise and fall in the equality case for $\Delta$ then leave only the two stated supports. If every count is positive, the two distance-two products incident with a peak already sum to at least $2M$, so equality $B=M$ is impossible. For the second and third bounds, the cases $M=0,1$ follow from $B\ge M$. Suppose some count is zero and choose a peak at index $0$. If $m_{-2}$ or $m_2$ is positive, its product with the peak adds at least $M$ to $\Delta\ge M$. Otherwise the four gradients from the peak to these two zero counts give \[ \Delta\ge\frac12\bigl((M-m_{-1})^2+m_{-1}^2 +(M-m_1)^2+m_1^2\bigr)\ge\frac{M^2}{2}. \] For $M\ge4$ this is at least $2M$. For $M=3$, equality in the integral lower bound $5$ forces both first neighbors to be positive, so their distance-two product adds at least one; otherwise the gradient contribution is already larger. For $M=2$, the only way those four gradients contribute just $2$ is $m_{-1}=m_1=1$, again adding their product; otherwise they contribute at least $3$. This proves \eqref{eq:common-twist-B-bounds}. \paragraph{Nonnegative depth bounds and their outer equality cases.} For an unanchored marked depth, combining \eqref{eq:common-twist-bare-orders} and \eqref{eq:common-twist-mark-orders} gives $B-m_k+m_{k\pm2}\ge0$. The unmarked bound is $B-1\ge0$. Since an unanchored cluster has at least two labels, equality in either bound requires a unit adjacent pair. For an anchored unmarked depth, $B-a-b'\ge0$ by the stage-$0$ bound in the residue proof. Apart from the matching first kernel, discarding the nonnegative numerator orders gives the marked lower bound \begin{equation}\label{eq:common-twist-anchored-bound} B-a-b'-m_k+1. \end{equation} If at most one of $a,b'$ is positive, this is nonnegative by $a+b'\le M$, $m_k\le M$ and $B\ge2M-1$. If both are positive, \[ B\ge M+ab',\qquad B-a-b'-m_k+1\ge (M-m_k)+(a-1)(b'-1)\ge0. \] The first inequality follows from $\Delta\ge M$ and the distance-two product $ab'$ when some count is zero. When all counts are positive it is strict, because another product incident with a peak supplies $M$ in addition to $ab'$ and further products remain. Equality requires no distance-two products other than $ab'$ and a contiguous monotone support with steps of size at most one. Hence the support is exactly $\{-1,0,1\}$, with $m_{-1}=m_1=1$ and $1\le m_0\le2$. At an outer anchored cluster, \eqref{eq:common-twist-outer-neutrality} rules out exactly one of $a,b'$ being positive: all nested counts at the absent end would still be zero, so the weighted sum would have the strict sign of the other end. If both vanish, the unmarked and marked bounds are strictly positive by $B\ge M$. Thus any outer equality in the bounds just considered has support $\{-1,0,1\}$ with unit end counts. The unmarked equality is the subcase $m_0=1$. \paragraph{The matching first kernel.} It remains to handle an evaluation mark with $s=1$ at its own site and offset $-1$. The anchored clusters inside an outer anchored cluster form a chain, because they all contain the site. Along the initial part of that chain containing the marked label, the exact bound is \[ (B-a-b')+(b'-a+1)-1=B-2a. \] After the marked label leaves the anchored chain, the bound is $B-a-b'$. If $b'>0$ at every marked depth, then $B\ge M+ab'\ge2a$ makes all these bounds nonnegative; equality at the outer depth again has the support and unit ends above. Otherwise some marked anchored depth has $b'=0$. Every deeper anchored depth then has $b'=0$ by nesting, including the depths after the mark has left. Add the zero quantity in \eqref{eq:common-twist-outer-neutrality} to the total anchored order for this outer cluster. The bounds become $B-a-b'$ while the mark is anchored and $B-2b'=B$ afterwards. They are all nonnegative. In a total equality, the outer bound is $B-a-b'=0$; outer neutrality and the marked label at $-1$ again force the same support with unit ends. \paragraph{Strict order on a neutral ray.} Suppose the total logarithmic order of one of the three forms were zero. The preceding nonnegative bounds would have to be equalities at every positive depth. If an outer cluster were unanchored, its unit adjacent pair would have pair charges $+1,-1$ and no strict subcluster, contradicting \eqref{eq:common-twist-neutrality}. Every remaining outer equality has support $\{-1,0,1\}$, unit counts at $-1,1$, and positive $m_0$. For its unique label at $-1$, the charge from this outer depth is \[ \underbrace{m_{-3}+m_0-m_{-2}-m_1}_{\text{pair charge}} +\underbrace{1}_{\text{site charge}}=m_0>0. \] At every deeper anchored depth containing that label the charge is $1+m_0-m_1\ge0$, since the count at $1$ is at most one. A deeper unanchored equality containing it must be an adjacent unit pair within the same support, necessarily using an offset-$0$ label, and contributes $+1$. Thus the total $e_i$ is positive, again contradicting neutrality. Every boundary ray meeting $\pi^{-1}(U)$ therefore has strictly positive logarithmic order for each form. Since these orders are integral, the pulled-back forms are holomorphic there, including at intersections of boundary divisors. \subsection{From holomorphic forms to scalar traces} It remains to justify that dividing by the Jacobian and summing over the fibers preserves holomorphy. The restriction $\pi:\pi^{-1}(U)\to U$ is proper. Let $\eta$ be any of the holomorphic top forms just obtained, and put $\Omega=\bigwedge_i d\log K_i$, a nonvanishing volume form on $U$. On regular fibers define \[ T_\eta(\mathbf K)= \sum_{x\in\pi^{-1}(\mathbf K)} \frac{\eta}{\pi^*\Omega}(x). \] For $\eta=\eta_f$, the denominator in root coordinates is $J$, so this is exactly $\tau(f)$. First, $T_\eta$ is a rational function on the generic target. If the target coordinates $-D_i$ are algebraically dependent, the generic fiber is empty and this trace is identically zero, so there is nothing further to prove. Otherwise $\mathbb C(p_1,\ldots,p_N)$ is a finite algebraic extension of $\mathbb C(-D_1,\ldots,-D_N)$, by equality of transcendence degrees. It is separable in characteristic zero. Adjoin the $p_i$ successively: outside the zeros of the minimal-polynomial discriminants and the relevant denominators, their distinct roots give precisely the generic inverse branches. A nonzero input denominator can be excluded by the nonzero norm of that field element on the target. Summing the rational function $\eta/\pi^*\Omega$ over the branches is symmetric at each field step, and is therefore rational. This also gives a finite generic degree $d_\pi$ bounding the number of branches. Now take a relatively compact coordinate polydisc $V\Subset U$. On its regular values, Cauchy--Schwarz and change of variables give \begin{equation}\label{eq:common-twist-L2} \int_{V_{\rm reg}}|T_\eta|^2\,dV_\Omega \ \le\ d_\pi\int_{\pi^{-1}(V)}dV_\eta \ <\ \infty, \end{equation} where $dV_\Omega=i^{N^2}\Omega\wedge\overline\Omega$ and $dV_\eta=i^{N^2}\eta\wedge\overline\eta$. Indeed $|\sum_{\ell=1}^{d_\pi}a_\ell|^2\le d_\pi\sum_\ell|a_\ell|^2$ for the branch values $a_\ell=\eta/\pi^*\Omega$, and the squared complex Jacobian in change of variables cancels their denominators. The last integral is finite because $\eta$ is holomorphic on the compact preimage of $\overline V$. Positive-dimensional exceptional fibers cause no problem, since the estimate uses only the regular fibers and the generic degree. A rational function with a pole along a divisor is not locally square-integrable: at a generic smooth point of that divisor its form is $z_1^{-m}$ times a holomorphic unit, $m\ge1$, and the integral in the transverse complex coordinate diverges. Thus \eqref{eq:common-twist-L2} excludes every polar divisor of $T_\eta$; cancellation of the irreducible factors of a local denominator then shows that it is holomorphic on $U$. The same argument works jointly with sites in a small neighborhood avoiding one another's shift orbits. The compactification is then made relatively in the site coordinates; wedging $\eta$ with their differentials and using the proper map $(\mathbf p,\mathbf b)\mapsto(-\mathbf D,\mathbf b)$ gives equal dimensions and the same estimate. In the punctured family just considered, these local proper-trace arguments cover the common product neighborhood $U\times\{0<|\epsilon|<\delta\}$. Consequently, the trace of a Laurent insertion whose coefficients are jointly holomorphic in the sites, twists and any row parameters is a finite holomorphic combination of the rational monomial traces, and is holomorphic there. When those coefficients are also rational in the parameter coordinates, the trace is rational and has no polar divisor in that product neighborhood. Any inverted scalar row denominators must stay regular and nonzero. This supplies the off-$\epsilon=0$ input to the later two-site proof, without asserting an order at the collision itself. Apply this result to the three families of forms. It gives the holomorphic continuation of $\tau(f)$. Telescoping a product of site evaluations expresses its difference from the corresponding product of $P(b)$ as a sum of the holomorphic marked traces multiplied by $K-K_i$ or $K^{-1}-K_i^{-1}$. Their restrictions to the diagonal vanish. After clearing the spectral denominators, the same argument for every Laurent-polynomial multiple of a coefficient gives its annihilation by the specialized trace. This proves all three assertions. \end{proof} \subsection{A finite-jet representation of the specialized trace} The collision limit need not be a sum of ordinary evaluations. The next corollary records the additional derivatives that can remain, and the relations satisfied at every point that contributes. It is used at generic $K$ satisfying the preceding hypotheses. Section~\ref{sec:physical-amplitude} reaches the physical values $K=i$ and $K=1$ by analytic continuation of the physical branch and its coefficient, after using this corollary at generic nearby twists. Write $\mathbf p=(p_1,\ldots,p_N)$, where $p_i=e^{2it_i}$, and let \[ \mathcal A=\mathbb C[p_1^{\pm1},\ldots,p_N^{\pm1}] \] be the Laurent polynomials regular on the entire root torus. \begin{corollary}[Finite jets at a common twist]\label{cor:common-twist-jets} Fix the sites and a generic common twist as in Proposition~\ref{prop:common-twist}, and denote the specialized trace by $\tau_0$. There are finitely many points $\mathbf p^{(\nu)}\in(\mathbb C^*)^N$, integers $m_\nu\ge0$, and constants $c_{\nu,\alpha}$ such that, for every $f\in\mathcal A$, \begin{equation}\label{eq:common-twist-jets} \tau_0(f)=\sum_\nu\ \sum_{|\alpha|\le m_\nu} c_{\nu,\alpha}\,\partial_{\mathbf p}^{\alpha}f(\mathbf p^{(\nu)}). \end{equation} Here $\alpha\in\mathbb Z_{\ge0}^N$, $|\alpha|=\sum_i\alpha_i$, and repeated points are combined and zero local functionals omitted. Every remaining point satisfies the coefficients of the cleared relation $Q(w)Q(w-1)(L-P)=0$; consequently $L=P$ there as a rational identity in the spectral argument. Such a point may have colliding roots. At a regular transverse solution in the original configuration locus, the contribution is the ordinary trace term $Wf/J$, with $J$ as in \eqref{eq:common-twist-trace}; its weight is nonzero whenever $W\ne0$. \end{corollary} \begin{proof} For $N=0$, the root torus is the single empty tuple, $\mathcal A=\mathbb C$, and $\tau_0(f)=f$. The representation is one order-zero jet with coefficient one, and the cleared relation vanishes because $L=P$. Assume $N\ge1$ below. If the generic fiber is empty, the trace is zero and the empty sum gives the representation, so we may also assume that the map is dominant. Choose a generic affine line through the fixed diagonal point, with parameter $\epsilon=0$ there. After shrinking the disk, its nonzero points avoid the exceptional algebraic locus of the finite generic fibers. The finitely many inverse branches have finite permutation monodromy, so after a power substitution $\epsilon=s^m$ they are single-valued on the punctured $s$-disk. Properness and the exclusion of toric ends keep every root coordinate bounded above and away from zero. Each branch therefore extends holomorphically to $s=0$, with a limit in $(\mathbb C^*)^N$; collisions are allowed. The branch weights $W/J$ are rational functions of the roots and parameters. On each extended branch they are meromorphic in $s$, with some finite pole order. Choose one integer $M$ bounding those orders for all branches. It is independent of $f$. Since the trace is holomorphic at $\epsilon=0$, $\tau_0(f)$ is the constant term in $s$ of the sum of the weighted branch values $f(\mathbf p(s))$. The Taylor expansion at a branch limit contributes to this constant term only through derivatives of order at most $M$, because $\mathbf p(s)-\mathbf p(0)=O(s)$. This gives \eqref{eq:common-twist-jets}; combine equal limits and discard zero local functionals. Let $c(\mathbf p)$ be any coefficient of the cleared relation. Proposition~\ref{prop:common-twist} gives $\tau_0(cf)=0$ for every $f\in\mathcal A$. The powers of the distinct maximal ideals at the finitely many support points are pairwise comaximal, so Laurent polynomials can prescribe their finite jets independently. If $c(\mathbf p^{(\nu)})\ne0$, multiplication by $c$ is invertible in the finite jet algebra at that point. We could then prescribe an arbitrary jet of $cf$ there and zero jets at the other points, forcing the local functional at $\nu$ to vanish. This contradicts its retention in the sum. Thus every cleared coefficient vanishes at every support point. Since $Q(w)Q(w-1)$ is a nonzero rational function of the spectral argument for torus roots, the cleared identity gives $L=P$ as a rational identity there. Finally, at a regular transverse solution in the original configuration locus the implicit function theorem gives a unique ordinary local branch. Its weight is holomorphic with value $W/J$, so its constant term is precisely $(W/J)f$ at that point, with no derivatives. It is nonzero when $W\ne0$, as asserted. \end{proof} ==== END Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//common-twist.tex ==== ==== BEGIN Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//counting.tex ==== \section{The real counting branch, including its edge roots}\label{sec:counting} The exact moment formula is now available. We locate its physical term by a real root system. The amplitude calculation will require an inverse for the counting Jacobian that is uniform in the number of roots and in the twist, including at the two extreme roots. At \(l=\pi/8\), take the real slice \(t_j=iy_j\), \(1\le j\le N\), and homogeneous sites zero. Define \[ p_k(y)=\frac{\sin(2kl)}{\pi(\cosh2y-\cos2kl)},\qquad p=p_1,\qquad r=p_1-p_2,\qquad R(y)=\int_{-\infty}^y r(v)\,dv, \] and, for an ordered configuration \(y_1<\cdots-0.012,\qquad \int_{\alpha_{\rm r}^{-1}\log(200/27)}^\infty |r(y)|\,dy<0.012. \end{equation} \paragraph{Boundary signs and containment.} The quantile intervals imply \[ \|E+s\|_\infty\le0.71. \] Indeed, between two successive atoms \(E+s=j-x+s\), and the two neighboring interval bounds put this value between \(-0.71\) and \(0.71\); the two exterior intervals obey the same bound. Since \(\delta=(1-\int r)s\), the residual has the useful form \[ F_i=x_i-(i-\tfrac12+s)+(r*(E+s))(y_i). \] On an untrimmed upper face \(x_i=i-\tfrac12+s+0.21\) it is positive, and on the corresponding lower face it is negative, because \[ \bigl|(r*(E+s))(y_i)\bigr| \le0.71\|r\|_1<0.20945<0.21. \] This also gives containment before the two trims. For any ordered solution of the consecutive equations, put \(M=\max_i|x_i-(i-\tfrac12+s)|\). The same interval argument now gives \(\|E+s\|_\infty\le\tfrac12+M\), including the exterior intervals since \(|s|\le1/3\). Therefore \[ M\le0.295(\tfrac12+M),\qquad M\le\frac{0.1475}{0.705}<0.21. \] Only the lower face for the first quantile and the upper face for the last quantile can be trimmed. We prove a slightly stronger statement for the first one: throughout the untrimmed box, \(x_1\le0.09\) forces \(T(y_1)<1/4\) for all sufficiently large \(N\), uniformly in \(s\). Put \(b=1+s\in[2/3,4/3]\). All atoms \(x_j\), \(j\ge2\), lie above \(b\), and \[ \eta_N(u)=\#\{j\ge2:x_j\le u\}-(u-b),\qquad b\le u\le N, \] satisfies \(\eta_N(b)=0\) and \(\|\eta_N\|_\infty\le0.71\). Comparing the atoms with the continuum identity for \(x(y_1)\), and separating the first atom, gives \[ T(y_1)=x_1+R(0)-\int_0^b R(y_1-y_N(u))\,du +\int_b^N R(y_1-y_N(u))\,d\eta_N(u). \] Here \(R(0)=1/8\), since \(r\) is even and has mass \(1/4\). The boundary terms in integration by parts in the last integral vanish: \(\eta_N(b)=0\), and \(R(y_1-y_N(u))\to0\) as \(u\uparrow N\). It follows that \[ T(y_1)\le x_1+\tfrac18+0.012b+ 0.71\int_{y_N(b)-y_1}^\infty |r(v)|\,dv. \] The exponential tail of \(\rho\) gives \[ e^{\alpha_{\rm r}(y_N(b)-y_N(x_1))} =\frac b{x_1}(1+o(1)) \] uniformly for \(00\) and \(N_0\) such that \begin{equation}\label{eq:counting-dominance} d_i-S_i\ge\epsilon \quad\hbox{for every }\boldsymbol x\in\mathcal B_{N,s}, \quad N\ge N_0,\quad |\delta|\le\tfrac14. \end{equation} Since \(r(0)>0\), this is exactly a lower bound for \(H_{ii}-\sum_{j\ne i}|H_{ij}|\). Also \(d_i\) is uniformly bounded above: \(r'*E=r'*(E+s)\) is bounded by \(0.71\|r'\|_1\), while \eqref{eq:quantile-density} and the trims bound \(a_i\) away from zero. We record the uniform control of the kernel sums used here and below. The quantile intervals \([x_j-0.04,x_j+0.04]\) are disjoint and lie in \((0,N)\). On each such interval \(a_N(u)\asymp a_j\), by \eqref{eq:quantile-density}; their rapidity lengths are therefore comparable to \(1/a_j\), and their rapidities stay a bounded distance from \(y_j\). The exponential decay of \(r\) and its derivatives then implies, for each of these kernels \(K\), \begin{equation}\label{eq:counting-kernel-tails} \sup_{\substack{N\ge2,\ |s|\le1/3,\ \boldsymbol x\in\mathcal B_{N,s}\\1\le i\le N}} \sum_{|y_j-y_i|>L}\frac{|K(y_i-y_j)|}{a_j} \longrightarrow0\qquad(L\longrightarrow\infty). \end{equation} For example, compare each summand with the integral of a constant multiple of \(e^{-c|y_i-z|}\) over its rapidity interval; the intervals are disjoint. This also bounds the full sums uniformly. Consider a sequence of rows for which \(\min(x_i,N-x_i)\to\infty\). Then \(a_i\to\infty\), so \(d_i=1+o(1)\). On any fixed rapidity window about \(y_i\), boundedness of \(\rho'/\rho\) gives \(N\rho\asymp a_i\to\infty\) uniformly. In quantile coordinates \(a_N'(u)=\rho'(y_N(u))/\rho(y_N(u))\) is bounded. Thus on the unit quantile cells containing the atoms, \(a_N\) changes by \(O(1)\), its relative change tends to zero, and the rapidity mesh tends to zero. Sampling at \(x_j\), whose offset in each cell is bounded by the box, gives \[ \sum_j\frac{|r(y_i-y_j)|}{a_j}\longrightarrow\int_{\mathbb R}|r(y)|\,dy. \] The passage from a fixed window to the full line follows from \eqref{eq:counting-kernel-tails}. Hence the limiting bulk margin is \(1-\|r\|_1>0.705\). \paragraph{The two limiting edge deficits.} It remains to examine sequences for which \(\min(x_i,N-x_i)\) stays bounded. Reflecting if necessary, take a subsequence with \(x_i\to x\ge0.09\), \(s\to s_*\), and every fixed left quantile \(x_j\to v_j\). The index \(i\) is fixed on this subsequence. Set \[ b=1+s_*\in[2/3,4/3],\qquad Q=\frac bx,\qquad h(u)=r\bigl(\alpha_{\rm r}^{-1}\log u\bigr)\quad(u>0). \] Write \(h_+=\max(h,0)\) and \(h_-=\max(-h,0)\). Let \(\mu=\sum_{j\ge1}\delta_{v_j}\). The discrepancy \[ D(v)=\mu((0,v])-v+s_* \] has absolute value at most \(0.71\). The tail quantile limits and \eqref{eq:counting-kernel-tails} give \begin{equation}\label{eq:limiting-edge-margin} \lim\alpha_{\rm r}(d_i-S_i) =\alpha_{\rm r} +\frac1x\int_0^\infty h(v/x)(\mu(dv)-dv) -\sum_{j\ge1}\frac{|h(v_j/x)|}{v_j}. \end{equation} The signed integral here is a paired difference, defined without subtracting two divergent integrals by \[ \int_0^\infty h(v/x)(\mu(dv)-dv) :=-\int_0^\infty D(v)\frac1x h'(v/x)\,dv. \] It converges because \(D\) is bounded and \(h\) has finite total variation. To justify \eqref{eq:limiting-edge-margin}, first restrict to bounded rapidity offsets, where \(a_N(x_i)\to\alpha_{\rm r}x\) and \(y_N(v)-y_N(x_i)\to\alpha_{\rm r}^{-1}\log(v/x)\). Then remove the restriction using \eqref{eq:counting-kernel-tails} for the absolute sum and the bound on \(E+s\) against the exponential tails of \(r'\) for the signed correction. Write \(\operatorname{TV}_I g\) for the total variation of \(g\) on \(I\). For the atoms \(j\ge2\), the discrepancy \[ D_b(v)=\#\{j\ge2:v_j\le v\}-(v-b),\qquad v\ge b, \] also has absolute value at most \(0.71\) and vanishes at \(b\). Integration by parts against \(D\) or \(D_b\) bounds a signed sum by \(0.71\) times the variation of its kernel. For a nonfirst atom, the exact range inherited from its quantile interval is \[ i\ge2,\qquad b+i-1.71\le x\le b+i-1.29,\qquad x\ge b+0.29\ge\frac{287}{300},\qquad 00\). Pairing the remaining signed atoms with Lebesgue measure only above \(b\) leaves the ordinary integral from \(0\) to \(b\). The deficit is therefore at most \begin{equation}\label{eq:first-edge-deficit} \begin{aligned} \mathcal D_1(Q,b)={}& \int_0^Q h(u)\,du +\frac{0.71}{b}Q\operatorname{TV}_{[Q,\infty)}h +\int_Q^\infty |h(u)|\,\frac{du}{u}\\ &\quad +\frac{0.71}{b}Q\operatorname{TV}_{[Q,\infty)}(|h|/u). \end{aligned} \end{equation} These are the two estimates that must be kept strictly below \(8/3\). \paragraph{Elementary bounds for the two deficits.} Put \(a=1/\sqrt2\), \(c=1-a\), \(z=u^{3/4}\), and \(v=(z+z^{-1})/2=\cosh(\tfrac34\log u)\). Direct substitution gives \begin{equation}\label{eq:edge-kernel-formulas} \begin{aligned} h(u)&=\frac1\pi\left(\frac a{v-a}-\frac1v\right) =\frac{a-cv}{\pi v(v-a)},\\ \mathcal P(u):=\int_0^u h(w)\,\frac{dw}{w} &=\frac8{3\pi}\left(\frac\pi4+ \arctan(z/a-1)-\arctan z\right). \end{aligned} \end{equation} In particular \(h(u)=h(1/u)\). Its positive interval is \((\zeta^{-1},\zeta)\), where \[ v_0=\frac a{1-a},\qquad \zeta=(v_0+\sqrt{v_0^2-1})^{4/3},\qquad 7.5<\zeta<8. \] The maximum is \(h(1)=\sqrt2/\pi<0.451\). The two negative minima occur at \(v=a/(1-\sqrt a)\), with common absolute value \[ m=\frac{(1-\sqrt a)^2}{\pi a}<0.012. \] Thus \(-0.012-0.012\). Since \(200/27<\zeta\), the absolute logarithmic tail there equals \[ \frac23-2\mathcal P(\zeta^{-1})-\mathcal P(200/27) <\frac23+2(0.0294)-0.6957<0.032 =0.012\alpha_{\rm r}. \] We next bound the variations of \(|h|/u\). Let \(G(u)=h_+(u)/u\). On the positive interval define \(\phi(t)=\log G(e^t)\). Differentiating the rational expression above yields \[ \phi''(t)=\frac9{16} \left(-\frac{c(av-c)}{(a-cv)^2}-\frac1{v^2} +\frac{av-1}{(v-a)^2}\right)<0, \qquad v=\cosh(3t/4). \] For completeness, if \(av\le1\), every displayed term is nonpositive and the first two are negative. If \(av>1\), then on the positive interval \(0c^2\ge(a-cv)^2\); also \(av-c>av-1\). The first negative term therefore dominates the last positive term. The logarithm of \(G\) is strictly concave. The following values locate its unique mode: \[ \begin{array}{c|ccc} t&-0.40&-0.38&-0.36\\ \hline \phi'(t)&(0.055,0.065)&(0.011,0.019)&(-0.036,-0.027). \end{array} \qquad G(e^{-0.38})<0.54. \] It follows that the mode \(u_*\) lies between \(e^{-0.40}\) and \(e^{-0.36}<0.70\). The tangent bound for the concave function \(\phi\) gives \[ \max_{u>0}G(u)<0.54e^{0.019(0.02)}<0.55. \] For \(00.24\) give the sharper bound \(1.10(0.70)-0.24=0.53\). Thus this positive contribution is always below \(0.55\). For the negative part on the left interval \((0,\zeta^{-1})\), \(h_-(u)/u\) decreases: \(v h_-(u)=(cv-a)/(\pi(v-a))\) increases with \(v\), while \(v\) decreases and \(uv\) increases with \(u\). Its scaled variation above \(Q\) is therefore at most \(h_-(Q)<0.012\). There is no left contribution for \(Q\ge1\). On the right interval \((\zeta,\infty)\), the variation of \(h_-\) is at most \(2m<0.024\). The product variation bound applied to \(h_-(u)/u\) gives, with \(q=\max(Q,\zeta)\), \[ Q\operatorname{TV}_{[q,\infty)}(h_-/u) \le Q\left(\frac{0.024}{q}+\frac{0.012}{q}\right) <0.036\min(1,Q/7.5). \] For \(00.60\). Inversion symmetry and \(h<0.451\) give \[ \int_Q^\infty |h(u)|\,\frac{du}{u} <0.393+0.451\log(5/3). \] Hence in this remaining domain \[ \mathcal D_{\ne1}(x,b) <\frac{0.675}{287/300}+\frac{0.55}{287/300} +0.393+0.451\log(5/3)+0.605 <2.510<\frac83. \] For \eqref{eq:first-edge-deficit}, recall \(11\), and \(0.71/b\le1.065\). Combining the correlated ranges in \eqref{eq:first-signed-variation}--\eqref{eq:first-ordinary-integral} gives a strict total in every case: \[ \begin{array}{ll} 10\). Then every derivative of \(\boldsymbol x\mapsto\boldsymbol x-\eta F(\boldsymbol x,\delta)\) has \(\ell^\infty\) row norm at most \(1-\eta\epsilon\). Coordinatewise projection onto \(\mathcal B_{N,s}\) is nonexpansive, so the projected map is a contraction of the rectangle. Its fixed point exists by iteration. The strict boundary signs put that point in the interior, where a projected fixed point satisfies \(F=0\). For any two points of the rectangle, the Jacobian averaged along their segment retains positive diagonal entries and the same strict row margin. It is invertible by the maximum-coordinate estimate, which also gives \(\|H^{-1}\|_\infty\le\epsilon^{-1}\). Applying the averaged Jacobian to the difference of two zeros proves uniqueness. The containment argument at the start places every ordered solution in this rectangle. Finally the analytic implicit function theorem, applied for fixed \(N\) to \eqref{eq:quantile-jacobian}, proves the stated local dependence and continuation. \end{proof} ==== END Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//counting.tex ==== ==== BEGIN Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//cylinder.tex ==== \section{Positive boundary mass and winding polygons}\label{sec:cylinder} We prove two estimates with constants uniform in the circumference: a boundary deficit of order \(N^{-1/4}\), and a mass of order \(N^{-1}\) for one winding polygon whose highest row is fixed. The first step is to prove that the boundary sums over all depths are finite. We then express the mass of one polygon as a polynomial in a tile weight, constrain that polynomial by fusion, and recover both estimates from a strip integral. Retain \(l=\pi/8\) and \(x=(2\cos l)^{-1}\), and put \(c=3/8\). Orient the unit triangular tiling with horizontal sides. A boundary port is the midpoint of a triangle side. A port-to-port path visits each triangle at most once and joins two distinct side midpoints through its centre; a path visiting \(m\) centres has weight \(x^m\). For an integer \(g>0\), let \(K_\infty(g)\) be the total weight of paths in the lower half-plane from horizontal boundary port \(0\) to port \(g\). In the quotient by horizontal translation by an integer \(N\ge1\), let \(K_N(g)\), \(1\le g0}K_\infty(g)\le1. \end{equation} This finite bound is the input for the annular cap construction. The equality and the endpoint tail estimate will follow after the deficit has been evaluated. \subsection{Annular diagram rows and the zero-fugacity cap} We first distinguish the two row operators that enter the argument. The spin row \(T^\zeta(u;\mathbf b)\) of Section~\ref{sec:local} acts on \(V^{\otimes N}\), with a twist \(\zeta^{s_a}\) in the auxiliary trace. The diagram row defined below acts instead on annular arch patterns. The two descriptions give the same scalars when the open diagram ports are closed with the specified caps and the orientations are summed; we do not identify their state spaces. A tile has two incoming legs, auxiliary from the right and quantum from below, and the opposite outgoing legs. Its scalar weights are the functions \(a,b,h,d,e\) of Section~\ref{sec:local}. An empty tile has weight \(1\). A single connection from one input to the other stream's output has weight \(a\), a single straight connection has weight \(h\), and a cup or cap joining the two inputs or the two outputs has weight \(b\). When all four legs are occupied, the double pairing of \(a\)-type has weight \(d\), and the other double pairing has weight \(e\). Contractible closed loops have weight zero, and simple winding loops have weight \(n\). At the physical argument \(t=l\), \[ (a,b,h,d,e)(l)=(x,x^2,x^2,x^2,0). \] These are exactly the weights of dual self-avoiding paths in the two triangles of a split rhombus. Here is the orientation comparison with the spin matrix. Choose \(\zeta\ne0\) with \(n=\zeta+\zeta^{-1}\) and orient each occupied arc. The spins \(+1,-1\) mean respectively along and against the positive stream direction. Straighten the positive stream tangents at the tile ports to parallel forward directions, with the auxiliary stream to the right of the quantum stream at entry and to its left at exit. A through connection then has net turn zero. A clockwise cup or cap has phase \(p=e^{-i\pi/4}\), and a counterclockwise one has phase \(p^{-1}\), exactly as in the terms involving \(C\) and \(p^{-m}\) in \(R(t)\). These spin phases use one quarter of the turn angle, separately from the boundary flux phase \(e^{icW}\) above. The changes of tangent used in straightening cancel at each join, so a closed curve of signed rotation \(\Theta\) has phase \(e^{i\Theta/4}\). Equivalently, the phase counts signed cups and caps and is unchanged by isotopy fixing the port tangents. A contractible simple loop has rotation \(\pm2\pi\); its two orientations therefore contribute \(i\) and \(-i\), with sum zero. A simple winding loop has rotation zero in the flat cylinder frame and seam flux \(\pm1\); its orientation sum is \(\zeta+\zeta^{-1}=n\). This remains true for a diagram closed by arches in earlier rows. For each such arch one chooses an actual representative below the row, with the prescribed port tangents and homotopy, and includes that representative's own turn and seam factors. Representatives above a row are handled in the same way. We also need the interchange and inversion identities with open connections. For interchange, place streams \(1,2,3\) from right to left at entry, all pointing forward, and in reverse order at exit. For fixed occupied boundary slots in a disk, the oriented tensors of distinct noncrossing matchings are linearly independent. Indeed, after straightening the port tangents, the nonzero arc phase depends only on the oriented connectivity. For a chosen matching, assign inward spins to its left endpoints and outward spins to its right endpoints in a linearized boundary order. The result is a Dyck word. At every gap, a matching compatible with this assignment has at least the Dyck height many open arcs; equality at every gap singles out the chosen matching. Ordering matchings by these height tests gives a triangular test matrix with nonzero diagonal. The spin identities of Section~\ref{sec:local} therefore imply the corresponding identities of open diagrams. Sliding a crossing of two auxiliary streams around the \(N\) sites now proves commutation of periodic diagram rows: insert the crossing and its inverse, use interchange at each site, and cancel them after the trace. Any loops formed during this operation are annular loops with the prescribed weights. The proof applies at generic regular parameters and extends as a meromorphic identity. Write \(\mathsf T_n(t)\) for the homogeneous diagram row with argument \(t\) and winding weight \(n\). A state specifies, at the bottom boundary circle, which of the \(N\) slots are vacant and how the occupied slots are joined by noncrossing arches below the circle, including their homotopy classes. No strand runs to infinity. A chosen lift of each arch has endpoint span of absolute value less than \(N\), since an arch disjoint from its translates cannot have interlacing lifted endpoints. There are consequently only finitely many states. To act by \(\mathsf T_n(t)\), attach a row above the input arches, evaluate closed loops, and read the output arch pattern. The preceding diagram interchange gives, for the same winding weight, \[ \mathsf T_n(u)\mathsf T_n(v)=\mathsf T_n(v)\mathsf T_n(u). \] Let \(|0\rangle\) be the all-vacant state and let \(\langle0|\) select its component. Set \(T(t):=\mathsf T_0(t)\), retaining this shorthand for the zero-fugacity diagram row. Every nonempty configuration with empty output contains a closed loop, so \[ \langle0|T(t)=\langle0|. \] We next prove, for each fixed \(N\), \begin{equation}\label{eq:cyl-cap-limit} T(l)^H\longrightarrow|\Psi\rangle\langle0| \qquad(H\to\infty),\qquad \langle0|\Psi\rangle=1. \end{equation} \paragraph{Empty input.} The column \(T(l)^H|0\rangle\) is the positive sum of physical arch systems of depth at most \(H\). For a fixed output pattern, lift each component path to the half-plane. Its endpoint separation has absolute value less than \(N\), and its total weight is bounded by the corresponding \(K_\infty(g)\). Ignoring mutual avoidance bounds the arch-system sum by products of these finite masses. Only finitely many endpoint and homotopy choices occur at fixed \(N\). Thus each component increases to a finite limit. These limits define the finite vector \(\Psi\); its empty component is \(1\). \paragraph{Nonempty input.} Fix a nonempty annular pattern \(\alpha\). It has a physical realization of positive weight in finitely many rows below the cut. To construct one, repeatedly choose a pair whose bounded disk side contains no endpoints still to be joined. Cap that pair across its auxiliary interval, using one \(a\), one \(b\), and the intervening \(h\)'s in bottom-to-top notation, while carrying all other occupations straight. This also permits an interval crossing the seam. The construction terminates, and all its weights at \(l\) are positive. Attach this one fixed realization beneath a diagram counted by the column \(T(l)^H|\alpha\rangle\). The union determines the original bulk diagram, so this attachment is injective. At zero fugacity no closed component survives. Hence every segment belongs to an output arch, and at least one such arch reaches the attached realization, at depth at least \(H-O_\alpha(1)\). The output weight multiplied by the fixed positive cap weight is therefore bounded by the tail of the convergent sum of completed arch systems. That tail tends to zero. Every nonempty input column thus tends to zero, proving \eqref{eq:cyl-cap-limit} in the finite state space. Commutation and the left empty identity now give the common right fixed vector: \[ T(t)|\Psi\rangle\langle0| =\lim_{H\to\infty}T(l)^H T(t) =|\Psi\rangle\langle0|T(t) =|\Psi\rangle\langle0|, \qquad T(t)\Psi=\Psi. \] The limit \eqref{eq:cyl-cap-limit} also says that \(1\) is a simple eigenvalue of \(T(l)\), with all remaining eigenvalues strictly inside the unit disk; the principal nonempty block has spectral radius less than \(1\). These are fixed-\(N\) conclusions. They give a locally analytic eigenpair under a small change of winding fugacity at each fixed \(N\), as used for the higher derivatives in Section~\ref{sec:projection}, but give no uniform spectral gap or uniform analyticity radius in \(N\). \subsection{The one-winding derivative, fusion, and its polynomial} Keep the cap \(\Psi\) at zero fugacity and define \begin{equation}\label{eq:cyl-bare-J} J_N(t):= \left\langle0\left| \left.\partial_n\mathsf T_n(t)\right|_{n=0} \right|\Psi\right\rangle . \end{equation} The derivative selects exactly one winding loop closed by the new row: \(\partial_n n|_{0}=1\), while higher powers of \(n\) vanish. The cap below contains ordinary open arches and no closed loops. At \(t=l\) the new row is the unique highest occupied row of the resulting polygon. Thus each actual unoriented quotient polygon is counted once, without a marked vertex or a division by horizontal translations, and \[ J_N(l)=P_N. \] The fugacity-\(2\) gas used later instead assigns \(2x^{|\gamma|}\) to each polygon. \paragraph{Fusion after closing the caps.} For any fixed annular bottom cap \(\alpha\), the local fusion identity implies a scalar relation before differentiation: \begin{equation}\label{eq:cyl-capped-fusion} \left\langle0\left| \mathsf T_n(t)\mathsf T_n(t+2l)-\mathsf T_n(t+l) \right|\alpha\right\rangle =t^N R_\alpha(t,\zeta), \qquad n=\zeta+\zeta^{-1}, \end{equation} where \(R_\alpha\) is analytic near \(t=0,\zeta=i\). To see the order \(N\), label the two auxiliary spaces \(1,2\). The local product \(R_{23}(t)R_{13}(t+2l)\) preserves \(E_2V\otimes V_3\) and acts there as the single tile \(R(t+l)\). Moreover, \(E_2\) intertwines \(\zeta^{s_1+s_2}\) with \(\zeta^s\), so the trace on that subspace is the fused row. On the quotient auxiliary space, each local product is \(O(t)\): at \(t=0\) it is \(P_{23}R_{13}(2l)\), whose range already lies in \(E_2V\otimes V_3\). Its product around \(N\) sites is therefore \(O(t^N)\). Closing with the chosen oriented representative of \(\alpha\), including its turn and seam factors, converts this spin trace relation to \eqref{eq:cyl-capped-fusion}. There are finitely many cap patterns, each with a finite coefficient in \(\Psi\), so \eqref{eq:cyl-capped-fusion} also holds with \(\alpha\) replaced by \(\Psi\). Near \(\zeta=i\), the quantity \(n=\zeta+\zeta^{-1}\) is a local analytic coordinate because \(dn/d\zeta=2\) there. Differentiate that scalar relation at \(n=0\), with the cap coefficients fixed. The identities \(T(t)\Psi=\Psi\) and \(\langle0|T(t)=\langle0|\) reduce the two derivatives of the product to \(J_N(t)\) and \(J_N(t+2l)\). Hence \begin{equation}\label{eq:cyl-differentiated-fusion} J_N(t)+J_N(t+2l)-J_N(t+l)=O(t^N) \qquad(t\to0). \end{equation} No derivative of \(\Psi\) is taken in this argument. \paragraph{Dependence on the straight weight.} At empty output, a full horizontal auxiliary loop contributes \(h(t)^N\) once. Otherwise the occupied auxiliary intervals have positive integer lengths \(s_1,\ldots,s_r\) and alternate with gaps of positive lengths \(g_1,\ldots,g_r\). Their bottom endpoints are paired by arches in \(\Psi\). We claim that each such arch crosses one consecutive gap. Indeed, in the lift an arch cuts off a bounded disk side. It encloses an even number of other endpoints, so its two endpoints are of opposite types in the alternating list of interval starts and ends. The incident row intervals therefore both point into that disk side or both point out. The first case closes a component inside the bounded lifted disk. In the second case, any enclosed endpoints are paired internally both by arches and by row intervals, again making a closed component inside that disk. Such components are contractible on the cylinder and have zero weight. The disk side must therefore contain no other endpoint, which is precisely a pairing across a consecutive gap. Put \(C_0=2\cos(3l)\). The local weight formulas give, for each occupied interval, \[ a(t)b(t)h(t)^{s_i-1}=C_0(1-h(t))h(t)^{s_i}. \] A labelled cyclic arrangement \(\mathcal A\) records its actual intervals on the \(N\)-site ring and the resulting arches across the gaps; write \(\Psi(\mathcal A)\) for that pattern's cap mass. Consequently \[ \begin{aligned} J_N(t)&=f_N(h(t)),\\ f_N(u)&=u^N+ \sum_{\mathcal A} C_0^{\,r}(1-u)^r u^{s_1+\cdots+s_r}\Psi(\mathcal A). \end{aligned} \] Each actual labelled arrangement is counted once, with no additional marked starting site. The lengths satisfy \(\sum_i(s_i+g_i)=N\), so \(r+\sum_i s_i\le N\). Thus \(f_N\) is a polynomial of degree at most \(N\) and \(f_N(1)=1\). At first derivative order at \(u=1\), only the full ring and the one-interval terms survive. For a gap of length \(g\in\{1,\ldots,N-1\}\), the one interval has \(N\) labelled starting positions, its cap mass is \(K_N(g)\), and the derivative of \((1-u)u^{N-g}\) at \(1\) is \(-1\). We have proved the exact normalization \begin{equation}\label{eq:cyl-polynomial} \deg f_N\le N,\qquad f_N(1)=1,\qquad f_N'(1)=N-C_0N\sum_{g=1}^{N-1}K_N(g)=ND_N. \end{equation} \subsection{The strip operator} Use the strip coordinate \(s=2t-3l\), and write \[ \alpha_{\rm s}:=\frac{\pi}{6l}=\frac43,\qquad \beta_{\rm s}:=\alpha_{\rm s}^{-1}=\frac34,\qquad S_a:=\{s\in\mathbb C:|\Re s|1\); integrating this bound and adding twice the negative mass to \(1/4\) gives \(p_g<1/2\). With these definitions, the degree bound gives the representation \begin{equation}\label{eq:cyl-strip-representation} \begin{gathered} G_f(s)=(I-\mathcal C_g)q(B(s))+\lambda \quad(s\in S_{3l}),\\ q(z)=\sum_{m=1}^N q_mz^m, \end{gathered} \end{equation} for some \(q\) and constant \(\lambda\). To prove \eqref{eq:cyl-strip-representation}, first Fourier transform \(j_f(iy)-f(1)\), which decays exponentially at both ends. Shifting the contour through one real period \(16l\) crosses the mirrored poles at \(7l,9l\), each of order at most \(N\). Their residues give, on \(S_{7l}\), \[ j_f(s)=f(1)+\int_{\mathbb R} \frac{H(k)}{\sinh(8lk)}e^{ks}\,dk, \qquad H(k)=P_e(k)\sinh(lk)+P_o(k)\cosh(lk), \] where \(P_e\) is even, \(P_o\) is odd, and both have degree at most \(N-1\). The displayed symbol therefore gives \[ G_f'=(I-\mathcal C_g)F,\qquad F(s):=\int_{\mathbb R}\frac{kH(k)}{\sinh(6lk)}e^{ks}\,dk. \] Here \(F\) lies in the holomorphic odd Bergman space on \(S_{3l}\), with area norm \(\int_{S_{3l}}|F|^2\,dA\). This space explains the degree bound in \eqref{eq:cyl-strip-representation}. Plancherel on the vertical lines \(s=r+iy\), followed by integration over \(-3l0\), \[ |z|<1+\epsilon,\qquad |\arg(1-z)|<\frac{\pi}{2}+\epsilon,\qquad z\ne1. \] Away from \(z=1\) this follows by compactness and the disk mapping. Near \(1\), inversion of \(\cos(\alpha_{\rm s}s)=d_0(1+z)/(1-z)\) gives both the continued domain and \(e^{-\alpha_{\rm s}|\Im s|}\asymp |1-z|\). Consequently, \[ |\Phi_m(z)| \le C\int_{\mathbb R}e^{-|y|} \exp\!\left(-c_2m|1-z|e^{-\alpha_{\rm s}|y|}\right)\,dy \le C'\min\{1,(m|1-z|)^{-\beta_{\rm s}}\}. \] For sufficiently large \(n\), use Cauchy's coefficient contour with outer radius \(1+\epsilon/2\), indented along \(1-z=r e^{\pm i(\pi/2+\epsilon/2)}\) to \(r=1/n\), then along the small circle of that radius. On the rays, \(|z|^{-n-1}\le C e^{-c n r}\), since their angle is greater than \(\pi/2\); the small circle has length \(O(1/n)\) and bounded \(|z|^{-n-1}\). The preceding bound on \(\Phi_m\) therefore gives \[ |M_{nm}|\le\frac Cn\min\{1,(n/m)^{\beta_{\rm s}}\}. \] Applying the same contour to \((1-z)\Phi_m(z)\) gives the difference bound in \eqref{eq:cyl-M-bounds}. For the finitely many smaller \(n\), use a fixed circle \(|z|=r_*<1\); there \(|1-z|\) is bounded away from zero, so the bound on \(\Phi_m\) is uniform in \(m\) and gives the same estimates after increasing \(C\). The symmetry \(nM_{nm}=mM_{mn}\) supplies the additional \((m/n)^{\beta_{\rm s}}\) factor, and the bound at \(z=0\) gives the stated estimate for \(M_{0m}\). \paragraph{The mass of a high column.} The entry and difference bounds also imply the sharper limit \begin{equation}\label{eq:cyl-column-limit} \sum_{n\ge1}|M_{nm}|\longrightarrow p_g \qquad(m\to\infty). \end{equation} To prove it, consider the finite signed measures \[ \mu_m:=\sum_{n\ge1}M_{nm}\,\delta_{n/m} \quad\text{on }(0,\infty). \] The entry bound makes their total variations uniformly bounded and tight in compact ratio ranges. More precisely, the total variation on \(n/m<\eta\) is \(O(\eta^{\beta_{\rm s}})\), and that on \(n/m>R\) is \(O(R^{-\beta_{\rm s}})\), uniformly in \(m\). For fixed \(a>0\) and all sufficiently large \(m\), choose \(Y_m\to\infty\) so that \(B(iY_m)=e^{-a/m}\). For each fixed \(y\), \[ B(i(Y_m+y))^m\longrightarrow \exp(-a e^{-\alpha_{\rm s}y}). \] Since \(0N\). We claim \begin{equation}\label{eq:cyl-c-margin} \begin{gathered} \sup_N\|c^{(N)}\|_1<\infty,\qquad c_m^{(N)}\longrightarrow0\ \text{for every fixed }m,\\ \limsup_{N\to\infty}\sum_{m0\), choose the remaining columns so that \(d_m\le p_g+\varepsilon\), let \(N\to\infty\), and then let \(\varepsilon\downarrow0\). With \(\|c\|_1=1+\sum_{m0. \end{gathered} \] Together with the fixed-\(N\) nonvanishing in \eqref{eq:cyl-normalized-G}, this also gives a uniform positive lower bound on \(|E_N|\). The positivity of \(E_N\) itself has so far been proved only for sufficiently large \(N\), which is all the asymptotic argument needs. The coefficient recurrence cancels every power below \(N\): \[ \widetilde G(s) =B(s)^N-\sum_{n\ge N}(Mc)_n B(s)^n. \] For real \(y\), \(z=B(iy)\) lies in \((0,1)\). The last strict bound on \(\|Mc\|_1\) and the positive bounds on \(E_N\) therefore give \begin{equation}\label{eq:cyl-positive-G} G(iy)\asymp B(iy)^N>0 \qquad(y\in\mathbb R,\ N\ge N_0), \end{equation} for some fixed \(N_0\), with constants uniform in \(y,N\). This is the two-sided positive estimate needed for both masses. \subsection{Recovering the two positive masses} The two conclusions of Proposition~\ref{prop:cylinder} come from different tails of the strip calculation. The inverse shift kernel at \(s=-l\) gives \(P_N\), while the slower tail of \(g\) determines \(f_N'(1)=ND_N\). The inverse of the shift symbol is \[ (2\cosh(2lk)-1)^{-1}=\frac{\cosh(lk)}{\cosh(3lk)}. \] Its Fourier kernel can be written \[ \mathcal K(y)= \frac{\operatorname{sech}(\alpha_{\rm s}y+i\pi/6) +\operatorname{sech}(\alpha_{\rm s}y-i\pi/6)}{12l}. \] It is analytic for \(|\Im y|<2l\). Uniformly for \(v\) in any closed inner subinterval of \((-2l,2l)\), its tails are \[ \begin{aligned} \mathcal K(u+iv)&\sim c_h e^{-\alpha_{\rm s}u}e^{-i\alpha_{\rm s}v} &&(u\to+\infty),\\ \mathcal K(u+iv)&\sim c_h e^{\alpha_{\rm s}u}e^{i\alpha_{\rm s}v} &&(u\to-\infty). \end{aligned} \qquad c_h=\frac{\sqrt3}{6l}>0. \] The endpoint behavior \(j(s)-1=O_N(e^{-|\Im s|})\) on closed inner substrips of \(S_{7l}\) justifies Fourier inversion of \(j-1\) and \(G-1\). The constant mode has multiplier \(1\), and \(\int\mathcal K=1\). Contour displacement to the imaginary axis, followed by taking real parts, therefore gives \[ P_N=j(-l)=\int_{\mathbb R} \operatorname{Re}\mathcal K(y-il)\,G(iy)\,dy. \] Here \(G(iy)\) is real, and evenness permits the opposite sign of the displacement as well. Since \(\alpha_{\rm s}l=\pi/6\), the real part of the displaced kernel is positively comparable to \(e^{-\alpha_{\rm s}|y|}\) for large \(|y|\). On the other hand \[ B(iy)=\frac{\cosh(\alpha_{\rm s}y)-d_0} {\cosh(\alpha_{\rm s}y)+d_0} =1-\Theta(e^{-\alpha_{\rm s}|y|})\in(0,1). \] On any bounded \(y\)-interval, \eqref{eq:cyl-positive-G} is exponentially small in \(N\). On either tail the same equation, followed by the substitution \(v=e^{-\alpha_{\rm s}|y|}\), gives \[ \int_{|y|>Y}e^{-\alpha_{\rm s}|y|}B(iy)^N\,dy \asymp \int_0^{e^{-\alpha_{\rm s}Y}}e^{-\Theta(Nv)}\,dv \asymp N^{-1} \] for fixed sufficiently large \(Y\). The positive tails dominate the bounded middle, proving \(P_N\asymp N^{-1}\). For the deficit, define \[ \kappa:=2(\cos l+\cos3l)(2\cos2l-1)>0,\qquad \gamma_g:=\frac{1-2^{-1/2}}{\pi}>0. \] At fixed \(N\), expansion of the rational expression for \(j\) at imaginary infinity gives \[ \lim_{y\to\infty}e^y(1-G(iy))=\kappa f_N'(1). \] For \(Q_N\), put \[ D_Q(u):=Q_N(1)-Q_N(B(iu)),\qquad I_m:=\int_{\mathbb R}e^u(1-B(iu)^m)\,du. \] The constant terms cancel in the endpoint difference of \eqref{eq:cyl-normalized-G}, leaving \[ 1-G(iy)=\frac1{E_N} \left(D_Q(y)-\int_{\mathbb R}g(t)D_Q(y+t)\,dt\right). \] For fixed \(N\), \(D_Q(u)=O_N(e^{-\alpha_{\rm s}|u|})\). Also \(g(v)\sim-\gamma_g e^{-|v|}\) and \(e^y|g(u-y)|\le Ce^u\). Since \(\alpha_{\rm s}>1\), the direct term tends to zero after multiplication by \(e^y\), and dominated convergence in the convolution gives the exact fixed-\(N\) limit \[ \kappa f_N'(1)=\frac{\gamma_g}{E_N} \sum_{m=1}^N c_m I_m. \] No limit in \(N\) has been interchanged with this endpoint limit. The integrals \(I_m\) are nonnegative and increase with \(m\), because \(00\). Also a half-plane path from \(0\) to \(N\) exists with positive weight. Together with \eqref{eq:cyl-half-plane} and \(K_N(g)\le K_\infty(g)\), this gives \[ D_N\ge 2\cos(3l)\sum_{g\ge N}K_\infty(g)>0. \] The finitely many remaining circumferences can therefore be absorbed into the comparison constants, completing Proposition~\ref{prop:cylinder}. Finally the deficit also saturates the half-plane bound. If \(C_0=2\cos3l\), then \[ 1-D_N=C_0\sum_{g0}K_\infty(g)\le1. \] Since \(D_N\to0\), the last sum equals \(1\), and the same comparison gives \[ C_0\sum_{g\ge N}K_\infty(g) =1-C_0\sum_{g0}K_\infty(g)=1,\qquad \sum_{g\ge N}K_\infty(g)=O(N^{-1/4}), \] as required by the later boundary comparisons. ==== END Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//cylinder.tex ==== ==== BEGIN Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//higher-jets.tex ==== \subsection{Uniform bounds for every fixed twist derivative} \label{subsec:v-higher-twist-jets} The first winding derivative controls the total mass of one winding polygon. We now bound every fixed derivative. These bounds will control the heights of the polygons before passage from the cylinder to the plane. The analyticity radius in the winding fugacity may depend on the circumference; no uniform radius will be needed. We use the following notation from Section~\ref{sec:cylinder}: \[ l=\frac{\pi}{8},\qquad \alpha_{\mathrm s}=\frac{\pi}{6 l}=\frac43,\qquad \beta_{\mathrm s}=\alpha_{\mathrm s}^{-1}=\frac34, \qquad S_a=\{s\in\mathbb C:|\Re s|J$. The approximate equation then gives the contradiction \[ 1\le\sum_{m\le J}d_m|v_m^{(r)}| +(p+\varepsilon)\|v^{(r)}\|_1+o(1) \longrightarrow p+\varepsilon<1. \] We also need the normalization of the one-dimensional homogeneous solution. Let \[ Q_N(z)=\sum_{m=1}^N c_mz^m,\qquad c_N=1,\qquad c_n=\sum_{m=1}^N M_{nm}c_m\quad(1\le n0. \end{equation} For completeness, $\mathcal E(Q_N)$ is nonzero at every remaining finite $N$ as well. Indeed the existing first jet has $P_1=0$; its order-$N$ zero makes the polynomial in \eqref{eq:v-jet-strip-representation} a multiple of $Q_N$, while its endpoint is $H_1(i\infty)=1$. Thus $\mathcal E(Q_N)=0$ is impossible. Together with \eqref{eq:v-jet-positive-endpoint}, this gives a uniform positive lower bound on $|\mathcal E(Q_N)|$. \paragraph{The forcing coefficients.} We proceed by induction on $k$. Assume the two strip bounds in the lemma for every lower order. When $k=1$ there are no lower-order terms and the forcing below is zero. For $k\ge2$, let \[ b_m=[z^m]P_k(s(z)). \] The function $P_k$ is even: replacing $s$ by $-s$ and $i$ by $k-i$ leaves its sum unchanged. Hence it descends holomorphically to the disk. We claim \begin{equation}\label{eq:v-jet-forcing-lone} \sum_{m=0}^{N-1}|b_m|\le C_k. \end{equation} Here is a contour proof, including the cancellation across its slit. Temporarily write \[ A=F_N(s-2 l,\nu),\quad C=F_N(s,\nu),\quad D=F_N(s+2 l,\nu),\quad E=AD-C, \qquad a_i=[\nu^i]A. \] Define the replacement \[ R_k(s)=\sum_{i=1}^{k-1}a_i[\nu^{k-i}](C/A). \] The exact identity $D-C/A=E/A$ gives \begin{equation}\label{eq:v-jet-replacement-identity} R_k-P_k=-\sum_{i=1}^{k-1}a_i[\nu^{k-i}](E/A). \end{equation} Only defect coefficients $G_r$ with $r0$ and $\theta\in(\pi/2,\pi/2+\alpha_{\mathrm s}\varepsilon_s)$ such that this branch is holomorphic on \[ \mathcal D=\{z:|z|<1+\eta,\quad |\arg(1-z)|<\theta\}\setminus[b_*,1), \] and satisfies \begin{equation}\label{eq:v-jet-inverse-domain} 0<\Re s(z)<3 l+\varepsilon_s, \qquad e^{\alpha_{\mathrm s}|\Im s(z)|}\le\frac C{|1-z|}. \end{equation} These facts follow directly from the inverse formula: away from $1$ one uses compactness and the disk mapping, while near $1$ the logarithmic inverse has real part at most $\alpha_{\mathrm s}^{-1}|\arg(1-z)|+o(1)$. The other branch point, $(1+d_0)/(1-d_0)$, is excluded by taking $\eta$ small. On the two banks of $b_*0$, the explicit formula for $B$ gives \[ |B(s)|\le\exp(-c_\eta e^{-\alpha_{\mathrm s}|\Im s|}) \quad(|\Re s|\le3 l-\eta). \] Since $e^{-cu}\le C_c\min\{1,u^{-2}\}$, the last bound on $G_k$ and the inductive sharp bounds for the two factors of $P_k$ imply \begin{equation}\label{eq:v-jet-sharp-forcing} |H_k(s)|\le C_{k,\eta} \min\left\{1,\frac{e^{2\alpha_{\mathrm s}|\Im s|}}{N^2}\right\} \quad(|\Re s|\le3 l-\eta). \end{equation} Convolution improves the squared scale to the desired scale. Indeed, with $a=e^{\alpha_{\mathrm s}|y|}/N$, the case $a\ge1$ uses $\int\mathcal K=1$. When $a<1$, split at $T=\alpha_{\mathrm s}^{-1}\log(1/a)$ and use \[ \int_{\mathbb R}e^{-\alpha_{\mathrm s}|t|} \min\{1,a^2e^{2\alpha_{\mathrm s}|t|}\}\,dt =\frac{4a-2a^2}{\alpha_{\mathrm s}} \le\frac{4a}{\alpha_{\mathrm s}}. \] Equations \eqref{eq:v-jet-inverse-formula} and \eqref{eq:v-jet-sharp-forcing} therefore prove \eqref{eq:v-jet-sharp-conclusion} on each closed inner substrip of $S_{3 l}$. To reach the asserted width, take $\eta=\min\{\delta, l\}$. For $3 l-\eta\le\Re s\le5 l-\eta$, all arguments on the right of \eqref{eq:v-jet-edge-rearrangement} lie in $|\Re s|\le3 l-\eta$. Apply the sharp bounds there, and use $\min\{1,a^2\}\le\min\{1,a\}$. The negative edge is identical. This proves \eqref{eq:v-jet-sharp-conclusion} and closes the induction. For $k=1$, the argument starts with $P_1=0$ and no contour forcing. \paragraph{The derivative at $h=1$.} No sign of a higher jet is needed for this step. Define \[ \kappa=2(\cos l+\cos3 l)(2\cos2 l-1)>0, \qquad \gamma_g=\frac{1-2^{-1/2}}\pi>0. \] For fixed $N,k$, expansion of the rational function $h$ at imaginary infinity gives \begin{equation}\label{eq:v-jet-endpoint-first-limit} \lim_{y\to\infty}e^y (\mathbf1_{k=1}-H_k(iy))=\kappa f_{N,k}'(1). \end{equation} For the polynomial $q$ found above, put \[ D_q(u)=q(1)-q(B(iu)),\qquad I_m=\int_{\mathbb R}e^u(1-B(iu)^m)\,du. \] The constant $\lambda$ cancels in the endpoint difference, so \[ \mathbf1_{k=1}-H_k(iy) =D_q(y)-\int_{\mathbb R}g(t)D_q(y+t)\,dt. \] For each fixed polynomial, $D_q(u)=O(e^{-\alpha_{\mathrm s}|u|})$. Also $g(v)\sim-\gamma_g e^{-|v|}$ and $e^y|g(u-y)|\le Ce^u$. Because $\alpha_{\mathrm s}>1$, dominated convergence gives \[ \lim_{y\to\infty}e^y (\mathbf1_{k=1}-H_k(iy)) =\gamma_g\sum_{m=1}^N q_m I_m. \] Here $00\), conjugate the trace over two auxiliary spaces by their \(R\), then use interchange at each site and conservation of the auxiliary spins. The same assertion holds for the empty list with the evident convention. Interchange also gives the adjacent exchange relation \[ T^\zeta(u;\ldots,A,B,\ldots)\check R_{j,j+1}(A-B) =\check R_{j,j+1}(A-B)T^\zeta(u;\ldots,B,A,\ldots). \] The analogous cyclic operation is as follows. Rotate the tensor factors one place to the left and multiply by \(\zeta^{s_1}\) on the site sent to the end. This map intertwines the row with the cyclically rotated list, by conservation in slot 1 and cyclicity of the auxiliary trace for entries at separate sites. Shifting a site parameter by \(\pi\) conjugates the row by the sign of occupation at that site. Finally, specialize adjacent entries to \(d_0+rl,d_0\). The row intertwines through \(E_r\) with the shorter list obtained by replacing this pair with \(d_0+l\) for \(r=2\), or by deleting it for \(r=3\). This follows by inverting the fusion identity with site 3 as the auxiliary space and \(t=d_0-u\). The row relations, like the local ones, are meromorphic identities in all parameters. ==== END Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//local.tex ==== ==== BEGIN Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//main.tex ==== \documentclass[11pt]{article} \pdfinfoomitdate=1 \pdftrailerid{} \pdfsuppressptexinfo=15 \usepackage[T1]{fontenc} \usepackage[utf8]{inputenc} \usepackage{lmodern} \usepackage{amsmath,amssymb,amsthm,mathtools,mathrsfs} \usepackage[margin=1in]{geometry} \usepackage{booktabs,longtable,array} \usepackage{graphicx,tikz} \usepackage{microtype} \usepackage[colorlinks=true,linkcolor=blue!50!black,citecolor=blue!50!black,urlcolor=blue!50!black]{hyperref} \usepackage{bookmark} \makeatletter \renewcommand*\l@subsection{\@dottedtocline{2}{1.5em}{3.6em}} \renewcommand*\l@subsubsection{\@dottedtocline{3}{3.8em}{4.6em}} \makeatother \hypersetup{pdfauthor={OpenAI},pdftitle={Cap-selected amplitudes and triangle chords for honeycomb walks}} \allowdisplaybreaks[2] \numberwithin{equation}{section} \newtheorem{theorem}{Theorem}[section] \newtheorem{lemma}[theorem]{Lemma} \newtheorem{proposition}[theorem]{Proposition} \newtheorem{corollary}[theorem]{Corollary} \newtheorem{claim}[theorem]{Claim} \theoremstyle{definition} \newtheorem{definition}[theorem]{Definition} \theoremstyle{remark} \newtheorem{remark}[theorem]{Remark} \DeclareMathOperator{\Pf}{Pf} \DeclareMathOperator{\Tr}{Tr} \DeclareMathOperator{\Res}{Res} \DeclareMathOperator{\diam}{diam} \DeclareMathOperator{\dist}{dist} \DeclareMathOperator{\supp}{supp} \DeclareMathOperator{\rank}{rank} \DeclareMathOperator{\Span}{span} \title{Cap-selected amplitudes and triangle chords for honeycomb walks} \author{OpenAI} \date{September 26, 2026} \begin{document} \maketitle \begin{abstract} At the critical honeycomb fugacity, self-avoiding port-to-port chords in an equilateral lattice triangle of side $R$, summed over both boundary endpoints and restricted to diameter at least $R/100$, have partition sum comparable to $R^{3/4}$ and mean length $R^{4/3+o(1)}$. We also determine the amplitude selected by two vacuum caps on a cylinder of circumference $N$ and prove that it grows as $N^{1/6+o(1)}$. \end{abstract} \setcounter{tocdepth}{1} \tableofcontents \clearpage \input{sections/introduction} \input{sections/local} \input{sections/cylinder} \input{sections/vacuum} \input{sections/residue} \input{sections/common-twist} \input{sections/site-moments} \input{sections/repeated-moments} \input{sections/counting} \input{sections/physical-amplitude} \input{sections/asymptotics} \input{sections/projection} \input{sections/chords} \begin{thebibliography}{99} \input{references} \end{thebibliography} \end{document} ==== END Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//main.tex ==== ==== BEGIN Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//physical-amplitude.tex ==== \section{Selecting the physical coefficient}\label{sec:physical-amplitude} We prove that the coefficient carried by the real counting roots is the coefficient selected by the two vacuum caps. First we identify the roots at the vacuum twist \(K=i\). At generic nearby parameters, the common-twist trace then isolates their ordinary \(W/J\) term. A separate estimate for the two ends of the cylinder makes the cap coefficient analytic at fixed \(N\); this permits normalization at \(K=i\) and continuation to the physical twist \(K=1\). The right cap keeps twist \(i\), and the left cap keeps twist \(-i\), throughout. The polynomial construction below works for every \(N\ge1\). The identification with the counting branch and the final coefficient formula use \(N\ge N_0\) from Proposition~\ref{prop:real-counting}. \subsection{The positive roots at the vacuum endpoint} Set \(l=\pi/8\) and \(q=e^{2il}\). For this subsection take additive sites \(b_a=i\eta_a\), so their multiplicative coordinates \(\xi_a=e^{2ib_a}=e^{-2\eta_a}\) are positive. Repeated sites are allowed. For \(z=e^{2iw}\), put \[ B(z)=\prod_{a=1}^N\prod_{j=3,4,5}(z-q^j\xi_a). \] A polynomial is said to have the missing modes if its coefficients of \(z^{4m+1}\) all vanish. \begin{proposition}[Vacuum endpoint]\label{prop:positive-vacuum-roots} For every \(N\ge1\) and every multiset \(\xi_1,\ldots,\xi_N>0\), there is a unique monic real polynomial \(P_0\) of degree \(N\) such that \(H=BP_0\) has the missing modes. The zeros of \(P_0\) are positive and simple. Write these zeros as \(e^{-2y_j}\), with \(y_1<\cdots0 \quad (x>0). \] The missing modes impose \(N\) real homogeneous equations on the \(N+1\) coefficients of a real polynomial of degree at most \(N\). Choose a nonzero solution \(P_0\). The positive-kernel argument belongs to Schoenberg's P\'olya-frequency and variation-diminishing theory~\cite{Schoenberg1950PFII}. We include the transform factorization and the strict confluent-minor argument because repeated sites require derivatives of the kernel. For a site value \(c=\xi_a^{4/3}\), integrate \[ H(z)(-z)^{1/3} \frac{dz}{z\prod_a(z^4+\xi_a^4)\bigl(c-(-z)^{4/3}\bigr)} \] around a keyhole about the positive ray, using the principal logarithm of \(-z\). If the site has multiplicity \(m\), also differentiate in \(c\) through order \(m-1\). The end circles vanish. At \(-\xi_a,q^3\xi_a,q^5\xi_a\) the poles are removable. At \(q\xi_a,q^7\xi_a\) the residues cancel to the full site multiplicity: the missing modes give \[ iH_j+H_{j+2}-iH_{j+4}-H_{j+6}=0, \] and hence \(H_7=iH_1\) modulo that multiplicity near \(\xi_a\). The local coordinates \((-z)^{4/3}\) at the two poles agree as germs, while the one-third powers give the cancelling phases. Writing \(r_0=x^{4/3}\), the jump therefore gives orthogonality of \(P_0(x)\), for the positive measure \[ d\mu(x)=\frac{B(x)x^{1/3}}{x\prod_a(x^4+\xi_a^4)}\,dx , \] against all the indicated \(c\)-derivatives of \[ \frac{c}{c^2+c r_0+r_0^2}. \] This family, with any initial collection of derivatives at each repeated parameter, has nonsingular interpolation matrices at distinct positive arguments. To see strictness, pass to logarithmic coordinates and remove positive factors. The kernel becomes \((1+2\cosh(u-v))^{-1}\). Its bilateral Laplace transform is proportional to \[ \frac{\sin(\pi t/3)}{\sin(\pi t)} \quad\text{and hence to}\quad \prod_{k=1,2}\Gamma((k+t)/3)\Gamma((k-t)/3), \] by the substitution to a positive half-line, the partial-fraction Beta integral, and the reflection formula. Thus the kernel is a positive constant times the convolution of the four kernels \[ f_{k,\pm}(s)=\exp\bigl(\pm ks-\exp(\pm3s)\bigr),\qquad k=1,2. \] Each matrix \(f_{k,\pm}(u-v)\) has strictly positive ordered minors, including confluent column minors with derivative orders starting at zero. After stripping positive factors, it is \(\exp(-XY)\) with \(X,Y\) oppositely ordered. The distinct row exponentials form an extended Chebyshev system in the column variable, by Rolle induction after division by the first exponential. Confluence gives the Vandermonde sign, which is positive in these orders. Integral determinant multiplication on ordered arguments preserves this strictness under convolution; differentiation under the integrals is valid because the differentiated kernels retain integrable exponential decay. If \(P_0\) had fewer than \(N\) sign changes on the positive ray, take one more of the preceding initial derivatives than the number of sign changes. The interpolation determinant vanishing at those change points supplies a function whose signs alternate strictly on their complementary intervals. Its product with \(P_0\) has one strict sign, contradicting orthogonality. Thus \(P_0\) has \(N\) positive simple zeros and degree \(N\). Two projectively different solutions would have a nonzero linear combination of degree below \(N\), which is impossible by the same argument. This proves uniqueness after making \(P_0\) monic. We next identify the vacuum TQ function. Formula~\eqref{eq:vacuum-TQ-quotient} follows term by term from the shifted factors in \(L/d_{\rm full}\); the constants are fixed, for example, at \(z=0\), where the three terms are \(K,1,K^{-1}\). The same quotient formula holds for any degree-\(N\) root polynomial with nonzero constant term. For the constructed \(H\), at every site and to its full multiplicity the eight germs have the form \[ (H_0,\ldots,H_7)=(a,b',c',0,0,0,c'+ia,ib'). \] This follows from the three zeros supplied by \(B\) and the missing-mode relation. Substitution, also after each cyclic shift, shows that \[ iH_1H_6+H_5H_0-iH_4H_7-H_7H_6 \] vanishes at all eight rotations of every site, to their multiplicities. These are \(8N\) nonzero zeros counted with multiplicity. The polynomial has degree at most \(8N\) and also vanishes at zero, so it is identically zero. Hence \(F_i=1\). The poles at the shifted roots in \(L\) must consequently cancel. Their first pole coefficients are nonzero units times \(i+D_j\), as in \eqref{eq:common-twist-pole-coefficients}. The units are nonzero because the roots and sites are positive and the roots are distinct; an unshifted equality between a root and a site still gives the nonzero factors \([l]\) and \([2l]\). Thus \(D_j=-i\). It remains to determine the integer labels of these phases. The integrated masses of \(p\) and \(r\) are \(3/4\) and \(1/4\). The product defining \(D_j\), without a self factor, starts at \(q\) as \(y\to-\infty\), while the self primitive at zero is \(1/8\). The phase derivative formulas in Section~\ref{sec:counting} therefore give \[ D_j=\exp\bigl(2\pi iT_{\boldsymbol\xi}(y_j)\bigr). \] We determine the labels first when the site scales are very widely separated. Let one site be \(\epsilon\to0\) while the others remain fixed. Normalize the coefficients of \(H\) to be bounded with a nonzero projective limit. That limit is divisible by \(z^3B_{\rm red}\), where \(B_{\rm red}\) omits this site. Write \(H=\sum h_mz^m\). Evaluation at \(q^{3,4,5}\epsilon\), using \(h_1=0\), shows \[ \max\{|h_0|,\epsilon^2|h_2|,\epsilon^3|h_3|\}=O(\epsilon^4). \] Indeed the matrix of \(1,z^2,z^3\) at \(q^{3,4,5}\) is invertible; its relevant symmetric sum is \(1+\sqrt2\). Hence the limit is divisible by \(z^4\). After dividing by \(z^4\), the missing modes and the uniqueness already proved for \(N-1\) sites identify the limit as \(H_{\rm red}\), up to scale; for the empty reduced list when \(N=1\), set \(H_{\rm red}=1\), and the assertion is immediate. It has full degree. Consequently exactly one positive zero of \(P_0\) tends to zero, and the others tend to the reduced zeros. Since \([z]H=0\), the sum of the inverses of all zeros of \(H\) is zero; the three new site zeros contribute \(-(1+\sqrt2)/\epsilon\). The small root is therefore asymptotic to \(\epsilon/(1+\sqrt2)\). Iterating this separation gives configurations with one root comparable to each site and arbitrarily large gaps between their scales. In increasing rapidity order, there \[ T_{\boldsymbol\xi}(y_j)\in j-1+[1/8,7/8]+o(1). \] Since \(D_j=-i\), the only possible phase label in this interval is \(j-1+3/4\). Finally, the monic solution varies continuously in the positive sites by projective uniqueness and full degree. Its ordered zeros remain simple and positive, so the integer phase labels cannot change. This proves the stated counting equations for all positive sites. \end{proof} At homogeneous sites \(\xi_a=1\), these phase labels form a contiguous ordered solution at \(\delta=1/4\). The containment estimate in the proof of Proposition~\ref{prop:real-counting} puts every such solution in its quantile rectangle, where uniqueness identifies it with that proposition's solution. For \(N\ge N_0\), the implicit function theorem then continues the same roots along \[ K=e^{i\theta},\qquad \delta=\theta/(2\pi),\qquad 0\le\theta\le\pi/2. \] Indeed \(T(y_j)=j-\tfrac12+\delta\) is equivalent to \(D_j=-K\). The roots remain distinct and positive in multiplicative coordinates. The branch is analytic locally in \(K\) and under small complex changes of the additive sites. \subsection{The isolated ordinary term at generic parameters} Return to additive site notation \(\mathbf b\). On the selected local root branch write \[ a(\mathbf b,K)=\frac{W}{J},\qquad J=\det\bigl(\partial_{t_j}\log D_i\bigr), \] and let \(F(u;\mathbf b,K)\) be its normalized TQ function. On a solution it can be evaluated as \(P(u)/d_{\rm full}(u)\), using the Laurent polynomial part from Section~\ref{sec:common-twist}; this expression has already cancelled the apparent poles at shifted roots. \begin{lemma}[Isolation of the ordinary trace term]\label{lem:ordinary-physical-term} Fix \(N\ge N_0\). In a sufficiently small neighborhood of the homogeneous vacuum parameters, take generic sites, a generic common twist \(K\) for which the toric ends in Proposition~\ref{prop:common-twist} are excluded, and a generic row argument \(u\) near \(l\). For all sufficiently large \(h\), the finite exponential-polynomial expansion of \[ G_h=\Lambda\bigl(T^{K^{-1}}(u;\mathbf b)\bigr)^h\Psi \] contains the term \[ \frac{N!a(\mathbf b,K)}{\tau(1)}\,F(u;\mathbf b,K)^h . \] Its coefficient is nonzero, and no other term has that exponential base. Here \(\tau(1)\) is nonzero generically and is independent of \(K\). \end{lemma} \begin{proof} Theorem~\ref{thm:row-moments} and Corollary~\ref{cor:common-twist-jets} apply on the stated generic locus. The jet functional is fixed as \(u,h\) vary; only the Laurent polynomial insertion \((P(u)/d_{\rm full}(u))^h\) changes. Differentiating its powers gives a finite sum of polynomials in \(h\) times exponential bases. Terms with zero base disappear for sufficiently large \(h\). At every jet support point the cleared TQ relation gives \(L=P\), so each remaining base is the value at \(u\) of that point's exact rational TQ function. A generic \(u\) separates nonidentical rational functions. We show that no other root polynomial has the same rational TQ function as the selected one at generic nearby parameters. In \eqref{eq:vacuum-TQ-quotient}, normalize the two corresponding polynomials \(H\) at zero. Multiplication of \(H\) by any invertible series in \(z^4\) leaves the quotient unchanged: every product in its numerator and denominator has one even and one odd shift. Suppose lower coefficients of two such series have been matched. At degree \(m\), equality of the quotients determines their relative coefficient with multiplier \[ \mu_m(K)= K(a_m-a_m^7)+(a_m^5+1-a_m^7-a_m^6) +K^{-1}(a_m^4-a_m^6),\qquad a_m=q^m . \] For \(4\nmid m\), its coefficient \(a_m-a_m^7=a_m-a_m^{-1}\) of \(K\) is nonzero. Thus, after excluding finitely many values of \(K\), all these multipliers are nonzero. At degrees divisible by four use a series in \(z^4\) to match the coefficients; at every other degree the triangular equation then matches them. The quotient of the two root polynomials is therefore invariant under \(z\mapsto iz\). The selected roots near the vacuum remain separated from their three quarter-turns. If that rational quotient had a pole at a selected root, quarter-turn invariance would require poles at all its rotations, where its denominator has no zero. It has no finite poles, and the degree bound makes it constant. The root polynomials are equal. At generic nearby sites the selected solution lies in the regular configuration locus and is transverse. Its local contribution to the continued trace is the ordinary \(W/J\) evaluation by Corollary~\ref{cor:common-twist-jets}. The \(N!\) labelings have the same evaluation and Jacobian, and the preceding uniqueness excludes every other support polynomial from their exponential base. The factors in \(W\) are nonzero there, and \(J\ne0\) by continuation from Proposition~\ref{prop:real-counting}. Finally, the independent-twist residue formula with the empty insertion makes \(\tau(1)\) independent of the twists and nonzero generically, as used in the distinct-site normalization. This proves the lemma. \end{proof} The selected quotient \(a=W/J\) also extends across an unshifted root--site coincidence on the positive slice. Such a point was excluded from the deliberately small regular configuration locus, but \(S(0)=-1\), \(V(0)=1/[l]^2\), and the marked-pole site numerators retain nonzero \([l][2l]\) factors there. With distinct roots and \(J\ne0\), the equations and \(W\) are holomorphic, and the implicit function theorem gives the same local inverse branch as the limits of generic independent-twist points. Along the whole homogeneous real branch \(W\ne0\) and \(J\ne0\), so \(a\) is analytic and nonzero locally at every point. The finite-jet representation in Lemma~\ref{lem:ordinary-physical-term} is used only at generic common twists where the toric ends have been excluded. The endpoint identity below will instead concern the analytic continuation of this selected quotient and of the independently defined cap coefficient. \subsection{The fixed-\(N\) annular gap and the cap ends} Let \(\mathsf M_N(\mathbf b,u,K)\) be the finite annular diagram row of Section~\ref{sec:cylinder}, with local arguments \(u-b_a\), contractible loop weight zero, and winding weight \(n=K+K^{-1}\). Let \(\mathsf A_N\) be its principal submatrix on nonempty annular states. At homogeneous sites and \(K=e^{i\theta}\), its weights are nonnegative for \(l\le u<2l\) and \(0\le\theta\le\pi/2\). We describe precisely the finite sums at the ends of \(G_h\). Expand the caps in their finite disk-matching bases and the rows in tile diagrams, summing the orientations of each connected component. A horizontal cut is empty when all its quantum ports are vacant. Let \(E^-_a\) be the sum from the bottom cap through \(a\) rows to an empty cut, with every earlier cut nonempty. Let \(E^+_b\) be the sum from an empty cut through \(b\) rows to the top cap, with every later cut nonempty. Set \(E^-_0=E^+_0=1\), the empty coefficients of the normalized caps. Let \(C_h\) be the sum of diagrams in \(G_h\) with every cut nonempty, including the two end cuts. These are finite sums whenever the cap entries are regular. Put \(Z_m^{00}=\langle0|\mathsf M_N^m|0\rangle\). At an empty cut the orientation sums split, and the portion between two empty cuts has exactly the ordinary \(0,n\) loop weights. Choosing the first and last empty cuts therefore gives the exact identity \begin{equation}\label{eq:cap-end-decomposition} G_h=C_h+\sum_{\substack{a,b\ge0\\a+b\le h}} E^+_b\,Z_{h-a-b}^{00}\,E^-_a . \end{equation} The case \(a+b=h\) has a single empty cut and \(Z_0^{00}=1\). At \(h=0\), null-loop contraction gives \(C_0=0\), consistently with \(\Lambda\Psi=1\). \begin{proposition}[Fixed-\(N\) annular gap and cap estimate]\label{prop:fixed-N-cap-estimate} Fix \(N\ge N_0\). There is \(\varepsilon_N>0\) with the following properties. For every \[ p_0=(\mathbf0,u_0,e^{i\theta_0}),\qquad l\le u_0\le l+\varepsilon_N,\quad 0\le\theta_0\le\pi/2, \] the matrix \(\mathsf M_N(p_0)\) has a simple dominant eigenvalue \(\lambda(p_0)>0\), and \[ \rho\bigl(\mathsf A_N(p_0)\bigr)<\lambda(p_0). \] The right cap \(\Psi\) at twist \(i\) and the left cap \(\Lambda\) at twist \(-i\) are holomorphic in the sites near \(\mathbf0\). At homogeneous sites they are limits of rows from vacuum-empty ends at their respective twists, with reversal and spin negation on the upper end as in the definition of \(\Lambda\). There is a complex neighborhood \(\mathcal U\) of \(p_0\), positive constants \(C,\gamma,\beta\), and holomorphic functions \(\lambda,c_{\rm phys}\) on \(\mathcal U\), with \[ \gamma<\beta<\inf_{p\in\mathcal U}|\lambda(p)|, \] such that, for \(p=(\mathbf b,u,K)\in\mathcal U\) and all \(h\ge0\), \begin{equation}\label{eq:fixed-N-cap-bounds} \begin{aligned} |E^-_h(p)|+|E^+_h(p)|+|C_h(p)|&\le C\gamma^h,\\ |G_h(p)-c_{\rm phys}(p)\lambda(p)^h|&\le C\beta^h. \end{aligned} \end{equation} Here \(\lambda(p)\) continues the stated annular eigenvalue. At \(K=i\) in a vacuum neighborhood, \(\lambda=1\) and \(c_{\rm phys}=1\). The gaps, constants and neighborhoods in this proposition may depend on \(N\). \end{proposition} \begin{proof} We first establish the annular gap. When \(n>0\), the nonnegative annular matrix is primitive. From a nonempty state, connect cyclically consecutive occupied ends of different arcs, carrying all other ends straight, until only one arc remains. Two single turns across the intervening empty interval perform each merger. Close the last arc across its complementary gap, producing a winding loop with positive weight \(n\). Conversely, reverse deletions of minimal disk-side arcs to build every state from the empty state, using the two output turns. These moves have positive weight also at \(u=l\), and the empty state has a positive self-step. For completeness, the finite matrix fact used here is as follows. If a nonnegative matrix \(M\) has a strictly positive power, its spectral radius \(r>0\) has a strictly positive eigenvector: normalize \(\sum_{h\ge0}t^{-h}M^h{\bf1}\) as \(t\downarrow r\). Its norm diverges because the entry sum of \(M^h\) dominates \(r^h\), and a limit is the required eigenvector. Conjugating by that vector and dividing by \(r\) makes \(M\) row-stochastic. Equality in the sup-norm averaging inequality for a positive power forces any unit-circle eigenvector to be constant, while bounded powers exclude Jordan growth. Thus \(r\) is simple and strictly dominant. Deleting a state strictly decreases the spectral radius, because a sufficiently high power of the restricted substochastic matrix loses mass in every row. Applied to the empty state, this gives \(\rho(\mathsf A_N)<\lambda\). At \(n=0,u=l\), Section~\ref{sec:cylinder} proves \[ (\mathsf M_N(\mathbf0,l,i))^h \longrightarrow |\Psi_{\rm ann}\rangle\langle0|. \] The empty functional is fixed. In finite dimension this convergence gives the simple eigenvalue \(1\) and spectral radius below \(1\) on the nonempty block. Continuity supplies the same gaps for \(u\) near \(l\) at \(n=0\). Compactness of \(0\le\theta\le\pi/2\) now supplies an \(\varepsilon_N>0\) for the stated interval. This compactness is for the fixed matrix size \(N\). We next obtain regular caps. At zero loop weights, forget the annular homotopy of every arc and send the resulting disk matching to its tensor in \(D\). This map is onto: every disk matching can be drawn in a cut cylinder. It intertwines the annular row with the spin row of twist \(i\). Indeed the auxiliary--spectator contraction of Section~\ref{sec:vacuum} straightens an arc or closes it with weight zero; the exterior contraction has the same alternatives and phases. Forgetting which side contains the annular hole is the only change in connectivity. The image of the nonempty annular subspace is \(D'\), the subspace of empty coefficient zero. Hence the restriction of the spin row to \(D'\) at homogeneous \(u=l\) is a quotient of the nonempty annular block and has spectral radius below one. The normalized fixed vector of this test row in \(D\) is therefore unique and holomorphic in nearby sites, by solving its finite linear equations on \(D'\). Commutation and empty normalization make it the simultaneous vacuum, so it agrees with \(\Psi\) at generic sites. Reversal, transposition and spin negation give the same conclusion for \(\Lambda\) at twist \(-i\). Projecting the preceding annular limit proves the asserted limits from the lower empty end; applying the same reversal and negation proves the upper limit. We now bound the end sums. Work first at a real point \(p_0\) in the proposition, and abbreviate \(\mathsf A=\mathsf A_N(p_0)\). Choose \[ \rho(\mathsf A)0\) and \(r\le N\), \begin{equation}\label{eq:cap-word-bound} \|\mathsf U_m\|\le C_N(m+1)^N a^m . \end{equation} The same estimate holds in a sufficiently small complex neighborhood with a slightly larger exponential rate. To see this, majorize the moduli of the local tile weights by nearby positive weights, and use \(|K+K^{-1}|\) for an ordinary winding loop. The resulting nonnegative block \(\widehat{\mathsf A}\) tends to \(\mathsf A\) as the neighborhood shrinks; the corresponding exceptional matrix remains bounded. A fixed resolvent circle gives a uniform power bound for \(\widehat{\mathsf A}\) with rate still below \(|\lambda|\). A boundary component can acquire an additional seam modulus under a complex perturbation. If the size of that perturbation is \(\epsilon\), its total cost through \(m\) rows is at most \(\exp(O(\epsilon Nm))\), since only the finitely many boundary strands are involved. Choose the neighborhood small enough that this cost, and then the polynomial in \eqref{eq:cap-word-bound}, can be absorbed into a rate \(\gamma<\inf|\lambda|\). The caps are bounded there by their holomorphy. This proves the first estimate in \eqref{eq:fixed-N-cap-bounds}. Finally define the convergent end series \[ \mathcal E^\pm(z)=\sum_{m\ge0}E^\pm_m z^m,\qquad \mathcal C(z)=\sum_{m\ge0}C_mz^m . \] They are holomorphic for \(|z|<1/\gamma\), jointly with the parameters after shrinking the neighborhood. The middle series is \[ \mathcal Z^{00}(z) =\langle0|(I-z\mathsf M_N)^{-1}|0\rangle =\frac{d}{1-z\lambda}+\mathcal R(z), \qquad d=\langle0|\Pi_\lambda|0\rangle, \] where \(\Pi_\lambda\) is the spectral projector. The annular gap makes \(\mathcal R\) holomorphic in a disk of radius strictly larger than \(1/|\lambda|\), uniformly on a smaller neighborhood. Summing \eqref{eq:cap-end-decomposition} gives \[ \sum_{h\ge0}G_hz^h =\mathcal C(z)+\mathcal E^+(z)\mathcal Z^{00}(z)\mathcal E^-(z). \] Its only possible pole in a suitable larger disk is at \(z=\lambda^{-1}\), of order at most one. The coefficient of \((1-z\lambda)^{-1}\) is \[ c_{\rm phys} =d\,\mathcal E^+(\lambda^{-1})\mathcal E^-(\lambda^{-1}), \] which is holomorphic in the parameters. Subtracting this polar term and applying Cauchy's coefficient bound on a fixed circle beyond \(z=\lambda^{-1}\) gives the second estimate in \eqref{eq:fixed-N-cap-bounds}, with \(\gamma<\beta<\inf|\lambda|\). At \(K=i\), the bulk spin twist \(K^{-1}=-i\) agrees with the left vacuum twist. Thus \(\Lambda T^{-i}(u)=\Lambda\) and \(G_h=\Lambda\Psi=1\). The empty annular functional is fixed at zero loop weight for all nearby local arguments, so the annular branch has \(\lambda=1\). The estimate forces \(c_{\rm phys}=1\). \end{proof} \subsection{Continuation to the homogeneous physical parameters} The generic trace term and the cap estimate now identify the same analytic coefficient. The following statement retains the distinct bulk and cap twists. \begin{proposition}[The cap-selected coefficient]\label{prop:physical-coefficient} Fix \(N\ge N_0\), homogeneous sites, and \(0\le\theta\le\pi/2\). Evaluate \(W,J\) at the real roots of Proposition~\ref{prop:real-counting} with \(\delta=\theta/(2\pi)\), and write the quotient as \((W/J)_\theta\). Put \(F_\theta(u)=F(u;\mathbf0,e^{i\theta})\). In the cap matrix element with bulk spin twist \(K^{-1}=e^{-i\theta}\), right cap twist \(i\), and left cap twist \(-i\), the eigenvalue selected in Proposition~\ref{prop:fixed-N-cap-estimate} is the TQ value \(F_\theta(l)\). Its coefficient is \begin{equation}\label{eq:physical-coefficient} c_{\rm phys}(\theta)=\frac{(W/J)_\theta}{(W/J)_{\pi/2}} . \end{equation} In particular \(c_{\rm phys}(\pi/2)=1\). \end{proposition} \begin{proof} Start near the vacuum with a row argument \(l\beta\). For generic parameters there, Lemma~\ref{lem:ordinary-physical-term} gives a nonzero isolated term with base \(F(u;\mathbf b,K)\). The generating function of an exponential polynomial has a pole at the reciprocal of each uncancelled nonzero base. The cap generating function has no such pole in this disk except possibly at \(\lambda^{-1}\). Therefore \(F(u;\mathbf b,K)=\lambda(\mathbf b,u,K)\), and comparison of the simple-pole coefficients gives \[ c_{\rm phys}(\mathbf b,u,K) =\frac{N!a(\mathbf b,K)}{\tau(1)} \] at these generic parameters. Take the limit of this equality to \(K=i\) at generic nearby sites. The cap coefficient tends to \(1\), and \(a(\mathbf b,K)\) tends to the nonzero value \(a(\mathbf b,i)\). Since \(\tau(1)\) is independent of \(K\), the common normalization is thereby identified, giving \[ c_{\rm phys}(\mathbf b,u,K) =\frac{a(\mathbf b,K)}{a(\mathbf b,i)} \] near the vacuum by analytic continuation. This limiting argument uses the generic trace identity only before taking the limit. No representation by jets supported only in the root torus is asserted at \(K=i\) or \(K=1\), which can be toric-end resonances. Now restrict the sites to the homogeneous value and continue along \(K=e^{i\theta}\) from \(i\) to \(1\). Proposition~\ref{prop:real-counting} supplies locally analytic roots and nonzero \(J\), and the preceding observation gives \(W\ne0\), all along this interval. Proposition~\ref{prop:fixed-N-cap-estimate} supplies locally analytic annular eigenvalues and cap coefficients there; the simple positive eigenvalue selects the same branch on overlaps. The identity theorem therefore continues both \(F=\lambda\) and the coefficient ratio through the interval. Continuity in the row argument includes \(u=l\), where \(P(u)/d_{\rm full}(u)\) is regular because \(d_{\rm full}(l)=([4][3])^N\ne0\). This proves \eqref{eq:physical-coefficient}. \end{proof} The next section estimates this quotient uniformly as \(N\to\infty\). The annular gaps and analytic neighborhoods used here were established separately for each fixed \(N\); they supply the exact coefficient identity and make no uniform claim about those gaps. ==== END Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//physical-amplitude.tex ==== ==== BEGIN Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//projection.tex ==== \section{Height moments and planar nesting}\label{sec:projection} The cap coefficient has growth $N^{1/6+o(1)}$ at winding fugacity $2$. To extract a planar nesting exponent, we must control the height of polygons before projecting the cylinder. Arbitrary fixed twist derivatives provide the required height moments. \input{sections/higher-jets} \subsection{From derivatives to polygon height moments} We now use the derivative bounds in positive path sums at the physical row $t=l$. Throughout this subsection a polygon $P$ has \emph{bare weight} $w(P)=x^{|P|}$, with no loop fugacity included. An open path still has $w(\gamma)=x^{m(\gamma)}$. Thus a winding polygon in the gas with fugacity $\nu$ has activity $\nu w(P)$; in particular, $P_N$ is the bare mass from Proposition~\ref{prop:cylinder}, not the mass at fugacity $2$. For a path or polygon let $r_-$ and $r_+$ be its lowest and highest occupied rhombus rows. Write \[ \operatorname{top}(P)=r_+(P),\qquad \operatorname{ht}(P)=r_+(P)-r_-(P)+1. \] For an arc below the cut at row $0$, put $\operatorname{dep}(\gamma)=-r_-(\gamma)$. Changing between row cuts and centre heights costs only an absolute number of rows. Constants denoted by $C_{\rm lat}$ below allow for this fixed lattice convention, and $u_+=\max\{u,0\}$. \paragraph{The winding pair coefficient.} Let $\mathscr W_{N,M}$ be the finite set of simple winding polygons whose occupied rows are among $-M,\ldots,-1$. Two polygons are incompatible if they share a dual vertex; two copies of the same polygon are also declared incompatible. Define the finite gas partition \[ \Xi_{N,M}(\nu)= \sum_{\substack{\mathcal B\subset\mathscr W_{N,M}\\ \mathcal B\text{ pairwise compatible}}} \prod_{\eta\in\mathcal B}\nu w(\eta). \] The empty set contributes one. Expanding its logarithm at $\nu=0$ gives the exact finite identity \begin{equation}\label{eq:winding-pair-finite} -2[\nu^2]\log\Xi_{N,M}(\nu) =\sum_{\gamma,\eta\in\mathscr W_{N,M}} w(\gamma)w(\eta)\mathbf1_{\{\gamma\not\sim\eta\}}. \end{equation} The sum is ordered and includes repetitions. In transfer notation, $\Xi_{N,M}(\nu)=\langle0|(\mathsf T_\nu(l))^M|0\rangle$. For each fixed $N$, the analytic gap at $\nu=0$ gives a nonzero analytic leading coefficient times $F_N(-l,\nu)^M$, with an analytic exponentially smaller remainder. Consequently \[ \lim_{M\to\infty}\frac1M[\nu^2]\log\Xi_{N,M}(\nu) =[\nu^2]\log F_N(-l,\nu). \] This limit is taken for fixed $N$; it requires no uniform analyticity radius in $\nu$. Let $\mathscr W_N$ be the winding polygons at all row positions on the infinite cylinder and let $\mathscr W_N^0$ be those with top row $0$. Horizontal translates are distinct members when their vertex sets are distinct, and $P_N=\sum_{\eta\in\mathscr W_N^0}w(\eta)$. For an anchored $\gamma\in\mathscr W_N^0$ and a placed $\eta\in\mathscr W_N$, let $r(\gamma,\eta)$ be the number of rows from the lowest to the highest occupied row of their union. Translating the pair vertically puts it in the $M$-row slab in exactly $(M-r(\gamma,\eta)+1)_+$ ways. Hence the right side of \eqref{eq:winding-pair-finite}, divided by $M$, is \[ \sum_{\gamma\in\mathscr W_N^0}\sum_{\eta\in\mathscr W_N} w(\gamma)w(\eta)\mathbf1_{\{\gamma\not\sim\eta\}} \frac{(M-r(\gamma,\eta)+1)_+}{M}. \] Each placement fraction increases to one. Monotone convergence and Lemma~\ref{lem:v-higher-twist-jets} therefore give the finite mass \begin{equation}\label{eq:winding-pair-mass} \begin{split} Q_N&:=\sum_{\gamma\in\mathscr W_N^0}w(\gamma) \sum_{\eta\in\mathscr W_N} w(\eta)\mathbf1_{\{\gamma\not\sim\eta\}}\\ &=-2[\nu^2]\log F_N(-l,\nu)\le \frac{C}{N}. \end{split} \end{equation} We relate this incompatibility mass to height. A simple winding curve separates the cylinder into an upper and a lower component. Write $u_\eta(z)$ for the indicator that a lattice centre $z$ lies strictly in the upper component of $\eta$. Let $\tau_{a,b}$ translate by $a+b e^{i\pi/3}$, where $a\in\mathbb Z/N\mathbb Z$ and $b\in\mathbb Z$. On a fixed $\gamma\in\mathscr W_N^0$, choose centres $z_+,z_-$ of the same triangle type within $C_{\rm lat}$ rows of its two extreme rows. Their row-height difference $\delta_\gamma$ satisfies $\delta_\gamma\ge\operatorname{ht}(\gamma)-C_{\rm lat}$. Since the two centres differ by one of the translations $\tau_{a,b}$, reindexing $a$ and telescoping in $b$ gives, for every winding shape $\eta$, \begin{equation}\label{eq:winding-indicator-telescope} \begin{split} N\delta_\gamma &=\sum_{a=0}^{N-1}\sum_{b\in\mathbb Z} \bigl(u_{\tau_{a,b}\eta}(z_+) -u_{\tau_{a,b}\eta}(z_-)\bigr)\\ &\le\sum_{a=0}^{N-1}\sum_{b\in\mathbb Z} \mathbf1_{\{\gamma\not\sim\tau_{a,b}\eta\}}. \end{split} \end{equation} The differences vanish for all sufficiently high or low translates, so the first sum is finite. For the inequality, if $\eta$ avoids the connected curve $\gamma$, the two centres on $\gamma$ lie in the same complementary component and the two indicators agree. Multiply \eqref{eq:winding-indicator-telescope} by $w(\eta)$ and sum over $\eta\in\mathscr W_N^0$. Every placed polygon on the right is counted exactly $N$ times: summing both the anchored shapes and all $N$ horizontal translations has this multiplicity even when a shape has a horizontal stabilizer. Dividing by $N$ yields \[ P_N\delta_\gamma\le \sum_{\eta\in\mathscr W_N}w(\eta) \mathbf1_{\{\gamma\not\sim\eta\}}. \] Summing over $\gamma$ and using $P_N\asymp N^{-1}$ in \eqref{eq:winding-pair-mass} proves \begin{equation}\label{eq:winding-height} \frac1{P_N}\sum_{\gamma\in\mathscr W_N^0} w(\gamma)\operatorname{ht}(\gamma) \le C_{\rm lat}+\frac{Q_N}{P_N^2}\le C N. \end{equation} \paragraph{Normalized arc entries.} For $1\le gL}}w(P)\\ &\quad\le \sum_{P\in\mathscr C_N^0}w(P) (\operatorname{ht}(P)+C_{\rm lat}) \mathbf1_{\{\operatorname{ht}(P)>L\}}\\ &\quad\le L^{-k}\sum_{P\in\mathscr C_N^0}w(P) (\operatorname{ht}(P)+C_{\rm lat})^{k+1} \le C_kN^{k+2}L^{-k}. \end{aligned} \end{equation} In particular the bare mass with height greater than $N^{1+\epsilon}$ is at most $C_kN^{2-k\epsilon}$. For clarity, a contractible cylinder polygon enclosing $v$ has a unique planar lift whose disk contains a specified lift of $v$. Indeed the boundaries of two translates of a lifted disk are disjoint. Their disks must then be disjoint or nested; strict nesting is impossible for two translates of the same bounded disk, since translation preserves area. Thus their interiors are disjoint. We call the lift containing the specified tip the centred lift and denote it by $\widetilde P$. Fix $\epsilon>0$ and put $M=\lfloor N^{1+2\epsilon}\rfloor$. Consider centred lifts with old height at most $N^{1+\epsilon}$ and diameter greater than $N^{1+3\epsilon}$. The translation $M e^{i\pi/3}$ changes the old row coordinate by $M$. For large $N$ this exceeds the old height, so the lift and every nonzero translate by this vector have disjoint old row ranges. It therefore projects to a simple contractible polygon on the cylinder of period $M$ in direction $\pi/3$. Write $\operatorname{ht}_0$ and $\operatorname{ht}_{\pi/3}$ for occupied-row heights in the original and rotated orientations. The two row coordinates are independent linear coordinates in the plane. There is an absolute $C$ such that \[ \diam(\widetilde P)\le C\bigl(\operatorname{ht}_0(\widetilde P) +\operatorname{ht}_{\pi/3}(\widetilde P)\bigr). \] The assumed diameter and old height thus force rotated height at least $c_1N^{1+3\epsilon}$ for large $N$. This projection is injective on the centred lifts: if two gave the same rotated-cylinder polygon, they would differ by a translate of its disk, while both disks contain the same specified tip lift. The same disjoint-disk argument on the rotated cylinder forces that translate to be zero. Applying \eqref{eq:tip-height-tail} after the $\pi/3$ lattice rotation gives \[ \begin{aligned} &\sum_{\substack{P\text{ enclosing }v,\ \operatorname{ht}_0(\widetilde P)\le N^{1+\epsilon}\\ \diam(\widetilde P)>N^{1+3\epsilon}}}w(P)\\ &\quad\le C_k M^{k+2}(c_1N^{1+3\epsilon})^{-k} \le C_kN^{2+4\epsilon-k\epsilon}. \end{aligned} \] Together with the old-height tail, this proves \begin{equation}\label{eq:tip-diameter-tail} \begin{split} q_{N,\epsilon}&:= \sup_v\sum_{\substack{P\text{ contractible, enclosing }v\\ \diam(\widetilde P)>N^{1+3\epsilon}}}w(P)\\ &\le C_k\bigl(N^{2-k\epsilon}+N^{2+4\epsilon-k\epsilon}\bigr)=o(1), \end{split} \end{equation} where one fixes $k>(2+4\epsilon)/\epsilon$ and then lets $N\to\infty$. This is a bound on total bare activity, not merely on the probability of a large polygon in a gas. \subsection{The two seam tips and the planar comparison} The height estimates permit the projection to the plane. We first identify the cap matrix element as a positive two-tip partition and control the cost of horizontal slab walls. At homogeneous sites let $G_H=\Lambda(T^1(l))^H\Psi$, where the superscript $1$ denotes trivial spin twist. The cap regularity argument in Section~\ref{sec:physical-amplitude} represents the lower and upper caps as limits of empty boundaries extended by physical rows with twists $i$ and $-i$, respectively; the upper construction uses transpose, reversal, and spin negation. In a finite such approximation the seam twist is $i$ below the tip at cut $0$, $1$ between cuts $0$ and $H$, and $-i$ above the tip at $H$, as indicated in Figure~\ref{fig:seam}. For an oriented loop let $c_t$ be its signed crossing number of the upward seam ray from tip $t$, with left-going auxiliary flow positive. For a counterclockwise contractible loop, $c_t$ is one if its disk contains $t$ and zero otherwise. Its turn phase is $i$, and its total full-seam crossing number is zero. For a winding loop oriented to the left, $c_t$ is one if $t$ lies in its lower complementary component and zero otherwise; its turn phase is one and its full-seam phase for the constant twist $i$ is $i$. In either case the two successive twist ratios are $1/i=(-i)/1=-i$. The sum of the two orientations is therefore \[ 2\Re\bigl(i(-i)^{c_0+c_H}\bigr) =\begin{cases} 2,&c_0+c_H=1,\\ 0,&c_0+c_H=0\text{ or }2. \end{cases} \] Thus precisely the contractible loops enclosing one tip and the windings separating the tips survive, each with activity $2x^{|P|}$. Every finite cap approximation is the corresponding positive gas in its finite cylinder. Exhaustion proves $G_H=\mathcal G_N(H)$. Let $\mathfrak b_N$ be the total bare weight of windings meeting a fixed horizontal cut. Counting their possible row placements and using \eqref{eq:winding-height} gives the uniform wall bound \begin{equation}\label{eq:winding-wall-mass} \mathfrak b_N\le\sum_{\eta\in\mathscr W_N^0}w(\eta) (\operatorname{ht}(\eta)+C_{\rm lat}) \le CNP_N\le C. \end{equation} The winding-only slab partition in the introduction is $L_N(h)=\Xi_{N,h}(2)$, after translating the slab. Stacking slabs and retaining only polygons internal to each gives the lower inequality below. For the upper inequality, keep the exclusions within the two slabs and discard all exclusions involving polygons that meet their common wall. Those polygons have total activity at most $2\mathfrak b_N$, so their unrestricted gas is at most $\prod_\eta(1+2w(\eta))\le e^{2\mathfrak b_N}$. Hence \begin{equation}\label{eq:slab-wall-inequalities} L_N(h)L_N(j)\le L_N(h+j) \le e^{2\mathfrak b_N}L_N(h)L_N(j). \end{equation} The positive annular transfer at fugacity $2$ is primitive by Section~\ref{sec:physical-amplitude}. Its dominant eigenvalue $\lambda_N$ is therefore the exponential growth rate of its empty-to-empty entry $L_N(h)$. Iterate \eqref{eq:slab-wall-inequalities} on multiples of $h$ and take the growth rate to obtain \begin{equation}\label{eq:slab-growth-bounds} e^{-2\mathfrak b_N}\lambda_N^h\le L_N(h)\le\lambda_N^h, \qquad 0\le\log\lambda_N\le 2P_N. \end{equation} For the last inequality, discarding all exclusions gives $\log L_N(h)\le2\sum_{\eta\in\mathscr W_{N,h}}w(\eta)\le2hP_N$. In particular, for $H>2N$, \begin{equation}\label{eq:inward-wall-cost} \frac{L_N(H-2N)}{\lambda_N^H} \ge e^{-2\mathfrak b_N}\lambda_N^{-2N} \ge e^{-2\mathfrak b_N-4NP_N}\ge e^{-C}. \end{equation} Write $Z(r)=\mathcal Z_v(r)$ for the planar nesting partition at either tip. Choose a small absolute $a_0>0$ so that a disk of radius $a_0N$ projects injectively to the period-$N$ cylinder and, when centred at the lower tip, lies below the cut at $N$; the analogous disk at the upper tip lies above the cut at $H-N$. For large $H$ the two disks are disjoint. Their nests coexist with every winding configuration internal to the reduced slab $[N,H-N]$, and all those windings separate the tips. Thus \[ \mathcal G_N(H)\ge Z(a_0N)^2L_N(H-2N). \] For the upper bound use the original slab $[0,H]$. A separating winding not internal to this slab must meet its lower or upper wall, so such windings have total activity at most $4\mathfrak b_N$. A contractible polygon with centred lift of diameter at most $D=N^{1+3\epsilon}$ lies within distance $D$ of its enclosed tip: the tip lies in the convex hull of the polygon. The centred lifts at each tip form a planar nest, whose partition is bounded by $Z(D)$. Finally, \eqref{eq:tip-diameter-tail} bounds the total activity of the remaining contractible polygons at both tips by $4q_{N,\epsilon}$. Retain the exclusions within each of the two small nests and within the slab-internal winding gas, and discard all other exclusions. In finite exhaustions this gives, and hence in the limit gives, \[ \mathcal G_N(H)\le L_N(H)Z(N^{1+3\epsilon})^2 \exp(4\mathfrak b_N+4q_{N,\epsilon}). \] These bounds in particular make the positive full-cylinder sum finite. Equations \eqref{eq:slab-growth-bounds}--\eqref{eq:inward-wall-cost} now give the comparison \begin{equation}\label{eq:planar-amplitude-squeeze} e^{-C}Z(a_0N)^2\le \lim_{H\to\infty}\lambda_N^{-H}\mathcal G_N(H) \le e^{C_\epsilon}Z(N^{1+3\epsilon})^2. \end{equation} The limit exists by the physical cap expansion and equals $c_{\rm phys}(0)=N^{1/6+o(1)}$. The left inequality in \eqref{eq:planar-amplitude-squeeze}, with monotonicity between successive $a_0N$, gives $\limsup_{r\to\infty}\log Z(r)/\log r\le1/12$. The right inequality, with $N=\lfloor r^{1/(1+3\epsilon)}\rfloor$, gives $\liminf_{r\to\infty}\log Z(r)/\log r\ge1/(12(1+3\epsilon))$. Letting $\epsilon\downarrow0$ proves $Z(r)=r^{1/12+o(1)}$ and Theorem~\ref{thm:amplitude}. ==== END Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//projection.tex ==== ==== BEGIN Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//references.tex ==== \bibitem{BatchelorYung1995} M.~T.~Batchelor and C.~M.~Yung. Exact results for the adsorption of a flexible self-avoiding polymer chain in two dimensions. \emph{Physical Review Letters} \textbf{74} (1995), 2026--2029. \href{https://doi.org/10.1103/PhysRevLett.74.2026} {doi:10.1103/PhysRevLett.74.2026}. \bibitem{DeConciniProcesi1995} C.~De Concini and C.~Procesi. Wonderful models of subspace arrangements. \emph{Selecta Mathematica, New Series} \textbf{1} (1995), 459--494. \href{https://doi.org/10.1007/BF01589496}{doi:10.1007/BF01589496}. \bibitem{DiFrancescoZinnJustin2005} P.~Di Francesco and P.~Zinn-Justin. Around the Razumov--Stroganov conjecture: proof of a multi-parameter sum rule. \emph{Electronic Journal of Combinatorics} \textbf{12} (2005), Research Paper~6. \href{https://doi.org/10.37236/1903}{doi:10.37236/1903}. \href{https://arxiv.org/abs/math-ph/0410061}{arXiv:math-ph/0410061}. \bibitem{DuminilCopinSmirnov2012} H.~Duminil-Copin and S.~Smirnov. The connective constant of the honeycomb lattice equals $\sqrt{2+\sqrt2}$. \emph{Annals of Mathematics} \textbf{175} (2012), 1653--1665. \href{https://doi.org/10.4007/annals.2012.175.3.14}{doi:10.4007/annals.2012.175.3.14}. \bibitem{FernandezProcacci2007} R.~Fern\'andez and A.~Procacci. Cluster expansion for abstract polymer models. New bounds from an old approach. \emph{Communications in Mathematical Physics} \textbf{274} (2007), 123--140. \href{https://doi.org/10.1007/s00220-007-0279-2}{doi:10.1007/s00220-007-0279-2}. \bibitem{FultonMacPherson1994} W.~Fulton and R.~MacPherson. A compactification of configuration spaces. \emph{Annals of Mathematics} \textbf{139} (1994), 183--225. \href{https://doi.org/10.2307/2946631}{doi:10.2307/2946631}. \bibitem{GarbaliNienhuis2017Ground} A.~Garbali and B.~Nienhuis. The dilute Temperley--Lieb $O(n=1)$ loop model on a semi infinite strip: the ground state. \emph{Journal of Statistical Mechanics: Theory and Experiment} (2017), 043108. \href{https://doi.org/10.1088/1742-5468/aa6a30}{doi:10.1088/1742-5468/aa6a30}. \href{https://arxiv.org/abs/1411.7020}{arXiv:1411.7020}. \bibitem{GarbaliNienhuis2017Sum} A.~Garbali and B.~Nienhuis. The dilute Temperley--Lieb $O(n=1)$ loop model on a semi infinite strip: the sum rule. \emph{Journal of Statistical Mechanics: Theory and Experiment} (2017), 053102. \href{https://doi.org/10.1088/1742-5468/aa6bc3}{doi:10.1088/1742-5468/aa6bc3}. \href{https://arxiv.org/abs/1411.7160}{arXiv:1411.7160}. \bibitem{GlazmanManolescu2017} A.~Glazman and I.~Manolescu. Self-avoiding walk on $\mathbb Z^2$ with Yang--Baxter weights: universality of critical fugacity and 2-point function. \emph{Annales de l'Institut Henri Poincar\'e, Probabilit\'es et Statistiques} \textbf{56} (2020), 2281--2300. \href{https://doi.org/10.1214/19-AIHP1024}{doi:10.1214/19-AIHP1024}. \bibitem{Griffiths1976Abel} P.~A.~Griffiths. Variations on a Theorem of Abel. \emph{Inventiones Mathematicae} \textbf{35} (1976), 321--390. \href{https://doi.org/10.1007/BF01390145}{doi:10.1007/BF01390145}. \bibitem{GrimmPearce1993} U.~Grimm and P.~A.~Pearce. Multi-colour braid--monoid algebras. \emph{Journal of Physics A: Mathematical and General} \textbf{26} (1993), 7435--7460. \href{https://doi.org/10.1088/0305-4470/26/24/018}{doi:10.1088/0305-4470/26/24/018}. \href{https://arxiv.org/abs/hep-th/9303161}{arXiv:hep-th/9303161}. \bibitem{HammersleyWelsh1962} J.~M.~Hammersley and D.~J.~A.~Welsh. Further results on the rate of convergence to the connective constant of the hypercubical lattice. \emph{Quarterly Journal of Mathematics} \textbf{13} (1962), 108--110. \href{https://doi.org/10.1093/qmath/13.1.108}{doi:10.1093/qmath/13.1.108}. \bibitem{KrachunPanagiotis2026} D.~Krachun and C.~Panagiotis. Quantitative sub-ballisticity of self-avoiding walk on the hexagonal lattice. \emph{Annals of Probability} \textbf{54} (2026), 1109--1125. \href{https://doi.org/10.1214/24-AOP1730}{doi:10.1214/24-AOP1730}. \href{https://arxiv.org/abs/2310.17299}{arXiv:2310.17299}. \bibitem{Li2009Wonderful} L.~Li. Wonderful compactification of an arrangement of subvarieties. \emph{Michigan Mathematical Journal} \textbf{58} (2009), 535--563. \href{https://doi.org/10.1307/mmj/1250169076}{doi:10.1307/mmj/1250169076}. \bibitem{Nienhuis1982} B.~Nienhuis. Exact critical point and critical exponents of $O(n)$ models in two dimensions. \emph{Physical Review Letters} \textbf{49} (1982), 1062--1065. \href{https://doi.org/10.1103/PhysRevLett.49.1062}{doi:10.1103/PhysRevLett.49.1062}. \bibitem{Penrose1967} O.~Penrose. Convergence of fugacity expansions for classical systems. In T.~A.~Bak (ed.), \emph{Statistical Mechanics: Foundations and Applications}, W.~A.~Benjamin, New York--Amsterdam, 1967, 101--109. \bibitem{Schoenberg1950PFII} I.~J.~Schoenberg. On P\'olya frequency functions. II: Variation-diminishing integral operators of the convolution type. \emph{Acta Scientiarum Mathematicarum (Szeged)} \textbf{12} (1950), 97--106. \href{https://acta.bibl.u-szeged.hu/13626/1/math_012_pars_b_097-106.pdf}{Original journal article}. \bibitem{ScottSokal2005} A.~D.~Scott and A.~D.~Sokal. The repulsive lattice gas, the independent-set polynomial, and the Lov\'asz local lemma. \emph{Journal of Statistical Physics} \textbf{118} (2005), 1151--1261. \href{https://doi.org/10.1007/s10955-004-2055-4}{doi:10.1007/s10955-004-2055-4}. \bibitem{ZhouBatchelor1997} Y.~K.~Zhou and M.~T.~Batchelor. Critical behaviour of the dilute $O(n)$, Izergin--Korepin and dilute $A_L$ face models: Bulk properties. \emph{Nuclear Physics B} \textbf{485} (1997), 646--664. \href{https://doi.org/10.1016/S0550-3213(96)00654-2}{doi:10.1016/S0550-3213(96)00654-2}. \href{https://arxiv.org/abs/cond-mat/9611156}{arXiv:cond-mat/9611156}. ==== END Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//references.tex ==== ==== BEGIN Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//repeated-moments.tex ==== \section{Arbitrary and repeated row arguments}\label{sec:repeated-moments} To obtain repeated rows at arbitrary arguments, add a fused pair of sites for each desired argument and remove the pairs after applying the distinct-site formula. The following collision calculation justifies that removal for the scalar trace. \begin{theorem}[All row moments between the vacuum caps]\label{thm:row-moments} Let $N,h\ge0$ be integers. With ordered site variables $\mathbf b=(b_1,\ldots,b_N)$, the right and left vacuum caps $\Psi,\Lambda$, the scalar trace $\tau$, and the Laurent polynomial part $P$ of Sections~\ref{sec:vacuum}--\ref{sec:site-moments}, the following is an identity of meromorphic functions of the site variables, nonzero twist $K$, and row variables $u_1,\ldots,u_h$: \begin{equation}\label{eq:all-row-moments} \Lambda\prod_{a=1}^h T^{K^{-1}}(u_a;\mathbf b)\Psi =\frac{\tau\!\left(\prod_{a=1}^h P(u_a)/d_{\rm full}(u_a)\right)}{\tau(1)}. \end{equation} At a common-twist collision, the trace is the holomorphic specialization near a generic diagonal supplied by Proposition~\ref{prop:common-twist}; colliding roots are retained by that continuation. Pointwise evaluation is at regular values of the meromorphic expressions, including repeated row arguments. For $N=0$ the empty site and root conventions apply, and for $h=0$ both sides are one. \end{theorem} We prove the theorem by a two-site reduction. The matrix side reduces by fusion. On the scalar side, we show that both traces have at most a simple pole, with residues equal to a common nonzero factor times the corresponding smaller traces. The denominator residue is nonzero generically, so that factor cancels in the normalized ratio. \subsection{A scalar two-site reduction} Let $\mathbf b_{\rm red}$ be the ordered list of $n\ge0$ sites obtained by deleting two neighboring sites $x,y$, in that order, from a larger list $\mathbf b_\epsilon$, where $x=y+3+\epsilon$. All other sites are generic. Use subscripts $\epsilon$ and ${\rm red}$ on $\tau,P,d_{\rm full},W,D_i,J$ to distinguish the two lists; a subscript $0$ on a regular coefficient denotes its specialization at $\epsilon=0$. \begin{lemma}[Two-site reduction of the normalized trace]\label{lem:scalar-two-site-fusion} Fix a generic common twist for both lists. For any finite list of row arguments $u_a(\epsilon)$ with limits $u_a(0)$, assume that for every $a$ both limiting denominator evaluations are regular and nonzero: \[ d_{{\rm full},0}(u_a(0))\ne0, \qquad d_{{\rm full},{\rm red}}(u_a(0))\ne0. \] After shrinking the $\epsilon$-disk, the nearby evaluations used below are then regular and nonzero. Put \[ f_\epsilon=\prod_a\frac{P_\epsilon(u_a(\epsilon))} {d_{{\rm full},\epsilon}(u_a(\epsilon))}, \qquad f_{\rm red}=\prod_a\frac{P_{\rm red}(u_a(0))} {d_{{\rm full},{\rm red}}(u_a(0))}. \] Then, at generic values of the remaining parameters, \begin{equation}\label{eq:scalar-two-site-fusion} \lim_{\epsilon\to0}\frac{\tau_\epsilon(f_\epsilon)}{\tau_\epsilon(1)} =\frac{\tau_{\rm red}(f_{\rm red})}{\tau_{\rm red}(1)}. \end{equation} The identity extends meromorphically in the remaining parameters. Empty products are allowed. The row arguments may vary with $\epsilon$ subject to the stated regularity. \end{lemma} \begin{proof} We first bound the order at the merging sites and identify the contributing root cluster. We then compute its residue and restrict the polynomial insertions. \paragraph{Preparing the site collision.} Keep the \(n+2\) twists \(K_i\) independent during the transverse calculation, while the polynomial insertion still uses the single variable \(K\). We calculate residues coefficientwise for a Laurent-polynomial insertion \(F_\epsilon\) in the roots whose coefficients are rational in the exponential parameter coordinates and jointly holomorphic in the site variables (including \(\epsilon\)), the twists \(K,K_i\), and the row variables near the chosen collision point and generic diagonal, after inverting the stated nonzero scalar denominators. Both \(F_\epsilon=1\) and \(F_\epsilon=f_\epsilon\) have this form. The calculation will show that \(\epsilon\tau_\epsilon(F_\epsilon)\) has a regular limit and will express that limit as a sum of smaller independent-twist traces. We keep the row variables independent in this calculation and substitute their prescribed limits at the end. The multiplicative site coordinates remain in a fixed compact annulus and have distinct shift orbits for all sufficiently small \(\epsilon\ne0\). We may therefore use the uniform punctured-neighborhood conclusion in the proof of Proposition~\ref{prop:common-twist}: one neighborhood \(U\) of a generic diagonal, avoiding the subset-product toric-end hypersurfaces for the large trace and for every smaller trace obtained by deleting two twist coordinates, works throughout this punctured site family. With the remaining site and row parameters near their chosen regular values, \(\tau_\epsilon(F_\epsilon)\) is holomorphic there, so no other polar divisor can approach the collision inside \(U\). Once at-most-simple order is proved for generic independent twists, \(\epsilon\tau_\epsilon(F_\epsilon)\) is therefore regular jointly near the generic diagonal: one may use holomorphic extension off codimension two, or cancellation of every possible denominator divisor in the rational function. This will carry the independent-twist residue identity to the diagonal. The transverse calculation itself is made off all other exceptional sets of twists on \(\epsilon=0\). \paragraph{The only contributing cluster.} Use the compactification argument of Lemma~\ref{lem:residue} over a small disk of site values, counting moving \(x\) as one more configuration coordinate relative to fixed anchors; end limits remain excluded. Compute order of \(W\prod dt_i\wedge d\epsilon\) relative to log poles. Single-site and unanchored depths give the earlier bounds for bare \(W\) (a cluster drifting with the moving site but not \(y\) counts as anchored), whereas until separation of \(x,y\) the cluster containing both gives \[ B-m_{-1}-m_1-m_2-m_4+1 \] with offsets relative to $y$, since the log discrepancy is now one larger. The following claim gives both the order bound and its possible equality case. \begin{claim}[Order at the merging sites]\label{claim:merging-site-order} For the offset counts of a cluster containing both sites, put \[ B=\frac12\sum_{k\bmod8}(m_k-m_{k+1})^2 +\sum_{k\bmod8}m_km_{k+2}. \] Then \[ B-m_{-1}-m_1-m_2-m_4+1\ge0, \] and equality requires unit support $\{1,2\}$. \end{claim} \begin{proof} If the minimum count is positive, the distance-two products already suffice; one can use the terms starting at each of the four anchored indices. Suppose henceforth that at least one count is zero, and write \[ a=m_{-1},\quad b=m_1,\quad c=m_2,\quad f=m_4, \qquad M=\max_k m_k. \] \emph{A doubly occupied anchored pair.} If either pair $(a,b)$ or $(c,f)$ has both counts positive, use $B\ge M+ab+cf$. When $a,b>0$, the inequality $(a-1)(b-1)\ge0$ gives $a+b\le ab+1$, and analogously $c+f\le cf+1$ when $c,f>0$. A singly occupied pair contributes at most $M$. If both pairs are doubly occupied, use also $2\le M+1$. In every case this gives \[ a+b+c+f\le M+ab+cf+1\le B+1; \] the first inequality is strict when the other pair is empty. Equality thus requires the other pair to be occupied too, equality in the variation bound, and no distance-two products beyond $ab,cf$. In particular the support must be contiguous. This is impossible: contiguous support containing the two indices of the first doubly positive pair and an index of the other either gives an additional unwanted product, or makes both pairs doubly occupied with unwanted products. \emph{Neither anchored pair is doubly occupied.} Use $B\ge2M-1$; the empty case is direct. The proof of that bound in Proposition~\ref{prop:common-twist} gives strictness for $M\ge3$. At $M=3$, the four gradients give an integral bound that could be sharp only if both first neighbors were positive, adding their product. If a second neighbor is positive, the variation and its product give at least $2M$ altogether. For $M=2$, equality only needs to be excluded when two peaks have height $2$. All their second neighbors must vanish. The variation-square mass on the four edges around one peak is already at least $2$, leaving at most $1$; hence no other height-$2$ peak can lie outside its first neighbors, since the remaining path starts and ends at zero. If one first neighbor equals $2$, the variation there is already at least $3$ and the other first neighbor must be unit, adding their product. Thus the two distinct anchored pairs cannot both supply height $2$. For $M=1$, the unit equality sets for $B=M$ give precisely the asserted pair $\{1,2\}$. This proves the claim. \end{proof} Over generic independent twists all levels are transverse on boundary strata and solutions over \(\epsilon=0\) lie on only one boundary component (codimension at most one). A nonzero log order cannot contribute to the limiting trace of \(\epsilon W f\) (residues divided by log Jacobian understood): locally use level values as coordinates with one further coordinate \(v\), \(\epsilon=v^m\) after absorbing a unit. Neutrality and log order zero are thus necessary. Clusters off the merging orbit are strictly excluded by the argument of common-twist regularity for \(W\). In the merging orbit equality requires exactly two distinct labels \(A,B\) at offsets 1,2 at every depth containing both sites. In the local boundary coordinate this implies \[ \operatorname{ord}(t_A-y-1),\ \operatorname{ord}(t_B-y-2) \ge \operatorname{ord}\epsilon, \] so \(t_A=y+1+\epsilon s,t_B=y+2+\epsilon u\) with bounded \(s,u\). The scaled equations will decide whether any further singular collision occurs in these coordinates. \paragraph{The two scaled root equations.} Since \(S(z)S(z-3)R(z-1)R(z-2)=1\), each surviving equation tends to its equation in the regular smaller system, independently of \(s,u\). After collecting their regular nonzero factors, the two other equations have the cross-multiplied limits \[ C_A(s-u)=-K_A s,\qquad C_B(u-1)=-K_B(s-u) \] with nonzero constants depending on the surviving regular solution and the sites, but not on \(s,u\). These equations exclude the further singular collisions \(s=0\), \(u=1\), and \(s=u\): the first would force \(u=0\) and contradict the second equation; the second would force \(s=1\) and contradict the first; the third would force \(s=u=0\) and again contradict the second. For generic independent twists the equations have a unique finite simple solution. Together with a regular smaller solution it gives an actual branch by the implicit function theorem. At that solution the limiting level maps can be written \[ D_A^{(0)}=C_A\frac{s-u}{s},\qquad D_B^{(0)}=C_B\frac{u-1}{s-u}. \] \paragraph{The common scalar residue.} We now compute the factor contributed by one ordered choice of the deleted root labels \((A,B)\). In the ordinary angle coordinate \(\epsilon\), substitution into the scalar functions \(V,G\) of Lemma~\ref{lem:residue} gives \[ \lim_{\epsilon\to0}\epsilon^3s(u-1)(s-u) V(t_A-x)V(t_A-y)V(t_B-x)V(t_B-y)G(t_A-t_B) =\frac{[l]}{8[2l]^3[3l]}. \] Each remaining site \(b\) supplies \(V(y+1-b)V(y+2-b)\), whereas every surviving root factor cancels by \(G(z-1)G(z-2)V(z)V(z-3)=1\). Define \[ \kappa(\mathbf b_{\rm red};y) :=\frac{[l]}{8[2l]^3[3l]} \prod_{b\in\mathbf b_{\rm red}}V(y+1-b)V(y+2-b). \] This is a nonzero meromorphic function at generic sites. It may depend on the collision center \(y\) and every remaining site, but is independent of the surviving roots, all twists, and the ordered choice of deleted labels. The full weight consequently satisfies \[ \epsilon^3W_\epsilon\longrightarrow \frac{\kappa(\mathbf b_{\rm red};y)W_{\rm red}} {s(u-1)(s-u)}. \] The scaled Jacobian has exactly the same factor in \(s,u\). Indeed \[ \partial_{(s,u)}\begin{pmatrix}\log D_A^{(0)}\\ \log D_B^{(0)}\end{pmatrix} =\begin{pmatrix} (s-u)^{-1}-s^{-1} & -(s-u)^{-1}\\ -(s-u)^{-1} & (u-1)^{-1}+(s-u)^{-1} \end{pmatrix}, \qquad \det=\frac1{s(u-1)(s-u)}. \] The surviving level maps are independent of \(s,u\) in the limit, so the full derivative matrix in the scaled coordinates is block triangular, even though \(C_A,C_B\) may depend on the surviving roots. Since \(t_A,t_B\) each contribute a factor \(\epsilon\) to the coordinate change, and simultaneous row and column reorderings introduce no sign, \[ \epsilon^2J_\epsilon\longrightarrow \frac{J_{\rm red}}{s(u-1)(s-u)}, \qquad \epsilon\frac{W_\epsilon}{J_\epsilon}\longrightarrow \kappa(\mathbf b_{\rm red};y)\frac{W_{\rm red}}{J_{\rm red}}. \] For a coefficientwise limit \(F_0\), write \[ \iota_{A,B}^*F_0:= \left.F_0\right|_{t_A=y+1,\ t_B=y+2} \] for its restriction to the contributing root cluster. Write \(\tau_{\rm red}^{\widehat A,\widehat B}\) for the smaller trace with the surviving independent twists. Summing the branch limits gives, for generic independent twists, \[ \operatorname{Res}_{\epsilon=0}\tau_\epsilon(F_\epsilon) :=\lim_{\epsilon\to0}\epsilon\tau_\epsilon(F_\epsilon) =\kappa(\mathbf b_{\rm red};y) \sum_{A\ne B}\tau_{\rm red}^{\widehat A,\widehat B} (\iota_{A,B}^*F_0). \] The sum is over the \((n+2)(n+1)\) ordered choices of \((A,B)\). The simple-order and joint-regularity argument above continues this rational identity to a generic diagonal as an equality of these residue limits. For \(F_\epsilon=1\) or \(F_\epsilon=f_\epsilon\), the restrictions are symmetric in the surviving roots, so all summands become the same smaller trace after relabeling. Thus the common nonzero trace residue factor is \[ \mathcal C_{\rm del}(\mathbf b_{\rm red};y) :=(n+2)(n+1)\kappa(\mathbf b_{\rm red};y), \] and at a generic common twist \begin{equation}\label{eq:scalar-fusion-residue} \lim_{\epsilon\to0}\epsilon\tau_\epsilon(F_\epsilon) =\mathcal C_{\rm del}(\mathbf b_{\rm red};y)\, \tau_{\rm red}(\iota^*F_0), \qquad F_\epsilon=1\ \text{or}\ f_\epsilon. \end{equation} Here \(\iota^*\) denotes any ordered choice after relabeling the surviving roots. In particular, \eqref{eq:trace-normalization} gives the nonzero generic denominator residue \(\mathcal C_{\rm del}(\mathbf b_{\rm red};y)\tau_{\rm red}(1)\). \paragraph{Restricting the polynomial insertions.} It remains to identify \(\tau_{\rm red}(\iota^*f_0)\) in \eqref{eq:scalar-fusion-residue}. At the collision put \[ d_{\rm pair}(w):=\frac{d_{{\rm full},0}(w)}{d_{{\rm full},{\rm red}}(w)} =[w-y][w-y-1][w-y+3][w-y+2], \qquad Q_{\rm pair}(w):=[w-y-2][w-y-3]. \] Then \(\iota^*Q_0=Q_{\rm red}Q_{\rm pair}\). In every summand of \(L/d_{\rm full}\), the added \(Q_{\rm pair}\) factors cancel the two added sites using \([z+8]=-[z]\). Equivalently, as rational functions of the spectral variable, \[ \iota^*L_0(w)=d_{\rm pair}(w)L_{\rm red}(w). \] For fixed generic remaining sites and common twist, let \(\mathcal A_{\rm red}=\mathbb C[p_i^{\pm1}:i\ne A,B]\), where \(p_i=e^{2it_i}\). Choose root and spectral monomials so that \(C_{\rm red}(z)\in\mathcal A_{\rm red}[z]\) is the polynomial form of \(Q_{\rm red}(w)Q_{\rm red}(w-1)\), with unit constant and leading coefficients. Set \[ \mathcal R_{\rm red}(z):=C_{\rm red}(z)(L_{\rm red}(w)-P_{\rm red}(w)) =\sum_j r_jz^j\in\mathcal A_{\rm red}[z^{\pm1}], \qquad I_{\rm red}:=(r_j:j\in\mathbb Z) \subset\mathcal A_{\rm red}. \] Multiplying \(\mathcal R_{\rm red}\) by a spectral monomial makes it a cleared polynomial of the kind used in Proposition~\ref{prop:common-twist} and does not change its coefficient ideal. Let \(\operatorname{pp}_z\) denote the Laurent polynomial part used to define \(P\). In the formal divisions defining \(P_\epsilon\), the leading and trailing coefficients of the root denominators are monomial units on the torus, so division and hence taking \(\operatorname{pp}_z\) commute with \(\iota^*\). Since \(d_{\rm pair}P_{\rm red}\) is already a Laurent polynomial, the exact restriction needed here is \begin{equation}\label{eq:scalar-fusion-polynomial-ideal} \iota^*P_0(w)-d_{\rm pair}(w)P_{\rm red}(w) =\operatorname{pp}_z\!\left(d_{\rm pair}(w) \frac{\mathcal R_{\rm red}(z)}{C_{\rm red}(z)}\right) \in I_{\rm red}\,\mathcal A_{\rm red}[z^{\pm1}]. \end{equation} Indeed formal division at zero and infinity expresses every coefficient of this finite polynomial part as an \(\mathcal A_{\rm red}\)-linear combination of the coefficients of \(\mathcal R_{\rm red}\). The divisions use the coefficients of \(C_{\rm red}\) and inverses only of its constant and leading root-monomial units. At each allowed row value, divide \eqref{eq:scalar-fusion-polynomial-ideal} by the nonzero scalar \(d_{\rm pair}(u_a(0))d_{{\rm full},{\rm red}}(u_a(0))\). Each normalized insertion difference then belongs to \(I_{\rm red}\), after adjoining the row-parameter scalars. Telescoping a product of such factors, including the empty product, gives \begin{equation}\label{eq:scalar-fusion-product-ideal} \iota^*f_0-f_{\rm red}\in I_{\rm red}. \end{equation} Proposition~\ref{prop:common-twist} says that \(\tau_{\rm red}\) annihilates every Laurent-polynomial multiple of every generator of this ideal. It therefore annihilates the product difference in \eqref{eq:scalar-fusion-product-ideal}, including its contributions from all collision jets. Substitution in \eqref{eq:scalar-fusion-residue} gives numerator residue \(\mathcal C_{\rm del}\tau_{\rm red}(f_{\rm red})\) and denominator residue \(\mathcal C_{\rm del}\tau_{\rm red}(1)\ne0\). Their common factor cancels. The assumed regularity lets the row variables approach any prescribed \(u_a(0)\), proving \eqref{eq:scalar-two-site-fusion}; the stated meromorphic continuation follows from the rational identities used throughout. \end{proof} \subsection{Removing the auxiliary pairs} \begin{proof}[Proof of Theorem~\ref{thm:row-moments}] First take the desired row arguments distinct and generic. Add one neighboring pair for each argument, with its first site $x$ equal to that argument. Apply the distinct-site identity \eqref{eq:site-moment} to the enlarged list, using the subset of all these first sites, and use Proposition~\ref{prop:common-twist} to replace $L$ by $P$ in the scalar trace. Remove the pairs successively at $x=y+3$. The denominator product at the first site satisfies $d_{\rm full}(x)\ne0$ on fusion. On the scalar side, Lemma~\ref{lem:scalar-two-site-fusion} gives the smaller normalized trace at each removal. On the transfer side, the vacuum injections and the $r=3$ transfer intertwining in Section~\ref{sec:local}, together with the vacuum fusion and left-cap contraction of Sections~\ref{sec:vacuum} and~\ref{sec:site-moments}, give the identical reduction. Thus the two sides pass through the same sequence: \[ \begin{aligned} \text{enlarged distinct-site moment} &\ \longrightarrow\ \text{successive two-site reductions}\\ &\ \longrightarrow\ \text{requested row moment on }\mathbf b. \end{aligned} \] This proves \eqref{eq:all-row-moments} for generic distinct arguments. Meromorphic continuation proves it for arbitrary row arguments at regular values, including repetitions. For the empty product, the identity is the normalization $\Lambda\Psi=1$ and $\tau(1)/\tau(1)=1$. \end{proof} ==== END Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//repeated-moments.tex ==== ==== BEGIN Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//residue.tex ==== \section{A residue formula with independent twists}\label{sec:residue} The next step is a scalar formula for the vacuum contraction. We keep the twists independent while taking residues so that the level intersections are transverse. The collision analysis below proves that the proposed residues exhaust the global residue theorem. Write \([u]=2\sin u\), shifts indicated by integers mean multiples of \(l=\pi/8\); all variables are taken modulo \(8l\) (binomial products below descend even when single brackets are antiperiodic). The variables lie in the complex torus $\mathbb C/(8l\mathbb Z)$, identified with $\mathbb C^*$ by $p=e^{2it}$. The regular configuration torus excludes every shifted root--root and root--site coincidence. For this calculation put \[ V(d)=-\frac1{[d-1][d+1]},\quad G(d)=-\frac{[d]^2[d-2][d+2]}{[d-1][d+1]},\quad S(d)=\frac{[d+1]}{[d-1]},\quad R(d)=\frac{[d-2][d+1]}{[d+2][d-1]} \] Here $V(d)$ and $R(d)$ denote scalar functions of a root difference, a convention local to the residue calculation. Take \(N\) sites (indexed below by \(b\)), split into lists \(X,Y\), with no two sites differing by a shift. In \(N\) root variables put \[ D_i=\prod_b S(t_i-b)\prod_{j\ne i}R(t_i-t_j),\qquad \Phi=\prod_{i,b}V(t_i-b)\prod_{i0$ put $H_i=D_i^F/(D_i^F+K_i)$. When several left multipliers occur, they may equivalently be inserted before the iterated residues: after insertion of its own multiplier an earlier variable has a simple pole, and every later multiplier is regular there for generic values of the other variables. The free pair factors remain $G$ in $\Phi_F$ and $R$ in $D_i^F$. The one-body factors at a left site $x$, up to nonzero constants independent of the free variables, are as follows. Here $d=t_i-x$ and $s$ is the number of slots already filled at $x$. \[ \begin{array}{c|cc} s&v_s(d)\text{ in }\Phi_F&d_s(d)\text{ in }D_i^F\\ 0&-1/[d+1]^2&S(d)\\ 1&[d+3]^2/[d+2]^2&[d+1][d+2]/([d][d+3])\\ 2&-[d+4]^2&[d+2]/[d+4] \end{array} \] Indeed freezing a left label at the additive slot $\xi$ multiplies a spectator's bare factor by \[ \frac{G(t_i-\xi)}{R(t_i-\xi)} =-\frac{[t_i-\xi]^2[t_i-\xi+2]^2}{[t_i-\xi+1]^2}, \] and its level function by $R(t_i-\xi)$. At a right site the two factors are $v_s(-d)d_s(-d)$ and $d_s(-d)^{-1}$, respectively. At the next left slot, $D_i^F$ has a simple zero and $H_i\sim D_i^F/K_i$; at the next right slot, $D_i^F$ has a simple pole and $H_i\sim1$. In each case the resulting form has a simple pole with exactly the scalar residue specified in the lemma. This also computes simultaneous next-slot residues when at most one label is frozen at each site and the remaining variables are generic. \subsection{Resolving collisions and coordinate boundaries} We need a boundary that records both relative rates of escape and collisions among variables at one scale. All bracket products are rational in $p_i=e^{2it_i}$; their finite zeros and poles lie on $p_i=q^k p_j$ or $p_i=q^k e^{2ib}$, where $q=e^{2il}$ and $k$ is taken modulo $8$. Thus the following boundary contains every divisor relevant to the form. \begin{claim}[A compact model adapted to the level equations]\label{cl:compactification} At each recursive stage, the regular configuration torus admits a smooth compactification, obtained by proper modifications of the projective ratio model, with a simple-normal-crossing boundary containing all coordinate and shifted-collision divisors. Every $D_i^F$ extends to a holomorphic map to $\mathbb P^1$. For generic independent nonzero twists, any collection of level divisors $D_i^F=-K_i$ meets every boundary stratum transversely whenever it meets it at all. \end{claim} Here and in the order calculation we suppress the superscript $F$ on the level functions. The construction has three parts. \paragraph{Relative scales.} Adjoin the fixed label $p_*=1$, and take the closure of the graph of all ratios $[p_a:p_b]\in\mathbb P^1$ for $a,b\in I\cup\{*\}$. For each ordering $\sigma(0),\ldots,\sigma(h)$ of these labels there is an affine chart with coordinates \[ u_a=\frac{p_{\sigma(a)}}{p_{\sigma(a+1)}},\qquad 0\le a0$ and zero if $e_i\le0$, away from its separate level divisor. This proves the last line. In particular, the maximum is taken after adding all depths in $e_i$. \paragraph{A lower bound valid for every final ray.} Let $A$ be the union of the free labels in the outermost positive-weight clusters anchored at left sites. Use \[ \sum_i\max\{0,e_i\}\ \ge\ \sum_{i\in A}e_i. \] The outermost clusters are disjoint. Every positive-weight cluster is contained in one of them or disjoint from it, by the cluster incidence rule including the site labels. Therefore every pair charge cancels in the sum on the right. Its only remaining terms are the site charges of the left-anchored clusters inside those outermost clusters. Combining them with the bare site orders in \eqref{eq:residue-bare-order} gives \begin{equation}\label{eq:residue-depth-bound} \operatorname{ord}_{\log,\rho}\omega_F \ \ge\ \mathfrak s_F(w)+\sum_C\lambda_C\,\ell(C), \end{equation} where the depth bounds are \[ \begin{array}{c|l} \text{cluster }C&\ell(C)\\\hline \text{unanchored}&B-1\\ \text{left site, stage }0&B-m_{-1}-m_1\\ \text{left site, stage }1&B-m_{-2}-m_0+m_{-3}+m_{-1}\\ \text{left site, stage }2&B+m_{-2}+m_{-4} \end{array} \] and right-site bounds are their reflections $k\mapsto-k$. All these bounds are nonnegative, with a short explicit list of zero cases. To see this, put $M=\max_km_k$ and $\Delta=\frac12\sum_k(m_k-m_{k+1})^2$. If every $m_k$ is positive, the distance-two sum is at least $\sum_km_k$ and strictly covers each subtracted pair of counts. Otherwise \[ \Delta\ge\frac12\sum_k|m_k-m_{k+1}|\ge M. \] Equality requires one contiguous support interval, rising and falling monotonically in steps of size at most one. For stage $0$, put $a=m_{-1}$ and $b'=m_1$. If only one is positive, equality in $B\ge M\ge a+b'$ requires that occupied offset to be a unit singleton or part of a unit adjacent pair. If both are positive, use $B\ge M+ab'\ge a+b'$; equality requires $a=b'=1$ and unit support $\{-1,0,1\}$. Thus the stage-$0$ zero supports are exactly \[ \{-1\},\quad\{-2,-1\},\quad\{-1,0\},\quad \{1\},\quad\{0,1\},\quad\{1,2\},\quad\{-1,0,1\}. \] Every listed count is one. The stage-$1$ expression is the stage-$0$ expression translated to the offsets $-2,0$, plus $m_{-3}+m_{-1}$. Its zero supports are therefore exactly $\{-2\},\{0\},\{0,1\}$. Stage $2$ has no zero case for a nonempty cluster. Finally $B=1$ has only a unit singleton or unit adjacent pair; an unanchored cluster has at least two labels, so its zero case is a unit adjacent pair. This proves Claim~\ref{cl:orders}. \paragraph{Which zero-order rays can occur.} For the total order in \eqref{eq:residue-final-order} to be zero, \eqref{eq:residue-scale-order} and \eqref{eq:residue-depth-bound} require $w=0$ and $\ell(C)=0$ at every positive depth. Equality in the bound on the level terms also requires \begin{equation}\label{eq:residue-equality-signs} e_i\ge0\quad(i\in A),\qquad e_i\le0\quad(i\notin A). \end{equation} These conditions, rather than the depth table alone, exclude most of its zero cases. For an outermost left-anchored cluster, the following table gives a label whose charge from that depth is negative. The charge includes both \eqref{eq:residue-pair-charge} and the site charge in \eqref{eq:residue-site-charges}. \[ \begin{array}{c|c|c|c} \text{stage}&\text{outer support}&\text{tested offset}&\text{charge}\\ \hline 0&\{-1,0\}&0&-1\\ 0&\{1\}&1&-1\\ 0&\{0,1\}&1&-2\\ 0&\{-1,0,1\}&1&-1\\ 0&\{1,2\}&2&-1\\ 1&\{0\}&0&-1\\ 1&\{0,1\}&1&-1 \end{array} \] Every strict positive-weight subcluster must again be one of the zero supports just listed if anchored, or a unit adjacent pair if unanchored. In the stage-$0$ rows tested at offset $1$, the only possible anchored subclusters containing that label are $\{1\}$ and $\{0,1\}$, with charges $-1$ and $-2$; in the other rows there is no anchored subcluster containing the label. Any unanchored subcluster containing a tested label is an adjacent pair with that label at its right endpoint, and contributes $-1$. Hence every eligible subcluster charge is nonpositive and the total $e_i$ of that label is negative, contradicting \eqref{eq:residue-equality-signs}. This leaves the next-slot singleton $\{-1\}$ at stage $0$, the next-slot singleton $\{-2\}$ at stage $1$, and one extra stage-$0$ support $\{-2,-1\}$. In that extra pair the charges at $-2,-1$ are respectively $+1,0$. The only strict zero-order subcluster containing the label at $-2$ is the unanchored pair, where its charge is again $+1$; an anchored $\{-1\}$ subcluster omits it. Thus the total charge of that extra label is strictly positive. Right-anchored clusters give the reflected list and opposite charges. An outermost unanchored zero support is an adjacent pair with charges $+1,-1$ and has no strict subcluster; its positive charge contradicts the second sign in \eqref{eq:residue-equality-signs}. We have now classified every possible boundary logarithmic pole, including rays combining several depths. \begin{claim}[The contributing boundary intersections]\label{cl:surviving-slots} Group each level pole by its free label. The boundary logarithmic poles can be assigned to those groups so that every intersection selecting one pole from each group consists of disjoint ordinary next-slot collisions and regular level intersections for the remaining labels. No additional collision or coordinate boundary contributes. \end{claim} Give group $i$ its level divisor $D_i=-K_i$. Every boundary logarithmic pole has, in at least one outermost cluster, a label projecting to an allowed next slot. Assign that component to the group of any such label; for an extra stage-$0$ pair choose its next-slot label, not the extra one. Claim~\ref{cl:compactification} implies that an intersection selecting one component from every group is transverse and finite. If it selects $r$ boundary components and $h-r$ levels, an extra boundary component would leave a stratum of dimension less than $h-r$, while an extra level would impose one more independent equation on the stratum of dimension $h-r$. Transversality excludes both. Suppose a selected component had an extra left pair, and let $i$ be its label at $x-2$. Its total charge is positive, so the regular map $D_i$ is identically zero on that component. Its finite nonzero level cannot be selected. Group $i$ would therefore have to select another boundary component, on which $i$ projects to an allowed next slot. At their intersection the two projections to the ratio model would agree, but $x-2$ is not the stage-$0$ next slot $x-1$, and no other site's shift orbit meets this one. This is impossible. At a right extra pair the same argument uses $D_i=\infty$. Every selected boundary ray therefore contains only next-slot singletons. Such a cluster has no strict subcluster, and two such clusters in one cone must be disjoint, including their site labels. The cone they span is already smooth. Each level charge on it is supported on its own label's slot ray, with one fixed sign, so none of the level hyperplanes cuts that cone; the refinement kept it unchanged. Thus its rays are the original codimension-one slot divisors, with at most one per site. A point with only those boundary components cannot project to any further collision or toric boundary: the inverse image of such a divisor is supported on rays with that incidence. The model is unchanged on this ordinary locus. The remaining levels are precisely regular transverse solutions of the restricted system, which proves Claim~\ref{cl:surviving-slots}. \subsection{Completing the residue sum} We apply the compact-manifold global residue theorem in the grouped-divisor form of Griffiths~\cite[Section~III(a), Equation~(3.7)]{Griffiths1976Abel}. Here compactness and the transverse finite intersections of the polar groups were proved above. The sum of the residues in group order is zero. One can see the sign directly: take residues along components of the first $h-1$ groups, leaving meromorphic differentials on smooth compact curves, and sum their ordinary residues. Poles at a second component of one of those first groups occur twice with opposite signs, by antisymmetry of the two residue orders. They cancel, leaving exactly the sum over one component of every group. For $h=1$ this is the ordinary residue theorem itself. Let $T_F$ denote the all-level residue of $\omega_F$, summed over its regular solutions; set $T_F=\Phi_F$ when $h=0$. At an all-level point, the local residue is $\Phi_F/\det(\partial_{t_j}\log D_i^F)$. By Claim~\ref{cl:surviving-slots}, every other grouped intersection freezes a nonempty subset $E$ of labels at distinct allowed next sites, then takes the all-level residue of the remaining system. Conversely each regular restricted solution occurs there. The local slot calculation therefore gives the recursion \begin{equation}\label{eq:residue-recursion} T_F=-\sum_{\substack{E\ne\varnothing\\ E\ {\rm an\ allowed\ simultaneous\ next\ freeze}}}T_{F\cup E}. \end{equation} Permuting the selected groups and the corresponding variables changes both orientations by the same sign, so no further sign occurs in this formula. Vanishing local residues may be included. Expand \eqref{eq:residue-recursion} to a full assignment $\mathcal L$. Residues at slots of distinct sites commute: all interaction factors between their disjoint shift orbits are regular, and the poles are ordinary transverse coordinate poles. Along one site the first slot must precede the second. Thus the scalar residue for $\mathcal L$ is independent of its division into successive nonempty layers, while a layering with $j$ layers has sign $(-1)^j$. For a fixed assignment with $c$ two-slot chains and $s$ single slots, $N=2c+s$. Let $a_j$ count its surjective layerings into $j$ ordered layers, strictly increasing along each chain. For nonnegative integers $p$, choosing the used layers gives the polynomial identity \[ \binom{p}{2}^{c}p^s=\sum_j a_j\binom{p}{j}. \] For the empty assignment the convention is $a_0=1$. Evaluating at $p=-1$ yields $\sum_j(-1)^j a_j=(-1)^s=(-1)^N$. Consequently every full assignment has the coefficient in the lemma. The $h=0$ convention supplies the terminal scalar residues, and the proof is complete. \end{proof} ==== END Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//residue.tex ==== ==== BEGIN Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//seam.tex ==== \begin{figure}[t] \centering \begin{tikzpicture}[x=1cm,y=.72cm,>=stealth] \draw[thick] (0,0)--(0,6) (4,0)--(4,6); \draw[dashed] (2,0)--(2,6); \fill (2,1) circle(2pt) node[right] {$0$}; \fill (2,5) circle(2pt) node[right] {$H$}; \node at (.35,.4) {$i$}; \node at (.35,2) {$1$}; \node at (.35,5.65) {$-i$}; \draw[blue!70!black,thick] (2,1) ellipse(.67 and .48); \draw[blue!70!black,thick] (2,1) ellipse(1 and .7); \draw[blue!70!black,thick] (2,5) ellipse(.8 and .53); \draw[orange!80!black,thick] (0,2.8) to[out=20,in=160] (2,3) to[out=-20,in=170] (4,2.8); \draw[orange!80!black,thick] (0,3.6) to[out=-15,in=170] (2,3.4) to[out=15,in=190] (4,3.6); \draw[<->] (0,-.4)--(4,-.4) node[midway,below] {circumference $N$}; \node[align=left,anchor=west] at (4.5,5) {contractible polygons\\around one tip}; \node[align=left,anchor=west] at (4.5,3) {winding polygons\\between the tips}; \end{tikzpicture} \caption{A schematic fundamental rectangle for the cylinder: its vertical sides are identified. The dashed seam has tips at cuts $0$ and $H$; the labels $i$, $1$, and $-i$ give its twist in the three row regions. The surviving contractible polygons enclose one tip, and the surviving winding polygons separate the tips. Dividing by the winding growth retains the two tip contributions. Curves are schematic and do not represent individual lattice edges.} \label{fig:seam} \end{figure} ==== END Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//seam.tex ==== ==== BEGIN Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//site-moments.tex ==== \section{Identifying the cap moments at distinct sites}\label{sec:site-moments} We compare the scalar trace with a contraction of the two disk spin caps when each row argument is a distinct site parameter. The comparison first determines the unnormalized scalar factor. Its empty-row case then fixes the denominator used for arbitrary repeated-row moments. Here and below an integer shift in a trigonometric argument denotes that integer times \(l\); brackets mean twice the sine. Use \(d,e,H,\Psi,Z\) from Sections~\ref{sec:local} and~\ref{sec:vacuum} and \(W,S,R_{\rm pair},V_{\rm pair},G,D_i,L,d_{\rm full}\) as in Lemma~\ref{lem:residue} and Proposition~\ref{prop:common-twist}. Thus \(R_{\rm pair},V_{\rm pair}\) here are the \emph{scalar} functions called \(R,V\) there, and \(d_{\rm full}(w)=\prod_b[w+3-b][w+2-b]\). Write \(\tau_N\) for the trace \(\tau\) on \(N\) sites and roots. We work first with ordered sites \(\mathbf b\) in distinct shift orbits and generic \(K\ne0\). A diagonal trace means the holomorphic specialization of the rational independent-twist trace near a generic diagonal \(K_i=K\), as in Proposition~\ref{prop:common-twist}; it includes contributions from colliding roots. For this \(N\)-site spin space write \(\Psi=\Psi_N\). Define \(\Lambda\) by taking the transpose of the right vacuum for the reversed site list, reversing its tensor slots back to the original order, and negating all spins. Since \(R^t\) equals simultaneous spin-negation conjugation of \(R\), this is the left vacuum at twist \(-i\). Reversal and negation together preserve matching tensors, so \(\Lambda\) lies in \(D^t\), has empty coefficient one, and is generically unique with that normalization. Moreover \(\Lambda\Psi=1\): two matching tensors pair only on identical occupancy, and an innermost arc then contracts the other matching by the pairing rules, eventually closing a zero loop unless the occupancy was empty. For a subset \(X\) of size \(k\), with complement \(Y\) of size \(m=N-k\), the comparison before normalizing the scalar trace is \begin{equation}\label{eq:site-moment-unnormalized} \Lambda \prod_{x\in X}T^{K^{-1}}(x)\Psi =\frac{(2\sqrt2)^N}{N!\,Z_N^2} \prod_{\{b,b'\}\subset\mathbf b}H(b-b')H(b'-b)\, \tau_N\!\left(\prod_{x\in X}\frac{L(x)}{d_{\rm full}(x)}\right). \end{equation} The product over \(\{b,b'\}\) takes each unordered pair of distinct sites once. The empty subset and \(\Lambda\Psi=1\) give the exact normalization \begin{equation}\label{eq:trace-normalization} \tau_N(1)= \frac{N!\,Z_N^2} {(2\sqrt2)^N\displaystyle\prod_{\{b,b'\}\subset\mathbf b} H(b-b')H(b'-b)}. \end{equation} This value is independent of the twist, including the independent twists before specialization. For \(N=0\), there is one empty root fiber, with \(W=J=1\); with \(Z_0=Z_1=1\), \(0!=1\), and empty pair products equal to one, the formula gives \(\tau_0(1)=1\) and \(\tau_1(1)=1/(2\sqrt2)\). In particular \(\tau_N(1)\ne0\) generically. Substituting this value in \eqref{eq:site-moment-unnormalized} gives \begin{equation}\label{eq:site-moment} \Lambda \prod_{x\in X}T^{K^{-1}}(x)\Psi = \frac{\tau_N\!\left(\prod_{x\in X}L(x)/d_{\rm full}(x)\right)}{\tau_N(1)}. \end{equation} These identities extend as identities of meromorphic expressions in the site and twist parameters, with pointwise evaluation at regular values. We use two further consequences of the cap construction. By the exchange and rotation identities of Section~\ref{sec:local}, \(\Lambda\) absorbs adjacent spectral swaps at exchanged arguments (the \(\check R\) are symmetric and preserve empty coefficient on matching tensors), and twisted rotations with twist \(-i\). On regular fusion divisors \(\Lambda E_r=\Lambda_{\rm red}\): uniqueness applies on the reduced list, the pairing rules give \(E_r^tD\subset D\) by contracting adjacent slots via \(C^t\) and squeezing slots with a single spin, and the empty normalization is preserved. The simultaneous disjoint specializations used here and for deletion fusions follow as well, since the symmetric \(Z\) satisfies the fusion recurrences with nonzero denominator at generic independent centers. Permute the listed sites to the front by exchange. At the first site itself the transfer equals a left rotation weighted by \(K^{-s_1}\) followed by adjacent swaps taking that rapidity back to the front (put \(R_{a1}=P_{a1}\) in the trace). Order the transfers with this evaluation rightmost by commutativity. Push those swaps from the rightmost transfer past remaining transfers into \(\Lambda\); now evaluate the next transfer at the new first site and repeat. Absorbing the accumulated rotations at twist \(-i\) gives \(\Lambda\eta^{\sum_{x\in X}s_x}\Psi\), \(\eta=i/K\). \subsection{Resolving the finite spin words} For an ordered block $I=(i_1,\ldots,i_r)$, define its diagonal auxiliary entries by \[ A^{I}_{\pm}(w)=\langle\pm|_a R_{a i_1}(w-b_{i_1})\cdots R_{a i_r}(w-b_{i_r})|\pm\rangle_a. \] Thus $A_+(w)$ below means $A^X_+(w)$ and $A_-(w)$ means $A^Y_-(w)$. We resolve the cap contraction by inserting the commuting spectral projectors of these auxiliary diagonal entries in the partial ordered row monodromies. Entries with a fixed extremal auxiliary spin commute at different arguments (fully polarized Yang--Baxter component, both as row and column). Within each block they conserve total spin and are triangular (rightmost changed label increases), with labels \(\sigma=s\) on \(X\), \(\sigma=-s\) on \(Y\); diagonal signatures are products over its sites of \(d,h,e\) at argument \(w-b\) for \(\sigma=+,0,-\), respectively. These signatures separate simple eigenlines by numerator zeros in each site's shift orbit on using a word-independent common denominator. The term of a joint word vanishes or both block label sums equal some \(r\), with multiplier \(\eta^r\), by conservation. Empty blocks have the obvious scalar conventions. It remains to evaluate \(\Lambda\Pi\Psi\) for each joint projector \(\Pi\). The next reduction removes zero labels while preserving this contraction. Split each zero site into the ordered pair \(b+1,b-1\), replacing its label by \((+,-)\). The monodromy fusion identity intertwines the block entries through the spin injection \(E_2\), and the signature matches there because \(d(t-l)e(t+l)=h(t)\). The signature remains simple on the fusion locus: in the common pair denominator representation the first numerator choices have zeros (relative to center) \((4,3),(4,1),(2,1)\), the second \((2,1),(2,-1),(0,-1)\) modulo 8. Separation at a generic test argument per block therefore makes the projector regular there. Applying vacuum fusion on both sides proves that \(\Lambda\Pi\Psi\) is unchanged. We can consequently compute the all-nonzero word on an enlarged generic list, performing all exchanges there before restricting the split pairs. For an all-nonzero word, the extremal orders identify the selected eigenlines. Within each block reorder rapidities to the minimum word (+ then -), or maximum word (- then +), keeping internal same-label orders fixed. Adjacent exchanges transport the projectors via Yang--Baxter. For fixed sum these words are extremal even allowing zeros, so the left eigenline in minimum order is the pure minimum covector, and the right eigenline in maximum order the pure maximum ket (simultaneously for both blocks). If \(U\) is the product of adjacent swaps taking maximum list to minimum, this gives the contraction \[ \frac{\Psi(\min)_{\min}\ \Lambda(\max)_{\max}} {\langle\min|U|\max\rangle},\qquad \langle\min|U|\max\rangle=\prod_{\text{within block }(+,-)} e(b_+-b_-). \] Indeed pass each + from first to last in their relative order across the minuses; to arrive still + it could not have donated any label spin, since each encountered minus is initially minimal; induction leaves the minuses unchanged. All these operations precede restriction to fused values. \subsection{Matching spin words with frozen slots} The projector calculation has reduced the matrix element to explicit products. We now assign a frozen-slot term to each word and check its phase and denominator. Recall that the label is \(\sigma=s\) on \(X\) and \(\sigma=-s\) on \(Y\), where \(s\) is the actual eigenword spin. Encode an original site's spin by the slot length \(j=1+s\): its slots are \(x-a\) on \(X\) or \(y+a\) on \(Y\), for \(a=1,\ldots,j\). Since both block label sums are \(r\), the numbers of occupied slots are \(k+r\) on the left and \(m-r\) on the right, hence \(N\) in total. To calculate the word's contribution, let its enlarged plus/minus \(\sigma\)-label counts on \(X,Y\) be \((A,B),(C,D_0)\); in particular a plus label on \(Y\) is an actual minus spin. Conservation gives \(M=A+D_0=B+C\). The minimum ket has actual spins \(+^A-^M+^{D_0}\), and the reversed negated maximum word has actual spins \(+^C-^M+^B\). Rotate the first plus segment to the end in each and use the two-block component formula \eqref{eq:two-block-component}; the phases are \(p^{2A-M+2C-M}=p^{2r}\), converting \(\eta^r\) to \(K^{-r}\). Besides \(K^{-r}\), the result has \(C_0=2\cos 3l\) per enlarged site divided collectively by \(Z_{\rm enlarged}^2\). Within each block equal labels contribute \(H(t)H(-t)\), \(t\) the difference; opposite labels contribute \(1/e(b_+-b_-)\). Across blocks, equal \emph{actual} spins \(n\) give \(H(-n(x-y))^2\); unequal spins give 1. Applying the fusion recurrence to \(Z_{\rm enlarged}^2/Z_N^2\) reorganizes this contribution as \[ \frac{K^{-r} (C_0\alpha)^N}{Z_N^2}\, \prod_{\{b,b'\}\subset\mathbf b}H(b-b')H(b'-b)\prod_{x,y}d(x-y) \ (-1)^{km+m}\operatorname{Res}_{\mathcal L} W\prod_{i,x}S(t_i-x)^{-1}\prod_{i\ {\rm left}}D_i , \] where $\alpha=4\cos l$ and \(\mathcal L\) uses any one labeling of the occupied slots. Here are details of the elementary reorganization. Besides the base pair products displayed, each occupied slot \(p'\) with another original site \(b\) has a factor \(V_{\rm pair}(p'-b)\) within a block and \([p'-b\mp1]^{-2}\) across (upper sign left slot). Slots at distinct sites interact by \(G(u)\) within block and \(([u][u-2]/[u-1])^2\) across, \(u=\) left minus right there. One verification uses offset/spin lists \(B_0=((0,-)), B_1=((1,+),(-1,-)), B_2=((0,+))\) (reflect offsets on Y). Explicitly put \(F_{ij}(t)=[t+3]^2[t+4]^2[t+5]^2\) for both indices 1, \([t+4]^2\) for just one, 1 otherwise (squared fusion factors). For \(z=t+a-b\) in the first identity, \(z=t+a+b\) in the second, \[ \begin{gathered} \prod_{(a,n)\in B_i,(b,n')\in B_j} \begin{cases}H(z)H(-z)&n=n'\\H(-nz)/([nz-1][nz])&n\ne n'\end{cases}\\ =F_{ij} H(t)H(-t)\prod_{c=1}^i V_{\rm pair}(t-c)\prod_{h=1}^j V_{\rm pair}(-t-h)\prod_{c,h}G(t-c+h),\\[4pt] \prod_{(a,n)\in B_i,(b,n)\in B_j} H(-nz)^2 =F_{ij} H(t)^2\prod_c[t-c-1]^{-2}\prod_h[t-h-1]^{-2}\\ {}\times\prod_{c,h}\left(\frac{[t-c-h][t-c-h-2]}{[t-c-h-1]}\right)^2 . \end{gathered} \] Both follow by inserting the three displayed lists, canceling with \([-v]=-[v], [v+8]=-[v]\). Use reflection of all rapidities for within Y. Single-site constants are \(C_0, C_0^2/((2\cos 2l)^2 e(2l))=1,C_0\); equivalently \(C_0\alpha^j\) times the successive left-slot residues of \(\prod V_{\rm pair}\prod G\) at that site alone (residues \(1,1/(2[2l]),[l]^2/(4[2l]^2)\) in ordinary angle coordinates), or their reflected right-slot residues multiplied by \((-1)^j\). Remaining factors match the displayed integrand since across blocks \(V_{\rm pair}S^{\pm1}\) and \(G R_{\rm pair}\) in oriented arguments give the stated cross corrections each with a minus. The total sign including reflection is \((-1)^{(k+r)m+(m-r)k+(m-r)(k+r)+m-r}=(-1)^{km+m}\). Finally \(K^{-r}=K^k\prod_{i\ {\rm left}} K^{-1}\) and \(\prod_x L(x)/d_{\rm full}(x)=(-1)^{Nk}K^k\prod_{xy}d(x-y)\prod_{i,x}S(t_i-x)^{-1}\). Each unlabelled frozen assignment has exactly \(N!\) labelings. At diagonal twists the residue lemma therefore turns the sum of the displayed word contributions into \eqref{eq:site-moment-unnormalized}: its sign agrees because \(N+Nk\equiv km+m\pmod2\), and its prefactor is \[ \frac{(C_0\alpha)^N}{N!\,Z_N^2} \prod_{\{b,b'\}\subset\mathbf b}H(b-b')H(b'-b), \qquad C_0\alpha=(2\cos3l)(4\cos l)=2\sqrt2. \] For \(X=\varnothing\), the cap contraction is one, proving \eqref{eq:trace-normalization}. In the independent-twist residue lemma this same empty subset has only right slots and no left twist multipliers. Hence \(\tau_N(1)\) is already independent of every \(K_i\) and equals its generic diagonal value. Substitution of the normalization now proves \eqref{eq:site-moment}. ==== END Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//site-moments.tex ==== ==== BEGIN Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//slit.tex ==== \begin{figure}[t] \centering \begin{tikzpicture}[x=1cm,y=1cm,>=stealth,font=\small, primal/.style={draw=black!50,line width=.45pt}, dual/.style={draw=blue!65!black,line width=.9pt}, arrival/.style={draw=orange!80!black,line width=1.1pt,->}] \begin{scope} \node at (0,1.8) {(a) The edge slit}; \draw[primal] (0,.9)--(-1.55885,0)--(0,-.9) (0,.9)--(1.55885,0)--(0,-.9); \draw[black,line width=.7pt] (0,.9)--(0,-.9); \draw[dual] (-1.03923,.9)--(-.519615,0)--(-1.03923,-.9) (1.03923,.9)--(.519615,0)--(1.03923,-.9); \draw[dual] (-.519615,0)--(-.10,0) (.10,0)--(.519615,0); \fill[blue!65!black] (-.519615,0) circle(1.6pt); \fill[blue!65!black] (.519615,0) circle(1.6pt); \draw[blue!65!black,fill=white,line width=.8pt] (-.10,0) circle(2pt) (.10,0) circle(2pt); \fill (0,.9) circle(1.5pt) node[above] {$v_1$}; \fill (0,-.9) circle(1.5pt) node[below] {$v_2$}; \node[above right] at (.12,.13) {$e$}; \node[align=center] at (0,-1.65) {open circles: the two slit ports}; \end{scope} \begin{scope}[xshift=4.8cm] \node at (0,1.8) {(b) Exterior arrival}; \path[fill=blue!4] (-.55,0)--(-.125,.736122)--(.725,.736122) --(1.15,0)--(.725,-.736122)--(-.125,-.736122)--cycle; \draw[dual] (-.55,0)--(-.125,.736122)--(.725,.736122) --(1.15,0)--(.725,-.736122)--(-.125,-.736122)--cycle; \draw[dual,->] (-.55,0)--(-.27,.484974); \draw[dual,->] (-.55,0)--(-.27,-.484974); \draw[arrival] (-1.55,0)--(-.61,0); \fill[orange!80!black] (-1.55,0) circle(1.7pt) node[above] {$s$}; \fill (-.55,0) circle(1.7pt) node[below left] {$q$}; \node[align=center] at (0,-1.65) {$\Delta W=\pm4\pi/3$\\$2\cos(c\Delta W)=0$}; \end{scope} \begin{scope}[xshift=9.6cm] \node at (0,1.8) {(c) Interior arrival}; \path[fill=blue!4] (-.35,0)--(-.675,.562917)--(-.35,1.125834) --(.30,1.125834)--(.625,.562917)--(.30,0) --(.625,-.562917)--(.30,-1.125834)--(-.35,-1.125834) --(-.675,-.562917)--cycle; \draw[dual] (-.35,0)--(-.675,.562917)--(-.35,1.125834) --(.30,1.125834)--(.625,.562917)--(.30,0) --(.625,-.562917)--(.30,-1.125834)--(-.35,-1.125834) --(-.675,-.562917)--cycle; \draw[dual,->] (-.35,0)--(-.57,.381051); \draw[dual,->] (-.35,0)--(-.57,-.381051); \draw[arrival] (-.025,0)--(-.30,0); \fill[orange!80!black] (-.025,0) circle(1.7pt) node[above] {$s$}; \fill (-.35,0) circle(1.7pt) node[left] {$q$}; \node[align=center] at (0,-1.65) {$\Delta W=\pm8\pi/3$\\$2\cos(c\Delta W)=-2$}; \end{scope} \end{tikzpicture} \caption{The slit and the two first-repeat cancellations, with $c=3/8$. In (a), black and grey segments are primal tiling sides and blue segments are honeycomb bonds. The bond crossing the side $v_1v_2$ is cut at its midpoint $e$; the two slit ports are separated slightly for visibility. Panels (b) and (c) are schematic representatives. The orange leg is the first prefix arrival at $q$ from the slit source $s$. The path then traverses the blue cycle and returns to $q$, its first repeated dual vertex. The two arrows leaving $q$ indicate the two cycle orientations. Their additional rotations are $\pm4\pi/3$ when $s$ is outside the cycle and $\pm8\pi/3$ when $s$ is inside. In (c), the arrival leg occupies the reflex interior sector at $q$.} \label{fig:slit} \end{figure} ==== END Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//slit.tex ==== ==== BEGIN Cap-selected-amplitudes-and-triangle-chords-for-honeycomb-walks//vacuum.tex ==== \section{The normalized disk vacuum}\label{sec:vacuum} The cylinder estimate used one differentiated row. For repeated rows with different boundary twists, we need normalized right and left vacuum caps. This section constructs the right cap as a simultaneous fixed vector of the spin rows, determines its exact polynomial denominator, and fixes the normalization under fusion. The left cap will then follow by transpose and spin reversal in Section~\ref{sec:site-moments}. Polynomial exchange equations and special-rapidity reductions have earlier counterparts in the dense $O(1)$ sum rule of Di Francesco and Zinn-Justin~\cite{DiFrancescoZinnJustin2005} and in the dilute $O(1)$ ground-state work of Garbali and Nienhuis~\cite{GarbaliNienhuis2017Ground,GarbaliNienhuis2017Sum}, which gives a first normalization recurrence in a different convention. Our null-loop vector has different boundary data, and we prove its normalization directly below. \subsection{The disk matching space} We use the spin transfers \(T^\zeta(u;\mathbf b)\) of Section~\ref{sec:local}, with their labelled and possibly inhomogeneous site parameters, and take \(\zeta=i=p^{-2}\). The annular diagram space of Section~\ref{sec:cylinder} retains the homotopy of each arc around the cylinder. The space used here is instead a subspace of \(V^{\otimes N}\) defined by ordinary disk matchings, with no annular winding data. For each noncrossing partial matching of the ordered sites in a disk, form a tensor by placing \(|0\rangle\) at each unmatched site and inserting \(C\) at the left and right endpoints of each arc, in that order. Let \(D\) be the span of these tensors. Every vector in \(D\) has total spin zero. The matching tensors are linearly independent: first separate them by their occupancy sets, then test a matching against the spin word that assigns \(+\) to each opening endpoint and \(-\) to each closing endpoint. Any matching compatible with that word has at least its word height many unclosed arcs at each cut, with simultaneous equality only for the chosen matching. Ordering by these heights gives a triangular test with nonzero diagonal. We call the coordinate on \(|00\ldots0\rangle\) the empty coefficient. The adjacent operators \(\check R\) preserve \(D\). Their nonempty contractions use \(C^t\): this kills an arc \(C\), and otherwise joins its two partners with the prescribed phases, since cup-cap straightening uses \(p^m p^{-m}=1\). New arcs are inserted as \(C\). The last arc therefore cannot be capped off with nonzero weight; because the empty-to-empty entry is 1, \(\check R\) also preserves the empty coefficient. Twisted left rotation preserves matching tensors as well, by \(p^m\zeta^m=p^{-m}\). The full row has the same two properties. To express its twisted trace using disk matchings, append an auxiliary and a spectator in \(|00\rangle+C\) at the right. Pass the auxiliary from last to first through the sites using adjacent \(\check R\)'s, and then contract the first and last exterior legs, namely the auxiliary and spectator, by \(\langle00|-C^t\). This gives exactly \(T^\zeta\), because \(-p^{2m}=\zeta^m\) for \(m=\pm1\). If the occupied exterior legs are joined to each other, their contraction is zero. If they are joined inward to partners, the contraction forms a new arc whose left spin is \(-m\), with phase \(-p^{3m}=p^{-m}\). Thus the output remains in \(D\). Again a last arc cannot disappear with nonzero weight, so only a completely empty input contributes to the empty output, with coefficient 1. Let \(D'\subset D\) be the invariant subspace of vectors with empty coefficient zero. \subsection{Construction of the generic fixed vector} We claim that for generic site parameters there is a unique \(\Psi_N(\mathbf b)\in D\) with empty coefficient 1 such that \(T^\zeta(u)\Psi_N=\Psi_N\) for every regular \(u\). Its components are rational in the half variables \(y_j=e^{ib_j}\). Set \(\Psi_0=1\) and \(z_j=y_j^2\). The proof finds finitely many row arguments for which a linear combination of \(T^\zeta-I\) is invertible on \(D'\); the empty normalization then gives the fixed vector by linear equations. For that purpose, let the sites have very separated scales in \(\mathbb C^*\): successive ratios of their moduli, in increasing order, tend to infinity. Evaluate the row at intermediate scales \(z=e^{2iu}\), one in each gap, with its ratios to the scales on either side tending to their respective extremes. For \(z/z_j\to\infty\), the local weights satisfy \(a,b\to0,\ h\to1,\ d\to q^3,\ e\to q^{-3}\). The limit preserves occupancy on both legs, and the auxiliary spin can only increase or stay fixed: on \(|1,-1\rangle\), the opposite-spin off-diagonal weight is \(q^3+p^{-2}q^{-3}=0\). Its diagonal entries are \(q^{3s_a s_j}\). For \(z/z_j\to0\), the limits of \(d,e\) are exchanged; the auxiliary spin can only decrease or stay fixed, and the diagonal entries are \(q^{-3s_a s_j}\). We will also need the same tests when the sites form nonempty blocks of comparable scales, with successive blocks very separated. Place the test scales between those blocks. The limiting rows preserve every site's occupancy and are simultaneously triangular on the full spin-zero sector, in any linear ordering compatible with the componentwise order of the block prefix sums. The blocks are ordered by increasing scale, which need not agree with their order around the disk. For a cut before the test gap, every spin in its prefix can only decrease. For a cut after the gap, every spin in the complementary suffix can only increase, so total-spin conservation again makes the prefix sum nonincreasing. If every block sum is preserved, all changes within a block have the same sign and hence every site spin is unchanged. At the cut after block \(k\), let \(S_k\) be the spin in the smaller blocks. A contribution on the diagonal has constant auxiliary spin, so the diagonal entry of the limiting row is \(1+i q^{6S_k}-i q^{-6S_k}\). All the tested diagonal entries can equal 1 only when every block spin is even. Consequently, on every occupancy sector containing an occupied singleton block, a generic linear combination of the limiting \(T^\zeta-I\) is invertible, with one choice working for all such sectors. This is the interblock exclusion we shall use below. For a fully separated list, every occupied sector has an occupied singleton block. Only the empty tensor remains, so the same combination is invertible on \(D'\), also before the limit when the scales are sufficiently separated. For \(N=1\), \(D'=0\) directly. Choose one regular configuration where invertibility holds and keep its test arguments and coefficients fixed. Rationality then makes invertibility generic in the site parameters. The combination maps \(D\) into \(D'\) and is invertible on \(D'\); its kernel in \(D\) is therefore one-dimensional and contains exactly one vector of empty coefficient 1. Since all rows commute with the combination and preserve the empty coefficient, this vector is fixed by every row. Solving the finite linear equations proves rationality in \(\mathbf y\) as well as the claimed uniqueness. The row identities of Section~\ref{sec:local} now transfer to the normalized vector by uniqueness. It is invariant under a common translation of the site parameters, covariant under a period shift at one site by the sign of occupation there, and covariant under twisted cyclic rotation. Its adjacent exchange equation is \[ \Psi_N(\ldots,A,B,\ldots)=\check R_{j,j+1}(A-B)\Psi_N(\ldots,B,A,\ldots). \] \subsection{Counting the poles} Existence and uniqueness determine the vacuum as a rational vector. To use its entries in the residue formula, we need its least common denominator, including the multiplicities of its factors. We first count poles while varying one site. The endpoint limits will then determine the leading coefficients, and symmetry and total degree will leave the fusion factors to be fixed. For each component, divide by \(\prod_{j:\,\text{occupied}}y_j\). Sign covariance makes the result rational in \(\mathbf z\). Let \(\mathscr D_N(\mathbf z)\) be the least polynomial common denominator of these rational components on \((\mathbb C^*)^N\), chosen with no coordinate factor. It is unique up to a nonzero constant. We will prove that \[ \deg_{z_j}\mathscr D_N=N-1\quad(1\le j\le N),\qquad \deg\mathscr D_N=\frac{N(N-1)}2, \] that \(\mathscr D_N\) is symmetric and homogeneous, and that none of the divisors \(z_B/z_A=q^{\pm2},q^{\pm3},1\) is a factor. These divisor assertions concern generic points of each divisor; their intersections require separate justification when used. Fix all sites except \(B\) and vary \(z_B\). Choose the fixed parameters generically, outside the loci where specialization creates cancellation, a degree drop, or a new coordinate factor, and then let their scales become very separated within this set. Such choices form a Zariski open set: a one-variable gcd identity keeps the reduced common fractions reduced after specialization. At a pole of the normalized vector, no combination of regular row transfers minus identity can be invertible on \(D'\), since the normalized linear equations would otherwise solve holomorphically there. Consider any sequence of such poles as the fixed scales become more separated. After passing to a subsequence, the scales separate into singleton blocks except possibly for one pair \(A,B\) with comparable modulus. If no such pair remained, the fully separated tests would give invertibility on \(D'\) and rule out the poles. We may therefore assume that \(z_B/z_A\) has a nonzero finite limit, and pass to further subsequences to fix the square-root ratios. The interblock exclusion rules out all occupancy sectors with an occupied site outside this pair. It remains to examine the pair. Evaluate additional rows at arguments for which \(e^{iu}/y_A\) is fixed and generic. In the limit these rows preserve the outside occupancies; on the subspace with empty spectators they reduce to the two-site transfer. The outside-occupancy decomposition also splits \(D\) and \(D'\). Hence degeneracy on \(D'\) requires a simultaneous eigenvalue 1 on the two-site \(D'\), which is the line spanned by \(C\). This follows either by commutation or by taking a generic combination of all the tests. \subsection{Localization and existence of the simple poles} Call the earlier site of the pair in disk order \(F\), and the other \(G\). In the two-site problem, evaluation at the first site \(u=b_F\) gives \(R_{FG}(b_F-b_G)\zeta^{s_F}\). On the disk matching space this acts by \[ C\mapsto d(b_F-b_G)C,\qquad |00\rangle\mapsto |00\rangle+b(b_F-b_G) C. \] At a regular first-site evaluation, subtracting the numerator and denominator of \(d(t)\) shows that \(d(t)=1\) requires \(\sin(2t)=0\). The case \(t=0\bmod\pi\) is excluded by evaluating instead at \(u=b_F+l\). Up to the sign gauges, the eigenvalue on \(C\) is then \(2x^4\ne1\): at \(l\) we have \(e=0,\ h=d=x^2\), and for outgoing pair spins \(m,-m\), the vacant traced auxiliary contributes \(p^m h^2\), while the occupied ones contribute \(\zeta^{-m}p^{-m}d^2=p^m d^2\). This test also excludes a limiting pair for which the first-site evaluation itself has a denominator pole. Multiply the two-site operator by its two scalar sine denominators. At \(u=b_F\), its restriction to \(D'\) minus the cleared identity is nonzero at any such pole, because the numerator of \(d\) is nonzero there. It is consequently nonzero at some nearby regular row argument for that limiting pair. The only remaining possibility is \(t=\pi/2\bmod\pi\), or \(z_B/z_A\to-1\). We now count the poles near each possible limit \(-1\). Take the interblock tests together with the first-site test \(u=b_F\), allowing that argument to move with the varying site. On a small disk in \(z_B/z_A\) about \(-1\), the transfers are analytic on each square-root sheet for sufficiently large separations and converge uniformly to the limiting operators. On \(D'\), those limits preserve the outside occupancies. Choose a generic linear combination of their differences from identity. On sectors with occupied spectators, the interblock exclusion makes it invertible at the center; the pair test can be given a sufficiently small nonzero coefficient to preserve this. On the remaining one-dimensional sector, the interblock differences vanish and the pair test has the simple zero of \(d-1\). The determinant of the combination on \(D'\) therefore has one simple zero at the center in the limit. The argument principle and the normalized linear equations show that, for sufficiently separated scales, the vector has at most one pole in the disk, and that pole is simple. There is in fact a pole in each disk. On the whole spectator-empty pair sector in \(D\), the interblock transfers are the identity: both pair spins lie on the same side of a test, so neither can change, and their sum is zero. On a sufficiently small circle about \(-1\), the normalized vector therefore converges to the unique normalized two-site solution. Its coefficient of \(C\) is \(b/(1-d)\), which has nonzero residue. Contour integration forces a pole of the vector inside the circle. The two square-root signs give equivalent gauges. Passing to subsequences makes the argument apply along any sufficiently separated sequence, since the limiting tests do not depend on the extreme ratios. There is thus exactly one simple pole near each antipodal value \(z_B=-z_A\), for the \(N-1\) fixed sites \(A\), and none away from those disks. Generic specialization gives \(\deg_{z_j}\mathscr D_N=N-1\) for every \(j\). The same tests show that \(\mathscr D_N\) has none of the divisors \(z_B/z_A=q^{\pm2},q^{\pm3},1\): make this pair comparable and the other sites separated, and use the pair tests away from \(-1\). They also show that a combination is invertible on \(D'\) at generic points of each such divisor. For \(N=1\), the degree and these assertions hold directly. The absence of the equal-site divisor is only a generic statement on that divisor; it does not by itself establish regularity when all sites coincide. Section~\ref{sec:physical-amplitude} proves that regularity for the sufficiently large sizes used below, from the projected annular row and its spectral gap. \subsection{Endpoint limits, symmetry, and degree} The pole count controls finite nonzero site values. We next let one site tend to zero or infinity while the others remain generic and fixed. In either limit, \(\Psi_N\) is bounded and converges to the embedded \(\Psi_{N-1}\) with that site vacant. To prove this, evaluate one row at an intermediate scale separating the single departing site. Its triangular limit preserves occupancy and excludes an occupied singleton. Other fixed generic row arguments also preserve that site's occupancy in the limit and reduce on its vacant block to the shorter transfers. Together they give an invertible combination on \(D'\). On the vacant-site spin-zero block, the intermediate-scale test is the identity: monotonicity prevents any other site from changing, and the relevant prefix spins vanish. Thus the normalized solution of the limiting equations is precisely the embedded shorter vacuum. Invertibility gives boundedness and convergence to that solution. Throughout these arguments, the rows are evaluated at regular arguments before and in the indicated limits, and degeneracy is tested by finite linear equations. Exchange and inversion now imply that \(\mathscr D_N\) is symmetric. Away from the exceptional ratio divisors, the exchange matrix is regular and invertible, so it preserves all polar divisors and their multiplicities. The exceptional divisors themselves are absent by the pole tests. Primitivity therefore makes a permutation act on \(\mathscr D_N\) by a scalar; an alternating sign is impossible because the equal-site divisor is absent. A common translation of the site parameters gives a common dilation of \(\mathbf z\). Translation invariance of the vector and primitivity similarly make \(\mathscr D_N\) homogeneous. Its total degree is \(M_N=N(N-1)/2\). Choose rationally independent real slopes \(w_j\) in any increasing order and let \(z_j=\exp(Lw_j)\), \(L\to\infty\). Nonzero exceptional polynomials arising in fixed-site specializations stay nonzero for large \(L\). The root count and locations then give \(z_j\partial_{z_j}\mathscr D_N/\mathscr D_N\) tending to the number of smaller sites. On the other hand, rational independence makes one monomial of \(\mathscr D_N\) uniquely dominant along this path, and these limits are its exponents. Their sum is \(0+1+\cdots+(N-1)=M_N\), which is the total degree by homogeneity. Finally, regard \(\mathscr D_N\) as a polynomial in one site and pass to the endpoint limit in the divided components that are vacant at that site. Its leading coefficient times each limiting reduced component is a Laurent polynomial. Thus the leading coefficient clears every reduced denominator, and minimality of \(\mathscr D_{N-1}\), with no coordinate factors, makes it divisible as a polynomial by \(\mathscr D_{N-1}\). The total degree makes the quotient a constant. Normalize \(\mathscr D_0=\mathscr D_1=1\), and recursively choose this leading coefficient to be \(\mathscr D_{N-1}\). Symmetry makes the choice independent of the site. \subsection{Fusion normalization} With this normalization, define \[ Z_N=\mathscr D_{N}(\mathbf z)\prod_j y_j^{-(N-1)},\qquad \mathcal U_N=Z_N\Psi_N . \] Both are Laurent polynomials in \(\mathbf y\). For adjacent slots \(d_0+rl,d_0\), \(r=2,3\), let \({\rm red}\) denote the shorter list in which the pair is replaced by \(d_0+l\) for \(r=2\), or deleted for \(r=3\). We claim the recurrences \begin{equation}\label{eq:vacuum-fusion} \Psi_N|_{\rm pair}=E_r\Psi_{\rm red},\quad Z_N|_{\rm pair}=2\cos(rl)\prod_{j\ {\rm other}}\prod_{h'=1}^{r-1} 2\cos(b_j-d_0-h'l)\ Z_{\rm red}. \end{equation} Here \(E_r\) is inserted in the indicated positions. The first formula follows at generic points of the pair divisor: \(E_r\) maps matching tensors into \(D\), preserves the empty normalization, and intertwines the transfers, so generic uniqueness applies there. To prove the second formula, first locate the extra zeros on a pair divisor. All specializations of Laurent polynomials in this argument take place on the torus, with a chosen square-root sheet for the half variables. An ordered adjacent triple \(c+6l,c+3l,c\) forces \(\mathcal U_N=0\): polynomial continuation from the two pair divisors puts it in both \(E_3\) images, and these images have zero intersection because the two-site vector has matrix rank three. This wheel zero persists in cyclic orders with generic spectators inserted, by moving those spectators away through regular exchanges and then rotating. It also persists modulo \(\pi\) by the sign gauges. On the specified adjacent pair divisor, the wheel zero therefore occurs at every spectator value \(b_j=d_0+h'l+4l\bmod\pi\). For \(r=3\), these values complete either end of the ordered wheel. For \(r=2\), the cyclic triple is \(d_0+2l,d_0,d_0+5l\), which modulo the period rotates to \(d_0+8l,d_0+5l,d_0+2l\). Since \(E_r\) is injective, the first fusion formula and polynomiality of \(\mathcal U_N\) show that \(\mathscr D_N|_{\rm pair}\Psi_{\rm red}\) is regular on the torus and vanishes on each of these extra divisors. Thus \(\mathscr D_N|_{\rm pair}\) is divisible as a polynomial by \[ \mathscr D_{\rm red}\prod_{j\ {\rm other}}\prod_{h'=1}^{r-1}(z_j+q^{h'}z_{d_0}),\qquad z_{d_0}=e^{2id_0}. \] This divisibility includes multiplicities even if one of the extra factors coincides with a divisor of \(\mathscr D_{\rm red}\). Indeed, at any divisor the minimum component order of \(\Psi_{\rm red}\) is minus the order of \(\mathscr D_{\rm red}\), by minimality of the denominator and the empty normalization. Adjoining a dummy variable or taking square roots on the torus does not change that assertion. Regularity supplies the order from \(\mathscr D_{\rm red}\), and the wheel zero supplies the additional order. The quotient after division has total degree 1 and no spectator degree, so it is proportional to \(z_{d_0}\). After centering, this gives the displayed recurrence up to a constant. Compare the highest terms in each spectator in succession. The leading coefficient of \(y_j^{N-1}\) in \(Z_N\) is \(Z_{N-1}\prod_{k\ne j}y_k^{-1}\), and the removed pair product \(e^{i(2d_0+rl)}\) equals the product of the half variables contributed by the replacement site, when present, and the \(h'\)-factors. The comparison reduces the constant to the two-site case, where symmetry and normalization give \(Z_2=2\cos(b_1-b_2)\). This proves the second recurrence in \eqref{eq:vacuum-fusion}. \subsection{An explicit component} The fusion recurrences also determine the component with two consecutive blocks of opposite spins. This component will fix the normalization of the cap moments in Section~\ref{sec:site-moments}. For \(N=2m\) and \(\epsilon=\pm1\), write \(\epsilon^m(-\epsilon)^m\) for the spin word consisting of \(m\) copies of \(\epsilon\) followed by \(m\) copies of \(-\epsilon\). Its component is \begin{equation}\label{eq:two-block-component} (\mathcal U_N)_{\epsilon^m(-\epsilon)^m} =(p^\epsilon\,2\cos 3l)^m\! \prod_{\substack{i