Plectis Atlas
Search

↑ and ↓ move, Enter opens, Esc closes.

Reading view

Almost-everywhere Fourier convergence in L log L

OpenAI, 23 September 2026

PDFTeX source

Part of family 075: The L log L Fourier-convergence conjecture. Read it with its family

This page is converted from the paper’s TeX source. Atlas reconstructs statement numbers from the paper’s typesetting source (TeX). Use the original PDF if they differ. 2 figures are shown as placeholders; see the PDF for them.

Abstract

We prove that the ordinary symmetric Fourier partial sums of every complex-valued function in L\log L on the circle converge almost everywhere to the function along the full sequence. This establishes the classical L\log L sufficiency conjecture.

1 Introduction

Let \mathbb T=\mathbb R/(2\pi\mathbb Z), with normalized Lebesgue measure d\mu(x)=dx/(2\pi). For an integrable complex-valued function f, write

\widehat f(k)=\int_{\mathbb T} f(x)e^{-ikx}\,d\mu(x), \qquad S_N f(x)=\sum_{k=-N}^{N}\widehat f(k)e^{ikx}.

The endpoint question concerns the Orlicz class

L\log L(\mathbb T) =\left\{f:\int_{\mathbb T}|f(x)|\log(2+|f(x)|)\,d\mu(x)<\infty\right\}.

Theorem 1.1 (Almost-everywhere convergence in L\log L). For every complex-valued f\in L\log L(\mathbb T), there is a measurable set N_f of measure zero such that

\lim_{N\to\infty}S_Nf(x)=f(x) \qquad (x\in\mathbb T\setminus N_f).

The limit is along the full sequence N=0,1,2,\ldots.

Carleson proved almost-everywhere convergence for square-integrable functions [3], and Hunt extended this conclusion to L^p for every p>1 [11]. Kolmogorov’s integrable counterexample shows that the conclusion cannot hold throughout L^1 [13]. Between these cases, Sjölin obtained convergence in L\log L\log\log L [19], and Soria developed larger convergence spaces through extrapolation [20, 21]. Antonov proved convergence in L\log L\log\log\log L [1]. Arias de Reyna subsequently introduced the quasi-Banach space \mathrm{QA}, which strictly contains both Antonov’s class and Soria’s space [2]. Carro, Mastyło and Rodríguez-Piazza further studied \mathrm{QA} and obtained a Lorentz convergence space strictly larger than Antonov’s space [4]. Here the iterated-log class names describe growth at infinity; the logarithms may be regularized near zero.

Lie gave a time-frequency proof of the full-sequence maximal bound from \mathrm{QA} to weak-L^1, avoiding extrapolation [16, Corollary 1.2(5)–(7)]. Di Plinio and Fragkos proved weak-L^p bounds of order (p-1)^{-1} as p\downarrow1, together with sparse estimates and weighted endpoint convergence results [7]. The L\log L conjecture is also formulated through a weak-L^1 maximal estimate in [17, Definition 1.1 and Conjecture 1]. Theorem 1.1 establishes the conjectured sufficiency of L\log L for the ordinary Fourier sums along the full sequence. Stein’s Banach-space theorem also yields the corresponding Luxemburg-norm weak-L^1 maximal bound; see Corollary 10.3. Our proof proceeds through the good-set integral estimate described next.

1.1 The estimate behind the convergence theorem

The main analytic estimate concerns one nonnegative input \rho of mass one and height at most e^d, where d is large. Its uncentered spatial maximal density is

M\rho(x)=\sup_{\substack{I\ni x\\I\text{ a circle arc}}} \frac{1}{\mu(I)}\int_I\rho\,d\mu, \qquad 0<\mu(I)\le1.

Set h=\log\log d and remove the points where M\rho>e^{4h}. The removed set has measure at most Ce^{-4h}, while on its complement we prove

\int_{\{M\rho\le e^{4h}\}}\sup_{N\ge0}|S_N\rho|\,d\mu \le Cd. \tag{1.1}

The two bounds have different roles. The integral bound is summable with the masses of the height layers of an L\log L function. The exceptional sets are summable when the layers have doubly exponential heights. Section 10 makes both statements precise and proves Theorem 1.1 from Equation (1.1).

Good-set integral estimates also occur in Lie’s near-endpoint argument [16, Proposition 3.1 and (3.13)]. Here the bounded-density estimate and exceptional-set bound are quantified for the summation at the L\log L endpoint.

To obtain Equation (1.1), we first use a finite binary tree of dyadic intervals. In unit-period coordinates, modulation multiplies the input by e^{2\pi iuy}, where u is an integer. At each interval, take half the difference between the averages of this modulated input on its two children. Along a uniformly chosen output path, this coefficient is known before the next left-or-right choice. Its output uses the sign of one further binary choice; a multiplier of modulus at most one may depend on the intervening choice. Stop the path when the average of the original input first exceeds e^{4h}. Proposition 6.1 bounds the integral of the supremum over all integer modulations of the resulting sum by Cd, uniformly in the tree depth. This is the finite estimate assembled by the combinatorial part of the proof.

Translation and scale averaging of these delayed Haar kernels produce a nonzero multiple of the Hilbert kernel, with an error controlled on the same good set. The periodic Hilbert transform then gives the ordinary symmetric Fourier projections. The averaging follows the dyadic representation method of Petermichl and Hytönen [18, 12]; Section 10 includes the kernel calculation and all limiting steps.

1.2 The structural argument

Fefferman’s proof organized the Carleson operator by a measurable frequency selector [8, 9]. Lacey and Thiele later gave a simplified proof of its weak-(2,2) bound using phase-plane tiles and trees [14].

The finite proof must control many modulation frequencies without a loss that grows with their number. It splits the input into pieces and, on each piece, uses a noisy observation of the input point to identify frequencies with nearly the same phase. The input point is sampled with probability proportional to the density, whereas the transform is estimated along a uniformly sampled output path. Information is counted under the first law and converted to path estimates under the second.

Three constructions make these estimates compatible.

  1. 1.

    One phase region for a large set of frequencies. The compression lemma starts with a finite set of integer frequencies and a uniform bound on noisy sums of independent samples from any probability law on that set. It produces a large subset whose phase defects, relative to one center frequency, are small on average over one common region. The conclusion remains valid for every probability law on the retained subset, so later spectral weights may be chosen after the subset is found. The proof introduces correlated samples with prescribed marginals and controls their relative entropy together with the information in the noisy sum. Two orders of refining their conditional laws reach the required phase concentration with a small loss of frequencies (Sections 3 and 4).

  2. 2.

    Two bounds for the cost of a new observation. For an input point Y, a label has the form Z=Y+\theta on the circle, where the noise density is proportional to e^{-A} and A is a nonnegative linear combination of functions 1-\cos(2\pi sy), s\in\mathbb Z. Nonnegative Fourier coefficients and the Griffiths–Ginibre positive-correlation method [10] control how this density changes when a new term is added to A. If too many large conditional Fourier coefficients are separated, compression supplies such an addition. It either produces an event unlikely under a comparison law, or enlarges the set of frequencies whose phases vary little under the noise. The first alternative is bounded by likelihood information, the second by a packing estimate. Together they bound the total cost of the labels (Sections 5 and 6).

  3. 3.

    Separate costs for information and for choosing frequencies. A part of the tree with fixed noise law and mass normalization is called a live plateau; its normalization is the mean of its input piece at entry. On each plateau a cumulative list approximates the relevant frequencies by signed sums. At selected spatial intervals we store a smaller list of representatives and the current mass measures divided by that normalization; these records are the snapshots. They are refreshed before the next output bit whenever those measures have changed enough.

    For a representative frequency, the path estimate controls a transform whose values are functions of the label. Integrating against a label character recovers the scalar transform at a shifted frequency, multiplied by the corresponding Fourier coefficient of the noise density. The choice of representative makes this coefficient close to one; the remaining frequency discrepancy contributes a small spatial phase error, which sums geometrically along the path.

    The delayed vector estimate charges label information with an absolute coefficient; the logarithms of the representative counts appear only in a separate probability tail. The distinction is needed when the local estimates are summed. The expected costs of choosing representatives and the remaining path terms, under uniform output sampling, are O(d/h) on each live plateau. The total starting weight of these plateaus is O(h), where a starting weight is the uniform probability of its root interval relative to the top interval, multiplied by its mass normalization. Their product is O(d). Information and errors from input near dyadic interval boundaries are instead charged directly across all plateaus (Sections 7, 8 and 9).

1.3 Organization

Sections 3 and 4 prove the entropy compression lemma. Section 5 constructs and estimates the label updates. Section 6 states the finite dyadic estimate, defines its splitting procedure, and proves the global budgets. Section 7 constructs frequency lists, controls exceptional edges, and fixes the timing of snapshots. Section 8 proves the scalar and vector path estimates, and Section 9 assembles the dyadic bound. Section 10 transfers that bound to Fourier sums and completes the layer summation.

2 Conventions and probability notation

After stating Theorem 1.1, we use unit-period coordinates y=x/(2\pi) and normalized length dy on \mathbb R/\mathbb Z. Characters are

\chi_u(y)=e^{2\pi iuy},\qquad u\in\mathbb Z.

For real u the same exponential notation is used on real intervals; only integer frequencies are regarded as characters of the circle. Periodic functions are extended to the real line when restricted to spatial intervals. All conclusions are unchanged by modifications on null sets.

All logarithms are natural. Write \log^-t=\max\{0,-\log t\} for t>0. The parameter d is sufficiently large; its lower threshold may depend on fixed auxiliary constants chosen earlier. Numerical constants denoted by c,C can change from line to line. Any other dependence is displayed or explicitly stated. No constant depends on the spatial depth or on a frequency cutoff. The notation \operatorname{poly}(x_1,\ldots,x_k) denotes a polynomial upper bound with fixed nonnegative coefficients and fixed finite exponents. Dyadic levels of a positive parameter are integer powers of two. Integer roundings that matter are stated explicitly.

For a probability law P on a finite set, its Shannon entropy is

\mathrm H(P)=-\sum_x P(x)\log P(x),

with 0\log0=0. Conditional entropy is averaged over the conditioning law. For general probability measures, relative entropy is

D(P\|Q)=\int\log\!\left(\frac{dP}{dQ}\right)dP

when P\ll Q, with value +\infty otherwise. The same convention applies when the integral diverges. Mutual information is \mathrm I(X;Y)=D(P_{X,Y}\|P_X\otimes P_Y), and conditional mutual information is the corresponding average over the conditioning variable. A random variable may stand for its law inside D or a distance between measures.

We use nonnegativity, the chain rules, and contraction of relative entropy under measurable maps. For background on these quantities and their chain rules, see [6, Chapter 2]. With

\|P-Q\|_{\mathrm{TV}}=\sup_A|P(A)-Q(A)|,

Pinsker’s inequality takes the form \|P-Q\|_{\mathrm{TV}}\le\sqrt{D(P\|Q)/2}. For probability densities this is one half of their L^1 distance. Numerical constants absorb the harmless choice between total variation and L^1 conventions when an estimate is stated with an unspecified constant.

In the entropy compression argument, all Shannon entropies concern finite-valued variables. For a finite set S\subset\mathbb Z, nS is the n-fold sumset with repetitions. Later, [D] denotes the signed span of a finite list D:

[D]=\left\{\sum_{v\in D}\varepsilon_v v: \varepsilon_v\in\{-1,0,1\}\right\}.

Lists retain their order of insertion; repeated coordinates, if present, are treated as separately indexed entries.

The dyadic construction uses two probability laws. The law \mathbf U follows a uniformly chosen spatial output point and therefore has fair binary bits. The true input law \mathbf P is weighted by the current nonnegative input. Additional noisy labels are observations conditional on the input point. The deterministic processing choices depend on the spatial and membership state, not on the realized label noise. These distinctions will be used explicitly in Section 6; they are essential when converting information under the input law into estimates on output paths.

3 Entropy foundations for compression

We first establish a finite probabilistic mechanism for replacing control of relatively short random sums by a common region of phase alignment. The letters in the intermediate construction are allowed to be correlated. The quantity that pays for this correlation is a relative entropy together with the information in a noisy sum; neither term will be used in isolation when coordinates are selected or conditioned.

We use the entropy conventions of Section 2. All Shannon entropies in this section concern finite-valued variables, and sums of letters are taken in \mathbb Z.

For the large parameter d, set

h=\log\log d,\qquad T=\lfloor d h^{-5}\rfloor,\qquad w=(\log d)^{1/4},\qquad L=d\exp(w),\qquad m_*=\lfloor\exp(w^3)\rfloor. \tag{3.1}

In particular, m_*=d^{o(1)}=o(T). Each assertion below is for all sufficiently large d; the threshold may depend on the fixed constants in the hypotheses, but never on the finite sets under consideration.

Lemma 3.1 (Entropy compression). Fix D_0,P_0\ge1. Let q be a strictly positive probability mass function on \mathbb Z. Suppose that, for a positive integer g and a probability measure \omega on \mathbb R/\mathbb Z,

H_q(x):=\frac{q(x)}{q(0)} =\int\chi_{gx}(\theta)\,d\omega(\theta),\qquad \chi_v(\theta)=e^{2\pi i v\theta},\qquad F_0:=-\log q(0)\le D_0d. \tag{3.2}

Let S\subset\mathbb Z be finite and nonempty and assume

\log|nS|\le D_0d h^{P_0}\log d \qquad\text{for every integer }1\le n\le d^2. \tag{3.3}

Suppose d/h^2\le E\le D_0d and, for every probability measure \lambda on S,

-\log\mathbb E_\lambda H_q(X_1+\cdots+X_T-X'_1-\cdots-X'_T)\le E, \tag{3.4}

where the 2T letters in this expectation are independent with common law \lambda. Then there exist x_0\in S and a nonempty subset S'\subset S such that

\log|S|-\log|S'|\le C(D_0)\frac{d^2}{Eh^2} \tag{3.5}

and, for every probability measure v on S',

