← SoftRank overview

SoftRank · How it works

Measure what softmax can see, in no particular coordinates

The rank an attention head needs is a property of the attention operator, not of the query and key vectors that happen to produce it. SoftRank defines that rank by the outputs it must preserve, bounds it with geometry that survives every change of coordinates, and brackets it directly on trained heads.

Quantitative Research Society, Singapore · Nanyang Technological University

Part 1 · The idea

1Why not measure the queries and keys?

An attention head computes logits H=βQK⊤ and normalises each row with softmax. The obvious way to say how complex the head is, is to measure the two point clouds that produce the logits: how far the queries and keys spread (their radii, combined into Γ=βρQρK), or the spectrum of Wq⊤Wk. A bound written in those terms is correct for the factorization it is given.

But the factorization is not unique. For any invertible matrix S, replacing Q by QS and K by KS−⊤ leaves every logit, and so the whole attention matrix, unchanged. The point clouds do change. In the paper’s rank-two example, four queries and keys at ±e1,±e2 give Γ=1; after the shear St with t=2,5,10 the same attention matrix has Γ=5,26,101. A measure that can be made as large as we like without touching the operator describes the coordinates, not the head.

Two ways to measure attention complexityTop: radii and spectra of one query-key factorization, which change under an invertible change of coordinates. Bottom: the attention operator, its softmax-visible interaction C, a retained part T with rank, residual and size, and a bound on the approximation rank that is the same for every factorization.Usual: measure one factorizationQ and Kone choice ofquery and key vectorsRadii, spectraρ_Q, ρ_K, spectraof W_qᵀW_kΓ = βρ_Qρ_Knative sizeof the interactionDepends on the frameQ→QS, K→KS⁻ᵀ moves Γwhile A stays fixedSoftRank: measure the operatorA = softmax(QKᵀ)the normalisedattention operatorC = P(log A)Pwhat row softmaxcan seeRetain T ≈ Crank r, residual τ(T),size κ(T)Bound on r_ε(A)the same for everyrefactorization
Two ways to measure attention complexityLeft: one factorization and its radius product, which depends on the coordinate frame. Right: the softmax-visible interaction C, a retained part T, and a bound that is the same for every factorization.Usual: one factorizationQ and Kone choice ofquery and keyΓ = βρ_Qρ_Kradii of thetwo cloudsDepends on the frameQ→QS, K→KS⁻ᵀmoves ΓSoftRank: the operatorC = P(log A)Pwhat softmaxcan seeRetain T ≈ Crank r, residual,size κ(T)Bound on r_ε(A)same for everyfactorization
Two ways to measure. Above, quantities of one query–key factorization, which move under a change of coordinates. Below, quantities of the normalised operator, which do not.

SoftRank changes both ends. The target is the operator’s own approximation rank. The geometry is read off the part of the logits that softmax can see, and its size is minimised over every factorization at once, so that no frame is privileged.

2Three questions, three objects

What must be preserved?

The approximation rank rε(A)

The least rank of any real matrix B within maximum-row-ℓ1 error ε of A. That error is exactly the worst output error over all value vectors of norm at most one, in any normed space, so rε(A) is the rank needed to preserve every bounded attention output.

What can softmax see?

The visible interaction C=Pm(logA)PN

Row softmax ignores any constant added to a row of logits, and terms that depend only on the key act as common base weights. Centring logA on both sides removes the first and keeps the second apart. C is a function of A alone, and R=rankC is the exact visible interaction dimension.

How large is it?

Residual τ(T) and size κ(T)

Keep a rank-r part T of C. The residual τ(T) is the output error that dropping the rest can cause; the size κ(T) is the smallest product of the two enclosing radii over every factorization T=XY⊤. Both are fixed by A and T, whatever coordinates produced them.

Put together, they give the central theorem: whenever τ(T)<ε,

rε(A)≤min{p,Cr(1+κ(T)(ε−τ(T))2)r/2},

Here p=min(m,N) is the trivial cap. The bound is never worse than the native bound with the same retained dimension, and it does not change under any invertible refactorization of the queries and keys.

3Coordinates against the operator

