Organizations: Department of Computer Science and Engineering Indraprastha Institute of Information Technology, Delhi, India. · Department of Mathematics Indraprastha Institute of Information Technology, Delhi, India.
Adversarial sharpness penalty per iteration of a base optimizer, known as sharpness-aware minimization (SAM), is an accepted technique to improve generalization in deep learning. Recently, Muon has emerged as a dominant optimization algorithm for training large language models (LLMs). It achieves lower sample complexity for heavy-tailed gradient updates - a phenomenon common in LLM training. However, sharpness penalties with Muon updates are yet largely unexplored in the literature. Moreover, the behavior of SAM under heavy-tailed noise is unknown both theoretically and empirically. A potential reason for this gap is the additional cost of extra gradient computation per iteration for penalty perturbations in SAM, which becomes prohibitive for LLMs. This paper introduces Sharpness-aware Orthogonalized Momentum Algorithm (SOMA) to combine a computationally efficient and theoretically grounded perturbation scheme with the Muon algorithm. The key idea is to use stale (purely) stochastic gradient noise for adversarial perturbations. Theoretically, we show that SOMA's sample complexity is invariant to the perturbation direction. We provide explorations and insights into SOMA's robustness to heavy-tailed gradient noise. We extensively explore SOMA's efficacy in (a) mitigating quality degradation due to post-training quantization, and (b) catastrophic forgetting during fine-tuning of LLMs.
Figures & tables
Figure 1: Performance under anisotropic heavy-tailed noise (Condition no., κε=10 ). Top: validation accuracy. Bottom: training loss (log scale). Red , Blue and Green lines indicate AdamW -based, SGD -based and Muon -based methods, respectively. Lines show the median over seeds, and shaded bands show the interquartile range.
Figure 2: Training curves of the fifteen optimizers (mean over 3 seeds) on FineWeb-Edu NanoGPT, FineWeb-Edu Pythia and the CIFAR-10 Airbench CNN: training loss, and validation perplexity (language models) or validation accuracy (CNN), against the optimizer step (100 steps per epoch). Red , Blue and Green lines indicate AdamW -based, SGD -based and Muon -based methods, respectively.
FineWeb-Edu NanoGPT
FineWeb-Edu Pythia
Optimizer
FP
INT4
Δ%
FP
INT4
Δ%
Muon
22.96 ± 0.01
23.87 ± 0.01
+3.98 ± 0.06
33.60 ± 0.01
37.28 ± 0.10
+10.96 ± 0.27
Ssm
23.09 ± 0.00
23.87 ± 0.03
+3.38 ± 0.11
33.65 ± 0.02
37.22 ± 0.01
+10.61 ± 0.08
FP-Soma
23.07 ± 0.00
23.86 ± 0.03
+3.42 ± 0.09
33.65 ± 0.02
37.08 ± 0.10
+10.20 ± 0.26
Soma
22.95 ± 0.01
23.86 ± 0.03
+3.96 ± 0.12
33.66 ± 0.01
37.21 ± 0.08
+10.55 ± 0.19
Rsm
22.96 ± 0.01
23.87 ± 0.01
+3.98 ± 0.02
33.61 ± 0.02
37.34 ± 0.10
+11.13 ± 0.26
Table 1: INT4 post-training quantization results on held-out FineWeb-Edu perplexity (lower is better), comparing full precision (FP) and 4-bit AWQ (INT4). Results show mean ± std over 3 pretraining seeds (INT4 averaged over 3 calibration seeds); bold indicates the lowest value per column. Row colour marks the base optimizer family: Muon -based ( Green ), AdamW -based ( Red ), SGD -based ( Blue ).
Appendix figures & tables22 assets
Supplementary material from the paper’s appendix.
Appendix
Algorithm
α=1.1
α=1.6
α=2.0
α=3.0
best acc.
ρ⋆
best acc.
ρ⋆
best acc.
ρ⋆
best acc.
ρ⋆
Muon
38.1
–
42.9
–
44
–
44.3
–
Ssm
38.6
0.5
43.7
0.5
46.3
0.5
47
0.5
FP-Soma
39.8
1
44.2
0.5
46.9
0.5
46.7
0.5
Soma
39.1
1
45
1
46.3
1
46.5
0.5
Rsm
38.3
0.01
43.3
0.002
44.5
0.1
44.8
1
Appendix
Table 2: Validation accuracy (%): best value over the perturbation radius ρ and the radius ρ⋆ where it is attained (median of 8 seeds; Frobenius-ball methods pool the global and per-layer radius modes). Row colour marks the base optimizer family: Muon -based ( Green ), AdamW -based ( Red ), SGD -based ( Blue ); baselines have no ρ⋆ . Best value per α in bold.
Algorithm
α=1.1
α=1.6
α=2.0
α=3.0
best loss
ρ⋆
best loss
ρ⋆
best loss
ρ⋆
best loss
ρ⋆
Muon
1.45
–
0.372
–
0.184
–
0.078
–
Ssm
1.45
0.02
0.372
0.002
0.183
0.001
0.0779
0.001
FP-Soma
1.44
0.2
0.372
0.01
0.183
0.001
0.0768
0.02
Soma
1.43
0.2
0.371
0.02
0.184
0.001
0.0783
0.01
Rsm
1.45
0.1
0.371
0.2
0.184
0.001
0.0779
0.05
Appendix
Table 3: Training loss: best value over the perturbation radius ρ and the radius ρ⋆ where it is attained (median of 8 seeds; Frobenius-ball methods pool the global and per-layer radius modes). Row colour marks the base optimizer family: Muon -based ( Green ), AdamW -based ( Red ), SGD -based ( Blue ); baselines have no ρ⋆ . Best value per α in bold.
Figure 3: Performance under anisotropic heavy-tailed noise (Condition no., κε=10 ). Optimality gap (log scale) vs ρ . Red , Blue and Green lines indicate AdamW -based, SGD -based and Muon -based methods, respectively.
Algorithm
α=1.1
α=1.6
α=2.0
α=3.0
gap
ρ⋆
gap
ρ⋆
gap
ρ⋆
gap
ρ⋆
Muon
0.413
–
0.062
–
0.0316
–
0.0165
–
Ssm
0.265
1
0.0473
0.5
0.0253
0.5
0.0141
0.5
FP-Soma
0.252
1
0.0457
0.5
0.0247
0.5
0.0143
0.2
Soma
0.263
1
0.0551
0.5
0.0299
0.2
0.0163
0.1
Rsm
0.408
1
0.0609
1
0.0311
1
0.0162
1
Appendix
Table 4: Optimality gap: best value over the perturbation radius ρ and the radius ρ⋆ where it is attained (median of 8 seeds; Frobenius-ball methods pool the global and per-layer radius modes). Row colour marks the base optimizer family: Muon -based ( Green ), AdamW -based ( Red ), SGD -based ( Blue ); baselines have no ρ⋆ . Best value per α in bold.
Figure 4: Performance under isotropic heavy-tailed noise (nonlinear setting)(Condition no., κϵ=1 ). Top: validation accuracy. Bottom: training loss (log scale). Red , Blue and Green lines indicate AdamW -based, SGD -based and Muon -based methods, respectively.
Algorithm
α=1.1
α=1.6
α=2.0
α=3.0
best acc.
ρ⋆
best acc.
ρ⋆
best acc.
ρ⋆
best acc.
ρ⋆
Muon
36.1
–
43.3
–
45.2
–
43.8
–
Ssm
36.6
0.2
44.2
0.5
45.5
0.005
46.2
0.5
FP-Soma
37.8
1
44.7
1
46.6
0.5
46.1
0.5
Soma
38.2
1
44.8
1
46.6
0.5
46.6
0.5
Rsm
36.3
0.05
43.6
0.1
45.2
0.2
44.8
0.2
Appendix
Table 5: Validation accuracy (%): best value over the perturbation radius ρ and the radius ρ⋆ where it is attained (median of 8 seeds; Frobenius-ball methods pool the global and per-layer radius modes). Row colour marks the base optimizer family: Muon -based ( Green ), AdamW -based ( Red ), SGD -based ( Blue ); baselines have no ρ⋆ . Best value per α in bold.
Algorithm
α=1.1
α=1.6
α=2.0
α=3.0
best loss
ρ⋆
best loss
ρ⋆
best loss
ρ⋆
best loss
ρ⋆
Muon
1.58
–
0.482
–
0.233
–
0.114
–
Ssm
1.58
0.02
0.482
0.001
0.236
0.001
0.113
0.01
FP-Soma
1.57
0.2
0.481
0.001
0.234
0.02
0.112
0.02
Soma
1.56
0.5
0.483
0.002
0.234
0.001
0.113
0.01
Rsm
1.58
0.005
0.482
0.05
0.233
0.5
0.113
0.2
Appendix
Table 6: Training loss: best value over the perturbation radius ρ and the radius ρ⋆ where it is attained (median of 8 seeds; Frobenius-ball methods pool the global and per-layer radius modes). Row colour marks the base optimizer family: Muon -based ( Green ), AdamW -based ( Red ), SGD -based ( Blue ); baselines have no ρ⋆ . Best value per α in bold.
Figure 5: Performance under isotropic heavy-tailed noise (linear setting)(Condition no., κϵ=1 ). Optimality gap (log scale) vs ρ . Red , Blue and Green lines indicate AdamW -based, SGD -based and Muon -based methods, respectively.
Algorithm
α=1.1
α=1.6
α=2.0
α=3.0
gap
ρ⋆
gap
ρ⋆
gap
ρ⋆
gap
ρ⋆
Muon
1.01
–
0.231
–
0.114
–
0.0552
–
Ssm
0.648
1
0.141
1
0.0733
1
0.0412
0.5
FP-Soma
0.654
1
0.147
1
0.0823
0.5
0.0413
0.5
Soma
0.63
1
0.176
0.5
0.086
0.5
0.0486
0.5
Rsm
1
1
0.224
1
0.108
1
0.0527
1
Appendix
Table 7: Optimality gap: best value over the perturbation radius ρ and the radius ρ⋆ where it is attained (median of 8 seeds; Frobenius-ball methods pool the global and per-layer radius modes). Row colour marks the base optimizer family: Muon -based ( Green ), AdamW -based ( Red ), SGD -based ( Blue ); baselines have no ρ⋆ . Best value per α in bold.
Tail index α
Direction
Geometry
1.1
1.6
2.0
3.0
Random
Frobenius
38.3
43.4
44.5
44.7
Spectral
38.5
44.7
46.5
47.2
Sam - Muon
Frobenius
38.3
43.2
45.0
45.6
Spectral
38.6
43.7
46.3
47.0
Soma
Frobenius
38.2
43.3
45.0
45.6
Appendix
Table 8: Nonlinear anisotropic setting: best validation accuracy (%, ↑ ) with a Frobenius or a spectral perturbation, for each tail index α . The better geometry for each method and α is in bold; ties are not marked.
Tail index α
Direction
Geometry
1.1
1.6
2.0
3.0
Random
Frobenius
0.410
0.0612
0.0312
0.0163
Spectral
0.347
0.0576
0.0301
0.0163
Sam - Muon
Frobenius
0.404
0.0581
0.0294
0.0149
Spectral
0.265
0.0473
0.0253
0.0141
Soma
Frobenius
0.399
0.0592
0.0303
0.0163
Appendix
Table 9: Linear anisotropic setting: best optimality gap ( ↓ ) with a Frobenius or a spectral perturbation, for each tail index α . The better geometry for each method and α is in bold; ties are not marked.
Tail index α
Geometry
Direction
Method
1.1
1.6
2.0
3.0
Frobenius
Fresh
Sam - Muon
0.3828
0.4316
0.4502
0.4561
Stale
S 2 -Sam - Muon
0.3818
0.4326
0.4502
0.4492
Spectral
Fresh
Ssm
0.3857
0.4365
0.4629
0.4697
Stale
S 2 -Sam - Muon
0.3887
0.4443
0.4678
0.4570
Spectral
Fresh
FP-Soma
0.3984
0.4424
0.4688
0.4668
Appendix
Table 10: Nonlinear anisotropic setting: best validation accuracy (%, ↑ ) with fresh or stale perturbation direction, for each tail index α . The better direction for each geometry and α is in bold; ties are not marked.
Tail index α
Geometry
Direction
Method
1.1
1.6
2.0
3.0
Frobenius
Fresh
Sam - Muon
0.4037
0.0581
0.0294
0.0149
Stale
S 2 -Sam - Muon
0.4068
0.0601
0.0312
0.0165
Spectral
Fresh
Ssm
0.2648
0.0473
0.0253
0.0141
Stale
S 2 -Sam - Muon
0.2846
0.0588
0.0309
0.0165
Spectral
Fresh
FP-Soma
0.2522
0.0457
0.0247
0.0143
Appendix
Table 11: Linear anisotropic setting: best optimality gap ( ↓ ) with a fresh or a stale perturbation direction, for each tail index α . The better direction for each geometry and α is in bold; ties are not marked.
FineWeb-Edu NanoGPT
FineWeb-Edu Pythia
CIFAR-10 Airbench CNN
Optimizer
Train loss
Val ppl
Train loss
Val ppl
Train loss
Val acc (%)
Muon
3.303 ± 0.004
24.41 ± 0.02
3.558 ± 0.003
35.72 ± 0.03
0.266 ± 0.002
87.27 ± 0.08
Ssm
3.305 ± 0.004
24.55 ± 0.03
3.559 ± 0.003
35.78 ± 0.03
0.271 ± 0.003
87.63 ± 0.34
FP-Soma
3.305 ± 0.004
24.54 ± 0.03
3.560 ± 0.003
35.79 ± 0.03
0.270 ± 0.001
87.37 ± 0.44
Soma
3.303 ± 0.004
24.40 ± 0.01
3.562 ± 0.003
35.79 ± 0.03
0.269 ± 0.002
87.51 ± 0.27
Rsm
3.303 ± 0.004
24.41 ± 0.02
3.558 ± 0.003
35.74 ± 0.03
0.266 ± 0.004
87.22 ± 0.07
Appendix
Table 12: Continued-pretraining results of the fifteen optimizers after the last epoch (FineWeb-Edu NanoGPT 15, FineWeb-Edu Pythia 14, CIFAR-10 Airbench CNN 7 epochs). Lower is better for train loss and validation perplexity (Val ppl), higher for validation accuracy (Val acc); the best validation value of each dataset is in bold. Mean ± std over 3 seeds. Row colour marks the base optimizer family: Muon -based ( Green ), AdamW -based ( Red ), SGD -based ( Blue ).
Figure 6: Catastrophic forgetting of nanoGPT (GPT-2 small) models pretrained with each optimizer, after fine-tuning on four datasets. Bars: rise Δ PT of the held-out FineWeb-Edu loss at a common fine-tuning loss (lower is better); diamond: average over datasets; error bars: standard deviation over 3 pretraining seeds. † starting model clearly weaker than AdamW’s (not comparable).
Figure 7: Catastrophic forgetting of Pythia-70M models pretrained with each optimizer, after fine-tuning on four datasets. Bars: rise Δ PT of the held-out FineWeb-Edu loss at a common fine-tuning loss (lower is better); diamond: average over datasets; error bars: standard deviation over 3 pretraining seeds. † starting model clearly weaker than AdamW’s (not comparable).
Lafter by fine-tuning dataset
Optimizer
Lbefore
codeparrot
stackmathqa
musicpile
tulu3
Average
Δ%adj
Muon
3.144 ± 0.000
3.516 ± 0.074
3.526 ± 0.004
3.259 ± 0.000
3.332 ± 0.001
3.409 ± 0.020
4.20 ± 0.31
Ssm
3.147 ± 0.000
3.386 ± 0.003
3.561 ± 0.006
3.251 ± 0.001
3.331 ± 0.040
3.382 ± 0.009
3.73 ± 0.15
FP-Soma
3.146 ± 0.000
3.393 ± 0.004
3.515 ± 0.084
3.251 ± 0.001
3.308 ± 0.002
3.367 ± 0.022
3.50 ± 0.35
Soma
3.144 ± 0.001
3.491 ± 0.086
3.511 ± 0.005
3.258 ± 0.002
3.330 ± 0.001
3.397 ± 0.023
4.03 ± 0.35
Rsm
3.144 ± 0.001
3.558 ± 0.009
3.526 ± 0.007
3.260 ± 0.000
3.333 ± 0.001
3.419 ± 0.004
4.36 ± 0.06
Appendix
Table 13: Catastrophic Forgetting of nanoGPT (GPT-2 small) models pretrained with different optimizers: held-out FineWeb-Edu loss before ( Lbefore ) and after ( Lafter ) fine-tuning on four datasets (lower is better). Δ%adj is the smoothed relative increase of the average Lafter . Mean ± std over 3 seeds. Row colour marks the base optimizer family: Muon -based ( Green ), AdamW -based ( Red ), SGD -based ( Blue ).
Lafter by fine-tuning dataset
Optimizer
Lbefore
codeparrot
stackmathqa
musicpile
tulu3
Average
Δ%adj
Muon
3.515 ± 0.000
4.361 ± 0.005
4.464 ± 0.015
3.732 ± 0.001
3.888 ± 0.056
4.111 ± 0.013
8.45 ± 0.19
Ssm
3.516 ± 0.001
4.360 ± 0.002
4.478 ± 0.026
3.733 ± 0.001
3.856 ± 0.004
4.107 ± 0.006
8.36 ± 0.09
FP-Soma
3.516 ± 0.001
4.361 ± 0.002
4.480 ± 0.016
3.733 ± 0.001
3.852 ± 0.001
4.107 ± 0.004
8.36 ± 0.07
Soma
3.517 ± 0.000
4.360 ± 0.004
4.481 ± 0.018
3.734 ± 0.001
3.855 ± 0.004
4.107 ± 0.003
8.36 ± 0.05
Rsm
3.514 ± 0.000
4.358 ± 0.007
4.468 ± 0.027
3.732 ± 0.001
3.857 ± 0.002
4.104 ± 0.008
8.35 ± 0.12
Appendix
Table 14: Catastrophic Forgetting of Pythia-70M models pretrained with different optimizers: held-out FineWeb-Edu loss before ( Lbefore ) and after ( Lafter ) fine-tuning on four datasets (lower is better). Δ%adj is the smoothed relative increase of the average Lafter . Mean ± std over 3 seeds. Row colour marks the base optimizer family: Muon -based ( Green ), AdamW -based ( Red ), SGD -based ( Blue ).
Optimizer
FineWeb-Edu NanoGPT
FineWeb-Edu Pythia
CIFAR-10 Airbench CNN
s / epoch
s / epoch
s / epoch
Muon
265.2 ± 3.6
287.2 ± 2.0
23.6 ± 0.3
Ssm
772.5 ± 6.7
236.8 ± 1.9
23.9 ± 0.4
FP-Soma
641.8 ± 129.0
240.8 ± 4.2
23.9 ± 0.6
Soma
272.1 ± 1.8
294.4 ± 4.8
23.2 ± 2.1
Rsm
267.2 ± 1.1
298.2 ± 2.6
23.1 ± 0.5
Appendix
Table 15: Pretraining runtime of the optimizers: wall-clock seconds per epoch (100 optimizer steps). Row colour marks the base optimizer family: Muon -based ( Green ), AdamW -based ( Red ), SGD -based ( Blue ).
Optimizer
FineWeb-Edu NanoGPT
FineWeb-Edu Pythia
s / run
s / run
Muon
48.6 ± 18.2
25.4 ± 7.0
Ssm
82.0 ± 33.7
58.6 ± 21.4
FP-Soma
84.2 ± 49.9
50.9 ± 20.2
Soma
36.0 ± 2.5
27.8 ± 6.8
Rsm
31.2 ± 0.1
25.7 ± 5.8
Appendix
Table 16: Fine-tuning runtime of the optimizers: wall-clock seconds per fine-tuning run, taken from the fixed-base fine-tuning study (per job the median over its successful runs, then mean ± std over 3 pretrained bases × 3 seeds). Runs shared the GPUs with other jobs, so only large differences are meaningful. Row colour marks the base optimizer family: Muon -based ( Green ), AdamW -based ( Red ), SGD -based ( Blue ).
nanoGPT
Pythia-70M
CIFAR-10 CNN
Optimizer
lr
wd
ρ
lr
wd
ρ
lr
wd
ρ
Muon
4.24×10−3
0.01
–
8.48×10−4
0.3
–
0.03
0
–
Ssm
4.24×10−3
0.01
0.01
8.48×10−4
0.3
3×10−4
0.03
0
0.01
FP-Soma
4.24×10−3
0.01
0.01
8.48×10−4
0.3
3×10−4
0.03
0
0.01
Soma
4.24×10−3
0.01
3×10−4
8.48×10−4
0.3
3×10−4
0.03
0
3×10−3
Rsm
4.24×10−3
0.01
0.01
8.48×10−4
0.3
9.5×10−4
0.03
0
0.1
Appendix
Table 17: Selected learning rate, weight decay and perturbation radius ρ of the continued-pretraining experiments.
nanoGPT
Pythia-70M
Optimizer
lr
ρ
lr
ρ
Muon
10−2
–
10−3
–
Ssm
10−2
10−2
10−3
3×10−4
FP-Soma
10−2
10−2
10−3
3×10−4
Soma
10−2
3×10−4
10−3
3×10−4
Rsm
10−2
10−2
10−3
9.5×10−4
Appendix
Table 18: Selected learning rate and perturbation radius ρ for fine-tuning.