π€ AI Summary
This work addresses the lack of a clear theoretical characterization of the βconceptsβ extracted by sparse autoencoders (SAEs), particularly regarding their applicability to complex language model representations. Bypassing assumptions about data-generating models, the study extends local optimality analysis for the first time to the non-negative SAE framework. By linking non-negative jointly optimized solutions to the underlying feature distribution, it reveals how L1 regularization and non-negativity jointly shape dictionary structure. This analysis yields precise constraints between SAE features and data distribution, successfully explaining empirical phenomena such as hierarchical splitting, feature absorption, residual structures, and dense antipodal features. The results provide a theoretical foundation for SAE interpretability and offer principled guidance for future architecture design.
π Abstract
Sparse Autoencoders (SAEs) have found success parsing neural representations into interpretable concepts, providing a basis for understanding and control. However, what exactly SAEs extract, and, correspondingly, the scientific conclusions we can draw from them, are not obvious. Empirically, the proof is in the pudding: SAEs learn interpretable features. Theoretically, we lack a clear account of what properties a 'concept' must satisfy for an SAE to extract it. There has been extensive identifiability work studying the conditions under which sparse coding recovers ground-truth features; however, these approaches tends to focus on simple data-generating models (e.g. sparse independent features) which poorly approximate the internet-swallowing language-model representations on which SAEs are trained. Here, avoiding data-generating models, we ask simply what properties any dictionary learning optimum must satisfy. Concretely, we extend local optimality analyses (Gribonval & Schnass, 2010) to the nonnegative joint-optimisation problem that vanilla SAEs approximate, and derive constraints relating optimal SAE features to their distributions. We use these constraints to explain a range of observed SAE behaviours - hierarchical splitting & absorption, the structure of residuals, and dense antipodal features - each reflecting how L1+nonnegativity interact with data to structure optimal dictionaries. Finally, we construct a novel large-dictionary convex problem and explore the wide atom-per-datapoint limit. In sum, we hope to tease model assumptions from unexpected observations, letting us learn more from SAEs' successes and provide principles for designing their successors.