The two sliders separate what changes the coordinates from what changes the head. The shear t moves the queries and keys; the matrix below is recomputed from the moved vectors each time and does not change. The logit scale β changes the head itself: the matrix sharpens, and the intrinsic size grows with it.

Coordinates explorer (static)Four queries and four keys at plus and minus e1 and e2, sheared by S_t with t = 2. The attention matrix is the same for every t. The native query-key radius product is 1 + t squared = 5; the intrinsic size kappa(C) stays 1.Queries and keysquerykeyt = 2softmaxAttention A0.530.070.200.200.070.530.200.200.200.200.530.070.200.200.070.53unchanged in every frameSize of the interaction (log scale)native radius product ΓΓ = 5intrinsic size κ(C)κ(C) = 1110100
Explorer: one operator, many frames. The paper’s rank-two example with β free: the native radius product is β(1+t2) and the intrinsic size is κ(C)=β.The paper’s rank-two example at t=2: the native radius product is 5, the intrinsic size κ(C) is 1, and the attention matrix is the one at t=0.

Here κ(C)=1 exactly at β=1: the original factorization reaches a radius product of 1, and the rectangular contrast C11−C12−C21+C22=4 shows that no factorization can do better.

4Where the exponent comes from

The bound rests on a cover of Gibbs rows: every softmax row is close, in ℓ1, to one of a small set of rows, and the size of that set is a question about convex bodies in fixed dimension. Its exponent is set by the geometry of the support. For self-attention on the unit sphere in ℝd, the worst-case rank grows like (1+β)(d−1)/2 with the temperature; on the full ball, with enough tokens, it grows like βd/2. One radial degree of freedom adds exactly one half to the exponent.

The intrinsic theorem inherits the same exponent, now with the visible dimension r in place of the ambient one, and it cannot be lowered. On a finite grid of Lr points with logits 2sa⊤b−s‖b‖2, the exact approximation rank is Lr while κ(C)=sr(L−1)2/2, so rε(A)=Θ(κ(C)r/2). Low interaction rank alone is not enough either: rank-one families with no residual still need growing rank, so a size term like κ cannot be dropped.

Finite-grid sharpness of the r/2 exponentExact approximation rank L^2 against intrinsic size kappa(C) for the grid family with r = 2 and epsilon = 0.25, L = 3, 5, 8, 12, 20, on log-log axes; the points follow the dashed slope-1 guide.101001000101001000intrinsic size κ(C)exact rank r_ε(A)slope r/2 = 1L = 3: κ(C) = 12.2411, r_ε(A) = 9L = 3L = 5: κ(C) = 48.9643, r_ε(A) = 25L = 5L = 8: κ(C) = 149.9533, r_ε(A) = 64L = 8L = 12: κ(C) = 370.2928, r_ε(A) = 144L = 12L = 20: κ(C) = 1104.7578, r_ε(A) = 400L = 20
Finite-grid sharpness of the r/2 exponentExact approximation rank L^2 against intrinsic size kappa(C) for the grid family with r = 2 and epsilon = 0.25, L = 3, 5, 8, 12, 20, on log-log axes; the points follow the dashed slope-1 guide.101001000101001000intrinsic size κ(C)exact rank r_ε(A)slope r/2 = 1L = 3: κ(C) = 12.2411, r_ε(A) = 9L = 3L = 5: κ(C) = 48.9643, r_ε(A) = 25L = 5L = 8: κ(C) = 149.9533, r_ε(A) = 64L = 8L = 12: κ(C) = 370.2928, r_ε(A) = 144L = 12L = 20: κ(C) = 1104.7578, r_ε(A) = 400L = 20
The r/2 exponent is reached. The grid family at r=2, ε=0.25: exact rank L2 against κ(C), for L=3,5,8,12,20. Exact values from the proposition.

5Measuring the rank directly

The constant Cr in the theorem is not known numerically, so on trained heads the rank is measured rather than predicted. Any explicit matrix of rank k that is checked to lie within ε proves rε(A)≤k; the upper search uses reweighted SVD and alternating-ℓ1 factorization. Lower certificates quantify over every real competitor: Frobenius tails, truncated pseudoinverses, residual-aware bounds and contractive decoders. A failed attempt to approximate is never used as a lower bound. Numerical attention matrices carry a slack χ: Lε+χ(A^)≤rε(A)≤Uε−χ(A^).

