Iterative Orthogonalization Scaling Laws
Devan Selvaraj
TL;DR
The paper investigates potential scaling limitations of muon's iterative orthogonalization as model sizes grow, showing that singular-value distributions shift with scale and may hamper large-scale deployments. It combines empirical experiments on randomly initialized matrices with Marchenko-Pastur-based theory to characterize how singular values shrink with dimension and how NS iterations affect this behavior. A key contribution is articulating a scaling law for singular values with dimension and discussing tuning requirements for NS iterations at scale, without proposing concrete remedies. The findings alert practitioners to possible scaling challenges in very large models and motivate further work on robust orthogonalization strategies and precision considerations in large-scale settings.
Abstract
The muon optimizer has picked up much attention as of late as a possible replacement to the seemingly omnipresent Adam optimizer. Recently, care has been taken to document the scaling laws of hyper-parameters under muon such as weight decay and learning rate. However, at much larger scales the iterative orthogonalization procedure present in muon may suffer a possible issue as the singular values of random matrices shrink with scale. This paper shows this scaling behavior theoretically and empirically on random matrices but does not suggest what to do about it.
