Alphabeta Math
Pipeline-generated
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

Overparameterized networks and training to interpolation

1 · Prerequisites

2 · Logical flowchart

3 · Definitions, theorems and proofs

LemmaStatement: Literature-sourcedProof: AI-generatedjudge pass (gpt-6-luna)Open item page →

The Gaussian activation normalizer is finite and positive

Statement

Assume AC. Let Z have the standard normal law on R. Suppose σ:R→R is analytic, nonpolynomial, satisfies ∣σ(0)∣≤B for some finite B, and is L-Lipschitz for a finite L. Then 0<E[σ(Z)2]<∞. Consequently the normalization constant

cσ:=(E[σ(Z)2])−1

is a well-defined positive finite real number.

Facts & Assumptions

Given: AC, the activation σ, its constants B,L, and a standard normal random variable Z.

[A1]

AC is assumed by the standard-normal-law and Gaussian-moment suppliers. The Axiom of Choice

[F1]

Proof

technique · bound the activation by an affine function of $|Z|$, then use continuity and nonpolynomiality for strict positivity
1.1A1F1givenalgebra

Lipschitz continuity gives ∣σ(z)∣≤B+L∣z∣, hence σ(z)2≤2B2+2L2z2. By [F1], Eσ(Z)2≤2B2+2L2<∞.

1.2A1F1givenchoose

A nonpolynomial analytic function is not identically zero, so choose z0 with σ(z0)≠0. Continuity gives an open interval around z0 on which ∣σ∣≥∣σ(z0)∣/2. The standard normal density is strictly positive on that interval, so its integral of σ2 is positive.

2.1step 1.1step 1.2algebra∎

The expectation is both finite and strictly positive by steps 1.1 and 1.2, so its reciprocal exists and is positive and finite.

LemmaStatement: Literature-sourcedProof: AI-generatedjudge pass (gpt-6-luna)Open item page →

Gram control and a Taylor remainder imply one-step loss contraction

Statement

Let r∈Rn, let G be a real symmetric matrix with eigenvalues in [λ,Λ] for some 0<λ≤Λ, and suppose 0<ηΛ≤1. If an update has the form r+=(I−ηG)r+R with ∥R∥2≤ηλ∥r∥2/4, then for the quadratic loss L=12∥r∥22 and L+=12∥r+∥22,

L+≤(1−ηλ/2)L.

Facts & Assumptions

Proof

technique · bound the spectral norm of the linear contraction and then square the residual-norm bound
1.1F1F2givenalgebra

By [F1], the eigenvalues μj of G lie in [λ,Λ], so those of I−ηG lie in [1−ηΛ,1−ηλ]⊆[0,1−ηλ]. Expanding r in an orthonormal eigenbasis gives ∥(I−ηG)r∥22=∑j(1−ημj)2rj2≤(1−ηλ)2∥r∥22. The norm triangle inequality [F2] with the remainder hypothesis then yields ∥r+∥2≤(1−3ηλ/4)∥r∥2.

2.1step 1.1algebra∎

Put t=ηλ∈(0,1]. Then (1−3t/4)2≤1−t/2, since 1−t/2−(1−3t/4)2=t(1−9t/16)≥7t/16>0. Squaring the estimate in step 1.1 and multiplying by 1/2 yields L+≤(1−ηλ/2)L.

DefinitionDefinition: Literature-sourcedProof: Not applicablejudge pass (gpt-6-luna)Open item page →

Fully connected and residual networks with the Du Gram kernel

Definition

Assume AC, DC and countable choice for the Gaussian probability-law constructions below. The displayed finite real network recurrences and empirical loss are deterministic definitions for fixed parameters. Fix positive integers d,n,H,m, training inputs x1,…,xn∈Rd with ∥xi∥2=1, labels y1,…,yn∈R, an analytic nonpolynomial activation σ satisfying ∣σ(0)∣<∞ and finite Lipschitz constants for σ and σ′, and the positive finite normalization cσ=(EZ∼N(0,1)[σ(Z)2])−1. A fully connected network of width m and depth H is the finite feedforward recurrence

x(0)=x,x(h)=cσ/m σ(W(h)x(h−1))(1≤h≤H),ffc(x,θ)=aTx(H),

where W(1)∈Rm×d, W(h)∈Rm×m for h≥2, a∈Rm, and σ acts coordinatewise. The residual network uses the same first layer and then

x(1)=cσ/m σ(W(1)x),x(h)=x(h−1)+cresHm σ(W(h)x(h−1))(2≤h≤H),fres(x,θ)=aTx(H),

with 0<cres<1. In either architecture, random initialization means all entries of all matrices W(h) and of a are independent N(0,1) random variables; all these entries are trained by full-batch gradient descent. The empirical quadratic loss is L(θ)=12∑i=1n(f(xi,θ)−yi)2.

For the fully connected architecture and the same fixed data, define the population matrices by

Kij(0)=⟨xi,xj⟩.

For 1≤h<H, let (Uij(h),Vij(h)) be a centered bivariate normal vector with covariance

Aij(h)=(Kii(h−1)Kij(h−1)Kji(h−1)Kjj(h−1)),

and put Kij(h)=cσE[σ(Uij(h))σ(Vij(h))]. For h=H, use the same bivariate-normal definition with Aij(H) formed from K(H−1), and set

Kij(H)=cσKij(H−1)E[σ′(Uij(H))σ′(Vij(H))].

When H=1, this last formula uses K(0) directly and there are no intermediate activation-covariance matrices. For the fully connected theorem, this terminal matrix is denoted Kfc(H).

For the residual architecture with H≥2, use a separate sequence S(h) and mean vector b(h). Set S(0)=K(0). For G(0)∼Nn(0,S(0)), put

Sij(1)=cσE[σ(Gi(0))σ(Gj(0))],bi(1)=cσE[σ(Gi(0))].

Write α=cres/H. For 2≤h<H, take a fresh centered Gaussian vector G(h−1) with covariance S(h−1), put μi(h−1)=E[σ(Gi(h−1))], and define

Sij(h)=Sij(h−1)+α(bi(h−1)μj(h−1)+bj(h−1)μi(h−1))+α2E[σ(Gi(h−1))σ(Gj(h−1))],bi(h)=bi(h−1)+αμi(h−1).

For G∼Nn(0,S(H−1)), define the residual terminal matrix by

Kres,ij(H)=cres2H2Sij(H−1)E[σ′(Gi)σ′(Gj)].

The fully connected and residual convergence statements use Kfc(H) and Kres(H), respectively. At H=1, the two architectures coincide; define Kres(1)=Kfc(1). The residual recursion above is the scalar specialization of the source's general mean/second-moment recursion, equation (16); in particular it includes the mean cross terms. Its terminal Gaussian covariance is S(H−1), as determined by the finite-width last-layer Gram formula (14).

The residual recurrence has one identity contribution from x(h−1) and one nonlinear contribution at every layer 2≤h≤H; the arrows below depict that convention for one generic layer and do not enter any proof.

diagram failed to render:
\begin{tikzpicture}[>=stealth,node distance=1.3cm]
  \node (in) {$x^{(h-1)}$};
  \node[draw,rectangle,right=1.7cm of in] (act) {$\frac{c_{\rm res}}{H\sqrt m}\sigma(W^{(h)}x^{(h-1)})$};
  \node[draw,circle,right=2.0cm of act] (sum) {$+$};
  \node[right=1.2cm of sum] (out) {$x^{(h)}$};
  \draw[->] (in) -- (act);
  \draw[->] (act) -- (sum);
  \draw[->] (in) to[bend right=28] (sum);
  \draw[->] (sum) -- (out);
\end{tikzpicture}
LemmaStatement: Literature-sourcedProof: AI-generatedjudge pass (gpt-6-luna)Open item page →