Every result therefore says one of three things about a rank budget k: feasible (a checked construction of rank k exists), impossible (a lower certificate exceeds k), or unresolved (neither). Certified here means a float64 computation without interval enclosure.

Part 2 · The evidence

Tested so far tested

  • 9,978 brackets on BERT, GPT-2, Qwen2.5-1.5B and two ViT-Base checkpoints
  • Six practical constructions against the lower certificate on a held-out panel
  • Fixed rank budgets across block widths 64 to 512, on Qwen2.5-1.5B and, blind, on Llama-3.2-1B
  • FLAN-T5 encoder, decoder and cross-attention; BLIP-2 query-to-image cross-attention
  • Native query–key geometry against a calibrated baseline, and under changes of coordinates

In progress In progress

  • Validated numerics: rounding and residual bounds for the stored certificates
  • Brackets for the intrinsic size κ(T) on learned heads

6Results: brackets on learned heads

The primary panel holds 2,302 operators: BERT-base at 256 × 256, the common-prefix blocks of GPT-2 and Qwen2.5-1.5B at 128 × 128, and two ViT-Base checkpoints at 197 × 197, from 64 fixed text sources and 64 ImageNet images. At three tolerances and with secondary decoder cells, the certificate pass produced 9,978 brackets. An independent recomputation found no failure, none of 458 targeted attacks on the lower ends succeeded, and 119 of 119 reconstructed upper witnesses passed the same check.

BERT-base

256 × 256 · 512 operators

42.8%

resolved; median upper ÷ lower 3.82 (lower 11, upper 57)

GPT-2

128 × 128 · 382 operators

99.7%

resolved; median upper ÷ lower 3.00 (lower 2, upper 6)

Qwen2.5-1.5B

128 × 128 · 384 operators

77.1%

resolved; median upper ÷ lower 2.63 (lower 3, upper 9)

ViT-Base (ImageNet-1k)

197 × 197 · 512 operators

59.0%

resolved; median upper ÷ lower 3.41 (lower 6, upper 25)

ViT-Base (ImageNet-21k)

197 × 197 · 512 operators

61.9%

resolved; median upper ÷ lower 3.08 (lower 5, upper 21)

Brackets on 2,302 learned attention operators at ε = 0.25Five panels, one per model family. Each operator is a vertical segment from its certified lower rank to its checked upper rank on a log scale; resolved operators are teal, unresolved grey.141664256BERT-baseresolved 42.8% · U/L 3.82GPT-2resolved 99.7% · U/L 3.00Qwen2.5-1.5Bresolved 77.1% · U/L 2.63ViT-Base (ImageNet-1k)resolved 59.0% · U/L 3.41ViT-Base (ImageNet-21k)resolved 61.9% · U/L 3.08rank (log scale)Each column is one operator, sorted by its upper end; the segment runs from the certified lower rank to the checked upper rank. Dashed: full matrix rank.
Brackets on 2,302 learned attention operators at ε = 0.25One panel per model family; each operator is a segment from its certified lower rank to its checked upper rank (log scale); resolved teal, unresolved grey.141664256BERT-baseresolved 42.8% · U/L 3.82141664256GPT-2resolved 99.7% · U/L 3.00141664256Qwen2.5-1.5Bresolved 77.1% · U/L 2.63141664256ViT-Base (ImageNet-1k)resolved 59.0% · U/L 3.41141664256ViT-Base (ImageNet-21k)resolved 61.9% · U/L 3.08
Every operator of the primary panel at ε = 0.25. Each segment runs from the certified lower rank to the checked upper rank. A bracket counts as resolved when the upper end is at most twice the lower end, or within two ranks or a tenth of full rank of it.

GPT-2, Qwen2.5 and both ViT arms meet the pre-specified identification criterion (resolved share, several disjoint bands, ordering coverage). BERT does not. At a common 128 × 128 shape the encoder and vision operators have larger checked upper budgets than the decoder blocks, and the two ViT checkpoints order difficulty almost identically while fine-tuning on ImageNet-1k raises the upper ends.

7Constructions against certificates

