Pith. sign in

REVIEW 2 major objections 17 references

Balancing Learning Rates Across Layers: Exact Two-Step Dynamics and Optimal Scaling in Linear Neural Networks

T0 review · 2 major / 0 minor · reviewed 2026-06-28 · grok-4.3

Pith's one-line read Test loss in linear neural networks is minimized by unequal learning rates across layers in the first gradient step, with equal rates optimal afterward.

desk verdict Exact closed forms for the first two GD steps are the real contribution; the optimality claim for unequal rates rests on an approximation whose effect on the argmin is not bounded. read the letter →

arxiv 2606.00340 v1 pith:OZ72RCI5 submitted 2026-05-29 cs.LG

classification cs.LG
keywords learningratescalinglinearneuralnetworksgradientdescentdynamicslayer-wiseratessurrogatelossearlytraining
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper derives exact closed-form expressions for the gradients and test loss after one and two steps of gradient descent in two-layer and three-layer linear networks trained on linear targets. It shows that a gradient approximation permits different per-layer learning rates in the first step to achieve lower loss, while equal rates become optimal in the second step. This approximation produces a surrogate loss whose error is provably small and tight. The result matters because it isolates a distinct early-training regime where balancing layer-wise rates improves outcomes, as the experiments confirm.

What carries the argument

Exact closed-form expressions for the gradients and test loss after one and two steps of gradient descent, which support the characterization of learning-rate scaling under the approximation.

What would settle it

Training a two-layer linear network on a linear target for two gradient steps and measuring whether the test loss after the first step is lower with unequal per-layer rates than with equal rates would confirm or refute the central claim.

Watch

Extended reading notes

Core claim

In two- and three-layer linear neural networks trained to learn linear target functions, the exact closed-form expressions for gradients and test loss after one and two steps of gradient descent show that optimal learning rates are unequal across layers at the initial step but equal in subsequent steps. Performing updates with the gradient approximation yields a tractable surrogate loss with a tight, small approximation error, enabling analysis of layer-wise scaling.

Load-bearing premise

The closed-form derivations and optimality claims hold only for linear networks with linear target functions and are restricted to the first two gradient steps using the gradient approximation.

Editorial extensions

If this is right

  • Unequal learning rates across layers reduce test loss more than equal rates do in the first step.
  • Equal learning rates become optimal from the second step onward.
  • The surrogate loss approximation has provably small error and supports further theoretical analysis of layer-wise rates.
  • Numerical experiments on two- and three-layer networks confirm the early-training regime where layer balancing matters.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The two-step exact dynamics could be checked for repetition in three or more steps to see if the unequal-to-equal transition pattern persists.
  • The initial unequal-rate regime might be tested as a practical heuristic in models that are approximately linear near initialization.
  • The surrogate loss construction could be applied to study scaling in wider linear networks without changing the core approximation.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 0 minor

Summary. The paper derives exact closed-form expressions for gradients and test loss after one and two steps of gradient descent in two- and three-layer linear neural networks learning linear targets. It characterizes optimal layer-wise learning-rate scaling under a gradient approximation for the first two steps, proves that updates with this approximation yield a tractable surrogate loss with tight small approximation error, and shows that unequal learning rates minimize test loss at the initial step while equal rates become optimal subsequently. Numerical experiments validate the theory; code is provided.

Significance. If the results hold, this supplies precise early-training dynamics for linear networks and highlights the value of unequal layer-wise rates at initialization. The reproducible code and numerical validation are explicit strengths that support the claims.