The recursively defined population Gram matrices are finite, symmetric, and positive semidefinite

Statement

Assume AC, DC and countable choice. Under the hypotheses and notation defining the Du population Gram recursion, every fully connected matrix K(h), every residual feature second-moment matrix S(h), and both terminal matrices are finite, symmetric, and positive semidefinite. The fully connected matrices satisfy Kii(h)=1 for 0≤h<H. For the residual sequence, S(h)−b(h)(b(h))T is positive semidefinite for 1≤h<H; its diagonal is not asserted to equal one.

Facts & Assumptions

Given: The finite unit-norm data, activation, normalization, Gaussian covariance convention, and matrices from the definition.

[F1]

The architecture and matrices are defined by the stated recurrences, with cσ>0 and finite. Fully connected and residual networks with the Du Gram kernel

[F5]

The normalizer satisfies cσE[σ(Z)2]=1 for standard normal Z. The Gaussian activation normalizer is finite and positive

[F6]

A property holding at zero and preserved by successor holds for every natural-number index. The principle of mathematical induction

[F2]

A finite-dimensional positive semidefinite covariance matrix defines a possibly singular centered multivariate normal law, with the listed covariance and finite second moments. Multivariate normal law, including singular covariance

[F4]

Countably many independent copies of a normal vector can be constructed, and expectations of products of integrable independent variables factor. Countably many independent copies of a prescribed law exist Expectations factor over finite products of independent random variables

[A1]

Assume AC as required by the multivariate normal-law supplier. The Axiom of Choice

[F7]

For two Borel probability laws on finite-dimensional real spaces, their product measure exists, has total mass one and satisfies the rectangle formula. Thus its two coordinate projections have the prescribed laws and are independent. For sigma-finite factors, the product measure exists, has the rectangle formula, is sigma-finite, and is unique

[F8]

If X,Y have finite second moments, then E∣XY∣≤(EX2)1/2(EY2)1/2. This ensures integrability of the Gaussian activation products in steps 2.1 and 4.1 and the terminal-coordinate products in steps 3.1 and 5.1. Cauchy-Schwarz for random variables

Proof

technique · induction on the layer index, using Gaussian feature vectors to express each quadratic form as an expectation of a square
1.1baseF1givenalgebra

The input matrix K(0) is symmetric, and for v∈Rn, vTK(0)v=∥∑ivixi∥22≥0. Its diagonal entries are ∥xi∥22=1, giving the base case.

2.1ihA1F1F2F3F5F8step 1.1algebra

Assume the induction hypothesis that K(h−1) is finite, symmetric positive semidefinite, and has diagonal 1, for some 1≤h<H. By [F2]–[F3], take G∼Nn(0,K(h−1)); each pair (Gi,Gj) has covariance Aij(h). The growth bound ∣σ(z)∣≤∣σ(0)∣+L∣z∣ and finite Gaussian second moments put each σ(Gi) in L2; [F8] then makes every product σ(Gi)σ(Gj) integrable. Thus K(h) is finite and symmetric, and for every v, vTK(h)v=cσE[(∑iviσ(Gi))2]≥0. Because Gi∼N(0,1), its diagonal is cσE[σ(Z)2]=1 by [F5]. This establishes the induction step.

3.1A1A2F2F3F4F6step 1.1step 2.1algebradischarge-induction

For the fixed H, let P(h) mean that K(h) is finite, symmetric positive semidefinite, and has diagonal 1, for 0≤h<H. Define Q(k) to mean P(j) for every j≤min⁡(k,H−1). Step 1.1 gives Q(0). If Q(k) and k+1<H, step 2.1 applied with h=k+1 gives P(k+1) and hence Q(k+1); if k+1≥H, Q(k+1) adds no new index. By [F6], Q(H−1) holds, so K(H−1) is finite, symmetric positive semidefinite, and has diagonal 1; when H=1, this is exactly the base matrix K(0). By [F4], take independent G,G′ both distributed as Nn(0,K(H−1)). The vector with coordinates Yi=cσ Gi′σ′(Gi) has finite second moments; σ′ is bounded because it is the derivative of a Lipschitz differentiable function. Factoring expectations of independent variables gives E[YiYj]=cσKij(H−1)E[σ′(Gi)σ′(Gj)]=Kij(H). Therefore K(H) is finite and symmetric, and vTK(H)v=E[(∑iviYi)2]≥0 for every v. This proves the fully connected terminal claim and discharges that induction.

3.2A1F1F2step 2.1constructalgebra

For the residual base layer let G0∼Nn(0,S(0)) and V(1)=cσ σ(G0), with activation applied coordinatewise. Then EV(1)=b(1) and E[V(1)(V(1))T]=S(1) by [F1]. The growth bound used in step 2.1 proves finite second moments. For every v, vTS(1)v=E[(vTV(1))2]≥0, while S(1)−b(1)(b(1))T is its covariance and has quadratic form E[(vT(V(1)−EV(1)))2]≥0. Thus both matrices are symmetric positive semidefinite.

4.1A1A2F1F2F4F6F7F8step 3.2constructalgebra

Suppose a finite-second-moment vector V(h−1) has mean b(h−1) and second-moment matrix S(h−1). Let μ be the Borel law of this generally non-Gaussian vector and let ν=Nn(0,S(h−1)), which exists by [F2]. Both are probability measures and hence sigma-finite. By [F7], on (Rn×Rn,μ×ν) the coordinate projections V and G have laws μ and ν and are independent: for Borel sets A,B, the probability of {V∈A,G∈B} is μ(A)ν(B). Put V(h)=V+ασ(G). Finite second moments follow from the growth bound and (u+v)2≤2u2+2v2. Independence gives EV(h)=b(h−1)+αμ(h−1) and E[Vi(h)Vj(h)]=Sij(h−1)+α(bi(h−1)μj(h−1)+bj(h−1)μi(h−1))+α2E[σ(Gi)σ(Gj)]=Sij(h). For every v, vTS(h)v=E[(vTV(h))2]≥0, and S(h)−b(h)(b(h))T is the covariance of V(h), so its quadratic form is also an expectation of a square. Thus both matrices are finite, symmetric positive semidefinite at this stage, which supplies the covariance hypothesis for the next stage. Iterating through the finitely many layers is justified by [F6].

5.1A1A2F1F2F4step 3.1step 4.1constructalgebradischarge-induction∎

For H≥2, take independent G,G′∼Nn(0,S(H−1)) and set Yi=(cres/H)Gi′σ′(Gi). The boundedness ∣σ′∣≤L gives finite second moments, and independence gives E[YiYj]=Kres,ij(H). Consequently this terminal matrix is symmetric positive semidefinite. At H=1 it equals the already proved fully connected matrix by definition. This proves all residual assertions.

LemmaStatement: Literature-sourcedProof: AI-generatedjudge pass (gpt-6-luna)Open item page →

Gradient descent induces an exact output Gram update with a Taylor remainder

Statement

For any fixed finite real network obeying either architecture recurrence defined on the A page, let u(θ)=(f(xi,θ))i=1n∈Rn, let y=(yi)i=1n, let J(θ)∈Rn×p have row i equal to the gradient of f(xi,θ) with respect to the vector of all trainable parameters, and set G(θ)=J(θ)J(θ)T. For the empirical quadratic loss L(θ)=12∥u(θ)−y∥22, one full-batch gradient-descent step with step size η>0 satisfies

u(θ+)=u(θ)−ηG(θ)(u(θ)−y)+R,

where θ+=θ−ηJ(θ)T(u(θ)−y) and R is the first-order Taylor remainder of the differentiable map u at θ in direction −ηJ(θ)T(u(θ)−y). No size bound on R is asserted here.

