Preprint · 2026

Beyond SAEs: Consistent Concept Discovery with Sparse Coding

Sparse autoencoders trained on the same activations can produce different dictionaries and coefficients across random seeds. We propose Consistent Sparse Dictionary Learning (ConSDL), which replaces the learned encoder with per-sample sparse coding to improve the consistency of both dictionaries and coefficients across training runs.

1Korea Advanced Institute of Science & Technology (KAIST) 2Helmholtz Munich
+17.6%
dictionary consistency
relative gain over baseline average
+35.9%
coefficient consistency
relative gain over baseline average
16
synthetic settings
near-perfect dictionary and coefficient recovery
4
model backbones
CLIP · DINOv2 · GPT-2 · Pythia

Consistency gains are averaged across four backbones relative to six SAE baselines, with nine dictionary-width and sparsity settings per backbone.

Cross-seed consistency and dictionary utilization across four backbones
More consistent dictionaries and coefficients across training runs. We compare ConSDL with six SAE baselines on CLIP, DINOv2, Pythia-160M, and GPT-2 Small, across 4×, 8×, and 16× dictionary expansion and sparsity levels ℓ0 ≈ 32, 64, and 128. ConSDL (red) improves cross-seed consistency on average while maintaining high dictionary utilization.

Method

ConSDL replaces the learned encoder with per-sample sparse coding, alternating between coefficient inference and dictionary updates. Coefficient refitting and co-active similarity regularization stabilize dictionary learning by correcting shrinkage bias from the sparsity penalty and reducing redundancy among co-active atoms.

ConSDL overview and training-step pseudocode
1

Per-sample sparse coding

For each model activation, FISTA iteratively optimizes nonnegative sparse coefficients under the current dictionary. Each iteration takes a gradient step on the reconstruction term, followed by nonnegative soft thresholding for the ℓ₁ penalty.

$$\mathbf{a}^{\mathrm{sc}}=\arg\min_{\mathbf{a}\ge0}\tfrac12\lVert\mathbf{x}-W_{\mathrm{dec}}\mathbf{a}-\mathbf{b}_{\mathrm{dec}}\rVert_2^2+\lambda\lVert\mathbf{a}\rVert_1$$
2

Coefficient refitting

The ℓ₁ penalty promotes sparsity but also shrinks nonzero coefficients, potentially underestimating selected atoms’ contributions during dictionary updates. ConSDL refits their magnitudes on the support selected by FISTA without this penalty, while retaining nonnegativity.

$$\mathbf{a}^{\mathrm{refit}}=\arg\min_{\mathbf{a}\ge0,\;a_j=0\;\forall j\notin S^{\mathrm{sc}}}\tfrac12\lVert\mathbf{x}-W_{\mathrm{dec}}\mathbf{a}-\mathbf{b}_{\mathrm{dec}}\rVert_2^2$$
3

Co-active similarity regularization

Although sparse coding improves coefficient estimation, atoms can still converge to similar directions during dictionary learning. ConSDL reduces this redundancy by penalizing high cosine similarity between each active atom and its most similar atom active for the same input.

$$\mathcal{L}_{\mathrm{sim}}=\frac{1}{|V_\mu|}\sum_{(i,j)\in V_\mu}(s_{i,j}-\mu)^2$$

Why refitting and regularization matter

Removing either coefficient refitting or co-active similarity regularization makes dictionary learning less stable. At 16× expansion and matched sparsity (ℓ₀ ≈ 128), refitting raises reconstruction R² from 0.87 to 0.92 on CLIP and eliminates the dead atoms observed with sparse coding alone on Pythia-160M. However, refitting alone still leaves thousands of duplicate atoms on both backbones. Adding similarity regularization reduces these duplicates to nearly zero, yielding high consistency with a fully usable dictionary.

Recovering ground-truth dictionaries and coefficients

We test whether ConSDL improves dictionary and coefficient recovery over SAE baselines as dictionary size and frequency imbalance increase. We generate 768-dimensional synthetic data with an expected 60 active atoms per sample and log-normal nonzero coefficients. Across 16 settings, dictionary sizes range from 3,072 to 12,288 atoms and Zipf exponents from β = 0 (uniform selection) to β = 1 (strong frequency imbalance).

Synthetic data generating process
Synthetic data with a known decomposition. Observations are sparse linear combinations of unit-norm dictionary atoms. The plots show the distributions of support size, nonzero coefficient magnitudes, and atom-selection frequencies.

SAE recovery is generally strong for smaller dictionaries and more balanced activation frequencies, but declines as the settings become more challenging. ConSDL achieves near-perfect dictionary and coefficient recovery across all 16 settings.