\int\exp\!\left( -L\sum_{x\in S'}v(x) (1-\operatorname{Re}\chi_{g(x-x_0)}(\theta)) \right)\,d\omega(\theta)\ge e^{-CE}. \tag{3.6}

The constant C in Equation (3.6) is numerical. The constant in Equation (3.5) depends only on D_0; the lower threshold for d may also depend on P_0.

The same subset and center work for every v in Equation (3.6). This uniformity will allow a later procedure to choose a distribution after the subset has been found. We prove Lemma 3.1 in stages. This section constructs the letter laws and the entropy quantities used in the argument; Section 4 performs the refinement and extracts the common region of phases.

3.1 A variational construction of correlated letters

Throughout the proof of Lemma 3.1, let U_S be the uniform law on S, and put

c(u)=-\log H_q(u),\qquad u\in\mathbb Z.

By Equation (3.2), 0<H_q(u)\le1, since it is a positive real characteristic coefficient. Thus c is finite and nonnegative, and c(0)=0.

The optimization below uses finite Gibbs variational duality: subtracting linear constraints from an entropy functional produces a normalized exponential law. Related entropy-minimization arguments appear in [15, Section 2.1]. Here the constraints prescribe the one-letter marginals, and the objective also retains the information in a noisy sum. We give the duality calculation for this precise setting.

Lemma 3.2 (A noisy sum with prescribed marginals). Under the hypotheses of Lemma 3.1, there is a law of X^T=(X_1,\ldots,X_T) and \Gamma\in TS-TS such that every X_i has law U_S, the joint law is invariant under permutations of the letters, and, with Z=\sum_{i=1}^T X_i+\Gamma,

D(X^T\|U_S^{\otimes T})+\mathrm I(X^T;Z) +\mathbb E c(\Gamma)\le E. \tag{3.7}

The construction may be made with Z\in TS.

Proof. Consider probability laws on S^T\times TS, with coordinates (X^T,Z), having uniform one-letter marginals. Minimize

T\mathrm H(U_S)-\mathrm H(X^T\mid Z) +\mathbb E c\!\left(Z-\sum_iX_i\right).

The feasible set is a nonempty compact polytope. Conditional entropy is continuous on finite probability simplices, including their boundaries, so the minimum is attained. It is a convex minimization problem because conditional entropy is concave in the joint law and the other terms are linear or constant.

We give the dual formula, including the marginal constraint. For a tuple \boldsymbol\mu=(\mu_1,\ldots,\mu_T) of probability laws on S, let f(\boldsymbol\mu) be the minimum of -\mathrm H(X^T\mid Z)+\mathbb E c(Z-\sum_iX_i) with these marginals. It is finite everywhere on the product of simplices: a product law for the letters and any fixed z\in TS is feasible and has finite cost. It is convex by mixing feasible laws. At the tuple of uniform laws, a relative interior point of this product, the finite convex function f has a supporting affine function. Equivalently, separating its epigraph from points strictly below its value yields finite multipliers for the marginal constraints. Therefore there is no duality gap, and the desired minimum equals

\begin{split} T\mathrm H(U_S)+\sup_{a_1,\ldots,a_T:S\to\mathbb R} \left\{\sum_{i=1}^T\mathbb E_{U_S}a_i -\log\max_{z\in TS} \sum_{\boldsymbol x\in S^T} e^{\sum_i a_i(x_i)} H_q\!\left(z-\sum_i x_i\right)\right\}. \end{split} \tag{3.8}

To check the expression inside this dual, set

B_z=\sum_{\boldsymbol x\in S^T} e^{\sum_i a_i(x_i)}H_q\!\left(z-\sum_i x_i\right)>0.

For a conditional law \pi of X^T given Z=z, its entropy and cost, after subtracting the multipliers, equal

\sum_{\boldsymbol x}\pi(\boldsymbol x) \log\frac{\pi(\boldsymbol x)}{ e^{\sum_i a_i(x_i)}H_q(z-\sum_i x_i)} =D(\pi\|B_z^{-1}e^{\sum_i a_i}H_q(z-\textstyle\sum_i x_i)) -\log B_z.

The minimum is -\log B_z, attained at the displayed normalized Gibbs law. The remaining optimization over the distribution of Z selects a maximizer of B_z, giving -\log\max_z B_z as in Equation (3.8).

The dual objective is concave in the tuple (a_1,\ldots,a_T): the negative logarithm in Equation (3.8) is the negative of a maximum of convex log-sum-exponential functions. The objective is also invariant under coordinate permutations. Averaging over all such permutations can only increase it. Hence it suffices to consider a_i=a for every i. Write

A_a=\sum_{x\in S}e^{a(x)},\qquad \lambda(x)=e^{a(x)}/A_a.

For these multipliers, the full expression in Equation (3.8) is

-T D(U_S\|\lambda) -\log\max_{z\in TS} \mathbb E_{\lambda^{\otimes T}} H_q\!\left(z-\sum_iX_i\right).

The maximum is at least its average when z is the sum of an independent T-tuple with law \lambda. That average is at least e^{-E} by Equation (3.4); the two independent tuples can be interchanged to obtain the sign convention there. Since relative entropy is nonnegative, the displayed dual value is at most E. This proves the same upper bound for the primal minimum.

Finally average a minimizing joint law of (X^T,Z) over all permutations of X^T. Convexity preserves the upper bound, and both the constraints and the cost are permutation invariant. Define \Gamma=Z-\sum_iX_i\in TS-TS. The identity

D(X^T\|U_S^{\otimes T})+\mathrm I(X^T;Z) =T\mathrm H(U_S)-\mathrm H(X^T\mid Z)

then gives Equation (3.7) and the required exchangeability, including the error variable.

3.2 Rules for keeping the same error

For a probability law p on a finite set of integer letters and a law of (X^s,\Gamma) supported on the support of p in every coordinate, define

\begin{split} \mathcal C_p(X^s,\Gamma) &:=D(X^s\|p^{\otimes s}) +\mathrm I(X^s;Z)+\mathbb E c(\Gamma)\\ &=\mathbb E[-\log p^{\otimes s}(X^s)] -\mathrm H(X^s\mid Z)+\mathbb E c(\Gamma), \qquad Z=\sum_iX_i+\Gamma. \end{split} \tag{3.9}

When p has zeros, all statements concern its support, on which the cross entropy is finite. The cost is nonnegative and is convex in the joint law of (X^s,\Gamma) for fixed p. For convexity, the map to (X^s,Z) is deterministic, the cross entropy and error cost are linear, and -\mathrm H(X^s\mid Z) is convex.

Lemma 3.3 (Projection and conditioning rules). Let the letters and error in Equation (3.9) be finite-valued, with p positive on every coordinate support. The following operations apply to that cost.

  1. (i)

    Keeping any fixed subset of the letters, in any fixed order, and adding the same error \Gamma to their sum does not increase \mathcal C_p.

  2. (ii)

    If J is any finite-valued variable, then the average cost of the laws conditional on J is at most the original cost plus \mathrm H(J). After this conditioning, the projection and reordering in (i) may depend on J.

  3. (iii)

    Conditioning on a fixed prefix of the letters and keeping the remaining k letters, with the same error, has no additional average cost. If those letters have a common conditional marginal \nu at a given prefix value, their comparison product may be changed from p^{\otimes k} to \nu^{\otimes k} without increasing the conditional cost.

Proof. For (i), write B and O for the kept and omitted vectors, and let Z_B=Z-\sum O. The cross entropy with respect to the product comparison separates over B,O. Also

\mathrm H(B\mid Z_B)\ge\mathrm H(B\mid Z,O).

It follows from the chain rule that

\begin{align*} \mathcal C_p(B,O;\Gamma)-\mathcal C_p(B,\Gamma) &=\mathbb E[-\log p^{\otimes |O|}(O)] -\mathrm H(B,O\mid Z)+\mathrm H(B\mid Z_B)\\ &\ge\mathbb E[-\log p^{\otimes |O|}(O)]-\mathrm H(O\mid Z) \ge0.\end{align*}

The last inequality follows by first removing the conditioning on Z and then using nonnegativity of D(O\|p^{\otimes |O|}). Reordering does not change either the sum or the product comparison.

For (ii), the cross entropy and error cost are unchanged after averaging conditional laws, whereas the conditional entropy becomes \mathrm H(X^s\mid Z,J). Thus the exact average increase is

\sum_j\mathbb P(J=j)\mathcal C_p(X^s,\Gamma\mid J=j) -\mathcal C_p(X^s,\Gamma) =\mathrm I(X^s;J\mid Z)\le\mathrm H(J). \tag{3.10}

This argument does not impose independence on J: it may depend on the letters, the error, and additional randomness. Apply (i) within each conditional law to obtain the final assertion.

For (iii), let B be the revealed prefix and O the remaining vector. Conditional on B, replacing Z by Z-\sum B is an invertible translation of the output. Consequently

\begin{align*} &\mathcal C_p(B,O;\Gamma) -\mathbb E_B\mathcal C_p(O,\Gamma\mid B)\\ &\hspace{1cm} =\mathbb E[-\log p^{\otimes |B|}(B)]-\mathrm H(B\mid Z) \ge0.\end{align*}

In a conditional law where each remaining letter has marginal \nu, changing the product comparison decreases its cross entropy by |O|D(\nu\|p). The other two terms are unchanged. This proves the claim, including when \nu vanishes at some letters.

3.3 An exchangeable law with nearly independent short blocks

Let T_0=\lfloor T/4\rfloor. Conditioning an exchangeable vector on a suitable prefix both preserves a long exchangeable remainder and makes short subblocks close to independent copies of its conditional marginal. The next lemma retains the noisy-sum cost at the same time.

Lemma 3.4 (Selection of a conditional marginal). Under the hypotheses of Lemma 3.1, there is a probability law \nu on S and a joint law \Psi of an exchangeable T_0-tuple with common marginal \nu and an error \Gamma\in TS-TS such that, writing Z=\sum_{i=1}^{T_0}X_i+\Gamma,

\begin{split} &\log|S|-\mathrm H(\nu)\le CE/T,\\ &D(X^{T_0}\|\nu^{\otimes T_0}) +\mathrm I(X^{T_0};Z)+\mathbb E c(\Gamma)\le CE,\\ &\|\mathcal L_\Psi(X^m)-\nu^{\otimes m}\|_{\mathrm{TV}} \le C m_*\sqrt E/T\qquad(1\le m\le m_*). \end{split} \tag{3.11}

All constants here are numerical.

Proof. Start with Lemma 3.2 and put

s_j=\mathrm H(X_{j+1}\mid X_1,\ldots,X_j),\qquad 0\le j<T.

These numbers are nonincreasing. Indeed conditioning decreases entropy, and exchangeability of the unrevealed coordinates gives

s_{j+1}\le\mathrm H(X_{j+2}\mid X_1,\ldots,X_j)=s_j.

The chain rule and Equation (3.7) also give

\sum_{j=0}^{T-1}(\log|S|-s_j) =D(X^T\|U_S^{\otimes T})\le E. \tag{3.12}

The nonnegative deficits in this sum are nondecreasing. Therefore, for j<T,

\log|S|-s_j\le\frac{E}{T-j}.

Let a=\lfloor T/4\rfloor, b=\lfloor T/2\rfloor, and choose a prefix length j uniformly in the integer interval [a,b]. Since m_*=o(T), we may assume b+m_*\le3T/4. For all indices in the enlarged interval [a,b+m_*], the entropy deficit, and hence the total variation of the sequence s_j on that interval, is at most CE/T.

For a fixed prefix realization x^j of positive probability, write \nu_{j,x^j} for the conditional law of X_{j+1}. It is the conditional marginal of every unrevealed coordinate. The expected divergence of the next m letters from its product is exactly

\begin{split} &\mathbb E_{X^j} D\!\left(\mathcal L(X_{j+1},\ldots,X_{j+m}\mid X^j) \,\middle\|\,\nu_{j,X^j}^{\otimes m}\right)\\ &\hspace{2cm}=m s_j-\sum_{i=0}^{m-1}s_{j+i} =\sum_{i=0}^{m-1}(s_j-s_{j+i}). \end{split} \tag{3.13}

For 0\le i<m_*, summing the differences over j cancels all interior terms:

\sum_{j=a}^{b}(s_j-s_{j+i}) =\sum_{j=a}^{a+i-1}s_j-\sum_{j=b+1}^{b+i}s_j \le i(s_a-s_{b+i})\le C iE/T.

The equality for i=0 means that both sides are zero; for i>0 it is valid because i\le b-a+1 for sufficiently large d. Averaging Equation (3.13) over j, with m=m_*, now gives

\mathbb E_{j,X^j} D\!\left(\mathcal L(X_{j+1},\ldots,X_{j+m_*}\mid X^j) \,\middle\|\,\nu_{j,X^j}^{\otimes m_*}\right) \le C m_*^2E/T^2. \tag{3.14}

The expectation of \log|S|-\mathrm H(\nu_{j,X^j}) is \log|S|-s_j\le CE/T. In addition, for each fixed j, Lemma 3.3 allows us to condition on the prefix, keep T_0 of the remaining coordinates, and change the comparison to \nu_{j,X^j}^{\otimes T_0} with cost at most E after averaging over X^j. Averaging these bounds over j retains the same upper bound.

All three quantities just bounded are nonnegative. Dividing each by its stated positive upper scale and averaging their sum shows that some common j and positive-probability prefix realization meet all three bounds, with a larger numerical constant. Freeze this choice, set its conditional marginal equal to \nu, and keep the first T_0 remaining letters. Conditional exchangeability of the unrevealed letters, including their joint law with \Gamma, gives the claimed exchangeable law \Psi. For the m_*-letter marginal, Pinsker’s inequality \|P-Q\|_{\mathrm{TV}}\le\sqrt{D(P\|Q)/2} gives the last line of Equation (3.11); projection contracts total variation and proves it for every smaller m.

3.4 Contexts and bounded exact-sum costs

Our remaining aim is to reach a longer block length with small noisy cost while losing little one-letter entropy under conditioning. The refinement will be controlled by a bounded exact-sum cost that decreases as information is revealed.

Fix the one-letter marginal X\sim\nu supplied by Lemma 3.4. A context is a finite-valued variable W jointly distributed with X. We start with a constant, denoted by \varnothing. A refinement replaces W by (W,A), where A is sampled by a finite channel conditional on the current pair (X,W); thus it keeps the old context and never changes the marginal of X.

For s independent copies of the current pair (X,W), define

G_s(W)=\mathrm H\!\left(\sum_{i=1}^{s}X_i\,\middle|\,W_1,\ldots,W_s\right), \qquad G_0(W)=0. \tag{3.15}

The independence in this definition is part of the notation. In contrast, the next quantity permits arbitrary dependence between its rows:

I_r(W)=\inf\left\{ r\mathrm H(X\mid W) -\mathrm H\!\left(X_1,\ldots,X_r\,\middle|\, \sum_{a=1}^{r}X_a,W_1,\ldots,W_r\right) \right\}. \tag{3.16}

The infimum is over joint laws of the r pairs (X_a,W_a), each with the prescribed one-row marginal (X,W). All these optimizations are over compact finite probability polytopes with continuous objectives, so their infima are attained.

We will use the following observation for several conditional costs. Suppose a trial has prescribed row marginals (B_a,V_a) and a possibly random output Z. Its conditional entropy cost is

\mathcal F=\sum_a\mathrm H(B_a\mid V_a) -\mathrm H(B^r\mid Z,V^r)\ge0.

Nonnegativity follows by conditioning and subadditivity; equivalently,

\mathcal F= \mathbb E_{V^r}D\!\left(\mathcal L(B^r\mid V^r) \,\middle\|\,\prod_a\mathcal L(B_a\mid V_a)\right) +\mathrm I(B^r;Z\mid V^r).

To extend this trial to new contexts, sample A_a independently across rows, conditional on their old pairs, and independently of Z given all the old pairs. Thus

\mathcal L(A^r\mid B^r,V^r,Z) =\prod_a Q_a(\,\cdot\mid B_a,V_a).

Let \mathcal F' be the cost after the contexts become (V_a,A_a). The chain rule gives

\begin{split} \mathcal F-\mathcal F' &=\sum_a\mathrm I(B_a;A_a\mid V_a) -\mathrm I(B^r;A^r\mid Z,V^r)\\ &=\sum_a\mathrm H(A_a\mid V_a) -\mathrm H(A^r\mid Z,V^r)\ge0. \end{split} \tag{3.17}

For the second equality, conditional independence cancels the terms \mathrm H(A_a\mid B_a,V_a) against \mathrm H(A^r\mid B^r,Z,V^r). The last inequality is again conditioning and subadditivity. This calculation applies even when the old rows are highly correlated.

Lemma 3.5 (Monotonicities and uniform entropy bounds). Let \nu,\Psi be as in Lemma 3.4, and define G_s,I_r by Equations (3.15) and (3.16). There is a constant C'_0(D_0) such that, with

D_*=C'_0(D_0)d, \tag{3.18}

every finite refinement of the initial constant context satisfies:

  1. (i)

    G_s(W) is nondecreasing in s, nonincreasing under refinement of W, and 0\le G_s(W)\le D_* for 0\le s\le m_*.

  2. (ii)

    I_r(W) is nonnegative, nondecreasing in r, and nonincreasing under refinement of W. Moreover I_r(W)\le D_* for 1\le r\le T_0.

The same D_* works for all the contexts in any finite sequence of refinements.

Proof. First consider sums under the correlated law \Psi, before introducing contexts. By Lemma 3.3, for every s\le T_0 the law of (X^s,\Gamma) has \mathcal C_\nu\le CE. Relative entropy to q gives

\mathrm H(\Gamma)\le\mathbb E[-\log q(\Gamma)] =F_0+\mathbb E c(\Gamma). \tag{3.19}

This comparison is legitimate although q lives on all of \mathbb Z: \Gamma is finite-valued and q is positive on its support. Write Z_s=\sum_{i=1}^sX_i+\Gamma. Given X^s, the two variables Z_s and \Gamma are translates of one another. Hence

\mathrm H(Z_s)=\mathrm I(X^s;Z_s)+\mathrm H(\Gamma\mid X^s) \le\mathrm I(X^s;Z_s)+\mathrm H(\Gamma).

Since \sum_iX_i=Z_s-\Gamma, it follows that

\mathrm H_\Psi\!\left(\sum_{i=1}^sX_i\right) \le\mathrm I(X^s;Z_s)+2\mathrm H(\Gamma) \le 2F_0+CE. \tag{3.20}

We transfer this bound to independent letters only for s\le m_*. The sum laws then have total variation distance at most \eta=Cm_*\sqrt E/T by Equation (3.11). For completeness, if two laws on a finite alphabet of size N have distance \eta', couple them as U,V with \mathbb P(U\ne V)=\eta'. Such a coupling uses their common submeasure, the pointwise minimum, for the equal values; the two residual submeasures have disjoint supports and total mass \eta'. Let J be the disagreement indicator. The chain rule gives

\mathrm H(U)\le\mathrm H(V)+\mathrm H(J)+\mathrm H(U\mid V,J) \le\mathrm H(V)+\log2+\eta'\log N.

Interchanging U,V yields the corresponding absolute difference bound. Apply this with the common alphabet sS. The added error is at most

\log2+C\frac{m_*\sqrt E}{T}\log|sS|=o(d/h^2). \tag{3.21}

Indeed m_*\sqrt E/T\le d^{-1/2+o(1)} for fixed D_0, while Equation (3.3) bounds the logarithm by D_0d h^{P_0}\log d; the product is d^{1/2+o(1)}. This estimate is uniform in s\le m_*, and s\le d^2 for large d. Combining Equations (3.20) and (3.21) shows G_s(\varnothing)\le C(D_0)d.

For independent current pairs, conditioning on all the contexts makes the letters independent. For any independent finite-valued A,B on \mathbb Z, one has \mathrm H(A+B)\ge\mathrm H(A+B\mid B)=\mathrm H(A). Applying this at each fixed context tuple and averaging proves that G_s(W) is nondecreasing in s. Adding the last context does not change the law of the first s pairs, so the averaged right-hand side is exactly G_s(W). Refinement decreases G_s by conditioning on more information in an iid extension of all pairs. The established initial bound therefore proves (i).

Nonnegativity of I_r follows from that of \mathcal F. For its monotonicity in r, take any (r+1)-row trial and put S_r=\sum_{a=1}^rX_a. The chain rule and removal of conditioning give

\begin{align*} &\mathrm H(X^{r+1}\mid S_{r+1},W^{r+1})\\ &\quad=\mathrm H(X_{r+1}\mid S_{r+1},W^{r+1}) +\mathrm H(X^r\mid X_{r+1},S_r,W^{r+1})\\ &\quad\le\mathrm H(X\mid W)+\mathrm H(X^r\mid S_r,W^r).\end{align*}

The projected first r pairs form an admissible trial. Subtracting this inequality from (r+1)\mathrm H(X\mid W) and taking infima proves I_{r+1}(W)\ge I_r(W). To refine a trial for I_r, use independent channels given each old pair as in Equation (3.17), with Z=S_r. They preserve the required new one-row marginal and do not increase the objective. Taking an optimal old trial proves refinement monotonicity.

Finally use the projection of \Psi to r\le T_0 letters as a trial at the constant context. Its exact-sum objective is

r\mathrm H(\nu)-\mathrm H(X^r\mid S_r) =D(X^r\|\nu^{\otimes r})+\mathrm H(S_r) \le 2F_0+CE,

using Equations (3.11) and (3.20) and Lemma 3.3. Choose C'_0(D_0) large enough for this and the G_s bound. Refinement monotonicity then proves (ii) with the same fixed D_*.

3.5 Noisy block costs and their full-vector interpretation

For m independent current pairs, put Y_m=\sum_{i=1}^mX_i and V_m=(W_1,\ldots,W_m). At a fixed positive integer r, define

\begin{split} K_m=K_m(W;r):=\inf\bigg\{ &r\mathrm H(Y_m\mid V_m) -\mathrm H\!\left(Y_{m,1},\ldots,Y_{m,r}\,\middle|\, \sum_{a=1}^rY_{m,a}+\Gamma,V_m^r\right)\\ &\hspace{3cm}+\mathbb E c(\Gamma)\bigg\}. \end{split} \tag{3.22}

Here each row pair (Y_{m,a},V_{m,a}) has the prescribed marginal (Y_m,V_m), the rows may be arbitrarily coupled, and the error may be arbitrarily coupled with them subject to \Gamma\in TS-TS. The notation K_m suppresses W and r only when both are fixed. As before, the infimum is attained and the objective is nonnegative.

Lemma 3.6 (Lifting block sums and refining contexts). Let m,r be positive integers and W a finite context of X\sim\nu. The value of Equation (3.22) is unchanged if each row sum Y_{m,a} is replaced in both entropy terms by its entire m-letter vector, with that row constrained to have the law of m iid current pairs. The output remains the sum of all mr letters plus \Gamma, and independence is required within each prescribed row marginal, but not between rows. Moreover K_m(W;r) is nonincreasing under finite context refinement.

Proof. Write B_a=(X_{a,1},\ldots,X_{a,m}), V_a=(W_{a,1},\ldots,W_{a,m}), Y_a=\sum_iX_{a,i}, and Z=\sum_aY_a+\Gamma. For any full-vector trial, the full-vector objective minus its projected sum objective equals

\sum_{a=1}^r\mathrm H(B_a\mid Y_a,V_a) -\mathrm H(B^r\mid Y^r,Z,V^r)\ge0. \tag{3.23}

The equality follows from the chain rule, since Y_a is a function of B_a; the inequality follows by subadditivity and conditioning. Thus projection to block sums does not increase the cost.

Conversely, start with any admissible joint law of (Y^r,V^r,\Gamma). In each row, independently of all the other new rows and of \Gamma conditional on (Y^r,V^r), sample the block from the genuine iid block law conditioned on its own sum and contexts:

\mathcal L(B^r\mid Y^r,V^r,\Gamma) =\prod_{a=1}^r \mathcal L_{\mathrm{iid}}(B_a\mid Y_a,V_a).

Each full row then has precisely the prescribed iid marginal. Conditioning further on Z, a function of (Y^r,\Gamma), does not change these independent conditional row laws. Consequently the last entropy in Equation (3.23) equals the sum of the first entropies after averaging. The two objectives are equal for this lift, which proves equality of their infima.

For refinement, use this full-vector formulation. At every one of its mr positions, run the prescribed one-letter channel independently given the old pairs, also independently of the old error given those pairs. The iid old-pair marginal within a row and the independent channels give exactly the iid refined-pair marginal required in that row. Apply Equation (3.17) with the whole row as B_a and its vector of new observations as A_a. The error law and its cost are unchanged, and the entropy objective cannot increase. Taking an optimal old trial proves the stated monotonicity.

A new one-letter channel may itself be specified by an experiment using auxiliary letters, as will occur in Section 4. In the construction just used, each channel is run with fresh auxiliary randomness conditional on its own old pair. One does not substitute the other actual letters of its row for those auxiliary samples. This distinction is what preserves the prescribed iid row marginal.

3.6 Thinning correlated blocks to exact iid marginals

The short blocks of \Psi are only approximately iid. A later refinement argument needs exact iid marginals inside each row. We obtain them by selecting blocks with a weighted acceptance rule and filling the missing row mass.

Lemma 3.7 (The low-end noisy cost). Under the hypotheses of Lemma 3.1, fix the law \nu from Lemma 3.4. Let m,r be positive integers satisfying

m\le m_*,\qquad 2mr\le T_0,\qquad r\le \tfrac12 T e^{-w^2}.

For the constant context and for every finite refinement of it,

K_m(W;r)\le CE, \tag{3.24}

where C is numerical, with the threshold for d allowed to depend on D_0,P_0.

Proof. It suffices to construct a trial for the constant context; the rest follows from Lemma 3.6. Take 2r disjoint blocks of m coordinates from \Psi. Their common marginal is a law \mu on S^m. Put Q=\nu^{\otimes m} and

\varepsilon=\|\mu-Q\|_{\mathrm{TV}} \le C m_*\sqrt E/T.

Define the common submeasure \gamma(\boldsymbol x) =\min(\mu(\boldsymbol x),Q(\boldsymbol x)), which has total mass 1-\varepsilon. Conditional on a block value \boldsymbol x with \mu(\boldsymbol x)>0, tag it with probability \gamma(\boldsymbol x)/\mu(\boldsymbol x). The choices of tags may be made independently given all block values. Each tagged block has unnormalized law \gamma; independence between the block values is not assumed.

Let B be the number of tagged blocks. If B<r, reject. If B\ge r, accept with probability B/(2r), using a fresh coin, and upon acceptance select r distinct tagged blocks uniformly in random order. Let the acceptance probability be 1-\delta. Each block fails to be tagged with probability \varepsilon, so

\mathbb E(2r-B)=2r\varepsilon, \qquad\mathbb P(B<r)\le2\varepsilon.

The rule consequently gives

\begin{split} \delta &=1-\mathbb E\left[\frac B{2r}\mathbf1_{\{B\ge r\}}\right]\\ &=\varepsilon+\mathbb E\left[\frac B{2r}\mathbf1_{\{B<r\}}\right] \le3\varepsilon. \end{split} \tag{3.25}

In particular \delta<1 for sufficiently large d.

We check the marginal of an accepted row before normalizing by the acceptance probability. Given the tags, a particular tagged block has probability 1/B of occupying any specified selected row. Multiplying this by the acceptance probability B/(2r) cancels B. Therefore all selected rows have the same unnormalized marginal \eta, and, for each \boldsymbol x\in S^m,

\begin{split} \eta(\boldsymbol x) &=\frac1{2r}\sum_{i=1}^{2r} \mathbb E\!\left[ \mathbf1_{\{B\ge r\}} \mathbf1_{\{i\text{ tagged}\}} \mathbf1_{\{\text{block }i=\boldsymbol x\}} \right]\\ &\le\gamma(\boldsymbol x)\le Q(\boldsymbol x). \end{split} \tag{3.26}

Its total mass is 1-\delta. This domination is the reason for the acceptance weight B/(2r).

On acceptance, retain the selected letters in their selected order and keep the same error \Gamma. To bound their cost, first project \Psi to the 2rm letters used above. Its cost with respect to \nu^{\otimes 2rm} is at most CE. Record a finite index J that is either a rejection symbol or the ordered tuple of the r selected block indices. The full tag pattern is not recorded. Thus

\mathrm H(J)\le\log\bigl(1+(2r)^r\bigr) \le\log2+r\log(2r).

By Lemma 3.3, conditioning on J and then projecting and reordering on accepted outcomes has average cost at most CE+\mathrm H(J). Nonnegativity permits us to drop the rejected outcomes from this upper bound. Finally, forgetting the ordered indices on the accepted outcomes uses convexity of the combined cost \mathcal C_\nu. If P_{\mathrm{acc}} is the normalized accepted joint law of the mr letters and the error, we obtain

(1-\delta)\mathcal C_\nu(P_{\mathrm{acc}}) \le CE+\log2+r\log(2r). \tag{3.27}

Keeping the acceptance mass on the left avoids any division by a small conditional acceptance probability for a particular ordered selection.

If \delta=0, the accepted rows already have marginal Q and the construction is complete. Otherwise let

R=\frac{Q-\eta}{\delta}.

This is a probability law on S^m by Equation (3.26). Complete the trial by mixing P_{\mathrm{acc}}, with weight 1-\delta, and the law of r independent blocks with marginal R, with weight \delta; in the latter law set \Gamma=0. Every row of the mixture now has exact marginal \eta+\delta R=Q. Thus it is an admissible full-vector trial for K_m. The error still lies in TS-TS, which contains zero.

For the filling law, R\le Q/\delta pointwise, and hence

D(R\|Q)\le\log(1/\delta).

Independence between its r blocks makes its divergence to \nu^{\otimes mr}=Q^{\otimes r} equal to rD(R\|Q). Its noisy sum is the exact sum, which lies in mrS, so the mutual information in the output is its entropy, at most \log|mrS|. The error cost is c(0)=0. Its weighted additional cost is therefore at most

\delta r\log(1/\delta)+\delta\log|mrS|. \tag{3.28}

Convexity of \mathcal C_\nu on mixing, followed by the full-vector identity of Lemma 3.6, bounds K_m by the sum of the right-hand sides of Equations (3.27) and (3.28).

It remains to verify the scale of all additional terms. Since w^2=\sqrt{\log d} and r\le Te^{-w^2}/2,

r\log(2r) \le Cdh^{-5}e^{-\sqrt{\log d}}\log d=o(d/h^2).

Also \delta\log(1/\delta)\le1/e, so the first term of Equation (3.28) is at most r/e=o(d/h^2). For its second term, Equation (3.25) gives \delta\le d^{-1/2+o(1)}, and mr\le T_0/2<d^2 allows the sumset hypothesis to give

\delta\log|mrS| \le d^{-1/2+o(1)}D_0d h^{P_0}\log d=o(d/h^2).

These estimates include all entropy paid for the selection index. They are uniform in admissible m,r. Since E\ge d/h^2, the additional terms can be absorbed into a numerical multiple of E once d is sufficiently large depending on D_0,P_0. This proves Equation (3.24).

We now have a fixed one-letter law \nu, uniformly bounded exact-sum quantities G_s and I_r, and the small noisy cost Equation (3.24) at the low block lengths. The next section uses context refinements to reach a substantially larger block length while retaining a comparable noisy cost, and then completes the proof of Lemma 3.1.

4 Refinement and common phase volume

We complete the proof of Lemma 3.1. All parameters and costs retain their meanings from Section 3; in particular, X\sim\nu, G_0=0, T_0=\lfloor T/4\rfloor, and D_*=C'_0(D_0)d bounds the relevant G_s and I_r. Choose C'_0(D_0) large enough that E\le D_* as well. We occasionally write K_m(W;r) to display the context and row count. A prime denotes a value after the specified refinement.

First, a block observation converts a noisy-cost difference into a decrease of an exact-sum cost. A finite search then finds a long block with small noisy cost while spending little one-letter entropy. Finally, we extract a common region of phases and a large subset of the original letters.

4.1 Gluing two couplings

Lemma 4.1 (A block observation pays a noisy-cost difference). Fix a current finite-valued context W, integers m\ge2 and r\ge1, and an independent block of m current pairs (X_i,W_i). Refine the context of X_1 by observing \bigl(\sum_{i=1}^mX_i,W_2,\ldots,W_m\bigr), retaining W_1. For the resulting context W',

G_1(W)-G_1(W')=G_m(W)-G_{m-1}(W), \tag{4.1}

and

I_r(W)-I_r(W')\ge K_m(W;r)-K_{m-1}(W;r). \tag{4.2}

If the observation is made only on an independent Bernoulli mark of probability \zeta\in[0,1], retaining the mark in the new context, its one-letter entropy loss is multiplied by \zeta, and its decrease of I_r is at least \zeta\bigl(I_r(W)-I_r(W')\bigr).

Proof. Independence of the pairs gives

\begin{split} &\mathrm I\left(X_1;\sum_{i=1}^mX_i,W_2,\ldots,W_m\mid W_1\right)\\ &\quad=\mathrm H\left(\sum_{i=1}^mX_i\mid W_1,\ldots,W_m\right) -\mathrm H\left(\sum_{i=2}^mX_i\mid W_2,\ldots,W_m\right), \end{split}

which proves Equation (4.1).

Take independent minimizing trials for I_r(W) and K_{m-1}(W;r). Their alphabets, including the error alphabet TS-TS, are finite, so minimizers exist. Write their row vectors and contexts as

U=(U_a)_{a=1}^r,\quad\mathcal W=(W_a)_{a=1}^r, \qquad Y=(Y_a)_{a=1}^r,\quad\mathcal V=(V_a)_{a=1}^r,

where each (Y_a,V_a) has the sum-and-context law of an independent (m-1)-letter block. Set

P=(\mathcal W,\mathcal V),\qquad Z_U=\sum_aU_a,\qquad Z_Y=\sum_aY_a+\Gamma.

The whole (U,\mathcal W) trial is independent of the whole (Y,\mathcal V,\Gamma) trial. Consequently (U+Y,P,\Gamma) is admissible for K_m(W;r), and (U,(U+Y,P)), with exact output Z_U, is admissible for I_r(W'). Each row has the prescribed marginal in both trials; independence across rows is not required.

The two positive one-row entropy terms add to

\begin{split} &\mathrm H(U_a+Y_a\mid W_a,V_a) +\mathrm H(U_a\mid U_a+Y_a,W_a,V_a)\\ &\quad=\mathrm H(U_a,Y_a\mid W_a,V_a) =\mathrm H(U_a\mid W_a)+\mathrm H(Y_a\mid V_a). \end{split}

For the subtracted joint entropies, condition both terms on both old outputs:

\begin{split} &\mathrm H(U+Y\mid Z_U+Z_Y,P)+\mathrm H(U\mid U+Y,Z_U,P)\\ &\quad\ge\mathrm H(U+Y\mid Z_U,Z_Y,P) +\mathrm H(U\mid U+Y,Z_U,Z_Y,P)\\ &\quad=\mathrm H(U,Y\mid Z_U,Z_Y,P)\\ &\quad=\mathrm H(U\mid Z_U,\mathcal W) +\mathrm H(Y\mid Z_Y,\mathcal V). \end{split} \tag{4.3}

The middle equality uses the invertible map (U,Y)\mapsto(U+Y,U); the last uses independence of the old trials. This argument permits random \Gamma. The error cost stays unchanged. Adding the new trial costs and taking infima proves K_m(W;r)+I_r(W')\le K_{m-1}(W;r)+I_r(W).

For a partial observation W^{(\zeta)}, conditioning on the independent mark gives its entropy loss. To bound its exact cost, synchronize the mark across rows, using a minimizing old coupling on mark zero and a minimizing fully refined coupling on mark one. This has the prescribed refined pair marginal in each row, so

I_r(W^{(\zeta)})\le(1-\zeta)I_r(W)+\zeta I_r(W').

The assertion follows.

4.2 A search with two orders of refinement

We choose block ranges whose lower endpoints decrease as their row counts increase. For sufficiently large d, put

Q=\left\lfloor\frac{w^3/2-w^2}{3w}\right\rfloor,\qquad M=Q+1,

and, for 1\le j\le M, set

\begin{aligned} m_{j,0}&=\left\lceil\exp\bigl(w^2+3w(M-j)\bigr)\right\rceil,\\ r_j&=\left\lfloor\frac{T_0}{2m_{j,0}}\right\rfloor, &m_{j,1}&=\left\lceil\frac L{r_j}\right\rceil . \end{aligned} \tag{4.4}

Here M\asymp w^2. Uniformly in j, T_0/(2m_{j,0})\to\infty, so

\frac{T_0}{4m_{j,0}}\le r_j\le\frac{T_0}{2m_{j,0}},\qquad 1<\frac{m_{j,1}}{m_{j,0}} \le\frac{4L}{T_0}+\frac1{m_{j,0}}\le e^{2w}.

Successive lower endpoints have ratio at least e^{3w}/(1+e^{-w^2})>e^{2w}. Thus the integer ranges [m_{j,0},m_{j,1}] are disjoint and descend as j increases. The floor bounds also show that r_j increases strictly. Moreover,

\begin{gathered} e^{w^2}\le m_{j,0}<m_{j,1}\le m_*, \qquad 2r_jm_{j,0}\le T_0,\qquad r_j\le Te^{-w^2}/8,\\ L\le r_jm_{j,1}<L+r_j<2L<d^2. \end{gathered} \tag{4.5}

For the upper endpoint, use m_{j,1}\le2\exp(w^3/2+2w)\le m_* for large d. The low endpoints therefore satisfy all hypotheses of Equation (3.24), and all the G_s,I_{r_j} used below are bounded by D_*.

Proposition 4.2 (Finite refinement search). Under the hypotheses of Lemma 3.1, there is a finite sequence of context refinements from empty context to W, and an index j, such that

\begin{split} G_1(\varnothing)-G_1(W)&\le\frac{CD_*^2}{Eh^2},\\ K_{m_{j,1}}(W;r_j)&\le K_{m_{j,0}}(W;r_j)+E. \end{split} \tag{4.6}

The constant C is numerical, and every context produced is finite-valued.

Proof. Normalize by

a_j=I_{r_j}/D_*,\qquad g_s=G_s/D_*,\qquad\tau=E/D_*\in(0,1].

These functions take values in [0,1]. At a fixed state a_j increases with j and g_s increases with s; all decrease under refinement. Measure refinement cost by \tau times the loss in g_1. Costs add, and total cost c means entropy loss cD_*^2/E.

Available updates. Suppose the desired K inequality fails at index j in a given state. Consider separately from that same state the full observations of Lemma 4.1 for m_{j,0}<m\le m_{j,1}. Let d_m\ge0 be their decreases in a_j, and c_m\ge0 their costs. Telescoping Equations (4.2) and (4.1) gives

\sum_m d_m>\tau,\qquad\sum_m c_m=\tau p_j,\qquad p_j=g_{m_{j,1}}-g_{m_{j,0}}.

For p_j>0, offers with c_m>2p_jd_m have total decrease less than \tau/2. The other offers have total decrease greater than \tau/2. There are at most m_* offers, so one satisfies

d_m\ge\eta:=\frac{\tau}{2m_*},\qquad c_m\le2p_jd_m. \tag{4.7}

If p_j=0, all costs are zero and the same conclusion follows. For any target x\in(0,d_m], use a partial observation with probability x/d_m. Its actual decrease is at least x and its cost at most 2p_jx.

Why two visit orders are needed. The available update decreases a_j at a cost controlled by the current width p_j. Although both endpoints g_{m_{j,0}} and g_{m_{j,1}} decrease under refinement, their difference need not decrease. We will therefore stop working at an index when the difference between the g-value at its upper endpoint on entry and the g-value at its lower endpoint in the current state is too large relative to the cost already spent there.

Two comparisons determine the order of the search. When we visit indices in increasing order, their block-length ranges descend: the g-value at the upper endpoint on entry to the next visited index is bounded by the g-value at the lower endpoint on exit from the earlier index. These endpoint differences telescope in [0,1] and will control the cost of skipped work. Completed decreases of a_j telescope in the opposite index order, because a_j increases with j and decreases under refinement. The nested partition below uses increasing indices within a subblock to control skips, decreasing subblocks within a macroblock to control completed decreases, and increasing macroblocks to control the remaining skips. The numerical thresholds make the total cost O(h^{-2}).

Partitions and visit orders. Put

\begin{gathered} N=\lfloor M^{1/10}\rfloor,\quad n=\lfloor M^{1/2}\rfloor,\quad u=(\log N)^{-2},\quad q_*=\lceil2/u\rceil,\\ \alpha_0=\frac{100u}{\log n},\quad\delta_0=n^{-1/2},\qquad \alpha_1=\frac{10^5}{\log n\log N},\quad\delta_1=N^{-1/2}. \end{gathered} \tag{4.8}

Take d large enough that N\ge4. Partition the first Nq_*n indices into N successive macroblocks, each containing q_* successive subblocks of n indices. There is room since

Nq_*n\le M^{3/5}\left(1+\frac{(\log M)^2}{50}\right)\le M.

Process macroblocks in ascending index order; within a macroblock, process subblocks in descending order; within a subblock, process indices in ascending order. A subblock finishes when a_j at one of its indices has decreased by at least u since entry to that index. Discovery of the desired K inequality stops the whole search. The skip rules below can end an index or a macroblock earlier.

One index and one subblock. At index j, let c_0 be the cost spent there so far and put

p=g_{m_{j,1}}(\text{entry to }j)-g_{m_{j,0}}(\text{current state}).

This nonnegative quantity dominates the current p_j. Skip this index if

p>c_0/\alpha_0+\delta_0. \tag{4.9}

Otherwise check the desired K inequality. If it fails, take an offer from Equation (4.7). Write \Delta_j=a_j(\text{entry to }j)-a_j(\text{current state})<u for the decrease already achieved at this index, and choose the partial observation with target

x=\min\{d_m,u-\Delta_j\}.

Repeat until a skip, completed target, or discovery. The sum of the targets x in an index is at most u, since actual decreases are at least their targets and a clipped last target finishes the attempt. Each unclipped target equals d_m\ge\eta, so the loop is finite.

At every update the skip test has failed, whence

\Delta c_0\le2(c_0/\alpha_0+\delta_0)x,\qquad c_{0,\mathrm{new}}+\alpha_0\delta_0 \le(c_{0,\mathrm{old}}+\alpha_0\delta_0)(1+2x/\alpha_0).

Starting at c_0=0, and using 1+t\le e^t, the cost at every stage satisfies

c_0\le\alpha_0\delta_0(e^{2u/\alpha_0}-1) \le\alpha_0 n^{-12/25}\le\alpha_0 n^{-2/5}. \tag{4.10}

Write H,\ell for the highest and lowest block-length endpoints in the current macroblock and set

P=g_H(\text{entry to the macroblock})-g_\ell(\text{current state})\in[0,1].

For successive skipped indices j<j' within a subblock, the block-length ranges descend, so m_{j',1}<m_{j,0}. Decrease over refinement time then gives

g_{m_{j',1}}(\text{entry to }j') \le g_{m_{j,0}}(\text{exit from }j). \tag{4.11}

Thus the values of p at skips sum to at most P, including at every intermediate state in this subblock. A skipped index has c_0<\alpha_0(p-\delta_0)\le\alpha_0p and p>\delta_0. There are therefore fewer than 1/\delta_0=\sqrt n skips; all n indices cannot be skipped. Skipped attempts cost at most \alpha_0P, and there is at most one other, active or terminal, index attempt. By Equation (4.10), every prefix or completed run of the subblock costs at most

\alpha_0(P+n^{-2/5}). \tag{4.12}

A completed subblock without discovery must consequently reach a target.

Figure 1: Visit orders in the refinement search. Macroblocks advance with j, subblocks inside one macroblock are visited in reverse order, and indices inside one subblock again advance with j. Block lengths m decrease with j. The endpoint telescope controls skipped work, while the exact-cost telescope controls completed targets. In the first sum, j_k runs through skipped indices in one subblock, U_k=g_{m_{j_k,1}} at entry to j_k, and V_k=g_{m_{j_k,0}} at exit from j_k. For a macroblock with highest and lowest block-length endpoints H,\ell, put P=g_H(\text{entry to the macroblock})-g_\ell(\text{current state}). Each P in the macroblock sum is evaluated at that macroblock’s own skip. Entry and exit values include changes caused by intervening refinements.

One macroblock. Let c_1 be the cumulative cost in the current macroblock. At boundaries before or after subblock runs, skip the rest of the macroblock if

P>c_1/\alpha_1+\delta_1. \tag{4.13}

Target indices decrease strictly between completed subblocks. For successive such indices j>j', monotonicity in index and time gives

a_{j'}(\text{entry to }j')\le a_j(\text{exit from }j). \tag{4.14}

Their actual targeted decreases telescope to at most one. At most 1/u subblocks can complete targets, whereas there are q_*>1/u subblocks. Thus this macroblock must be skipped or discover the desired inequality before it can be exhausted.

For a completed subblock after which no macroblock skip is made and processing continues, apply Equation (4.12) with the new P, and then the failure of Equation (4.13):

c_{1,\mathrm{new}}-c_{1,\mathrm{old}} \le\alpha_0(c_{1,\mathrm{new}}/\alpha_1+\delta_1+n^{-2/5}).

Put \beta_*=\alpha_0/\alpha_1=1/(1000\log N) and \delta_*=\delta_1+n^{-2/5}. Iteration over at most 1/u such runs yields

c_1\le\alpha_1\delta_*\bigl((1-\beta_*)^{-1/u}-1\bigr)\le\alpha_1. \tag{4.15}

Indeed n\ge N^5, so n^{-2/5}\le N^{-2}. Also -\log(1-\beta_*)\le2\beta_*, and hence (1-\beta_*)^{-1/u}\le N^{1/500}. The factor multiplying \alpha_1 is at most N^{-249/500}+N^{-999/500}<1 for N\ge4. One last or partial subblock costs at most 2\alpha_0 by Equation (4.12) and P\le1.

Across macroblocks and termination. Macroblocks advance in index order, so their block-length ranges descend as in Equation (4.11). The values of P at macroblock skips sum to at most one. Their total cost is at most \alpha_1, and fewer than 1/\delta_1=\sqrt N macroblocks can be skipped. Since a macroblock cannot be exhausted without a skip or discovery, the whole search must discover the desired K inequality. Its loops are finite by the positive minimum offer before clipping. Total cost through discovery is bounded by

\underbrace{\alpha_1}_{\text{skipped macroblocks}} +\underbrace{\alpha_1}_{\text{continuing subblocks}} +\underbrace{2\alpha_0}_{\text{last subblock}}\le4\alpha_1.

Both telescopes allow decreases caused by intervening refinements; see Figure 1.

Since h=4\log w and M\asymp w^2, both \log n and \log N are of order h. Thus 4\alpha_1\le C/h^2. Undoing the cost normalization proves Equation (4.6). Each new observation and mark has a finite alphabet, so the finite sequence leaves a finite-valued context.

4.3 Freezing contexts and extracting one phase region

Completion of the proof of Lemma 3.1. Use the context W, index j, and high endpoint m=m_{j,1} supplied by Proposition 4.2. Put r=r_j and s=mr\ge L. Combining Equation (4.6) with the low-end bound Equation (3.24), preserved under refinement, gives

K_m(W;r)\le CE \tag{4.16}

with numerical C. Lift a minimizing coupling to all s letters using the full-vector interpretation of the cost. Write \boldsymbol W=(W_i)_{i=1}^s for their whole context stack, and \nu_z for the prescribed one-letter conditional law given context z.

At any fixed stack \boldsymbol w of positive probability, compute

\begin{split} C(\boldsymbol w)&= D\left(X^s\,\middle\|\,\prod_i\nu_{w_i}\right) +\mathrm I(X^s;Z)+\E c(\Gamma),\\ &\hspace{35mm}Z=\sum_iX_i+\Gamma, \end{split} \tag{4.17}

under the actual conditional law at that stack. All three terms are nonnegative. The actual conditional letter marginals after freezing the whole stack may differ from the comparator laws \nu_{w_i}. Nevertheless the comparison has finite divergence: an event X_i=x,W_i=w_i with \nu_{w_i}(x)=0 has zero probability even before conditioning, so it also has zero probability at any positive-probability stack. Before freezing, however, each row has its prescribed independent block marginal, and consequently

\E\left[-\log\prod_i\nu_{W_i}(X_i)\right]=s\mathrm H(X\mid W).

Subtracting \mathrm H(X^s\mid Z,\boldsymbol W) and adding the error cost therefore shows that \E C(\boldsymbol W)=K_m(W;r)\le CE.

For the separate quantity B(\boldsymbol w)=s^{-1}\sum_iD(\nu_{w_i}\|U_S), the prescribed pair marginals give

\begin{split} \E B(\boldsymbol W) &=\log|S|-\mathrm H(X\mid W)\\ &=\log|S|-\mathrm H(\nu)+G_1(\varnothing)-G_1(W)\\ &\le CE/T+\frac{CD_*^2}{Eh^2} \le\frac{C(D_0)d^2}{Eh^2}. \end{split} \tag{4.18}

This uses Equations (3.11) and (4.6). The E/T term is absorbed because T\asymp d/h^5, E\le D_0d, and E^2h^7/d^3\to0 for fixed D_0. Markov’s inequality, with factor four for each quantity, selects one stack such that

C(\boldsymbol w)\le CE,\qquad B(\boldsymbol w)\le\frac{C(D_0)d^2}{Eh^2}. \tag{4.19}

A zero-mean quantity is zero almost surely, which covers that case as well. Freeze this stack.

A large independent noisy-sum atom. Let p=\prod_i\nu_{w_i}, and compare the actual conditional law of (X^s,\Gamma) to the independent law p\times q. The map (X^s,\Gamma)\mapsto(X^s,Z) is invertible, and -\log q(\Gamma)=F_0+c(\Gamma). Thus

C(\boldsymbol w)=D\bigl((X^s,\Gamma)\|p\times q\bigr) +\mathrm H(Z)-F_0. \tag{4.20}

Let Q_Z be the sum law under p\times q. Data processing gives

C(\boldsymbol w)+F_0 \ge D(\mathcal L(Z)\|Q_Z)+\mathrm H(Z) =\E[-\log Q_Z(Z)].

Some integer z_* in the support of the actual output therefore has Q_Z(z_*)\ge e^{-F_0-CE}=q(0)e^{-CE}. All Q_Z masses are positive since q is strictly positive. Divide the convolution formula by q(0) and use the characteristic formula:

\begin{split} e^{-CE} &\le Q_Z(z_*)/q(0)\\ &=\int\chi_{gz_*}(\theta) \prod_i\E_{\nu_{w_i}}\chi_{-gX}(\theta)\,d\omega(\theta)\\ &\le\int\prod_i \left|\E_{\nu_{w_i}}\chi_{gX}(\theta)\right|\,d\omega(\theta). \end{split} \tag{4.21}

The last step takes the absolute value of the preceding integral, whose value is a positive real mass.

One region and one coordinate. Put \phi_i(\theta)=\E_{\nu_{w_i}}\chi_{gX}(\theta) and A(\theta)=\sum_i(1-|\phi_i(\theta)|^2). The bound 0\le|\phi_i|\le1 implies \prod_i|\phi_i|\le e^{-A/2}. Choose a numerical C_R larger than twice the constant in Equation (4.21), with a fixed additional margin, and set \mathcal R=\{A\le C_RE\}. Since E\ge d/h^2\ge1, the contribution of \mathcal R^c to the product integral is at most half its lower bound. The product is at most one, so

\omega(\mathcal R)\ge e^{-CE},\qquad \sum_i\int(1-|\phi_i|^2)\,d\omega_{\mathcal R}\le CE, \tag{4.22}

where \omega_{\mathcal R} is normalized restriction to \mathcal R. Apply Markov’s inequality to a uniform choice of coordinate, using Equations (4.19) and (4.22). One coordinate satisfies both

\int(1-|\phi_i|^2)\,d\omega_{\mathcal R}\le CE/s,\qquad D(\nu_{w_i}\|U_S)\le\frac{C(D_0)d^2}{Eh^2}. \tag{4.23}

Only numerical factors are lost in the first bound. Write \mu=\nu_{w_i}.

A center and a subset valid for every weight. For independent X,\breve X\sim\mu,

1-|\phi_i(\theta)|^2 =\E[1-\operatorname{Re}\chi_{g(X-\breve X)}(\theta)].

Averaging over \omega_{\mathcal R} and then over \breve X selects x_0\in\operatorname{supp}\mu\subset S with

\sum_x\mu(x)\int(1-\operatorname{Re}\chi_{g(x-x_0)}) \,d\omega_{\mathcal R}\le CE/s.

Prune by Markov’s inequality to a subset S'\subset\operatorname{supp}\mu of \mu-mass at least 1/2, on which every x satisfies

\int(1-\operatorname{Re}\chi_{g(x-x_0)})\,d\omega_{\mathcal R}\le CE/s. \tag{4.24}

If the average defect is zero, its zero-defect set has full \mu-mass and can be used for S'. In all cases S' is nonempty.

Write t=\mu(S')\ge1/2 and v_0=|S'|/|S|. Binary data processing gives

D(\mu\|U_S)\ge t\log\frac t{v_0}+(1-t)\log\frac{1-t}{1-v_0} \ge t\log\frac1{v_0}-\log2.

Together with Equation (4.23), this implies

\log|S|-\log|S'|\le2D(\mu\|U_S)+2\log2 \le\frac{C(D_0)d^2}{Eh^2}.

The constant is absorbed because d^2/(Eh^2)\to\infty.

Finally, let v be any probability measure on this same S'. By Equation (4.24), its averaged phase defect on \mathcal R is at most CE/s. Jensen’s inequality and L\le s yield

\begin{split} &\int\exp\left[-L\sum_{x\in S'}v(x) (1-\operatorname{Re}\chi_{g(x-x_0)}(\theta))\right]\,d\omega(\theta)\\ &\quad\ge\omega(\mathcal R) \exp\left[-L\sum_{x\in S'}v(x) \int(1-\operatorname{Re}\chi_{g(x-x_0)})\,d\omega_{\mathcal R}\right]\\ &\quad\ge e^{-CE}e^{-CEL/s}\ge e^{-CE}. \end{split}

The region, center, and subset were fixed before v was chosen, establishing the conclusion for every probability on S'. The exponential constant is numerical: dependence on D_0 occurs in the entropy loss and in the sufficiently large threshold for d, whereas Equation (4.16) and all subsequent constant-factor selections have numerical bounds. This completes the proof.

5 Thermal labels and the spectral update

This section constructs the channels that will be attached to the spatial pieces. It also supplies two ways to pay for a change of channel: a deterministic increase in a metric ball, and an observable event that is unlikely under a comparison law. The choices of channels and of input membership sets will always be made without inspecting a sampled label.

5.1 Positive Fourier kernels and form order

Definition 5.1 (Thermal forms). A thermal form on the unit circle is a finite sum

A(y)=\sum_s\lambda_s(1-\operatorname{Re}\chi_s(y)), \qquad \lambda_s\ge0,\quad s\in\mathbb Z.

We write A\preceq A' if A'=A+D for a thermal form D. Define

\begin{gathered} q_A(u)=\int_{\mathbb T}e^{-A(y)}\overline{\chi_u(y)}\,dy, \qquad F(A)=-\log q_A(0),\qquad K_A=e^{F(A)-A},\\ \Delta_A(u)=\left(\int_{\mathbb T}|1-\chi_u(y)|^2K_A(y)\,dy\right)^{1/2}, \qquad B_A(r)=\{u\in\mathbb Z:\Delta_A(u)\le r\}. \end{gathered} \tag{5.1}

Thus K_A is a probability density. All circle integrals in this section use normalized measure of total mass one.

The covariance argument below belongs to the Griffiths–Ginibre positive-correlation method [10]. We include the circle calculation, together with the form-order consequences needed for the label updates.

Lemma 5.2 (Thermal positivity and monotonicity). For every thermal form, q_A is a nonnegative, even probability mass function on \mathbb Z. The coefficients \int\chi_uK_A=q_A(u)/q_A(0) are real and nonnegative. For all integers u,v,

\operatorname{Cov}_{K_A}(\operatorname{Re}\chi_u, \operatorname{Re}\chi_v)\ge0.

Consequently these characteristic coefficients increase in form order, \Delta_A(u) decreases in form order, and \Delta_A(u+v)\le\Delta_A(u)+\Delta_A(v). The function F increases in form order. For every thermal form D, the increment F(A+D)-F(A) decreases when A increases in form order.

If P=\sum_jP_j is a finite sum of thermal forms and P_{\rm lo}\preceq P\preceq P_{\rm hi}, then, for arbitrary shifts z_j and every integer u,

\left|\int\chi_u(y)\prod_j e^{-P_j(z_j-y)}\,dy\right| \le q_P(u) \le e^{-F(P_{\rm lo})}\int\chi_u K_{P_{\rm hi}}. \tag{5.2}

Proof. Expand each factor as

e^{-\lambda(1-\operatorname{Re}\chi_s)} =e^{-\lambda} \exp\bigl((\lambda/2)\chi_s\bigr) \exp\bigl((\lambda/2)\chi_{-s}\bigr).

Its Fourier coefficients are nonnegative, with sum one. The expansions are absolutely summable, as is their finite product. Evaluation at zero therefore gives \sum_uq_A(u)=1; evenness gives q_A(-u)=q_A(u).

For completeness, write c_s(y)=\cos(2\pi sy) and Z_A=\int\exp(\sum_s\lambda_sc_s). Express the covariance using two independent copies as half the expectation of the product of their differences. The map (x,z)\mapsto(x+z,x-z) pushes Haar probability on the two-dimensional torus to Haar probability. The two differences are -2\sin(2\pi ux)\sin(2\pi uz) and the analogous expression for v, whereas the unnormalized product density is \exp(2\sum_s\lambda_sc_s(x)c_s(z)). Expanding this last exponential gives the identity

\operatorname{Cov}_{K_A}(c_u,c_v) =\frac{2}{Z_A^2} \sum_{(n_s)\ge0}\left(\prod_s\frac{(2\lambda_s)^{n_s}}{n_s!}\right) \left(\int\sin(2\pi ux)\sin(2\pi vx) \prod_sc_s(x)^{n_s}\,dx\right)^2.

Absolute convergence justifies the expansion and proves nonnegativity. Differentiation under the integral yields

\frac{\partial}{\partial\lambda_v}\int c_uK_A =\operatorname{Cov}_{K_A}(c_u,c_v),\qquad \frac{\partial F}{\partial\lambda_u}=\int(1-c_u)K_A.

The latter derivatives are nonnegative and decrease in every nonnegative coefficient direction. Integrating these derivative statements proves the assertions about F and its increments. Since \Delta_A(u)^2=2(1-\int c_uK_A), the metric monotonicity follows as well. Subadditivity follows from 1-\chi_{u+v}=(1-\chi_u)+\chi_u(1-\chi_v) and the triangle inequality in L^2(K_A).

Translation multiplies each Fourier coefficient of e^{-P_j} by a unit complex number. Absolute values in the convolution of these coefficients are therefore bounded by the corresponding convolution of the centered coefficients, which is q_P(u). Finally q_P(u)=e^{-F(P)}\int\chi_uK_P, and both required comparisons follow from the monotonicities already proved.

For a real interval p of length at most one, set

R_p=|p|^{-1},\qquad S_p(y)=R_p^2(1-\operatorname{Re}\chi_1(y)).

The interval may be viewed on the circle when integrating. We will use the following bounds, also with any larger spatial scaling in place of R_p:

\begin{gathered} F(S_p)=\log R_p+O(1),\qquad \int(1-\operatorname{Re}\chi_1)K_{A+S_p}\le C R_p^{-2},\\ q_{S_p}(u)\ge cR_p^{-1}\quad (|u|\le cR_p),\qquad \Delta_{A+S_p}(u)\le C|u|/R_p. \end{gathered} \tag{5.3}

Indeed, if \|y\| denotes distance to the nearest integer, then 8\|y\|^2\le1-\cos(2\pi y)\le2\pi^2\|y\|^2. Gaussian comparison gives \int e^{-S_p}\asymp R_p^{-1} and \int(1-\operatorname{Re}\chi_1)K_{S_p}\le CR_p^{-2}. Monotonicity gives the second bound for A+S_p. The elementary inequality |1-\chi_u|\le|u||1-\chi_1| gives the last bound and, for sufficiently small |u|/R_p, gives \int\chi_uK_{S_p}\ge1/2, proving the remaining lower bound. In particular \int S_pK_{A+S_p}\le C.

5.2 Channels, scales, and the cost of tightening a form

Take b_0 dyadic with h^{-16}/2<b_0\le h^{-16}. The computation levels are the dyadic b from b_0/8 to 1, and

k_b=A_0\frac d h(1+\log(1/b))^2,\qquad r_{\rm sep}=h^{-30},\qquad \sigma=h^{-40}. \tag{5.4}

The constants will be fixed in the order stated below. Every estimate holds for sufficiently large d, uniformly in the spatial depth and in any frequency cutoff.

A label of type A is Z=Y+\theta modulo one, where the error \theta has density K_A. Initially A=0, so a dummy independent uniform label may be drawn. Only the current main channel is used in the spectra and transforms. A successful update adds a spatial form and a metric form, A^+=A+S_p+D_i. The errors used for different label draws are independent given Y. Their parameters are selected from spatial prefixes and earlier input membership events. In particular neither a sampled old label nor a sampled fresh label selects a membership set, channel, or test. Conditional on Y and its discrete parameter path, the errors thus retain their prescribed independent laws.

We will need to tighten the forms of past labels at small free-energy cost and to keep their own-error energies \int A K_A small. Both estimates follow by choosing each metric addition on a sufficiently flat interval of free energy. We require, with A_{\rm base} the form immediately before the addition and including its new spatial term,

F(A_{\rm base}+h^{140}D_i)-F(A_{\rm base}+D_i/2) \le d h^{-400}. \tag{5.5}

Before a requested update assume that fewer than C_0h main additions have occurred and that

F(A)\le\log R_*+C_*d. \tag{5.6}

Here p_* is the interval of the most recent addition, R_*=R_{p_*}, and the current p is contained in p_*. Before the first addition p_* is the top interval, whose length lies in [1/2,1]. The procedure in Section 6 will establish these hypotheses inductively.

Fix a sufficiently large numerical C_s and put \widetilde A=h^{130}A+C_sS_p. Then

\begin{split} F(\widetilde A)-F(A+S_p) &\le C+CC_0h\bigl(dh^{-400}+C\log h\bigr),\\ \int A K_A&\le CC_0h\bigl(dh^{-400}+1\bigr). \end{split} \tag{5.7}

To prove the first line, enlarge the constituent forms one at a time. For an old metric term D_i, decreasing increments bound its cost of tightening by the left side of Equation (5.5): the other terms already contain its original A_{\rm base}. Tightening an old spatial term costs at most F(h^{130}S_{p_i})-F(S_{p_i})\le C(1+\log h) by Equation (5.3). Tightening the current newly inserted S_p to C_sS_p costs O(1). For the energy bound, concavity along the D_i direction and Equation (5.5) give

\int D_iK_A\le 2\{F(A_{{\rm base},i}+D_i)-F(A_{{\rm base},i}+D_i/2)\} \le2dh^{-400}.

Each spatial constituent contributes O(1) by Equation (5.3) and monotonicity. This energy argument uses the count bound and Equation (5.5), but does not use Equation (5.6); it remains valid, with the same type of bound, just after a final allowed addition that breaches a later stopping cap.

The spatial increment itself obeys

F(A+S_p)\le F(A)+\log(R_p/R_*)+C. \tag{5.8}

After an addition, A contains S_{p_*}, so q_A(u)\ge ce^{-F(A)} for |u|\le cR_*. Convolution with q_{S_p} over those integers gives q_{A+S_p}(0)\ge ce^{-F(A)}R_*/R_p. Initially use the term u=0 alone and 1\le R_*\le2.

5.3 From a separated spectrum to a finite compression problem

Consider a subdensity a at an interval p, with

0\le a\le e^d,\qquad e^{-d}\le m=\mathbb E_pa\le2\alpha, \qquad \alpha>0.

Here \mathbb E_p is uniform spatial average, and \mathcal P is the law of Y with density a/m relative to that average. For either prospective child p' define

\begin{split} M_{p'}(z)&=\mathbb E_{p'}[aK_A(z-Y)]/\alpha,\\ f_{p'}^u(z)&=\mathbb E_{p'}[a\chi_u(Y)K_A(z-Y)]/\alpha,\\ G_{p'}^u&=\int |f_{p'}^u(z)|^2/M_{p'}(z)\,dz. \end{split} \tag{5.9}

The integrand is zero at zero mass. An integer is high at b when some G_{p'}^u\ge b^2. A test requests an update if there is a collection of \lceil e^{k_b}\rceil high integers whose pairwise differences avoid

B_A(r_{\rm sep})+\{j\in\mathbb Z:|j|\le\sigma R_p\}. \tag{5.10}

One offending level and one collection are used at a time.

Keep at least half this collection having a common high child, and snap each retained integer to a nearest multiple gx, where

g=\max(1,\lfloor\sigma R_p/20\rfloor).

The snapping error is at most \sigma R_p/20, and is zero when g=1. The snapped indices are distinct, since otherwise the original difference would violate Equation (5.10). Denote by S the resulting set of integer indices x. For each x\in S choose a measurable complex function v_x, with |v_x|\le1, so that

\operatorname{Re}\mathbb E_{\mathcal P} [\mathbf1_{\rm child}v_x(Z)\chi_{gx}(Y)]\ge\beta, \qquad \beta=b^2/100. \tag{5.11}

The expectation here includes the type-A channel. To see the bias, |f_{p'}^u|\le M_{p'} implies \int|f_{p'}^u|\ge G_{p'}^u. Choosing the polarizing sign gives bias (\alpha/(2m))\int|f_{p'}^u|\ge b^2/4 before snapping. A constant phase rotation at the center of p makes the error of replacing u by gx at most C\sigma. This is smaller than the required slack because \sigma=o(b^2) uniformly over the computation levels.

We apply Lemma 3.1 with

d\omega=K_{\widetilde A}\,dy,\qquad q(x)=\frac{q_{\widetilde A}(gx)}{N_g},\qquad N_g=\sum_{x\in\mathbb Z}q_{\widetilde A}(gx).

The spatial frequency-one term makes every Fourier coefficient strictly positive: in its exponential expansion, every integer frequency occurs with positive coefficient. Therefore q is strictly positive. The grid identity and Gaussian comparison give

N_g=\frac1g\sum_{l=0}^{g-1}e^{-\widetilde A(l/g)},\qquad \frac1g\le N_g\le\frac Cg. \tag{5.12}

Indeed g\le R_p, and the summands are bounded by \exp(-cR_p^2\|l/g\|^2). Moreover R_p/g\le C/\sigma, also when g=1 by its definition. Thus Equations (5.6), (5.7) and (5.8) imply

F_0=-\log q(0)=F(\widetilde A)-\log g+O(1)\le D_0d

for a fixed sufficiently large D_0 depending only on the fixed caps. Evenness gives precisely q(x)/q(0)=\int\chi_{gx}\,d\omega.

We verify both remaining uniform hypotheses of compression.

Sumset size. We use the dissociated-spectrum method underlying Chang’s lemma [5, Section 3], in its entropy and positive-product form [15, Section 3.2]. The following calculation includes the side information carried by the label: the comparison law makes that information independent of the uniform circle coordinate.

The law of gY modulo one has density at most e^{2d}(1+R_p/g): its preimage count in p is at most g/R_p+1, and the density of Y relative to circle length is at most R_pe^{2d}. Also

\mathrm I(Y;Z)\le F(A)-\log R_*+C. \tag{5.13}

Compare each K_A(z-Y) to the fixed density K_{S_{p_*}}(z-\operatorname{center}(p_*)). The average relative entropy is

F(A)-F(S_{p_*})-\int AK_A +\mathbb E S_{p_*}(\theta+Y-\operatorname{center}(p_*)).

The last expectation is O(1): use S(v+w)\le2S(v)+2S(w), the spatial energy bound, and |Y-\operatorname{center}(p_*)|\le|p_*|/2. In the initial case A=0 and R_*\le2, the same bound is immediate. Dropping the nonpositive term proves Equation (5.13).

Let Q denote the side information consisting of Z and the child bit. Data processing and its one-bit entropy imply

D\bigl((gY,Q)\,\|\,\text{uniform}\otimes\mathcal L(Q)\bigr) =D(gY\|\text{uniform})+\mathrm I(gY;Q)\le C(C_*)d.

If \Lambda\subset S is dissociated, meaning that no nonzero coefficient choice from \{0,\pm1\} sums to zero, form the positive product

\prod_{x\in\Lambda} \left(1+c\beta\operatorname{Re} [\mathbf1_{\rm child}v_x(Z)\chi_x(gY)]\right).

For sufficiently small numerical c, the factors lie in [1/2,3/2]. Under the reference law the first coordinate is uniform and independent of Q; expansion and dissociation make the expectation of the product exactly one. Relative entropy is at least the true expected logarithm of this product. The inequality \log(1+s)\ge s-Cs^2, together with Equation (5.11), bounds the expected logarithm of each factor below by c'\beta^2. Hence |\Lambda|\le C(C_*)d/\beta^2. A maximal dissociated subset spans S by coefficients 0,\pm1: otherwise adjoining an unspanned element preserves dissociation. Consequently

|nS|\le(2n+1)^{|\Lambda|},\qquad \log|nS|\le D_0d h^{P_0}\log d\quad(1\le n\le d^2)

for fixed D_0,P_0, because \beta\ge c h^{-32}.

The uniform moment cap. For any probability \nu on S set

\begin{split} E_T(\nu)&=-\log\mathbb E_\nu H_q\left(\sum_{i=1}^TX_i-\sum_{i=1}^TX_i'\right),\\ e^{-E_T(\nu)}&= \int\left|\sum_x\nu(x)\chi_{gx}(y)\right|^{2T}K_{\widetilde A}(y)\,dy, \end{split} \tag{5.14}

where the 2T letters are independent with common law \nu. Then

0\le E_T(\nu)\le D_0d\qquad\hbox{for every such }\nu. \tag{5.15}

Here is a direct proof of the upper bound. Restrict Z to the domain \Omega at circle distance at most h^{40}/R_* from p. The channel spatial variance bounds the discarded probability by Ch^{-80} (the assertion is trivial if \Omega is the full circle). This loses o(\beta) in Equation (5.11). Summing the biases with weights \nu and applying Jensen on a subprobability measure gives

\mathbb E_{\mathcal P} \left[\mathbf1_{\rm child}\mathbf1_{\{Z\in\Omega\}} \left|\sum_x\nu(x)v_x(Z)\chi_{gx}(Y)\right|^{2T}\right] \ge(\beta/2)^{2T}.

For an upper bound discard the child indicator. At fixed z, first majorize the nonnegative moment integrand using the spatial density bound R_pe^{2d} and

\mathbf1_p(y)\le \exp\bigl(C-S_p(y-\operatorname{center}(p))\bigr).

Only after this majorization expand the even power and the shifted thermal weights. By Equation (5.2), the Fourier coefficients of those weights are bounded absolutely by

Cq_{A+S_p}(u)\le Ce^{-F(A+S_p)}\int\chi_u K_{\widetilde A}.

The coefficients v_x(z) have modulus at most one. The absolute coefficient sum in the even-power expansion therefore averages this envelope at u=g(\sum X_i-\sum X_i') over the independent letters in Equation (5.14). Using |\Omega|\le Ch^{40}/R_*, the joint moment is at most

C R_p e^{2d+F(A)}\frac{h^{40}}{R_*} e^{-F(A+S_p)}e^{-E_T(\nu)}.

Now F(A+S_p)\ge\log R_p-C, F(A)\le\log R_*+C_*d, and T\log(2/\beta)=o(d). Comparison with the lower moment proves Equation (5.15) after increasing the fixed D_0. This argument uses a density bound only in a nonnegative integral; it does not assert a Fourier coefficient bound for a.

5.4 The iterative update and the deterministic alternative

Lemma 5.3 (Spectral update). Fix C_0,C_*\ge1. Let the current form A arise from the initial form zero by fewer than C_0h additions of spatial and metric forms, each metric satisfying Equation (5.5). Let p\subset p_*, where p_* is the latest addition interval, or the top interval of length in [1/2,1] if there has been no addition, and assume F(A)\le\log R_{p_*}+C_*d. On p suppose 0\le a\le e^d and e^{-d}\le m=\mathbb E_pa\le2\alpha, with \alpha>0. Suppose there is a request at a computation level b in the sense of Equation (5.10), using the prospective child moments in Equation (5.9) and the type-A channel. There are numerical constants C_1,C_2 and a sufficiently small numerical \varepsilon>0 with the following property. After the caps are fixed, choose fixed sufficiently large D_0,P_0, then A_0, and finally take d sufficiently large. A finite deterministic construction supplies a parameter

d/h^2\le t\le D_0d,\qquad t\ge\varepsilon k_b,

a thermal metric form D_i=\lambda B, and an input-only success set. On success set A^+=A+S_p+D_i. The new metric satisfies Equation (5.5) and

\begin{gathered} A^+=A+S_p+D_i,\qquad \int B K_{A^+}\le C_1t/\lambda,\\ \lambda\ge dh^{200},\qquad h^{140}\lambda\le L, \end{gathered} \tag{5.16}

and F(A^+)-F(A+S_p)\le C_1t. There are two alternatives:

  1. (i)

    A regular test has success probability \delta\ge\beta/2 under the current input law, where \beta=b^2/100. It also supplies a probability \nu on a set S' and a center x_0\in S for which E_T(\nu)\ge t/4.

  2. (ii)

    A cheap update succeeds deterministically, has t=\varepsilon k_b, and uses the uniform probability on S'.

In both alternatives \log|S'|\ge0.8k_b and B=\sum_x\nu(x)(1-\operatorname{Re}\chi_{g(x-x_0)}). Track a radius starting at r_{\rm sep}/4 and increasing by \eta=r_{\rm sep}/(10C_0h) at every main addition. As long as the count cap applies, the gap F(A)-\log|B_A(r)| starts at zero, stays bounded below by a numerical negative constant, increases by at most C_2t at a regular success, and decreases by at least k_b/2 at a cheap success. On failure there is no main addition and no label draw. All these conclusions are independent of spatial depth and frequency cutoff.

Proof. Selecting the subset and its weight. On S and then on its successive restrictions set

t=\max\left(\varepsilon k_b, \max_{\nu\text{ supported on the current set}}E_T(\nu) \right).

The maximum exists: the simplex is compact, its moment integral is continuous and strictly positive, and therefore its negative logarithm is continuous. For large d, the floor is between d/h^2 and D_0d; the upper bound follows also for the maximum from Equation (5.15). Apply Lemma 3.1 with E=t to obtain a subset S' and center x_0. Its conclusion holds for every probability on this fixed S', so the weight may now be chosen by the following test. If the maximum of E_T on S' is at least t/4, keep a maximizing witness and terminate with a regular test. If that maximum is smaller than t/4 and t=\varepsilon k_b, take uniform \nu on S' and terminate with a cheap update. Otherwise replace the current set by S' and repeat. Every nonterminal iteration reduces t by a factor of at least four, except possibly for the last transition to the fixed floor. Thus the process terminates in finitely many iterations, and the sum of the logarithmic cardinality losses is at most

C(D_0)\frac{d^2}{\varepsilon k_bh^2}.

The initial child restriction lost at most \log2 in logarithmic cardinality. Since k_b\ge A_0d/h, choosing A_0 sufficiently large in terms of D_0,P_0,\varepsilon makes the final loss less than 0.2k_b, proving \log|S'|\ge0.8k_b.

Choosing a flat coefficient range. For the last center and chosen probability, the compression conclusion is

F(\widetilde A+LB)-F(\widetilde A)\le Ct.

Put A_{\rm base}=A+S_p. Monotonicity of F and A_{\rm base}\preceq\widetilde A give

F(A_{\rm base}+LB)-F(A_{\rm base}) \le \bigl[F(\widetilde A+LB)-F(\widetilde A)\bigr] +\bigl[F(\widetilde A)-F(A_{\rm base})\bigr].

The first bracket is controlled by compression and the second by the tightening bound in Equation (5.7). Hence

F(A_{\rm base}+LB)-F(A_{\rm base})\le C_1t. \tag{5.17}

The constant C_1 is numerical: the volume constant of compression is numerical, and the cap-dependent error in Equation (5.7) is o(t) after the fixed caps and A_0 have been selected.

To find a flat range, set

\lambda_j=dh^{200}(2h^{140})^j, \qquad N=\left\lfloor \frac{\log(2L/(dh^{200}))}{\log(2h^{140})}\right\rfloor.

For 0\le j<N, the intervals [\lambda_j/2,h^{140}\lambda_j] are adjacent and lie in [dh^{200}/2,L]. Since \log(L/d)=w=e^{h/4}, N\asymp e^{h/4}/\log h>C_1D_0h^{400} for sufficiently large d. The sum of the free-energy increments on these intervals is bounded by Equation (5.17), hence by C_1D_0d. At least one increment is at most dh^{-400}. Take its \lambda and D_i=\lambda B; this is Equation (5.5). Concavity along the B direction gives

\int BK_{A_{\rm base}+\lambda B} \le\frac{F(A_{\rm base}+\lambda B)-F(A_{\rm base})}{\lambda} \le C_1t/\lambda,

which proves Equation (5.16) and the metric increment bound.

The input-only success set. For a regular test define the success set by

\left|\sum_x\nu(x)\chi_{gx}(Y) \mathbb E[v_x(Z)\mid Y]\right|\ge\beta/2. \tag{5.18}

The expectation uses the specified type-A channel, rather than the sampled old observation. The expression has modulus at most one, and Equation (5.11), averaged with \nu, gives its expected modulus at least \beta. Thus its probability of exceeding \beta/2 is at least \beta/2. In the cheap alternative take the entire input as success. The data a,p,\alpha,A determine all these sets and choices.

The change in the gap potential. It remains to prove the gap assertions. Since r\le r_{\rm sep}/2 through the count cap, summability of q_A gives

|B_A(r)|e^{-F(A)}(1-r^2/2)\le1, \qquad F(A)-\log|B_A(r)|\ge\log(1-r^2/2)\ge-C.

The balls are finite by the same inequality. Initially A=0 and B_0(r)=\{0\}, so the gap is zero. The spatial addition raises F, but also allows short integer translates of the old ball inside a slightly larger new ball. We compare these changes using the same overlap multiplicity. Write

U=B_A(r),\qquad D_p=\mathbb Z\cap[-\sigma R_p/10,\sigma R_p/10],\qquad l=|(U-U)\cap(D_p-D_p)|.

Here l\ge1. On this intersection \Delta_A\le2r and |u|\le\sigma R_p/5, so convolution and Equation (5.3) imply

F(A+S_p)-F(A)\le\log R_p-\log l+C, \qquad |U+D_p|\ge |U||D_p|/l. \tag{5.19}

For the second inequality, each representation of a fixed element of U+D_p differs from one fixed representation by an element of the intersection defining l, so its multiplicity is at most l.

After either addition, U+D_p lies in B_{A^+}(r+\eta/2), since monotonicity preserves the old radius and \Delta_{A^+}(j)\le C|j|/R_p\le C\sigma=o(\eta) for j\in D_p. On a cheap success, uniformity of \nu and Equation (5.16) give

\frac1{|S'|}\sum_{x\in S'}\Delta_{A^+}(g(x-x_0))^2 \le 2C_1t/\lambda.

At least half the x therefore have \Delta_{A^+}(g(x-x_0))\le\sqrt{4C_1t/\lambda} \le C(D_0)h^{-100}=o(\eta). Their translates g(x-x_0)+U+D_p lie in the new tracked ball B_{A^+}(r+\eta). They are disjoint. A collision would put g(x-x') in U-U+D_p-D_p; after restoring the two snapping errors, the difference of the original frequencies would lie in

B_A(2r)+\{j\in\mathbb Z:|j|\le3\sigma R_p/10\},

contrary to Equation (5.10).

Since R_p/|D_p|\le C/\sigma, including the case D_p=\{0\}, Equation (5.19) cancels the same \log l in the free-energy and cardinality costs. On a regular success the gap increase is at most C_1t+O(\log h)\le C_2t for a numerical C_2. On a cheap success it is at most

C_1\varepsilon k_b-\log(|S'|/2)+O(\log h) \le C_1\varepsilon k_b-0.8k_b+O(\log h)\le-k_b/2,

on choosing \varepsilon sufficiently small numerically and then taking d sufficiently large. These choices precede the cap choices; only the size-loss bound and the small-error thresholds require the subsequent choices of A_0 and d.

5.5 Evidence under the comparison law and under truth

At every success draw a new main label V of type A^+. On a regular success also draw N_h=\lceil h^{120}\rceil fresh, independent type-A clones Z'. They serve only the evidence calculation. These draws are made even when this very success will breach a cap and transfer the piece out of live processing. A failure causes neither a new main draw nor an evidence draw.

The comparison experiment starts with Y uniform on the top interval. Spatial bits continue to observe Y, and all label draws use the same channels prescribed on their recorded branches. At a membership test, however, success is an independent coin with probabilities \delta,1-\delta computed from the current true input law before conditioning on sampled labels. For a fixed recorded spatial and membership branch, the computed state is fixed. Its membership coin probabilities therefore contribute no y-dependent likelihood. Given the full past observations at a test in p, the comparison posterior has density on p proportional to

\exp\left(-\sum_jA_j(z_j-y)\right). \tag{5.20}

The sum includes past main labels and past evidence clones. Each A_j\preceq A, and their number is at most 1+C_0h(1+N_h). This statement conditions on the observed branch, not on the actual truth-membership event within the reference law.

Set s_d=dh^{-100}. A past is admissible when some y_0\in\overline p satisfies

\sum_jA_j(z_j-y_0)\le s_d.

This is an observable condition: the continuous function in the display attains its infimum on the compact interval closure.

Proposition 5.4 (Regular evidence). Consider any finite procedure starting with the type-zero label, making at most \lfloor C_0h\rfloor successful main additions and using the updates of Lemma 5.3 with its hypotheses at every attempt. Require all discrete choices to depend only on spatial prefixes and input memberships, with independent label errors given Y and its parameter path, and include all the success draws specified above. Then an admissible regular attempt has the following conditional reference bound, given its past and its success coin. The observable event

\begin{gathered} \text{for some fresh clone }Z',\qquad A(Z'-Z)\le4s_d,\\ \left|\sum_x\nu(x)v_x(Z')\chi_{gx}(V)\right|\ge\beta/4. \end{gathered} \tag{5.21}

has probability at most e^{-t/8}, where Z is the old main label. Under the true input law, all attempted pasts are admissible and Equation (5.21) occurs at every regular success simultaneously, except on a set of probability o(h^{-20}). This includes the label draws at a final success breaching a cap. The bounds are uniform in depth, cutoff, and the deterministic choices made at the states.

Proof. First bound the posterior normalizing integral from below at an admissible past. The phase triangle inequality gives A_j(v+u)\le2A_j(v)+2A_j(u), and the observation count gives 2\sum_jA_j\preceq h^{130}A for sufficiently large d. Thus

\int_p\exp\left(-\sum_jA_j(z_j-y)\right)dy \ge c\exp(-2s_d-F(\widetilde A)). \tag{5.22}

Indeed in one direction from y_0 a segment of length |p|/2 remains inside \overline p. By symmetry and the spatial variance estimate, the density K_{\widetilde A}(y-y_0) assigns this segment at least a numerical positive mass if C_s is large enough: the circular second-moment bound puts a fixed positive mass within distance |p|/2 of zero, and symmetry gives half of that mass in either direction. On this segment the integrand is at least e^{-2s_d}\exp(-h^{130}A(y-y_0)), which is at least e^{-2s_d}\exp(-\widetilde A(y-y_0)).

Fix one fresh clone and integrate over its value z' in the consistency set \{z':A(z'-Z)\le4s_d\}. Also use the variable \theta^+=V-Y, whose independent conditional density is K_{A^+}. At fixed z',\theta^+ expand the 2T-moment of

\left|\sum_x\nu(x)v_x(z')\chi_{gx}(Y+\theta^+)\right|

in the posterior integral. Its weights include the extra clone factor e^{-A(z'-y)} and the spatial majorant used above. The centered summed form is

P=\sum_jA_j+A+S_p, \qquad A+S_p\preceq P\preceq\widetilde A.

The upper comparison follows from the observation count and large d. By Equation (5.2), the y integral before the clone normalization and the posterior division is at most

Ce^{-F(A+S_p)}e^{-E_T(\nu)} \le Ce^{-F(A+S_p)}e^{-t/4}.

Here the phase factors involving \theta^+ have modulus one, and |v_x(z')|\le1, so the same coefficientwise moment bound applies. The consistency set has length at most e^{4s_d-F(A)}, since e^{-A}\ge e^{-4s_d} there. This cancels the clone normalization e^{F(A)} up to e^{4s_d}. The density in \theta^+ integrates to one. Dividing by Equation (5.22), applying the moment tail bound at \beta/4, and summing over the clones gives

C N_h\exp\left( 6s_d+F(\widetilde A)-F(A+S_p) +2T\log(4/\beta)-t/4\right)\le e^{-t/8}.

Indeed Equation (5.7), t\ge\varepsilon k_b, and 2T\log(4/\beta)=o(t) absorb every positive term. The reference success coin is independent of Y given the past and so has not altered its posterior.

For the true-law assertion, condition on Y and its entire discrete parameter path. At a regular success the conditional mean in Equation (5.18) has modulus at least \beta/2. The empirical mean of the N_h clone polynomials evaluated at Y approximates it to error \beta/16 except with probability C/(\beta^2N_h), by the second-moment bound for independent bounded complex variables. Replacing Y by V, and removing the common phase \chi_{gx_0}(V-Y), changes each of these polynomials by at most

\sum_x\nu(x)|1-\chi_{g(x-x_0)}(V-Y)|.

Its conditional expectation is at most \sqrt{2\int BK_{A^+}}\le C\sqrt{t/\lambda}=O(h^{-100}) by Equation (5.16). Markov’s inequality makes this error at most \beta/16 except with probability C\sqrt{t/\lambda}/\beta=o(h^{-30}), since \beta\ge ch^{-32}. Outside these errors the clone average after replacement has modulus at least 3\beta/8, so at least one clone satisfies the large-polynomial condition in Equation (5.21). The empirical-mean error probability is O(h^{-56}), also o(h^{-30}).

Finally sum the own-error energies of every observation on the capped true path, each measured in its own form. Conditional independence and the input-only selection rule let us apply the energy estimate in Equation (5.7) at each draw, including final draws. There are at most O(h^{121}) observations, each of conditional expected energy O(h(dh^{-400}+1)). Their total expected energy is therefore

O(dh^{-278}+h^{122})=o(s_dh^{-30}).

With exception probability o(h^{-30}), their total energy is at most s_d. This single event makes every past admissible by taking y_0=Y. It also gives every consistency condition, since the old main and its clone both have type A and

A(Z'-Z)\le2A(Z'-Y)+2A(Z-Y)\le4s_d.

There are only O(h) regular successes. Unioning their conditional large-polynomial failure bounds and adding the total-energy exception proves the simultaneous o(h^{-20}) assertion. The bound is over successes rather than over all attempts; past admissibility at all attempts followed from the single total-energy event.

The numerical constants C_1,C_2 and the permitted small \varepsilon have now been fixed independently of the later caps. The next section chooses those caps, then fixes D_0,P_0,A_0, and finally increases the threshold for d. The evidence clones enter only the likelihood estimate for regular charges. The main labels also enter the conditional-information bounds for the spatial budgets. The realized values of both kinds of labels do not control the spatial splitting procedure.

6 The finite dyadic procedure and its budgets

We seek a bound for modulated sums whose coefficients are formed from one input bit and emitted with a later output bit. The bound must be uniform in the tree depth and in the number of frequencies. We state it first; the construction in this section supplies the pieces and budgets used in its proof in Section 9.

Proposition 6.1 (Modulated delayed Haar estimate). There are absolute constants C and d_0 with the following property. Let d\ge d_0, put h=\log\log d, and let I\subset\R be a half-open interval with 1/2\le |I|\le1. Let \rho:I\to[0,\infty) be measurable, with \rho\le e^d almost everywhere and

m_{\rm top}=\frac1{|I|}\int_I\rho\le2.

Give I its normalized uniform probability \mathbf U. Write \mathcal F_j for its depth-j dyadic sigma field, p_j(x) for the depth-j interval containing x, and r_{j+1} for the sign that is +1 on the left child and -1 on the right child of each such interval. Let J\ge1 and Q\ge0 be integers. For 0\le j<J, let \epsilon_j:I\to\C be \mathcal F_{j+1}-measurable with |\epsilon_j|\le1, and use the same \epsilon_j for all modulations. The bits through depth J+1 are retained, so the last output sign is r_{J+1}.

Set

\tau(x)=\min\bigl(\{0\le j<J:\ \E_{p_j(x)}\rho>e^{4h}\}\cup\{J\}\bigr), \qquad \E_p v=\frac1{|p|}\int_p v(y)\,dy,

and, with \chi_u(y)=e^{2\pi iuy}, define

\mathcal H_{I,J}^{u}\rho(x) =\sum_{0\le j<\tau(x)} \E_{p_j(x)}[\rho(Y)\chi_u(Y)r_{j+1}(Y)]\, \epsilon_j(x)r_{j+2}(x).

Then

\int_I\max_{\substack{u\in\Z\\|u|\le Q}} |\mathcal H_{I,J}^{u}\rho(x)|\,\frac{dx}{|I|} \le Cd. \tag{6.1}

The constants are independent of I,\rho,J,Q, and the multipliers. Consequently the same estimate holds with \sup_{u\in\Z} in place of the finite maximum.

6.1 Input pieces and the order of processing

Fix the data of Proposition 6.1. At an output point x in a depth-j interval p, the summand in that proposition is

\mathbb E_p[\rho(Y)\chi_u(Y)r_{j+1}(Y)]\, \epsilon_j(x)r_{j+2}(x),\qquad u\in\mathbb Z. \tag{6.2}

The input variable Y and output point x play different roles. The coefficient is known before the output bit r_{j+1}(x), whereas its multiplier \epsilon_j(x) need only be known before r_{j+2}(x). The bits retained through depth J+1 supply this delayed sign even for the last contributing jump. The same multipliers are used for every piece and every label coordinate.

Every piece and frequency uses the common stopping rule \tau of Proposition 6.1, defined from the original density \rho. No processing is needed at a node where this rule has already stopped the output path. At all other nodes, the processing below is completed before the jump.

Discard the zero input. If 0<m_{\rm top}<1/h, put the whole input into static handling at the top. Otherwise define the true input law

d\mathbf P(Y)=\frac{\rho(Y)}{m_{\rm top}}\,d\mathbf U(Y).

The word “true” will always refer to this input law, or to its conditional laws, and not to uniform output sampling.

A processing state at a spatial interval p consists of a subdensity a, its mean m=\mathbb E_p a, and its recorded spatial and membership history. The subdensity is the restriction of \rho by the membership sets on that history. A zero-mass state is discarded. A live state also carries a main thermal form A, a count n of successful main additions, and the charge

\mathcal Q=\sum_{\text{regular successes so far}}t.

Initially A=0, n=0, and \mathcal Q=0. These three quantities carry along the entire branch and are never reset at a new epoch or plateau.

Fix constants R_0,C_0,C_* in the order specified below. The processing rules are as follows.

  1. 1.

    At every incoming live state, transfer its entire current subdensity to static handling if

    m<e^{-d},\qquad \mathcal Q>R_0d, \qquad\text{or}\qquad n\ge\lfloor C_0h\rfloor.

    This transfer ends its live processing. Its density is subsequently only restricted to spatial descendants; it is never split by another membership test.

  2. 2.

    The initial live piece starts an epoch. Every successful main addition starts a new epoch on its success piece unless that piece transfers immediately. Within an epoch, a plateau starts with normalization \alpha=m at entry. At a subsequent incoming state it ends and a new plateau starts if m>2\alpha. The new normalization is the incoming mean. No new main label is drawn at such a reset.

  3. 3.

    At a live state that has passed the transfer and reset rules, apply the spectral test of Lemma 5.3 to (a,p,\alpha,A), using the prospective child moments. Thus e^{-d}\le m\le2\alpha. The tests range over all computation levels b\in[b_0/8,1] that are powers of two; stabilization means that at each such level there is no separated high collection of size \lceil e^{k_b}\rceil. If a test requests an update, its success set E is a measurable set of inputs, determined from this state. Replace the current density by the two subdensities

    a\mathbf1_E\quad\hbox{and}\quad a\mathbf1_{E^c}.

    Only the success piece receives the new form A^+ and the increment of the count; a regular success also adds t to \mathcal Q. Its mean is \delta m, where \delta is its conditional success probability. The failure piece, if it has positive mass, has mean (1-\delta)m and retains the old epoch, plateau, form, and charge. For a cheap update \delta=1.

  4. 4.

    Process both resulting pieces at the same spatial interval p, enforcing transfer and reset rules anew. Continue until every remaining live piece at p passes all spectral tests. Then, and only then, take the spatial jump for each surviving live piece and each static piece, and proceed to its children subject to the common stop.

Static pieces also have plateaus: start with normalization equal to their mean at transfer (or at initial static entry), and restart whenever the incoming mean exceeds twice that normalization. They undergo no spectral tests or membership splits. Their estimates will use only the scalar construction.

Every test, witness, and membership set is chosen from the spatial prefix and the preceding membership outcomes. In particular, none of these choices uses a realized noisy observation or the next spatial bit. The label experiments are attached to these already specified discrete branches. Start with a type-0 main observation, independent of the input. Draw a fresh main observation of type A^+ on every success. On a regular success also draw the old-type evidence clones specified in Proposition 5.4. Make these draws even when that success immediately causes a transfer. Moving to a smaller interval, following a failure, or resetting \alpha does not redraw the main observation. Given the input and its discrete branch, every draw uses fresh independent error with its specified channel law.

At a jumping node the density belonging to a plateau is its final post-processing survivor. Its prospective child densities are the restrictions of that survivor, before any processing at the children. Removed mass has instead been assigned to another plateau or to static handling. Therefore the densities of all jumping pieces at p sum exactly to \rho|_p, and their normalized contributions, after multiplication by their respective normalizations, sum exactly to Equation (6.2). A true input follows just one of these pieces; uniform output estimates must sum all of their starting weights.

6.2 Statement of the budgets

For a live plateau e, write p_e for its root and \alpha_e for its initial mean, and put \mathbf U_e=\mathbf U(\,\cdot\mid p_e). At a final live jumping state p belonging to e, use its post density a and its fixed main channel to define

M_\pm(z)=\frac{\mathbb E_{p_\pm}[a(Y)K_A(z-Y)]}{\alpha_e},\qquad M=\frac{M_++M_-}{2},\qquad D=\frac{M_+-M_-}{2},\qquad J_p=\int\frac{|D|^2}{M}\,dz.

The piece and plateau are implicit in J_p; zero divided by zero is taken to be zero. Notice that \int M=m_p/\alpha_e\le2, where m_p=\mathbb E_p a, whereas \int M_\pm\le4. In particular J_p\le2.

Proposition 6.2 (Live budgets). There are fixed absolute choices of the constants in the procedure and in Lemma 5.3 such that, for all sufficiently large d, the procedure is finite in every finite spatial tree and preserves the exact splitting identity above. When 1/h\le m_{\rm top}\le2, it has the following bounds, with constants independent of J, I, the input, and any frequency cutoff.

  1. (i)

    Before every requested live update the count and form hypotheses of Lemma 5.3 hold. On every input branch, including its final success if a cap is crossed,

    \sum_{\rm successes}k_b+\sum_{\rm successes}t +\sum_{\rm successes}\bigl(F(A^+)-F(A+S_p)\bigr)\le Cd, \qquad n\le C h.

    The count flag cannot cause transfer unless the charge flag has already been crossed at that success.

  2. (ii)

    Let B_{\rm mem} be the sum of the negative logarithms of the conditional probabilities of all observed membership outcomes along the true path through live termination. Let X_* =\max\log^-m over its processing states, including the state at transfer or at the common stop, but excluding later static descendants. Then

    \|B_{\rm mem}\|_{L^2(\mathbf P)} +\|X_*\|_{L^2(\mathbf P)}\le C h, \tag{6.3}

    and \mathbf P(\mathrm{transfer})\le C h^{-10}.

  3. (iii)

    The live plateau weights satisfy

    \sum_{e\text{ live}}\mathbf U(p_e)\alpha_e\le C h.

    If \xi_p is the local average of mass removed from the survivor of one plateau at p, divided by its \alpha_e, then \mathbb E_{\mathbf U_e}\sum_p\xi_p\le1. Full removal at an exit may be included. Starting from a final jumping post state p_0, the corresponding future bound under uniform probability on p_0 is at most m_{p_0}/\alpha_e\le2.

  4. (iv)

    The global label Fisher sum obeys

    \sum_{\text{all live jumping pieces }p} \mathbf U(p)\alpha_e J_p\le C d. \tag{6.4}
  5. (v)

    Ignore the labels and use the filtration of spatial bits and memberships. After live termination one may keep the final piece fixed and reveal any finite number of additional spatial bits, making no further membership observations. The expected sum of the spatial bit relative entropies to a fair bit is O(d). If this spatial extension is stopped at the common density cutoff, the bound is O(h).

Proof.

Caps, constant order, and finite termination

Let C_1 be the numerical constant bounding a relative metric cost by C_1t in Lemma 5.3, and let C_2 be its numerical regular gap-increase constant. Choose the numerical \varepsilon>0 there small enough for the cheap gap decrease, and take R_0=1000. Consider a prefix of successes for which the update hypotheses have been justified. The gap potential

F(A)-\log|B_A(r)|

starts at zero and is bounded below by an absolute constant. Here r begins at r_{\rm sep}/4 and increases by r_{\rm sep}/(10C_0h) at each addition, so r\le r_{\rm sep}/2 under the count cap. Telescoping its regular increases and cheap decreases gives

\sum_{\rm cheap}k_b\le 2C_2\sum_{\rm regular}t+C, \qquad \sum_{\rm successes}t\le (1+2\varepsilon C_2)\sum_{\rm regular}t+C. \tag{6.5}

In the second inequality we used t=\varepsilon k_b on cheap successes. Regular successes also satisfy \varepsilon k_b\le t. Consequently, while \mathcal Q\le R_0d, the sum of all success k_b’s is at most C(R_0,\varepsilon,C_2)d. Since every k_b\ge d/h, the number of those successes is at most C'h. Choose C_0>C'+2. A first count breach without a charge breach would obey this same prefix bound, contradicting that choice for large d.

The form cap is also inductive. If the preceding main addition occurred at p_*, the spatial and metric estimates give

F(A^+)-F(A) \le\log(R_p/R_*)+C+C_1t.

The spatial scale ratios telescope along a branch. Before a charge breach, Equation (6.5) and the count estimate therefore give F(A)\le\log R_*+C_*d for a sufficiently large fixed C_*. The initial case is A=0 with p_*=I. At the next update, the flat-increment condition Equation (5.5) is supplied by the update itself. Thus the assumptions needed to construct each successive update have been established from its predecessors.

Fix D_0,P_0 in the compression application using these caps, then fix A_0 sufficiently large for the cardinality loss, and finally take d sufficiently large for all required asymptotic inequalities. The constants C_1,C_2 and the choice of \varepsilon precede the caps; their numerical nature is what permits this order. A regular success that first crosses the charge cap adds at most D_0d and is followed immediately by transfer. Hence that final overshoot only changes the absolute constants in the claimed path budgets. No subsequent update uses a form cap for this terminal form. The energy estimate in Equation (5.7) remains valid for the labels drawn at that success, as it uses the flat-increment and count bounds. This proves all the deterministic budget assertions.

At a fixed spatial node, the number of successes along any recurrence path is bounded by the count cap. Between successes, every failure multiplies m by at most 1-c(b_0/8)^2. Since m\le e^d, only finitely many consecutive failures are possible before m<e^{-d}. There is a uniform finite bound, depending on d but not the node, on the length of this local binary recursion. Its branching is finite; induction over the finite spatial depth proves termination of the entire procedure. Integer frequencies and finite offending collections can be chosen in a fixed enumeration. Polarization uses measurable bounded signs of the moment functions, and all membership probabilities are computed from their channel laws. Thus all resulting input sets are measurable. No joint choice over different dyadic grids is required.

The spatial likelihood and membership costs

At a positive-mass discrete state with density a on p,

\mathbf P(\text{state})= \frac{\mathbf U(p)m}{m_{\rm top}},\qquad \mathcal L(Y\mid\text{state})= \frac{a}{m}\,\mathcal L(\mathbf U\mid p). \tag{6.6}

This follows directly by restricting the initial density to the recorded spatial interval and membership sets. Define, at each such state,

L_{\rm sp}=\log(m/m_{\rm top})+B_{\rm mem}.

A membership outcome with probability q>0 replaces m by qm and increases B_{\rm mem} by -\log q, so L_{\rm sp} is unchanged. At a spatial step its increment is \log(2q), where q is the conditional probability of the actual child under the unlabelled true law. It follows that

\mathbb E_{\mathbf P}\!\left[ e^{-L_{\rm sp,new}}\mid\text{past}\right] \le e^{-L_{\rm sp,old}}.

Indeed the conditional sum over positive-probability children is \sum_q q/(2q)\le1; equality need not hold if a child has zero mass. Thus e^{-L_{\rm sp}} is a nonnegative supermartingale starting at one. It remains so after stopping or after keeping the piece fixed and revealing further spatial bits.

A success has surprisal at most C(1+\log(1/b)) because \delta\ge cb^2; a cheap success has zero surprisal. The deterministic k_b budget and

1+\log(1/b)\le C(h/d)k_b

bound the total success surprisal by Ch on every path. Between one success and the next, let F accumulate only failure surprisal. Until that stretch ends, e^F is unchanged at spatial steps; at a membership test keep it multiplied by (1-\delta)^{-1} on failure and send it to zero on success. The resulting process is a nonnegative supermartingale starting at one. Therefore its maximal inequality implies, conditionally on the stretch’s start,

\mathbf P(F_{\rm final}>s\mid\text{start})\le e^{-s}, \qquad s\ge0.

Here F_{\rm final} includes the failures preceding a terminal success, even though the auxiliary process is then sent to zero. There are at most Ch+1 stretches. Their conditional second moments are at most 2; pad with zero stretches and use the triangle inequality in L^2. This proves \|B_{\rm mem}\|_2\le Ch. The same maximal inequality applied to e^{-L_{\rm sp}} bounds the L^2 norm of \max(-L_{\rm sp})_+ by an absolute constant. Since m_{\rm top}\ge1/h,

X_*\le\log h+B_{\rm mem}+\max(-L_{\rm sp})_+,

which proves Equation (6.3), including the mean at transfer.

We shall also need the expected number N_{\rm att} of regular attempts. Each has conditional success probability at least c(b_0/8)^2\ge c'h^{-32}, whereas there are at most Ch successes. Summing conditional success probabilities over the finite processing tree gives

\mathbb E_{\mathbf P}N_{\rm att}\le C h^{33}. \tag{6.7}

Likelihood evidence and the transfer probability

The small-mass flag implies X_*>d, so its probability is at most Ch^2/d^2=o(h^{-10}). The count flag is already accounted for by the charge flag. To control the latter, attach all the main and evidence observations specified above. Compare the true experiment with a reference law \mathbf Q in which the prior for Y is \mathbf U, the spatial bits and all channel draws still observe Y, but each membership outcome is an independent toss with the state’s prescribed probabilities \delta,1-\delta. On a recorded branch those probabilities are constants in Y. All subsequent states and channel parameters are the same predetermined functions of its spatial and membership history as in the true experiment. One may stop any zero-true-mass branch in the reference as well.

Let \mathcal O denote the terminal observations through live termination, including draws at its final success, and let R=d\mathbf P_{\mathcal O}/d\mathbf Q_{\mathcal O}. On true histories,

R\le\exp(d+\log h+B_{\rm mem}). \tag{6.8}

To verify this even for histories of different lengths, fix one complete discrete history and all its label values, using densities for the latter. The true likelihood integrates the same channel weights over inputs compatible with both the spatial bits and the actual membership sets, with prior density at most e^d/m_{\rm top}. The reference integrates over all inputs compatible with the spatial bits, and its additional coin factors have product e^{-B_{\rm mem}}. The domain of the true integral is a subset of that reference domain. This proves Equation (6.8) on each terminal branch separately.

At every regular attempt whose observed past is admissible, multiply a score by the following factor, writing E_{\rm ev} for the evidence event in Equation (5.21):

\frac{\exp\bigl((t/16)\mathbf1_{ \{\mathrm{success}\}\cap E_{\rm ev}}\bigr)} {1+e^{-t/16}};

otherwise leave the score unchanged. Group the fresh observations of this test into that step. The conditional reference probability of the evidence given success is at most e^{-t/8} by Proposition 5.4. Consequently the conditional mean of the numerator is at most 1+(e^{t/16}-1)e^{-t/8}\le1+e^{-t/16}. The score \mathcal S is therefore a nonnegative reference supermartingale, with \mathbb E_{\mathbf Q}\mathcal S\le1 at termination. Change of measure gives the particularly useful bound

\mathbf P(\mathcal S>Re^d) =\int_{\{\mathcal S>Re^d\}}R\,d\mathbf Q \le e^{-d}\mathbb E_{\mathbf Q}\mathcal S\le e^{-d}. \tag{6.9}

The simultaneous true admissibility and evidence assertion of Proposition 5.4 fails with probability o(h^{-20}) under the count cap; it includes evidence at a success that immediately transfers. Also \mathbf P(B_{\rm mem}>d)\le Ch^2/d^2 by Equation (6.3). Finally, the total logarithmic denominator cost in the score is at most N_{\rm att}e^{-c d/h}. Its expectation tends to zero faster than h^{-10} by Equation (6.7), so with that allowed exception it is at most one. Outside these exceptions and Equation (6.9),

\frac{\mathcal Q}{16}-1 \le\log\mathcal S \le\log R+d \le2d+\log h+B_{\rm mem} \le3d+\log h.

This excludes \mathcal Q>1000d for large d. It proves the stated transfer probability. The score and evidence draws have only served to estimate this probability; they never select a piece or an update.

The sum of live plateau weights

At a plateau’s root its mean is \alpha_e, so Equation (6.6) gives

\sum_{e\text{ live}}\mathbf U(p_e)\alpha_e =m_{\rm top}\, \mathbb E_{\mathbf P}[\text{number of live plateau starts}].

Epoch starts number at most Ch+1 on each path. It remains to bound upward resets without multiplying that bound by another h.

Consider one epoch, with start and end stopping times \sigma,\tau in the unlabelled filtration. Include its ending membership success in \tau if there is one. Write \Delta=L_{\rm sp}(\tau)-L_{\rm sp}(\sigma), and let N be its number of upward resets. At its jth reset time \rho_j,

L_{\rm sp}(\rho_j)\ge L_{\rm sp}(\sigma)+j\log2:

the normalizations have more than doubled at each reset and the membership surprisal is nondecreasing. For any fixed r>0, optional sampling of the inverse likelihood gives

\mathbf P\bigl(L_{\rm sp}(\tau)<L_{\rm sp}(\rho_j)-r \mid\mathcal F_{\rho_j}^{\rm disc}\bigr)\le e^{-r}, \qquad \mathbb E[(-\Delta)_+\mid\mathcal F_\sigma^{\rm disc}]\le1.

Here \mathcal F^{\rm disc} denotes the spatial and membership history, not the uniform-output filtration used earlier. The second bound follows by integrating the corresponding exponential tail. Every reset not in the exceptional event in the first bound has j\log2\le\Delta+r. Summing reset indicators and then conditioning at the epoch start yields

(1-e^{-r})\mathbb E[N\mid\mathcal F_\sigma^{\rm disc}] \le\frac{\mathbb E[(\Delta+r)_+ \mid\mathcal F_\sigma^{\rm disc}]}{\log2} \le\frac{\mathbb E[\Delta\mid\mathcal F_\sigma^{\rm disc}]+r+1} {\log2}.

Sum this inequality over epochs, weighting by the probability that each epoch starts, and set r=2. From one epoch’s end to the next epoch’s start L_{\rm sp} does not change, because only a membership step separates them. Thus all \Delta terms telescope to the terminal L_{\rm sp}. At termination,

m\le2e^{4h},\qquad L_{\rm sp}\le4h+\log2+\log h+B_{\rm mem}.

The factor two allows the last spatial step that first crosses the common cutoff; membership restrictions cannot increase the mean. The expected total number of resets is therefore O(h) by Equation (6.3). Since m_{\rm top}\le2, this proves the plateau weight bound.

For a fixed plateau, the portions removed from its survivor before its end are disjoint subsets of its original input mass. Multiplying each local removed mean by \mathbf U_e(p) and summing counts each removed input at most once. Division by the entry mean \alpha_e gives \mathbb E_{\mathbf U_e}\sum\xi_p\le1. The identical argument from a final post state p_0, excluding its past removals, starts with mass m_{p_0}/\alpha_e. This proves the remaining assertions in part (iii).

The information cost of main labels

We first recall the finite entropy accounting used here. For any finite remaining spatial bit string, track its conditional relative entropy to uniform future bits. Revealing the next bit consumes, in conditional expectation, precisely that bit’s relative entropy to fair, by the chain rule. An intervening observation increases the expected deficit by its conditional information about the remaining string, at most its conditional information about Y. The remaining deficit is nonnegative, so the same inequality holds when the process is stopped. The initial deficit is at most

D(\mathbf P\|\mathbf U) =\mathbb E_{\mathbf P}\log(\rho/m_{\rm top}) \le d+\log h.

For a membership bit its information cost is at most its conditional entropy. Averaging and removing all observed-label conditioning only increases this entropy, so the total membership cost is at most \mathbb E B_{\rm mem}=O(h).

For the main-label filtration we retain all main observations, but not the evidence clones. The initial type-0 observation has zero information. At a success at p, let Z be the old main observation and V the new observation of type A^+. Given the full main-observation history through that success, compare the conditional channel of V with the probability density

Q_Z(v)=\frac{\exp[-A(v-Z)-S_p(v-c_p)]} {\int\exp[-A(w-Z)-S_p(w-c_p)]\,dw}, \qquad c_p=\operatorname{center}(p).

For every conditional prior on Y, its information about V is bounded above by the expected relative entropy of the channel to any fixed comparison density. The centered envelope bounds the logarithm of the denominator above by -F(A+S_p). Expanding that relative entropy and dropping the nonpositive new-form energy gives

\begin{aligned} \mathrm I(Y;V\mid\text{history through success}) \le{}& F(A^+)-F(A+S_p)\\ &+\mathbb E[A(V-Z)+S_p(V-c_p) \mid\text{history through success}]. \end{aligned} \tag{6.10}

The notation denotes the usual averaged conditional-information bound; it may first be read for each fixed value of the history.

Now average this inequality over the main-observation histories at the fixed discrete success state. Conditional on Y and that discrete state, the old and new errors still have their prescribed independent laws, because every membership decision is input-only. The phase triangle inequality and form monotonicity imply

\mathbb E A(V-Z) \le2\int A K_{A^+}+2\int A K_A \le4\int A K_A \le C h(dh^{-400}+1).

This averaging step is essential: conditioning on an observed old label alone would not give its unconditional error law. Also Y\in p and A^+ contains S_p, so

\mathbb E S_p(V-c_p) \le2\mathbb E S_p(V-Y)+2S_p(Y-c_p)\le C.

Multiply each such estimate by the true probability of its success state and sum. Those probabilities sum to the expected number of successes, at most Ch. Their total old-error energy is therefore

O\bigl(h^2(dh^{-400}+1)\bigr)=o(d),

and their spatial energy is O(h). The relative metric terms in Equation (6.10) sum to O(d) on every path by part (i). Thus all incremental main-channel informations together cost O(d), including a terminal cap-crossing success.

The Fisher sum and the unlabelled deficits

At a final live state p, the conditional density of the current main label, ignoring the older labels, is \alpha_e M/m_p. Given this label, the next spatial bit has bias D/M. Hence its expected squared bias at this state is

\frac{\alpha_e}{m_p}\int\frac{|D|^2}{M} =\frac{\alpha_e}{m_p}J_p.

Multiplication by the state’s probability in Equation (6.6) gives exactly \mathbf U(p)\alpha_e J_p/m_{\rm top}. Conditioning also on all older main labels can only increase the average squared conditional bias, by conditional Jensen. For a binary law of bias b, its relative entropy to a fair bit is at least b^2/2. The finite entropy accounting just proved therefore bounds the sum of these weighted J_p’s by a constant times m_{\rm top}(d+\log h+\mathbb E B_{\rm mem}+O(d))\le Cd. This proves Equation (6.4) with no lower bound on the surviving mass needed in this conversion.

Finally ignore all labels. The conditional expected increment of L_{\rm sp} at a spatial bit is exactly its relative entropy to fair, while memberships leave L_{\rm sp} unchanged. Starting from L_{\rm sp}=0, the expected sum of bit deficits equals its expected terminal value. On any finite spatial extension with the last piece fixed, m\le e^d, so

L_{\rm sp,final}\le d+\log h+B_{\rm mem}.

If that extension stops at the common density cutoff, instead m\le2e^{4h} and L_{\rm sp,final}\le4h+\log2+\log h+B_{\rm mem}. Taking expectations and using Equation (6.3) gives respectively O(d) and O(h). This completes the proof.

7 Capture, edge charges, and snapshots

The frequency lists in this section serve two different purposes. A cumulative list gives a deterministic approximation by signed sums. At selected spatial prefixes we compress the relevant signed sums into a smaller list of representatives. The first list controls how long approximation can fail; only the smaller lists will enter the probability estimates of Section 8. All constructions take place in the finite procedure of Section 6, uniformly in its depth and in any frequency cutoff.

7.1 The two laws on a fixed plateau

For static pieces, start a plateau at entry to static handling and restart whenever the incoming mean exceeds twice its current normalization. There are no membership tests on these plateaus. They obey the same common stopping rule as live pieces. Discard zero pieces.

Fix a live or static plateau e, with root p_e and normalization \alpha=\alpha_e>0. Under \mathbf U_e, spatial uniform probability on p_e, follow only this plateau’s surviving density, through its failures, until its exit or the common stop. Its jumps on an output path occur at consecutive depths. At a jumping node p, write a for the final post density, after all removals there, and m_p=\mathbb E_p a. The two prospective child measures use this same density restricted to the children, before any processing in a child. Thus

m_p/\alpha\le2, \qquad \mathbb E_{p'}a/\alpha\le4 \quad(p'\text{ a child of }p). \tag{7.1}

Let \xi_p be the mass removed at p, measured as its local spatial average divided by \alpha. Full removal at an exit can be included. The removed input sets are disjoint, so

\mathbb E_{\mathbf U_e}\sum_p\xi_p\le1, \qquad \mathbb E_{\mathbf U_{p_0}}\sum_{p\text{ after the post state }p_0}\xi_p \le m_{p_0}/\alpha\le2. \tag{7.2}

The second sum excludes removals already made at p_0. Here and below a sum along a path includes only states of the specified plateau. Put

d_1=d\quad\text{on a live plateau},\qquad d_1=d+2+\log^-\alpha\quad\text{on a static plateau}.

For comparison, let \mathbf P_e be the input experiment starting with the entry piece normalized to probability on p_e. It observes the spatial bits of its input Y and its actual memberships; hence it can leave the plateau on a success. It may subsequently reveal more spatial bits, with no further observations, when a finite extension is useful. A final jumping post state has probability

\mathbf P_e(\text{post state at }p) =\mathbf U_e(p)\frac{m_p}{\alpha}. \tag{7.3}

Conditional on that state, the input density is a/m_p relative to uniform on p. If the experiment instead begins at a jumping post state p_0, the corresponding probability of a later post state p\subseteq p_0 is \mathbf U_{p_0}(p)m_p/m_{p_0}. These identities follow directly because each survivor is the entry density restricted by its spatial prefix and membership events.

The uniform output law does not condition on a random input membership. At each reached output prefix, the survivor and all removals are already determined functions. Both prospective children can therefore be examined before the next fair output bit. For a live plateau the side information will be its fixed main label Z, with conditional density K_A(z-Y) given Y. The discrete decisions do not inspect its realized value. For a scalar construction there is no side information.

Lemma 7.1 (Local entropy and Fisher budgets). On the plateau just described, the expected sum of spatial bit relative entropies to a fair bit is O(d_1), with or without the fixed side label. This includes every finite extension after exit. If finite-valued observations are made at predictable states and their total conditional entropy cost has expectation at most B, the bound becomes C d_1+B.

Let M be the normalized joint post mass measure of the side information and the finite observations so far, and let D_{\rm jump} be the half-difference of its prospective child measures. Then

\mathbb E_{\mathbf U_e}\sum_p \int\frac{|D_{\rm jump}|^2}{M}\le C(d_1+B). \tag{7.4}

For a forward estimate from a live jumping post state p_0, restart the input experiment with its post density, retaining the fixed main label when present but discarding all earlier auxiliary observations. The same bounds then hold with d_1=d, with B charging the auxiliary observations made in this restarted experiment, and with an extra factor m_{p_0}/\alpha\le2 in the mass conversion. Already completed removals are excluded. In particular the future main-label Fisher sum has conditional expectation at most Cd. If earlier auxiliary information W is retained instead, its conditional information \mathrm I(Y;W\mid Z) under the restarted input law must also be added to the initial budget; omit Z when there is no side label. All quotients are zero at zero mass; for measures, densities to a common dominating measure are understood.

Proof. The initial density deficit relative to spatial uniform is at most d-\log\alpha\le C d_1. On a live plateau with jumps, \alpha\ge e^{-d}. To pay for its fixed label, let p_* be the interval of the last main addition and c_* its center. Compare its channel with the fixed density K_{S_{p_*}}(z-c_*). Writing Z=Y+\theta gives

\begin{split} \mathrm I(Y;Z)&\le\mathbb E_Y D\bigl(K_A(z-Y)\,dz\,\|\,K_{S_{p_*}}(z-c_*)\,dz\bigr)\\ &=F(A)-F(S_{p_*})-\mathbb E A(\theta) +\mathbb E S_{p_*}(Y+\theta-c_*)\\ &\le F(A)-\log R_{p_*}+C\le Cd. \end{split}

Indeed p_e\subseteq p_*, the channel includes S_{p_*}, and the spatial variance and form monotonicity in Section 5 bound the last energy by a constant: use the phase triangle inequality and |Y-c_*|\le |p_*|/2. The negative energy can be dropped. The initial dummy channel has zero information. This estimate holds for every input law supported on p_e, so also for the post law at a later live jumping state.

Within one live plateau a true path has at most one success, which is terminal for that plateau. Its success probability at a test is at least c b_0^2, so its success surprisal is O(1+\log h). The accumulated failure surprisal has an exponential tail by the failure-stretch supermartingale used in Section 6. Thus the expected total membership entropy cost is O(1+\log h). Static plateaus have no such cost. Forward from a live post state, m_p\ge e^{-d}, so the initial density deficit is again at most 2d, and the same membership argument applies. This forward argument begins without earlier auxiliary observations. Retaining such information W instead increases the initial information term from \mathrm I(Y;Z) to \mathrm I(Y;Z,W)=\mathrm I(Y;Z)+\mathrm I(Y;W\mid Z).

For completeness, apply the entropy chain rule to any finite remaining spatial bit string. Its expected deficit relative to the uniform string starts below the density deficit, plus the information of any label already observed. Revealing a bit consumes its conditional deficit relative to a fair bit. An additional observation increases the expected remaining deficit by its conditional information about the remaining string, at most its conditional Shannon entropy if finite-valued. The remaining deficit is nonnegative at stopping. Membership entropy can be bounded before conditioning on side and auxiliary observations, using the discrete state history that prescribes the tests. This proves the first assertion, including finite extensions.

At a post state p, the joint observations have conditional density (\alpha/m_p)M under truth and the next bit bias is D_{\rm jump}/M. Their averaged squared bias is therefore

\frac{\alpha}{m_p}\int\frac{|D_{\rm jump}|^2}{M}.

Observing further past information only increases the mean squared conditional bias, by conditional Jensen. For a bit with bias \theta, its relative entropy to fair is at least \theta^2/2. Multiply by Equation (7.3) and sum: the factors m_p/\alpha and \alpha/m_p cancel. This proves Equation (7.4) without any lower bound on surviving mass. Starting at p_0 instead leaves the factor m_{p_0}/\alpha\le2, as asserted.

7.2 The capture statement

There are two separate constructions at each allowed dyadic level b. For the main criterion, used on live plateaus for b_0\le b\le1, set \kappa=2 and say that u is high if \max_{p'}G_{p'}^u\ge b^2, using the prospective labeled child moments. For the scalar criterion, set \kappa=1 and say that u is high if \max_{p'}|f_{p'}^u|\ge b, where f_{p'}^u=\mathbb E_{p'}[a\chi_u]/\alpha has no label. Scalar levels are 0<b\le b_0 on live plateaus and 0<b\le1 on static plateaus. In particular the main and scalar constructions at b_0 are distinct.

Set

l_b=\left\lceil20\log_2\frac{d_1+2}{b}\right\rceil, \quad \delta_p=2^{-l_b}|p|, \quad e_{p,b}=\frac{\mathbb E_p[a\mathbf1_{\mathrm{edge}(p,b)}]}{\alpha}, \tag{7.5}

where \mathrm{edge}(p,b) consists of points inside p within 4\delta_p of an endpoint or the midpoint. The jump is bad for this criterion and level if e_{p,b}\ge c b^\kappa, for a sufficiently small fixed positive c. Everything in this definition is known at the post state before the next output bit.

Proposition 7.2 (Capture on a plateau). In the finite procedure, for every plateau, allowed level, and criterion above, there is a cumulative finite list D of integers, initially empty, whose appends are prescribed at jumping post states, such that, after the appends at p,

\begin{gathered} u\text{ high and the jump not bad}\quad\Longrightarrow\quad |u-v||p|\le\poly(1+|D|,d_1,1/b) \quad\text{for some }v\in[D],\\ [D]=\left\{\sum_{n\in D}\eta_n n:\eta_n\in\{0,1,-1\}\right\}, \qquad \mathbb E_{\mathbf U_e}|D|_{\rm final}\le C d_1 b^{-C}. \end{gathered} \tag{7.6}

The construction is independent of the caps used below.

The bad-edge charge at p and level b is e_{p,b}\mathbf1_{\{p\text{ bad at level }b\}}. Its global live sum, with weights \mathbf U(p)\alpha_e, is O(d) when one scalar level is selected predictably at each piece, and O(\poly(h)(h+\log d)) when all main levels are summed. The local one-scalar-level charge is O(d_1) under \mathbf U_e. Forward from a live post state the conditional expected number of bad jumps at a fixed main level is at most d\poly(h).

For each deterministic dyadic D_{\max}\ge1, the scheme stopped when |D|>D_{\max} admits snapshots with expected number \poly(D_{\max},d_1,1/b). A snapshot at level b stores at most \lceil e^{k_{b/2}}\rceil main representatives, or \exp(C(d_1+\log(1/b))) scalar representatives. Each birth and its list are determined by the spatial prefix before its birth jump.

At a jump p, call a high integer u activated when some v\in[D] satisfies |u-v||p|\le\sigma in the main case, or |u-v||p|\le c_8b^4 in the scalar case, for a sufficiently small fixed c_8>0. At activation its chosen representative s has, for the remainder of that plateau, a fixed decomposition

\begin{array}{ll} u=s+n+a',\quad \Delta_A(n)\le r_{\rm sep},\quad |a'||p_{\rm act}|\le C\sigma,&\text{main},\\ u=s+a',\quad |a'||p_{\rm act}|\le c_1b^4,&\text{scalar}, \end{array}

where p_{\rm act} is the activation node and c_1 can be made a sufficiently small absolute constant. Before first activation, all non-bad high jumps for this u lie in a single block of at most

H_b=C\log(CD_{\max}d_1/b)+C\log(1/\sigma)

consecutive depths, as long as the cap is not breached. The construction does not use sampled side or auxiliary values on an output path.

We prove the assertions in turn. The intermediate lemmas also specify the conditioning needed for their later use.

7.3 Broad capture by trimmed masks and polynomial products

Suppose u is high at a jumping post piece. Polarize on one high child: there is a measurable complex v(z) with |v(z)|\le1 (a constant in the scalar case) such that the normalized real bias

\operatorname{Re}\frac1\alpha\mathbb E_p [a\mathbf1_{p'}(Y)v(Z)\chi_u(Y)]\ge b^\kappa/2. \tag{7.7}

For the main criterion, |f_{p'}^u|\le M_{p'} implies \int|f_{p'}^u|\ge G_{p'}^u; choosing its pointwise conjugate phase proves the assertion. In particular every append state below satisfies m_p/\alpha\ge b^\kappa/2.

We first replace the child indicator by a trigonometric polynomial. Write n=|D| at this stage. Trim a distance \delta_p from both ends of p', and convolve its indicator on the unit circle with the positive normalized kernel proportional to

\left(\frac{\sin(\pi N x)}{N\sin(\pi x)}\right)^{2s}, \qquad N=\left\lceil C2^{l_b}/|p|\right\rceil, \qquad s=\left\lceil C(1+n+\log(1/b))\right\rceil.

Its continuous values at integers are understood. The result P_p is a real trigonometric polynomial, 0\le P_p\le1, of degree at most sN. The unnormalized kernel integral is at least c/(N\sqrt{s}): on circle distance at most c/(N\sqrt{s}) from zero, the ratio has 2s-th power bounded below by a positive numerical constant. The bound |\sin\pi x|\ge2\operatorname{dist}(x,\mathbb Z) gives its unnormalized tail outside distance \delta_p at most Ce^{-c's}/N, after choosing the constant in N large. Consequently the normalized tail is at most C\sqrt{s}e^{-c's}. Outside the edge inside p, the mask approximates \mathbf1_{p'} uniformly with this error. Integrating the convolution in the other order also gives

\deg P_p\le \frac{C2^{l_b}(1+n+\log(1/b))}{|p|}, \qquad \int_{\mathbb T\setminus p}P_p \le |p|\exp[-C'(1+n+\log(1/b))]. \tag{7.8}

Here C' can be chosen as large as needed by increasing the fixed kernel constants. The edge error in Equation (7.7) is at most e_{p,b}, and the remaining error is at most 2C\sqrt{s}e^{-c's}. Thus, at a non-bad jump, the replacement retains bias \ge c_2b^\kappa. This argument also covers intervals that cross the circle coordinate boundary.

For each fixed side value z, keep a full-circle polynomial product Q_z, initially one. Its restriction to the current p, normalized by its integral there, is the reference probability density for Y conditional on that side value. After n appends its factors lie between 1/2 and 2, so

2^{-n}\le Q_z\le2^n,\qquad \int_pQ_z\ge2^{-n}|p|.

If a non-bad high u is outside [D] enlarged by the sum of all previous mask bandwidths and the prospective bandwidth in Equation (7.8), append it and multiply Q_z by

1+\eta b^\kappa\operatorname{Re}[v(z)\chi_u P_p], \tag{7.9}

where \eta>0 is a sufficiently small fixed constant. The old product has Fourier support inside [D] enlarged by its accumulated mask bandwidths. Thus the full-circle integral of Q_z\chi_u P_p vanishes. Restriction to p changes its normalized integral by at most

\frac{2^n\int_{\mathbb T\setminus p}P_p}{2^{-n}|p|} \le4^n\exp[-C'(1+n+\log(1/b))].

Choose C' so this is a suitably small multiple of b^\kappa, uniformly in n and z.

The relevant entropy is the true average, over the current conditional side law, of the relative entropy of Y conditional on that side value to this reference. Its drop at an append equals the true expectation of the logarithm of Equation (7.9), minus the average logarithm of the reference expectation of that factor. By Equation (7.1), the retained normalized bias becomes at least c_2b^\kappa/2 when divided by m_p/\alpha to obtain a probability-law bias. The inequality \log(1+x)\ge x-Cx^2 for |x|\le1/2 therefore bounds the first logarithm below by c\eta b^{2\kappa}, after fixing \eta small. The leakage estimate bounds the normalization logarithm by a smaller multiple of the same quantity. Every append hence consumes at least c_3b^{2\kappa} of this conditional entropy.

Here is the budget across states, including the change in the side law. Denote a discrete state history by H, and the current reference by Q_{H,z}, restricted and normalized on its interval. Set

\mathcal B(H)=\mathbb E_{Z\mid H} D(\mathcal L(Y\mid H,Z)\|Q_{H,Z}).

At a spatial bit X, restrict the reference to the selected child. The relative entropy chain rule gives exactly

\sum_x\mathbb P(x\mid H)\mathcal B(H,x) =\mathcal B(H)-\mathbb E_{Z\mid H} D(\mathcal L(X\mid H,Z)\|\mathcal L_{Q_{H,Z}}(X)) \le\mathcal B(H).

At a membership observation C, keep the old reference on both outcomes. Its average entropy increment is \mathrm I(Y;C\mid H,Z)\le\mathrm H(C\mid H), since the current side posterior is used on each outcome. Nonnegative entropy on exited branches can be discarded. Initially the budget is D(\mathcal L(Y)\|\mathbf U_e)+\mathrm I(Y;Z)\le Cd_1; the membership costs are those in Lemma 7.1. Thus

\mathbb E_{\mathbf P_e}(\text{number of appends}) \le Cd_1b^{-2\kappa},\qquad \mathbb E_{\mathbf U_e}|D|_{\rm final} \le Cd_1b^{-3\kappa}. \tag{7.10}

The second estimate uses Equation (7.3) and m_p/\alpha\ge b^\kappa/2 at every append. At a fixed positive-mass state the entropy is finite and drops strictly, so the append loop is finite. This finiteness propagates down every positive-probability branch of the finite processing tree.

If no further append is possible, every non-bad high u is within the accumulated bandwidth of [D]. Previous append intervals are ancestors of p, and previous list lengths are at most the current length. Hence that bandwidth, multiplied by |p|, is bounded by

C2^{l_b}(n+1)(n+1+\log(1/b)).

For example C(d_1+2)^{20}b^{-21}(n+1)^2 is a permissible fixed polynomial upper bound. This proves Equation (7.6). The signs, products, and appended frequencies are chosen from the laws at each state; no realized side value is used to decide an append.

7.4 Charging thin spatial edges

Lemma 7.3 (Edge charges). With the bad-jump convention Equation (7.5), the following bounds hold.

  1. (i)

    Across all live jumping pieces, the sum of \mathbf U(p)\alpha_e e_{p,b}\mathbf1_{\rm bad} is at most Cd if a single scalar level b is selected at each piece before its next output bit, without using noisy label values.

  2. (ii)

    The corresponding sum over all main levels and live pieces is at most C\poly(h)(h+\log d).

  3. (iii)

    Within one plateau, the expected one-scalar-level sum under \mathbf U_e is at most Cd_1.

  4. (iv)

    From any live jumping post state, under spatial uniform conditional on that prefix, the expected number of future bad jumps of the same plateau at main level b is at most Cd b^{-2}, and hence d\poly(h).

Proof. For (i) and (ii), fix the global unlabelled input experiment of Proposition 6.2, following memberships through all live epochs and freezing the terminal piece only at live termination. For (iii) and (iv), first fix one plateau and use the corresponding local input experiment of Lemma 7.1, freezing the terminal piece at that plateau’s exit. In each application, all windows below belong to one common finite spatial extension of this fixed experiment and use the same predictable spatial entropy sequence.

At a final post state the conditional input edge probability is

q=\frac{\alpha e_{p,b}}{m_p}\ge c b^\kappa/2 \quad\text{if the jump is bad}.

From this state compare the next l=l_b spatial bits, starting with its jump bit, to a reference experiment with fair spatial bits. At intervening membership tests use the very same state-conditional outcome probabilities as in the true unlabelled experiment. Thus membership likelihood factors cancel. This reference can be generated by sampling each bit fairly, regardless of earlier memberships, and by tossing the stated membership coins on their recorded branches. In particular its spatial path remains uniform. Continue bits after exits or stops as needed, keeping the piece fixed and taking no new observations. Histories of zero true mass can be extended arbitrarily under the reference while leaving its bits fair.

The edge is determined by the next l bits, up to null endpoints, and its reference probability is at most 16\,2^{-l}. If L is the true-to-reference likelihood ratio of this segment, then

\mathbb P_{\rm true}(\mathrm{edge}\cap\{L<2^{l/2}\}\mid\text{post state}) \le16\,2^{-l/2}\le q/2.

The last inequality follows uniformly from the choice of l_b, the badness lower bound on q, and sufficiently large d. Call the starting depth marked if the edge occurs and \log L\ge(l/2)\log2. Conditional on a bad post state, its probability is at least q/2.

Let s_j be the predictable relative entropy to fair of the spatial bit at depth j, in this unlabelled true history. If its bias is \theta_j, then s_j\ge\theta_j^2/2, whereas the positive part of its realized log likelihood increment satisfies

\bigl[\log(1\pm\theta_j)\bigr]_+ \le |\theta_j|\le\sqrt{2s_j}.

Memberships contribute no log increment to L. On a mark, Cauchy–Schwarz on the length-l window consequently yields

\sum_{j\text{ in its window}}s_j\ge c_4 l. \tag{7.11}

There is at most one final jumping piece at each spatial depth along an actual input path, although membership splitting may have preceded it. For any finite family of marked starting indices, select disjoint windows whose triples cover the union, for example by successively taking a longest remaining window. The number of starting indices is no larger than a constant times the sum of the selected window lengths. By Equation (7.11), this is at most C\sum_j s_j, since the selected windows are disjoint. This proves the interval charging bound for variable window lengths as well as for fixed l.

Globally, the unlabelled full spatial entropy budget of Proposition 6.2 is O(d) when live pieces occur. Apply the interval bound with the one selected scalar level per index. The marking probability and the state law imply

\sum_{\rm live\ pieces}\mathbf U(p)\alpha_e e_{p,b}\mathbf1_{\rm bad} =m_{\rm top}\mathbb E_{\mathbf P} \sum_{\rm live\ post\ states}q\mathbf1_{\rm bad} \le2m_{\rm top}\mathbb E_{\mathbf P}(\text{number of marks}) \le Cd.

Only a finite spatial extension is needed for any finite selection of windows; the bound for the full sum follows by nonnegativity. This also handles arbitrarily small selected scalar levels.

For a fixed main level, l_b is constant along the live path. Windows wholly before the common cutoff use only its O(h) truncated unlabelled bit budget. At most l_b starting indices straddle that cutoff. The truncated budget remains valid on paths that transferred earlier: keep their membership state fixed and continue their spatial sampling. Therefore the main-level weighted edge sum is O(h+l_b). There are O(\log h) main levels and l_b=O(\log d+\log(1/b_0)) on them, proving (ii).

For (iii), use \mathbf P_e and Lemma 7.1, with no side label, in the same proof. The state conversion now has root mass normalization one, giving Cd_1. Starting instead at a live post state p_0, the conversion leaves a factor m_{p_0}/\alpha\le2, and the forward entropy budget is Cd. Thus the expected forward sum of e_{p,b}\mathbf1_{\rm bad} under \mathbf U_{p_0} is at most Cd. At a bad main jump e_{p,b}\ge cb^2; dividing gives (iv).

7.5 Snapshots in fixed rounded coordinates

Fix a deterministic dyadic cap D_{\max}\ge1 and stop this level’s snapshot scheme on a branch as soon as its independently constructed broad list has |D|>D_{\max}. Use a fixed partition of the unit circle into O(D_{\max}/b^6) arcs, rounding each \chi_n(Y), n\in D, to a chosen point of its arc with error at most c_5b^6/D_{\max}. Track the joint finite mass measure of these rounded coordinates and the side information, normalized by \alpha; do not renormalize it to probability. For a finite signed measure \mu on \Omega, write \|\mu\|_{\rm var}=|\mu|(\Omega) for its full variation norm, without a factor 1/2.

After the broad-list step at each jumping post node, refresh the snapshot once if it is the first snapshot, if any coordinate was appended, or if some prospective child pre measure differs in total variation by more than c_6b^6 from the stored post measure of the current snapshot. The last comparison uses the same coordinates on both measures; an append itself causes replacement. On refresh store the present post measure. All these decisions inspect deterministic measures at the current prefix.

At a snapshot take a maximal separated subcollection of the elements of [D] high at threshold b/2 at this node. In the main case separation means that differences avoid

B_A(r_{\rm sep})+\{j\in\mathbb Z:|j||p|\le\sigma\};

in the scalar case it means that differences have magnitude greater than c_7b^4/|p|. Take all the constants c_5,c_6,c_7 sufficiently small.

Lemma 7.4 (Snapshot covers and birth count). The main and scalar snapshot lists have respectively the cardinality bounds in Proposition 7.2. While a snapshot is current, it covers every v\in[D] that is high at threshold 0.9b at a jumping node, with its stored separation neighborhoods measured at the snapshot’s birth interval. Moreover, if N_{\rm snap} is the number of snapshots on the capped uniform output path, then

\mathbb E_{\mathbf U_e}N_{\rm snap} \le 1+D_{\max}+C b^{-6} +C b^{-12}\bigl(d_1+D_{\max}\log(CD_{\max}/b)\bigr). \tag{7.12}

Every snapshot is born before its node’s next fair bit. Its stored representatives may subsequently be used on the same plateau even after replacement of that snapshot.

Proof of the cover assertions. At a stabilized live node, a separated high collection of size \lceil e^{k_{b/2}}\rceil would request another update. The computation levels include b/2. Thus a maximal collection has at most this size and covers all high elements of [D] by the indicated symmetric neighborhoods.

For scalar lists, Parseval applied to a\mathbf1_{p'} on each prospective child gives

\sum_{u\in\mathbb Z}|f_{p'}^u|^2 =\frac{1}{\alpha^2|p'|^2}\int_{p'}a^2 \le\frac{C e^d}{\alpha|p|}.

After rotating by the phase at a point of p', its moments vary by at most C|u-u'||p|. Around each separated representative take an integer neighborhood of radius c_7b^4/(3|p|). These neighborhoods are disjoint; each contains at least c b^4/|p| integers, with the bound also true when that number is less than one. On one of the two children their moments are still at least c b. Summing their squared moments over the two children bounds the list length by C e^d/(\alpha b^6)\le\exp(C(d_1+\log(1/b))). For the final inequality use \alpha\ge e^{-d} on a live plateau and the definition of d_1 on a static plateau.

For the transfer from one node to another, the phase of any signed sum v\in[D] is approximated by the corresponding product of rounded phases with uniform error at most c_5b^6. Scalar moments consequently change by O((c_5+c_6)b^6) between a prospective child measure and the stored post measure if no refresh occurs. For the main functional, if \mu is a joint mass measure of input phase and side value, write M for its side marginal and f for its phase moment. Then

\int\frac{|f|^2}{M} =\sup_{|\gamma(z)|\le1} \int\bigl(2\operatorname{Re}(\gamma(z)\chi_v)-|\gamma(z)|^2\bigr) \,d\mu. \tag{7.13}

The tests have absolute value at most three. Phase rounding contributes at most twice the phase error times the total mass, which is at most four; total variation contributes at most three times the distance. Hence the same O((c_5+c_6)b^6) transfer applies to the main functional. A prospective child high at 0.9b therefore makes the stored post law high with ample slack above b/2. The scalar triangle inequality and the convexity of (M,f)\mapsto |f|^2/M imply that at least one prospective child of that stored post law is high at b/2. Its signed sum belongs to the stored cover. If a fresh snapshot is born at the current node, use its prospective cover directly. This proves the transfer assertion in every case.

Proof of the birth count. Each rounded coordinate is observed only once, when appended, and costs at most C\log(CD_{\max}/b) in conditional entropy. The successive coordinate spaces increase by retaining all previous coordinates. Before cap breach there are at most D_{\max} new coordinates on a path. Thus Lemma 7.1 gives

\mathbb E_{\mathbf U_e}\sum_p \int\frac{|D_{\rm jump}|^2}{M} \le C d_1+C D_{\max}\log(CD_{\max}/b). \tag{7.14}

Here M has the full coordinate list available at that jump. A Fisher term computed after forgetting some coordinates is smaller, by conditional Jensen. This permits comparisons in old coordinate spaces while charging all of them to Equation (7.14).

Start a comparison interval at each snapshot’s post state in its fixed coordinates. If its next replacement is caused by an append, or if the cap is breached, end the old comparison before that node’s jump, in the old coordinates. Intervening removals at that node can be included. If the next replacement is caused by distance with no append, include its node’s jump in the old comparison and end at the selected child pre measure. The new snapshot nevertheless starts at that trigger node’s post measure, before this same bit. A terminal interval ends at the plateau stop. Figure 2 records this order. Inclusion of a trigger bit is predictable, because the trigger is a test of both prospective children before either is selected.

Figure 2: A distance trigger without an append. At the post state p, both child pre measures are compared with the old snapshot M_*. The snapshot is replaced by M_p before r_{j+1} is read. The old comparison includes that bit and ends at the actual child pre measure; the new comparison starts at M_p, so this bit is used twice at most. Any child removals occur afterward. The coefficient of the jump from p is emitted using the delayed sign r_{j+2}, even if no child jump is made. Here c_p denotes that fixed normalized transform coefficient. An append instead ends the old comparison before r_{j+1}, in its old coordinates.

On one comparison interval let k be each removed positive joint mass measure, expressed in its fixed coordinate space. Retain k unchanged after its removal and let K be the sum retained so far. Then

N=M+K

is a positive measure martingale under the uniform output bits. At a kill, M loses precisely the measure added to K. At a spatial bit its two values are N\pm D_{\rm jump}, with |D_{\rm jump}|\le M\le N. In particular the martingale starts from the snapshot measure M_*, and all later measures are absolutely continuous with respect to M_*: a zero coordinate set remains zero under positive child averaging and removals.

Use generalized relative entropy for finite measures,

\mathcal E(N\mid M_*)= \int(n\log n-n+1)\,dM_*,\qquad n=dN/dM_*.

At a spatial bit, writing t=D_{\rm jump}/N, its average increment is

\int \frac N2\bigl((1+t)\log(1+t)+(1-t)\log(1-t)\bigr) \le C\int\frac{|D_{\rm jump}|^2}{N} \le C\int\frac{|D_{\rm jump}|^2}{M}.

The scalar ratio in this inequality is bounded on -1\le t\le1, with its limit at zero understood. Kills leave N unchanged. Bounded stopping in the finite-depth tree therefore bounds the expected terminal entropy by the expected Fisher sum on the comparison interval.

At a distance trigger, one of the two prospective children differs from M_* by more than c_6b^6, so the selected child has this property with conditional probability at least one half. On that event either \|K\|\ge c_6b^6/2, or \|N-M_*\|_{\rm var}>c_6b^6/2. On the complementary event \|K\|<c_6b^6/2, the latter alternative holds, N has mass at most 4+c_6b^6/2, and M_* has mass at most two. The scalar inequality

\frac{(x-1)^2}{1+x}\le C(x\log x-x+1),\qquad x\ge0,

followed by Cauchy–Schwarz gives

\|N-M_*\|_{\rm var}^2 \le C\bigl(N(\Omega)+M_*(\Omega)\bigr)\mathcal E(N\mid M_*),

where \Omega is the fixed observation space of the comparison interval. Thus the latter event forces \mathcal E(N\mid M_*)\ge c b^{12}. Conditional on the comparison start, the probability that its next replacement is a distance trigger is at most

C\mathbb E\left[ b^{-6}\sum_{\rm interval}\xi_p +b^{-12}\sum_{\rm interval} \int\frac{|D_{\rm jump}|^2}{M} \,\middle|\,\text{comparison start}\right]. \tag{7.15}

Indeed the far-child event has at least half the trigger probability, and the preceding deterministic alternatives charge it to removals or to terminal entropy, whose expectation was just bounded.

Sum Equation (7.15) with the probabilities of the starts. Every jump is in at most two comparison intervals: at a distance trigger it is the final bit of the old interval and the first bit of the new one. Appends stop old comparisons before their jump. Removals likewise have bounded overlap. By Equations (7.2) and (7.14), distance triggers contribute the last two terms in Equation (7.12). Initialization contributes at most one, and append replacements at most D_{\max}. This proves the count. No coordinate is re-observed merely because a snapshot is replaced, and the conditional argument did not condition on a realized side or rounded value.

7.6 Activation and the one block of missed jumps

Lemma 7.5 (Activation and localization before activation). On a capped scheme, declare a high u activated at a jump p if there is v\in[D] such that |u-v||p|\le\sigma in the main case, or |u-v||p|\le c_8b^4 in the scalar case. A representative from the current snapshot then gives the fixed decompositions and bounds in Proposition 7.2. All non-bad high jumps before first activation lie in one block of at most H_b depths as stated there.

Proof. Freeze the spatial phase \chi_{u-v} at any point of p. The normalized scalar error on a prospective child is at most C|u-v||p|, by its mass bound. In the main case conditional Cauchy–Schwarz gives

\left|\sqrt{G_{p'}^u}-\sqrt{G_{p'}^v}\right| \le\left(\frac1\alpha\mathbb E_{p'} [a|\chi_{u-v}(Y)-\chi_{u-v}(y_p)|^2]\right)^{1/2} \le C|u-v||p|.

Thus v is high at 0.9b, since \sigma\ll b_0 and the scalar prefactor is small. By Lemma 7.4, it is covered by a stored s. If that snapshot was born at q\supseteq p, its main cover gives v-s=n+j, with \Delta_A(n)\le r_{\rm sep} and |j||q|\le\sigma. Set a'=u-v+j. Then |a'||p|\le2\sigma. The channel A is fixed throughout the plateau, so n retains its metric bound. The scalar cover instead gives |v-s||q|\le c_7b^4, and hence |(u-v)+(v-s)||p|\le(c_7+c_8)b^4. Fix this s,n,a' at activation; all their spatial error bounds only improve on descendants.

Now let p_0 be the first non-bad high jump before any activation, if one exists. By Equation (7.6) there is a signed combination v\in[D] at that node with

|u-v||p_0|\le P(D_{\max},d_1,1/b)

for a fixed polynomial P. The same signed combination remains available as the list grows. After r further depths the left side with the new interval is at most 2^{-r}P(D_{\max},d_1,1/b). For

r\ge C\log(CD_{\max}d_1/b)+C\log(1/\sigma)

this is at most the relevant activation tolerance, \sigma or c_8b^4. Any subsequent high jump then activates, whether or not it is bad. Hence every non-bad high jump strictly before first activation lies in the one depth block beginning at p_0 of the claimed length. If no such p_0 exists there is no missed non-bad high jump. This proves the assertion and completes Proposition 7.2.

The cap has been deterministic throughout the probability estimates. In Section 8 the estimates will first be made simultaneous over all dyadic caps, after which a cap comparable to 1+|D|_{\rm final} may be selected on each path. In particular Equation (7.10) and Jensen give

\mathbb E_{\mathbf U_e}\log(C(1+|D|_{\rm final})) \le C\log(Cd_1/b),

so this later choice incurs a logarithmic expected cost and excludes no large-list paths.

8 Delayed transforms along a plateau path

The representatives from Section 7 will be used on intervals of depths chosen while estimating an actual modulation. We therefore prove interval bounds under explicit amplitude conditions, then make the bounds simultaneous over all snapshot representatives and dyadic list caps. Scalar moments handle small amplitudes and static plateaus; moments that retain the main label handle the main levels on live plateaus.

Fix a plateau e, with normalization \alpha=\alpha_e, and use uniform probability \mathbf U_e on its root. Write \mathcal F_j for the spatial prefix through depth j, together with the deterministic processing data at that prefix. In particular, every post state, its two prospective pre-child states, and the decision to retain its jump are \mathcal F_j-measurable. The sign r_{j+1} is independent and fair conditional on \mathcal F_j. The multiplier \epsilon_j satisfies |\epsilon_j|\le1 and is \mathcal F_{j+1}-measurable. It is the same for every label coordinate. Thus a coefficient formed at depth j has its output multiplied by the later sign r_{j+2}. The additional bit allowed in the finite model supplies this sign even for the final contributing jump.

All moments below are divided by \alpha. At a jumping post state the total mass is at most 2, and at either prospective pre-child it is at most 4. Intervening removals have normalized masses \xi_p; full removal at an exit can be included. We write K_e=\sum_p\xi_p for the sum along the entire plateau path. These quantities and the local budget statements are those of Lemma 7.1. No realized input membership or noisy label is sampled on this uniform output path. All statements use the sufficiently large d convention already fixed; their constants are uniform in the finite spatial depth and in every integer-frequency cutoff.

8.1 The predictable sign inequality

We use the standard exponential-supermartingale proof of a maximal inequality, based on Ville’s stopping argument [23, Chapter V, first section, Theorem 1, pp. 84–85]. The finite-depth version below also records the conditional and predictable-masking forms required later.

Lemma 8.1 (Predictable signs). Let (r_{i+1}) be fair signs in a filtration (\mathcal F_i), and let c_i be real \mathcal F_i-measurable coefficients. For every \theta>0 and t\ge0, outside an event of probability at most e^{-t},

\sum_{i<n}c_i r_{i+1} \le \frac{\theta}{2}\sum_{i<n}c_i^2+\frac{t}{\theta} \qquad\text{for every prefix }n.

The same statement holds conditionally at a stopping-time origin and with predictable stopping or masking of the coefficients. For complex sums one may apply the assertion to both signs of each real component, at a total failure probability at most 4e^{-t}.

Proof. The elementary inequality \cosh x\le e^{x^2/2} shows that

\exp\left(\theta\sum_{i<n}c_i r_{i+1} -\frac{\theta^2}{2}\sum_{i<n}c_i^2\right)

is a nonnegative supermartingale starting at one. Stop at the first crossing of e^t to obtain the asserted maximal inequality. Everything is in finite depth, so optional sampling needs no limiting argument. Starting the same calculation with conditional expectation given the origin proves the conditional assertion. Predictable masks are absorbed in c_i.

8.2 Scalar moments and their masked transforms

Proposition 8.2 (Scalar plateau paths). Let 0<b\le1 be a scalar level, and let s be a representative selected at a predictable snapshot birth on a live or static plateau. At a jumping post state put

f_p=\frac{\mathbb E_p[a\chi_s]}{\alpha},\qquad u_p=\frac{f_{p,+}-f_{p,-}}2,\qquad \eta_p=\mathbf1_{\{\max(|f_{p,+}|,|f_{p,-}|)\le C_3b\}},

where C_3 is a fixed constant and f_{p,\pm}=f_p\pm u_p are prospective pre-child moments. For every t\ge0, with conditional failure probability at most Ce^{-t} at the birth, every contiguous interval I of subsequent plateau jumps satisfies

\left|\sum_{p\in I}\eta_p u_p\epsilon_{j(p)}r_{j(p)+2}\right| \le C b^{1/4}(1+K_e+t). \tag{8.1}

Constants are independent of depth, s, and b; they may depend on the fixed mask constant C_3.

Proof. Identify a complex moment with a point in \mathbb R^2 and set

\psi(x)=b^{1/2}\big((|x|^2+b^2)^{3/4}-b^{3/2}\big).

On |x|\le4 one has 0\le\psi(x)\le Cb^{1/2} and |\nabla\psi(x)|\le Cb^{1/2}. More precisely,

\nabla^2\psi(x) =\frac32 b^{1/2}(|x|^2+b^2)^{-5/4} \left((|x|^2+b^2)\operatorname{Id}-\frac12 x x^{\mathsf T}\right).

Its least eigenvalue is at least \frac34 b^{1/2}(|x|^2+b^2)^{-1/4}. Consequently it is at least c b^{1/2} throughout that ball, and at least c(C_3)>0 on |x|\le C_3b.

Define the nonnegative averaging drift and the signed fluctuation by

\delta_p=\frac{\psi(f_p+u_p)+\psi(f_p-u_p)}2-\psi(f_p),\qquad n_p=\frac{\psi(f_p+u_p)-\psi(f_p-u_p)}2.

Taylor’s formula on the two segments f_p\pm t u_p, 0\le t\le1, gives

\delta_p\ge c b^{1/2}|u_p|^2, \qquad \delta_p\ge c(C_3)\eta_p|u_p|^2, \qquad n_p^2\le Cb|u_p|^2\le Cb^{1/2}\delta_p. \tag{8.2}

For the second inequality the masked segment lies in the ball of radius C_3b, since its endpoints do. The last inequality follows from the gradient bound and the first inequality.

At the actual pre-child the potential is exactly \psi(f_p)+\delta_p+r_{j+1}n_p. A removal of normalized mass \xi changes the complex moment by at most \xi, and hence changes \psi by at most Cb^{1/2}\xi. Telescoping from birth to any prefix end therefore gives

\sum\delta_p\le Cb^{1/2}(1+K_e)-\sum r_{j+1}n_p.

Apply Lemma 8.1 to the negative fluctuation with \theta=\gamma b^{-1/2}, where \gamma>0 is a sufficiently small numerical constant. By Equation (8.2) its quadratic term is at most \frac12\sum\delta_p. Outside probability e^{-t} this proves, on every prefix,

\sum\eta_p|u_p|^2\le Cb^{1/2}(1+K_e+t).

Now \eta_pu_p\epsilon_j is known in \mathcal F_{j+1} before r_{j+2}. Apply Lemma 8.1 to the real and imaginary components of this delayed transform with \theta=b^{-1/4}. Their quadratic sums are bounded by the preceding display. Thus every prefix transform has absolute value at most Cb^{1/4}(1+K_e+t) except on probability Ce^{-t}. Subtracting the two prefixes bounding any contiguous interval proves Equation (8.1). This subtraction retains precisely the fixed mask \eta_p within the interval.

8.3 The energy of label-valued moments

For the main levels we need an L^1(dz) estimate with an absolute coefficient on the Fisher sum. This will let the assembly charge information over disjoint intervals while paying separately for the choice of a representative. The label norm controls the scalar transform because integration against a label character recovers a shifted scalar modulation multiplied by a thermal Fourier coefficient. The main representative is chosen so that the coefficient needed in this recovery is close to one; Section 9 gives the exact decoding and the remaining spatial phase error. We begin with the energy identities for the label-valued moments.

Fix a main representative s on a live plateau, so its main channel A is fixed. At a jumping post state p with post density a_p, define

M(z)=\frac{\mathbb E_p[a_p(Y)K_A(z-Y)]}{\alpha},\qquad f(z)=\frac{\mathbb E_p[a_p(Y)\chi_s(Y)K_A(z-Y)]}{\alpha}.

Here the observation consists only of the main label: the rounded coordinates used to construct snapshots have been integrated out. Define M_\pm,f_\pm by the same formulas on the two prospective children, using the same post density a_p. Write

f=Mv,\qquad f_\pm=M_\pm v_\pm,\qquad |v|,|v_\pm|\le1, \qquad D=\frac{M_+-M_-}2,\qquad \mathfrak a=\frac{f_+-f_-}2.

The functions v are set to zero at zero mass. Define

G=\int M|v|^2\,dz,\qquad J=\int\frac{D^2}{M}\,dz, \qquad G_\pm=\int M_\pm|v_\pm|^2\,dz.

In particular J agrees with Equation (6.4). Quotients with zero denominator are assigned zero; nonnegativity and child averaging make their numerators zero as well.

Lemma 8.3 (Energy drift and innovations). Put

g_z=\frac12\big(M_+|v_+-v|^2+M_-|v_--v|^2\big),\qquad g=\int g_z\,dz,\qquad U_p=\frac{G_+-G_-}2.

Then

\begin{align*} &g_z=\frac{M_+|v_+|^2+M_-|v_-|^2}{2}-M|v|^2\ge0, &|\mathfrak a-vD|^2&\le M g_z,\tag{8.3}\\ &G_{\mathrm{actual\ pre}}=G+g+r_{j+1}U_p, &|U_p|^2&\le C\max(G_+,G_-)(J+g),\tag{8.4}\\ &&\left(\int|\mathfrak a|\,dz\right)^2&\le C(GJ+g). \tag{8.5}\end{align*}

A removal of normalized mass \xi changes G in absolute value by at most 3\xi. From any live jumping post state,

\mathbb E\sum_{\mathrm{future}}g\le C, \qquad \mathbb E\sum_{\mathrm{future}}J\le Cd, \qquad \mathbb E\sum_{\mathrm{future}}\xi\le2, \tag{8.6}

where the conditional law is spatially uniform on that state.

Proof. Child averaging gives M_+(v_+-v)=\mathfrak a-vD and M_-(v_--v)=-(\mathfrak a-vD). Expansion around v proves the first identity. If both child masses are positive, it also gives

g_z=\frac{M|\mathfrak a-vD|^2}{M^2-D^2}.

This implies the innovation inequality. If a child mass is zero, the corresponding moment is zero and \mathfrak a-vD=0, giving the same inequality. Write e=\mathfrak a-vD. Expanding the difference of the two child energies yields

U_p=\int D|v|^2\,dz+2\operatorname{Re}\int\overline v e\,dz+R, \qquad |R|\le g.

Cauchy–Schwarz bounds the first two integrals in absolute value by \sqrt{GJ} and \sqrt{Gg}, respectively. Since G+g=(G_++G_-)/2\le\max(G_+,G_-), squaring proves Equation (8.4). The decomposition \mathfrak a=vD+e gives

\int|\mathfrak a|\,dz\le\sqrt{GJ}+\sqrt{\left(\int M\,dz\right)g} \le\sqrt{GJ}+\sqrt{2g},

which proves Equation (8.5).

For the removal assertion, the energy has the variational expression

\int\frac{|f|^2}{M}\,dz =\sup_{|\gamma(z)|\le1} \int\big(2\operatorname{Re}(\overline\gamma f)-|\gamma|^2M\big)\,dz.

Removing k_M\ge0 and k_f with |k_f|\le k_M changes every test in this supremum by at most 3\int k_M=3\xi. Both the old and new moment-to-mass ratios have modulus at most one, so the same test class represents both energies. Finally telescope Equation (8.4) over the entire future stopped plateau path. Its martingale terms have zero conditional expectation, its terminal pre energy is at most 4, and the removals have expected total at most 2. This gives the first bound in Equation (8.6). The other two are Lemma 7.1.

8.4 A vector estimate on an outer-valid run

For a main dyadic level b\in[b_0,1], declare a jump exceptional when b<1 and it is bad for the main criterion at level 2b. Put

e'_p=C b_0^{-2}e_{p,2b}\mathbf1_{\{p\text{ exceptional}\}} \quad(b<1),\qquad e'_p=0\quad(b=1).

Choose C so that e'_p\ge1 on every exceptional jump; this is possible because a bad main jump has edge mass at least c(2b)^2. Let E_{e,b}=\sum_p e'_p on the whole plateau path. Fix a sufficiently large numerical C_4. For the specified s,b, call a jump inner-valid if it is exceptional or \max(G_+,G_-)\le C_4b^2, and outer-valid if it is exceptional or \max(G_+,G_-)\le2C_4b^2. At b=1 choose C_4 large enough that every jump is inner-valid.

Lemma 8.4 (Drift on a run). Start predictably at an outer-valid jump, and stop before the first outer-invalid jump or at plateau end. For t\ge0, outside conditional probability e^{-t}, every prefix P of this run satisfies

\sum_{p\in P}g_p \le C\left(1+\sum_{p\in P}\xi_p+\sum_{p\in P}e'_p\right) +Cb\sum_{p\in P}J_p+Cb t. \tag{8.7}

Removals here include those between consecutive jumps and may also be enlarged to the whole plateau sum.

Proof. Telescope the exact energy identity of Equation (8.4). Initial and terminal energies are bounded by 4 and removals cost at most 3\xi. On an exceptional jump, |U_p|\le4\le4e'_p, so pay its noise directly. On the remaining jumps, Equation (8.4) gives |U_p|^2\le Cb^2(J_p+g_p). Apply Lemma 8.1 to minus this remaining noise with parameter \theta=\gamma/b. Its quadratic bound is C\gamma b\sum(J_p+g_p). Choose \gamma small enough to absorb the g term, using b\le1. Its other term is Cb t, which proves the simultaneous prefix assertion.

Lemma 8.5 (Delayed norm estimate on a run). Under the hypotheses of Lemma 8.4, set \lambda=b^{1/2} and W_P(z)=\sum_{p\in P}\mathfrak a_p(z)\epsilon_{j(p)}r_{j(p)+2}. With conditional failure probability at most Ce^{-t}, simultaneously for every run prefix,

\int|W_P(z)|\,dz \le C\sum_{p\in P}J_p +C\lambda^{-1}\left(1+\sum_{p\in P}\xi_p+\sum_{p\in P}e'_p\right) +C\lambda t. \tag{8.8}

In particular, the coefficient of the Fisher sum is absolute, independent of b and of the number of representatives.

Proof. We prove the norm estimate using both signs of a two-step block. Separate contributing depths into their two parities. For one parity start at its first depth in the run and group successive spatial jumps into blocks j,j+1, with j contributing to this parity. Only complete blocks lying in the prefix will be used in the norm telescope. Denote the accumulated transform for this parity at the block start by W^{(0)}, and put

X_0=(W^{(0)},f_0,\lambda M_0),\qquad \Phi=|X_0|=\big(|W^{(0)}|^2+|f_0|^2+\lambda^2M_0^2\big)^{1/2}.

The complex coordinates are viewed as pairs of real coordinates. Subscript 0 denotes the first post state. After its sign r=r_{j+1}, write (\bar M,\bar f) for the actual pre-child, and (M_1,f_1) for its post state after removals (k_M,k_f), where

\bar M=M_0+rD_0,\quad \bar f=f_0+r\mathfrak a_0,\quad M_1=\bar M-k_M,\quad f_1=\bar f-k_f, \quad k_M\ge0,\quad |k_f|\le k_M.

Here and in this proof these are pointwise identities in z; in particular \int k_M=\xi_{j+1}. At the end of the second jump, before its child’s processing, the change in X_0 is exactly

\Delta=\big(\mathfrak a_0\epsilon_jr',\;r\mathfrak a_0-k_f+r'\mathfrak a_1, \lambda(rD_0-k_M+r'D_1)\big),\qquad r'=r_{j+2}. \tag{8.9}

We first prove the required pointwise quadratic estimate. Put

q=(|v_0|^2+\lambda^2)^{1/2},\qquad H=M_0q\le\Phi, \qquad j_i=\frac{D_i^2}{M_i},\qquad e_i=\mathfrak a_i-v_iD_i.

For M_0>0, q\ge\lambda and |v_0|^2/q\le1. Positivity of the child masses and of the removed density gives

0\le M_1\le\bar M\le2M_0, \quad 0\le k_M\le\bar M\le2M_0, \quad |D_i|\le M_i, \quad |e_i|^2\le M_i g_{i,z}. \tag{8.10}

There is also a comparison of the child’s post mean with the starting mean:

M_1|v_1-v_0|^2\le4g_{0,z}+8k_M. \tag{8.11}

To verify it, let \bar v=\bar f/\bar M when \bar M>0, and zero otherwise. The actual child’s term in the definition of g_{0,z} gives \bar M|\bar v-v_0|^2\le2g_{0,z}. If M_1=0, the left side of Equation (8.11) is zero. If 0<M_1\le\bar M/2, then M_1\le k_M and the bound |v_1-\bar v|\le2 gives M_1|v_1-\bar v|^2\le4k_M. If M_1>\bar M/2, the identity

M_1(v_1-\bar v)=k_M\bar v-k_f

instead gives M_1|v_1-\bar v|^2\le4k_M^2/M_1\le4k_M, since k_M<M_1. Thus the same bound covers every intermediate or arbitrarily small remaining mass. The squared triangle inequality, and M_1\le\bar M, now prove Equation (8.11).

Splitting \mathfrak a_0=v_0D_0+e_0 and retaining the factor q yields

\frac{|\mathfrak a_0|^2}{H}\le2j_0+\frac{2g_{0,z}}\lambda. \tag{8.12}

For the child jump use the three terms \mathfrak a_1=v_0D_1+(v_1-v_0)D_1+e_1. Their separate squared ratios obey

\begin{align*} \frac{|v_0D_1|^2}{H}&\le2j_1,\\ \frac{|v_1-v_0|^2D_1^2}{H} &\le\frac{2}{\lambda}M_1|v_1-v_0|^2 \le\frac{8g_{0,z}+16k_M}{\lambda},\\ \frac{|e_1|^2}{H}&\le\frac{2g_{1,z}}\lambda.\end{align*}

The factor M_1/M_0\le2 supplies the first and third bounds; the second uses |D_1|\le M_1 and Equation (8.11). In particular no lower bound on M_1/M_0 is used. Hence

\frac{|\mathfrak a_1|^2}{H} \le6j_1+\frac{24g_{0,z}+6g_{1,z}+48k_M}{\lambda}. \tag{8.13}

The mass-coordinate and removal terms satisfy

\frac{\lambda^2D_0^2}{H}\le j_0, \qquad \frac{\lambda^2D_1^2}{H}\le2j_1, \qquad \frac{k_M^2+|k_f|^2}{H}\le\frac{4k_M}{\lambda}.

These follow again from Equation (8.10) and 0<\lambda\le1. If M_0=0, all the masses, half-differences, and removals in Equation (8.9) are zero by positivity, so \Delta=0 whether or not W^{(0)} is zero. This deals with every zero denominator without an approximation argument.

For \Phi>0, concavity of the square root gives the global norm bound

|X_0+\Delta|\le\Phi+\left\langle\frac{X_0}{\Phi},\Delta\right\rangle +\frac{|\Delta|^2}{2\Phi}.

Combining the preceding pointwise inequalities and integrating proves

\int\frac{|\Delta|^2}{2\Phi}\,dz \le C(J_0+J_1)+C\lambda^{-1}(g_0+g_1+\xi_{j+1}). \tag{8.14}

The absolute Fisher coefficient results from the estimates of v_0D_i with q, rather than replacing q by \lambda in those terms.

It remains to control the integrated linear term and its timing. Write \Gamma_0=X_0/\Phi, with value zero where \Phi=0. Its norm is at most one. The linear term is

r L_0+r'L_1+L_{\rm kill},

where L_{\rm kill} has absolute value at most 2\xi_{j+1}, and L_0,L_1 are real integrated coefficients. Specifically L_0 pairs \Gamma_0 with (0,\mathfrak a_0,\lambda D_0), whereas L_1 pairs it with (\mathfrak a_0\epsilon_j,\mathfrak a_1,\lambda D_1). By Equation (8.5) and (\int|D_i|)^2\le(\int M_i)J_i\le2J_i, their squares satisfy

|L_0|^2+|L_1|^2 \le C\sum_{i=0,1}\big((\lambda^2+G_i)J_i+g_i\big). \tag{8.15}

At a block start W^{(0)} is \mathcal F_j-measurable: the previous output of the same parity, if present, used r_j. Thus L_0 is known before r_{j+1}. After that bit the removals, second post state, and \epsilon_j are all known, so L_1 is known before r_{j+2}. Form one scalar sign scheme by including rL_0 whenever the first jump is in the run, and including r'L_1 only if the second jump is also in the run. Both inclusion decisions are known before their respective signs. The linear sum through the last complete block of any prefix is itself a prefix of this scheme. If the run stops after a first step, no first coefficient is retrospectively conditioned on the next validity decision.

On a nonexceptional included state G_i\le2C_4b^2. At an exceptional state both G_i and J_i are at most 2, and e'_i\ge1. Therefore

G_iJ_i\le Cb^2J_i+Ce'_i.

Apply Lemma 8.1 to the linear scheme with parameter 1/\lambda. Its quadratic term, by Equation (8.15), is at most

C\sum J_p+C\lambda^{-1}\sum(g_p+e'_p),

because \lambda=b^{1/2} and b\le1. The tail term is \lambda t. All sums here can be restricted to the complete blocks of the prefix. At the first block W^{(0)}=0 and \int\Phi\le C. Between complete blocks, passing from a pre state to the next post state changes the integrated norm by at most 2\xi, by the norm’s Lipschitz property. Telescoping, using Equation (8.14), proves for the completed blocks the bound

C\sum J_p+C\lambda^{-1}\left(1+\sum\xi_p+\sum e'_p+\sum g_p\right) +C\lambda t.

At most one contributing output of this parity is unmatched at a prefix end. It costs at most \int|\mathfrak a_p|\le\int M_p\le2, even if its delayed sign occurs after the run has stopped. Add this output directly. The other parity is started conditionally on reaching its first jump, one step later if necessary, and obeys the same bound with the same conditional failure estimate. Each prefix jump is used in at most two of the completed-block estimates. Finally apply Equation (8.7) on this prefix. Its g contribution becomes C\lambda^{-1}(1+\sum\xi+\sum e')+C(b/\lambda)\sum J+C(b/\lambda)t. Since b/\lambda=\lambda\le1, adding the two parity bounds gives Equation (8.8).

8.5 Predictable runs and Fisher checkpoints

The run estimate charges the Fisher sum from its origin. Subtracting two prefixes to estimate a later interval would therefore charge the earlier part of the run twice. We introduce predictable checkpoints so that this extra Fisher charge is bounded by an absolute constant.

Proposition 8.6 (Vector estimates on arbitrary inner-valid intervals). Fix a live plateau, a main level b\in[b_0,1], and a representative s at its snapshot birth. There is a family of predictable estimate origins, containing run starts and the checkpoints described below, whose conditional expected size at the birth is at most P(d,b_0^{-1}) for an absolute fixed polynomial P. For every t\ge0, outside a conditional event of probability at most C P(d,b_0^{-1})e^{-t}, every contiguous interval I of inner-valid jumps after this birth satisfies

\int\left|\sum_{p\in I}\mathfrak a_p(z)\epsilon_{j(p)}r_{j(p)+2}\right|\,dz \le C\sum_{p\in I}J_p +Cb^{-1/2}(1+K_e+E_{e,b})+Cb^{1/2}t. \tag{8.16}

The Fisher sum is taken only over I. The representative remains eligible for these estimates until the plateau ends, even if another snapshot has meanwhile replaced the one at which it was born.

Proof. Start a run at the first inner-valid opportunity after birth, possibly at the birth jump itself. Follow consecutive outer-valid jumps, stopping before the first outer-invalid jump or at plateau end. After a break, start the next run at the first later inner-valid opportunity. These are predictable rules: both prospective child energies and the exceptional indicator are known before the jump sign. Every contiguous inner-valid interval lies in one of these outer-valid runs.

We first bound the number of run starts conditionally on the birth prefix. A start on an exceptional jump can be charged to that jump. Lemma 7.3 bounds the expected number of these by d\poly(h); there are none when b=1. At a nonexceptional start, G_{\rm start}\le C_4b^2. If its run breaks, include the first invalid jump in a comparison interval, although the transform estimate stops before it. At that jump \max(G_+,G_-)>2C_4b^2. Conditional on the prefix there, with probability at least 1/2 its actual pre-child energy exceeds 2C_4b^2. Mark such a break. The comparison intervals for different runs are disjoint, including their invalid jumps.

Let Q denote one such start-to-break comparison interval, and put \mathcal M_Q=\sum_{p\in Q}r_{j(p)+1}U_p. By the exact energy identity and the removal bound, a marked break implies

C_4b^2\le\sum_{p\in Q}g_p+3\sum_{p\in Q}\xi_p+|\mathcal M_Q|.

Consequently, with a numerical constant depending only on C_4,

\mathbf1_{\{Q\text{ marked}\}} \le Cb^{-2}\sum_{p\in Q}(g_p+\xi_p)+Cb^{-4}|\mathcal M_Q|^2. \tag{8.17}

To use orthogonality, first include every start-to-break interval of this type and the possible unfinished terminal interval, without conditioning on a mark. Its starting rule and inclusion through the invalid bit are predictable. Thus stopped martingale orthogonality gives

\mathbb E\sum_Q|\mathcal M_Q|^2 =\mathbb E\sum_{p\in\bigcup Q}U_p^2 \le C\mathbb E\sum_p(J_p+g_p),

using Equation (8.4) and \max(G_+,G_-)\le4. All expectations here are conditional at birth. Each nonexceptional break is marked with conditional probability at least 1/2, so Equations (8.17) and (8.6) imply

\mathbb E\#\{\text{nonexceptional-start runs that break}\} \le Cb^{-2}+Cdb^{-4}.

At most one final run fails to break. Including exceptional starts proves a bound P(d,b_0^{-1}) for all starts, uniformly in s. The factor \poly(h) is absorbed in a fixed power of b_0^{-1}.

Within each outer-valid run, also start an estimate whenever the accumulated J since the preceding origin reaches one. More precisely, after including the jump that reaches this threshold, reset the counter and use the next jump as a new origin if it remains in the run. Each such origin’s estimate continues to that run’s same break or end. The counter decision is known at this next origin; no future bit is used. The number of additional origins is at most the total J of the path, and J_p\le\int M_p\le2 bounds a threshold overshoot. Therefore their conditional expected number is at most Cd by Equation (8.6).

Apply Lemma 8.5 conditionally at every origin, with the same parameter t. Summing its conditional failure probabilities weights each origin by its probability of being reached. The expected origin count just proved bounds the resulting probability by C P(d,b_0^{-1})e^{-t}. On the event of simultaneous success, fix any desired interval I and choose the latest origin q in its outer-valid run at or before its first jump. The total J between q and that first jump is at most 3; indeed the counter is below one before a threshold crossing, and one jump adds at most two. The interval transform is the difference of two prefixes from q. Their Fisher sums add to

\sum_{p\in I}J_p+2\sum_{q\le p<\min I}J_p \le\sum_{p\in I}J_p+6.

Their removal and exceptional sums are bounded by twice K_e and E_{e,b}, respectively. The extra constant is absorbed by b^{-1/2}(1+K_e+E_{e,b}). The triangle inequality in L^1(dz) then proves Equation (8.16). This argument also permits the interval to begin long after the snapshot’s replacement.

8.6 Simultaneous representatives and random list caps

Lemma 8.7 (Parameters for all representatives at a level). On each plateau and at each permitted level there are nonnegative random parameters t_b,H_b, independent of the actual frequency later being estimated, with the following properties. Almost surely the scalar bound of Proposition 8.2, or the vector interval bound of Proposition 8.6, holds with t=t_b for every representative in the snapshot construction used on the path. The activation conclusion of Lemma 7.5 holds with miss-block length H_b. Moreover

\begin{align*} \mathbb E_{\mathbf U_e}t_b &\le Ck_{b/4}+C\log(Cd/b) &&\text{for main levels}, \tag{8.18}\\ \mathbb E_{\mathbf U_e}t_b &\le C\big(d_1+\log(1/b)\big) &&\text{for scalar levels}, \tag{8.19}\\ \mathbb E_{\mathbf U_e}H_b &\le C\log(Cd_1/b)+C\log(1/\sigma) &&\text{for either criterion}. \tag{8.20}\end{align*}

These assertions hold simultaneously over the countably many levels.

Proof. Begin with a deterministic dyadic cap D_{\max}=2^m, m\ge0. The broad list is constructed independently of this cap. The capped snapshot scheme creates snapshots only while its list length is at most D_{\max}; estimates for representatives already born continue to the plateau end. By Lemma 7.4, the expected number of births is bounded by a fixed polynomial in (D_{\max},d_1,1/b). The number of representatives in each snapshot is at most

R_b^{\rm main}=\lceil e^{k_{b/2}}\rceil, \qquad R_b^{\rm sc}=\exp\big(C(d_1+\log(1/b))\big),

as established in that lemma. Scalar estimates originate at these births. For each main representative, Proposition 8.6 bounds the conditional expected number of run and checkpoint origins by a fixed polynomial in (d,b_0^{-1}). For sufficiently large d this factor is also bounded by a fixed power of d, because b_0^{-1}\le2h^{16}\le d.

Here is the conditional union argument, which also applies when the number of origins is random. Enumerate predictable origins by their finite tree addresses and, within an address, by their representative indices. If B_o is the event that origin o is reached and an estimate starting there fails with deterministic parameter t, then

\mathbb P(B_o) \le Ce^{-t}\mathbb P(o\text{ is reached}).

Sum this inequality over origins and then use their expected count. Thus the failure probability for a capped scheme is at most Ce^{-t} times its stated representative cardinality and polynomial count. No union over the signed span of the broad list is taken.

Choose sufficiently large fixed constants and define baseline parameters

T_b^{\rm main}=Ck_{b/4}+C\log(Cd/b),\qquad T_b^{\rm sc}=C(d_1+\log(1/b)).

They absorb the logarithms of the representative cardinalities and all polynomial factors except the cap. For a sufficiently large further fixed C_{\rm cap}, use

t=T_b+C_{\rm cap}\log(CD_{\max})+z

in each deterministic capped scheme. Its failure probability is at most C2^{-2m}e^{-z}, after increasing C_{\rm cap} if necessary. Summing over all dyadic caps gives Ce^{-z}. These events are nested as z increases, since every right-hand side in the path estimates increases. Let Z_b be the least nonnegative integer z for which all these capped schemes succeed. Integer tail summation gives

\mathbb P(Z_b>n)\le Ce^{-n},\qquad \mathbb E Z_b\le C,

and Z_b is finite almost surely.

Only now choose a cap from the realized output path: let \widehat D be the least dyadic integer at least 1+|D|_{\rm final} for the broad list at this level. Thus 1+|D|_{\rm final}\le\widehat D<2(1+|D|_{\rm final}), and its capped scheme never breaches the cap on this path. The simultaneous assertion over deterministic caps permits this selection without a stopping-time claim about \widehat D. Set

t_b=T_b+C_{\rm cap}\log(C\widehat D)+Z_b.

Broad capture gives \mathbb E|D|_{\rm final}\le Cd_1b^{-C}, so Jensen’s inequality yields

\mathbb E\log(C\widehat D) \le C+\log\big(C(1+\mathbb E|D|_{\rm final})\big) \le C\log(Cd_1/b).

This proves Equations (8.18) and (8.19), using d_1=d for main levels and \log d_1\le d_1 for scalar levels. Use Lemma 7.5 with this same cap and put

H_b=C\log(C\widehat Dd_1/b)+C\log(1/\sigma).

Its expectation satisfies Equation (8.20). In the static case the optional term \log(1/\sigma)=40\log h is at most Cd\le Cd_1; thus it causes no loss in the scalar plateau bound used below.

For each level the construction has a probability-one success event. Their countable intersection still has probability one. The resulting parameters depend on the plateau path and the level, but not on the integer frequency whose transform will be compared to these representatives.

9 Assembly of the dyadic estimate

We prove the modulated delayed Haar estimate stated in Proposition 6.1. First we bound the contribution of a single plateau, using its representatives to recover each actual modulation. We then sum these contributions with the plateau starting weights, keeping the live Fisher and edge charges global.

Proof of Proposition 6.1. The zero input is immediate. For positive input, apply the finite procedure of Section 6, with the common stopping rule in the statement. In particular, the cutoff always refers to the original \rho, including when a subpiece is being processed. The parameter constants are fixed in the order specified there, and d_0 is increased to satisfy their finitely many large-d requirements.

The contribution of one plateau

For a plateau e, write p_e for its root, \alpha_e for its normalization, and \mathbf U_e for uniform probability on p_e. On an output path x\in p_e, let \mathcal P_e(x) be its consecutive jumping nodes belonging to this plateau. The post density at a node p\in\mathcal P_e(x) is denoted by a_p. Thus all removals at p have already occurred, and a_p is the density whose two prospective child moments are used in that jump. In sums over nodes p, write j=\operatorname{depth}(p). Define

c_p^u=\frac1{\alpha_e}\E_p[a_p(Y)\chi_u(Y)r_{j+1}(Y)], \qquad \mathcal T_e^u(x)=\sum_{p\in\mathcal P_e(x)} c_p^u\epsilon_j(x)r_{j+2}(x).

At each jumping node the post pieces partition the original density. Consequently their transforms satisfy the exact identity

\mathcal H_{I,J}^{u}\rho(x) =\sum_{e:\,x\in p_e}\alpha_e\mathcal T_e^u(x).

Put w_e=\mathbf U(p_e)\alpha_e. The triangle inequality and integration with respect to \mathbf U give

\int_I\max_{\substack{u\in\mathbb Z\\|u|\le Q}} |\mathcal H_{I,J}^{u}\rho|\,d\mathbf U \le\sum_e w_e\mathbb E_{\mathbf U_e} \sup_{u\in\mathbb Z}|\mathcal T_e^u|. \tag{9.1}

We will bound this weighted sum by Cd, charging the live Fisher and edge terms across all pieces and the remaining terms one plateau at a time.

The local mass bound gives |c_p^u|\le2 and prospective child masses at most 4, in this normalization. Let K_e(x) denote the sum of all normalized removal masses \xi_p on the plateau’s processing path, including the intervening and terminal removals allowed in Proposition 7.2; then

\E_{\mathbf U_e}K_e\le1. \tag{9.2}

For live plateaus set d_1=d, and for static plateaus set d_1=d+2+\log^-\alpha_e.

Use the main levels \mathcal B_{\rm m}=\{2^{-j}:b_0\le2^{-j}\le1\} on live plateaus. Use scalar levels \mathcal B_{\rm s}(B)=\{B2^{-j}:j\ge0\}, where B=b_0 on a live plateau and B=1 on a static plateau. At the shared level b_0 the main and scalar constructions are distinct. By Lemma 8.7, they have nonnegative, frequency-independent path parameters t_b^{\rm m},t_b^{\rm s} and miss-block lengths H_b^{\rm m},H_b^{\rm s}, satisfying

\begin{align*} \E_{\mathbf U_e}t_b^{\rm m} &\le C k_{b/4}+C\log(Cd/b), &&b\in\mathcal B_{\rm m},\tag{9.3}\\ \E_{\mathbf U_e}t_b^{\rm s} &\le C(d_1+\log(1/b)), &&b\in\mathcal B_{\rm s}(B),\tag{9.4}\\ \E_{\mathbf U_e}H_b^{\rm m},\quad \E_{\mathbf U_e}H_b^{\rm s} &\le C\log(Cd_1/b)+C\log(1/\sigma) &&\text{for the levels where they are used.} \tag{9.5}\end{align*}

These statements already include all snapshot births, representative choices, starts, checkpoints, and dyadic list caps. In particular we may use the same parameters for every actual integer u.

At a live jumping piece write

E_{p,b}^{\rm m}=e_{p,b}\mathbf1_{\{p\text{ bad at main level }b\}}.

For the scalar construction let \bar b(p) be its largest bad level, if there is one, and put

E_p^{\rm s}=e_{p,\bar b(p)} \quad\text{if a bad scalar level exists},\qquad E_p^{\rm s}=0 \quad\text{otherwise}.

A nonempty subset of the dyadic levels bounded above by B has a largest member. This choice is determined by the jumping piece before its bit and is independent of u. The scalar edge regions are nested as b increases: their defining l_b decreases. Hence

b\le\bar b(p)\quad\Longrightarrow\quad e_{p,b}\le e_{p,\bar b(p)}. \tag{9.6}

Main amplitudes and exact decoding

Fix an integer u and a live plateau path. Starting with its first main activation, including that jump, retain the largest main level activated so far. This level is nondecreasing. Its positive intervals therefore use each level b at most once. On the interval assigned to b, keep the representative s and decomposition supplied at its first activation:

u=s+n+a',\qquad \Delta_A(n)\le r_{\rm sep},\qquad |a'|\,|p_{\rm act}|\le C\sigma. \tag{9.7}

The form A is fixed on this plateau. All subsequent nodes in the interval lie below p_{\rm act}.

For clarity, the comparison of criteria uses no loss depending on the number of representatives. For a prospective child q of p, write f_q^v(z)=\alpha_e^{-1}\E_q[a_p(Y)\chi_v(Y)K_A(z-Y)] and M_q(z)=\alpha_e^{-1}\E_q[a_p(Y)K_A(z-Y)]. For the fixed output point x\in p, conditional Cauchy–Schwarz gives

\left\|f_q^u-\chi_n(\cdot)\chi_{a'}(x)f_q^s \right\|_{L^2(M_q^{-1}dz)} \le C\Delta_A(n)+C|a'|\,|p|.

Indeed the square of the channel part is at most

\frac1{\alpha_e}\E_q\left[ a_p(Y)\int |\chi_n(Y)-\chi_n(z)|^2K_A(z-Y)\,dz\right] =\frac{\E_q a_p}{\alpha_e}\Delta_A(n)^2,

and the spatial part uses |\chi_{a'}(Y)-\chi_{a'}(x)|\le2\pi|a'|\,|p| and \E_q a_p/\alpha_e\le4. The triangle inequality in this weighted Hilbert space proves

\left|\sqrt{G_q^u}-\sqrt{G_q^s}\right| \le C(r_{\rm sep}+\sigma)\ll b_0\le b. \tag{9.8}

If b<1, level 2b has not activated at any jump of this interval. Its non-bad high jumps before activation lie in one block of at most H_{2b}^{\rm m} depths, by Proposition 7.2. Remove the intersection of that entire block with the current interval and pay C H_{2b}^{\rm m} directly, using |c_p^u|\le2. The remaining jumps form at most two contiguous intervals. On each, a nonexceptional jump is not high for u at 2b; thus Equation (9.8) makes it inner-valid for s,b in Proposition 8.6. Here the exceptional jumps are precisely the bad main jumps at 2b. When b=1, every jump is inner-valid after the fixed constant in that definition is chosen large enough, since all G_q^s\le4; no block is removed and no exception is needed.

We next verify the decoding of the vector estimate back to the actual scalar modulation. Put

\mathfrak a_p^s(z) =\frac1{\alpha_e}\E_p[a_p(Y)\chi_s(Y)r_{j+1}(Y)K_A(z-Y)], \qquad c_A(n)=\int\chi_n(\theta)K_A(\theta)\,d\theta.

Evenness and positivity of the thermal Fourier coefficients imply

c_A(n)=1-\tfrac12\Delta_A(n)^2\ge\tfrac12, \qquad \int\chi_n(z)\mathfrak a_p^s(z)\,dz=c_A(n)c_p^{s+n}. \tag{9.9}

This is an exact identity, obtained by putting z=Y+\theta in the inner integral. On any subset of the nodes below p_{\rm act},

\sum_p|c_p^u-\chi_{a'}(x)c_p^{s+n}| \le C|a'|\sum_p|p| \le C|a'|\,|p_{\rm act}|\le C\sigma.

The sum of lengths is geometric because the nodes lie on one descending path. The bound applies in particular to each retained subinterval. Since the multipliers are common to all label coordinates, Equation (9.9) now gives

\left|\sum_{p\in L}c_p^u\epsilon_jr_{j+2}\right| \le2\int\left|\sum_{p\in L}\mathfrak a_p^s(z) \epsilon_jr_{j+2}\right|\,dz+C\sigma \tag{9.10}

for every retained subinterval L.

Apply Proposition 8.6 to L. Its Fisher sum is charged only on L. The intervals at distinct retained levels are disjoint, and there are at most two retained subintervals at each level. Thus these Fisher charges have an absolute bounded multiplicity at every jumping piece. All kill and exceptional sums may use the whole plateau path. Writing the exceptional proxy as e'_p=C b_0^{-2}E_{p,2b}^{\rm m}, and using b^{-1/2}b_0^{-2}\le b_0^{-3} for b\ge b_0, we have bounded the whole main part by a constant times

\sum_pJ_p+b_0^{-3}\sum_p\sum_{b\in\mathcal B_{\rm m}}E_{p,b}^{\rm m} +\sum_{b\in\mathcal B_{\rm m}} \left[b^{-1/2}(1+K_e)+b^{1/2}t_b^{\rm m}+H_b^{\rm m}\right]. \tag{9.11}

The constants in the vector estimates and the geometric decoding errors are absorbed by the displayed b^{-1/2} terms.

The scalar prefix and all static paths

On a live plateau consider the prefix strictly before the first main activation, or the whole path if no main activation occurs. On a static plateau consider the entire path. Use the scalar levels through B specified above, again retaining the largest activated level so far, including activations at the current jump, and using zero if none has activated. There are finitely many changes on this finite path even though the available scalar levels extend arbitrarily close to zero. Each positive level b is retained on at most one contiguous interval.

For its fixed representative s, scalar activation supplies u=s+a' with |a'|\,|p_{\rm act}|\le Cb^4, with the sufficiently small constant in Proposition 7.2. On this level interval sum first over the fixed scalar mask

\max\{|f_{p,+}^s|,|f_{p,-}^s|\}\le C_3b.

This is precisely the masked transform in Proposition 8.2; that proposition applies to the masked sum over the whole contiguous level interval. Freezing \chi_{a'} at x changes this masked sum by at most Cb^4, by the same geometric sum of interval lengths as above. Hence all such masked sums cost at most

C\sum_{b\in\mathcal B_{\rm s}(B)}b^{1/4} (1+K_e+t_b^{\rm s}). \tag{9.12}

Treat every remaining jump directly. This includes all jumps while the retained level is zero. Put A'_p=\max\{|f_{p,+}^u|,|f_{p,-}^u|\}\le4 in the unlabelled scalar notation. The direct jump costs at most A'_p, and a zero value needs no charge. If a positive level b is currently retained, the jump is outside its mask; the spatial phase comparison gives A'_p\ge C_3b-Cb^4>4b, on taking C_3 large enough.

For a live prefix, any direct jump with A'_p>4b_0 is high at main level b_0. In fact f_{p,\pm}^u=\int f_{p,\pm}^u(z)\,dz and Cauchy–Schwarz gives |f_{p,\pm}^u|^2\le(\int M_{p,\pm})G_{p,\pm}^u\le4G_{p,\pm}^u. No main level has yet activated in this prefix. All its non-bad main high jumps at b_0 are therefore paid by its one miss block, and its bad jumps are paid by C b_0^{-2}E_{p,b_0}^{\rm m}, since each direct coefficient is at most 4. These costs are already included, up to constants, in Equation (9.11).

At every other nonzero direct jump we have A'_p\le4B; this holds automatically on a static path with B=1. Choose the largest available scalar level \widetilde b\le A'_p. Then

A'_p/4\le\widetilde b\le A'_p,

and the jump is high at \widetilde b. This level has not activated: if the current level is positive, A'_p>4b implies \widetilde b>b, and otherwise none has activated. All non-bad jumps assigned to a fixed \widetilde b consequently lie in its single preactivation miss block and cost at most C\widetilde b H_{\widetilde b}^{\rm s} together. If the assigned level is bad, then e_{p,\widetilde b}\ge c\widetilde b and Equation (9.6) gives

A'_p\le4\widetilde b\le C e_{p,\widetilde b}\le C E_p^{\rm s}.

Thus only one predictably selected scalar edge level per jump is charged. This establishes the scalar accounting for every positive amplitude, including arbitrarily small amplitudes.

Combining these conclusions, the same frequency-independent quantity bounds \sup_u|\mathcal T_e^u| on a live plateau:

\begin{align*} C^{-1}\sup_{u\in\Z}|\mathcal T_e^u| \le{}&\sum_{p\in\mathcal P_e(x)} \left(J_p+b_0^{-3}\sum_{b\in\mathcal B_{\rm m}}E_{p,b}^{\rm m} +E_p^{\rm s}\right)\\ &+\sum_{b\in\mathcal B_{\rm m}} \left[b^{-1/2}(1+K_e)+b^{1/2}t_b^{\rm m}+H_b^{\rm m}\right]\\ &+\sum_{b\in\mathcal B_{\rm s}(b_0)} \left[b^{1/4}(1+K_e+t_b^{\rm s})+bH_b^{\rm s}\right]. \tag{9.13}\end{align*}

For a static plateau the corresponding bound is

\sup_{u\in\Z}|\mathcal T_e^u| \le C\sum_pE_p^{\rm s} +C\sum_{b\in\mathcal B_{\rm s}(1)} \left[b^{1/4}(1+K_e+t_b^{\rm s})+bH_b^{\rm s}\right]. \tag{9.14}

The path constructions and bounds used here hold simultaneously outside one \mathbf U_e-null set. The countable scalar-level sums have finite expectation by the estimates below. In particular no exceptional set depending on the actual modulation is introduced in these majorants.

Summing levels and live weights

Write b_0=2^{-M}, so M=O(1+\log h). Since k_b=A_0(d/h)(1+\log(1/b))^2, the main cardinality tails satisfy

\sum_{b\in\mathcal B_{\rm m}}b^{1/2}k_{b/4} =A_0\frac dh\sum_{j=0}^{M}2^{-j/2} \bigl(1+(j+2)\log2\bigr)^2 \le C\frac dh. \tag{9.15}

The additional logarithms in Equation (9.3) sum to O(\log(Cd)) with the same b^{1/2} weights. Also,

\sum_{b\in\mathcal B_{\rm m}}b^{-1/2}\E(1+K_e) \le Cb_0^{-1/2}\le Ch^8, \qquad \sum_{b\in\mathcal B_{\rm m}}\E H_b^{\rm m} \le C(1+\log h)\bigl(\log(Cd)+\log h\bigr).

Every fixed polynomial in h, multiplied by 1+\log d, is o(d/h) because h=\log\log d. Thus the second line of Equation (9.13) has expectation O(d/h).

For the scalar series write b=B2^{-j}. The convergent sums \sum_{j\ge0}2^{-j/4} and \sum_{j\ge0}j2^{-j/4}, together with Equations (9.2), (9.4) and (9.5), give

\begin{align*} &\sum_{b\in\mathcal B_{\rm s}(B)} \E\left[b^{1/4}(1+K_e+t_b^{\rm s})+bH_b^{\rm s}\right]\\ &\hspace{1cm}\le C B^{1/4}\bigl(d_1+1+\log(1/B)+\log(1/\sigma)\bigr). \tag{9.16}\end{align*}

For example the miss-block part alone is bounded by CB(\log(Cd_1/B)+\log(1/\sigma)+1), which is no larger than the displayed right side. On live plateaus B=b_0\le h^{-16} and d_1=d, so this is O(d/h^4) and hence O(d/h). On static plateaus B=1 and d_1\ge d, so \log(1/\sigma)=40\log h\le Cd_1 makes it O(d_1). The local scalar edge bound in Proposition 7.2 therefore proves

\E_{\mathbf U_e}\sup_u|\mathcal T_e^u|\le C(d+2+\log^-\alpha_e) \quad\text{for every static plateau.} \tag{9.17}

For every nonnegative jump quantity V_p of a plateau, change of the uniform root measure gives exactly

\sum_e w_e\E_{\mathbf U_e}\sum_{p\in\mathcal P_e(x)}V_p =\sum_{\text{jumping pieces }p}\mathbf U(p)\alpha_e V_p. \tag{9.18}

The Fisher bound in Proposition 6.2 bounds the right side by Cd when V_p=J_p and e is live. The global scalar edge assertion of Proposition 7.2 does the same for V_p=E_p^{\rm s}. Its global main edge assertion yields

b_0^{-3}\sum_{\text{live jumping pieces }p} \mathbf U(p)\alpha_e\sum_{b\in\mathcal B_{\rm m}}E_{p,b}^{\rm m} \le \poly(h)(h+\log d)=o(d).

Finally \sum_{e\text{ live}}w_e\le Ch by Proposition 6.2. Multiplying the O(d/h) per-plateau costs in Equation (9.13) by these weights therefore proves

\sum_{e\text{ live}}w_e\E_{\mathbf U_e}\sup_u|\mathcal T_e^u|\le Cd. \tag{9.19}

Static roots, small mass, and rare transfers

Consider static handling starting with a positive piece of mean m at an interval q. Its later plateaus restart only when their incoming mean exceeds twice the previous normalization. Along any input path each new normalization is therefore more than twice its predecessor, and every normalization is at least m. A plateau with a jump has normalization at most e^{4h}; allowing a final non-jumping overshoot changes this bound only by a factor of two. Consequently the number N_{\rm st} of static starts on the path obeys

N_{\rm st}\le C(1+h+\log^-m),\qquad d+2+\log^-\alpha_e\le d+2+\log^-m. \tag{9.20}

Starts with no subsequent jump may also simply be omitted.

Let \mathbf P_q be the probability with density equal to this entry piece divided by m relative to uniform probability on q. Since the piece is never split again, the probability of reaching a later static plateau root p_e is \mathbf U_q(p_e)\alpha_e/m. Thus the sum of its starting weights is the original weight \mathbf U(q)m multiplied by \E_{\mathbf P_q}N_{\rm st}. Applying Equations (9.17) and (9.20), the entire static descendant family costs at most

C\mathbf U(q)m(1+h+\log^-m)(d+2+\log^-m). \tag{9.21}

This uses the input path only to sum weights; each individual transform estimate remains an integral over its uniform output path.

If 0<m_{\rm top}<1/h, the procedure takes the top static immediately. Writing \ell=\log(1/m_{\rm top}), its total cost is bounded by

Cm_{\rm top}(1+h+\ell)(d+2+\ell)\le Cd.

For an explicit verification, m_{\rm top}h\le1, \sup_{0<t\le1}t\log(1/t)=1/e, and \sup_{0<t\le1}t\log^2(1/t)=4/e^2. Expanding the product and using h\le d proves the inequality.

It remains to sum static families transferred out of an initially live top, so m_{\rm top}\ge1/h. Every true input path undergoes at most one such transfer. Let A_{\rm tr} be the transfer event under its input probability \mathbf P, and let X_* be the maximum negative logarithm of the mean through transfer or stopping, including the transfer state. Proposition 6.2 gives

\mathbf P(A_{\rm tr})\le Ch^{-10},\qquad \|X_*\|_{L^2(\mathbf P)}\le Ch. \tag{9.22}

This includes both charge/count transfers and the small-mass flag. In particular the small-mass event has X_*>d, hence probability at most Ch^2/d^2\le Ch^{-10}. The cap accounting and the evidence argument in that proposition bound the probability of the other transfers; there is no further sequence of live restarts after transfer.

For a transfer piece of mean m, its true reaching probability is \mathbf U(q)m/m_{\rm top} and \log^-m\le X_* on that event. Summing Equation (9.21) over all transfer roots therefore bounds their entire cost by

Cm_{\rm top}\E_{\mathbf P} [\mathbf1_{A_{\rm tr}}(1+h+X_*)(d+2+X_*)].

Only the second moment in Equation (9.22) is needed. Indeed, with \delta=\mathbf P(A_{\rm tr}),

\begin{align*} &\E[\mathbf1_{A_{\rm tr}}(1+h+X_*)(d+2+X_*)]\\ &\quad=(1+h)(d+2)\delta +(d+3+h)\E[\mathbf1_{A_{\rm tr}}X_*] +\E[\mathbf1_{A_{\rm tr}}X_*^2]\\ &\quad\le Cdh^{-9}+C(d+h)h^{-4}+Ch^2\le Cd.\end{align*}

Here Cauchy–Schwarz gives \E[\mathbf1_{A_{\rm tr}}X_*]\le\sqrt\delta\,\|X_*\|_2\le Ch^{-4}, whereas the quadratic term uses its full bound \E X_*^2\le Ch^2. Multiplying by m_{\rm top}\le2 preserves Cd.

Completion

The live estimate Equation (9.19) and the static-family estimates bound the right side of Equation (9.1) by Cd. This proves Equation (6.1). Every construction was in fixed finite depth and every estimate was independent of Q. Increasing Q and applying monotone convergence proves the assertion for all integer modulations as well.

10 Fourier sums and the endpoint summation

We now return to periodic functions. In this section the circle has period one and measure dy, as in Section 2. The Fourier coefficients use \chi_k(y)=e^{2\pi iky}; changing coordinates at the end recovers precisely the normalization of Theorem 1.1.

10.1 A good-set integral bound

For a nonnegative periodic function \rho, let

M\rho(x)=\sup_{\substack{I\ni x\\0<|I|\le1}} \frac1{|I|}\int_I\rho(y)\,dy

be its uncentered interval maximal function, with intervals interpreted on the periodic extension. The elementary interval covering argument gives

\bigl|\{M\rho>\lambda\}\bigr|\le \frac{C}{\lambda}\int_0^1\rho(y)\,dy. \tag{10.1}

Indeed, for a compact subset of the level set, choose finitely many witnessing intervals and a disjoint subfamily whose intervals enlarged by a fixed factor cover that subset. The sum of the lengths of the disjoint intervals is at most \lambda^{-1}\int\rho, up to an absolute periodic overlap constant. Inner regularity proves the stated inequality. In particular, for a mass-one input and h=\log\log d,

G_\rho=\{M\rho\le e^{4h}\}, \qquad |G_\rho^c|\le Ce^{-4h}=C(\log d)^{-4}.

Theorem 10.1 (Truncated integral estimate). There are absolute constants C and d_* such that, whenever d\ge d_* and \rho is a nonnegative periodic function with \int_0^1\rho=1 and \rho\le e^d almost everywhere,

\int_{G_\rho}\sup_{N\ge0}|S_N\rho(x)|\,dx\le Cd. \tag{10.2}

Proof. We first pass from Proposition 6.1 to modulated convolutions, then identify the periodic Hilbert transform.

The two Haar profiles. Extend the following functions by zero outside [0,1):

\mathfrak h =\mathbf1_{[0,1/2)}-\mathbf1_{[1/2,1)},\qquad \mathfrak g =\mathbf1_{[0,1/4)}-\mathbf1_{[1/4,1/2)} -\mathbf1_{[1/2,3/4)}+\mathbf1_{[3/4,1)}.

For a dyadic interval I=[a,a+r) put \mathfrak h_I(y)=\mathfrak h((y-a)/r) and \mathfrak g_I(x)=\mathfrak g((x-a)/r). The shift summand

\left(\frac1{|I|}\int_I \rho(y)\chi_n(y)\mathfrak h_I(y)\,dy\right) \mathfrak g_I(x)

is an admissible summand in Proposition 6.1: on the output side choose \epsilon_j=r_{j+1}, so that \mathfrak g_I=r_{j+1}r_{j+2}.

Tile the line by top intervals of length s\in[1/2,1], translate the tiling by a\in[0,s), and use depths j=0,\ldots,J-1 below every top, where J\ge1. Denote the resulting finite-depth shift by \mathcal S_{a,s,J}. This shift need not be periodic, so in the next estimate G_\rho denotes its representative in [0,1). At most a fixed number of top intervals meet [0,1). Each has input average at most 2, since \rho has mass one in a period and s\ge1/2. At every x\in G_\rho, the density cutoff in Proposition 6.1 is inactive on all intervals of its spatial path. Thus the stopped and unstopped shifts agree there, and

\int_{G_\rho}\sup_{n\in\mathbb Z} |\mathcal S_{a,s,J}(\rho\chi_n)(x)|\,dx\le Cd \tag{10.3}

uniformly in a,s,J. Here a uniform-probability integral on a top interval is multiplied by its length to obtain the ordinary spatial integral. Only a bounded number of top intervals is involved.

Translation and scale averaging. Define

\kappa(u)=\int_{\mathbb R}\mathfrak g(v)\mathfrak h(v-u)\,dv.

For a fixed length r=s2^{-j}, uniform translation modulo s is also uniform translation modulo r. Changing variables from an interval’s left endpoint to the output’s position v in that interval shows that the translation average of its full tiling is convolution with r^{-1}\kappa(t/r). This identity is just Fubini’s theorem applied to the bounded, compactly supported profiles. For F=\rho\chi_n, integrating the translation average in ds/s gives

\int_{1/2}^1\int_0^s\mathcal S_{a,s,J}F(x) \,\frac{da}{s}\,\frac{ds}{s}=(K_J*F)(x),

where, because the intervals [2^{-j-1},2^{-j}] partition [2^{-J},1],

K_J(t)=\int_{2^{-J}}^1\kappa(t/r)\,\frac{dr}{r^2}. \tag{10.4}

The triangle inequality, followed by Equation (10.3), therefore yields

\int_{G_\rho}\sup_{n\in\mathbb Z} |K_J*(\rho\chi_n)(x)|\,dx\le Cd. \tag{10.5}

The measure ds/s on [1/2,1] has mass \log2; it has not been normalized to a probability measure. This fixed factor is absorbed in C. The functions being averaged are finite shift outputs, measurable in the translation and scale parameters. No measurability of the auxiliary processing choices used to prove their uniform bound is required.

The nonzero singular kernel. Reflection about 1/2 fixes \mathfrak g and reverses the sign of \mathfrak h, so \kappa is odd. It is supported in [-1,1] and is Lipschitz, since it is a finite linear combination of overlaps of intervals. Moreover

\begin{align*} c:=\int_0^1\kappa(u)\,du &=\int_0^1\mathfrak g(v) \left(\int_0^v\mathfrak h(t)\,dt\right)\,dv \\ &=\frac1{32}-\frac3{32}-\frac3{32}+\frac1{32} =-\frac18. \tag{10.6}\end{align*}

The inner primitive is the tent function \min(v,1-v). Thus the average in Equation (10.4) has a nonzero Hilbert-transform component. This calculation agrees with the Haar-shift kernel evaluation in [12, Section 3, pp. 1135–1136].

Put \delta=2^{-J}. When \delta<t\le1, the substitution u=t/r and the support of \kappa give

K_J(t)=\frac1t\int_t^1\kappa(u)\,du =\frac ct-\frac1t\int_0^t\kappa(u)\,du.

The second term is bounded by an absolute constant; oddness gives the same conclusion for negative t. On |t|\le\delta,

|K_J(t)|\le \|\kappa\|_\infty \int_\delta^1\frac{dr}{r^2}\le C\delta^{-1}.

Compare K_J to

\pi c\,\cot(\pi t)\, \mathbf1_{\{\delta<|t|\le1/2\}}.

On the central annulus the difference is bounded, because \pi\cot(\pi t)-1/t is bounded on [-1/2,1/2]. On 1/2<|t|\le1 the original kernel is bounded as well. Those bounded errors have convolution at most C in absolute value against any \rho\chi_n, using periodicity and mass one. The inner window contributes at most

C\delta^{-1}\int_{|t|\le\delta}\rho(x-t)\,dt \le C M\rho(x)\le Ce^{4h}\qquad(x\in G_\rho).

It follows from Equation (10.5) and e^{4h}=o(d) that

\int_{G_\rho}\sup_{n\in\mathbb Z} |H_\delta(\rho\chi_n)(x)|\,dx\le Cd, \qquad H_\delta F(x)= \int_{\delta<|t|\le1/2}F(x-t)\cot(\pi t)\,dt. \tag{10.7}

Division by the fixed nonzero number \pi c is harmless.

The Hilbert multiplier and Fourier projections. For completeness, H_\delta converges in L^2 to the periodic Hilbert transform H whose multiplier is -i\,\operatorname{sgn}(k) for k\ne0 and zero at k=0. To see this directly, oddness reduces its multiplier to a sine integral. These multipliers are uniformly bounded: write \cot(\pi t)=1/(\pi t)+O(1), split the sine integral at inverse frequency, and integrate by parts on the remaining interval. For a fixed positive integer k,

\int_{-1/2}^{1/2}\sin(2\pi kt)\cot(\pi t)\,dt=1. \tag{10.8}

Indeed, setting x=\pi t and using

\sin(2kx)=2\sin x\sum_{a=1}^k\cos((2a-1)x)

reduces the integral to the constant term in 2\cos x\sum_{a=1}^k\cos((2a-1)x); all its nonconstant even-frequency cosines have integral zero on [-\pi/2,\pi/2]. The constant term occurs only for a=1 and gives one after division by \pi. Negative frequencies have the opposite sign. Pointwise multiplier convergence and the uniform bound now imply L^2 convergence by Parseval’s identity.

Every \rho\chi_n is bounded and hence belongs to L^2. For a fixed finite set of integers, the maximum of the differences |H_\delta(\rho\chi_n)-H(\rho\chi_n)| tends to zero in L^1, by the finite sum of their L^2 norms. Pass to the limit in Equation (10.7) for that finite set and then increase the set to \mathbb Z. Monotone convergence gives

\int_{G_\rho}\sup_{n\in\mathbb Z}|H(\rho\chi_n)(x)|\,dx\le Cd. \tag{10.9}

Let P_0 be projection onto constants and let \Pi_{\ge0} be projection onto the nonnegative Fourier frequencies. As L^2 multiplier identities,

\Pi_{\ge0}=\frac{I+iH+P_0}{2}, \qquad S_N\rho= \chi_{-N}\Pi_{\ge0}(\rho\chi_N) -\chi_{N+1}\Pi_{\ge0}(\rho\chi_{-(N+1)}).

They hold almost everywhere simultaneously for all the indicated integer indices, after removing one null set. Since |\widehat\rho(k)|\le1, their identity and constant terms give the pointwise upper bound

\sup_{N\ge0}|S_N\rho| \le \rho+1+\sup_{n\in\mathbb Z}|H(\rho\chi_n)|.

Integrating and applying Equation (10.9) proves Equation (10.2).

10.2 Bounded inputs

Proposition 10.2. Every bounded complex-valued periodic function has almost-everywhere convergence of its symmetric Fourier partial sums along the full sequence.

Proof. Let \|f\|_\infty\le B. The case B=0 is immediate. Take the Fejér means P_\ell of f. They are trigonometric polynomials satisfying \|P_\ell\|_\infty\le B and \|f-P_\ell\|_1\to0. The latter fact follows from the approximate-identity property of the nonnegative, mass-one Fejér kernels and continuity of translations in L^1.

Split f-P_\ell into the positive and negative parts of its real and imaginary components. Each nonnegative part q_\ell is at most 2B and has mass m_\ell\to0. If m_\ell=0 its partial sums vanish. Otherwise normalize to \rho_\ell=q_\ell/m_\ell and choose

d_\ell=d_*+\max\{0,\log(2B/m_\ell)\}.

As m_\ell\to0 through positive values, d_\ell\to\infty, m_\ell d_\ell\to0, and \rho_\ell\le e^{d_\ell}. Theorem 10.1 and Markov’s inequality show, for each \eta>0,

\left|\left\{\sup_N|S_Nq_\ell|>\eta\right\}\right| \le C(\log d_\ell)^{-4} +\frac{Cm_\ell d_\ell}{\eta}\longrightarrow0.

Applying this to the four component parts proves \sup_N|S_N(f-P_\ell)|\to0 in measure.

For each \ell, the partial sums of P_\ell are eventually equal to P_\ell itself. Consequently the measurable function

A(x)=\limsup_{N\to\infty}|S_Nf(x)-f(x)|

satisfies

A(x)\le \sup_N|S_N(f-P_\ell)(x)|+|f(x)-P_\ell(x)|

almost everywhere. Both terms on the right tend to zero in measure. For each \eta>0 this bounds |\{A>\eta\}| by a quantity tending to zero. Thus A=0 almost everywhere, as required.

10.3 Summation of the height layers

Proof of Theorem 1.1. It suffices to work in unit-period coordinates. Write

f=f_{\mathrm{bd}}+\sum_{j\ge1}f_j,\qquad f_{\mathrm{bd}}=f\,\mathbf1_{\{|f|<e^2\}},\qquad f_j=f\,\mathbf1_{\{e^{2^j}\le|f|<e^{2^{j+1}}\}}.

For each j, split f_j into four nonnegative parts with coefficients 1,-1,i,-i, using the positive and negative parts of its real and imaginary components. Consider one such sequence q_j, and write m_j=\int q_j. Discard its zero-mass members. Put

\rho_j=q_j/m_j,\qquad d_j=2^{j+1}+\log^-m_j.

Then \int\rho_j=1 and \rho_j\le e^{d_j}. For all sufficiently large j, Theorem 10.1 applies.

We claim

\sum_j m_jd_j<\infty. \tag{10.10}

On the support of q_j, \log|f|\ge2^j, and q_j\le|f|. Disjointness of the layers therefore gives

\sum_j2^{j+1}m_j \le2\int_{\{|f|\ge e^2\}}|f|\log|f|<\infty.

For the additional entropy term, if m_j\ge e^{-2^j}, then m_j\log^-m_j\le2^jm_j. If m_j<e^{-2^j}, the monotonicity of t\log(1/t) on (0,e^{-1}) gives m_j\log(1/m_j)\le2^je^{-2^j}. Both bounds are summable. This proves Equation (10.10).

Let G_j be the good set associated to \rho_j. Since \log d_j\ge(j+1)\log2,

\sum_j|G_j^c|\le C\sum_j(\log d_j)^{-4}<\infty.

Borel–Cantelli implies that almost every point belongs to all but finitely many G_j. On the other hand, Tonelli’s theorem and Equations (10.2) and (10.10) give

\int\sum_j\mathbf1_{G_j}(x)\sup_N|S_Nq_j(x)|\,dx \le C\sum_jm_jd_j<\infty,

where a finite number of initial layers may be omitted. Therefore, at almost every point, the partial-sum maxima of the tail layers are absolutely summable. Every individual q_j is bounded, so Proposition 10.2 also gives its full-sequence limit and the finiteness of its partial-sum maximum outside a null set. Restoring finitely many omitted layers causes no difficulty. Intersect these full-measure sets for the four component sequences and for f_{\mathrm{bd}}.

The layers are summable in L^1. For each fixed N, Fourier coefficients can consequently be summed term by term; equivalently, the finite-rank operator S_N is bounded from L^1 to L^\infty. Thus

S_Nf=S_Nf_{\mathrm{bd}}+ \sum_{j\ge1}S_Nf_j

with the four component parts understood in the last sum. At every point of the full-measure set just constructed, the absolute sum of their partial-sum maxima is finite. Its tail tends to zero uniformly in N. Take the limit first on a finite collection of layers, using Proposition 10.2, and then let the collection increase. The result is S_Nf(x)\to f_{\mathrm{bd}}(x)+\sum_jf_j(x)=f(x). All excluded sets are measurable, and only countably many are used. Changing back to x=2\pi y proves the stated theorem with normalized measure dx/(2\pi).

10.4 The maximal consequence

Corollary 10.3 (Weak-L^1 maximal bound). On \mathbb T=\mathbb R/(2\pi\mathbb Z) with d\mu=dx/(2\pi), set \Phi(t)=t\log(2+t) and define the Luxemburg norm by

\|f\|_\Phi =\inf\left\{b>0:\int_{\mathbb T}\Phi(|f|/b)\,d\mu\le1\right\}.

For the ordinary symmetric partial sums of Theorem 1.1, put S^*f=\sup_{N\ge0}|S_Nf|. There is an absolute constant C such that

\sup_{\lambda>0}\lambda\,\mu\{S^*f>\lambda\} \le C\|f\|_\Phi \qquad(f\in L\log L(\mathbb T;\mathbb C)). \tag{10.11}

Proof. The real Luxemburg space B=L^\Phi(\mathbb T;\mathbb R) is Banach. Since \Phi(2t)\le4\Phi(t), its underlying set is exactly the real L\log L class used in Theorem 1.1; moreover \|f\|_1\le\|f\|_\Phi/\log2. Translations are isometries of B by invariance of \mu. Each S_N is a real linear operator commuting with translations and is bounded from B into itself: indeed, \|S_Nf\|_\infty\le(2N+1)\|f\|_1, and L^\infty embeds continuously in B. These bounds need not be uniform in N. By Theorem 1.1, S^*f is finite almost everywhere for every f\in B. Stein’s theorem [22, Appendix, Theorem 10], applied to the sequence S_0,-S_0,S_1,-S_1,\ldots, therefore gives Equation (10.11) for real inputs; the signed supremum in that theorem is precisely S^*f for this sequence. For f=u+iv, use S^*f\le S^*u+S^*v and \|u\|_\Phi,\|v\|_\Phi\le\|f\|_\Phi, splitting the level at \lambda/2. This proves the complex-valued assertion after enlarging C by a fixed factor.

References

  1. [1]N. Yu. Antonov. Convergence of Fourier series. East Journal on Approximations, 2(2):187–196, 1996.
  2. [2]Juan Arias-de-Reyna. Pointwise convergence of Fourier series. Journal of the London Mathematical Society, 65(1):139–153, 2002. doi:10.1112/S0024610701002824.
  3. [3]Lennart Carleson. On convergence and growth of partial sums of Fourier series. Acta Mathematica, 116:135–157, 1966. doi:10.1007/BF02392815.
  4. [4]María J. Carro, Mieczysław Mastyło, and Luis Rodríguez-Piazza. Almost everywhere convergent Fourier series. Journal of Fourier Analysis and Applications, 18(2):266–286, 2012. doi:10.1007/s00041-011-9199-9.
  5. [5]Mei-Chu Chang. A polynomial bound in Freiman’s theorem. Duke Mathematical Journal, 113(3):399–419, 2002. doi:10.1215/S0012-7094-02-11331-3.
  6. [6]Thomas M. Cover and Joy A. Thomas. Elements of Information Theory. John Wiley & Sons, Second edition, 2006. doi:10.1002/047174882X.
  7. [7]Francesco Di Plinio and Anastasios Fragkos. The weak-type Carleson theorem via wave packet estimates. Transactions of the American Mathematical Society, 378(3):1551–1592, 2025. doi:10.1090/tran/9303, arXiv:2204.08051.
  8. [8]Charles Fefferman. Pointwise convergence of Fourier series. Annals of Mathematics. Second Series, 98(3):551–571, 1973. doi:10.2307/1970917.
  9. [9]Charles L. Fefferman. Erratum: Pointwise convergence of Fourier series. Annals of Mathematics. Second Series, 146(1):239, 1997. doi:10.2307/2951836.
  10. [10]Jean Ginibre. General formulation of Griffiths’ inequalities. Communications in Mathematical Physics, 16(4):310–328, 1970. doi:10.1007/BF01646537.
  11. [11]Richard A. Hunt. On the convergence of Fourier series. In Deborah Tepper Haimo, editor, Orthogonal Expansions and their Continuous Analogues, pages 235–255, Southern Illinois University Press, 1968. Proceedings of the conference held in Edwardsville, Illinois, April 27–29, 1967.
  12. [12]Tuomas Hytönen. On Petermichl’s dyadic shift and the Hilbert transform. Comptes Rendus Mathématique, 346(21–22):1133–1136, 2008. doi:10.1016/j.crma.2008.09.021.
  13. [13]Andrey Kolmogoroff. Une série de Fourier–Lebesgue divergente presque partout. Fundamenta Mathematicae, 4:324–328, 1923. doi:10.4064/fm-4-1-324-328.
  14. [14]Michael Lacey and Christoph Thiele. A proof of boundedness of the Carleson operator. Mathematical Research Letters, 7(4):361–370, 2000. https://archive.intlpress.com/site/pub/files/_fulltext/journals/mrl/2000/0007/0004/MRL-2000-0007-0004-a001.pdf.
  15. [15]James R. Lee. Covering the large spectrum and generalized Riesz products. SIAM Journal on Discrete Mathematics, 31(1):562–572, 2017. doi:10.1137/15M1048604, arXiv:1508.07109.
  16. [16]Victor Lie. On the boundedness of the Carleson operator near L^1. Revista Matemática Iberoamericana, 29(4):1239–1262, 2013. doi:10.4171/rmi/755.
  17. [17]Victor Lie. Pointwise convergence of Fourier series (I). On a conjecture of Konyagin. Journal of the European Mathematical Society, 19(6):1655–1728, 2017. doi:10.4171/JEMS/703, arXiv:1408.4783.
  18. [18]Stefanie Petermichl. Dyadic shifts and a logarithmic estimate for Hankel operators with matrix symbol. Comptes Rendus de l’Académie des Sciences, Série I, Mathématique, 330(6):455–460, 2000. doi:10.1016/S0764-4442(00)00162-2.
  19. [19]Per Sjölin. An inequality of Paley and convergence a.e. of Walsh–Fourier series. Arkiv för Matematik, 7(6):551–570, 1969. doi:10.1007/BF02590894.
  20. [20]Fernando Soria. Note on differentiation of integrals and the halo conjecture. Studia Mathematica, 81(1):29–36, 1985. doi:10.4064/sm-81-1-29-36.
  21. [21]Fernando Soria. On an extrapolation theorem of Carleson–Sjölin type with applications to a.e. convergence of Fourier series. Studia Mathematica, 94(3):235–244, 1989. doi:10.4064/sm-94-3-235-244.
  22. [22]Elias M. Stein. On limits of sequences of operators. Annals of Mathematics. Second Series, 74(1):140–170, 1961. doi:10.2307/1970308.
  23. [23]Jean Ville. Étude critique de la notion de collectif. PhD thesis, Faculté des sciences de Paris, 1939. https://numdam.org/item/THESE_1939__218__1_0/.