With a fixed sample budget T, a larger batch size B reduces gradient noise but also reduces the number of updates. The paper separates these effects using a noisy quadratic model.[1]
Model
f(w)=21i∑hiwi2,gi=hiwi+Bcizi,K=BT.
hi>0: curvature; ci>0: per-sample gradient-noise variance; zi∼N(0,1). The direction-dependent curvature-to-noise ratio is CNRi=hi2/ci.
Ideal Muon reduces to a sign operation on a 1×1 matrix. Start with SignSGD with momentum:
Directional intervention. Starting from the same checkpoint, increase the batch from 128K to 2M, retaining small-batch updates only in selected directions. The recovered fraction of the local loss gap is
Metric
R=Llarge−LsmallLlarge−Lmixed.
Directions retaining small-batch updates
Local gap recovered
The 768 sharpest directions, about 0.001%
Up to 59.5%
Random subspace of the same dimension
Change of at most about 3%
Local curvature diagnostic on the same data stream. Sharp directions come from a Muon-preconditioned curvature approximation.[2]
Related updateAbstract only
Noise can smooth the sign—and change the order of feature learning
The authors report that small-batch noise restores sequential feature learning absent from full-batch spectral descent, and formalize the mechanism in deep linear models.[3][4]
The following illustrative scalar derivation explains noise smoothing. For z∼N(0,1) and σ>0,
Chinchilla (2022) used small-scale training runs to study compute allocation between model size and data. SLDBench (2025; ICLR 2026) then asked agents to discover predictive laws automatically. The new question: obtaining the fitting data itself costs compute.[6][7]
Candidate law
L(N,D)=E+aN−α+bD−β.
The expression above is only a candidate form. The agent neither knows the correct functional form nor gets all observations for free; it receives a catalog of priced experiments.[5]
Experiment catalog
U={(xi,ci)}i=1n,ci>0.
xi is a training configuration; ci is usually its training FLOPs. The outcome yi is revealed only after purchase. Omitting group indices, the agent buys At each round and stops at time τ:
Spending
St+1=St∪At,ρ=∑i∈Uci∑i∈Sτci.
Evaluation records (extrapolation error at larger scales,ρ). The agent must jointly decide which experiments to buy, when to stop, and which law to fit. Test configurations are larger than candidate experiments and cannot be purchased directly.[5]
Setting
Range of median R2 across systems
Experiment cost
Passive · all observations
0.723–0.753
100%
Active · chosen purchases
0.235–0.570
About 8%–44%, averaged within each system
30 selected difficult tasks and 8 agent systems. The R² range is across systems, not a confidence interval.[5]
BriefAbstract only
Monarch-Muon: break a large orthogonalization into smaller blocks
Monarch factors impose structure on the weights, reducing weight, gradient, and optimizer-state storage while producing smaller update subproblems. The abstract reports these measurements for a 6.89B system:[8][9]
Measurement
Original setup
Monarch-Muon
Peak memory
68.9 GB
52.9 GB
Muon optimizer step
1288 ms
304 ms
Daily classic · 001Relevant method, theory, and experiment sections read
Long sequences require both memory and efficient training. HiPPO (2020) designed continuous memory through polynomial projections; S4 (2022) sought a state-space parameterization that is fast to compute and numerically stable.[10][13]
Consider one channel with scalar input and output, and hidden state x(t)∈RN. A governs memory evolution, B writes the input, and C reads the state. Omit the separately computable skip term.[12]
Continuous system
x˙(t)=Ax(t)+Bu(t),y(t)=Cx(t).
Apply a bilinear discretization with step size Δ>0, assuming the inverse exists. Δ is the sampling timescale and can be learned as a model parameter.[12]
Discretization
Aˉ=(I−2ΔA)−1(I+2ΔA),Bˉ=(I−2ΔA)−1ΔB.
With zero initial state x−1=0 and fixed parameters, stepwise recurrence and full-sequence convolution are exactly equivalent:[12]
Computing the kernel is the real bottleneck. Direct eigendecomposition of a HiPPO matrix can be severely ill-conditioned. S4 uses a normal matrix plus a low-rank correction, then a unitary change of basis gives DPLR:[10]
Structured parameterization
A=V(Λ−PQ∗)V∗,V∗V=I,P,Q∈CN×r,r≪N.
Λ is diagonal. Write A=Λ−PQ∗ and R(ω)=(ωI−Λ)−1. When the relevant inverses exist, Woodbury’s identity gives:[10]
Low-rank correction
(ωI−A)−1=R−RP(I+Q∗RP)−1Q∗R.
A high-dimensional inverse becomes a diagonal inverse plus an r×r correction. The paper reduces the finite-length kernel’s generating function to Cauchy computations and an inverse FFT. For fixed low rank, theoretical kernel construction costs O(N+L), and each recurrent step costs O(N). These are per-SSM costs and assume a fast Cauchy algorithm.[10]
Path-X: connectivity classification over length 16,384
Accuracy
Random guessing
50.00%
S4 original
88.10%
The “original” row in v3 Appendix Table 10; later improvements are listed separately as “updated.” Original Path-X setup: 6 layers, 256 channels, state dimension 64 per channel, one A100 40GB.[10]