Ground-truth recovery and cross-seed consistency
Ground-truth recovery and cross-seed consistency. SAE recovery declines as dictionary size and frequency imbalance increase, while ConSDL maintains near-perfect dictionary and coefficient recovery. Cross-seed consistency strongly correlates with ground-truth recovery (Pearson r = 0.90 for dictionaries and r = 0.99 for coefficients). These correlations support its use when ground truth is unavailable, although consistency alone does not establish recovery.

Results on vision and language models

We train ConSDL and six SAE baselines on penultimate-layer activations from CLIP ViT-B/16 and DINOv2-B/14 using ImageNet, and on layer-9 activations from GPT-2 Small and Pythia-160M using 500M FineWeb tokens. We evaluate cross-seed consistency at dictionary expansions of 4×, 8×, and 16× and sparsity levels ℓ₀ ≈ 32, 64, and 128, with three training seeds per configuration.

Dictionary consistency (DCS) measures the mean decoder cosine similarity after Hungarian matching. Coefficient consistency (MCC) measures the mean coefficient correlation under the same matching. Cross-seed scores are averaged over all seed pairs. The usable atom fraction measures how much of the dictionary is neither dead nor duplicated.

ConSDL achieves the highest average DCS and MCC across the four backbones. MP-SAE has comparable average dictionary consistency but substantially lower coefficient consistency. ConSDL improves both metrics while maintaining high dictionary utilization.

Analysis

Finer-grained concepts

On CLIP ViT-B/16 (16×, ℓ0 ≈ 128), ConSDL learns about 9K consistent atoms (cross-seed cosine > 0.7), compared with 4K for ReLU. Where ReLU learns one broad atom, ConSDL often learns several matched atoms that each capture a more specific feature.

A ReLU beer-bottle atom and its matched ConSDL atoms in CLIP ViT-B/16
One ReLU atom, several ConSDL atoms. A ReLU “beer bottle” atom matches five ConSDL atoms (decoder cosine > 0.7) for bottle necks, bubbles, bottle shapes, and beer in glasses. The ConSDL atoms differ from each other (pairwise cosine ~0.5) and respond to different image regions.

Larger dictionaries, more detailed concepts

Going from 4× to 16× expansion, a broad ConSDL atom splits into several more specific atoms, while ReLU and TopK keep nearly the same broad atom. Across all four backbones, consistent ConSDL atoms at 4× more often have multiple counterparts at 16× than ReLU or TopK atoms.

A 4× ConSDL retriever atom and its 16× counterparts in CLIP ViT-B/16
4× vs. 16× dictionaries. A 4× ConSDL “Golden Retriever” atom corresponds to six 16× atoms for fur, tails, muzzle shapes, and puppies. ReLU and TopK keep a nearly identical retriever atom (cosine 0.99 and 0.98).

Consistency beyond individual atoms

We examine whether atom-level consistency extends to the overall structure of the dictionary. For CLIP ViT-B/16, we cluster atoms based on similarity between decoder directions and compare cluster sizes and centroids across training runs.

Dictionary cluster consistency across training runs
Dictionary cluster consistency across training runs. (A) Joint UMAP of ConSDL decoder atoms at 16×, ℓ0 ≈ 128, colored by top 30 clusters from a reference run. (B) Cluster-size profiles for ConSDL and Matryoshka. ρ denotes the Spearman correlation of cluster sizes across runs. (C) Normalized cluster-centroid agreement across granularities. Lines show the mean over nine CLIP settings, and shading shows the interquartile range.

Consistent activation patterns

Qualitative examples in CLIP and Pythia-160M show how matched atoms activate across training runs. In CLIP, matched atoms highlight similar image regions and retrieve similar visual examples. In Pythia, matched atoms highlight similar tokens and contexts, including animal nouns and sentence-final punctuation.

Region decomposition and cross-seed latents in CLIP ViT-B/16
Region decomposition and cross-seed latent analysis in CLIP ViT-B/16. At 16× expansion and ℓ0 ≈ 128, (A) image patches are independently clustered by coefficient vectors across three seeds. (B–C) Activation maps and top-activating images compare selected latents with their nearest cross-seed matches, showing animals in close contact (B) and a heart/romance motif (C).

Interactive Demo

Explore how matched dictionary atoms activate across training seeds and methods.

Coming soon

BibTeX

@misc{consdl2026,
  title  = {Beyond SAEs: Consistent Concept Discovery with Sparse Coding},
  author = {Choi, Jinho and Lim, Hyesu and Choo, Jaegul and Schneider, Steffen},
  year   = {2026},
  eprint = {2610.23801},
  archivePrefix = {arXiv},
  url    = {https://arxiv.org/abs/2610.23801}
}