A held-out panel of 1,152 operators (32 text sources and 32 images) asks how close practical low-rank constructions come to the lower certificate. Six fixed implementations were run at ε = 0.25; alternating ℓ1 reached the lowest rank in every cell.

Practical constructions against the lower certificate (ε = 0.25)Left: per family, the rank of the best of six constructions divided by the certified lower rank, with 95% intervals. Right: rank needed by each construction relative to alternating l1. Bottom: where failed attempts fall.Best construction ÷ lower certificate1×2×3×4×BERT-baseBERT-base: 3.73 [3.48, 3.89]; L 11.25, best 57.53.73GPT-2GPT-2: 2.58 [2.50, 3.00]; L 2, best 62.58Qwen2.5-1.5BQwen2.5-1.5B: 2.58 [2.42, 2.78]; L 3, best 9.252.58ViT-Base (ImageNet-1k)ViT-Base (ImageNet-1k): 3.47 [3.00, 3.67]; L 5, best 24.53.47Rank needed ÷ rank for alternating ℓ₁1×2×4×8×Alternating ℓ₁1.00×Row-reweighted SVD1.49×Truncated SVD1.62×Column pivoted-QR1.80×Row pivoted-QR2.09×Representative rows6.45×Where the 435,436 failed construction attempts fallbelow the lower certificate: no matrix of that rank works: 21.1%21.1%between the certificate and the best construction: unresolved: 33.4%33.4%at or above a rank another construction reaches: 45.5%45.5%below the lower certificate: no matrix of that rank worksbetween the certificate and the best construction: unresolvedat or above a rank another construction reaches
Practical constructions against the lower certificate (ε = 0.25)Best construction over certificate per family; each construction relative to alternating l1; where failed attempts fall.Best construction ÷ lower certificate1×2×3×4×BERT-baseBERT-base: 3.73 [3.48, 3.89]; L 11.25, best 57.53.73GPT-2GPT-2: 2.58 [2.50, 3.00]; L 2, best 62.58Qwen2.5-1.5BQwen2.5-1.5B: 2.58 [2.42, 2.78]; L 3, best 9.252.58ViT-Base (ImageNet-1k)ViT-Base (ImageNet-1k): 3.47 [3.00, 3.67]; L 5, best 24.53.47Rank needed ÷ rank for alternating ℓ₁1×2×4×8×Alternating ℓ₁1.00×Row-reweighted SVD1.49×Truncated SVD1.62×Column pivoted-QR1.80×Row pivoted-QR2.09×Representative rows6.45×Where the 435,436 failed construction attempts fallbelow the lower certificate: no matrix of that rank works: 21.1%21.1%between the certificate and the best construction: unresolved: 33.4%33.4%at or above a rank another construction reaches: 45.5%45.5%below the lower certificate: no matrix of that rank worksbetween the certificate and the best construction: unresolvedat or above a rank another construction reaches
The best construction needs 2.6–3.7 times the certified lower rank. Family medians with 95% intervals; truncated SVD needs 1.62 times the rank of alternating ℓ1 under this error. The bottom bar sorts the 435,436 failed attempts into what the certificate can say about them.

The gap is between a construction and a certificate; which end is loose is not identified. The bracket closes exactly in 0%, 17.7%, 19.3%, 0.4% of BERT, GPT-2, Qwen2.5 and ViT-1k operators.

8Rank budgets as the block widens

For Qwen2.5-1.5B the operator is the renormalised N × N common-prefix block of one 1,024-token pass, for 6 fixed heads and 32 arXiv documents, at widths N = 64 to 512. The question is a fixed budget: is rank 8, 16 or 32 enough at ε = 0.25?