Facts & Assumptions

Given: One finite network, its trainable parameter vector, a fixed data set and labels, and a step size η>0.

[F1]

The network is a finite composition of linear maps and an analytic coordinatewise activation, so its output map is continuously differentiable in all trainable parameters. Fully connected and residual networks with the Du Gram kernel

[F2]

Reverse accumulation through a finite C1 parameterized recurrence gives the exact parameter gradient, including all parameter uses. Reverse accumulation sums gradients over shared parameter uses

[F3]

Fixed-step gradient descent updates the parameter by subtracting η times the Euclidean gradient. Gradient descent

[F4]

Differentiability in finite-dimensional real parameter spaces means a linear first-order approximation with a remainder o(∥h∥2). The real Frechet derivative on real and complex matrix spaces with the Frobenius norm

Proof

technique · express the loss gradient by the output Jacobian and apply the defining first-order expansion of the output map
1.1F2F3givenalgebra

Write r=u(θ)−y. By the scalar chain rule applied to 12∑i(ui−yi)2, the parameter gradient is J(θ)Tr. By [F3], h:=θ+−θ=−ηJ(θ)Tr.

2.1F1F4step 1.1algebra∎

Since u is differentiable, [F4] gives u(θ+h)=u(θ)+J(θ)h+R, where R:=u(θ+h)−u(θ)−J(θ)h is the first-order remainder. Substitute h=−ηJ(θ)Tr to get J(θ)h=−ηJ(θ)J(θ)Tr=−ηG(θ)r, which is the claimed identity.

LemmaStatement: Literature-sourcedProof: AI-generatedjudge pass (gpt-6-luna)Open item page →

Gaussian output initialization concentrates conditionally on the hidden features

Statement

Assume AC, DC and countable choice. Fix hidden-layer feature vectors ϕ1,…,ϕn∈Rm, form the matrix Φ∈Rn×m with rows ϕiT, and let the independent output-weight vector a∼Nm(0,Im). Put u=Φa and G=ΦΦT. Then for every 0<δ<1,

P{∥u∥2≤tr⁡(G)+∥Φ∥op2log⁡(1/δ)}≥1−δ,

where ∥Φ∥op=sup⁡∥v∥2=1∥Φv∥2. The probability is conditional on the fixed hidden-layer features and is over a only.

Facts & Assumptions

Given: AC, DC, countable choice, fixed feature vectors, and an independent standard Gaussian output vector.

[F1]

For a standard Gaussian vector Z and an L-Lipschitz real function g, P{g(Z)≥Eg(Z)+r}≤e−r2/(2L2) for L>0,r>0; for L=0 the upper-tail event is empty. Lipschitz concentration for standard Gaussian measure

[F2]

A standard multivariate normal vector has covariance identity, and a fixed linear functional has variance given by its covariance quadratic form. Multivariate normal law, including singular covariance

[F3]

Cauchy–Schwarz bounds E∣XY∣ by the product of the L2 norms. Cauchy-Schwarz for random variables

Proof

technique · apply one-sided Gaussian concentration to the norm of a fixed linear image and bound its mean by its second moment
1.1F2F3F4givenalgebra

[F4] gives ∣∥Φa∥2−∥Φb∥2∣≤∥Φ(a−b)∥2≤∥Φ∥op∥a−b∥2, so a↦∥Φa∥2 is ∥Φ∥op-Lipschitz. Since E[aaT]=Im by [F2], E∥Φa∥22=tr⁡(ΦΦT)=tr⁡(G); [F3] with X=∥Φa∥2 and Y=1 yields E∥Φa∥2≤tr⁡(G).

2.1A1F1F2step 1.1algebra∎

If ∥Φ∥op=0, then u=0 and the asserted event has probability one. Otherwise apply [F1] to g(a)=∥Φa∥2 with r=∥Φ∥op2log⁡(1/δ)>0. The upper-tail probability is at most exp⁡(−log⁡(1/δ))=δ, and the mean bound in step 1.1 gives the displayed threshold.

LemmaStatement: AI-adaptedProof: AI-generatedjudge pass (gpt-6-luna)Open item page →

Local Gram positivity and loss smoothness give confined gradient descent

Statement

Let u:Rp→Rn be C2, L=12∥u−y∥2, and r0=u(θ0)−y. On the closed ball B(θ0,R) suppose a block J∗ of the output Jacobian satisfies J∗J∗T⪰(λ/2)In and ∥∇2L∥≤S, with λ,S>0. If R>22∥r0∥/λ and 0<η≤1/S, every fixed-step gradient-descent iterate stays in the ball,

L(θk)≤(1−ηλ/2)kL(θ0),∑k≥0∥θk+1−θk∥≤22∥r0∥/λ.

The assertion is trivial for r0=0; otherwise the hypotheses imply 0<ηλ/2≤1.

Facts & Assumptions

Given: The output map, quadratic loss, ball, spectral and Hessian bounds, and step size.

[F1]

Gradient descent uses θ+=θ−η∇L(θ). Gradient descent

[F2]

For a scalar C2 function on a segment, the first-order Taylor remainder is the integral of its second derivative against the linear remainder kernel. Taylor's theorem with an integral remainder

[F3]

The Hessian norm is the induced operator norm on the finite-dimensional Euclidean parameter space. The operator norm is zero on the zero domain and otherwise is max_{||v||=1} ||Tv||

[F4]

Euclidean Cauchy–Schwarz gives ∣aTb∣≤∥a∥2∥b∥2 for all real vectors of the same finite dimension. Cauchy-Schwarz ∣⟨x,y⟩∣≤∥x∥2∥y∥2 with its equality case, the triangle inequality for ∥⋅∥2, the parallelogram law and polarisation

Proof

1.1F1F2F3F4givenalgebra

Write g=∇L=JTr. Since JJT⪰J∗J∗T, ∥g∥2=rTJJTr≥(λ/2)∥r∥2=λL in the ball. For any step θ+=θ−hg, 0<h≤η, whose segment stays in the ball, put ψ(s)=L(θ−sg) for 0≤s≤h. Then ψ′(0)=−∥g∥2 and ψ′′(s)=gTHsg, where Hs=∇2L(θ−sg). By Euclidean Cauchy–Schwarz [F4] and the induced norm [F3], ∣gTHsg∣≤∥g∥∥Hsg∥≤∥g∥∥Hs∥∥g∥≤S∥g∥2; for g=0 the same bound is immediate. By [F2], L+=L−h∥g∥2+∫0h(h−s)ψ′′(s) ds≤L−h(1−Sh/2)∥g∥2≤(1−hλ/2)L. If L>0, nonnegativity of L+ gives hλ/2≤1.

2.1step 1.1algebra

For such a segment with L>0, L−L+=(L−L+)/(L+L+)≥h∥g∥2/(4L)≥(λ/4)h∥g∥. If L=0, r=0 and hence g=0, so the step is stationary. Telescoping any finite sequence of admissible full or truncated steps gives total parameter path length at most 4L0/λ=22∥r0∥/λ.

3.1step 1.1step 2.1algebra∎

Suppose there is a first proposed full step exiting the ball. Truncate that step at the first intersection with its boundary; this is a segment of length h∥g∥ with h≤η, wholly in the ball. Step 2.1 applies to all previous steps and this truncated step, giving total path length strictly less than R. Reaching the boundary from θ0 requires length at least R, a contradiction. Thus no exit occurs. Step 1.1 applied at every full step gives the loss rate, and taking the supremum of finite partial path lengths in step 2.1 gives the asserted series bound.

LemmaStatement: AI-adaptedProof: AI-generatedjudge pass (gpt-6-luna)Open item page →

Parameter-ball control of the last-layer Gram and loss Hessian

Statement

