Summary of "Tokens-per-Parameter Coverage Is Critical for Robust LLM Scaling Law Extrapolation"

Summary (Overview)

  • Core finding: Scaling-law fits trained on collinear designs (fixed tokens-per-parameter ratio, D=kND = kN) are inherently ill-conditioned when the exponents governing model size (NN) and data size (DD) are nearly equal, making scale coefficients practically unidentifiable.

  • Key theoretical result: The condition number of the Gauss-Newton least-squares problem grows as Θ(ε2)\Theta(\varepsilon^{-2}) where ε\varepsilon is the gap between the NN- and DD-exponents. This is proven for four scaling-law formalisms: Chinchilla, repeated-data, Kaplan, and Droppo-Elibol.

  • Practical consequence: Non-collinear designs outperform collinear ones on held-out data with a 97.3% win rate across four laws, five corpora, and multiple floating-point precision modes (including BF16).

  • Constructive contribution: A closed-form TPP-diversity threshold (Proposition 2) that is necessary and sufficient for well-conditioned estimation, plus an identified reduced model that recovers what single-ray designs can constrain.

  • Root cause: The degeneracy is in Jacobian geometry, not the loss function—any smooth estimation objective whose curvature involves the Jacobian inherits the same ill-conditioning.

Introduction and Theoretical Foundation

Background and Motivation

Scaling laws approximate language model loss as a power-law function of parameter count NN and token count DD. The seminal Chinchilla study (Hoffmann et al., 2022) proposed D20ND \approx 20N for compute-optimal training. However, the authors identify a critical methodological flaw:

"This ratio is a prescription for how to allocate compute when training a single model, not a design for the experimental grid from which scaling laws are derived."

Fixed-TPP grids are widely used for scaling-law estimation despite this distinction, appearing in Cerebras-GPT, OLMo, and other major efforts. When all training runs lie on a single ray D=kND = kN in the (N,D)(N, D) plane, the predictor variables become collinear.

Theoretical Foundation

The standard Chinchilla scaling law is:

L(N,D)=ANα+BDβ+E,(1)L(N, D) = A N^{-\alpha} + B D^{-\beta} + E, \tag{1}

where A,B,α,β,ER>0A, B, \alpha, \beta, E \in \mathbb{R}_{>0} and N,DR>0N, D \in \mathbb{R}_{>0}.

On a single ray (K=1K=1), substituting D=kND = kN collapses the two power-law terms when β=α\beta = \alpha:

Definition 2 (Reduced Chinchilla model under fixed TPP):

L(N;ψ,α,E)=(A+Bkα)Nα+E=ψNα+E,(2)L(N; \psi, \alpha, E) = (A + B k^{-\alpha}) N^{-\alpha} + E = \psi N^{-\alpha} + E, \tag{2}

where ψ:=A+Bkα\psi := A + B k^{-\alpha}.

Definition 3 (Exponent gap): The law-specific exponent gap ε\varepsilon is defined as:

  • Chinchilla/repeated-data: ε:=αβ\varepsilon := |\alpha - \beta|
  • Kaplan: ε:=αDαN\varepsilon := |\alpha_D - \alpha_N|
  • Droppo-Elibol: ε:=γNγD\varepsilon := |\gamma_N - \gamma_D| where γN:=αN/α\gamma_N := \alpha_N/\alpha, γD:=αD/α\gamma_D := \alpha_D/\alpha

Methodology

Gauss-Newton (GN) Framework

The authors fit by nonlinear least squares, minimizing S(θ):=12r(θ)2S(\boldsymbol{\theta}) := \frac{1}{2}\|\mathbf{r}(\boldsymbol{\theta})\|^2 where residuals ri:=LiL^(Ni,Di;θ)r_i := L_i - \hat{L}(N_i, D_i; \boldsymbol{\theta}) stack into r(θ)\mathbf{r}(\boldsymbol{\theta}), with Jacobian JRm×pJ \in \mathbb{R}^{m \times p} given by Jij=ri/θjJ_{ij} = \partial r_i / \partial \theta_j.

Table 1: Jacobian columns under D = kN (showing collinearity structure)