major comments (2)
  1. [Abstract] Abstract: the claim that the surrogate loss has a 'tight, small approximation error' underpins the optimality conclusion for unequal rates at step 1, yet no Lipschitz or sensitivity bound is supplied showing that the error does not shift the argmin over the learning-rate vector.
  2. [Abstract] Abstract: the exact closed-forms and optimality statements are derived only for the first two steps under linear networks and linear targets; the paper provides no argument that the identified early-training regime (unequal then equal rates) survives beyond these restrictions.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The comments highlight two areas where the abstract claims can be strengthened with additional rigor and clearer scoping. We address each point below and will revise the manuscript accordingly.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the claim that the surrogate loss has a 'tight, small approximation error' underpins the optimality conclusion for unequal rates at step 1, yet no Lipschitz or sensitivity bound is supplied showing that the error does not shift the argmin over the learning-rate vector.

    Authors: We agree that a formal sensitivity analysis would strengthen the link between the approximation error and the optimality of unequal rates. The current manuscript demonstrates small error numerically and shows that the surrogate preserves the qualitative ordering of test loss, but does not supply an explicit Lipschitz or perturbation bound on the argmin. In the revision we will add a short sensitivity lemma bounding the change in the optimal learning-rate vector as a function of the approximation error, using the fact that the surrogate loss is quadratic in the rates under the linear-network setting. revision: yes

  2. Referee: [Abstract] Abstract: the exact closed-forms and optimality statements are derived only for the first two steps under linear networks and linear targets; the paper provides no argument that the identified early-training regime (unequal then equal rates) survives beyond these restrictions.

    Authors: The paper deliberately restricts attention to the first two gradient steps in linear networks with linear targets precisely because this regime admits exact closed forms. We do not claim that the unequal-then-equal pattern extends to deeper networks, nonlinear activations, or later training phases; the contribution is the exact characterization and the resulting insight that layer-wise rates should be balanced after the initial step. In the revision we will modify the abstract and introduction to state the scope more explicitly and add a brief paragraph on the limitations and possible extensions. revision: partial

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity; derivations are self-contained exact closed forms plus independent error bound

full rationale

The paper derives exact closed-form expressions for gradients and test loss after one and two gradient steps on linear networks, then introduces a gradient approximation to obtain a surrogate loss whose approximation error is separately bounded as tight and small. No step reduces a claimed prediction or optimality result to a fitted parameter by construction, nor does any load-bearing premise rest on a self-citation chain, imported uniqueness theorem, or ansatz smuggled from prior work. The central early-training regime claim follows directly from the closed forms and the stated error bound without circular reduction to inputs.

Assumptions & free parameters 0 free parameters · 2 assumptions · 0 invented entities

Analysis rests on standard gradient descent update rules and linearity of both model and target; no free parameters or new entities are introduced in the abstract.

assumptions (2)
  • standard math Gradient descent follows the standard parameter update rule using the gradient of the loss.
    Invoked implicitly when deriving closed-form expressions after one and two steps.
  • domain assumption Both the neural network and the target function are linear.
    Stated in the problem setup; required for closed-form tractability.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Balancing Learning Rates Across Layers: Exact Two-Step Dynamics and Optimal Scaling in Linear Neural Networks." pith.science (2026). https://pith.science/paper/OZ72RCI5

@misc{pith2026260600340,
  author       = {Pith},
  title        = {Pith review of: Balancing Learning Rates Across Layers: Exact Two-Step Dynamics and Optimal Scaling in Linear Neural Networks},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OZ72RCI5}},
  note         = {Machine review of arXiv:2606.00340}
}
read the original abstract

We study optimal learning-rate selection in two-layer and three-layer linear neural networks trained to learn linear target functions. In particular, we derive the exact closed-form expressions for the gradients and test loss after one and two steps of gradient descent, enabling a precise characterization of early training dynamics. We characterize how learning rates should scale under the gradient approximation in the first two steps, and prove that performing updates with this approximation yields a tractable surrogate loss with a tight, small approximation error. This formulation enables the theoretical analysis of layer-wise learning rates and reveals a distinct early-training regime: test loss can be minimized by unequal learning rates at the initial step, while equal learning rates become optimal in subsequent steps. Our numerical experiments validate the theory and demonstrate the importance of balancing layer-wise learning rates early during training. The code is available at: https://github.com/TDCSZ327/Layer-Balancing.

Figures

Figures reproduced from arXiv: 2606.00340 by the authors.