Share of operators whose rank budget is shown feasible, by block width (ε = 0.25)Two panels, Qwen2.5-1.5B and Llama-3.2-1B; lines for rank budgets 8, 16 and 32 across block widths 64 to 512.Qwen2.5-1.5B6 heads, 32 documents0%25%50%75%100%64128256512k = 8, N = 64: 63.5% feasiblek = 8, N = 128: 50.0% feasiblek = 8, N = 256: 44.7% feasiblek = 8, N = 512: 29.8% feasiblek = 16, N = 64: 98.4% feasiblek = 16, N = 128: 80.2% feasiblek = 16, N = 256: 49.4% feasiblek = 16, N = 512: 46.6% feasiblek = 32, N = 64: 100.0% feasiblek = 32, N = 128: 100.0% feasiblek = 32, N = 256: 90.4% feasiblek = 32, N = 512: 50.7% feasible51%47%30%Llama-3.2-1B8 heads, 32 fresh documents, blind64128256512k = 8, N = 64: 100.0% feasiblek = 8, N = 128: 97.3% feasiblek = 8, N = 256: 86.7% feasiblek = 8, N = 512: 80.1% feasiblek = 16, N = 64: 100.0% feasiblek = 16, N = 128: 100.0% feasiblek = 16, N = 256: 99.2% feasiblek = 16, N = 512: 87.9% feasiblek = 32, N = 64: 100.0% feasiblek = 32, N = 128: 100.0% feasiblek = 32, N = 256: 100.0% feasiblek = 32, N = 512: 99.6% feasible100%88%80%common-prefix block width N (keys)shown feasiblerank 8rank 16rank 32
Share of operators whose rank budget is shown feasible, by block widthQwen2.5-1.5B above, Llama-3.2-1B below; rank budgets 8, 16 and 32.Qwen2.5-1.5B6 heads, 32 documents0%25%50%75%100%64128256512k = 8, N = 64: 63.5% feasiblek = 8, N = 128: 50.0% feasiblek = 8, N = 256: 44.7% feasiblek = 8, N = 512: 29.8% feasiblek = 16, N = 64: 98.4% feasiblek = 16, N = 128: 80.2% feasiblek = 16, N = 256: 49.4% feasiblek = 16, N = 512: 46.6% feasiblek = 32, N = 64: 100.0% feasiblek = 32, N = 128: 100.0% feasiblek = 32, N = 256: 90.4% feasiblek = 32, N = 512: 50.7% feasible51%47%30%Llama-3.2-1B8 heads, 32 fresh documents, blind0%25%50%75%100%64128256512k = 8, N = 64: 100.0% feasiblek = 8, N = 128: 97.3% feasiblek = 8, N = 256: 86.7% feasiblek = 8, N = 512: 80.1% feasiblek = 16, N = 64: 100.0% feasiblek = 16, N = 128: 100.0% feasiblek = 16, N = 256: 99.2% feasiblek = 16, N = 512: 87.9% feasiblek = 32, N = 64: 100.0% feasiblek = 32, N = 128: 100.0% feasiblek = 32, N = 256: 100.0% feasiblek = 32, N = 512: 99.6% feasible100%88%80%block width N (keys)rank 8rank 16rank 32
Fixed budgets are shown feasible less often as the block widens. Share of operators with a checked construction at the budget. Llama-3.2-1B was run blind, on 32 fresh documents, by agents that had not seen the Qwen results.

On Qwen2.5 the pre-specified class is a decline at all three budgets: rank 16 drops from 98% to 47% of operators and rank 32 from 100% to 51%. The bracket widens with N, the lower ends nearly flat (median 2.5 → 3.5) and the first SVD-type upper ends growing (8.5 → 72), so the operators lose their feasible verdict to unresolved rather than to impossible.

On Llama-3.2-1B the class is a decline at budgets 8 and 16 (feasible share down by 19.9 points and 12.1) and unresolved at 32. Compared with Qwen2.5, budgets 8 and 16 replicate and 32 does not resolve. Two records of that comparison exist and are reported together in the limitations.

9Other attention regimes

The same certificates run unchanged on rectangular and encoder–decoder attention.

FLAN-T5 encoder self-attention

256 × 256

3.81

median upper ÷ lower

FLAN-T5 cross-attention

64 × 256

2.53

median upper ÷ lower

FLAN-T5 decoder self-attention block

32 × 32

1.50

median upper ÷ lower

BLIP-2 query → image

32 × 257

11 of 32

median upper rank; median lower 4.5

10Native coordinates as a predictor

The earlier, native-coordinate version of the geometry (native Γ, SVD-selected interaction dimensions and shape diagnostics) was tested as a predictor of the checked upper budget on held-out sources. Against a weak structural baseline it carries coarse information: held-out error falls by 0.0279 [0.0254, 0.0300] cap units. Against a calibrated baseline of logit, concentration and weight-space features, it adds nothing detectable: the aggregate gain in held-out error is −0.00005 cap units and the relative gain is −0.0035 [−0.0422, +0.0179].

