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
- Absolute and Conditional Convergence; Rearrangement; Products
- Absolute Continuity and the Sharp Fundamental Theorem of Calculus
- Algebraic Closure, Embeddings, and Separability
- Algebraic Extensions, Extension Degree, and Finite Fields
- Areas of Elementary Plane Figures
- Automatic differentiation and backpropagation
- Automatic differentiation and computational graphs
- Banach Gaussian priors, RKHS support, and concentration
- Binary Operations, Monoids, Groups and Subgroups
- Bounded Variation and the Riemann–Stieltjes Integral
- Brownian Motion Construction and Continuity
- Central Limit Theorems
- Characteristic Functions Inversion and Continuity
- Compactness
- Compactness in Metric Spaces
- Complete Metrizability, Čech-Completeness, and Baire Category
- Completeness, Completion, and Uniform Continuity
- Complex Lp Spaces and Test-Function Conventions
- Composition Series, the Jordan–Hölder Theorem and Solvable Groups
- Computational foundations: factorization, KL, and optimization
- Concentration
- Congruences, the Integers Modulo n and the Chinese Remainder Theorem
- Construction of the Natural Numbers
- Construction of the Real Numbers via Cauchy Sequences
- Construction of the Real Numbers via Dedekind Cuts
- Continuity, IVT, EVT, and Uniform Continuity
- Cosets, Index and Lagrange's Theorem
- Countability and Uncountability
- Countability Axioms and Cardinal Functions
- Cyclic Groups and Direct Products
- Density Separability and Convolution in Lᵖ
- Determinants of Matrices over a Commutative Ring
- Diagonalisation and the Minimal Polynomial
- Differentiation of Monotone Functions and the Vitali Covering Theorem
- Direct Matrix Factorisations: LU, Cholesky and QR
- Divisibility, Euclidean Domains, Principal Ideal Domains and Unique Factorisation
- Divisibility, Greatest Common Divisors and Bézout's Identity
- Dual Spaces, Bilinear and Quadratic Forms, and Sylvester's Law of Inertia
- Eigenvalues, Eigenvectors and the Characteristic Polynomial
- Feedforward networks and single-hidden-layer universality
- Filters and Ultrafilters
- Finite Counting, Factorials and Binomial Coefficients
- Finite Probability and the Probabilistic Method
- Finite Probability Spaces and Random Variables
- Foundations of the Real Numbers for Analysis
- Fourier Transform Convolution and Approximate Identities
- Fubini and Change of Variables
- Fundamental Trigonometric Identities
- Gaussian Elimination, Elementary Matrices and Reduced Row Echelon Form
- Group Actions, Orbits, Stabilisers and Cayley's Theorem
- Group Homomorphisms and the Isomorphism Theorems
- Hausdorff via the Diagonal
- Ideals, Quotient Rings and the Isomorphism Theorems for Rings
- Improper and Parameter-Dependent Multiple Integrals
- Improper Integrals
- Independence Borel Cantelli and Zero One Laws
- Infinite Product Measures and Kolmogorov Extension
- Inner Product Spaces, Gram-Schmidt, Projections and Adjoints
- Lebesgue Measure on Euclidean Space
- Limits of Real Functions
- limsup, liminf, and Subsequential Limits
- Linear Independence, Bases and Dimension
- Linear Transformations, Rank-Nullity and Quotient Spaces
- Matrices, the Matrix of a Linear Map, and Change of Basis
- Matrix Differentiation and First-order Spectral Perturbation
- Matrix Norms, Condition Numbers and Numerical Stability
- Measurable Functions and Simple Approximation
- Measures and Their Basic Properties
- Metric Spaces
- Mixed effects and longitudinal models
- Mixed Partials, Taylor Formulae, and Extrema
- Modes of Convergence Egorov and Lusin
- Monotone Functions, Discontinuities, and Continuity Sets
- Monotone Sequences, Bolzano-Weierstrass, and Cauchy Completeness
- Normal Subgroups and Quotient Groups
- Order, Zorn's Lemma, and the Axiom of Choice
- Outer Measure and the Caratheodory Extension Theorem
- Outward rounded interval and ball arithmetic
- Partitions of Unity and Paracompactness
- Pointwise Extended-Real Function Products
- Polynomial Rings, the Division Algorithm and Roots
- Power Series and Real-Analytic Functions
- Primes, Euclid's Lemma and the Fundamental Theorem of Arithmetic
- Probability Spaces Random Variables and Expectation
- Product Measures and the Fubini Tonelli Theorems
- Properties of the Integral and the Working FTC
- Radon Measures and the Riesz Markov Kakutani Theorem
- Relations, Functions, and Quotients
- Rings, Subrings, Integral Domains and Fields
- Rⁿ as a Normed Space; Vector-Valued Functions
- Roots, Rational Powers, and Classical Inequalities
- Separation Axioms: the Hierarchy
- Sequences and Limits
- Sequences and Series of Functions; Uniform Convergence
- Series: Convergence and the Nonnegative Tests
- Sigma Algebras and Borel Sets
- Signed and Complex Measures Hahn and Jordan
- Simple Field Extensions and the Construction of the Complex Numbers
- Sine, Cosine, and the Definition of Pi
- Splitting Fields
- State-space model laws
- Subspaces, Products, and Quotients
- Suprema and Infima
- Sylow's Theorems, p-Groups and Nilpotent Groups
- Symmetric Groups, Cycle Decomposition and the Sign Homomorphism
- The Cantor Set, Baire Category, and Measure Zero in ℝ
- The Complex Exponential and Euler's Formula
- The Derivative and the Mean Value Theorems
- The Determinant of a Linear Operator, Cofactors and Cramer's Rule
- The Exponential Function
- The Fundamental Theorem of Algebra
- The Fundamental Theorem of Finite Abelian Groups
- The Galois Correspondence
- The Inverse and Implicit Function Theorems
- The Lebesgue and Riemann Integrals Compared
- The Lebesgue Integral and the Convergence Theorems
- The Logarithm and General Powers
- The Lᵖ Spaces Holder Minkowski and Riesz Fischer
- The Maximal Function and Lebesgue Differentiation
- The Radon Nikodym Theorem and Lebesgue Decomposition
- The Riemann Integral in Rᵐ and Jordan Content
- The Riemann Integral: Definition and Integrability
- The Spectral Theorem, Positive Operators and Singular Value Decomposition
- The Topology of Euclidean Space
- The Total Derivative in ℝᵐ → ℝⁿ
- The ZFC Axioms and the Basic Set Constructions
- Topological Spaces and Continuity
- Topology of ℝ
- Triangularisation, Generalised Eigenspaces and Jordan Canonical Form
- Urysohn's Lemma and the Tietze Extension Theorem
- Vector Spaces, Linear Subspaces, Span and Direct Sums
- Volumes of Elementary Solids and Solids of Revolution
- Weak Convergence Tightness and Representation
2 · Logical flowchart
3 · Definitions, theorems and proofs
The Gaussian activation normalizer is finite and positive
Statement
Assume AC. Let have the standard normal law on . Suppose is analytic, nonpolynomial, satisfies for some finite , and is -Lipschitz for a finite . Then . Consequently the normalization constant
is a well-defined positive finite real number.
Statement references
Facts & Assumptions
Given: AC, the activation , its constants , and a standard normal random variable .
AC is assumed by the standard-normal-law and Gaussian-moment suppliers. The Axiom of Choice
The standard normal law has second moment . Standard normal and normal laws Gaussian even moments for Brownian increments
Proof
Lipschitz continuity gives , hence . By [F1], .
A nonpolynomial analytic function is not identically zero, so choose with . Continuity gives an open interval around on which . The standard normal density is strictly positive on that interval, so its integral of is positive.
The expectation is both finite and strictly positive by steps 1.1 and 1.2, so its reciprocal exists and is positive and finite.
Gram control and a Taylor remainder imply one-step loss contraction
Statement
Let , let be a real symmetric matrix with eigenvalues in for some , and suppose . If an update has the form with , then for the quadratic loss and ,
Statement references
Facts & Assumptions
Given: The residual vector, symmetric Gram matrix, eigenvalue bounds, step size, and remainder bound.
A real symmetric matrix has an orthonormal eigenbasis. Real spectral theorem: a self-adjoint endomorphism of a finite-dimensional real inner product space has an orthonormal eigenbasis
The Euclidean norm satisfies the triangle inequality. Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation
Proof
By [F1], the eigenvalues of lie in , so those of lie in . Expanding in an orthonormal eigenbasis gives . The norm triangle inequality [F2] with the remainder hypothesis then yields .
Put . Then , since . Squaring the estimate in step 1.1 and multiplying by yields .
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 , training inputs with , labels , an analytic nonpolynomial activation satisfying and finite Lipschitz constants for and , and the positive finite normalization . A fully connected network of width and depth is the finite feedforward recurrence
where , for , , and acts coordinatewise. The residual network uses the same first layer and then
with . In either architecture, random initialization means all entries of all matrices and of are independent random variables; all these entries are trained by full-batch gradient descent. The empirical quadratic loss is .
For the fully connected architecture and the same fixed data, define the population matrices by
For , let be a centered bivariate normal vector with covariance
and put . For , use the same bivariate-normal definition with formed from , and set
When , this last formula uses directly and there are no intermediate activation-covariance matrices. For the fully connected theorem, this terminal matrix is denoted .
For the residual architecture with , use a separate sequence and mean vector . Set . For , put
Write . For , take a fresh centered Gaussian vector with covariance , put , and define
For , define the residual terminal matrix by
The fully connected and residual convergence statements use and , respectively. At , the two architectures coincide; define . 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 , as determined by the finite-width last-layer Gram formula (14).
The residual recurrence has one identity contribution from and one nonlinear contribution at every layer ; 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}Definition references
Du et al., fully connected and residual architectures, equations (2)–(3)
Du et al., population Gram recursions, Definitions 5.1 and 6.1, equations (7) and (8)
Du et al., Appendix C, last-layer Gram formula (14)
Du et al., Appendix E.1, general population mean and second-moment recursion (16)–(17)
Finite feedforward computation graphs and layered perceptrons
Princeton Theory of Deep Learning, Chapter 8, §§8.1–8.5, infinite-width context
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 , every residual feature second-moment matrix , and both terminal matrices are finite, symmetric, and positive semidefinite. The fully connected matrices satisfy for . For the residual sequence, is positive semidefinite for ; its diagonal is not asserted to equal one.
Statement references
Facts & Assumptions
Given: The finite unit-norm data, activation, normalization, Gaussian covariance convention, and matrices from the definition.
The architecture and matrices are defined by the stated recurrences, with and finite. Fully connected and residual networks with the Du Gram kernel
The normalizer satisfies for standard normal . The Gaussian activation normalizer is finite and positive
A property holding at zero and preserved by successor holds for every natural-number index. The principle of mathematical induction
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
A real self-adjoint matrix has an orthonormal eigenbasis. Real spectral theorem: a self-adjoint endomorphism of a finite-dimensional real inner product space has an orthonormal eigenbasis
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
Assume AC as required by the multivariate normal-law supplier. The Axiom of Choice
Assume DC and countable choice as required by the independent-copy supplier. The axiom of dependent choice: a relation in which every element is related to something admits an -indexed chain The Axiom of Countable Choice ()
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
If have finite second moments, then . 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
The input matrix is symmetric, and for , . Its diagonal entries are , giving the base case.
Assume the induction hypothesis that is finite, symmetric positive semidefinite, and has diagonal , for some . By [F2]–[F3], take ; each pair has covariance . The growth bound and finite Gaussian second moments put each in ; [F8] then makes every product integrable. Thus is finite and symmetric, and for every , . Because , its diagonal is by [F5]. This establishes the induction step.
For the fixed , let mean that is finite, symmetric positive semidefinite, and has diagonal , for . Define to mean for every . Step 1.1 gives . If and , step 2.1 applied with gives and hence ; if , adds no new index. By [F6], holds, so is finite, symmetric positive semidefinite, and has diagonal ; when , this is exactly the base matrix . By [F4], take independent both distributed as . The vector with coordinates has finite second moments; is bounded because it is the derivative of a Lipschitz differentiable function. Factoring expectations of independent variables gives . Therefore is finite and symmetric, and for every . This proves the fully connected terminal claim and discharges that induction.
For the residual base layer let and , with activation applied coordinatewise. Then and by [F1]. The growth bound used in step 2.1 proves finite second moments. For every , , while is its covariance and has quadratic form . Thus both matrices are symmetric positive semidefinite.
Suppose a finite-second-moment vector has mean and second-moment matrix . Let be the Borel law of this generally non-Gaussian vector and let , which exists by [F2]. Both are probability measures and hence sigma-finite. By [F7], on the coordinate projections and have laws and and are independent: for Borel sets , the probability of is . Put . Finite second moments follow from the growth bound and . Independence gives and . For every , , and is the covariance of , 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].
For , take independent and set . The boundedness gives finite second moments, and independence gives . Consequently this terminal matrix is symmetric positive semidefinite. At it equals the already proved fully connected matrix by definition. This proves all residual assertions.
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 , let , let have row equal to the gradient of with respect to the vector of all trainable parameters, and set . For the empirical quadratic loss , one full-batch gradient-descent step with step size satisfies
where and is the first-order Taylor remainder of the differentiable map at in direction . No size bound on is asserted here.
Statement references
Facts & Assumptions
Given: One finite network, its trainable parameter vector, a fixed data set and labels, and a step size .
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
Reverse accumulation through a finite parameterized recurrence gives the exact parameter gradient, including all parameter uses. Reverse accumulation sums gradients over shared parameter uses
Fixed-step gradient descent updates the parameter by subtracting times the Euclidean gradient. Gradient descent
Differentiability in finite-dimensional real parameter spaces means a linear first-order approximation with a remainder . The real Frechet derivative on real and complex matrix spaces with the Frobenius norm
Proof
Write . By the scalar chain rule applied to , the parameter gradient is . By [F3], .
Since is differentiable, [F4] gives , where is the first-order remainder. Substitute to get , which is the claimed identity.
Gaussian output initialization concentrates conditionally on the hidden features
Statement
Assume AC, DC and countable choice. Fix hidden-layer feature vectors , form the matrix with rows , and let the independent output-weight vector . Put and . Then for every ,
where . The probability is conditional on the fixed hidden-layer features and is over only.
Statement references
Facts & Assumptions
Given: AC, DC, countable choice, fixed feature vectors, and an independent standard Gaussian output vector.
AC, DC and countable choice are assumed by the Gaussian concentration supplier. The Axiom of Choice The axiom of dependent choice: a relation in which every element is related to something admits an -indexed chain The Axiom of Countable Choice ()
For a standard Gaussian vector and an -Lipschitz real function , for ; for the upper-tail event is empty. Lipschitz concentration for standard Gaussian measure
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
Cauchy–Schwarz bounds by the product of the norms. Cauchy-Schwarz for random variables
The Euclidean norm satisfies the triangle inequality and hence the reverse triangle inequality. Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation
Proof
[F4] gives , so is -Lipschitz. Since by [F2], ; [F3] with and yields .
If , then and the asserted event has probability one. Otherwise apply [F1] to with . The upper-tail probability is at most , and the mean bound in step 1.1 gives the displayed threshold.
Local Gram positivity and loss smoothness give confined gradient descent
Statement
Let be , , and . On the closed ball suppose a block of the output Jacobian satisfies and , with . If and , every fixed-step gradient-descent iterate stays in the ball,
The assertion is trivial for ; otherwise the hypotheses imply .
Statement references
Facts & Assumptions
Given: The output map, quadratic loss, ball, spectral and Hessian bounds, and step size.
Gradient descent uses . Gradient descent
For a scalar 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
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||
Euclidean Cauchy–Schwarz gives for all real vectors of the same finite dimension. Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation
Proof
Write . Since , in the ball. For any step , , whose segment stays in the ball, put for . Then and , where . By Euclidean Cauchy–Schwarz [F4] and the induced norm [F3], ; for the same bound is immediate. By [F2], . If , nonnegativity of gives .
For such a segment with , . If , and hence , so the step is stationary. Telescoping any finite sequence of admissible full or truncated steps gives total parameter path length at most .
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 with , wholly in the ball. Step 2.1 applies to all previous steps and this truncated step, giving total path length strictly less than . Reaching the boundary from requires length at least , 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.
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 , all matrices with have operator norm at most , , , and all backward feature derivatives have coordinate maxima at most , where . For fully connected networks allow for a fixed , while are fixed independent of . There are depending only on , these fixed bounds and the activation, with the following bounds; write . In the residual case the constants depend also on fixed , and every below is replaced by a fixed independent of . Uniformly on any parameter ball with (fully connected) or (residual),
In the residual last inequality the upper bound is for . The residual first-layer and output-vector derivatives are bounded without a factor . Consequently, if and , the loss Hessian on the ball satisfies
The constants in this last bound can again be fixed independent of for residual networks. For any segment from to contained in the ball, its output Taylor remainder satisfies
with replaced by fixed in the residual case.
Statement references
Facts & Assumptions
Given: The deterministic initialization bounds, the architecture, the parameter ball, and fixed activation constants , , .
The architectures and scalings are defined locally. Fully connected and residual networks with the Du Gram kernel
The chain rule gives the forward variations and backward feature derivatives of the finite network. Reverse accumulation sums gradients over shared parameter uses
Euclidean Cauchy–Schwarz bounds vector inner products, coordinate products, and sums over parameter blocks in the estimates below. Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation
Restricting a scalar 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
A vector-valued map whose derivative norm is bounded by on a segment is -Lipschitz along that segment. If is differentiable with integrable then ; and a bounded derivative makes Lipschitz
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 on real or complex matrices
For , . The -norms for rational , and
Proof
Write and in both architectures, for fully connected , and for residual . A layer has feature Jacobian in the fully connected case, and in residual layers . Within the stated ball, . Products of consecutive , , therefore have norm at most in the fully connected case and in the residual case. All estimates below use only for matrix operator norms; the input is fixed.
Let . Lipschitz continuity gives . For fully connected layers, ; for residual layers it gives , 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 in the first case and in the second. The imposed smallness of keeps all feature norms in the required bounds and justifies the recurrences successively. Also or , since it is bounded by for , and by for .
For a unit parameter direction , let . At the first layer . Subsequent layers obey . The product bounds in step 1.1, the feature bounds in step 2.1, and Cauchy–Schwarz over layer blocks give or . For the residual case use , retaining the separate first-layer term. Thus or , uniformly in . Since and , each output gradient has norm at most or , proving or . This explicitly retains the output and first-layer terms with no bound.
Put . The backward recurrence is , with . Decompose its difference as . The latter is the sum of a weight-difference term bounded by (respectively ) and a gate-difference term bounded by , because . Summing the backward recurrence with the product bounds gives or . Hence or . The crude bound or follows separately from step 1.1 and .
For the last-hidden-layer Gram, put , and or . Direct differentiation gives . Subtract its initialized version in three parts: feature-inner-product change, derivative-gate change with , and output-weight-square change. Step 2.1 bounds the first part by or before multiplying by . For the second part, Cauchy–Schwarz gives ; multiply by , and use step 2.1. For the third part, . Feature norms and bound all remaining factors. The entry error is consequently at most or . The operator norm is at most times this entry bound.
For two unit parameter directions , twice differentiating the output gives the mixed output-weight terms , and the sum over layers of 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 , steps 3.1–3.2, and or bound the first part by or . The terms involving use the backward bound and or ; summing with Cauchy–Schwarz over matrix blocks gives the same bound. The output-weight mixed terms are at most or . This proves the output-Hessian bound, with polynomial factors in absorbed into only for the fully connected case.
The identity gives . Integrating the Jacobian along the straight segment to initialization gives or . Substitution of steps 3.1 and 4.1 proves the claimed loss-Hessian bound, after enlarging to absorb squares of fixed exponential envelopes. For each coordinate the integral Taylor formula is . Step 4.1 bounds its absolute value by or the corresponding residual constant. Summing coordinate squares gives the stated vector Taylor bound.
Simultaneous initialization estimates for the Du architectures
Statement
Use either of the two architectures, its own terminal population kernel , and independent Gaussian initialization. Assume AC, DC and countable choice and the activation and bounded-label hypotheses of the convergence theorems. Fix . Write , and . There are constants and , depending only on the fixed activation bounds, its normalizer, the label bound, and, in the residual case, , with the following property. For the fully connected case put ; if
then with probability at least , simultaneously:
Here is the Jacobian block of the last hidden-layer matrix, and is the derivative of the scalar output with respect to the hidden feature vector , so . In the residual case , the same conclusions hold with replaced everywhere in the conclusions by , provided
For residual depth one use the fully connected case. Constants may be enlarged, and may be increased, simultaneously in these conclusions.
Statement references
Facts & Assumptions
Given: The architecture, its population recursion, the data, Gaussian initialization, width and confidence parameters in the Statement.
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
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
Gaussian vectors have the stated covariance and even moments . Multivariate normal law, including singular covariance Gaussian even moments for Brownian increments
Conditional on hidden features, the Gaussian output norm is at most the square root of their total squared norm plus their operator norm times , with failure probability at most . Gaussian output initialization concentrates conditionally on the hidden features
If have finite second moments, then ; expanding a square and applying this bound gives the triangle estimate used for residual feature variances. Cauchy-Schwarz for random variables
Euclidean Cauchy–Schwarz bounds inner products and the Euclidean norm satisfies the triangle inequality. Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation
A nondegenerate bivariate normal vector has its Gaussian Lebesgue density. A positive-definite multivariate normal law has its Gaussian density
Differentiation may pass through an integral when parameter derivatives have a common integrable majorant. Differentiation under the integral sign
Integration by parts is valid on a compact interval when the derivatives are integrable. If are differentiable on with integrable, then
Dominated convergence passes limits through integrals under a common integrable majorant. Dominated convergence
For nonnegative and , . Markov's inequality for random variables
For square-integrable and , . Chebyshev's inequality for random variables
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 of Cauchy's: for continuous on with and differentiable on there is with
Analyticity gives , and the -Lipschitz property of gives wherever needed in step 2.1. [given]
A -net of the Euclidean unit sphere in exists with at most points. Finite epsilon nets of the Euclidean sphere
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||
For , . The -norms for rational , and
The Frobenius norm is the Euclidean norm of the matrix entry list, ; for any unit , Euclidean Cauchy–Schwarz gives , hence . The Frobenius norm on real or complex matrices Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation
In the residual coupling, the feature vector and fresh Gaussian vector are independent, so . 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.
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
The probability of a finite union of events is at most the sum of their probabilities. The finite union bound
Proof
We give the scalar estimates used below explicitly. If are centered Gaussian variables with variances at most a fixed , then and [F3] imply ; [F5] then gives for integers . The same last bound holds for when is an independent standard normal, because . After centering, each such variable satisfies , enlarging using . Expanding its exponential series, for one has . For independent copies, [F20] factors the exponential moment of their sum. Applying [F11] to with , and then to , to obtain for . All constants are uniform over the covariances with these variance bounds.
The Gaussian expectation maps needed in the population recurrences are Lipschitz functions of the three entries of their bivariate covariance on any region , where is fixed. This applies to for and , and to the residual marginal map . At fixed , changing changes these expectations by at most : use the Lipschitz constants , , and [F5] to bound the resulting products. For fixed and , the covariance matrix is positive definite, so [F7] supplies the density . On any compact interval , its eigenvalues lie between and . Direct differentiation of the Gaussian density gives . Uniform eigenvalue bounds imply that is bounded by an integrable polynomial times , uniformly on this interval; has at most linear growth and is bounded. Thus [F8] justifies differentiating the expectation under the integral. On finite rectangles apply [F9] successively in and ; 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 . Its absolute value is at most for and for , since and . At , the representation , is dominated after multiplication by either pair of activations by , which is integrable; [F10] gives continuity at both endpoints. Finally is Lipschitz in the covariance entries when .
Write and, for residual layers, . Conditional on previous layers, the preactivation pairs in distinct rows are independent centered Gaussians with covariance given by . In the fully connected case, step 1.1 bounds each difference of from its conditional mean by , except with probability , where . Its conditional mean is exactly the fully connected population expectation evaluated at . Since the population diagonal is one, step 2.1 gives, while the error is at most a fixed small constant, , where and . If is a fixed base for this recurrence, then after enlarging its prefactor. Choose the statement constant and enlarge the width constant; its term , with , then makes this error smaller than both the fixed stopping threshold and . Conditional bounds and the finite union bound [F21] over therefore justify the estimates through layers ; 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 .
In the residual case the population feature variances lie in a fixed compact interval independent of . Indeed the vector construction in [F1] gives , whose diagonal is . Positivity of follows from nonconstancy, continuity and the positive Gaussian density. Also with , so iteration gives a fixed upper bound. The empirical conditional mean satisfies exactly the same residual second-moment/mean recursion with replaced by . The nonlinear quadratic part has conditional fluctuation at most by step 1.1. Each cross part is ; its centered value has absolute size at most by [F2], since its Lipschitz constant in the independent row Gaussians is at most . The empirical mean has fluctuation at most by [F2]. Step 2.1 now yields , taking the maximum of the covariance and mean errors. The first layer has by step 1.1 and [F2]. Hence uniformly in . The width assumption and the same stopping argument, now in , justify the conditional estimates and give the feature bound through ; the finite union bound [F21] combines their layerwise failure probabilities. One further conditional step gives , using the same upper-norm recurrence, so .
For either architecture, conditional on layers before the last, the last-hidden-layer Jacobian Gram entries are , where for the fully connected case and for residual . This follows by differentiating the last layer. The summands are independent across . Steps 1.1–2.1, the feature bounds and errors just proved give entry error at most in the fully connected case, for a fixed base depending only on the fixed activation bounds, and at most in the residual case. Choose the statement constant also with ; the displayed logarithmic width term and a sufficiently large fixed multiplier then make both errors at most . A finite union bound [F21] covers the terminal-kernel entry events. If an matrix has entries bounded by , then for any , by [F6], hence . Thus . Both matrices are symmetric, and for every , , where . The Rayleigh-quotient characterization of the smallest eigenvalue now gives .
By [F16], choose a -net of the unit sphere with at most points. The Euclidean norm duality follows from [F6], so . For unit , choose within of them. Then [F6] gives . Taking suprema yields . Each is standard normal, whose tail is at most by [F2]. The finite union bound [F21] gives for all with failure at most when , which follows from either displayed width because in the fully connected case and 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 and with failure at most when , also implied by the widths. These arguments never bound on the whole input space, so no input-dimension condition is used.
Conditional on hidden weights, , where is the product of downstream feature Jacobians. Their operator norms have fixed envelopes in the fully connected case and in the residual case, using for residual layers. Increase the conclusion constants so these envelopes are at most and , respectively. Each coordinate is centered Gaussian of variance at most or . Since , its two-sided tail is at most ; the finite union bound [F21] over at most coordinates gives the asserted coordinate maxima with failure at most . Finally [F4] and give or with failure at most .
By the finite union bound [F21], allocate at most to each fully connected feature-covariance or terminal-feature-norm event in step 3.1; there are at most such events. Allocate at most to each residual covariance, cross-term, mean, or terminal-feature-norm event in step 4.1; there are at most such events. Allocate at most to each terminal-kernel entry event. The weight-matrix and readout norm events in step 5.2 each use ; backward coordinates use and the output norm uses in step 6.1. Even summing both architecture allocations, these budgets total at most . The thresholds are constant multiples of or ; 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.
Wide fully connected networks converge geometrically to zero empirical training loss
Statement
Assume AC, DC and countable choice. Fix positive integers , training data and labels with and for a fixed finite constant . Let the activation satisfy the Lipschitz and smoothness conditions , , and for finite constants, and also be analytic and nonpolynomial. Use the fully connected architecture and population matrix 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 . Write and put and assume . For , if the width and fixed step size satisfy
then with probability at least over initialization, every gradient-descent iterate satisfies
The constants hidden in and may depend on the fixed activation and label bounds but not on ; the step-size constant is small enough that . In particular, geometrically as .
Statement references
Facts & Assumptions
Given: The data, architecture, activation, Gaussian initialization, positive terminal population eigenvalue , width, confidence and step-size hypotheses of the Statement.
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
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
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
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
Put and . Fix an exponential envelope for all initialization constants in [F2] and another for [F3] with those initialization bounds. These bases depend only on the permitted fixed bounds. The source-width form permits , where is chosen sufficiently large in terms of the fixed bases of ; in particular take . All these choices are fixed before selecting the radius. On the event from [F2], of probability at least , satisfies , the initialized last-layer Gram is at least , and all other hypotheses of [F3] hold. Also by the diagonal of [F1].
The displayed width also implies , after fixing sufficiently large. For the logarithmic term this follows without an extra logarithmic width factor: writing and , one has , since for . The constant part of is absorbed by using and . Choose the fixed radius . These inequalities give and . Hence [F3] supplies on the entire ball of radius .
The same width implies , so and . The loss-Hessian bound in [F3], , and yield , where after enlarging its fixed base. Indeed expand in that bound; because and is bounded above, every term is at most an allowed exponential constant times , before the exterior factor . The function is bounded for , so the Hessian is bounded by on the ball, with . The allowed sufficiently small step-size constant ensures .
The radius satisfies . 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 for all . Choosing the step-size constant also ensures , so geometric convergence to zero follows. The initialization event has probability at least , proving the claim. The same application bounds the total parameter path by ; [F3] therefore supplies last-layer Gram stability and the displayed output Taylor-remainder bound on every accepted update segment, for all iterates simultaneously.
Wide residual networks converge geometrically to zero empirical loss with polynomial depth dependence
Statement
Assume AC, DC and countable choice. Fix positive integers , training data and labels with and for a fixed finite constant . Let the activation satisfy the Lipschitz and smoothness conditions , , and for finite constants, and also be analytic and nonpolynomial. Use the residual architecture with fixed and the population matrix 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 . Write and put and assume . For , if
then with probability at least over initialization, every gradient-descent iterate satisfies
The constants hidden in and may depend on the fixed activation, label, and residual-scaling bounds but not on ; the step-size constant is small enough that . In particular, geometrically as .
Statement references
Facts & Assumptions
Given: The data, architecture, activation, Gaussian initialization, positive terminal population eigenvalue , width, confidence and step-size hypotheses of the Statement.
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
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
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
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
For square-integrable real random variables, . Cauchy-Schwarz for random variables
Proof
First take and put , . At the stated width, [F2] supplies with probability at least a simultaneous event with all constants independent of , on which and . By [F1], and . For , put and . The residual recursion gives . By [F1] and [F5], and ; Lipschitz growth and [F5] also give . Hence . Starting from and iterating with gives for every , uniformly in . The terminal diagonal formula and now give ; positive semidefiniteness in [F1] implies . The residual neighborhood bounds [F3] retain separate first-layer and readout gradients without an erroneous factor.
Choose , with fixed sufficiently large for the initial-residual bound. The width implies and . To see the logarithmic absorption, use , and , which follows from exactly as in the fully connected argument. Dividing by only decreases the requirement, while the approved width contains and . The constant part is absorbed using and . Consequently and , so the last-layer Gram is at least throughout this ball.
All constants in the residual Hessian estimate are independent of . Substitution of and the initial-residual bound in [F3] yields , because and . The approved width and imply . Thus , and : explicitly, for all , so the ratio is at most . Also . Therefore on the ball. The permitted sufficiently small fixed step-size constant makes .
Since , [F4] confines every iterate to the ball and gives . The step-size constant can ensure , because . This proves geometric convergence on the probability- event. For 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 ; [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 , whereas the displayed term is . 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
- Simon S. Du, Jason D. Lee, Haochuan Li, Liwei Wang, and Xiyu Zhai, Gradient Descent Finds Global Minima of Deep Neural Networks, Condition 3.1 and equation (2)
- Simon S. Du, Jason D. Lee, Haochuan Li, Liwei Wang, and Xiyu Zhai, Gradient Descent Finds Global Minima of Deep Neural Networks, §§4–6
- Simon S. Du, Jason D. Lee, Haochuan Li, Liwei Wang, and Xiyu Zhai, Gradient Descent Finds Global Minima of Deep Neural Networks, §§3.3, 5, and 6
- Simon S. Du, Jason D. Lee, Haochuan Li, Liwei Wang, and Xiyu Zhai, Gradient Descent Finds Global Minima of Deep Neural Networks, Supplement, Appendix C equation (14) and Appendix E.1 equations (16)–(17)
- Princeton Theory of Deep Learning, Chapter 8, §§8.1–8.5
- Simon S. Du, Jason D. Lee, Haochuan Li, Liwei Wang, and Xiyu Zhai, Gradient Descent Finds Global Minima of Deep Neural Networks, Definition 5.1 and Appendix E
- Simon S. Du, Jason D. Lee, Haochuan Li, Liwei Wang, and Xiyu Zhai, Gradient Descent Finds Global Minima of Deep Neural Networks, §§4 and B.1
- Simon S. Du, Jason D. Lee, Haochuan Li, Liwei Wang, and Xiyu Zhai, Gradient Descent Finds Global Minima of Deep Neural Networks, Appendix B
- Du, Lee, Li, Wang and Zhai, Gradient Descent Finds Global Minima of Deep Neural Networks, Appendices B, C, E and G
- Simon S. Du, Jason D. Lee, Haochuan Li, Liwei Wang, and Xiyu Zhai, Gradient Descent Finds Global Minima of Deep Neural Networks, Theorem 5.1
- Simon S. Du, Jason D. Lee, Haochuan Li, Liwei Wang, and Xiyu Zhai, Gradient Descent Finds Global Minima of Deep Neural Networks, Theorem 6.1
- Simon S. Du, Jason D. Lee, Haochuan Li, Liwei Wang, and Xiyu Zhai, Gradient Descent Finds Global Minima of Deep Neural Networks, Appendix C