Fix finite real parameters obeying either Du architecture recurrence, with all parameters equipped with the Euclidean norm (Frobenius norm on each matrix block). At a fixed initialization assume all hidden feature norms are at most Bx, all matrices W(h)(0) with h≥2 have operator norm at most Bwm, ∥a(0)∥2≤Bam, ∥a(0)∥4≤B4m1/4, and all backward feature derivatives bi,h(0)=∇xi(h)fi have coordinate maxima at most BbM, where M≥1. For fully connected networks allow Bx,Bb≤A0H for a fixed A0≥2, while Bw,Ba,B4 are fixed independent of H. There are A≥2,C≥1 depending only on A0, these fixed bounds and the activation, with the following bounds; write T=AH. In the residual case the constants depend also on fixed cres, and every T below is replaced by a fixed C independent of H. Uniformly on any parameter ball ∥θ−θ0∥≤D with D/m≤c/T (fully connected) or D/m≤c (residual),

max⁡i,h∥xi(h)(θ)−xi(h)(θ0)∥≤TD/m,∥J(θ)∥op≤Tn,

max⁡i∥∇2fi(θ)∥op≤TM(1+D)m,∥JH(θ)JH(θ)T−JH(θ0)JH(θ0)T∥op≤TnDm.

In the residual last inequality the upper bound is CnD/(H2m) for H≥2. The residual first-layer and output-vector derivatives are bounded without a factor 1/H. Consequently, if r0=u(θ0)−y and L=12∥u−y∥2, the loss Hessian on the ball satisfies

∥∇2L∥op≤Tn+Tn (∥r0∥+TnD)M(1+D)m.

The constants in this last bound can again be fixed independent of H for residual networks. For any segment from θ to θ+h contained in the ball, its output Taylor remainder satisfies

∥u(θ+h)−u(θ)−J(θ)h∥≤TnM(1+D)2m∥h∥2,

with T replaced by fixed C in the residual case.

Facts & Assumptions

Given: The deterministic initialization bounds, the architecture, the parameter ball, and fixed activation constants ∣σ(0)∣≤B, ∣σ′∣≤L, ∣σ′′∣≤β.

[F1]

The architectures and scalings are defined locally. Fully connected and residual networks with the Du Gram kernel

[F2]

The chain rule gives the forward variations and backward feature derivatives of the finite network. Reverse accumulation sums gradients over shared parameter uses

[F3]
[F4]

Restricting a scalar C2 output to a parameter line segment gives the first-order Taylor formula with integral remainder used in step 5.1. Taylor's theorem with an integral remainder

[F5]

A C1 vector-valued map whose derivative norm is bounded by M on a segment is M-Lipschitz along that segment. If f:[a,b]→Rm is differentiable with integrable f′ then ∫abf′=f(b)−f(a); and a bounded derivative makes f Lipschitz

[F6]

Matrix operator norms are the induced norms on finite-dimensional inner-product spaces, and the parameter norm uses the Frobenius norm on matrix blocks. The operator norm is zero on the zero domain and otherwise is max_{||v||=1} ||Tv|| The Frobenius norm ∥A∥F=(∑i,j∣aij∣2)1/2 on real or complex matrices

[F7]

For p=4, ∥a∥4=(∑r∣ar∣4)1/4. The p-norms ∥x∥p for rational p≥1, and ∥x∥∞

Proof

1.1F1F2givenalgebraF6

Write zh=Whxh−1 and s1=cσ in both architectures, sh=cσ for fully connected h≥2, and sh=cres/H for residual h≥2. A layer has feature Jacobian Ph=(sh/m)diag⁡(σ′(zh))Wh in the fully connected case, and Ph=I+(sh/m)diag⁡(σ′(zh))Wh in residual layers h≥2. Within the stated ball, ∥Wh∥≤(Bw+D/m)m. Products of consecutive Ph, h≥2, therefore have norm at most AH in the fully connected case and ∏h(1+C/H)≤eC in the residual case. All estimates below use only h≥2 for matrix operator norms; the input is fixed.

2.1F1step 1.1algebra

Let eh=∥xh(θ)−xh(θ0)∥. Lipschitz continuity gives e1≤cσLD/m. For fully connected layers, eh≤Ceh−1+CA0HD/m; for residual layers it gives eh≤(1+C/H)eh−1+CD/(Hm), using the fixed upper feature bound in the residual case and the initial exponential envelope in the fully connected case. Summing the geometric recurrences gives eh≤TD/m in the first case and eh≤CD/m in the second. The imposed smallness of D/m keeps all feature norms in the required bounds and justifies the recurrences successively. Also ∥zh(θ)−zh(θ0)∥≤TD or CD, since it is bounded by ∥Wh(θ)∥eh−1+D∥xh−1(θ0)∥ for h≥2, and by D for h=1.

3.1F1F2step 1.1step 2.1algebraF3

For a unit parameter direction v, let ξh=Dxh[v]. At the first layer ∥ξ1∥≤C∥vW1∥/m. Subsequent layers obey ξh=Phξh−1+(sh/m)diag⁡(σ′(zh))vWhxh−1. The product bounds in step 1.1, the feature bounds in step 2.1, and Cauchy–Schwarz over layer blocks give max⁡h∥ξh∥≤T/m or C/m. For the residual case use ∑h≥2sh∥vWh∥≤cres(∑h≥2∥vWh∥2)1/2/H≤C, retaining the separate first-layer term. Thus ∥Dzh[v]∥≤T or C, uniformly in h. Since Dfi[v]=vaTxH+aTξH and ∥a∥≤Cm, each output gradient has norm at most T or C, proving ∥J∥≤Tn or Cn. This explicitly retains the output and first-layer terms with no 1/H bound.

3.2F2step 1.1step 2.1algebra

Put dh=bh(θ)−bh(θ0). The backward recurrence is bh−1=PhTbh, with dH=a−a0. Decompose its difference as Ph(θ)Tdh+(Ph(θ)−Ph(θ0))Tbh(θ0). The latter is the sum of a weight-difference term bounded by TshD (respectively CshD) and a gate-difference term bounded by CshBbM∥zh(θ)−zh(θ0)∥, because ∥diag⁡(σ′(zh)−σ′(zh,0))bh,0∥≤βBbM∥zh−zh,0∥. Summing the backward recurrence with the product bounds gives ∥dh∥≤TMD or CMD. Hence ∥bh(θ)∥∞≤TM(1+D) or CM(1+D). The crude ℓ2 bound ∥bh(θ)∥2≤Tm or Cm follows separately from step 1.1 and ∥a∥2≤Cm.

3.3F1F2step 2.1algebraF3F7

For the last-hidden-layer Gram, put Xi=xi(H−1), Zi=WHXi and γ=cσ or cres2/H2. Direct differentiation gives GH,ij=γ⟨Xi,Xj⟩m−1∑rar2σ′(Zi,r)σ′(Zj,r). Subtract its initialized version in three parts: feature-inner-product change, derivative-gate change with a02, and output-weight-square change. Step 2.1 bounds the first part by TD/m or CD/m before multiplying by γ. For the second part, Cauchy–Schwarz gives m−1∑ra0,r2∣ΔZi,r∣≤m−1∥a0∥42∥ΔZi∥2≤C∥ΔZi∥2/m; multiply by βL, and use step 2.1. For the third part, m−1∑r∣ar2−a0,r2∣≤m−1∥a−a0∥2(∥a∥2+∥a0∥2)≤CD/m. Feature norms and ∣σ′∣≤L bound all remaining factors. The entry error is consequently at most TD/m or CD/(H2m). The operator norm is at most n times this entry bound.

4.1F1F2step 3.1step 3.2algebra