Figure 1
Figure 1. 2-layer NN under orthogonal initialization. Here we set η1 + η2 ≤ O(h 3 2 ) and h = 1000. We observe that the theoretical losses closely track the empirical test losses measured after either one or two updates. Moreover, a clear qualitative shift emerges: after a single update, symmetric learning rates across layers are suboptimal, whereas after two updates they become locally optimal. We discuss this in Section 6. … view at source ↗
Figure 2
Figure 2. 3-layer NN under orthogonal initialization. Here we set η1 +η2 ≤ O(h 2 3 ) and h = 1000. We observe conclusions and results similar to those in [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. We visualize the ESDs of gradient matrices and weight matrices {A0 l } 2 l=1, {B0 l } 2 l=1, {G0 l } 2 l=1, {Af1 l } 2 l=1, {Bf1 l } 2 l=1, {Gf1 l } 2 l=1, {Wg1 l } 2 l=1 and {Wg2 l } 2 l=1 and norm gap with η1 = η2 = h 3 2 and h = 1000. Orthogonal initialization. We also ran the same set of experiments under Gaussian initialization In [PITH_FULL_IMAGE:figures/full_fig_p070_3.png] view at source ↗
Figures from the paper (14 more)
Figure 4
Figure 4. Figure 4: Test loss of a 2-layer NN under orthogonal initialization with width h = 1000. Here we set η1 + η2 ≤ O(h 3 2 ). Our theory accurately predicts the test loss after 1-step and 2-step gradient descent updates with varying learning rates. In particular, we highlight the ro…
Figure 5
Figure 5. Figure 5: More-steps-empirical-loss for 2-layer NN under Orthogonal initialization. Here we set η1 + η2 ≤ O(h 3 2 ) and h = 1000. 0 10000 η1 0.965 0.970 0.975 0.980 Test Loss 10000 12000 14000 16000 18000 20000 η1 + η 2 (a) Step=1 0 10000 η1 0.90 0.92 0.94 0.96 Test Loss 10000 1…
Figure 6
Figure 6. Figure 6: More-steps-empirical-loss for 2-layer NN under Gaussian initialization. Here we set η1 + η2 ≤ O(h 3 2 ) and h = 1000. 0 20 η1 0.945 0.950 0.955 0.960 Test Loss 30 32 34 36 38 40 η1 + η 2 (a) Step=1 0 20 η1 0.82 0.84 0.86 0.88 0.90 Test Loss 30 32 34 36 38 40 η1 + η 2 (…
Figure 7
Figure 7. Figure 7: More-steps-empirical-loss for 3-layer NN under Orthogonal initialization. Here we set η1 + η2 ≤ O(h 2 3 ) and h = 1000. 0 20 η1 0.94 0.95 0.96 Test Loss 30 32 34 36 38 40 η1 + η 2 (a) Step=1 0 20 η1 0.80 0.85 0.90 Test Loss 30 32 34 36 38 40 η1 + η 2 (b) Step=2 0 20 η1…
Figure 8
Figure 8. Figure 8: More-steps-empirical-loss for 3-layer NN under Gaussian initialization. Here we set η1 + η2 ≤ O(h 2 3 ) and h = 1000. 72 [PITH_FULL_IMAGE:figures/full_fig_p072_8.png]
Figure 9
Figure 9. Figure 9: More-steps-empirical-loss for 3-layer NN under Orthogonal initialization up to 512 steps. Here we set η1 + η2 ≤ O(h 2 3 ) and h = 1000. (a) 2-layer NN Step=1 (b) 2-layer NN Step=2 (c) 3-layer NN Step=1 (d) 3-layer NN Step=2 [PITH_FULL_IMAGE:figures/full_fig_p073_9.png]
Figure 10
Figure 10. Figure 10: (a)(b)2-layer NN with label noise ξ ∈ N (0, ρ) under orthogonal initialization. Here we set η1+η2 ≤ O(h 3 2 ) with h = 1000 and ρ = 0.001. (c)(d)3-layer NN with label noise ξ ∈ N (0, ρ) under orthogonal initialization. Here we set η1 + η2 ≤ O(h 2 3 ) with h = 1000 and…
Figure 11
Figure 11. Figure 11: 2-layer NN under orthogonal initialization. Here we set η1 + η2 ≤ O(h 3 2 ) and h = 5000. 73 [PITH_FULL_IMAGE:figures/full_fig_p073_11.png]
Figure 12
Figure 12. Figure 12: 3-layer NN under orthogonal initialization. Here we set η1 + η2 ≤ O(h 2 3 ) and h = 5000. 0 1000 η1 0.65 0.70 0.75 0.80 Test Loss 1000 1200 1400 1600 1800 2000 η1 + η 2 (a) 1-step (theory) 0 1000 η1 0.65 0.70 0.75 0.80 Test Loss 1000 1200 1400 1600 1800 2000 η1 + η 2 …
Figure 13
Figure 13. Figure 13: 2-layer NN under orthogonal initialization. Here we set η1 + η2 ≤ O(h 3 2 ) and h = 100. We can see since h does not satisfy the condition on h in Corollary 5.4, the balanced learning-rate allocation is not locally optimal. (a) 4-Layer 1-step (b) 4-Layer 2-step (c) 8-…
Figure 14
Figure 14. Figure 14: 4-layer and 8-layer NN under orthogonal initialization for 1 and 2-step updates. For 4-NN, we set η1 + η2 + η3 = C ≤ O(h 2 3 ) with h = 1000 and we set η2 = η3 = C−η1 2 . For 8-NN, we set η1 + η2 + η3 + η4 + η5 + η6 + η7 = C ≤ O(h 2 3 ) with h = 1000 and we set η2 = η…
Figure 15
Figure 15. Figure 15: 3-NN nonlinear under orthogonal initialization for 1 and 8-step updates. Here we consider student model is f(xi) = √1 h σ(σ(x ⊤ i W1)W2)a, and the teacher model is yi = σ(β ∗⊤xi), with σ being the ReLU activation. 75 [PITH_FULL_IMAGE:figures/full_fig_p075_15.png]
Figure 16
Figure 16. Figure 16: (a) 3-layer NN under orthogonal initialization for special 2-step update. Here we set η1 + η2 = C ≤ O(h 2 3 ) with h = 1000. We set the first step update step to use an asymmetric learning-rate allocation: the first layer is updated with learning rate C, while the sec…
Figure 17
Figure 17. Figure 17: Insights for designing layer-wise lr scheduler. Here we condier W1 ∈ R 60×100 ,W2 ∈ R 100×60 and M ∈ R 60×60, with ∥W1∥F = 1 and ∥W2∥F = 6 at initialization. We use 100 training samples and 20 test samples, base lr = 0.0001, adopt the MSE loss, and train the model usi…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

17 extracted references · 1 canonical work pages

  1. [1]

    SGDR: Stochastic Gradient Descent with Warm Restarts

    PMLR, 2015. Loshchilov, I. and Hutter, F. Sgdr: Stochastic gra- dient descent with warm restarts.arXiv preprint arXiv:1608.03983, 2016. Lu, H., Zhou, Y ., Liu, S., Wang, Z., Mahoney, M. W., and Yang, Y . Alphapruning: Using heavy-tailed self regular- ization theory for improved layer-wise pruning of large language models.Advances in Neural Information Pro...

  2. [2]

    = 1√ h (24) A0 2 = A0 2 F = q tr(A0 2 ⊤A0

  3. [3]

    We have B0 1 = B0 1 F = q tr(B 0 1 ⊤B0

    = 1√ h .(25) We also haveB 0 1 = 1 h W 0 1 W 0 2 aa⊤W 0⊤ 2 ,B 0 2 = 1 h W 0⊤ 1 W 0 1 W 0 2 aa⊤ are both rank-1 matrices. We have B0 1 = B0 1 F = q tr(B 0 1 ⊤B0

  4. [4]

    = 1 h (26) B0 2 = B0 2 F = q tr(B 0 2 ⊤B0

  5. [5]

    Thus, we get that G0 1 −A 0 1 ≤ 1√ h−1 G0 1 , G0 2 −A 0 2 ≤ 1√ h−1 G0 2

    = 1 h .(27) SinceG 0 1 =B 0 1 −A 0 1,G 0 2 =B 0 2 −A 0 2, we obtain that G0 1 −A 0 1 ≤ 1√ h A0 1 ≤ 1√ h ( G0 1 + G0 1 −A 0 1 ) G0 2 −A 0 2 ≤ 1√ h A0 2 ≤ 1√ h ( G0 1 + G0 2 −A 0 2 ). Thus, we get that G0 1 −A 0 1 ≤ 1√ h−1 G0 1 , G0 2 −A 0 2 ≤ 1√ h−1 G0 2 . (28) Based on this, we can get √ h G0 1 = Θh,P(1), √ h G0 1 F = Θh,P(1), √ h G0 2 = Θh,P(1), √ h G0 2...

  6. [6]

    = 1 h (32) A0 2 = 1 h √ h , A0 2 F = q tr(A0 2 ⊤A0

  7. [7]

    = 1 h .(33) We also haveB 0 1 = 1 h2 W 0 1 W 0 2 W 0⊤ 2 ,B 0 2 = 1 h W 0⊤ 1 W 0 1 W 0 2 , so we can get that B0 1 = 1 h2 , B0 1 F = q tr(B 0 1 ⊤B0

  8. [8]

    = 1 h √ h (34) B0 2 = 1 h2 , B0 2 F = q tr(B 0 2 ⊤B0

Show all 17 references
  1. [9]

    Thus, we get that G0 1 −A 0 1 ≤ 1√ h−1 G0 1 , G0 2 −A 0 2 ≤ 1√ h−1 G0 2

    = 1 h √ h .(35) SinceG 0 1 =B 0 1 −A 0 1,G 0 2 =B 0 2 −A 0 2, we obtain that G0 1 −A 0 1 ≤ 1√ h A0 1 ≤ 1√ h ( G0 1 + G0 1 −A 0 1 ) G0 2 −A 0 2 ≤ 1√ h A0 2 ≤ 1√ h ( G0 1 + G0 2 −A 0 2 ). Thus, we get that G0 1 −A 0 1 ≤ 1√ h−1 G0 1 , G0 2 −A 0 2 ≤ 1√ h−1 G0 2 . (36) Based on thi...

  2. [10]

    1 h W 1 1 W 1 2 −M ⊤ ˜x⊤ 0 ˜x0 1 h W 1 1 W 1 2 −M #! =tr EW 0 1 ,W 0 2 ,ξ, ˜x0,X

    Also consider each row row (or column) of Q is a random vector uniformly distributed on the unit sphere in Rh. Hence by the definitions of orthogonal group, we have hX a=1 Q2 ia = 1⇒E Q2 ia = 1 h , furthermore, if we consider flipping the sign of one row or one column like lef...

  3. [11]

    Then, for any α in this range, the point η1 =η 2 =h α is not a local minimum of the lossLtwo-layer(W 1 1 ,W 1 2 ). Moreover, for1< α≤ 3 2, if h >max{h ∗,256} , then η1 =η 2 =h α is a local minimum of the lossL two-layer(W 2 1 ,W 2 2 ), whereh ∗ is the root of the following equ...

  4. [12]

    Given the fixed 1< α < 8 5, we will give how large h is to ensure that η1 =η 2 =h α will are local minima, Case 1.Ifα= 3 2, we need 8h 1 2 −32−32−64−o(1)>0, which meansh >256 +o(1). 29 Balancing Learning Rates Across Layers: Exact Two-Step Dynamics and Optimal Scaling in Linea...

  5. [13]

    +η 1η2A0 1(A1 2 −fA1 2)−η 1η2B0 1A1 2 +η 1η2A1 1A1 2 −η 1η2B1 1A1 2 −η 1η2fA1 1 fA1 2 +η 1η2fB1 1 fA1 2 +η 2W 0 1 (B1 2 − fB1

  6. [14]

    +η 1η2A0 1(B1 2 − fB1

  7. [15]

    1√ h W 1 1 W 1 2 a−β ∗ ⊤ ˜x⊤ 0 ˜x0 1√ h W 1 1 W 1 2 a−β ∗ # =tr EW 0 1 ,W 0 2 ,a,ξ, ˜x0,X

    +η 1η2B0 1B1 2 −η 1η2A1 1B1 2 +η 1η2B1 1B1 2 +η 1η2fA1 1 fB1 2 −η 1η2fB1 1 fB1 2 −η 2W 0 1 B0 2 −η 1η2A0 1B0 2 +η 1η2B0 1B0 2 −η 1η2A1 1B0 2 +η 1η2B1 1B0 2 34 Balancing Learning Rates Across Layers: Exact Two-Step Dynamics and Optimal Scaling in Linear Neural Networks We know ...

  8. [16]

    Given the fixed 0< α < 2 3, we will give how large h is to ensure that η1 =η 2 =h α will be the local minima, We need 32h3α−2 + 33hα−1 + 74hα−2 + 2h−α + 10h−α−1 + 36h3α−3 + 4h5α−4 −8<0. □ E. Gaussian Initialization In this section, to obtain more general and practical results,...

  9. [17]

    1 h W 1 1 W 1 2 −M ⊤ ˜x⊤ 0 ˜x0 1 h W 1 1 W 1 2 −M #! =tr EW 0 1 ,W 0 2 ,ξ, ˜x0,X

    Then, for any α in this range, the point η1 =η 2 =h α is not a local minimum of the lossL two-layer(W 1 1 ′ ,W 1 2 ′ ). We do simulations in Figure 6 in Appendix F to support Corollary E.5. E.2.2. THREE-LAYERNEURALNETWORKS Given test data ˜x0 ∼ N(0,I d), we consider the test l...

Pith tools

Reviewed June 28, 2026 · model on record in the stance chip above.