Law (L^\hat{L})Collinear pairj1j_1j2j_2ε\varepsilon
ANα+BDβ+EAN^{-\alpha} + BD^{-\beta} + E(A,B)NαN^{-\alpha}kβNαNεk^{-\beta}N^{-\alpha}N^{\varepsilon}$
ANα+BDβ+EAN'^{-\alpha} + BD'^{-\beta} + E(A,B)NαN'^{-\alpha}kβ(D/N)βNαNεk^{-\beta}(D'/N')^{-\beta}N'^{-\alpha}N^{\varepsilon}$
[(Nc/N)αN/αD+Dc/D]αD[(N_c/N)^{\alpha_N/\alpha_D} + D_c/D]^{\alpha_D}(NcN_c, DcD_c)NαD\propto N^{-\alpha_D}NαDNε\propto N^{-\alpha_D}N^{\varepsilon}$
[(NC/N)αN+(DC/D)αD+E1/α]α[(N_C/N)^{\alpha_N} + (D_C/D)^{\alpha_D} + E^{1/\alpha}]^\alpha(NCN_C, DCD_C)NγN\propto N^{-\gamma_N}NγNNε\propto N^{-\gamma_N}N^{\varepsilon}$

Key Theoretical Results

Proposition 1 (Full-matrix conditioning): If two Jacobian columns satisfy jb=cja+δ\mathbf{j}_b = c\mathbf{j}_a + \boldsymbol{\delta} with δ=O(ε)\|\boldsymbol{\delta}\| = O(\varepsilon), then as ε0\varepsilon \to 0:

κ(JTJ)=Θ(ε2).(3)\kappa(J^T J) = \Theta(\varepsilon^{-2}). \tag{3}

Corollary 1 (Confidence interval inflation): Under i.i.d. Gaussian noise on a single ray D=kND = kN:

CI0.95(A)CI0.95(ψ)=Θ(ε1).(4)\frac{\mathrm{CI}_{0.95}(A)}{\mathrm{CI}_{0.95}(\psi)} = \Theta(\varepsilon^{-1}). \tag{4}

Proposition 2 (TPP diversity threshold): For any κtarget>0\kappa_{\text{target}} > 0 and ordered training ratios k1kKk_1 \leq \cdots \leq k_K:

κ(JTJ)κtarget    VKτK,(5)\kappa(J^T J) \leq \kappa_{\text{target}} \iff V_K \geq \tau_K, \tag{5}

where VKV_K is the second central moment of {kβeff}=1K\{k_\ell^{-\beta_{\text{eff}}}\}_{\ell=1}^K and βeff\beta_{\text{eff}} is the law-specific data-size exponent.

Table 2: Statistical identifiability under single-TPP designs

LawExponentsε\varepsilonε1\varepsilon^{-1}ε2\varepsilon^{-2}κA,B\kappa_{A,B}CI inflation
Chinchillaα0.34,β0.28\alpha \approx 0.34, \beta \approx 0.280.061727810310410^3-10^417×\gtrsim 17\times
KaplanαN0.076,αD0.095\alpha_N \approx 0.076, \alpha_D \approx 0.0950.019532.8×103\sim 2.8\times 10^310410510^4-10^553×\gtrsim 53\times
Droppo-ElibolγNγD\gamma_N \approx \gamma_D0.050.10\sim 0.05-0.1010-20100-40010310410^3-10^41020×\gtrsim 10-20\times

Empirical Validation / Results

Experimental Setup

  • Models: LLaMA-style transformers, 5.04–76.5M parameters across 14 log-spaced configurations
  • Corpora: C4, RedPajama, Wikipedia, Cosmopedia, peS2o
  • Designs: Collinear (CO) with k{1,1.5,1.9,2,2.5,2.7,3,3.3,3.5,4,4.5,5}k \in \{1, 1.5, 1.9, 2, 2.5, 2.7, 3, 3.3, 3.5, 4, 4.5, 5\} vs. Non-collinear (NC) with 14×1214 \times 12 grid spanning a two-dimensional region
  • Fitting: Nonlinear least squares with 100 random restarts and differential-evolution polish, 30 independent optimizer seeds

Aggregate Results

Table 3: R2R^2 and RMSE summary across training designs

SplitDesignR2R^2 Mean95% CIRMSE Mean95% CI
Train (D\mathcal{D})CO0.9848[0.985, 0.985]0.1737[0.173, 0.174]
Train (D\mathcal{D})NC0.9526[0.952, 0.953]0.2023[0.202, 0.203]
Holdout (H\mathcal{H})CO0.8370[0.832, 0.842]0.2373[0.233, 0.241]
Holdout (H\mathcal{H})NC0.9319[0.930, 0.934]0.1561[0.154, 0.158]