Native-coordinate geometry beyond the calibrated baselineRelative gain of adding the native geometry block to the calibrated baseline, with 95% intervals, per family and aggregated; every interval crosses zero and stays within the dashed ±0.10 criterion.-0.10-0.050+0.05+0.10BERT-baseGPT-2Qwen2.5-1.5BViT-Base (ImageNet-1k)ViT-Base (ImageNet-21k)Equal-weight aggregateshare of the baseline’s remaining error removed by the native block
Native-coordinate geometry beyond the calibrated baselineRelative gain of adding the native geometry block to the calibrated baseline, with 95% intervals, per family and aggregated; every interval crosses zero and stays within the dashed ±0.10 criterion.-0.10-0.050+0.05+0.10BERT-baseGPT-2Qwen2.5-1.5BViT-Base (ImageNet-1k)ViT-Base (ImageNet-21k)Equal-weight aggregateshare of the remaining error removed
No family excludes zero. Relative gain from adding the native block to the calibrated baseline, with 95% source-cluster bootstrap intervals; dashed, the pre-specified ±0.10 criterion.

The same native quantities move when the coordinates move. In a gauge audit of 144 operators, all 576 operator–gauge records pass the measurement’s own invariance check, while native Γ changes by more than 10% in every operator under the diagonal and dense gauges, by up to 21×.

11Limitations and what didn’t work

Scope

  • The panels are fixed checkpoints with fixed head rosters: one BERT encoder, GPT-2, Qwen2.5-1.5B, two ViT-Base checkpoints, FLAN-T5-base, Llama-3.2-1B and the BLIP-2 Q-Former. The GPT-2 and Qwen2.5 panels use six of eight pre-specified non-local heads. Nothing is claimed for other models, larger language models or long context.
  • The theorem assumes one common positive support. Causal masks need more theory, so decoder results come from common-prefix blocks, and the native-causal brackets are measurements, not uses of the theorem.
  • “Block width” is the width of a common-prefix block, not a context length. The query tokens change with N, so there is no causal length intervention, and each model was run on one dataset.
  • The results concern representation rank. No runtime, memory, sparsity, accuracy or hardness conclusion follows.

What the evidence does not show

  • Exact optimum ranks. Where the two ends differ the rank is an interval, and no midpoint is reported. The formal identification gate of the bracketing study was not met: BERT resolves 219 of 512 operators (42.8%), with layer 0 at 0/128.
  • How much rank a trained head needs, read from its geometry. The computable intrinsic score agrees with certified rank order in BERT (0.996) and ViT-1k (0.995) but not in GPT-2 (0.159) or Qwen2.5 (0.391); FLAN-T5 and BLIP-2 align, though in FLAN-T5 the plain logit scale does better in the encoder (−0.068) and decoder (−0.076). Added to the calibrated baseline, the intrinsic block adds nothing (gain −0.000116 cap units) and makes leave-one-head-out prediction 47.2% worse. That test was retrospective: its targets and baseline were known before its protocol was frozen.
  • A numerical value for the theorem. Cr is not estimated; with the explicit constant the bound returns the trivial cap at the tested dimensions for r≥2. The experiments measure the rank instead, and their size proxy is the balanced-factor surrogate κbal, which can exceed κ by an unbounded factor. The computable specialisations (balanced-factor size, canonical logit-SVD family) passed one internal proof check, the second with conditions on how they are computed; they do not solve the optimisation that defines κ.
  • Impossibility at budget. Rank-k impossibility certificates are rare: on Qwen2.5 they reach 6.9% of operators at rank 8 and N = 512 and never appear at 16 or 32, and Llama has none. At N = 512 63.3%, 53.4%, 49.3% of Qwen2.5 operators stay unresolved at ranks 8, 16 and 32. “Degrades” is a statement about these certificates and constructions, not about the true rank.
  • The native null as an equivalence. The calibrated baseline’s error was already below the pre-set resolution in several families, so the test rules out a large relative gain, not every effect.
  • Claims of the NeurReps abstract that the newer work dropped. The per-instance law there used the native Γ=βρQρK, which the intrinsic κ(T) replaces. Its 84-head BERT calibration (a smaller effective dimension in 41.9% of head–temperature cells; Spearman 0.574 and 0.606 with constructive upper certificates) is no longer offered as support for the geometry, and the reading that effective geometry tracks low-rank demand did not survive the held-out test above.

