Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Song–Zhang v1: polynomial coefficients control the spectral gap

Part of the first version of Song–Zhang, Chapter Song–Zhang, first version: polynomial estimates and curvature; the reading order is on the full proofs page.

Overview. This reconstructs Section 5 of the version-pinned Song & Zhang, 2026. The argument starts from a first eigenfunction, iterates a normalized inverse square root followed by differentiation, and bounds the mass removed by centering. A combinatorial recovery inequality converts polynomial testing into a bound on this loss. A convolution estimate absorbs all normalization defects, uniformly in the number of generations and the terminal polynomial degree. The only quantitative input beyond the profile in the statement is the universal coefficient estimate Theorem 7.1.

Analytic setup and normalized families

The analytic preparation Lemma 7.1, proved in Lemma 120.3, establishes the following facts for the present regular measure. The operator H=−Δ+∇W⋅∇H=-\Delta+\nabla W\cdot\nabla is the self-adjoint operator associated with the closed gradient form, its kernel is the constants, and its first positive eigenvalue λ=CP(ν)−1\lambda=C_P(\nu)^{-1} is attained. Its form domain is the weighted Sobolev space. On its operator domain,

∥Hg∥22=E∥D2g∥HS2+E⟨D2W∇g,∇g⟩.\|Hg\|_2^2=\mathbb E\|D^2g\|_{\rm HS}^2+ \mathbb E\langle D^2W\nabla g,\nabla g\rangle.

In particular all second weak derivatives exist in L2L^2 and commute. Inverse powers act on centered functions. For centered form-domain hh, H−1/2h∈Dom⁡(H)H^{-1/2}h\in\operatorname{Dom}(H), ∥HH−1/2h∥22=∥H1/2h∥22\|H H^{-1/2}h\|_2^2=\|H^{1/2}h\|_2^2, and ∥∇H−1/2h∥22=∥h∥22\|\nabla H^{-1/2}h\|_2^2=\|h\|_2^2. Polynomial tests belong to the form domain. These facts apply componentwise to finite families, with summed squared norms.

Put P+h=h−EhP_+h=h-\mathbb Eh, D=P+∇H−1/2D=P_+\nabla H^{-1/2}, and Lh=E[Xh]Lh=\mathbb E[Xh]. Both DD and LL are contractions: for LL use E(z⋅X)2≤∣z∣2\mathbb E(z\cdot X)^2\le|z|^2 and duality. Testing the form against xix_i gives E∇H−1/2h=LH1/2h\mathbb E\nabla H^{-1/2}h=LH^{1/2}h.

For a nonzero centered form-domain family uu, define

v=∥u∥22,b=∥H−1/2u∥22,β=v/b,h=βH−1/2u,χ=∥H1/2u∥22−βv.v=\|u\|_2^2,\quad b=\|H^{-1/2}u\|_2^2,\quad \beta=v/b,\quad h=\sqrt\beta H^{-1/2}u,\quad \chi=\|H^{1/2}u\|_2^2-\beta v.

The scalar β\beta belongs to the entire family, not individual components. Cauchy--Schwarz gives v2≤b∥H1/2u∥22v^2\le b\|H^{1/2}u\|_2^2, whence λ≤β≤∥H1/2u∥22/v\lambda\le\beta\le\|H^{1/2}u\|_2^2/v and χ≥0\chi\ge0. Furthermore

∥h∥22=v,H1/2h=βu,z:=H−1/2h−β−1/2u\|h\|_2^2=v,\qquad H^{1/2}h=\sqrt\beta u, \qquad z:=H^{-1/2}h-\beta^{-1/2}u

satisfies the exact identity

∥∇z∥22=βb−2v+β−1∥H1/2u∥22=χ/β.\|\nabla z\|_2^2 =\beta b-2v+\beta^{-1}\|H^{1/2}u\|_2^2=\chi/\beta.

This uses spectral calculus and the form identity; it does not commute a spatial derivative with an inverse spectral power.

Start with a centered unit first eigenfunction F0F_0. Set uj=DFju^j=DF_j and normalize uju^j as above to obtain Fj+1F_{j+1}, normalizer βj+1\beta_{j+1}, and defect χj\chi_j. Write

vj=∥Fj∥22,ej=∥H1/2Fj∥22,pj=∥LH1/2Fj∥2,Ej=ej−λvj.v_j=\|F_j\|_2^2,\quad e_j=\|H^{1/2}F_j\|_2^2,\quad p_j=\|LH^{1/2}F_j\|^2,\quad E_j=e_j-\lambda v_j.