For two unit parameter directions v,w, twice differentiating the output gives the mixed output-weight terms vaTDxH[w]+waTDxH[v], and the sum over layers of shmbhT(diag⁡(σ′′(zh))(Dzh[v]⊙Dzh[w])+diag⁡(σ′(zh))(vWhDxh−1[w]+wWhDxh−1[v])). This identity follows by differentiating every primitive twice and propagating each resulting second-order term through its downstream Jacobians; the linear skip has zero second derivative. The coordinate-product inequality ∥p⊙q∥1≤∥p∥2∥q∥2, steps 3.1–3.2, and ∑hsh≤CH or C bound the first part by TM(1+D)/m or CM(1+D)/m. The terms involving vWhDxh−1[w] use the ℓ2 backward bound and ∥Dxh−1[w]∥≤T/m or C/m; summing with Cauchy–Schwarz over matrix blocks gives the same bound. The output-weight mixed terms are at most T/m or C/m. This proves the output-Hessian bound, with polynomial factors in H absorbed into AH only for the fully connected case.

5.1step 3.1step 4.1algebraF3F4F5∎

The identity ∇2L=JTJ+∑iri∇2fi gives ∥∇2L∥≤∥J∥2+n∥r∥max⁡i∥∇2fi∥. Integrating the Jacobian along the straight segment to initialization gives ∥r∥≤∥r0∥+TnD or ∥r0∥+CnD. Substitution of steps 3.1 and 4.1 proves the claimed loss-Hessian bound, after enlarging A,C to absorb squares of fixed exponential envelopes. For each coordinate the integral Taylor formula is fi(θ+h)−fi(θ)−Dfi(θ)h=∫01(1−t)D2fi(θ+th)[h,h] dt. Step 4.1 bounds its absolute value by TM(1+D)∥h∥2/(2m) or the corresponding residual constant. Summing coordinate squares gives the stated vector Taylor bound.

LemmaStatement: AI-adaptedProof: AI-generatedjudge pass (gpt-6-luna)Open item page →

Simultaneous initialization estimates for the Du architectures

Statement

Use either of the two architectures, its own terminal population kernel K, and independent Gaussian initialization. Assume AC, DC and countable choice and the activation and bounded-label hypotheses of the convergence theorems. Fix 0<δ<1. Write λ=λmin⁡(K)>0, q=1+log⁡(16/δ) and M=1+log⁡(16mnH/δ). There are constants A≥2 and C≥1, depending only on the fixed activation bounds, its normalizer, the label bound, and, in the residual case, cres, with the following property. For the fully connected case put T=AH; if

m≥Cmax⁡{Tn4/λ4,Tn/δ,Tn2log⁡(16Hn/δ)/λ2},

then with probability at least 1−δ, simultaneously:

max⁡i,h∥xi(h)(0)∥≤T,max⁡h≥2∥W(h)(0)∥op≤Cm,∥a(0)∥2≤Cm,∥a(0)∥4≤Cm1/4,

max⁡i,h∥bi,h(0)∥∞≤TM,λmin⁡(JH(0)JH(0)T)≥3λ/4,∥u(0)−y∥2≤Tnq.

Here JH is the Jacobian block of the last hidden-layer matrix, and bi,h is the derivative of the scalar output with respect to the hidden feature vector xi(h), so bi,H=a. In the residual case H≥2, the same conclusions hold with T replaced everywhere in the conclusions by C, provided

m≥Cmax⁡{n4H6/λ4,n2/(λ2H2),n/δ,n2log⁡(16Hn/δ)/λ2}.

For residual depth one use the fully connected case. Constants may be enlarged, and A may be increased, simultaneously in these conclusions.

Facts & Assumptions

Given: The architecture, its population recursion, the data, Gaussian initialization, width and confidence parameters in the Statement.

[F1]

The architecture-specific population recurrences, including residual means, and their positive semidefiniteness are supplied locally. Fully connected and residual networks with the Du Gram kernel The recursively defined population Gram matrices are finite, symmetric, and positive semidefinite

[F2]

Lipschitz functions of a standard Gaussian vector obey the stated one-sided concentration inequality; applying it to both signs gives a two-sided bound. Lipschitz concentration for standard Gaussian measure

[F3]

Gaussian vectors have the stated covariance and even moments (2p−1)!!. Multivariate normal law, including singular covariance Gaussian even moments for Brownian increments

[F4]

Conditional on hidden features, the Gaussian output norm is at most the square root of their total squared norm plus their operator norm times 2log⁡(1/ε), with failure probability at most ε. Gaussian output initialization concentrates conditionally on the hidden features

[F5]

If X,Y have finite second moments, then E∣XY∣≤(EX2)1/2(EY2)1/2; expanding a square and applying this bound gives the L2 triangle estimate used for residual feature variances. Cauchy-Schwarz for random variables

[F7]

A nondegenerate bivariate normal vector has its Gaussian Lebesgue density. A positive-definite multivariate normal law has its Gaussian density

[F8]

Differentiation may pass through an integral when parameter derivatives have a common integrable majorant. Differentiation under the integral sign

[F10]

Dominated convergence passes limits through integrals under a common integrable majorant. Dominated convergence

[F11]

For nonnegative Y and a>0, P(Y≥a)≤EY/a. Markov's inequality for random variables

[F12]

For square-integrable Y and a>0, P(∣Y−EY∣≥a)≤Var⁡(Y)/a2. Chebyshev's inequality for random variables

[F13]

A real function continuous on a compact interval and differentiable in its interior has some derivative equal to its average rate of change; this turns the bounded ρ-derivatives in step 2.1 into Lipschitz bounds. The mean value theorem, as the case g(x)=x of Cauchy's: for f continuous on [a,b] with a<b and differentiable on (a,b) there is c∈(a,b) with f(b)−f(a)=f′(c)(b−a)

[F14]

Analyticity gives σ′′, and the β-Lipschitz property of σ′ gives ∣σ′′∣≤β wherever needed in step 2.1. [given]

[F16]

A 1/4-net of the Euclidean unit sphere in Rm exists with at most 9m points. Finite epsilon nets of the Euclidean sphere

[F17]

Matrix operator norms are induced by the Euclidean vector norm. The operator norm is zero on the zero domain and otherwise is max_{||v||=1} ||Tv||

[F18]

For p=4, ∥a∥4=(∑r∣ar∣4)1/4. The p-norms ∥x∥p for rational p≥1, and ∥x∥∞

[F19]

The Frobenius norm is the Euclidean norm of the matrix entry list, ∥Φ∥F2=∑i,rΦir2; for any unit v, Euclidean Cauchy–Schwarz gives ∥Φv∥22=∑i(∑rΦirvr)2≤∑i,rΦir2, hence ∥Φ∥op≤∥Φ∥F. The Frobenius norm ∥A∥F=(∑i,j∣aij∣2)1/2 on real or complex matrices Cauchy-Schwarz ∣⟨x,y⟩∣≤∥x∥2∥y∥2 with its equality case, the triangle inequality for ∥⋅∥2, the parallelogram law and polarisation

[F15]

In the residual coupling, the feature vector and fresh Gaussian vector are independent, so Cov⁡(V+ασ(G))=Cov⁡(V)+α2Cov⁡(σ(G))⪰Cov⁡(V). This is the construction in The recursively defined population Gram matrices are finite, symmetric, and positive semidefinite and yields the monotone variance lower bound in step 4.1.

[F20]

For finitely many independent real random variables, expectations of products of nonnegative Borel functions factor, including when the expectations are extended nonnegative values. Expectations factor over finite products of independent random variables

[F21]

The probability of a finite union of events is at most the sum of their probabilities. The finite union bound

Proof

1.1F2F3F5F11F20givenalgebra

