1. Introduction
Let (X,\mathcal F,\mu) be a probability space and let T:X\to X be an invertible measure-preserving transformation with measurable inverse. Given bounded measurable functions f_1,f_2,f_3 and distinct nonzero integers \lambda_1,\lambda_2,\lambda_3, consider
The issue is convergence at almost every point, rather than only in norm. We assume that T is mixing, in the sense that
Theorem 1.1. For every system satisfying (1.2), every triple of pairwise distinct nonzero integers \lambda_1,\lambda_2,\lambda_3, and every bounded measurable f_1,f_2,f_3:X\to\C, the averages in (1.1) satisfy
as N tends to infinity through all positive integers. The exceptional null set may depend on the slopes and the functions.
The probability space need not be standard or countably generated. Both positive and negative slopes are allowed. The proof depends on the analytic results of the companion manuscripts [H, V]; their role and the extension required here are described below.
1.1. Context and predecessors
Birkhoff’s pointwise ergodic theorem [Bir] treats one function. Products along several iterates arise naturally in Furstenberg’s ergodic proof of the multiple recurrence theorem and its connection with arithmetic progressions [Fur]. Multiple recurrence, norm convergence, and pointwise convergence are different questions. In particular, convergence in L^2 does not by itself establish convergence at almost every point along the full sequence.
For weakly mixing transformations, the product-of-integrals limit in L^2 is classical. Bergelson’s weakly mixing PET theorem [Ber, Theorem 1.2 and p. 340] includes distinct nonzero integer slopes of either sign. For general systems, Host and Kra [HK] and Ziegler [Zie] developed the norm-convergence theory using characteristic factors with nilpotent structure. These results explain the limiting structure without resolving the corresponding pointwise question.
Bourgain’s double recurrence theorem [Bou] establishes almost-sure convergence for two bounded functions evaluated at two powers of one transformation. Pointwise theorems for longer products also have substantial precedents under additional assumptions on the system. For example, Derrien and Lesigne [DL, Theorem 4.2] prove polynomial pointwise theorems for K-systems and exact endomorphisms; their invertible case includes the linear configurations considered here. Other structural hypotheses yielding pointwise convergence for all finite progression lengths include singular spectrum on the Pinsker factor of a weakly mixing system, by Assani [Ass], and the pairwise-independent-joining property for weakly mixing systems, by Gutman, Huang, Shao, and Ye [GHSY]. The significance of Theorem 1.1 is that its dynamical assumption is mixing alone, while the three slopes are arbitrary distinct nonzero integers.
The harmonic-analytic method separates convergence from identification of the limit. A uniform oscillation estimate gives pointwise convergence; mixing identifies that limit by an L^2 argument. Oscillation methods in ergodic theory are developed in [Bou, JKRW], and the passage from translation estimates to orbit estimates follows the transference principle of Calderón [Cal]. We give that passage directly, so it imposes no additional structure on the probability space.
1.2. Analytic input and the coefficient extension
The manuscript [H] proves an L^3 bound for the trilinear Hilbert transform through local four-linear forms. Its proof includes chart detection, localization by quadratic data, comparison on short depth blocks, and a separated-frequency counting estimate. The quantitative inverse theorem of Leng, Sah, and Sawhney [LSS] enters its chart detection argument. Arithmetic regularity and inverse-theorem methods also occur in Tao’s sublogarithmic cancellation theorem for truncated multilinear Hilbert transforms [Tao]; that cancellation theorem is distinct from the uniform local estimate used here.
The companion [V] adapts the local estimates to inputs that vary with depth and deduces pointwise convergence for the slopes 1,2,3. The local forms in both companions use the four coordinates x,x+t,x+2t,x+3t. A power of T and a permutation of the functions do not turn an arbitrary triple into this configuration. Thus their stated theorems cannot simply be invoked for the slopes in Theorem 1.1.
The new work is the extension of their local proofs to x+b_jt, 0\le j\le3, for fixed distinct integers with b_0=0. Three features are particularly useful beyond this application. First, an integer left inverse of the whole primitive configuration replaces a unimodular pair of coordinates; this removes extra lattice or torus configurations in discretization. Second, a Vandermonde relation provides the quadratic cancellation and the Fourier randomization for every set of four distinct slopes. Third, a height base divisible by all pair differences preserves the divisibility required in the separated-frequency argument. The last step must preserve a bound polynomial in the dimension: an exponential loss there would not fit the proof’s parameter choices.
These are extensions of the proofs in [H, V], not new proofs of their coefficient-independent projection and row estimates. We state the local conclusion completely, prove the coefficient-sensitive identities and estimates, and specify the results imported from the companions at their points of use. In particular, the technical part is a proof by reference to those manuscripts, whereas the ergodic deduction is presented in full.
1.3. Proof structure
Section 2 defines the admissible kernels and states the local estimate for inputs with controlled variation on blocks of depths. Section 3 proves Theorem 1.1 from that estimate: smooth scale differences are realized by cancellative kernels, then transferred to integer sequences and finite orbits. The resulting oscillation bound gives convergence for smooth profiles; moment-corrected approximations to interval indicators recover all integer lengths.
The remaining sections prove the local estimate. Section 4 establishes the coefficient modifications in detection and localization. Section 5 proves the short-block, height, and counting modifications. Section 6 closes the parameter argument, first for fixed inputs and then for depth-dependent inputs. All analytic constants may depend on the fixed slopes; no bound is required uniformly over triples of integers.
2. The local estimate and the coefficient algebra
The analytic estimate concerns four inputs on the real line. Its constant may depend on the slopes, but must be independent of the number of scales and of the number of changes of the inputs. This uniformity will give the oscillation estimate in Section 3.
2.1. Normalization and local forms
Set
Then S is invertible, probability preserving, and mixing, and T^{\lambda_j n}=S^{b_jn}. Henceforth the four distinct integers b_0,b_1,b_2,b_3 are fixed, and \gcd(\abs{b_1},\abs{b_2},\abs{b_3})=1. We make no assertion of uniformity as these integers vary.
Dyadic intervals are half-open. The depth of a descendant I of a dyadic interval J is s(I)=\log_2(\abs{J}/\abs{I}). A finite rooted tree contains every ancestor between any of its intervals and its root. The dyadic lattice may be translated.
Definition 2.1 (Local kernels). Fix \delta_0>0 and C_0\ge1. For an interval I of length r, a kernel w_I\in C_c^\infty(\R^2;\C) is admissible if
- 1.
\dist(x+b_jt,\partial I)\ge\delta_0r and x+b_jt\in I for every (x,t)\in\supp w_I and 0\le j\le3;
- 2.
for every v\in\R and 0\le j\le3,
\int_\R w_I(v-b_jt,t)\,dt=0; \tag{2.1} - 3.
for nonnegative integers a,d and n=a+d,
\norm{\partial_x^a\partial_t^d w_I}_\infty \le r^{-n}\bigl(C_0(n+2)\bigr)^{C_0(n+2)}. \tag{2.2}
For scalar measurable inputs define
whenever the integral is absolutely convergent.
For a nonempty integer block B=\{u,\ldots,v\} and functions z_s on a common domain, put
This is a pointwise bound on both size and variation within the block; no variation across block boundaries is included.
Theorem 2.2 (Local estimate with variation in depth). Fix the primitive distinct integers b_0=0,b_1,b_2,b_3 and the parameters \delta_0,C_0. Let \mathcal T be a finite dyadic tree with root J and largest depth d, and let \mathcal I\subseteq\mathcal T. Assign an admissible kernel to every I\in\mathcal I. For each 0\le j\le3, let z_{j,s}:J\to\C, 0\le s\le d, be measurable, and partition \{0,\ldots,d\} into W_j\ge1 deterministic nonempty consecutive blocks. Suppose that
for each block in the jth partition. Extend every input by zero outside J. Then
The constant is independent of the lattice, tree, depth range, subcollection, kernels, and the four partitions and their sizes.
The schedules are prescribed at every depth, even depths with no summation intervals. This is exactly the hypothesis of [V, Theorem 3.2], except for the coefficients in the local forms. Sections 4–6 prove the theorem by extending that proof and its analytic inputs from [H]. Before turning to those technical arguments, Section 3 proves the implication from Theorem 2.2 to the ergodic theorem.
2.2. Integer reconstruction and quadratic cancellation
Two features of the coefficients recur in both the analytic proof and transference. Pairs of coordinates determine (x,t) over \R, but need not determine it over \Z. All four coordinates together do determine the integral step. Quadratic cancellation requires a separate identity, furnished by Vandermonde interpolation.
Lemma 2.3 (Coefficient identities). There exist integers e_0,e_1,e_2,e_3 and nonzero integers c_0,c_1,c_2,c_3 such that
and
The nullspace over \R of the two rows (1)_j,(b_j)_j is spanned by (c_j)_j and (b_jc_j)_j. Consequently, for y_j=x+b_jt,
For any distinct i,j, the pair map (x,t)\mapsto(x+b_it,x+b_jt) is injective on \Z^2 and pushes normalized Haar measure on (\R/L\Z)^2 to normalized Haar measure, for every L>0.
Proof. Bezout’s identity gives integers e_1,e_2,e_3 with \sum_{j=1}^3e_jb_j=1; set e_0=-\sum_{j=1}^3e_j. Let P(z)=\prod_{j=0}^3(z-b_j). Lagrange interpolation of z^a, 0\le a\le2, followed by comparison of the coefficient of z^3, gives \sum_jb_j^a/P'(b_j)=0. Multiply by a common positive integer multiple of the denominators to obtain (2.6) with every c_j\ne0. The two displayed null vectors are independent because the b_j are distinct, and the nullspace has dimension two. The identities for y_j follow by expansion. Finally the pair matrix has determinant b_j-b_i\ne0, so it is injective on the lattice and surjective on the real torus. A surjective continuous homomorphism of compact groups pushes normalized Haar measure to normalized Haar measure.
Remark 2.4. Primitivity does not mean that some pair is unimodular. For example, the slopes (b_0,b_1,b_2,b_3)=(0,6,10,15) are primitive, but no pair difference is 1 in absolute value. This is why the integer inverse in (2.5), rather than division by a pair determinant, is used in the lattice and torus arguments below.
2.3. The fixed-input constant
The proof of Theorem 2.2 first establishes a fixed-input estimate. For 2<q<3, let M_q f=(M(\abs f^q))^{1/q}, where M is the uncentered Hardy–Littlewood maximal operator. Let C_N(q) be the least constant for which
for all finite collections \mathcal A of intervals in one translated dyadic lattice, whose lengths belong to at most N consecutive scales, admissible kernels, and compactly supported scalar L^q inputs. This is the definition in [H, Section 2], with our coefficients. Pair-coordinate estimates imply C_N(q)\lesssim_bN, so this constant is finite before any absorption argument. The analytic proof will show \sup_NC_N(q)<\infty for a suitable fixed q>2, and then use that bound in the scheduled-input argument. This order avoids using Theorem 2.2 to prove its own analytic input.
3. From the local estimate to pointwise convergence
We now deduce Theorem 1.1 from Theorem 2.2. This section is independent of the internal notation of the analytic proof. It follows the transference scheme of [V, Sections 2–3], including its use of signed smooth profiles. We give all the details, both to handle arbitrary integer slopes and to separate the two roles of the hypotheses: the local estimate gives almost everywhere convergence of smooth averages on every invertible probability preserving system; mixing identifies their limit.
Throughout, b_0=0,b_1,b_2,b_3 are the fixed primitive slopes, S is an invertible measure-preserving transformation, and e_j are the integers in (2.5). All constants in this section may depend on the slopes.
3.1. Smooth profiles and cancellation kernels
Definition 3.1 (Averaging profile). An averaging profile is a real \varphi\in C_c^\infty(\R) with
and, for some G\geq1,
Write \varphi_a(t)=a^{-1}\varphi(t/a) for a>0. Profiles are necessarily signed: a nonnegative smooth function of mass one cannot have zero second moment.
The derivative class in (3.2) is closed under fixed translations, dilations, differentiation, and finite linear combinations. It contains a nonnegative bump \rho of mass one supported in [-1,1]. For completeness, normalize
For |t|<1, put d=1-t^2. On a complex disk of radius cd about t, with sufficiently small absolute c>0, the real part of (1-z^2)^{-1} is at least c'/d. Cauchy’s estimate gives
For each n this tends to zero at the endpoints, and its maximum over 0<d\leq1 is at most C^{n+1}n!(n+1)^n. Thus the zero extension is smooth and satisfies (3.2) after enlarging G. Lemma 3.9 will construct averaging profiles arbitrarily close to interval indicators.
Fix a profile \varphi, and put
Lemma 3.2 (Realizing a scale difference). There exist a dyadic integer D\geq1, a constant C_0\geq1, and a real function \mathcal W\in C_c^\infty(\R^2), depending only on \varphi and the slopes, with the following property. For an interval I of length r=D2^k and midpoint m_I, the kernel
is admissible for Theorem 2.2, with interior margin \delta_0=1/8 and derivative constant C_0. Moreover,
Proof. The moments of \psi through order three vanish: its jth moment is (2^j-1)\int t^j\varphi(t)\,dt. Choose B\geq1 with \operatorname{supp}\psi\subseteq[-B,B], and define
For t>B, expansion of the cubic and the four zero moments show that g(t)=0; it also vanishes for t<-B. Hence g\in C_c^\infty(\R), g^{(4)}=\psi, and g satisfies (3.2) with a possibly larger constant. Indeed, its first four derivative bounds follow from integration over [-B,B], and its subsequent derivatives are derivatives of \psi.
Take u(z)=8\rho(8z), with the bump constructed above. Choose a dyadic integer D>4B\max_j|b_j|, and set
The constant-coefficient operators commute. On the support of \mathcal W,
After rescaling, every x+b_jt therefore lies in I, at distance greater than r/8 from its endpoints. Put G_j=\prod_{i\ne j}L_iF. The chain rule gives
The function differentiated is compactly supported in t, so its integral is zero. This proves all four line cancellations. The product rule and (3.2) give the required derivative bounds, with fixed D,b_j absorbed into C_0.
Finally, integration in z kills every term of \prod_jL_jF containing a z-derivative. Since \int u=1,
Division by r=D2^k yields (3.4).
We use the block size and variation \mathcal V_B from Section 2. Block endpoints are always deterministic; the input values may depend on the spatial variable.
Lemma 3.3 (A real-line test). Let 0\leq k_-\leq k_+, and partition \{k_-,\ldots,k_+\} into w nonempty consecutive blocks. Suppose that bounded Borel functions h_1,h_2,h_3 and h_{0,k} are supported in a common interval of length L>0, that |h_j|\leq1 for 1\leq j\leq3, and that \mathcal V_{B_0}(h_0)\leq w^{-1/2} on every block B_0. Then
The constant is independent of the scale range, partition, and inputs.
Proof. Put r_k=D2^k and R=r_{k_+}. For 0\leq a<R, use the translated dyadic lattice
At scale r_k, use the kernels from Lemma 3.2. Only roots of length R meeting the common support interval can contribute: the zeroth coordinate of a supported kernel belongs to its interval. These roots have total length at most L+2R.
Restrict all inputs to each such root. This changes no local form, since all four coordinates on its kernel support belong to the summation interval. In the root tree, depth is s=k_+-k. Reversing the scale order preserves the block size and variation. Thus slot zero satisfies the schedule hypothesis of Theorem 2.2, and the three fixed slots satisfy it with one block each. Summing that theorem over the roots gives
For r=r_k, the following sum is r-periodic in a. Since r divides R, unfolding one period yields
As rH_I is the integral against r^{-1}w_I, averaging (3.6) in a gives (3.5). All exchanges of integrals and sums involve finitely many scales, bounded supports, and uniformly finitely many relevant intervals.
3.2. Integer sequences and finite orbit segments
The next step is where primitivity is useful. We do not assume that any pair of slopes has difference one.
Lemma 3.4 (An integer-sequence test). Under the scale and block hypotheses of Lemma 3.3, let h_1,h_2,h_3 and h_{0,k} instead be sequences on \Z, supported in a common integer interval of cardinality L\geq1. Assume the same pointwise size and block variation bounds. Then
Proof. Write E_b=\sum_j|e_j|, and fix 0<\rho<1/4 so small that
Plant every sequence value on the half-open interval [m-\rho,m+\rho), and set it equal to zero elsewhere. The resulting Borel functions have the same size and variation bounds, and their common real support lies in an interval of length at most L-1+2\rho\leq L.
Suppose a configuration lies in four planted intervals, with integer centers p_j, and write \delta_j=x+b_jt-p_j, so |\delta_j|<\rho away from null boundaries. Set
By (2.5), \eta=\sum_je_j\delta_j. Consequently
The expression inside the absolute value is an integer, so p_j=m+b_jn for all four slots. Conversely, a configuration lies in the intervals with these centers precisely when
up to boundaries. This polygon has positive finite area A_\rho: it contains a neighborhood of zero and is bounded because two distinct coordinate forms are linearly independent. Distinct integer pairs (m,n) give disjoint regions up to boundaries, by uniqueness of the centers and (2.5). On P_\rho, |\eta|\leq\rho E_b.
It follows that the real-line integral at scale k is exactly
Let D_k be the discrete sum at scale k in (3.7). The derivative bound \|\Delta_k'\|_\infty=2^{-2k}\|\psi'\|_\infty gives
If this difference is nonzero, then |n|\leq B2^k+\rho E_b. Since k\geq0, there are O_{\varphi,b}(2^k) such integers n, and at most L contributing integers m. Every coefficient product has absolute value at most one. Hence
Sum over k, use \sum_{k\geq0}2^{-k}=2, and apply Lemma 3.3 to the planted inputs. Division by A_\rho>0 proves (3.7).
For measurable functions |f_j|\leq1 on the probability space, put
All sums are finite. If \operatorname{supp}\varphi\subseteq[-B_\varphi,B_\varphi], with B_\varphi\geq1, then
Proposition 3.5 (Finite orbit oscillation). For every deterministic list 0\leq k_1<\cdots<k_{w+1},
This holds on every invertible probability preserving system, with a constant independent of that system, the functions, and the endpoints.
Proof. First take sequences h_j, bounded by one and supported in an integer interval \mathcal K of cardinality L, and set
For m in a subset \mathcal E\subseteq\mathcal K, let \ell_i(m) be the smallest maximizing index of |\mathsf V_\ell(m)-\mathsf V_{k_i}(m)| in the finite range k_i\leq\ell\leq k_{i+1}. Write d_i(m)=\mathsf V_{\ell_i(m)}(m)-\mathsf V_{k_i}(m), and take \alpha_i(m)=\overline{d_i(m)}/|d_i(m)| when d_i(m)\ne0, and zero otherwise. In the block k_i\leq k<k_{i+1}, set
On each block this sequence has size at most (4\sqrt w)^{-1} and at most one jump of that size. Its block size and variation is therefore at most w^{-1/2}. The block endpoints remain deterministic, although the prefix and phase depend on m.
The identity \Delta_k=\varphi_{2^{k+1}}-\varphi_{2^k} telescopes, so Lemma 3.4 gives
Indeed, the sum tested in that lemma is exactly (4\sqrt w)^{-1} times the nonnegative left side, since \alpha_i(m)d_i(m)=|d_i(m)|.
We now use finite orbit segments, the elementary form of Calderón’s transference principle [Cal]. Fix the endpoints, and choose an integer
For each M\geq1 and x\in X, use the sequences
Their common support has cardinality M+2Q. For 1\leq m\leq M and all scales in the endpoint list, no term of the average is cut off, so \mathsf V_k(m)=V_k^\varphi(S^mx). Apply (3.12) with \mathcal E=\{1,\ldots,M\}, integrate over x, and use invariance under S^m. This gives
The maxima are measurable finite maxima. Divide by M and let M\to\infty, keeping the endpoints and hence Q fixed. This proves (3.11).
Lemma 3.6 (Oscillation implies convergence). Let (Z_k)_{k\geq0} be a uniformly bounded sequence of measurable complex functions on a probability space. Suppose that, for every deterministic 0\leq k_1<\cdots<k_{w+1},
Then Z_k converges on a measurable conull set.
Proof. We include the standard finite-block argument; compare the oscillation methods in [Bou, JKRW]. The measurable function
vanishes exactly where the sequence is Cauchy. If its positive set has positive measure, there are \delta>0 and E=\{\omega>2\delta\} with \mu(E)=d>0. For every K,
By monotone convergence of the finite maxima, a deterministic L>K satisfies \int\max_{K\leq\ell\leq L}|Z_\ell-Z_K|\,d\mu\geq\delta d/2. Starting with k_1=0, choose successive endpoints this way. The assumed estimate would give w\delta d/2\leq C\sqrt w for every w, a contradiction. Thus \omega=0 almost everywhere, and completeness of \C proves the claim.
Proposition 3.5, (3.10), and Lemma 3.6 already give almost everywhere convergence of V_k^\varphi. Mixing has not yet been used.
3.3. Identifying the limit by mixing
We include the norm argument to make explicit that it accommodates arbitrary distinct nonzero slopes and both signs of time. The set-mixing assumption on S implies, by uniform approximation by simple functions, that for bounded measurable g,h,
Lemma 3.7 (The Hilbert space van der Corput criterion). Let (u_n)_{n\geq1} be bounded vectors in a complex Hilbert space, with inner product linear in the first variable. If the limits
exist for each h\geq1 and tend to zero as h\to\infty, then N^{-1}\sum_{n=1}^Nu_n\to0 in norm.
Proof. We record the finite estimate underlying the criterion; compare [HK, Lemma D.1]. Extend u_n by zero outside \{1,\ldots,N\}. The identity
and Cauchy’s inequality give
If \|u_n\|\leq M, first let N\to\infty with H fixed, obtaining
The right side tends to zero with H, since the Cesàro means of (|\gamma_h|) tend to zero.
Lemma 3.8 (Norm convergence at distinct slopes). If S is mixing, a_1,\ldots,a_r are distinct nonzero integers, and g_1,\ldots,g_r are bounded measurable functions, then
Consequently, for F_n in (3.9) and p=\prod_{j=1}^3\int_Xf_j\,d\mu, every \chi\in C_c^1(\R) of integral one satisfies
Proof. Write U^mf=f\circ S^m. These operators are unitary on L^2(\mu), with inner product \langle f,g\rangle=\int f\overline g\,d\mu. We induct on r, starting with the empty product r=0. Suppose first that one of the g_j has mean zero. Put
By invariance,
For fixed h, induction applies to the distinct nonzero slopes a_j-a_1, 2\leq j\leq r. The average of this correlation therefore tends to
Here removing the final h terms and dividing by N instead of N-h does not change the limit. By (3.13), the jth factor tends to |\int g_j\,d\mu|^2 as h\to\infty. At least one limit is zero, so \gamma_h\to0. Lemma 3.7 proves (3.14) in the centered case. For general inputs, write each g_j as its mean plus a mean-zero function and expand. The constant term is the asserted limit, and every other term contains a centered factor, so the centered case applies. This completes the induction. Only fixed h is used in the induction; no estimate uniform in h is needed.
Apply (3.14) to b_j and -b_j. The partial sums
Choose B\geq1 containing the support of \chi, and put M_a=\lfloor Ba\rfloor+1, w_n^\pm=a^{-1}\chi(\pm n/a). Then w_{M_a}^\pm=0, and summation by parts gives
Its norm is at most a^{-2}\|\chi'\|_\infty\sum_{n<M_a}\|R_n^\pm\|_2=o(1). To see the last assertion, split at a fixed N_0, and use \|R_n^\pm\|_2\leq\epsilon n beyond N_0; the tail is O_B(\epsilon), and the finite initial part tends to zero. The n=0 term is O(a^{-1}) in norm, while the Riemann sum a^{-1}\sum_n\chi(n/a) tends to one. Adding the two signs and the zero term proves (3.15).
For each averaging profile, Lemma 3.8 gives V_k^\varphi\to p in L^2. Its almost everywhere limit is therefore p. More explicitly, if E_\varphi is the measurable conull convergence set, Fatou’s lemma gives
3.4. Removing the smooth weights
The moment restrictions need not obstruct approximation in L^1. The correction below, from [V, Lemma 2.8], moves a small amount of signed mass to a large scale. Uniform control of the supports or derivative constants is neither asserted nor needed.
Lemma 3.9 (Profiles approximate interval indicators). For every c\in[1,2] and \epsilon>0, there is an averaging profile \varphi with
Proof. Let \rho be the mass-one bump constructed above, and set \rho_\delta(t)=\delta^{-1}\rho(t/\delta). For sufficiently small \delta>0, the real smooth function v=h_c*\rho_\delta has mass one and \|v-h_c\|_1<\epsilon/2. Indeed,
and averaging this identity against \rho_\delta proves the claim. Also \|v^{(n)}\|_\infty\leq\delta^{-n-1}\|\rho^{(n)}\|_\infty, so v satisfies (3.2). Fix v, and write m_j=\int t^jv(t)\,dt for 1\leq j\leq3.
Choose four distinct real numbers z_0,z_1,z_2,z_3, and define
The jth moment of \rho_{i,A} is A^jP_j(z_i). Since P_j is monic of degree j, the fixed matrix \mathsf M=(P_j(z_i))_{j,i=0}^3 is an invertible triangular change of the Vandermonde matrix (z_i^j)_{j,i=0}^3. Solve
The inverse matrix is fixed, hence \sum_i|\beta_i|=O_v(A^{-1}). Thus
has mass one, its first three moments vanish exactly, and its added L^1 error is O_v(A^{-1}). Take A large enough that this error is less than \epsilon/2. Finite addition and fixed dilations preserve (3.2), so this real compactly supported function is the required profile.
Completion of the proof of Theorem 1.1. Continue with |f_j|\leq1, and write
For every rational c\in[1,2] and integer h\geq1, choose an averaging profile \varphi_{c,h} with \|h_c-\varphi_{c,h}\|_1<2^{-h}. The preceding results give a measurable conull set on which V_k^{\varphi_{c,h}}\to p. Intersect these sets over the countable family, obtaining a conull set E.
For N_k=\lfloor c2^k\rfloor, the choice of endpoints in h_c gives the exact identity
Since |F_n|\leq1, for every profile \varphi,
The second term tends to \|h_c-\varphi\|_1: its summand is a bounded compactly supported Riemann-integrable function, with possible discontinuities only at 0,c. Apply (3.16) with \varphi_{c,h}, then let h\to\infty. For every x\in E,
Fix an integer s\geq1, and use the finite rational grid c_i=1+i/s, 0\leq i\leq s. If 2^k\leq N<2^{k+1}, choose the largest i such that n=\lfloor c_i2^k\rfloor\leq N. Then 0\leq N-n\leq2^k/s+1. Averages of terms bounded by one satisfy, whenever 1\leq n\leq N,
By (3.17), the maximum error over the finite grid tends to zero for x\in E. Consequently
Let s\to\infty to obtain convergence along all positive integers.
Finally undo the normalization of the functions and use S=T^g, b_j=\lambda_j/g, as in the initial reduction. If the functions are specified only as essentially bounded classes, choose bounded representatives and remove the union of all integer translates of their exceptional sets. This is still a measurable null set. Every other exceptional set used above is measurable and every intersection is countable. Thus the proof needs no standardness or completeness assumption on the probability space.
4. Quadratic structure and localization for the new slopes
We begin the proof of Theorem 2.2 by establishing the quadratic detection and localization estimates for the configuration y_j=x+b_jt. These estimates supply the input to the short-block and frequency arguments of Section 5. The relevant constructions are those of [H, Sections 3–6] and [V, Sections 4–6]. We retain their interval projections and stopping constructions, and prove below the coefficient-dependent identities and estimates on which those constructions rely.
Throughout this section, the slopes are fixed. Constants may depend on them and on the admissibility parameters. A rational number has height at most H if its numerator and positive denominator in reduced form are at most H in absolute value; matrix heights are entrywise. We use the quantitative convention of [H, Section 3.1]: a quasipolynomial bound in parameters H_1,\ldots,H_m\ge2 is a bound of the form
The degree C is chosen independently of q, the number of depths, and the input schedules. A fixed number of such constructions may be composed. Constructions repeated a variable number of times retain the explicit count estimates in H and V.
4.1. Local bounds and the chart class
For a scalar or Hilbert-valued function, write \norm z_{p,I}=(|I|^{-1}\int_I\norm{z(y)}^p\,dy)^{1/p}, with the essential-supremum convention for p=\infty. Norms without an interval subscript use ordinary Lebesgue measure. The elementary bound used throughout the argument is unchanged except for its constant.
Lemma 4.1 (Local absolute bound). If 1\le p_j\le\infty and \sum_jp_j^{-1}\le2, then
The assertion also holds for Hilbert-valued inputs, using their norms in the integrand.
Proof. For two distinct roles i,j, the map (x,t)\mapsto(y_i,y_j) has determinant b_j-b_i\ne0. This proves the endpoint with p_i=p_j=1 and the other exponents infinite. The support and kernel bounds give the all-infinite endpoint, and normalized norm monotonicity gives the endpoints with only one exponent equal to one. Log-convexity of the positive integral fills the stated polytope, as in [H, Lemma 2.1].
On the support of a local kernel, all four coordinates lie in an interval of length r=|I|. Consequently
In any pair of coordinate variables, the rescaled kernel thus has support in a fixed box and the same type of derivative bounds. Its integral in either coordinate with the other fixed is zero, by the corresponding one-coordinate marginal hypothesis. These facts also give C_N(q)\lesssim_b N and preserve the normalization-removal argument of [H, Lemma 2.4].
We recall the concrete class in which the quadratic structure is detected. Write e(u)=\exp(2\pi i u). A quadratic chart atom is a function
where \theta\in\R^d, \lambda,\beta\in\R, C\in\mathbb Q^{d\times d}, and \chi is smooth and compactly supported. An overall constant of modulus one is allowed. We use exactly the complexity convention of [H, Definition 4.1]: for complexity H\ge e^2, d\le\log H, the height of C is at most H, \operatorname{supp}\chi\subset[-H,H]^d, and
with a fixed chart constant c_{\rm ch}. The real oscillation parameters \lambda,\beta,\theta are unrestricted. A bounded atom has supremum norm at most one; any normalization factor is put in its external scalar coefficient.
The elementary operations, Weyl alternative, and chart expansion of [H, Lemmas 4.2–4.4] concern this class itself. Their proofs therefore apply without a coefficient change. We will need the precise splitting supplied by the Weyl alternative. In its multivariable version replace \theta y by Vz, and the ordinary quadratic phase by z^T\Lambda z, where \Lambda is real symmetric and z has at most three coordinates. Put
Let Q(z) be this multivariable chart. Suppose 0<\delta<1/2 and a Lipschitz function \vartheta, supported in a fixed box with supremum and Lipschitz norms at most H, satisfies
The Weyl alternative then supplies a rational subspace U, the orthogonal projection P_f onto U, P_s=1-P_f, and the bounds
The projection heights are at most H', and H' is quasipolynomial in H,\delta^{-1}. The subsequent chart expansion has the absolute coefficient and pooled Fourier-tail bounds of [H, Lemma 4.4]. These are the bounds used below, including when expansions must be summed over depths.
4.2. Separation along the configuration and local detection
The first coefficient-dependent step moves two structured inputs into the other two coordinates and the step variable. This is the form of separation needed for both detection and the later predictor estimates.
Lemma 4.2 (Progression separation). Let i,k,a,l be the four distinct roles and let \phi_i,\phi_k be bounded chart atoms of complexity at most H. On y_j=x+b_jt, their product has an absolutely convergent expansion
with bounded chart atoms as factors and the same shell bounds as [H, Lemma 4.5], with constants depending on b. In particular, for a quasipolynomial D=D_b(H), a dyadic Fourier shell of size T\ge1 has absolute coefficient sum at most DT^{-100}, and its factors have complexity at most (DT)^C.
Proof. First separate \phi_i, temporarily suppressing its linear phase. Choose the fixed rationals
They satisfy \sum_j\alpha_j=\sum_jb_j\alpha_j=0. Insert the smooth integer partitions of unity from H at \theta y_a,\theta y_l,\theta t, and write the lifts as d_j=\theta y_j-m_j and d'=\theta t-m. For j=i,l, set
These carries are integer vectors. The support restrictions bound them by the chart budget, with a fixed slope factor. Conversely, k_i,m_a,m determine m_i uniquely. Each remaining specified carry equality is imposed exactly by a smooth cutoff equal to one at the specified integer vector and zero at all other integer vectors.
For S_\alpha=\sum_j\alpha_jb_j^2, direct expansion gives
Indeed, substituting d_j=d_a+(b_j-b_a)d'+k_j cancels the common-lift and x-terms. Here S_\alpha=(b_i-b_a)(b_i-b_l). Solving for the phase in role i leaves chart matrices -\alpha_aC,-\alpha_lC,S_\alpha C in the three desired variables. For fixed carries, the last sum is linear in x,t, and hence splits between y_a,y_l; so does the suppressed linear phase.
The coupled amplitude is smooth in d_a,d_l,d'. Its Fourier expansion on a larger fixed box separates those variables, exactly as in the proof of [H, Lemma 4.5]. The new substitutions cost at most C_b^n in derivatives of order n, and their rational heights grow by fixed factors. The possible integer choices increase by at most C_b^{O(d)}. Since d\le\log H, these changes are within the quasipolynomial chart budget. H’s integration by parts, with derivative order exceeding the Fourier dimension by at least 110, gives the stated absolute shell bound. Separate the second atom with \alpha_i=0 and combine the factors by direct sums of horizontal coordinates. This completes the expansion.
For one separated term, let
where \zeta is smooth, nonnegative, equal to one on [-1/3,1/3], and has \norm\zeta_1\le1, as in [V, Section 5.4]. The two-input form in the remaining roles has norm at most C_bA_r on normalized L^2(I) inputs. To see this, use
Fourier-expand the smooth normalized kernel in those two coordinates. Each term is a modulated convolution with a fixed dilation of Q(t)\zeta(t/r), so Plancherel gives the claimed bound. If A_r\ge\delta, the Weyl alternative and chart expansion further give the finite smooth separation in [H, Lemma 4.6]: for every 0<\varepsilon<1/2, the two-coordinate kernel has a separated approximation with uniform error at most \varepsilon, whose number of terms, atom complexities, and absolute coefficient sum are quasipolynomial in H,\delta^{-1},\varepsilon^{-1}. The retained phases are \lambda_*t^2+\gamma t, with |\lambda_*|r^2 bounded by the stated budget, and the linear phase splits as
Thus the fixed rational-height description of the effective quadratic and linear frequencies in that lemma is preserved.
Lemma 4.3 (Local detection). Fix q>2 and C_1\ge1. Suppose that the four inputs R,h,u,v, in any slot order, satisfy
where 0<\delta<1/2. There is a bounded chart atom \phi of complexity at most D such that
The degree C is independent of q.
Proof. We verify the progression estimate in [H, Lemma 4.11] and then apply its detection proof. Normalize I=[0,1). Extend the kernel to the same fixed circle used there. By (4.1), its real support fits in the periodic (x,t)-box, and all real y_j on that support already belong to I. Hence zero-extended inputs on the circle give the correct original integral, even when some slopes are large or negative. The kernel’s absolute Fourier coefficient sum is uniformly bounded.
The integer left inverse in Lemma 2.3 distributes a circle character with indices (p,q') into slot characters with indices p\mathbf1_{j=0}+q'e_j. To a character-weighted progression average apply H’s three Cauchy–Schwarz steps, taking Haar differences in the directions (-b_j,1). The last cube has three independent increments multiplied by nonzero differences b_i-b_j. Each such integer multiplication pushes circle Haar measure to itself. The last integral is therefore the eighth power of the same circle U^3 norm as in H. Character multiplication preserves that norm, the L^2 norm, and the supremum norm. This proves H’s local generalized von Neumann estimate; its accompanying two-L^1 bound follows from Lemma 4.1.
The continuous inverse theorem and bounded-input decomposition in [H, Sections 4.5–4.6] concern functions on this circle and its chart atoms, and require no change. In the proof of [H, Lemma 4.12], clip the three opposing inputs at height B. The local absolute bound gives total error at most C_bC_1B^{-(q-2)/q}. Taking B of order (C_bC_1/\delta)^{q/(q-2)}, apply that bounded-input decomposition twice. The error estimates use the generalized von Neumann and two-L^1 bounds just proved. Two bounded atoms remain opposite R and an input of normalized L^2 norm at most one. Lemma 4.2 and the preceding single-scale separation give a finite sum of products of one-coordinate moments, with controlled absolute coefficient sum. One moment of R consequently has the asserted size. These are a fixed number of quasipolynomial constructions. The clipping exponent contributes only a q-dependent multiplier to a logarithm, so the final logarithmic degree is independent of q, as required.
4.3. The base reduction and inputs that vary with depth
Detection is used in H and V to construct finite approximation spaces on interval trees. We recall their meaning before verifying the scale estimates. A column with label v, born on an interval J_v, is an augmented function
where \phi_v is a scalar chart atom, \epsilon_v>0, and the \mathbf u_v are distinct orthonormal private coordinates. Only the first, physical coordinate enters H_I. A local projection on I is the orthogonal projection onto constant-coefficient combinations of the available columns restricted to I. Lists are inherited on descendants. Accordingly the direct-sum projections on successively finer interval partitions have nested ranges. The private-coordinate, compression, and height-bin estimates of [H, Section 3] depend only on these definitions and the column budgets, and apply unchanged.
There are two coefficient-dependent scale estimates in the base construction. For two fixed chart atoms, the separation above reduces the other two inputs to a kernel of the form w_I(x,t)Q(t), with fixed bounded one-coordinate multipliers. The proof of [H, Lemma 5.1 and Corollary 5.2] expands Q, groups nearby linear frequencies, and uses the zero marginals of w_I. In coordinates (y_a,y_l), its double primitive has compact support in an interior rectangle. Fourier expansion followed by one derivative in each coordinate expresses it as a sum of tensor products of mean-zero smooth bumps. In role p\in\{a,l\}, the resulting scalar rows have fixed bounded carrier functions
where \psi_p is the fixed separated multiplier in that role. These are the only coefficient replacements in the row estimate. They alter frequency separations and rational heights by fixed factors. H’s height-bin argument, ordinary Bessel estimate, and compression step therefore give its absolute scale sum and finite predictor list with the same quantitative orders, including the pooled-tail bound needed for summation over all depths. The finer partition comparison in [H, Lemma 5.3] uses the resulting one-coordinate predictor moments, so its projection-increment proof applies with the same replacement.
For the second estimate, fix three bounded chart columns and take the transpose in role l at y+h. Its other coordinates are
After squaring the normalized L^2 norm, the six atoms form a multivariable chart in (h,t,t'). At a fixed path point y, all parameters before multiplication by r or r^2 remain fixed as the depth changes. H’s Weyl alternative, chart expansion, and height-bin argument therefore apply. Away from the bounded number of intermediate scales, the bounded effective quadratic parameters are tiny and each retained t-frequency is tiny or large. Smoothness handles the large frequencies. Removing the tiny oscillations leaves exactly
This proves the active-scale count and the small-norm scale sum in [H, Lemma 5.7], with constants depending on b.
For scheduled inputs we use [V, Section 4]. Its row estimates concern Hilbert-space operators on refining interval partitions, evaluated on inputs at prescribed nondecreasing read depths. Their hypotheses are support on the input partition, representation of the adjoint ranges on the observation partition by the inherited columns, the stated order of successive partitions, and an ordinary Bessel bound before evaluation on the varying inputs. None involves the maps y_j. Restriction to a new spatial root preserves the pointwise block bound in Theorem 2.2. The original W_j and original depth schedules are retained; subsequences may have empty intersections with some blocks, as allowed in V’s row argument. Its nearest-endpoint extension is retained as well.
The smooth rows in [V, Lemma 4.6] require the ordinary Bessel bound just verified, with the fixed carriers above. Thus the scheduled scale-sum and predictor results [V, Lemma 5.3 and Corollary 5.4], and its delayed replacement Lemma 5.5, hold with the new coefficients. The remaining steps of its Section 5 use the same local detection, fixed-cell fits, projection energies, and staggered pair comparisons. A finer projection in a comparison is still applied to the input at the prescribed read depth, and its row estimate has no factor depending on W_j.
We now describe the output of the base construction. For the fixed-input argument, embed the scalar input as F_{j,s}=(f_j,0), independent of s, on a subtree satisfying \norm{f_j}_{q,I}\le1 on every cell, as in [H, Proposition 2.3]. For the scheduled argument, put F_{j,s}=(z_{j,s},0) with the inputs of Theorem 2.2. The adaptive fits and two-earliest-increments expansion of [H, Lemmas 5.4–5.5], or [V, Sections 5.1–5.2], express the original form, at any finite accuracy cutoff, in tuples with two differences of fitted projections and two original-input residuals, together with a remainder whose absolute length-weighted sum tends to zero on the fixed finite tree. Each tuple has an earlier and a later accuracy index, and a choice of the four roles. Their exact staggered ordering is retained from those results. High-position truncation, using the predictor estimates above, is [H, Proposition 5.6] or [V, Sections 5.5–5.6]. It assigns each original summation cell to a segment root R with inherited projection families; write \mathcal I_R for the cells assigned to R. Auxiliary projection cells are not additional summands.
Fix one tuple and let k be its later accuracy logarithm. The two positions called low in this construction are denoted i,j, and the other two by a,l. On a cell of depth s, the resulting arguments satisfy
All expressions are restricted to I; a projection indexed by a finer depth is the direct sum on that finer partition. The Q^{p,h} and P^{\rm new}_p are the inherited basic families just constructed, and K_0 is their positive integer delay. Their column sizes, inverse private weights, scalar complexities, coefficient sums, and path counts have a joint budget e^k\le D_{\rm base}\le\exp(C_{q,b}k^C). On every segment they give
One may take \Lambda_k to dominate C_{q,b}(1+K_0)\log(2D_{\rm base}). For the fixed-input argument the smaller bound in H also suffices.
For noninitial accuracy levels, every pair of full Y arguments has form norm at most \eta=\exp(-k^\theta), and |H_I(Y_I)|\le\eta. Here the pair norm is the supremum against normalized L^2(I) inputs in the complementary two positions; 0<\theta<1 is fixed before q. Define \mathcal A_R to consist of those cells for which some available triple of basic scalar columns has transpose norm greater than D_{\rm base}^{-C_{\rm act}}, with the fixed sufficiently large exponent in H or V. The active-scale estimate gives \sum_{I\in\mathcal A_R}\mathbf1_I\le D_{\rm base}^C.
There are three remaining patterns: z^\varnothing=B, z^a replaces only role a of B by U_a, and z^l replaces only role l by U_l. Put
The output of [H, Proposition 5.8] and [V, Proposition 5.6] is
Indeed the term with both U factors uses the full complementary pair bound and (4.6). The inactive triple scale sum bounds the omitted cells. On the others, apply \min(\eta,\sum x_m)\le\sum\min(\eta,x_m) to the full tuple before separating its patterns. No individual pattern is asserted to have the full-pair smallness.
For each fixed tuple the segment roots have total length bounded by a fixed multiple of the original root length, and the sum of the high-truncation replacement errors is at most C_{q,b}e^{-k} times that length. These are the segment-packing assertions in the cited base constructions. The finitely many initial accuracy levels use the active path bound and local sizes, without the factor \eta, and hence have depth-independent bounds. All these estimates are valid anew on a normalized subroot with its inherited lists and original input schedule. We have therefore reduced both input classes to estimating (4.7) uniformly over these concrete base patterns.
4.4. Curvature tags and separation of distinct tags
We now localize the base arguments according to their quadratic structure. Recall the nonzero integers c_j of Lemma 2.3. A column in role i has normalized curvature tag v=(\theta_v,C_v,\lambda_v) when its quadratic chart phase is written as
Its linear modulation is not part of its tag. This is H’s normalization with (-1,3,-3,1) replaced by c. Division by the fixed nonzero c_i preserves the chart budgets.
At length r, the difference data of two tags v,v' are
The ambient graph in [H, Section 6.1] joins two tags when these data admit (4.3) with bounds at most the chosen edge height H_{\rm ed}; equal tags have distance zero. This graph includes all hypothetical tags within the basic budget. Its distances depend on tags and depth, and do not increase when the interval length is halved.
Lemma 4.4 (Separation of distinct tags). For single-column complexities at most \exp(C_{q,b}k^C) and \log(1/\tau)\le C_{q,b}k^C, the edge height can be chosen with \log H_{\rm ed}\le C_{q,b}k^C so that tags at graph distance greater than one have both of the following properties. Their same-slot normalized physical inner product is at most \tau in absolute value. In distinct slots, their form norm against two arbitrary normalized L^2(I) spectators is at most \tau.
Proof. For the first assertion use the Weyl alternative on the conjugate product, after the smoothing of the interval indicator in [H, Lemma 6.1]. Its fixed nonzero factor c_i changes the allowed rational heights by a fixed amount, and multiplying back gives the splitting for the tag difference.
For distinct source slots i,j, let the spectators be at a,l, and write the centered source coordinates in terms of the normalized spectator coordinates z,w. Their linear parts are r(\xi_pz+\zeta_pw), where
These are nonzero fixed rationals. The cancellation of polynomials of degree at most two gives
Use the difference data (4.9). For a small return V'h=m+d_h, choose an integer period \Pi and require m\in\Pi\Z^{d_v+d_{v'}}, where d_v,d_{v'} are the horizontal dimensions of the two tags. Choose \Pi sufficiently divisible that \xi_i\Pi,\xi_j\Pi are integers and \Pi C_v,\Pi C_{v'},c_i\xi_i\Pi C_v,c_j\xi_j\Pi C_{v'} have integral entries. This clears the lattice-index and phase denominators together. Translate z by h, and the source lift indices by \xi_im_v and \xi_jm_{v'}, respectively. The phase-increment calculation in the proof of [H, Lemma 4.3], together with (4.10), gives the modulation in w
Indeed each source contribution has the multiplier c_p\xi_p\zeta_p; the two contributions therefore combine into these difference data. The remaining linear terms in z are absorbed in the first spectator. The amplitude and kernel errors are bounded by a quasipolynomial times |h|+\norm{d_h}, as in H.
Call the two spectators f,g, and clip them at a sufficiently large inverse polynomial in \tau. Their discarded L^1 tails are small by their L^2 bounds and Lemma 4.1. If the form is larger than \tau, applying the return calculation to sufficiently many multiples of (h,m) gives large correlations of one fixed bounded kernel with tests f_\ell(z)g_\ell(w), where the f_\ell have bounded L^2 norms and
up to a common bounded multiplier. The small return threshold is chosen before these multiples so that all translation errors remain small. No regularity of the first spectator is required: it is translated by a change of variables. If |\Omega(h,m)| exceeds a fixed constant, the second tests are Bessel on the fixed box, with bound controlled by the clipping height. Cauchy–Schwarz in z then contradicts the number of large correlations. Thus \Omega is bounded on every sufficiently small return.
The concluding return-to-splitting argument in the proof of [H, Lemma 4.3] uses precisely this bound: its alternating return identity yields an isotropic lattice space, successive minima bound the rational fast projection and slow part, and coordinate-direction returns bound the effective quadratic form. The enlarged denominator-clearing period changes only its quasipolynomial constants. It supplies (4.3), contradicting the assumed graph distance. This is the distinct-slot argument of [H, Lemma 6.1] for the present coefficients.
4.5. Same-tag tests and the localization output
The complementary fact is that a column of one tag can detect a residual in another role using a template of that same tag. This keeps the approximation spaces finite without repeatedly increasing their nonlinear complexity. A pure source column means one scalar column label, rather than a linear combination of labels.
Lemma 4.5 (Uniform same-tag test). Fix an original tag v and its original source columns within the complexity budget \exp(C_{q,b}k^C). There are fixed template families in all four slots, consisting of bounded columns of tag v with arbitrary external linear modulations. Suppose a pure source column from the original or template family in role i has form norm greater than t_0>0 against a physical residual in L^2(I) in role j and two normalized L^2(I) spectators. The residual then has a normalized moment of magnitude at least t_0/H_0 against a template in role j, where H_0\le\exp(C_{q,b}k^C). This bound is independent of the external linear modulations and of the number of earlier template insertions.
Proof. Use the nonnegative partition function \rho and template amplitudes \rho(d)e(\ell^Td/P_0) of [H, Lemma 6.6], with \ell ranging over residue vectors modulo a sufficiently large integer period P_0. For the four lifts d_p=\theta_vy_p-n_p, set
The Bézout identities show that d_\diamond is a lift of \theta_vt modulo the integer lattice and every k_p is integral. All are bounded on the branch boxes. The common quadratic phases cancel by Lemma 2.3; for fixed carries the remaining phase is
which is linear in the coordinates. Smooth equality cutoffs in the four lifts impose (4.11) exactly, because their arguments are integer valued on the configuration. Their Fourier expansion on a sufficiently large box therefore produces the same template amplitudes. Increasing the period and the possible carry count costs only the allowed chart budget.
As in H, write each Fourier index as \ell=P_0q'+r'. The factor e((q')^Td_p) is an external linear modulation, so only r' enters the template amplitude. This gives the same bound even when the source was itself a previously inserted template. For an original source, the remaining smooth horizontal amplitude, after phase cancellation, is periodized and expanded in ordinary torus characters. Its remaining dependence on the source coordinate is a linear modulation, which is absorbed into the spectator tests. All these absolute coefficient sums are bounded by the same quasipolynomial budget. After separating the smooth kernel, express the target coordinate as y_j=\alpha y_a+\beta y_l, with fixed nonzero rational \alpha,\beta. Fourier-expand the demodulated residual on a fixed larger interval. If M is the supremum of its normalized template moments, one separated term is bounded by
The last inequality is Bessel’s inequality for spectator frequencies with nonzero spacings proportional to \alpha,\beta. Summing the absolute Fourier coefficients gives H_0M, and proves the assertion.
We record how these two tests enter the remainder of the imported construction. This also fixes the meaning of its output used in Section 5. Put P=\lceil\eta^{-\alpha}\rceil, where \alpha>0 will be chosen in the closing argument. The vertices are occurrences of the original columns and are ordered by admission. Added main and helper columns attach to those vertices and have their normalized tags; they create no new graph vertices. H’s short-path, color, and historical-mask constructions [H, Lemmas 6.2–6.5] form groups of these vertices and assign current groups to each summation cell. An actual use is such an occurrence of a current group at a cell. A mask is the subset of available vertices allowed in a specified local projection. A pad enlarges that subset by the prescribed graph radius; a positive smeared filter is the average of P nested pad projections. A history is the ordered product of such filters or their complements used in the binary localization. The number \mathcal L of colors is polynomial in k, and these main histories have O(\log\mathcal L) factors. Each group records a distinguished tag z, the primary winner in H’s color construction. Every vertex admitted to its masks has, already at admission, a graph path to z of length polynomial in P,k. The short-path lemma therefore gives the endpoint splitting (4.3) relative to z, with logarithmic height polynomial in P,k. This is the splitting used for the first matching of a column in Section 5.
We use the precise masks and historical-use convention of those lemmas. Its relevant properties are as follows. Masks are inherited on descendants; at an actual use, enlarging a pad by one includes every neighbor of the preceding mask. An endpoint is usable when its mask is nonempty and its cell has a downstream actual use of its group. Distinct groups of one color have tags separated by graph distance greater than one on overlapping usable endpoints. The endpoint bound along short paths follows from the rational splitting (4.3). The usable-space Bessel bound follows from same-slot separation in Lemma 4.4. These proofs involve only rational linear algebra, inherited interval projections, and the indicated splitting and separation bounds; thus their path, rank, and pad counts keep the quantitative orders in H.
For completeness, the closure condition imposed on these spaces has an exact form. At rank r, let n_{r-1} bound the number of labels through the preceding rank along a path, let D=D_{\rm base}, and let \tau_3=D^{-C_\tau}. For a basic argument a, each prescribed history \mathcal H, target mask \mathcal M, and pure lower-rank source whose tag belongs to that mask, the form norm against the residual (1-P_{r,\mathcal M})\mathcal Ha and two normalized L^2 spectators is required to be at most
On failure, Lemma 4.5 supplies a new column and hence a projection increment. The products-of-prefixes estimate of H Section 3 bounds the sum of those increments, exactly as in its Lemma 6.7. Its rank recurrence and main/helper catalog counts are therefore unchanged. This count depends on the names of masks and histories, including anticipated helper bases, but not on the numerical edge height or on later mode precision. As in [H, Section 6.4] and [V, Section 6.1], first fix these counts, then the mismatch tolerance, and finally the edge height needed in Lemma 4.4. The lower-rank count in (4.12) is retained, so no full source count is lost when these errors are summed.
It follows that the binary localization and outgoing-filter identities of [H, Proposition 6.8 and Lemma 6.9] hold for the present forms: they use (4.12), distinct-tag form separation, the local absolute bound, and the same nested projection algebra. In the color construction of [H, Lemma 6.3], the first failure of its marking conditions may remove a subtree. Increasing the fixed thresholds makes the sum of these roots’ lengths an arbitrarily small fixed fraction of |R|. Section 6 restarts the original summands on those roots. After these marking-failure cuts, let \mathcal A_R^{\rm main} denote the remaining active cells. For each base pattern and I\in\mathcal A_R^{\rm main} of depth s, the localization is
Here \lambda ranges over polynomially many fixed binary search leaves. The inner history \mathcal B^\lambda_{p,s} is independent of G, and \overline L^\lambda_{p,G,s} is the outer ball average. The minimum in (4.7) is distributed across the fixed leaves, while the group sum remains inside the original cell contribution. The filters act on the whole basic arguments in this identity. We next record their endpoint bounds, then the algebra that moves the complementary filters onto a side operand.
For varying inputs, this localization acts at each cell on its one actual tuple at that cell’s original depth. The insertion tests freeze that operand while the column list grows. The products-of-prefixes estimate therefore remains a fixed-vector estimate, and gives the same closure count; this is the argument of [V, Section 6.1]. After the operator families have been fixed, its Lemma 6.1 and Corollary 6.2 control their changes between endpoint depths. More explicitly, for one slot, search leaf, and color, write the full localized endpoint output as Z_{G,s}=\overline L_{G,s}\mathcal B_s B_s, where \overline L is the outer positive ball average and \mathcal B its group-independent inner history. If depth intervals [r,t] have edge-overlap multiplicity at most a, and E_{G,r,t} is a union of t-cells with downstream actual uses certifying the indicated group, those results give
The shared histories without an outer ball satisfy the corresponding estimate counted once across groups. Endpoints before the attempt root are zero, and all other endpoints use inherited continuations and local refits. The proof first establishes the ordinary Bessel bounds for operator changes, using usable endpoint ranges and compression, and then evaluates the rows at their prescribed depth inputs. Input changes use only the original block variation bound. Thus the constants in (4.14) are independent of the partitions and W_j.
4.6. Transferring the complementary filters
The side estimate will encounter two distinguished operands while the other two retain their whole localized histories. We state the projection identity for this situation explicitly. Fix an actual group use at a cell I, one search leaf, and four distinct roles a,a',c,d. Let S_a,S_{a'} be fixed Hilbert-valued L^2(I) operands in the first two roles. They are called side operands here because they are retained while the complementary histories are simplified; their particular construction and cross-depth bounds are given in Section 5. The other two roles have the whole localized arguments with basic inputs B_c,B_d from (4.13). Choose the side operands before closing the helper catalog. The helper version of (4.12) is imposed with these fixed states as bases and lower-rank main columns as sources, as in [H, Lemma 6.7].
Expand each complementary factor 1-\overline Q as the difference of the identity and a positive average. There are only 2^{O(\log\mathcal L)} resulting products. In one such product, remove the largest-rank outermost positive filter in roles c,d. The positive main ranks are distinct in the source construction, so the remaining source in the other complementary role has strictly lower rank. Suppose the filter being removed is on role c. Let Q_{c,p}, 1\le p\le P, be its individual pad projections, and Q_{a,p} the helper-rank projections in role a with the same masks and pads. For j=a,c, set
The A_j^p are mutually orthogonal projections and sum to the identity. Their final band is included, with weight zero in the positive average.
Let x_c be the inner history on B_c after deleting the outgoing filter, and let x_d be the entire remaining history on B_d. Write x_a=\mathcal D_aS_a, where \mathcal D_a is the product of helper averages already transferred, and x_{a'}=S_{a'}. In the next displays arguments are written in role order a,a',c,d; each still occupies its original slot in H_I. Multilinearity and (4.15) give
Indeed the coefficient of the pair of bands is the weight of band q in role c minus the weight of band p in role a. This proves H’s Equation (15), including its terminal bands, for the present form without a coefficient change.
The estimate accompanying [H, Lemma 6.9] permits only |p-q|\le2 to be retained: the other terms at this cell and group have total absolute value at most
for any prescribed fixed C_3, after the numerical accuracies are chosen. Its hypotheses are precisely the closure and separation estimates proved above. If p>q+2, expand the lower c-band into pure main sources and use the helper closure test on the higher a-band, written as a difference of residual projections. If q>p+2, expand the lower-rank source in role d. A source tag within one edge of the lower helper mask belongs, by neighbor closure, to the next outgoing mask; the main closure test then controls the higher c-band. For a tag farther away, expand the helper primal and use Lemma 4.4. The coefficient counts are canceled by the lower-rank denominator in (4.12), except for fixed count powers already included when the mismatch tolerance was chosen. The O(P^2) band pairs and polynomially many transfer positions fit those accuracy choices because \log P=O(k).
For each fixed offset h=q-p\in\{-2,-1,0,1,2\}, retain the coupled sum
The sum over p stays inside the cell contribution. When the current group is also summed, its sum stays there as well. This is the expression to which the short-block estimate is applied; taking absolute values of its individual bands would lose the matched-band summation used in that estimate.
Forward only the transferred main term in (4.16) to the next step. Removing the remaining positive filters in decreasing rank ends at
Restoring B_c,B_d to the full Y_c,Y_d uses their pair bound from the base reduction. Each correction has an additional basic U factor and is therefore a three-side term.
For the retained defects, the endpoint convention matters. If the outgoing filter is a ball filter, split A_c^{P+1}=1-Q_{c,P} into the shared identity and its positive ball endpoint. The other complementary band endpoints have usable ball ranges. If the outgoing filter is cumulative, its outer ball has already been removed and the entire complementary band expression is shared across groups. Role d likewise either has one positive outer ball filter or is shared. In the endpoint energy estimates, a shared identity is evaluated on the union of required cells and counted once; the original form still retains every actual use. These are exactly the endpoint classes controlled by (4.14) and [V, Corollary 6.2]; they give the side estimate in Section 5 its coupled band operands.
We have obtained the base minimum, the localized group expressions, their outgoing-filter identity, and the endpoint estimates with the same main-column budgets as H and V. The remaining work is to estimate the coupled near bands and the remaining localized quadratic contributions. Section 5 establishes the coefficient-dependent comparisons used for those two estimates.
5. Short blocks, rational shifts, and the core estimate
Starting from the localization proved in the preceding section, we now split the arguments into side terms and core terms. Side terms are controlled by square estimates for projection increments, together with a comparison on a bounded block of depths. Core terms are controlled by separation of the frequencies left after removing a common quadratic phase. We prove the coefficient-dependent parts of both arguments here. The projection, graph, and amplitude summations that we use from [H, V] are identified at their points of application.
Throughout this section, the primitive integer vector b=(0,b_1,b_2,b_3), the integers e_j, and the nonzero integers c_j are those of Lemma 2.3. Constants may depend on these fixed integers and on the admissible kernel bounds. We write e(t)=\exp(2\pi i t). For an interval I, the notation \|f\|_{2,I} denotes (|I|^{-1}\int_I|f|^2)^{1/2}, also for Hilbert-space-valued functions.
5.1. A comparison on a short block
We first establish the version of [H, Lemma 7.2] needed by the side argument. This estimate compares randomized inputs with the fixed-input constant C_N(q); in particular, the exponent of the block length must tend to zero as q tends to two.
Proposition 5.1 (Short-block comparison). Let 2<q<3. Let a finite dyadic tree piece have root R, with internal summation intervals in at most K consecutive depths and within the original range of at most N depths. Its terminal intervals, called leaves, partition R and occur in this depth block or one depth below it. In role j, let \mathbf z_j be a deterministic Hilbert-space-valued function supported on R, with \|\mathbf z_j\|_{2,L}\le s_j on every leaf L. Suppose scalar randomizations z_j^\sigma satisfy
for a=2,4,8, with constants independent of the tree and K. For measurable signs |\epsilon_I(\sigma)|\le1 and any subcollection of internal intervals,
The constant C_{q,b} may also depend on the fixed moment bounds C_a and the kernel bounds.
Proof. Homogeneity reduces to s_j=1; a zero size makes the form zero. Multiplying all signs by one random scalar of modulus one reduces the expected absolute value to the absolute value of an expected signed sum.
We use the partition of the four leaf depths from the proof of [H, Lemma 7.2]. Here is the information from that partition that will be needed. If the least depth occurs more than once, designate two roles attaining it. If it occurs once, designate that role and separate the least and second least depths by the smallest dyadic interval of depth indices containing them in opposite halves. Specify the nonempty set of roles whose depths belong to its upper half, and designate one of those roles by a fixed ordering. Every role that is neither the minimum role nor in that specified set has depth greater than the upper endpoint of that dyadic interval. This partitions the depth tuples into rectangles of masks E_j, each a union of leaves. There are O(1+\log(K+2)) length classes and a bounded number of role assignments. Within each class, there are O(K+1) positions; the masks in each of the two designated roles, say a,a', are disjoint over those positions. For a fixed rectangle there is an index t such that every role has leaf depth at least t, with the possible exception of role a. Only summation depths s<t contribute. Indeed a leaf meeting an internal interval must be its strict descendant. The same assertion holds after the small coordinate enlargements below, by the kernel’s fixed interior support margin.
Write r_s for the length at depth s, and set \delta=\beta r_t, where \beta>0 is a sufficiently small fixed dyadic number depending on b and the support margin. Decompose the physical coordinates as
The slow coordinates n,m are integers. Choose the fast representatives from the parallelogram
This is a fundamental domain for translation by \mathbb Z^2, of area one. To decompose a point, first choose n so that its first coordinate minus n belongs to [0,1), and then choose m so that u+b_1v belongs to [0,|b_1|). The latter choices differ by b_1, so existence and uniqueness hold for either sign of b_1.
Let \mathbb T_0=\mathbb R/(T_0\mathbb Z), with Haar probability measure, where T_0 will be a sufficiently large fixed number. Define F_j^\sigma(p,w) by sampling z_j^\sigma(\delta(p+w))\mathbf1_{E_j}(\delta(p+w)) on a bounded interval of real lifts and setting it to zero elsewhere on \mathbb T_0. Use lift intervals [0,1) and [0,|b_1|) in roles zero and one, and bounded intervals containing the ranges of w_2,w_3 on \mathcal P_b in the other roles.
The fast Haar integral reproduces integration over \mathcal P_b, up to its fixed normalization. To check the potentially noninjective pair coordinates, suppose a torus progression has all four bounded lifts \widetilde w_j in these windows. Put
The integer inverse of Lemma 2.3 gives \widetilde w_j=\widetilde u+b_j\widetilde v\pmod {T_0}. All the discrepancies are bounded in terms of the fixed windows and b,e. Choose T_0 larger than these bounds. Each discrepancy then vanishes, so these are equalities of real numbers and (\widetilde u,\widetilde v)\in\mathcal P_b. Conversely, every point of the tile gives exactly these lifts. This verifies the sampling identity without asserting that any coordinate pair is unimodular.
Replace the kernel by its value at (\delta n,\delta m). At depth s, the error is O_b(\delta/r_s), with a fixed support enlargement. The normalized sampled form has factor (\delta/r_s)^2. On the product of local counting measure and the fast torus, any two distinct coordinate maps are injective on the lattice and Haar-surjective on the torus. Consequently the two- L^1, two-L^\infty bounds hold, and interpolation gives the absolute form estimate whenever the four reciprocal exponents sum to at most two.
For clarity, let M_{j,I} be the squared mass of \mathbf z_j in all masked sampling windows relevant to I, including their bounded overlap multiplicity. Choosing \beta small enough puts these windows inside I. The leaf normalization therefore gives
The absolute L^2 estimate and the moment hypothesis bound the sampling errors at depth s by
The sum over s<t is bounded, and summing rectangle positions uses the disjointness of both designated masks.
We next improve the integrability of the fast inputs. On their fast Fourier coefficients define
The fast progression integral is unchanged under simultaneous application of these unitaries. Its Fourier support obeys \sum k_j=\sum b_jk_j=0, so Lemma 2.3 gives k_j=c_j(h+b_jh') for real h,h'. Thus
For a single role, integrating the fourth power first in the fast variable and then in \gamma equates the sums and the sums of squares of two pairs of Fourier indices. Those pairs agree as multisets, whence
Fourier polynomial approximation is justified by the pairwise L^2 form estimate and this same fourth-moment estimate for differences.
Let \mathfrak e_j(p) now denote the fast L^2 size of the deterministic stack in the sampling window. Window overlap gives \delta\sum_p \mathfrak e_j(p)^2\le C_b|E_j|. If all leaves of a role have depth at least t, then \mathfrak e_j(p)\le C_b: disjoint leaves meeting a window of length O_b(\delta) have total length at most O_b(\delta)+2r_t. In every role, including the possible remaining role a,
Minkowski’s inequality transfers the assumed sign moments to these window sizes. Intermediate moment orders below eight are bounded by the eighth moment using Hölder’s inequality.
Write G_j=U_{j,\gamma}F_j. If role a has leaves above depth t, clip it radially at A=(K+2)^{10}. The clipped tail has local L^1 norm at most A^{-1}\|F_a\|_2^2, whose L^4 norm in the random parameters is O_b(A^{-1}). In the other roles the local L^3 norms have bounded fourth moments, by (5.2) and the preceding window bounds. The absolute form estimate with exponents (1,3,3,3) therefore gives expected tail error O_b(A^{-1}) per normalized cell. Even summing O((K+1)^2) positions and depths is harmless. For the rest of the proof, replace G_a by this clipped function when that role is present.
To return to the original continuous form, fix the random parameters and (u,v)\in\mathbb T_0^2, and plant bursts
Choose c'>0 so small that c'(2+\max_j|b_j|\sum_l|e_l|)<1. If a real configuration lies in bursts centered at p_j, set n=p_0, m=\sum_j e_jp_j. Then
The last quantity is an integer and hence is zero. For this slow pair the allowed offsets form the polygon \{|u'+b_jv'|<c'\ (0\le j\le3)\}, of fixed positive finite area. The center-valued kernel integral on the polygon therefore gives exactly the sampled product, times that area and \delta^2. The second kernel replacement has the same error estimate as the first: unitarity and clipping preserve fast L^2 bounds and do not enlarge the set of slow indices. We can now apply C_N(q) to these fixed compactly supported burst inputs, using the original admissible kernels.
For completeness, the exponent calculation is unchanged. Put
Use fast exponent p_* in the clipped role and exponent four in the other roles. Fast pairwise Hölder, random Hölder of exponent four, spatial Hölder with the p_j, and the mixed-norm maximal theorem [H, Lemma 7.1] apply: all exponents in that maximal theorem exceed q. The resulting regular-role factors are C_b|E_j|^{1/p_j}; the clipped factor is C_{q,b}A^{1-2/p_*}|E_a|^{1/p_*}. Summing O(K+1) positions in the two disjoint designated masks costs (K+2)^{1-2/p_*}. Hence the total power is
The logarithmic number of length classes finishes the proof.
This proposition supplies exactly the input used by the side estimates [H, Lemmas 7.3–7.4 and Proposition 7.5] and [V, Section 6.5]. Their remaining arguments concern nested projections, frozen stacks, and variation in depth. In addition to the local absolute form bound, they assume the square estimate H(17) for the linear sides. We identify its source below; no coordinate substitution is needed in those projection arguments.
5.2. A common gauge and a lattice of rational shifts
We describe the data from the mode construction of [H, Section 8] before checking its arithmetic. Fix a localized group G, with central tag z=(\theta_z,C_z,\lambda_z), and a role l. A primary is an available main column of this role’s outer projection. Its first matching cell is the cell where the finite expansion described below is chosen; that choice and its labels are then inherited on descendants. Thus later freezing does not resample a primary’s expansion.
Choose an integer Q_0 clearing the denominators of C_z, and put \vartheta=\theta_z/Q_0. For an offset \rho in the unit cube of horizontal coordinates, define
Choose \delta_c=\exp(-k^{C_{\rm cut}}) sufficiently small after the main budgets, with fixed exponent C_{\rm cut}. Let \zeta^\rho be the product Gevrey cutoff of [H, Section 8], vanishing within distance \delta_c of each branch cut and equal to one outside distance 2\delta_c. It is bounded between zero and one. For each role and group, one offset is kept fixed throughout its history. Its expectation is at least 1/2, and \mathbb E_\rho|1-(\zeta^\rho)^2|^2\le C d_{\max}\delta_c, where d_{\max} denotes an upper bound for the dimensions.
For a primary \Phi=(\varphi_\Phi,\epsilon_\Phi^{1/2}\mathbf e_\Phi) with normalized tag v, the endpoint splitting from the localization section and the chart expansion of [H, Lemma 4.4] give, at its first match,
The chart expansion uses this same error tolerance \delta_c, and all modes of one primary have the same quadratic coefficient. As in [H, Lemma 8.1], use mode atoms
The coefficient sum, mode count, normalized curvature, and rational heights have logarithms polynomial in the accuracy and localization parameters k,P. In particular, for modes of one primary,
Multiplication of the splitting by the fixed integer c_l changes only fixed constants in these bounds.
The resulting lift \mathcal R_{l,G} is one fixed pointwise bounded operator on augmented inputs. Indeed, set q_\Phi=\sum_m a_{\Phi,m}(\psi_m)_{\rm phys} -(\zeta^\rho)^2\varphi_\Phi, so that |q_\Phi|\le\delta_c. For an augmented input F, its lifted physical component is
and its mode coordinate belonging to (\Phi,m) is \epsilon_{\rm m}^{1/2}a_{\Phi,m} \epsilon_\Phi^{-1/2}F_\Phi. Choose \delta_c small and \epsilon_{\rm m} small after the main column budgets. Cauchy–Schwarz then gives a uniform operator norm. The private coordinate of an outer projection output records its column coefficient, so this operator gives exactly its mode synthesis. It is cell-local and independent of the freezing block.
There is also a common periodic multiplier for the four gauges. For four offsets, put d_j=d^{\rho_j}(y_j) and
The k_j are bounded integer vectors. Between branch cuts the lifts are affine and their affine parts cancel in this expression; hence the carries are locally constant. The three identities \sum c_jb_j^a=0, 0\le a\le2, give
For fixed carries the right side is an ordinary character on the torus with coordinates (\vartheta y_{l'},\vartheta t), for any anchor l': Q_0 clears C_z, and changing the anchor uses integer coefficients b_j-b_{l'}. Smooth cutoffs slightly larger than the original ones extend these character pieces by zero at the cuts. They give the smooth periodic multiplier, bounded by one, required in [H, Section 8.3].
The number of carries and the derivative bounds acquire only fixed factors per coordinate or derivative. Their logarithmic costs remain polynomial in P,k. Choose a mode budget D_*\ge e^k, \log D_*\le\operatorname{poly}_b(P,k), covering the mode lists, coefficient sums, heights, inverse private weights, and absolute Fourier sums, and large enough to dominate the dimensions and all fixed slope factors. Only then choose a Fourier cutoff B_* with tails at most D_*^{-300}. This preserves the order of choices in H; a coefficient-sum budget is not defined through its own truncation radius.
To control rational divisions at later frequency tests, fix an integer
The mode graphs will use shifts in the lattice generated by \vartheta/L_h, where
Here E_h bounds integer shift indices. The choice of D_{\rm s} is fixed independently of all depth and accuracy parameters. From this point, constants indexed by b may also depend on this fixed choice of D_{\rm s}.
We recall the scale geometry to make the graph conditions explicit. In a depth block [s_0,s_0+K), freeze coefficients at u_{s_0}=s_0-4K, compare on a cell J at depth v_{s_0}=s_0-2K, and put t_{s_0}=s_0+K. A mode has instantaneous frequency \omega_m(x)=2\lambda_mx+\gamma_m. On the list frozen at u_{s_0}, join two modes in the PRE graph if
The CURRENT graph has the same definition with threshold \Xi/r_{t_{s_0}}. Write \mathcal E_h(m), \mathcal F_h(m) for their components. Both graphs increase with height, because a witness n at height h becomes D_{\rm s}n at height h+1, provided A_0\ge D_{\rm s}. They also increase along successive depth blocks on inherited lists.
Here are parameter choices sufficient for the arithmetic. Let d_{\max}\ge2 bound all relevant dimensions, set \Delta=10, and put
Choose the largest height J_*, with 0\le h\le J_*, large compared with \delta_*^{-2}d_{\max}R'\log(2D_{\rm m}), where D_{\rm m} is the main-column budget and \delta_*>0 is a sufficiently small fixed number. Take A_0 to be a sufficiently large fixed power of C_*D_{\rm s}^{R'}B_*D_*, and \Xi=D_*^{1000}. Choose a rational height reserve and a frequency gap with
Finally choose K, polynomial in P,k, after these logarithmic budgets, so the height-bin packets of [H, Lemma 3.5] and the gap \mathcal T fit inside a block. Increase it so the pre-lag Taylor errors are as small as any prescribed fixed inverse power of D_*. All degrees can be fixed independently of q.
For a primary with horizontal parameters \Theta=(\theta_v,\vartheta), call a block nonexceptional when every rational vector a of height at most \mathcal H satisfies
The height-bin proof of [H, Lemma 8.2] gives O(d_{\max}) exceptional blocks per primary. At a nonexceptional block let K_{\rm rat} be the rational span of the tiny vectors, and set K_B=\{w:(0,w)\in K_{\rm rat}\}. Cramer’s rule and the choice of \mathcal T show that every height-\mathcal H vector in K_{\rm rat} is itself tiny.
Lemma 5.2 (Height stabilization). For a nonexceptional primary, all but O(d_{\max}R') heights h\in\{0,\ldots,J_*\} have the following property: two of its modes in one CURRENT component at height h are directly adjacent in PRE at height h-\Delta.
Proof. This is [H, Lemma 8.3] with the shift lattice (5.6); we give its divisibility and size argument. A simple CURRENT path has at most D_* edges. The equal quadratic coefficients of its endpoint modes and (5.3) yield
Indeed the path bounds its evaluation by D_*\Xi/r_{t_{s_0}}, and the gap forces the tiny alternative.
Let V_i be the rational span of the height-D_* vectors q for which q+(0,m/L_i)\in K_{\rm rat} with |m|_\infty\le D_*E_i. Let U_i be the span of the integer vectors in K_B of size at most C_*^3D_*E_i. These spaces increase. Discard heights below R' and those where either rank changes between h-R' and h; there are O(d_{\max}R') such heights.
At a remaining height, express a\in V_{h-R'} in a defining basis a=\sum_j\beta_ja_j, with witnesses m_j. Cramer’s rule bounds the coefficients and their common denominator by C_*. Subtracting the relations gives
After clearing denominators its size is at most C_*^3D_*E_h, so it lies in U_h=U_{h-R'}=:U. The quotient \mathbb Z^{\dim\theta_z}/(U\cap\mathbb Z^{\dim\theta_z}) is free, since the intersection with the rational space U is saturated. In an integer basis of this quotient,
For each prime p\mid D_{\rm s}, every integral coordinate on the right has valuation at least R'v_p(D_{\rm s})-\log_2C_*\ge\Delta v_p(D_{\rm s}). Thus [m_h] is divisible by D_{\rm s}^\Delta.
Choose an integral lift m' of this quotient class divided by D_{\rm s}^\Delta. Then
Round its coefficients in an integer spanning set for U whose vectors have size at most C_*^3D_*E_{h-R'}. The adjusted lift remains integral and satisfies
The last inequality follows from the stated fixed-power choice of A_0. No lattice basis or covolume estimate is involved. Since m_h-D_{\rm s}^\Delta m'\in K_B, the vector a+(0,m'/L_{h-\Delta}) belongs to K_{\rm rat} and remains within the height reserve. It is tiny by (5.7), which gives the required PRE edge.
We now specify the remaining output imported from [H, Section 8]. Write \pi(\mathcal D) for the primary labels represented in a set \mathcal D of modes. A mode is good at height h\ge\Delta when
and at least nine tenths of the primaries in \mathcal F_h(m) are nonexceptional and satisfy Lemma 5.2 at h. Heights below \Delta are bad. Call a mode core if it is bad at at most \delta_*(J_*+1) heights; call every other mode a linear-side mode. The increasing graphs and height stabilization give the coherence of these predicates, the bad-density count, and hereditary representatives by exactly [H, Lemmas 8.4–8.8, Proposition 8.9, and Lemma 8.10]. Those arguments use component inclusion, cardinality, projections, and the rational gap; their coordinate inputs have all been established above. A minimal PRE or CURRENT component means a component at height zero. Good modes in a CURRENT component have one PRE predecessor at height h-\Delta. Moreover, core modes in one minimal CURRENT component have identical good-height predicates throughout the upper range
Each core mode has bad fraction at most 2\delta_* in this range.
We also record the decomposition of the original localized outputs. Let Z_{l,G,s} be the current augmented output in role l. For a basic-B role in the block beginning at s_0, put
Let A_{l,G,s_0} and V_{l,G,s_0} be its exact core and linear-side mode subtotals, with all private coordinates retained. Define the full side by
In a basic-U role set S_{l,G,s}=Z_{l,G,s} and A_{l,G,s_0}=\mathcal E_{l,G,s_0}=0. On physical components these definitions give
For a basic-B role, the physical error is
The offsets can be fixed so that every term containing this physical error is negligible. The explicit lift, applied to Z_{l,G,u_{s_0}}=\sum_\Phi\beta_\Phi\Phi, gives
The main coefficients and original actual-use cells were fixed before the offsets. Their coefficient, size, and path-count bounds, together with the cutoff expectation, imply the error estimate from [H, Section 8, subsection on fixing offsets]:
With \delta_c chosen this small in the main budget, fix offsets satisfying this bound. All later selections only remove nonnegative terms; they need not be independent of the offsets. The linear-side square estimate is uniform in the offsets. The local absolute L^2 form bound and Cauchy–Schwarz therefore make the terms containing \mathcal E negligible, after the same size stops as in H. This explains why the subsequent estimates may use the exact core and full-side arguments.
Here are the quantitative conclusions used below. A successful CURRENT component contains a good core mode, whose PRE component at h-\Delta is its uniquely determined good predecessor. For fixed height, role, group, and the pair of floor logarithms of the predecessor’s primary and mode counts, choose a representative mode inherited under predecessor inclusion. On comparable actual use cells I,I', distinct representatives satisfy H(25), namely
for x_*\in I\cap I'. On an actual cell, a successful component \mathcal F has amplitude
Here \mathcal D runs over minimal CURRENT components, A_{\mathcal D} is the full core subtotal in the augmented mode space, and \alpha_m is its frozen mode coefficient. This is H(26); the added coefficient terms absorb uniform Taylor and Fourier errors.
The Bessel estimates and bad-density count give H(17):
This is the square hypothesis needed for the side estimates following Proposition 5.1. After the stopping construction of [H, Proposition 8.12], the representatives occurring in a dyadic amplitude bin \sigma_j along a retained path satisfy H(28):
Here \mathcal A_j^G is the set of distinct representative labels in that bin. All local side and component sizes are at most P^C.
For scheduled inputs, these statements are used in the form proved in [V, Section 6.3–6.4, especially Lemma 6.3]. The rows are cell-local projections composed with the fixed bounded lift \mathcal R_{l,G}. Their ordinary Bessel estimates come from the separated mode construction, and their adjoint main ranges lie in the usable endpoint spaces established by localization. V evaluates the depth input after the final main basic projection. Its row composition result therefore applies to the same size, variation, reading-depth, and endpoint hypotheses. In particular, the possibly large mode budget D_* is not inserted into the strong square constant G_0. This is the distinction required by V’s side estimate and by the final choice of q.
5.3. The side estimate and the near-band gain
We now estimate the contribution of the full sides to the cell minimum in (4.7). The input is the localized expression (4.13), its outgoing-filter identity, the endpoint variation estimate (4.14), and the linear-side square bound H(17). These are precisely the hypotheses of [H, Proposition 7.5]; for scheduled inputs their counterparts are assembled in [V, Section 6.5].
Expand each localized argument into its core, full side, and negligible physical error. Terms with at least three side factors have a length-weighted p-power bound P^C|R|, for any fixed 2/3<p<1. Indeed Cauchy–Schwarz first sums the group index in two side roles; writing the three resulting square sizes as a_s,b_s,c_s, their contribution is bounded using
This uses 3p>2 and the retained path caps; the fourth argument uses its P^C size bound. The group sum is formed before the power is taken.
In a term with exactly two sides, restore its two core complements to the whole current arguments Z_{c,G,s},Z_{d,G,s}. The corrections have a third side. The outgoing-filter identity then removes the complementary projections one at a time, transferring helper averages to one side, say role a. Its main term is
where \mathcal D_{a,s} is the resulting product of helper averages. Restoring the full base pair Y_c,Y_d introduces another third-side error. The pair-smallness from the base reduction bounds this main term by \eta times the two side sizes. Summing costs at most \eta P^C|R|.
At one transfer step, use the outgoing and helper bands defined in (4.15), with these two full sides as its fixed side operands. The exact outgoing-minus-helper defect is (4.16). Its far bands have negligible error. For each remaining offset h, |h|\le2, the retained expression (4.17) has coefficient -h/P=O(P^{-1}). Its paired sum over the band index must remain inside the cell contribution.
Freeze the complementary band stack and the other complementary argument at the block’s first depth. Their differences from current values have square budgets by (4.14). Freeze errors therefore supply a third strong factor and have the preceding p-power bound. The same argument handles cells where a frozen complement has large size: its difference from its bounded current value is then large. On the remaining ancestor hulls inside a block, both frozen complements have bounded local L^2 size.
The binary-prefix decomposition in [H, Lemma 7.4] now writes the two varying sides, at their uses, as selections from fixed vector stacks. The terminal complementary band 1-Q_{c,P} is first split into the identity and a positive ball term. Positive outer-ball ranges have joint group budgets, whereas group-independent identity and cumulative histories are counted once rather than once per group. This convention is needed for the following stack energies. One stack includes the band index. There are only G_0 overhead factors. Individual pad choices stay inside the actual positive averaged operators: a choice is coupled across depths, and Jensen’s inequality restores the average after estimating it. They are not additional stack configurations. The stack energies satisfy H(18):
Here \nu records the group, block, and fixed stack configuration. The functions remain fixed after subsequent cuts. Independent Rademacher signs select the two position indices, and a third family matches the band in \mathbf A_\nu with the frozen complementary band stack. Iterated Khintchine gives the moment hypothesis of Proposition 5.1.
For precision, the passage from these energies to a sum of local stack sizes uses the common hull cuts of [H, Section 7]. A hull consists of all ancestors between a use and its maximal use cell, including their full intervals. Write \Omega_\nu for the union of the original use cells. The main path and catalog counts give \sum_\nu|\Omega_\nu|\le M_0|R|, with \log M_0\le\operatorname{poly}_{q,b}(k). For a family of total stack energy at most E|R|, write \widetilde m_\nu for its full hull maxima. They obey
Stopping at level C_1E and integrating up to that level costs at most CE(1+\log(2+C_1M_0))|R|. Thus only a main-budget logarithm enters G_0. For the retained hull define
with value zero for an empty supremum; define m_{B,\nu} in the same way. The weak L^1 estimate for the dyadic maximal function, followed by truncation and integration, gives common individual cuts; Markov’s inequality for the truncated sum gives the joint cuts for which these squared hull maxima satisfy
The removed roots have arbitrarily small fixed total relative length. They are first crossings of partial ancestor quantities, so retained caps hold on whole intervals even in branches containing a later cut. A sparse decomposition of each retained hull into tree pieces gives disjoint major subsets of the piece roots. If a_Q,b_Q are the two normalized stack sizes at such a root, then
Complete each piece by its boundary children. The two side stacks then have the required leaf-size bounds, and the frozen complements have bounded leaf sizes. Proposition 5.1 applies to at most K depths of each piece. Multiplying by the transfer coefficient O(P^{-1}) yields the near-band bound
The logarithmic factor of the short-block estimate is included in G_0. Arbitrary signs at the original cells are allowed, so this is an absolute estimate for their coupled band contributions. No absolute sum over individual paired band indices has been used.
Thus [H, Proposition 7.5], with the coefficient-dependent comparison now proved, gives a complete conditional side bound. Third-side and freeze errors have p-power cost P^C|R|; transferred main terms have absolute cost \eta P^C|R|; the remaining near bands satisfy (5.10). Applying \min(\eta,x)\le\eta^{1-p}x^p only to the powered terms gives their contribution C_{q,b}\eta^\varepsilon P^C|R|, with \varepsilon=1-p>0. All cuts remove original summands at first crossings and the data are reconstructed on each normalized subroot. This is the root-uniform side estimate used in the closing argument. The scheduled version follows from the same fixed stacks and the row estimates already stated, as in [V, Section 6.5].
5.4. Frequency tests for three or four core arguments
Fix a retained path point x_*, an actual cell I=I_s with r=|I|, a group, and a common successful height h\in\mathcal H_{\rm up} for the core arguments. Put R_s=r^{-1} and
The representative a_j belongs to the good predecessor of the component in role j. Summing a PRE path and using the pre-lag Taylor bound expresses its selected physical argument, up to the controlled uniform error, as
The function g_j is a cutoff times a trigonometric polynomial. After collecting the common multiplier (5.4), all Fourier indices in any anchor coordinates lie in the box of radius B'_h, enlarged by fixed factors depending on b. This is H(29) with its complete coefficient-sum and error bounds.
To compare physical and torus amplitudes, take the horizontal tiny space K_B of a nonexceptional primary in a successful component, and let
This is a connected torus because the integer lattice being annihilated is saturated. The tiny/huge alternative makes physical interval averages of the relevant characters agree with Haar averages on \mathbb H, to the prescribed small error. The component floors above absorb those errors, giving H(30):
The selected sets are whole minimal CURRENT subtotals, by coherence of the height predicates on \mathcal H_{\rm up}. For distinct roles, u+(b_j-b_{l'})v and u+(b_i-b_{l'})v have independent Haar distributions on \mathbb H^2. Interpolating the six two-L^1 bounds gives
A bounded progression multiplier changes only the constant.
For four core roles put
Choose, by a fixed rule independent of the mode masks, a reference q'=Q+n\cdot\xi, t'=T+m\cdot\xi in the full Fourier box with \max(|q'|,|t'|)\le U_0R_s, where U_0=D_*^{250}. Introduce dyadic parameters by
Discard absent references and inner minima below U_0^{-1}. The four marginals of the normalized kernel make its Fourier transform vanish on t'=b_jq'. The mean value theorem near these lines and rapid Fourier decay give C_{L,b}\mu\nu^{-L} for every fixed L. The tiny/huge dichotomy puts all moderate frequencies in one pair of character cosets modulo K_B; summing that pair of cosets is a Haar integral bounded by (5.11). Thus its contribution is at most C_{L,b}\mu\nu^{-L}\prod_j a'_j.
For three core roles with side role l, anchor at y_l and put
Choose a moderate reference T_l+m\cdot\xi and write \beta=r(T_l+m\cdot\xi). Set \mu\asymp\min(1,|\beta|), \nu\asymp\max(1,|\beta|), with the same small and large discards. Partial Fourier integration of the kernel in the step variable gives a function W_\beta whose fixed-order derivatives are bounded by C_{L,b}\mu\nu^{-L}. At zero this uses the l-marginal; at large frequency it uses integration by parts. Duality in (5.11) bounds the remaining horizontal coefficient sequence in \ell^2.
More explicitly, average a fresh offset cutoff \zeta^{\rho_l}/c_\zeta in the side, where c_\zeta=\mathbb E_{\rho_l}\zeta^{\rho_l}\ge1/2 is constant. This insertion is an exact identity in expectation. Keep the bounded factor c_\zeta^{-1} outside the estimate and set S'=\zeta^{\rho_l}e(-c_lp^{\rho_l})S_{\rm phys} on I. Its normalized L^2 norm is at most that of the side. Write \widehat S'_I(k') for its normalized Fourier coefficients. The horizontal pairings decay in |k'+r(Q+n\cdot\xi)|. Schur’s estimate on a dyadic shell of radius \Lambda, followed by averaging the fresh offset, gives H(32) with
The contribution is bounded by C_{L,b}\mu\nu^{-L}\prod_{j\ne l}a'_j \sum_{1\le\Lambda\le U_0}\Lambda^{-L}b_\Lambda, where the sum is dyadic. Notice that the windows use the whole Fourier box and are independent of the mode masks.
The error estimates in [H, Section 9.1–9.2] now apply with these explicit tests. They use the absolute local form estimate, Fourier tails, the pre-lag, and the rational gap, all with reserves already chosen above. Per tuple the error is bounded by a fixed multiple of D_*^{-100} times the absolute core coefficient sums, and by the side size if a side is present. The common-height argument samples heights independently and uniformly from \mathcal H_{\rm up} until all roles are good; its first-success weights sum to one on each pure tuple. A core mode has bad fraction at most 2\delta_*, so the sum over failure assignments converges. Consequently these errors sum absolutely to O(e^{-k}|R|), and the principal cell terms have the form H(33):
Here the \sigma_j are dyadic component-amplitude bins, with the side bin belonging to b_\Lambda; N_s^G counts the representative tuples passing the tests. For four cores, \Lambda=1. We next bound these counts across scales.
5.5. Sparse mass for separated tests
The finite counting statement is useful independently of the mode construction. Its integer divisibility is the reason for the choice (5.5).
Lemma 5.3 (Sparse mass). Let J\subseteq\{0,1,2,3\} consist of three or four core roles. In each role let a finite label set \mathcal A_i be identified injectively with real numbers, with |\mathcal A_i|\le M, M\ge1. Let d\ge0 be an integer. Fix \zeta\in\mathbb R^d, B,L\ge1, Q,\Sigma,R_0>0, and 0<\mu\le1\le\nu. Every shift vector below belongs to \mathbb Z^d. For finitely many distinct integers s, at dyadic scales R_s=R_02^s, let \mathcal E_s\subseteq\prod_{i\in J}\mathcal A_i and N_s=|\mathcal E_s|. Assume the following five conditions.
- 1.
Two labels in different core roles determine at most one full tuple, and that pair occurs at at most L scales.
- 2.
If there is a side role, there are w_s\ge0, \sum_s w_s\le M, such that every fixed core label is incident to at most w_s tuples at scale s.
- 3.
Distinct same-role labels occurring at scales s,t satisfy
|a-a'+n\cdot\zeta|>\Sigma\max(R_s,R_t) \quad(|n|_\infty\le Q). - 4.
For every triple T\subseteq J, with excluded role l'\in\{0,1,2,3\}\setminus T and \alpha_i=b_i-b_{l'}, each tuple has a witness |n|_\infty\le B such that
\mu R_s\le\left|\sum_{i\in T}\alpha_i a_i +D_{\rm s}^2 n\cdot\zeta\right|\le\nu R_s. - 5.
At every occupied scale, |n\cdot\zeta|\notin[\mu R_s/1000,1000\nu R_s] for all |n|_\infty\le Q.
There are constants C_b,C_{{\rm res},b} such that, if
then, with \ell=1+\log_2(2+(d+1)^2\nu/\mu),
The exponents and constants are independent of the dimension and the number of labels and scales, apart from the indicated factors.
Proof. We follow the finite argument of [H, Lemma 9.1], spelling out the arithmetic and the dimension dependence. The proof first covers most label differences by neighborhoods associated with later scales. Three such future covers then place sums of triples in a thin slab. Separated anchor labels reduce that slab’s cover from three factors of M to two, while the pure-shift gap makes many translated triples distinct inside the cover. This comparison forces small mass at sparsely occupied scales.
Approximate all finite real data on a sufficiently fine common grid, taking approximants on a fixed divisible multiple of its step so divisions by the \alpha_i remain on the grid. For example, a multiple of D_{\rm s}^2 suffices. Choosing the grid after the data preserves separation and annular inequalities with fixed relaxed margins and leaves a narrower pure-shift gap. Round scalar widths upward by at most one grid step. The estimates below are uniform in that step.
Restrict to N_s\le\delta M^2. Pair determination gives total mass at most LM^2. For parameters \delta\le\alpha/2 and 0<\kappa<1, greedily partition the scales in decreasing order of R_s into counting blocks of mass at most \alpha M^2. There are at most J_0=3L/\alpha blocks. For a block \mathfrak b, let R_{\mathfrak b}^+ be its largest scale, take N_0=\lceil C_bB\rceil with C_b>32D_{\rm s}, and choose a grid width \epsilon_{\mathfrak b}\asymp\nu R_{\mathfrak b}^+. Use the finite neighborhoods
They are symmetric and satisfy U(t)+U(t')\subseteq U(t+t').
Fix a role and a counting block with uncovered mass exceeding \kappa M^2. Project onto this role and another core role, retaining one scale witness per pair, and color edges by a third core label. The graph has at least (\kappa/L)M^2 edges and a proper coloring with at most M colors. The elementary four-step-walk observation in the proof of [H, Lemma 9.1] gives a subset B_0 of at least c\eta M vertices, \eta=\kappa/L, such that any two have at least c\eta^5M^3 four-step walks with distinct ordered color tuples. Repeated vertices are allowed. That observation is purely a graph statement: choose a random neighborhood, remove vertices with many low-codegree partners, and count the two choices of a common neighbor on either side of a retained midpoint.
Alternating the four annular relations along a walk cancels the intermediate labels. Dividing by the nonzero integer \alpha_i gives
Here q depends only on the ordered four edge colors; these c_i are graph colors, unrelated to the fixed quadratic weights. The scalar error after division is at most 4\nu R_{\mathfrak b}^+, and D_{\rm s}^2/\alpha_i is an integer multiple of D_{\rm s}, with its index covered by N_0. Packed differences outside U_{\mathfrak b}(2) cannot use the same color tuple. The M^4 possible color tuples therefore bound their number by C\eta^{-5}M. Repeating after removing B_0 covers all but \kappa M^2 occurrences in each role and block by at most K_1=C(L/\kappa)^6 subsets, each with its difference set covered by K_1M translates of U_{\mathfrak b}(2).
Discard uncovered occurrences. Call a core label future when it belongs to a cover subset in a later counting block. Nonfuture label sets are disjoint over blocks. Let m_i^{\mathfrak b} be their cardinalities in core role i. Treat the side as nonfuture, with m_l^{\mathfrak b}=\sum_{s\in\mathfrak b}w_s. Occurrences with fewer than three future roles have two nonfuture roles. Their total mass is at most C\sqrt{L\alpha}M^2, by summing \min(\alpha M^2,Lm_i^{\mathfrak b}m_j^{\mathfrak b}) and Cauchy–Schwarz. The discarded uncovered mass is at most 4J_0\kappa M^2.
Assign each remaining occurrence three future roles and one future cover subset in each. A class fixes its current counting block and its three future block/subset choices. There are at most 4J_0^4K_1^3 classes. Fix one of mass \rho M^2, with future subsets B_1,B_2,B_3. The subscripts now enumerate those roles; their coefficients remain the corresponding \alpha_i. Set \Gamma=C_{{\rm gap},b}(d+1)^2\nu/\mu. If \rho\le C_b\delta\ell, retain this bound. Otherwise split the mass into four consecutive ranges by quartiles, omitting the at most three scale bins crossing a cut. Each range retains at least \rho M^2/8 mass. The first provides anchor sets B_{i0}\subset B_i of size at least c\rho M/L. Let R^\# be the largest scale of the second range; these anchors are separated modulo reserved shifts at scale R^\#. The fourth range gives a set F of at least c\rho M^2/L distinct triples projected onto these three roles. The third range supplies the scale margin, and a greedy thinning of the second gives q\ge c_b\rho/(\delta\ell) triples t_v, also projected onto these roles, at scales separated by factors at least \Gamma. All future-block scales and fourth-range scales are smaller than these by at least the required factor \Gamma.
Write U_i for the neighborhoods of the three future subsets. Subdividing index and scalar intervals covers U_i(t) by at most (2t+1)^{d+1} translates of U_i(1). Hence there is an odd H_i\le100(d+1)^2 with |U_i(H_i+2)|\le2|U_i(H_i)|. The preliminary exponential cover is used only to find this small-growth index, and is never paid in the cardinality estimate. The translates a+U_i(H_i), a\in B_{i0}, are disjoint: an intersection would contradict the anchor separation with a shift index O_b((d+1)^2B) and scalar error O_b((d+1)^2\nu R^\#).
Set Z_i=B_{i0}+U_i(H_i). The difference covers give |Z_i-B_i|\le C(K_1L/\rho)|Z_i|. Choose a nonempty X_i\subset Z_i minimizing r_i=|X_i-B_i|/|X_i|. Petridis’s minimal-growth subset argument, in the form reproduced in [H, Lemma 9.1], gives |X_i-mB_i|\le r_i^m|X_i|. One anchor translate contains a subset Y_i\subseteq X_i of at least |X_i|/M points. If \mathcal P\subseteq mB_i has all distinct differences outside U_i(2H_i), the sets Y_i-p, p\in\mathcal P, are disjoint and lie in X_i-mB_i. Thus |\mathcal P|\le r_i^mM. A maximal packing covers mB_i by that many translates of U_i(2H_i). Only m=2,3 is used.
The outputs t_v+F lie in the slab
Put r_0=CK_1L/\rho. Cover the first two sumsets by r_0^2M boxes each. Above a fixed pair, dividing the slab equation by \alpha_3 places the third coordinate in
Every shift remains integral: both D_{\rm s}^2/\alpha_3 and D_{\rm s}\alpha_i/\alpha_3 are integers. If \mathcal P is a packing of the local third-coordinate set outside U_3(2H_3), then B_{30}+\mathcal P\subset3B_3 is still a packing. Collisions with different anchors are excluded by the same shift and scalar bounds. Thus |B_{30}|\,|\mathcal P|\le r_0^3M. The entire slab is covered by at most CK_1^7(L/\rho)^8M^2 product boxes.
Each such box contains at most one of the indexed outputs t_v+F. Within one translation this follows from label separation. Across translations, let R_{\rm big} be the larger selected scale. Its annular residual has magnitude at least \mu R_{\rm big}; the other three residuals together are at most \mu R_{\rm big}/10. Membership in one box adds scalar error at most C_b(d+1)^2\nu\max_i R_{{\rm future},i}^+ \le\mu R_{\rm big}/10, by the choice of C_{{\rm gap},b}. The resulting pure shift has magnitude between 4\mu R_{\rm big}/5 and 3\nu R_{\rm big}, with index O_b((d+1)^2B)\le Q, contradicting the forbidden gap.
Comparing q|F| with the slab cover yields \rho^{10}\le C_b\delta\ell K_1^7L^9, and hence \rho\le C_b\ell K_1L\delta^{1/10}. Including all choices and losses gives
Take \alpha=\delta^{1/1040}, \kappa=\alpha^2 for sufficiently small \delta, and use the crude mass bound for the remaining values. This proves the stated exponent 1/2080.
5.6. Application of the count and the resulting core bound
We verify the hypotheses rather than apply the original arithmetic form of [H, Proposition 9.2]. Fix a path, height, role types, frequency bands, and amplitude bins in (5.14). Identify each representative label with its value \omega_a(x_*), which is injective in a fixed role and type by (5.8). Define
Groups with a nonzero count have M'\ge1.
Two fixed core labels determine the others. For four cores, subtract the equations for Q,T of two candidate tuples. The two unknown label differences solve a system with determinant b_j-b_i, modulo bounded rational shifts and errors O_b(\nu\max(R_s,R_t)). For three cores the equation for T_l determines the last label, dividing by b_j-b_l. The shift reserve and separation force these differences to vanish. For every core triple with excluded role l', the tests also give the annulus
For four cores this is t'-b_{l'}q'; for three it is the T_l test. Once the labels are fixed, the height-bin lemma [H, Lemma 3.5] bounds their scale occurrences by L_0\le P^C(1+\log(\nu/\mu))^C.
If one role is a side, windows (5.13) for different triples incident to a fixed core label are disjoint. An overlap gives both a relation for their Q values with error O(\Lambda R_s) and their two T_l relations. Solving the resulting two-by-two system again divides by a b_j-b_i; separation then identifies the tuples. Couple the fresh offset across these windows. Parseval shows their incidence count is at most w_s=\sigma_l^{-2}\|S_{l,G,s}\|_{2,I_s}^2, giving the side-budget hypothesis.
Here are the promised reserve checks for all these divisions. In Lemma 5.3 take
Every divided \xi-shift is integral in this finer lattice, because each determinant divides D_{\rm s}. The fixed-power choice of A_0 ensures
Also \nu,\Lambda\le C_bD_*^{250}, while \Xi=D_*^{1000}, so the separation reserve holds. All rational shifts remain below \mathcal H. Every successful component contains a nonexceptional primary, even if its representative is exceptional. Applying (5.7) to its pure horizontal shifts excludes the entire intermediate interval required by the pure-shift hypothesis, since \mu\ge U_0^{-1} and the interval length |I_s| lies between r_{t_{s_0}} and r_{v_{s_0}}. This checks all five hypotheses.
Lemma 5.3, pair determination, and the side incidence bound therefore give the two counts H(34): for some fixed p_0<1,
Indeed the sparse-mass exponent c'=1/2080 allows any 1-p_0<c'. On the bins 2^{-j-1}(M')^2<N_s^G\le2^{-j}(M')^2, apply the sparse bound and sum the geometric series with ratio 2^{1-p_0-c'}<1.
Finally the amplitude and group summation of [H, Proposition 9.3] uses precisely (5.14), (5.9), and (5.16)–(5.17). Interpolate the two counts with a sufficiently small fixed weight \theta'>0, giving power p=(1-\theta')+\theta'p_0. The path square budgets sum the group factors, and the core count cap D_* controls the common amplitude shift. The three amplitude-gap sums are geometric, as in that proposition. The frequency losses after interpolation are at most \mu^{-C\theta'}\nu^{C\theta'}\Lambda^{C\theta'}; choose \theta' small and then the derivative order L large so they are summable against \mu^p\nu^{-Lp}\Lambda^{-Lp}. Applying the original cell minimum before subadditivity yields
for the terms with at least three core arguments, with fixed \varepsilon>0 and C_4 independent of q. The same calculation is the core step of [V, Section 6.6]: its frozen scheduled inputs provide exactly the amplitudes, representatives, and same-cell side incidence verified above. Together with Proposition 5.1 and the side square bound, this supplies the side and core estimates required to close the proof of Theorem 2.2.
6. Parameter choice and completion of the local estimate
The preceding sections provide the coefficient-dependent inputs to the completion arguments in [H, Section 10] and [V, Section 7]. We now explain the common parameter choice and the two distinct conclusions obtained from it. First we bound the fixed-input constant C_N(q) independently of N. Only after this absorption do we complete the proof for depth-dependent inputs.
6.1. The retained-cell estimate
We use the base reductions of [H, Proposition 5.8] and [V, Proposition 5.6], with the modifications proved above. At one noninitial accuracy level their target, on a normalized root R, is a sum of terms \abs I\min\{\eta,\abs{H_I(z_I)}\} over a specified subset of descendant cells. Here z_I is the four-tuple supplied by the base expansion at I, and \eta=\exp(-k^\theta) is its quality parameter, with k the accuracy log and 0<\theta<1 fixed. The constructions then decompose each of these four-tuples into the localized terms estimated in the preceding sections.
The stopping constructions in [H, Section 10.1] and [V, Section 7.1] are used in their stated order. Main columns and their uses are fixed first; mode expansions and the whole side functions are fixed next; helper projections are constructed after those side functions. Size and counting stops are then imposed on these fixed objects. A stop removes original summands from an attempt; it does not redefine a previously fixed stack function. In particular, the maximal functions in the short-block estimates remain defined on whole retained cells, including the portions below later stopping cells.
Let \mathcal A_R^{\mathrm{ret}} denote the cells not removed by these stops in the attempt rooted at R. The conclusions of the preceding sections give the same estimate as the retained-cell calculations in the two companions:
Here \varepsilon>0 is fixed, P is the transfer parameter, K is the depth-block length, and G_0 is the row and stack overhead. The constants may depend on b,\delta_0,C_0.
For clarity, the estimate has three sources. Terms with at least three core factors use the powered core estimate, and terms with at least three side factors use their square-energy bounds. With p=1-\varepsilon>2/3, both are converted to the first term of (6.1) by \min\{\eta,a\}\le\eta^{1-p}a^p. Terms with exactly two sides use the transfer estimate for the whole complementary arguments. Its main term is again included in the first term. The remaining coupled near-band sum is estimated at power one by the short-block comparison; it gives the second term, including its factor 1+C_N(q), by (5.10). Physical and numerical errors are estimated at power one as well. These alternatives exhaust the four-position expansion. The minimum is applied to the whole contribution of a cell before subadditivity is used; no smallness of a separately selected group is being assumed.
6.2. Restarts and the complexity hierarchy
The integrated size, path-count, and square-energy bounds allow the stopping thresholds to be enlarged by fixed factors so that the first stopped cells occupy at most \abs R/4. At each such cell send the original summands there and below to a new attempt. The local normalizations and column bounds are available on every new root. Linear-side energy is proved anew by the construction on that root, not by asserting a proportional restriction of a global energy bound. The resulting attempt roots satisfy
All the trees are finite, so the restart argument can equivalently be read as induction on their remaining depths.
In the scheduled-input construction, the original depths, functions, and block partitions remain attached to the summands. Restriction to a new root preserves the pointwise size and variation bounds with the same W_j. The row estimates of [V, Section 4] therefore apply at each restart with no factor depending on W_j.
The distinction between the two complexity scales is essential. The bounds supplied by the constructions are
Here D_{\mathrm m} bounds the main and helper catalogs, and D_* bounds the larger mode, Fourier, and rational-height data. The latter data enter the counting and mode estimates, but are not fed back into the main catalog construction. Strong stack estimates compress main projection families only; mode projections have already been included in the Bessel rows. This is precisely the separation in [H, Section 10.3] and [V, Section 7.3].
The coefficient modifications preserve this separation. Chart substitutions multiply heights and derivative bounds by fixed slope-dependent factors; the change of height base occurs only in the mode construction. The sparse-count estimate retains its dimension-polynomial reserves and its power saving. Thus all the degrees C,C_4,C_5,C_{19} can be fixed before choosing q. Their numerical values may be enlarged to cover both the fixed-input and scheduled-input constructions. Only coefficients such as C_q are allowed to depend on the later choice of q.
Choose q\in(2,3) sufficiently close to 2 that
for both constructions. This common choice is legitimate because local detection is available for every q>2 and short-block comparison for every 2<q<3. Next put P=\lceil\eta^{-\alpha}\rceil, choosing a fixed \alpha>0 so small that \alpha C_4<\varepsilon/2, as well as the analogous finitely many upper bounds already included in the powered costs. Finally choose the starting accuracy sufficiently large. Since P dominates every fixed power of k and G_0=P^{o(1)}, we obtain
Also P^{-1/4}\le\eta^{\alpha/4}, and e^{-k} is smaller than these quality powers for large k. Consequently, for some c_*>0, (6.1) and (6.2) give the full base-segment bound
6.3. Fixed inputs: accuracy summation and absorption
The base accuracy levels satisfy k_{t+1}=k_t^A for a fixed A>1. For tuples whose later level is t, there are O(t) choices of the earlier level and a bounded number of slot patterns. The segment-root packing in [H, Section 5], followed by (6.2), gives a bounded total root length for each tuple. The noninitial terms in (6.5) therefore have total coefficient controlled by
This series is finite, and its value is arbitrarily small when the starting accuracy is sufficiently large. The replacement and two-increment errors have the same summable multiplicity. The finitely many initial levels have depth-independent bounds from the active-scale reduction, as in [H, Section 10.4]. Thus the modified construction proves the following normalized conclusion.
Proposition 6.1 (Normalized absorption). For the fixed q chosen above and every sufficiently small \tau>0, there is A_\tau<\infty, independent of N, such that the following holds. Suppose a finite rooted dyadic subtree has root I_*, its summation cells lie in at most N consecutive depths, and its fixed scalar inputs satisfy
Then, for any subset of summation cells and admissible kernels,
The argument just given is the proof: take the coefficient in (6.6) below \tau and include its constant part, the initial levels, and the numerical errors in A_\tau. The base expansion is first taken at a finite accuracy cutoff. At each fixed cell the terminal residual tends to zero by the pair-comparison estimates of the base reduction; since there are only finitely many cells, the passage to all accuracy levels is legitimate. This is the concluding expansion step in [H, Section 10.4].
We include the outer normalization argument to make the absorption explicit. Start with disjoint roots at the largest scale in a finite collection from (2.7). At a root J, write a_{j,J}=(\abs J^{-1}\int_J\abs{f_j}^q)^{1/q}. If some a_{j,J} vanishes, all forms below J vanish. Otherwise stop at the first strict descendants J' on which an average of \abs{f_j}^q exceeds 16a_{j,J}^q. The union of first stops has measure at most \abs J/4. Dividing f_j by 16^{1/q}a_{j,J} normalizes the intervening tree, so Proposition 6.1 bounds its sum by
Repeat at the stopped roots. The sets E_J left after deleting their first stops are disjoint, have measure at least 3\abs J/4, and satisfy M_qf_j\ge a_{j,J} on E_J. Hence
Taking the supremum in the definition of C_N(q) gives
where C is independent of N and bounded for 2<q<3. Choose C\tau<1/2 and use the already established finiteness of C_N(q) to absorb the last term. We have proved
6.4. Depth-dependent inputs
Now perform the scheduled-input construction with the same q, using its own subsequent choices of \alpha and starting accuracy. Insert (6.7) into (6.5). The factor 1+C_N(q) is now a fixed constant. The scheduled base expansion of [V, Section 5] has the same accuracy multiplicity and segment-root packing, with the depth-row estimates supplying the uniform bounds. Its noninitial terms therefore sum by (6.6). Its finitely many initial levels are bounded by its active-scale reduction independently of the number of depths and of every W_j.
Finally, at a fixed original cell, all accuracy updates are evaluated on its same actual scheduled input z_{j,s(I)}. The per-cell pair comparisons therefore make the terminal residual tend to zero. There are finitely many cells, so the finite-cell sum commutes with this limit. The triangle inequality and the uniform summed bounds give (2.4), exactly as in [V, Section 7.4]. This proves Theorem 2.2, and hence, by Section 3, Theorem 1.1.
References
- [Ass]I. Assani, Multiple recurrence and almost sure convergence for weakly mixing dynamical systems, Israel J. Math. 103 (1998), 111–124.
- [Ber]V. Bergelson, Weakly mixing PET, Ergodic Theory Dynam. Systems 7 (1987), no. 3, 337–349.
- [Bir]G. D. Birkhoff, Proof of the ergodic theorem, Proc. Natl. Acad. Sci. USA 17 (1931), no. 12, 656–660.
- [Bou]J. Bourgain, Double recurrence and almost sure convergence, J. Reine Angew. Math. 404 (1990), 140–161.
- [Cal]A. P. Calderón, Ergodic theory and translation-invariant operators, Proc. Natl. Acad. Sci. USA 59 (1968), no. 2, 349–353.
- [DL]J.-M. Derrien and E. Lesigne, Un théorème ergodique polynômial ponctuel pour les endomorphismes exacts et les K-systèmes, Ann. Inst. H. Poincaré Probab. Statist. 32 (1996), no. 6, 765–778.
- [Fur]H. Furstenberg, Ergodic behavior of diagonal measures and a theorem of Szemerédi on arithmetic progressions, J. Analyse Math. 31 (1977), 204–256.
- [GHSY]Y. Gutman, W. Huang, S. Shao, and X. Ye, Almost sure convergence of the multiple ergodic average for certain weakly mixing systems, Acta Math. Sin. (Engl. Ser.) 34 (2018), no. 1, 79–90.
- [HK]B. Host and B. Kra, Nonconventional ergodic averages and nilmanifolds, Ann. of Math. (2) 161 (2005), no. 1, 397–488.
- [JKRW]R. L. Jones, R. Kaufman, J. M. Rosenblatt, and M. Wierdl, Oscillation in ergodic theory, Ergodic Theory Dynam. Systems 18 (1998), no. 4, 889–935.
- [LSS]J. Leng, A. Sah, and M. Sawhney, Quasipolynomial bounds on the inverse theorem for the Gowers U^{s+1}[N]-norm, preprint, arXiv:2402.17994v3, August 10, 2024.
- [H]OpenAI, An L^3 bound for the trilinear Hilbert transform, OpenAI Math Release preprint PDF, 2026. Reading view on this site: An L3 bound for the trilinear Hilbert transform.
- [V]OpenAI, Pointwise convergence of triple ergodic averages for mixing transformations, OpenAI Math Release preprint PDF, 2026. Reading view on this site: Pointwise convergence of triple ergodic averages for mixing transformations.
- [Tao]T. Tao, Cancellation for the multilinear Hilbert transform, Collect. Math. 67 (2016), 191–206.
- [Zie]T. Ziegler, Universal characteristic factors and Furstenberg averages, J. Amer. Math. Soc. 20 (2007), no. 1, 53–97.