The gradient before centering has norm squared vjv_j and mean norm squared pjp_j. Bochner applied to H−1/2FjH^{-1/2}F_j bounds its derivative energy by ej−avje_j-av_j. Normalization removes exactly χj\chi_j. Thus

vj+1=vj−pj,ej+1≤ej−avj−χj,0≤Ej+1≤Ej+λpj−avj−χj.(1)v_{j+1}=v_j-p_j,\qquad e_{j+1}\le e_j-av_j-\chi_j, \qquad 0\le E_{j+1}\le E_j+\lambda p_j-av_j-\chi_j. \tag{1}

In particular vj≤1v_j\le1, ej≤λe_j\le\lambda, and pj≤min⁡(vj,ej)≤λp_j\le\min(v_j,e_j)\le\lambda. At the first step, even if u0=0u^0=0, λ(1−p0)≤∥H1/2u0∥22≤λ−a\lambda(1-p_0)\le\|H^{1/2}u^0\|_2^2\le\lambda-a, so a≤λp0≤λ2a\le\lambda p_0\le\lambda^2. The same inequalities hold with zero successors after a vanishing family, but normalizers will only be used on the nonvanishing stopped prefix constructed below. Each uju^j is in the form domain by Bochner; hence every operation on that prefix has the required domain.

Recovering a partially symmetric tensor

We will also use, for an ll-slot array YY and its adjacent transpositions sis_i,

∥Y−Sym⁡lY∥≤l∑i=1l−1∥Y−siY∥.(2)\|Y-\operatorname{Sym}_lY\|\le l\sum_{i=1}^{l-1}\|Y-s_iY\|. \tag{2}

Indeed insertion sort writes each permutation with each generator at most ll times. For a product of isometries, telescope I−U1⋯UmI-U_1\cdots U_m as ∑i=1mU1⋯Ui−1(I−Ui)\sum_{i=1}^mU_1\cdots U_{i-1}(I-U_i); all differences are thereby applied to the original array. Averaging the resulting bounds proves (2).

Polynomial tests, adjacent swaps, and the dyadic loss

Let Ak\mathcal A_k be the tensor-valued Appell polynomial, so Pk[T]=⟨Ak,T⟩P_k[T]=\langle\mathcal A_k,T\rangle, and put Qkh=E[Akh]Q_kh=\mathbb E[\mathcal A_kh]. Duality gives ∥Qk∥≤Kk\|Q_k\|\le\sqrt{K_k}. If Pk+1\mathsf P_{k+1} symmetrizes the kk polynomial indices and the newest derivative index, then for j≥1j\ge1

(k+1)Pk+1QkDFj=Qk+1H1/2Fj=βjQk+1uj−1.(3)(k+1)\mathsf P_{k+1}Q_kDF_j =Q_{k+1}H^{1/2}F_j =\sqrt{\beta_j}Q_{k+1}u^{j-1}. \tag{3}

To check this, pair against an arbitrary symmetric (k+1)(k+1)-tensor, use ∇Pk+1=(k+1)Pk\nabla P_{k+1}=(k+1)P_k with one free index, and test the form against H−1/2FjH^{-1/2}F_j. Centering contributes zero because EAk=0\mathbb E\mathcal A_k=0.

Suppose the normalizers in question satisfy βj≤bλ\beta_j\le b\lambda, b≥1b\ge1. The normalization defect gives

uj+1=βj+1−1/2P+∇uj+P+∇zj,∥∇zj∥22=χj/βj+1.u^{j+1}=\beta_{j+1}^{-1/2}P_+\nabla u^j+P_+\nabla z^j, \qquad \|\nabla z^j\|_2^2=\chi_j/\beta_{j+1}.

The first term is symmetric in its newest two indices: uju^j is a centered gradient and taking another weak derivative gives a Hessian. Thus the newest swap has defect at most 2χj/λ2\sqrt{\chi_j/\lambda}. Each subsequent map DβH−1/2D\sqrt{\beta}H^{-1/2} has norm at most b\sqrt b and commutes with permutations of old indices. Its scalar is kept fixed at the value for the original entire family, even when the map is applied to a difference. Numbering derivative slots newest first, the swap of slots a,a+1a,a+1 in uJu^J was created at uJ−a+1u^{J-a+1} and underwent exactly a−1a-1 subsequent maps. Hence

∥uJ−sauJ∥2≤2b(a−1)/2χJ−a/λ.(4)\|u^J-s_au^J\|_2\le2b^{(a-1)/2}\sqrt{\chi_{J-a}/\lambda}. \tag{4}