We give the scalar estimates used below explicitly. If U,V are centered Gaussian variables with variances at most a fixed Q, then ∣σ(z)∣≤B+L∣z∣ and [F3] imply ∥σ(U)∥2p≤Cp; [F5] then gives ∥σ(U)σ(V)∥p≤Cp for integers p≥2. The same last bound holds for a2σ′(U)σ′(V) when a is an independent standard normal, because ∣σ′∣≤L. After centering, each such variable X satisfies E∣X∣p≤p!Cp, enlarging C using pp≤epp!. Expanding its exponential series, for ∣t∣≤1/(2C) one has EetX≤1+∑p≥2(C∣t∣)p≤e2C2t2. For m independent copies, [F20] factors the exponential moment of their sum. Applying [F11] to exp⁡(t∑rXr) with t=s/(4C2), and then to exp⁡(−t∑rXr), to obtain P(∣m−1∑rXr∣>s)≤2e−cms2 for 0<s≤c. All constants are uniform over the covariances with these variance bounds.

2.1F3F5F7F8F9F10step 1.1algebraF13

The Gaussian expectation maps needed in the population recurrences are Lipschitz functions of the three entries of their bivariate covariance on any region v≤Qii,Qjj≤Q, where v>0 is fixed. This applies to E[g(U)g(V)] for g=σ and g=σ′, and to the residual marginal map s↦Eσ(sZ). At fixed ρ, changing s,t changes these expectations by at most C(∣Δs∣+∣Δt∣): use the Lipschitz constants L,β, s,t≥v, and [F5] to bound the resulting products. For fixed s,t and ∣ρ∣<1, the covariance matrix is positive definite, so [F7] supplies the density ϕs,t,ρ. On any compact interval ∣ρ∣≤r<1, its eigenvalues lie between v(1−r) and Q(1+r). Direct differentiation of the Gaussian density gives ∂ρϕs,t,ρ=st ∂u∂vϕs,t,ρ. Uniform eigenvalue bounds imply that ∣g(u)g(v)∂ρϕs,t,ρ(u,v)∣ is bounded by an integrable polynomial times e−c(u2+v2), uniformly on this interval; g has at most linear growth and g′ is bounded. Thus [F8] justifies differentiating the expectation under the integral. On finite rectangles apply [F9] successively in u and v; the boundary terms tend to zero as the rectangles expand because of Gaussian decay and the growth bound. The same integrable majorant and [F10] pass the rectangle limits, giving ∂ρE[g(U)g(V)]=stE[g′(U)g′(V)]. Its absolute value is at most QL2 for g=σ and Qβ2 for g=σ′, since ∣σ′∣≤L and ∣σ′′∣≤β. At ρ=±1, the representation U=sZ, V=t(ρZ+1−ρ2Z′) is dominated after multiplication by either pair of activations by C(1+∣Z∣+∣Z′∣)2, which is integrable; [F10] gives continuity at both endpoints. Finally ρ=Qij/st is Lipschitz in the covariance entries when s,t≥v.

3.1F1F3F21step 1.1step 2.1algebra

Write S^ij(h)=⟨xi(h)(0),xj(h)(0)⟩ and, for residual layers, b^i(h)=m−1/2∑rxi,r(h)(0). Conditional on previous layers, the preactivation pairs in distinct rows are independent centered Gaussians with covariance given by S^(h−1). In the fully connected case, step 1.1 bounds each difference of S^ij(h) from its conditional mean by Ct0, except with probability δ/(32Hn2), where t0=log⁡(64Hn2/δ)/m. Its conditional mean is exactly the fully connected population expectation evaluated at S^(h−1). Since the population diagonal is one, step 2.1 gives, while the error is at most a fixed small constant, eh≤Ceh−1+Ct0, where eh=max⁡ij∣S^ij(h)−Kij(h)∣ and e0=0. If B>1 is a fixed base for this recurrence, then eh≤BHt0 after enlarging its prefactor. Choose the statement constant A≥B2 and enlarge the width constant; its term Tn2log⁡(16Hn/δ)/λ2, with T=AH, then makes this error smaller than both the fixed stopping threshold and cλ/n. Conditional bounds and the finite union bound [F21] over h,i,j therefore justify the estimates through layers h<H; the last hidden-feature norm is bounded by one further conditional application of the activation moment estimate, without identifying its covariance with the terminal derivative kernel. A first layer crossing the stopping threshold would contradict this recurrence. In particular all feature norms are bounded by a constant, hence by T.

4.1F1F2F3F21step 1.1step 2.1step 3.1algebraF15

In the residual case the population feature variances lie in a fixed compact interval [v,Q] independent of H. Indeed the vector construction in [F1] gives S(h)−b(h)(b(h))T⪰S(1)−b(1)(b(1))T, whose diagonal is v=cσVar⁡(σ(Z))>0. Positivity of v follows from nonconstancy, continuity and the positive Gaussian density. Also Sii(h)≤Sii(h−1)+α(B+LSii(h−1)) with α=cres/H, so iteration gives a fixed upper bound. The empirical conditional mean satisfies exactly the same residual second-moment/mean recursion with S,b replaced by S^,b^. The nonlinear quadratic part has conditional fluctuation at most Cα2t0 by step 1.1. Each cross part is αm−1/2∑rxi,r(h−1)σ(Uj,r); its centered value has absolute size at most Cαt0 by [F2], since its Lipschitz constant in the independent row Gaussians is at most Cα∥xi(h−1)∥/m. The empirical mean has fluctuation at most Cαt0 by [F2]. Step 2.1 now yields eh≤(1+C/H)eh−1+Ct0/H, taking the maximum of the covariance and mean errors. The first layer has e1≤Ct0 by step 1.1 and [F2]. Hence eh≤Ct0 uniformly in H. The width assumption and the same stopping argument, now in [v/2,2Q], justify the conditional estimates and give the feature bound through h<H; the finite union bound [F21] combines their layerwise failure probabilities. One further conditional step gives ∥xi(H)∥≤C, using the same upper-norm recurrence, so max⁡i,h∥xi(h)∥≤C.

5.1F1F3F6F21step 1.1step 2.1step 3.1step 4.1algebra

For either architecture, conditional on layers before the last, the last-hidden-layer Jacobian Gram entries are sH2S^ij(H−1)m−1∑rar2σ′(zi,r)σ′(zj,r), where sH2=cσ for the fully connected case and sH2=cres2/H2 for residual H≥2. This follows by differentiating the last layer. The summands are independent across r. Steps 1.1–2.1, the feature bounds and errors just proved give entry error at most B2Ht0 in the fully connected case, for a fixed base B2 depending only on the fixed activation bounds, and at most Ct0/H2 in the residual case. Choose the statement constant A also with A≥B22; the displayed logarithmic width term and a sufficiently large fixed multiplier then make both errors at most λ/(4n). A finite union bound [F21] covers the terminal-kernel entry events. If an n×n matrix E has entries bounded by ε, then for any z, ∥Ez∥2≤n ε∑j∣zj∣≤nε∥z∥2 by [F6], hence ∥E∥op≤nε. Thus ∥JHJHT−K∥op≤λ/4. Both matrices are symmetric, and for every z, zTJHJHTz≥zTKz−∣zTEz∣≥(λ−∥E∥op)∥z∥22≥3λ∥z∥22/4, where E=JHJHT−K. The Rayleigh-quotient characterization of the smallest eigenvalue now gives λmin⁡(JHJHT)≥3λ/4.

5.2F2F3F6F12F16F21step 3.1step 4.1constructalgebraF17F18

