REVIEW 2 major objections 2 minor 50 references
Mildly Overparameterized ReLU Networks on Orthogonal Data: Incremental Learning and Implicit Bias
T0 review · 2 major / 2 minor · reviewed 2026-06-29 · grok-4.3
Pith's one-line read Gradient flow on two-layer ReLU networks with orthogonal data converges in the small-initialization limit to a saddle-to-saddle jump process that activates one neuron per saddle.
desk verdict This paper derives a saddle-to-saddle jump process for gradient flow on orthogonal data, giving incremental neuron activation plus a new sqrt(n) scaling for the learned interpolator norm. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The saddle-to-saddle jump process, in which the flow lingers near saddles before a new neuron activates and the trajectory jumps to the next saddle.
What would settle it
Numerical integration of the gradient flow on orthogonal data that fails to exhibit one-neuron-per-saddle activation or produces a final squared norm far from order √n.
Extended reading notes
Core claim
Gradient flow on mildly overparameterized two-layer ReLU networks with orthogonal inputs converges in the small-initialization limit to a saddle-to-saddle jump process in which a new neuron activates at each saddle; this process produces an interpolating solution whose squared ℓ₂-norm scales as √n.
Load-bearing premise
The training data consists of exactly orthogonal vectors.
Editorial extensions
If this is right
- The network reaches exact interpolation of the training set with high probability once its width satisfies m ≳ log(n).
- The interpolating solution obtained by the flow has squared ℓ₂-norm of order √n and is therefore within a constant factor of the minimal-norm interpolator.
- Mild overparameterization suffices for the network to reach an interpolator whose complexity is of the same order as the optimal one.
Reading between the lines
- The same saddle-to-saddle structure may persist under small perturbations of orthogonality, allowing the incremental phenomenon to hold for nearly orthogonal data.
- Analogous limiting jump processes could govern training dynamics in deeper ReLU networks or with other piecewise-linear activations when data satisfy suitable separation conditions.
- The √n-norm scaling indicates that the implicit bias of gradient flow from tiny initialization already selects near-minimal solutions without needing explicit regularization.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes gradient flow dynamics of two-layer ReLU networks trained on orthogonal data from small initialization. It proves that as the initialization scale tends to zero the flow converges to a saddle-to-saddle jump process in which neurons activate incrementally one at a time. This recovers the known high-probability interpolation result for width m ≳ log(n) and derives a new implicit-bias statement that the squared ℓ₂-norm of the learned interpolator scales as √n (within a constant factor of the minimal-norm interpolator).
Significance. If the central claims hold, the work supplies the first rigorous characterization of an incremental learning process for ReLU networks and a novel quantitative implicit-bias result in the mildly overparameterized regime. The orthogonal-data assumption is used explicitly to decouple neuron contributions, enabling an explicit limiting ODE and jump process; this is a clear methodological strength. Recovery of the Dana et al. interpolation threshold together with the new √n scaling provides concrete, falsifiable predictions about both convergence and solution complexity.
major comments (2)
- [§4.2] §4.2, the passage from the finite-scale gradient-flow ODE to the saddle-to-saddle jump process: the argument that the initialization-scale limit produces well-defined jump times at which a new neuron activates while previous neurons remain at equilibrium relies on orthogonality decoupling; the verification that the resulting jump process indeed yields the claimed √n norm scaling (rather than a different constant or scaling) is load-bearing for the implicit-bias claim and requires an explicit check against the minimal-norm interpolator on orthogonal data.
- [§5.1] §5.1, Eq. (27) (or the corresponding statement of the norm result): the derivation that the incremental activation produces ||w||² ~ √n assumes a specific ordering of activation times; if the ordering is random or data-dependent, the constant factor relative to the minimal-norm solution must be shown to remain O(1) uniformly, otherwise the “within a constant factor” claim is not fully supported.
minor comments (2)
- The citation to Dana et al. (2025) appears in the abstract and introduction; the reference list should give the full arXiv identifier and title for consistency with journal style.
- Notation for the network width m and sample size n is introduced early but the dependence on the initialization scale ε is sometimes written as o(1) without an explicit ε→0 statement in every theorem; a uniform convention would improve readability.
Simulated Author's Rebuttal
We thank the referee for their thoughtful review and valuable suggestions. We address each major comment below, providing clarifications and indicating where revisions will be made to strengthen the manuscript.
read point-by-point responses
-
Referee: [§4.2] §4.2, the passage from the finite-scale gradient-flow ODE to the saddle-to-saddle jump process: the argument that the initialization-scale limit produces well-defined jump times at which a new neuron activates while previous neurons remain at equilibrium relies on orthogonality decoupling; the verification that the resulting jump process indeed yields the claimed √n norm scaling (rather than a different constant or scaling) is load-bearing for the implicit-bias claim and requires an explicit check against the minimal-norm interpolator on orthogonal data.
Authors: The orthogonality of the data is explicitly used to decouple the neuron dynamics, enabling the derivation of the limiting jump process as described in §4.2. For the norm scaling, the proof establishes that the incremental process produces an interpolator with ||w||² = Θ(√n). On orthogonal data, the minimal-norm interpolator also satisfies ||w_min||² = Θ(√n). We will include an explicit comparison in the revised version to verify the constant factor is bounded. revision: partial
-
Referee: [§5.1] §5.1, Eq. (27) (or the corresponding statement of the norm result): the derivation that the incremental activation produces ||w||² ~ √n assumes a specific ordering of activation times; if the ordering is random or data-dependent, the constant factor relative to the minimal-norm solution must be shown to remain O(1) uniformly, otherwise the “within a constant factor” claim is not fully supported.
Authors: The activation ordering in the jump process is data-dependent but, under the orthogonal data assumption, the total squared norm is the sum of contributions from each activated neuron, each adding a term of order 1/√n. This sum is invariant to the order of summation, yielding the same Θ(√n) scaling uniformly. We will add a remark in §5.1 clarifying that the constant factor remains O(1) independent of the activation sequence. revision: yes
Circularity Check
No significant circularity; derivation self-contained on orthogonal data
full rationale
The paper derives the saddle-to-saddle jump process from the gradient-flow ODE limit as initialization scale tends to zero, under the explicit orthogonal-data assumption that decouples neuron contributions. From this process it recovers (via external citation to Dana et al.) the interpolation probability m ≳ log n and separately derives the novel ||·||² ~ √n scaling. No equation is shown to equal its own input by construction, no fitted parameter is relabeled as a prediction, and the sole citation is to non-overlapping prior work. The central claims therefore remain independent of the target results.
Assumptions & free parameters
assumptions (2)
- domain assumption Continuous-time gradient flow governs the training dynamics
- domain assumption Training inputs are pairwise orthogonal
Cite this review
Pith. "Pith review of Mildly Overparameterized ReLU Networks on Orthogonal Data: Incremental Learning and Implicit Bias." pith.science (2026). https://pith.science/paper/IB2775ON
@misc{pith2026260527097,
author = {Pith},
title = {Pith review of: Mildly Overparameterized ReLU Networks on Orthogonal Data: Incremental Learning and Implicit Bias},
year = {2026},
howpublished = {\url{https://pith.science/paper/IB2775ON}},
note = {Machine review of arXiv:2605.27097}
}
abstract
The successful training of neural networks hinges on the use of first order optimization methods, yet the theoretical characterization of these methods remains incomplete. This is especially true in settings with mild overparameterization. In this work, we study the gradient flow dynamics of two-layer ReLU networks from small initialization with orthogonal training data. We prove the limiting flow converges to a saddle-to-saddle jump process as the initialization scale tends to zero, revealing an incremental learning phenomenon in which a new neuron activates at each saddle. This analysis recovers the known result of Dana et al. (2025, arXiv:2502.16977) that the network interpolates the training data with high probability as soon as $m \gtrsim \log(n)$, where $m$ is the network width and $n$ is the number of training samples. This incremental process characterization also allows us to derive a novel implicit bias result: the learned interpolator has a squared $\ell_2$-norm scaling as $\sqrt{n}$, which is within a constant factor of the minimal $\ell_2$-norm interpolator. More broadly, our work provides the first rigorous proof of an incremental learning process for ReLU networks, whilst suggesting mildly overparameterized networks can converge to interpolating solutions whose complexity is of the same order as that of the optimal interpolator.
Figures
Figures from the paper (3 more)
Reference graph
Works this paper leans on
-
[1]
Gradient Descent on Two-layer Nets: Margin Maximization and Simplicity Bias
2, 20 Kaifeng Lyu, Zhiyuan Li, Runzhe Wang, and Sanjeev Arora. Gradient Descent on Two-layer Nets: Margin Maximization and Simplicity Bias. In NeurIPS, pages 12978–12991, 2021. 3 Wesley J. Maddox, Gregory W. Benton, and Andrew Gordon Wilson. Rethinking Parameter Counting in Deep Models: Effective Dimensionality Revisited. CoRR, abs/2003.02139, 2020. 4 Har...
-
[2]
On the Convergence of Gradient Descent Training for Two-layer ReLU-networks in the Mean Field Regime
4 Stephan Wojtowytsch. On the Convergence of Gradient Descent Training for Two-layer ReLU-networks in the Mean Field Regime. CoRR, abs/2005.13530, 2020. 2 Blake E. Woodworth, Suriya Gunasekar, Jason D. Lee, Edward Moroshko, Pedro Savarese, Itay Golan, Daniel Soudry, and Nathan Srebro. Kernel and Rich Regimes in Overparametrized Models. In COLT, pages 3635–3673,
-
[3]
Wide Residual Networks
2 Sergey Zagoruyko and Nikos Komodakis. Wide Residual Networks. In BMVC, 2016. 1 14 Appendix Table of contents A Experiments 15 A.1 Experimental details . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 A.2 Additional experiments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16 B Technical backgro...
2016
-
[4]
∀j ∈ N (k) U , ℓ◦ j(κ) + logα(∥ewj(κ)∥) < ξ
-
[5]
∀j ∈ N (k) F , n D(k) j − ∥ewj(κ)∥2 < α δ
-
[6]
∀i ∈ S(k) F , |heθ(κ)(xi) − yi| < α δ
-
[7]
In the above, and throughout the proof, we use ∥·∥ to denote the ℓ2-norm
∀j ∈ [m], ew ⊤ j (κ)D (k) j > 1 − αδ. In the above, and throughout the proof, we use ∥·∥ to denote the ℓ2-norm. The reader should note that δ depends on κ in Hypothesis 1; in particular we have limκ→0 δ = 0. This does not present a problem, since the proof is of convergence on compact subsets of R which exclude (tk)p k=1. Therefore, κ may be viewed as a q...
-
[8]
∀j ∈ Js, s jew ⊤ j (κ)D (0) j > 1 − αδ
Show all 50 references
-
[9]
∀j ∈ Js, ∀t ∈ (0, κ), ℓ◦ j(t) + logα(∥ewj∥) < ξ . 23
-
[10]
The monotonicity of the gradient flow means there exists a maximum rate of growth it is possible for any neuron to achieve
∀i ∈ I−s, ∀j ∈ Js, ew⊤ j (κ)xi ≤ 0. The monotonicity of the gradient flow means there exists a maximum rate of growth it is possible for any neuron to achieve. Therefore, the upper bound on κ appears as a technicality to ensure that even if a neuron were to have grown at this ...
-
[11]
∀i ∈ S(k) F , ∀t ∈ [κ, τ1], |heθ(xi) − yi| < α ε
-
[12]
∀j ∈ [m], ∀t ∈ [κ, τ1], ew ⊤ j D (k) j > 1 − αε/2
-
[13]
∀j ∈ N (k) F , ∀t ∈ [κ, τ1], r n D(k) j − αε/4 < ∥ewj∥ < r n D(k) j + αε/4
-
[14]
tk+1 − 3ε 2 D(k) j⋆ < τ 1 < t k+1 − ε 2 D(k) j⋆
-
[15]
∀j ∈ N (k) U , ∀t ∈ [κ, τ1], ℓ◦ j(t) + logα(∥ewj∥) < ε 2. Proof. Proof of Item 1 Lemma 4. Let δ be such that the conditions of Hypothesis 1 hold. In order to check that |heθ(xi) − yi| < α ε, we will bound the growth in norm of the quantity E(k) = − 1 n X i∈S(k) F heθ(xi) − yi ...
-
[16]
τ2 − τ1 < 3ε D(k) j⋆
-
[17]
∀j ∈ N (k) F , ∀t ∈ [τ1, τ2], r n D(k) j − αε/8 < ∥ewj∥ < r n D(k) j + αε/8
-
[18]
∀i ∈ S(k) F , ∀t ∈ [τ1, τ2], maxi∈S(k) F ew ⊤ j⋆ xi + < α ε/6
-
[19]
∀t ∈ [τ1, τ2], ew ⊤ j⋆ D (k) j⋆ > 1 − αε/12
-
[20]
∀i ∈ S(k) F , ∀t ∈ [τ1, τ2], |heθ(xi) − yi| < α ε/6
-
[21]
∥ewj⋆(τ2)∥ = η. We observe at this point that although Item 4 Lemma 5 would seem to imply a statement on realignment with fitted data, we are required in the proof to obtain Item 3 Lemma 5 first in order to prove Item 4 Lemma 5. This is an interesting result of the fact that a...
-
[22]
τ3 − τ2 < ε′ D(k) j⋆
-
[23]
∀j ∈ N (k) F , ∀t ∈ [τ2, τ3], r n D(k) j − αε′/2 < ∥ewj∥ < r n D(k) j + αε′/2
-
[24]
∀i ∈ S(k) F , ∀t ∈ [τ2, τ3], |heθ(xi) − yi| < α ε′
-
[25]
∀i ∈ S(k) F , ∀t ∈ [τ2, τ3], maxi∈S(k) F ew ⊤ j⋆ xi + < α 5ε′/4
-
[26]
∀t ∈ [τ2, τ3], ew ⊤ j⋆ D (k) j⋆ > 1 − αε′
-
[27]
∥ewj⋆(τ3)∥2 = n D(k) j⋆ − αε′/4. Proof. We work out a useful preliminary bound to begin this section. We upper bound the network output on i ∈ S(k) j⋆ : heθ(xi) (1) < ew ⊤ j⋆ xi∥ewj⋆ ∥2 + mα2cU (2) < yi n D(k) j⋆ + √ 2αε′/2 n D(k) j⋆ − αε′/4 + mα2cU (3) < y i − αε′/3. ...
-
[28]
tk+1 < τ 4 < t k+1 + 4ε D(k) j⋆ + ε′ D(k) j⋆
-
[29]
∀i ∈ S(k+1) F , ∀t ∈ [τ3, τ4], |heθ(xi) − yi| < α ε′/6
-
[30]
∀j ∈ N (k+1) F , ∀t ∈ [τ3, τ4], r n D(k) j − αε′/8 < ∥ewj∥ < r n D(k) j + αε′/8
-
[31]
∀j ∈ N (k+1) U , ew ⊤ j (τ4)D (k+1) j > 1 − αε′/15
-
[32]
∀i ∈ N (k+1) U , ∀t ∈ [τ1, τ4], ℓ◦ j(t) + logα(∥ewj∥) < ζ j(ε, ε′), where ζj(ε, ε′) := ε 2 + D(k) j + q 2L(0) n 9ε 2 D(k) j⋆ + ε′ D(k) j⋆ . Proof. The reader should note it will be convenient to derive the bounds in Item 1 Lemma 7 as conse- quences of analyses for the other po...
-
[33]
∀i ∈ S(k+1) F , ∀t ∈ [τ4, τ5], |heθ(xi) − yi| < α ε′/8. 52
-
[34]
∀j ∈ [m], ∀t ∈ [τ4, τ5], ew ⊤ j D (k+1) j > 1 − αε′/32
-
[35]
∀j ∈ N (k+1) F , ∀t ∈ [τ4, τ5], r n D(k) j − αε′/10 < ∥ewj∥ < r n D(k) j + αε′/10
-
[36]
∀j ∈ N (k+1) U , ∀t ∈ [τ4, τ5], ℓ◦ j(t) + logα(∥ewj∥) < ζ j(ε, ε′) + ε′, where ζj(ε, ε′) is as defined in Lemma 7. Proof. Given any κ′′ ∈ (tk+1, tk+2), Item 1 Lemma 7 guarantees that we may select ε and ε′ sufficiently small that τ4 < κ ′′. Assume that this has been done and t...
-
[37]
∀j ∈ N (p) U , ℓ◦ j(κ) + logα(∥ewj(κ)∥) < ξ
-
[38]
∀j ∈ N (p) F , n D(p) j − ∥ewj(κ)∥2 < α δ
-
[39]
∀i ∈ [n], |heθ(κ)(xi) − yi| < α δ
-
[40]
Then given any κ′ ∈ [κ, ∞),eθα converges uniformly to θ◦ on [κ, κ′]
∀j ∈ [m], ew ⊤ j (κ)D (p) j > 1 − αδ. Then given any κ′ ∈ [κ, ∞),eθα converges uniformly to θ◦ on [κ, κ′]. Proof. We note that Item 3 Lemma 9 and the monotonicity of the gradient flow imply that for allj ∈ [m] and for all t ∈ [κ, ∞), ∥Dj∥ < αδ √n. To estimate the change in the...
-
[41]
∀j ∈ N (k+1) U , ℓ◦ j(κ′′) + logα(∥ewj(κ′′)∥) < ξ ′′
-
[42]
∀j ∈ N (k+1) F , n D(k+1) j − ∥ewj(κ′′)∥2 < α δ′′
-
[43]
∀i ∈ S(k+1) F , |heθ(κ′′)(xi) − yi| < α δ′′ . 57
-
[44]
∀j ∈ [m], ew ⊤ j (κ′′)D (k+1) j > 1 − αδ′′ . The proof consists of two applications of Lemma 8: one choice of ε, ε′ in order to get the existence of δ′′ such that the alignment, norm, and loss conditions all hold at time κ′′; the other choice of ε, ε′ in order to get the arbit...
-
[45]
■ Finally, we are able to provide a proof of Theorem 2
+ ε′| < ξ ′′ for all j ∈ N (k+1) U gives the claim. ■ Finally, we are able to provide a proof of Theorem 2. This proof consists of iteratively concluding uniform convergence on arbitrary closed intervals contained in (tk, tk+1) using Hypothesis 1, and then using Lemma 2 to app...
2017
-
[46]
Since |Isj | ≥ min(n−, n+) almost surely by definition, we then have that a column j is zero with probability at most 1 2 min(n−,n+)
Conditionally on aj, (Ai,j)i∈[n] are independent variables, each of them being either the deterministic variable 0 if i ̸∈ Isj , or a Bernoulli 1 2 i ∈ Isj . Since |Isj | ≥ min(n−, n+) almost surely by definition, we then have that a column j is zero with probability at most 1...
-
[47]
Note that two non-zero columns j and j′ are equal if and only if both sgn(aj) = sgn( aj′) and 1(w⊤ j xi) = 1(w⊤ j′ xi) for all i ∈ Isj . All these variables are independent Bernoulli 1 2 , so that a union bound argument again yields that the probability for two non-zero column...
-
[48]
F ∈ C 1([T1, T2]) and F ′ is Lipschitz continuous
-
[49]
F ≥ 0 and F (T1) = 0
-
[50]
Then for every t ∈ [T1, T2], F ′(t) ≤ q√p + F ′(T1) exp(√p(t − T1))
F ′′(t) ≤ pF (t) + q for almost every t ∈ [T1, T2], where p, q > 0. Then for every t ∈ [T1, T2], F ′(t) ≤ q√p + F ′(T1) exp(√p(t − T1)). Proof. Let G = F + q/p. Then G obeys for almost every t ∈ [T1, T2] the differential inequality G′′ < pG. We compute using this inequality th...
Reviewed June 29, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.