Implementation and reproducibility

  • No validated numerics yet. Certificates are float64 without interval enclosure: computed singular values carry no accuracy bounds, the pseudoinverse witness’s residual is not charged, the common-prefix slack omits the rounding of renormalisation, and upper checks have no evaluation-error bounds. The validated re-check and the κ(T) brackets were started and are not finished; no result from them exists.
  • Replication records. The blind Llama comparison has a pre-registered record in which every label is a defect, caused by a line-ending fault in the comparing code, and a corrected post-hoc reading (replicated at 8 and 16, unresolved at 32) that changes no statistic. Its experiment audit failed on closing-record wording while reproducing every statistic. An earlier Llama replication was not blind: its outcomes were seen before its design was pinned.
  • Verification. Independent recomputations ran in isolated contexts from the same model vendor; there is no cross-vendor replication, and some verifiers covered a sample of operators rather than all of them. The intrinsic finite-grid sharpness result has one proof pass and one independent re-check; the native form of the r/2 sharpness was verified independently.
  • Release. Code and per-operator records are to be released on publication and are not public yet.

12Where this is going proposed

  1. Validated certificates

    A post-processing pass with rounding-error and residual bounds over the stored witnesses and constructions, so impossible, unresolved and feasible hold without the float64 qualifier.

  2. κ(T) on learned heads

    Primal and dual brackets for the intrinsic size on a representative subset, to test the theorem’s own quantity rather than its balanced-factor surrogate.

  3. Beyond one common support

    Theory for causal masks; and the open questions of the paper: equality of the intrinsic and native optima, joint sharpness in κ, residual and ε, and an efficient global algorithm for the profile.

Part 3 · Technical details

13Definitions and the theorem

For an attention matrix A∈ℝm×N and 0<ε<1,

rε(A)=min{rankB:maxi∑j|Aij−Bij|≤ε}, ‖M‖row,1=supmaxj‖vj‖E≤1maxi‖∑jMijvj‖E.

On one common positive support, with Pm=I−1m11⊤, the visible interaction and the projective residual of a centred T are

C=Pm(logA)PN, q(E)=infbmaxiosc(Ei·−b⊤), τ(T)=2tanh(q(C−T)4).

The residual rests on the sharp inequality, with σ the softmax map,

‖σ(x+h)−σ(x)‖1≤2tanh(osc(h)/4).

The size is a minimum over factorizations, attained in rank dimension:

κ(T)=infT=XY⊤rad(X)rad(Y), rad(X)=infamaxi‖Xi·−a‖2.

The proof replaces A by the softmax of T plus a common column term, which costs at most τ(T) in row-ℓ1; translating an optimal factorization to the centres of its enclosing balls turns that matrix into a rank-r weighted Gibbs family of radius product κ(T); the weighted Gibbs-row cover at budget ε−τ(T) then gives the bound. For comparison, the native bound uses the radii of one factorization, Γ=βρQρK, in the same place; at full interaction, minimising Γ over invertible gauges recovers κ(C).

14Support-geometry laws

For self-attention Aβ(X)ij∝eβ⟨xi,xj⟩ with n tokens, write ℜ for the worst case of rε over configurations. For fixed d≥2 and ε, on the sphere, and on the ball when β≥β0(d,ε) and n≥Cdeβ/8:

ℜεsph(n,d,β)≍min{n,(1+β)(d−1)/2}, ℜε(n,d,β)≍βd/2.

The lower bounds come from decoding constructions (well-separated sphere centres; a discrete Gaussian on a lattice in the ball) whose attention is close to a block projection, so any good approximant must be invertible on the messages. The ball law needs exponentially many tokens, so it is a worst-case capacity statement, not a prediction for ordinary context lengths. The finite-grid family of section 4, Hab=2sa⊤b−s‖b‖2 on {0,…,L−1}r with s=log8r1−ε, has rankC=r, κ(C)=sr(L−1)2/2 and rε(A)=Lr.