By [F16], choose a 1/4-net N of the unit sphere with at most 9m points. The Euclidean norm duality ∥z∥2=sup⁡∥x∥2=1∣xTz∣ follows from [F6], so ∥W∥op=sup⁡∥x∥=∥y∥=1∣xTWy∣. For unit x,y, choose v,w∈N within 1/4 of them. Then [F6] gives ∣xTWy−vTWw∣≤∥x−v∥2∥Wy∥2+∥v∥2∥W(y−w)∥2≤12∥W∥op. Taking suprema yields ∥W∥op≤2max⁡v,w∈N∣vTWw∣. Each vTWw is standard normal, whose tail is at most 2e−t2/2 by [F2]. The finite union bound [F21] gives ∥W(h)∥≤Cm for all h≥2 with failure at most δ/16 when m≥Clog⁡(16H/δ), which follows from either displayed width because λ≤cσL2 in the fully connected case and λ≤C/H2 in the residual case. The latter bound follows by bounding the diagonal of its terminal kernel with step 4.1. Gaussian second and eighth moments and [F12] give ∥a∥2≤Cm and ∥a∥4≤Cm1/4 with failure at most δ/16 when m≥C/δ, also implied by the widths. These arguments never bound ∥W(1)∥ on the whole input space, so no input-dimension condition is used.

6.1F2F3F4F21step 3.1step 4.1step 5.1step 5.2F19

Conditional on hidden weights, bi,h(0)=Pi,hTa, where Pi,h is the product of downstream feature Jacobians. Their operator norms have fixed envelopes BH in the fully connected case and B in the residual case, using ∏l>h(1+C/H)≤eC for residual layers. Increase the conclusion constants so these envelopes are at most T/2 and C/2, respectively. Each coordinate is centered Gaussian of variance at most (T/2)2 or (C/2)2. Since M2=1+log⁡(16mnH/δ), its two-sided tail is at most 2e−2M2≤δ/(8mnH); the finite union bound [F21] over at most mnH coordinates gives the asserted coordinate maxima with failure at most δ/8. Finally [F4] and ∥Φ∥op≤∥Φ∥F≤nmax⁡i∥xi(H)∥ give ∥u(0)∥≤CTnq or Cnq with failure at most δ/16.

7.1F2F21step 3.1step 4.1step 5.1step 5.2step 6.1algebra∎

By the finite union bound [F21], allocate at most δ/(32Hn2) to each fully connected feature-covariance or terminal-feature-norm event in step 3.1; there are at most 2Hn2 such events. Allocate at most δ/(128Hn2) to each residual covariance, cross-term, mean, or terminal-feature-norm event in step 4.1; there are at most 4Hn2 such events. Allocate at most δ/(64n2) to each terminal-kernel entry event. The weight-matrix and readout norm events in step 5.2 each use δ/16; backward coordinates use δ/8 and the output norm uses δ/16 in step 6.1. Even summing both architecture allocations, these budgets total at most 27δ/64<δ. The thresholds are constant multiples of t0 or M; increasing the fixed constants in the width and conclusions makes each displayed concentration bound meet its allocation without changing their forms. Adding the deterministic bounded-label contribution proves the simultaneous event.

TheoremStatement: Literature-sourcedProof: AI-generatedjudge pass (gpt-6-luna)Open item page →

Wide fully connected networks converge geometrically to zero empirical training loss

Statement

Assume AC, DC and countable choice. Fix positive integers d,n,H,m, training data x1,…,xn∈Rd and labels y1,…,yn∈R with ∥xi∥2=1 and ∣yi∣≤Cy for a fixed finite constant Cy. Let the activation satisfy the Lipschitz and smoothness conditions ∣σ(0)∣≤Cσ, ∣σ(z)−σ(z′)∣≤L∣z−z′∣, and ∣σ′(z)−σ′(z′)∣≤β∣z−z′∣ for finite constants, and also be analytic and nonpolynomial. Use the fully connected architecture and population matrix Kfc(H) defined on this page, initialized with independent standard normal entries in all weight matrices and the output vector, and trained in all layers by full-batch gradient descent on L(θ)=12∑i=1n(f(xi,θ)−yi)2. Write K(H)=Kfc(H) and put λ0=λmin⁡(K(H)) and assume λ0>0. For 0<δ<1, if the width and fixed step size satisfy

m=Ω ⁣(2O(H)max⁡{n4λ04,nδ,n2log⁡(Hn/δ)λ02}),0<η=O ⁣(λ0n22O(H)),

then with probability at least 1−δ over initialization, every gradient-descent iterate k≥1 satisfies

L(θ(k))≤(1−ηλ02)kL(θ(0)).

The constants hidden in O and Ω may depend on the fixed activation and label bounds but not on n,H,m,δ; the step-size constant is small enough that 0<ηλ0/2<1. In particular, L(θ(k))→0 geometrically as k→∞.

Facts & Assumptions

Given: The data, architecture, activation, Gaussian initialization, positive terminal population eigenvalue λ0, width, confidence and step-size hypotheses of the Statement.

[F1]

The architecture-specific population kernel, including residual means and scaling, is defined locally and is positive semidefinite. Fully connected and residual networks with the Du Gram kernel The recursively defined population Gram matrices are finite, symmetric, and positive semidefinite

[F2]

The simultaneous Gaussian initialization event gives feature, parameter, backward-coordinate, last-layer Gram, and initial-residual bounds at the stated width. Simultaneous initialization estimates for the Du architectures

[F3]

On parameter balls, the last-layer Gram perturbation and output/loss Hessian have explicit width-dependent bounds. Parameter-ball control of the last-layer Gram and loss Hessian

[F4]

Local Gram positivity and loss-Hessian control confine all GD iterates by a first-exit argument and give geometric loss contraction. Local Gram positivity and loss smoothness give confined gradient descent

Proof

1.1F1F2givenalgebra

Put λ=λ0 and q=1+log⁡(16/δ). Fix an exponential envelope E=A0H for all initialization constants in [F2] and another F=A1H for [F3] with those initialization bounds. These bases depend only on the permitted fixed bounds. The source-width form permits m≥Wmax⁡{n4/λ4,n/δ,n2log⁡(16Hn/δ)/λ2}, where W=A2H is chosen sufficiently large in terms of the fixed bases of E,F; in particular take W≥CEF2. All these choices are fixed before selecting the radius. On the event from [F2], of probability at least 1−δ, r0=u(0)−y satisfies ∥r0∥2≤Enq, the initialized last-layer Gram is at least 3λI/4, and all other hypotheses of [F3] hold. Also λ≤cσL2 by the diagonal of [F1].

2.1F1F3step 1.1algebra

The displayed width also implies m≥CEF2n3q/λ3, after fixing W sufficiently large. For the logarithmic term this follows without an extra logarithmic width factor: writing A∗=n4/λ4 and B∗=n/δ, one has n3λ−3log⁡(1/δ)=n−1/4A∗3/4B∗1/4δ1/4log⁡(1/δ)≤(4/e)max⁡(A∗,B∗), since te−t/4≤4/e for t≥0. The constant part of q is absorbed by A∗ using λ≤cσL2 and n≥1. Choose the fixed radius D=4Enq/λ. These inequalities give D/m≤c/F and FnD/m≤λ/4. Hence [F3] supplies JH(θ)JH(θ)T⪰λI/2 on the entire ball of radius D.

3.1F3step 1.1step 2.1algebra

The same width implies m≥Cmax⁡(n,H,1/δ), so M=1+log⁡(16mnH/δ)≤Clog⁡(em) and q≤Clog⁡(em). The loss-Hessian bound in [F3], D=4Enq/λ, and λ≤cσL2 yield ∥∇2L∥≤T′n2λ−1(1+qM/m), where T′=A′H after enlarging its fixed base. Indeed expand (Enq+FnD)(1+D) in that bound; because n,q≥1 and λ is bounded above, every term is at most an allowed exponential constant times n3/2q/λ, before the exterior factor nM/m. The function log⁡(em)3/2/m is bounded for m≥1, so the Hessian is bounded by S=T′′n2/λ on the ball, with T′′=A′′H. The allowed sufficiently small step-size constant ensures η≤1/S.