No old defect has been differentiated.

Fix an integer L≥3L\ge3 and put

B=2L/(L−2),A=bB2,J=B4bL+1,md=(L+1)d/2−1.B=2\sqrt{L/(L-2)},\qquad A=bB^2,\qquad J=B^4b^{L+1},\qquad m_d=(L+1)d/2-1.

For j≥mdj\ge m_d and dyadic d≥2d\ge2, we claim

pj≤(Aλ)d/2cd+∑k<dk dyadictk∑a=1Lk−1χj−k−a,tk=4LkBJk/2ckλ(k−1)/2.(5)\sqrt{p_j}\le(A\lambda)^{d/2}c_d+ \sum_{\substack{k<d\\k\text{ dyadic}}}t_k \sum_{a=1}^{Lk-1}\sqrt{\chi_{j-k-a}}, \quad t_k=\frac{4Lk}{B}J^{k/2}c_k\lambda^{(k-1)/2}. \tag{5}

Here and below dyadic indices start at one. To prove (5), set Tk=Qkuj−kT_k=Q_ku^{j-k} at successive dyadic degrees. There are kk symmetric polynomial slots and j−k+1≥Lkj-k+1\ge Lk derivative slots for k<dk<d. Symmetrize the newest LkLk derivative slots to obtain SkS_k. The block lemma with (s,l,q)=(k,Lk,2k)(s,l,q)=(k,Lk,2k) has recovery constant

Ck2=(2kk)(Lk)k‾((L−1)k)k‾≤4k(L/(L−2))k=B2k.C_k^2=\binom{2k}k\frac{(Lk)_{\underline k}}{((L-1)k)_{\underline k}} \le4^k(L/(L-2))^k=B^{2k}.

Every denominator factor exceeds (L−2)k(L-2)k, which proves this bound. Repeated application of (3), whose inner projections are absorbed by the outer symmetrization, yields

∥P2kTk∥≤(bλ)k/2k!(2k)!∥T2k∥.\|\mathsf P_{2k}T_k\|\le(b\lambda)^{k/2}\frac{k!}{(2k)!}\|T_{2k}\|.

Equations (2), (4) and ∥Qk∥≤Kk\|Q_k\|\le\sqrt{K_k} give

Ek:=∥Tk−Sk∥≤2LkbLk/2Kk/λ∑a=1Lk−1χj−k−a.E_k:=\|T_k-S_k\| \le2Lk b^{Lk/2}\sqrt{K_k/\lambda} \sum_{a=1}^{Lk-1}\sqrt{\chi_{j-k-a}}.

Thus ∥Tk∥≤Bk(bλ)k/2k!/(2k)! ∥T2k∥+2BkEk\|T_k\|\le B^k(b\lambda)^{k/2}k!/(2k)!\,\|T_{2k}\|+2B^kE_k. The initial factor is pj=βj∥T1∥\sqrt{p_j}=\sqrt{\beta_j}\|T_1\|. The preceding dyadic degrees sum to k−1k-1, and the factorials telescope, so the coefficient before stage kk is at most

(bλ)k/2Bk−1/k!.(b\lambda)^{k/2}B^{k-1}/k!.

At the terminal stage ∥Td∥≤Kd\|T_d\|\le\sqrt{K_d} because ∥uj−d∥2≤1\|u^{j-d}\|_2\le1. Its contribution is (Aλ)d/2cd/B(A\lambda)^{d/2}c_d/B. Multiplying the stage coefficient by 2BkEk2B^kE_k gives exactly tkt_k above. Every defect index is nonnegative: j≥(L+1)k−1j\ge(L+1)k-1. This proves (5), including the boundary degree d=2d=2.

Uniform absorption and stopping

Fences respected. The proposed node has no bounded_by edges. The argument proves a curvature-dependent general-test estimate, not CMH, a sharp third-moment constant, or a universal-time occupation estimate. In particular it does not invert any sufficient-condition implication in the brief. Its constants require R≥240ε−2R\ge2^{40}\varepsilon^{-2} uniformly in the degree; retaining this threshold is essential when the theorem is iterated with depth-dependent ε\varepsilon. All polynomial assumptions are hypotheses of the displayed implication; the universal small-degree estimate is an actual dependency, not an assumed open antecedent. Independent review of this proof and that dependency is required.

References
  1. Song, Z., & Zhang, X. (2026). An O(4\log^* n) Bound for the KLS Constant. https://arxiv.org/abs/2610.01447v1