Increasing the batch size during training is a common practice in large language model (LLM) pretraining, yet the theoretical justification behind its success is not well understood. Analyses of stochastic optimization often assume uniformly bounded stochastic gradient variance, yet recent evidence suggests that this assumption fails in many practical nonconvex problems. The Blum--Gladyshev (BG-0) noise model relaxes this assumption by allowing the variance to grow quadratically with the distance from initialization, suggesting that batch size schedulers can help by controlling the variance growth during training. However, this growth can be overly conservative in practice. We empirically investigate variance growth in LLM pretraining and observe that a generalized BG model with a tunable growth exponent provides a tighter description of practical noise behavior. Motivated by this observation, we introduce the generalized BG-a noise model, which interpolates between bounded variance (a=0) and BG-0 noise (a=2). Under L-smoothness, we derive an information-theoretic lower bound with growth-dependent oracle complexity Ω(ε−(4+a)) and establish a matching upper bound in ε-dependence by increasing the batch size as the iterates move away from initialization. Finally, we propose an adaptive batch scheduler that controls variance growth through dynamic batch size adjustments during training. In pretraining OLMo2 models of up to 1B parameters on C4, our scheduler achieves a lower validation loss than both small and large batch training under matched token budgets, while using less than 10% of the iterations of small batch training.
Figures & tables
Figure 1: Empirical stochastic-gradient noise of OLMo2 models on C4 as a function of displacement from initialization. For each radius r , we evaluate x=x0+ru along ten independently sampled, normalized Gaussian directions and plot the corresponding minibatch variance estimates. The dashed red curve in each panel is the fitted BG- a model, ( B2∥x−x0∥a+G2 ). The fitted parameters are: (a) OLMo2-100M: B=0.379,G=5.41,a=0.98 ; (b) OLMo2-600M: B=0.132,G=8.64,a=1.38 ; and (c) OLMo2-1B: B=0.0916,G=9.77,a=1.49 . Across all three model sizes, the fitted curve captures both the initial noise floor governed by ( G2 ) and the distance-dependent growth governed by B2∥x−x0∥a , with the fitted exponent a increasing with model size. The corresponding results for the vision models are provided in Appendix E.2.1 .
Figure 2: Estimated BG- a noise-model parameters across OLMo2 model sizes. We observe that as model size increases, (a) B tends to decrease, while (b) G and (c) a tend to increase.
Figure 3: Results for OLMo2-100M (top), OLMo2-600M (middle) and OLMo2-1B (bottom) pretrained on C4. (a) Validation loss for a fixed small batch ( b=32 for 100M and 600M, b=64 for 1B), a fixed large batch ( b=512 for 100M and 600M, b=2048 for 1B), and BG- a , with batch size across iterations shown in the inset. (b) Modeled variance ratio for BG- a , with parameters B=0.379,G=5.41,a=0.98 (100M), B=0.132,G=8.64,a=1.38 (600M) and B=0.0916,G=9.77,a=1.49 (1B).
Model
Batch Size
Validation Loss
Iterations ( T )
Iterations (%)
OLMo2-100M
Small ( b=32 )
3.4609±0.0037
30,518
100.00%
Large ( b=512 )
3.5139±0.0071
1,908
6.25%
BG- a
3.4347±0.0085
2,673
8.76%
OLMo2-600M
Small ( b=32 )
2.9188±0.0062
183,106
100.00%
Large ( b=512 )
2.8710±0.0079
11,444
6.25%
BG- a
2.8482±0.0098
14,824
8.10%
Table 1: Validation performance and optimization cost of OLMo2 models pretrained on C4. Validation loss is reported as mean ± std across 3 seeds, with iteration percentages relative to small batch runs.
Appendix figures & tables10 assets
Supplementary material from the paper’s appendix.
Appendix
Model
Batch Size
HellaSwag
CommonsenseQA
OLMo2-100M
Small ( b=32 )
28.98±0.08
29.03±0.06
Large ( b=512 )
27.97±0.25
27.27±0.01
BG- a
29.15±0.30
29.07±0.46
OLMo2-600M
Small ( b=32 )
40.43±0.11
32.43±0.04
Large ( b=512 )
43.05±0.15
32.92±0.11
BG- a
43.11±0.32
33.33±0.34
Appendix
Table 2: Downstream performance of OLMo2 models pretrained on C4. Entries are accuracy (%), reported as mean ± std across 3 seeds.
Figure 4: Empirical stochastic-gradient noise of ResNet50 as a function of displacement from initialization. For each radius r , we evaluate x=x0+ru along ten independently sampled, isotropic Gaussian directions and plot the corresponding minibatch variance estimates. The dashed red curve in each panel is the fitted BG- a model, ( B2∥x−x0∥a+G2 ). The fitted parameters are: (a) CIFAR-10: B=24.9 , G=305 , a=1.25 ; (b) CIFAR-100: B=29 , G=337 , a=1.19 ; and (c) ImageNet-1K: B=7.01 , G=56.3 , a=0.96 .
Figure 5: Results for ResNet50 trained on CIFAR-10. (a) Training loss for fixed batches b=64 , b=512 , and BG- a , with batch size across iterations shown in the inset. (b) Modeled variance ratio for BG- a , with parameters a=1.25,B=24.9,G=305
Figure 6: Results for ResNet50 trained on CIFAR-100. (a) Training loss for fixed batches b=64 , b=512 , and BG- a , with batch size across iterations shown in the inset. (b) Modeled variance ratio for BG- a , with parameters a=1.19,B=29,G=337
Figure 7: Results for ResNet50 trained on ImageNet-1K. (a) Training loss for fixed batches b=128 , b=1024 , and BG- a , with batch size across iterations shown in the inset. (b) Modeled variance ratio for BG- a , with parameters a=0.96,B=7.01,G=56.3
Dataset
Batch Size
Train Loss
Test Loss
Iterations ( T )
Iterations (%)
CIFAR-10
Small ( b=64 )
0.000262±0.000020
0.2253±0.0086
156,200
100.00
Large ( b=512 )
0.000204±0.000021
0.2508±0.0028
19,400
12.42
BG- a
0.000200±0.000009
0.2412±0.0028
23,800
15.24
CIFAR-100
Small ( b=64 )
0.002215±0.000128
0.9522±0.0052
156,200
100.00
Large ( b=512 )
0.001496±0.000009
0.9951±0.0244
19,400
12.42
BG- a
0.001746±0.000193
0.9609±0.0160
23,440
15.00
Appendix
Table 3: Training and test performance of the ResNet50 model trained on CIFAR-10, CIFAR-100 and ImageNet-1K. Entries are loss, reported as mean ± std across 3 seeds. Iterations are averaged across seeds and given as total iterations and percentage relative to the small-batch baseline within each dataset.
Dataset
Method
Epochs
Batch size
Base LR
Time (min)
CIFAR-10
Small batch
200
64
0.05
338
Large batch
200
512
0.4
206
BG- a
200
64 ∗
0.05 ∗
240
CIFAR-100
Small batch
200
64
0.05
337
Large batch
200
512
0.4
206
BG- a
200
64 ∗
0.05 ∗
240
Appendix
Table 4: Vision training configurations and mean wall-clock times, using four NVIDIA A100 GPUs. For BG- a , the batch size and base learning rate refer to the initial batch size and its associated base learning rate. ( ∗ ) means that the given values are the initial starting values.
Model
Method
Tokens
Batch size
Base LR
Time
OLMo2-100M
Small batch
2B
32
0.004
28m 37s
Large batch
2B
512
0.016
16m 13s
BG- a
2B
32 ∗
0.004 ∗
17m 20s
OLMo2-600M
Small batch
12B
32
0.001
12h 58m
Large batch
12B
512
0.004
9h 03m
BG- a
12B
32 ∗
0.001 ∗
9h 17m
Appendix
Table 5: Language model pretraining configurations on C4 and mean wall-clock times, using four NVIDIA H100 GPUs. For BG- a , the batch size and base learning rate refer to the initial batch size and its associated base learning rate. ( ∗ ) means that the given values are the initial starting values.
Model
Dataset
H100 GPUs
Time
ResNet50
CIFAR-10
1
2m 44s
ResNet50
CIFAR-100
1
2m 38s
ResNet50
ImageNet-1K
1
11m 49s
OLMo2-100M
C4
1
50m 08s
OLMo2-600M
C4
4
2h 34m 28s
OLMo2-1B
C4
4
2h 56m 39s
Appendix
Table 6: Computational resources and wall-clock times for the BG- a noise growth ablation studies.
OLMo2-100M
OLMo2-600M
OLMo2-1B
Parameters
100M
600M
1B
Layers
12
16
16
Hidden dimension
512
1344
2048
Attention heads
8
16
16
Head dimension
64
84
128
FFN hidden dimension
2048
5376
8192
Appendix
Table 7: Architectural details of the OLMo2 models used in our experiments. All models use a decoder-only Transformer architecture with SwiGLU feed-forward layers, RMSNorm, QK-normalized attention, and rotary positional embeddings.