4.1F3F4step 1.1step 2.1step 3.1algebra∎

The radius satisfies D>22∥r0∥/λ. Apply [F4] with the last-hidden-layer Jacobian block and the ball and smoothness bounds from steps 2.1 and 3.1. Every iterate stays in that ball, and L(θ(k))≤(1−ηλ/2)kL(θ(0)) for all k. Choosing the step-size constant also ensures 0<ηλ/2<1, so geometric convergence to zero follows. The initialization event has probability at least 1−δ, proving the claim. The same application bounds the total parameter path by 22∥r0∥/λ; [F3] therefore supplies last-layer Gram stability and the displayed output Taylor-remainder bound on every accepted update segment, for all iterates simultaneously.

TheoremStatement: Literature-sourcedProof: AI-generatedjudge pass (gpt-6-luna)Open item page →

Wide residual networks converge geometrically to zero empirical loss with polynomial depth dependence

Statement

Assume AC, DC and countable choice. Fix positive integers d,n,H,m, training data x1,…,xn∈Rd and labels y1,…,yn∈R with ∥xi∥2=1 and ∣yi∣≤Cy for a fixed finite constant Cy. Let the activation satisfy the Lipschitz and smoothness conditions ∣σ(0)∣≤Cσ, ∣σ(z)−σ(z′)∣≤L∣z−z′∣, and ∣σ′(z)−σ′(z′)∣≤β∣z−z′∣ for finite constants, and also be analytic and nonpolynomial. Use the residual architecture with fixed 0<cres<1 and the population matrix Kres(H) defined on this page, initialized with independent standard normal entries in all weight matrices and the output vector, and trained in all layers by full-batch gradient descent on L(θ)=12∑i=1n(fres(xi,θ)−yi)2. Write K(H)=Kres(H) and put λ0=λmin⁡(K(H)) and assume λ0>0. For 0<δ<1, if

m=Ω ⁣(max⁡{n4H6λ04,n2λ02H2,nδ,n2log⁡(Hn/δ)λ02}),0<η=O ⁣(λ0H2n2),

then with probability at least 1−δ over initialization, every gradient-descent iterate k≥1 satisfies

L(θ(k))≤(1−ηλ02)kL(θ(0)).

The constants hidden in O and Ω may depend on the fixed activation, label, and residual-scaling bounds but not on n,H,m,δ; the step-size constant is small enough that 0<ηλ0/2<1. In particular, L(θ(k))→0 geometrically as k→∞.

Facts & Assumptions

Given: The data, architecture, activation, Gaussian initialization, positive terminal population eigenvalue λ0, width, confidence and step-size hypotheses of the Statement.

[F1]

The architecture-specific population kernel, including residual means and scaling, is defined locally and is positive semidefinite. Fully connected and residual networks with the Du Gram kernel The recursively defined population Gram matrices are finite, symmetric, and positive semidefinite

[F2]

The simultaneous Gaussian initialization event gives feature, parameter, backward-coordinate, last-layer Gram, and initial-residual bounds at the stated width. Simultaneous initialization estimates for the Du architectures

[F3]

On parameter balls, the last-layer Gram perturbation and output/loss Hessian have explicit width-dependent bounds. Parameter-ball control of the last-layer Gram and loss Hessian

[F4]

Local Gram positivity and loss-Hessian control confine all GD iterates by a first-exit argument and give geometric loss contraction. Local Gram positivity and loss smoothness give confined gradient descent

[F5]

For square-integrable real random variables, ∣EX∣≤(EX2)1/2. Cauchy-Schwarz for random variables

Proof

1.1F1F2F3F5givenalgebra

First take H≥2 and put λ=λ0, q=1+log⁡(16/δ). At the stated width, [F2] supplies with probability at least 1−δ a simultaneous event with all constants independent of H, on which ∥r0∥2≤Cnq and JH(0)JH(0)T⪰3λI/4. By [F1], Sii(1)=1 and Sii(h)−(bi(h))2≥0. For 2≤h<H, put qh=Sii(h) and Yh=σ(Gi(h−1)). The residual recursion gives qh2=Sii(h−1)+2αbi(h−1)EYh+α2EYh2. By [F1] and [F5], ∣bi(h−1)∣≤qh−1 and ∣EYh∣≤(EYh2)1/2; Lipschitz growth and [F5] also give (EYh2)1/2≤Cσ+Lqh−1. Hence qh≤(1+αL)qh−1+αCσ. Starting from q1=1 and iterating with α=cres/H gives qh≤eLcres(1+Cσcres) for every h<H, uniformly in H. The terminal diagonal formula and ∣σ′∣≤L now give Kii(H)≤C/H2; positive semidefiniteness in [F1] implies λ≤Kii(H)≤C/H2. The residual neighborhood bounds [F3] retain separate first-layer and readout gradients without an erroneous 1/H factor.

2.1F1F3step 1.1algebra

Choose D=4Cnq/λ, with C fixed sufficiently large for the initial-residual bound. The width implies m≥Cn3q/(λ3H4) and m≥CD2. To see the logarithmic absorption, use A∗=n4/λ4, B∗=n/δ and n3λ−3log⁡(1/δ)≤(4/e)max⁡(A∗,B∗), which follows from δ1/4log⁡(1/δ)≤4/e exactly as in the fully connected argument. Dividing by H4 only decreases the requirement, while the approved width contains H6A∗ and B∗. The constant part is absorbed using λH2≤C and n,H≥1. Consequently D/m≤c and CnD/(H2m)≤λ/4, so the last-layer Gram is at least λI/2 throughout this ball.

3.1F3step 1.1step 2.1algebra

All constants in the residual Hessian estimate are independent of H. Substitution of D and the initial-residual bound in [F3] yields ∥∇2L∥≤Cn+Cn2qM/(λm), because λ≤C and n,q≥1. The approved width and λ≤C/H2 imply m≥Cmax⁡(n,H6,1/δ). Thus M≤Clog⁡(em), q≤Clog⁡(em) and qM/m≤C/H2: explicitly, log⁡(em)3/2≤Cm1/6 for all m≥1, so the ratio is at most Cm−1/3≤C/H2. Also n≤Cn2/(λH2). Therefore ∥∇2L∥≤S=Cn2/(λH2) on the ball. The permitted sufficiently small fixed step-size constant makes η≤1/S.

4.1F1F2F3F4step 1.1step 2.1step 3.1algebra∎

Since D>22∥r0∥/λ, [F4] confines every iterate to the ball and gives L(θ(k))≤(1−ηλ/2)kL(θ(0)). The step-size constant can ensure 0<ηλ/2<1, because λH2≤C. This proves geometric convergence on the probability-1−δ event. For H=1 the architectures and population kernels coincide; the width is exactly the fully connected depth-one restriction up to fixed constants, and its local proof applies with the same contraction and step-size form. The same application bounds the total parameter path by 22∥r0∥/λ; [F3] therefore supplies last-layer Gram stability and the displayed output Taylor-remainder bound on every accepted update segment, for all iterates simultaneously.

Remarks

The displayed sufficient width retains the approved stronger depth restriction. The first term in the official Theorem 6.1 is n4/(λ04H6), whereas the displayed term is n4H6/λ04. The proof here establishes the displayed approved restriction by a parameter-ball argument, retaining all trainable layers. It does not claim the sharper source width from this argument. The residual population kernel includes its mean cross terms and differs from the fully connected kernel.

4 · Examples, counterexamples and false statements

None yet.

Sources