15Full tables

Brackets on the primary panel at ε = 0.25
FamilyShapeOperatorsResolvedMedian LMedian UMedian U/L
BERT-base256 × 25651242.8%11573.82
GPT-2128 × 12838299.7%263.00
Qwen2.5-1.5B128 × 12838477.1%392.63
ViT-Base (ImageNet-1k)197 × 19751259.0%6253.41
ViT-Base (ImageNet-21k)197 × 19751261.9%5213.08
Held-out construction benchmark at ε = 0.25
FamilyMedian lower certificateMedian best constructionRatio (95% interval)
BERT-base11.2557.53.73 [3.48, 3.89]
GPT-2262.58 [2.50, 3.00]
Qwen2.5-1.5B39.252.58 [2.42, 2.78]
ViT-Base (ImageNet-1k)524.53.47 [3.00, 3.67]
Share of operators shown feasible at fixed rank budgets, ε = 0.25
Model · budgetN = 64128256512Class
Qwen2.5-1.5B · rank 863.5%50.0%44.7%29.8%declines with width
Qwen2.5-1.5B · rank 1698.4%80.2%49.4%46.6%declines with width
Qwen2.5-1.5B · rank 32100.0%100.0%90.4%50.7%declines with width
Llama-3.2-1B (blind) · rank 8100.0%97.3%86.7%80.1%declines with width
Llama-3.2-1B (blind) · rank 16100.0%100.0%99.2%87.9%declines with width
Llama-3.2-1B (blind) · rank 32100.0%100.0%100.0%99.6%trend unresolved

16Reproducibility notes

Certificate pass
3,326 operators at tolerances 0.15, 0.25 and 0.35: 9,978 brackets, every upper end an explicit matrix checked in maximum-row-ℓ1; independent recomputation with no failure; 0 of 458 targeted lower-end attacks succeeded; 119/119 reconstructed upper witnesses pass.
Numerics
Certificates and constructions in float64 without interval enclosure, with the slack χ of each numerical attention matrix charged to both ends.
Protocols
Each study fixed its panel, tolerances and decision rules in a frozen protocol before its run of record, and a separate verifier rechecked its numbers against the stored per-operator records.
Exact checks
The gauge example and the grid family are computed exactly from their definitions; the small exact oracle closes six rational cells, one with unrestricted rank one and row-stochastic rank two.

Summary

17In brief

Design

The complexity of a head is the rank its normalised operator needs to keep every bounded output within ε. The quantities that control it are taken from the doubly centred log-attention and minimised over factorizations, so no query–key frame is preferred.

Results

The bound carries the dimension exponent of the support geometry, and that exponent is attained. On trained encoders, decoders and vision transformers, the target can be bracketed directly, and practical constructions sit a factor of a few above what the certificates exclude.

Open questions

Whether the theorem’s own size, rather than a surrogate, tracks trained-head difficulty; how to close the wide brackets of BERT and of large blocks; and whether the verdicts survive validated arithmetic. Limitations →

18Sources

Paper: Intrinsic Interaction Geometry Controls the Low-Rank Complexity of Softmax Attention, Yuhe Sui, Jianing Zhang, Yingzhi Tang, arXiv 2608.28150 (v1 28 August 2026, v2 22 September 2026). [arXiv]

Extended abstract: The Approximation Rank of Softmax Attention: Sharp Geometric Laws and Robust Interaction Dimension, Yuhe Sui, Jianing Zhang, NeurReps Workshop (Symmetry and Geometry in Neural Representations), NeurIPS 2026, accepted. [OpenReview]

Accepted Research

  • The Approximation Rank of Softmax Attention: Sharp Geometric Laws and Robust Interaction Dimension

    Yuhe Sui, Jianing Zhang

    NeurReps Workshop, NeurIPS 2026 (extended abstract) · 2026 · Accepted [OpenReview]

Numbers on this page come from the project’s study records: the bracketing study, the held-out construction benchmark, the block-width studies on Qwen2.5 and Llama, the FLAN-T5 and BLIP-2 studies, and the native and intrinsic predictor tests. Where a newer study or the current text of the paper supersedes the extended abstract, this page follows the newer one.

Back to top ↑