Key observation: CO overfits training data (higher R2R^2), but NC generalizes dramatically better on holdout data—improving R2R^2 from 0.837 to 0.932 and cutting RMSE from 0.237 to 0.156.

Table 4: Win rate breakdown (NC vs. CO), overall 97.3% (1,460/1,500)

DatasetWin RateScaling LawWin RateEpochWin Rate
C493.0%Chinchilla97.6%first99.8%
Cosmopedia94.7%Droppo-Elibol98.0%second97.6%
peS2o99.3%Kaplan95.6%final93.8%
RedPajama100.0%Repeated-Data100.0%
Wikipedia99.7%

Precision Invariance

Table 5: BF16 mixed-precision results

SplitDesignR2R^2 MeanRMSE Mean
Train (D\mathcal{D})CO0.98960.1854
Train (D\mathcal{D})NC0.99070.1457
Holdout (H\mathcal{H})CO0.94120.2546
Holdout (H\mathcal{H})NC0.96570.1949

The BF16 win rate (98.7%, CI: [96.6%, 99.5%]) is comparable to full-precision results, confirming the effect is precision-invariant.

Budget-Matched Subset Enumeration

Table 6: Regime A NC-win rate (%) by law and dataset

LawC4RedPajamaWikipediaCosmopediapeS2oDataset Avg.
Chinchilla77.178.688.051.847.168.5 [67.3, 69.8]
Droppo-Elibol62.267.265.945.546.257.4 [56.1, 58.7]
Kaplan76.179.972.950.950.366.0 [64.1, 67.9]

The effect is strongest on web-crawled corpora (77–88% for Chinchilla) and systematically smaller on specialized corpora (Cosmopedia, peS2o).

Theoretical and Practical Implications

Theoretical Implications

  1. Jacobian geometry is the root cause: The degeneracy is not an artifact of the loss function—any smooth estimation objective whose curvature involves the Jacobian inherits the same ill-conditioning.

  2. Explains replication gaps: The Kaplan-Chinchilla discrepancy (Besiroglu et al., 2024; Porian et al., 2025) can be reframed as a predictable consequence of Θ(ε2)\Theta(\varepsilon^{-2}) ill-conditioning: AA's confidence intervals inflate 17×\geq 17\times at Chinchilla exponents and 53×\geq 53\times at Kaplan exponents.

  3. Transferability: Because the degeneracy is purely Jacobian-geometric, the analysis transfers to generative vision, time-series, and scientific ML applications.

Practical Implications

  1. Design diagnostic: Proposition 2's condition (Eq. 5) can be evaluated before training using a literature estimate of βeff\beta_{\text{eff}}, serving as an a-priori design diagnostic.

  2. Reduced model reporting: When the design is collinear, the reduced parameterization (Definition 2) is the appropriate reporting target—along D=kND = kN rays, only the combined effect of NN and DD is well constrained.

  3. IsoFLOP prediction: IsoFLOP curves cut across the training ray and probe the sloppy direction that collinear fitting leaves unconstrained, making them an informative diagnostic.

  4. Scaling 12 TPP ratios doesn't rescue collinearity: The condition number κε2\kappa \propto \varepsilon^{-2} depends on the exponent gap, not the number of rays.

Conclusion

Main Takeaways

  • Single-TPP fits do not support coefficient-level interpretation; claims about the relative importance of model size vs. data volume require training that spans a two-dimensional region of (N,D)(N, D).

  • Collinear designs produce "deceptively confident fits that fail to generalize"—overfitting training data while degrading sharply off the training ray.

  • Non-collinear designs are strictly preferred for scaling-law extrapolation, with a 97.3% win rate across all experimental conditions.

Limitations

  • Theorem 1 gives an ordering, not a quantitative gap
  • Corollary 1 assumes i.i.d. Gaussian residuals (temporal correlation would understate, not invert, the pathology)
  • Small-scale models (5.04–76.5M parameters), single architecture, English-only text
  • Budget-matched validation softens to near-50% on specialized corpora (Cosmopedia, peS2o)

Future Directions

  • Downstream-task and cross-domain extensions that inherit the Jacobian argument
  • Testing the framework on larger-scale models and additional architectures
  • Extending the analysis to data-mixture weights and interaction terms
  • Exploring isoFLOP-grid designs as a natural way to avoid the collinearity pitfall

The paper concludes that the degeneracy is fundamental and design-level: any fixed-TPP experimental design will produce ill-conditioned scaling-law estimates, and researchers should either use non-collinear designs or report the reduced model (Definition 2) when collinearity is unavoidable.

Related papers