---
title: Low-Rank Stabilizer Decomposition
url: https://www.emergentmind.com/topics/low-rank-stabilizer-decomposition
type: topic
---

# Low-Rank Stabilizer Decomposition

Low-rank stabilizer decomposition denotes the representation of a quantum state as a superposition of a small number of stabilizer states. For an $n$-qubit state $|\psi\rangle$, the associated minimal number of terms is the stabilizer rank $\chi(\psi)$, and the corresponding approximate notion is the $\delta$-approximate stabilizer rank $\chi_\delta(\psi)$. In the literature on classical simulation of Clifford-dominated quantum circuits with non-Clifford resources, these quantities are central because the running time of several simulation methods is determined by the stabilizer rank of tensor powers of magic states; consequently, the subject has developed along two tightly linked directions: constructive low-rank decompositions and lower bounds showing when such decompositions must fail [1808.00128; 2106.03214].

## 1. Formal framework

The exact stabilizer rank of a pure state is defined as the smallest integer $r$ such that
\[
|\psi\rangle=\sum_{j=1}^r c_j|\varphi_j\rangle,
\]
where each $c_j\in\mathbb{C}$ and each $|\varphi_j\rangle$ is a stabilizer state. A stabilizer state is a state of the form $U|0^n\rangle$ for a Clifford unitary $U$ [2106.03214]. The approximate version, $\chi_\delta(\psi)$, is the minimal stabilizer rank among all states $|\varphi\rangle$ such that $\|\psi-\varphi\|_2\leq \delta$ [2106.03214].

A closely related quantity is the stabilizer extent. For a normalized state $|\psi\rangle$,
\[
\xi(\psi)=\min \|c\|_1^2,
\]
where the minimum is taken over all stabilizer decompositions $|\psi\rangle=\sum_\alpha c_\alpha |\phi_\alpha\rangle$ [1808.00128]. The stabilizer fidelity
\[
F(\psi)=\max_{|\phi\rangle\ \text{stabilizer}} |\langle \phi|\psi\rangle|^2
\]
measures maximal overlap with a stabilizer state [1808.00128].

These notions separate several different questions. Exact stabilizer rank concerns exact synthesis in the stabilizer basis. Approximate stabilizer rank concerns simulation and approximation under norm error. Stabilizer extent is tailored to sparsification and simulation cost bounds. Stabilizer fidelity connects decomposition complexity to proximity to a single stabilizer state. A recurrent theme in the literature is that these measures are related but not interchangeable.

## 2. Role in classical simulation of quantum circuits

The algorithmic importance of low-rank stabilizer decomposition was systematized in the simulation framework for circuits containing Clifford gates and non-Clifford diagonal gates [1808.00128]. One route is a sum-over-Cliffords method: if a circuit is written as
\[
U=D_m V_m D_{m-1}V_{m-1}\dots D_1V_1D_0,
\]
with $D_j$ Clifford and $V_j$ diagonal, then each $V_j$ can be decomposed as
\[
V_j=\sum_\alpha c^{(j)}_\alpha K^{(j)}_\alpha
\]
over Clifford operators. The full circuit is then expressed as a superposition of Clifford circuits, and the output state becomes a superposition of stabilizer states [1808.00128].

A second route is gadget-based simulation. In Clifford+$T$ and Clifford+CCZ settings, non-Clifford gates are replaced by state-injection gadgets using magic states, and the magic-resource state is itself decomposed into stabilizer states. Two variants were described: a fixed-sample approach that requires exact decompositions, and a random-sample approach that can use approximate decompositions [1808.00128].

The resulting complexity bounds are expressed in terms of stabilizer extent. For Clifford+diagonal circuits, the strong-simulation cost scales as
\[
\tilde O\left(\prod_j \xi(V_j)\right),
\]
and the weak-simulation cost as
\[
\tilde O\left(\delta^{-2}\prod_j \xi(V_j)\right)
\]
[1808.00128]. This framework extended earlier Clifford+$T$ methods to arbitrary diagonal gates and yielded simulations of 40–50 qubits with over 60 non-Clifford gates without high-performance computers. The same work reported a QAOA simulation processing superpositions of $\chi\sim 10^6$ stabilizer states and sampling from the full $n$-bit output distribution, and Hidden Shift simulations with circuits including up to 64 $T$ gates or 16 CCZ gates [1808.00128].

The significance of these results is operational rather than merely representational: low-rank stabilizer decompositions are not just compact descriptions of states, but algorithmic surrogates for the non-Clifford content of a circuit.

## 3. Approximate decompositions, sparsification, and extent

Approximate low-rank stabilizer decomposition is governed by a sparsification principle. If
\[
|\psi\rangle=\sum_j c_j |\phi_j\rangle
\]
with normalized stabilizer states $|\phi_j\rangle$, then for integer $k$ there exists a state $|\Omega\rangle$ that is a sum of $k$ stabilizer states sampled according to $|c_j|$ such that
\[
\mathbb{E}\| |\psi\rangle-|\Omega\rangle\|^2=\frac{\|c\|_1^2}{k}.
\]
This yields the bound
\[
\chi_\delta(\psi)\le 1+\frac{\|c\|_1^2}{\delta^2},
\]
and hence
\[
\chi_\delta(\psi)\le 1+\frac{\xi(\psi)}{\delta^2}
\]
[1808.00128].

For Clifford magic states, the extent is especially tractable: $\xi(\psi)=1/F(\psi)$ [1808.00128]. For the single-qubit $T$ state, since $F(|T\rangle)=\cos^2(\pi/8)$, one obtains
\[
\chi_\delta(T^{\otimes m})\le O\left(\delta^{-2}\cos^{-2m}(\pi/8)\right)
\]
[1808.00128]. The same work also established multiplicativity of stabilizer extent for tensor products of single- or few-qubit states, up to three qubits per factor:
\[
\xi\left(\bigotimes_j |\psi_j\rangle\right)=\prod_j \xi(|\psi_j\rangle)
\]
[1808.00128].

These results clarify why approximate decompositions became central to simulation practice. Exact rank can be too rigid for algorithm design, whereas extent and approximate rank admit probabilistic sparsification with explicit error control. At the same time, the literature later showed that approximate rank can behave very differently from exact rank, so any interpretation of “low rank” must specify which notion is meant.

## 4. Lower bounds and limitations

The lower-bound theory developed around tensor powers of single-qubit magic states. For $|B\rangle$ equal to $|H\rangle$ or $|R\rangle$, a 2021 result proved
\[
\chi(|B^{\otimes n}\rangle)=\Omega(n),
\]
improving the previous $\Omega(\sqrt n)$ lower bound, and also established that for a sufficiently small constant $\delta$,
\[
\chi_\delta(|B^{\otimes n}\rangle)=\Omega(\sqrt n/\log n)
\]
[2106.03214]. The same work represented stabilizer states as quadratic functions over affine subspaces of $\mathbb{F}_2^n$:
\[
|\varphi\rangle=\sum_{x\in A} i^{\ell(x)}(-1)^{q(x)}|x\rangle,
\]
and recast the amplitudes of $|H\rangle^{\otimes n}$ through the function
\[
H_n(x)=\cos(\pi/8)^{n-|x|}\sin(\pi/8)^{|x|}.
\]
The proofs combined directional derivatives of quadratic polynomials, combinatorial arguments using Kleitman’s theorem, and complexity-theoretic tools including Razborov–Smolensky low-degree approximations and correlation bounds against majority [2106.03214].

A complementary 2021 paper introduced number-theoretic and algebraic-geometric techniques. It showed that if the computational-basis coordinates of a state contain an exponentially increasing subsequence of length $p$, then
\[
\chi(\psi)\ge \frac{p}{4\log_2 p},
\]
and in particular
\[
\chi(T^{\otimes n})\ge \frac{n+1}{4\log_2(n+1)}.
\]
More generally, for any non-stabilizer qubit state $\psi$,
\[
\chi(\psi^{\otimes n})=\Omega(n/\log n),
\]
and for some constant $\delta>0$,
\[
\chi_\delta(\psi^{\otimes n})=\Omega(\sqrt n/\log n)
\]
[2110.07781].

A later probabilistic approach strengthened approximate lower bounds further. It stated that for $|T\rangle^{\otimes n}$ and a wide range of approximation parameters, the approximate stabilizer rank satisfies a $\tilde \Omega(n^2)$ lower bound, using a combination of Haar-random-state arguments, a magic-state teleportation analysis, and a result about trading Clifford operations with $T$ gates [2305.10277].

| Family or setting | Lower bound | Source |
|---|---|---|
| $|H\rangle^{\otimes n}$ or $|R\rangle^{\otimes n}$ | $\chi=\Omega(n)$ | [2106.03214] |
| $|H\rangle^{\otimes n}$ or $|R\rangle^{\otimes n}$, constant-error approximation | $\chi_\delta=\Omega(\sqrt n/\log n)$ | [2106.03214] |
| Any non-stabilizer qubit $\psi^{\otimes n}$ | $\chi=\Omega(n/\log n)$ | [2110.07781] |
| Any non-stabilizer qubit $\psi^{\otimes n}$, some constant $\delta>0$ | $\chi_\delta=\Omega(\sqrt n/\log n)$ | [2110.07781] |
| $|T\rangle^{\otimes n}$, wide range of approximation parameters | $\tilde \Omega(n^2)$ for approximate rank | [2305.10277] |

The conceptual consequence is that low-rank stabilizer decomposition is a powerful simulation primitive only up to a point. The lower-bound papers explicitly connect these limitations to classical simulation hardness and to representation problems for sums of quadratic phases over $\mathbb{F}_2$ [2106.03214; 2305.10277].

## 5. Structural phenomena: separations, multiplicativity, and genericity

Beyond asymptotic lower bounds, the theory uncovered several structural effects. One paper refined a theorem of Moulton to exhibit an explicit sequence of product states with exponential stabilizer rank but constant approximate stabilizer rank, and then gave a stronger field-theoretic construction in which an $n$-qubit product state has exact stabilizer rank $2^n$ while $\chi_\delta(\psi_n^\theta)=1$ for sufficiently large $|\theta|$ [2110.07781]. This established an extreme separation between exact and approximate notions.

The same work gave the first non-trivial examples of multiplicative stabilizer rank under tensor product. For the explicitly constructed two-qubit family
\[
\psi_\alpha=\frac{1}{\sqrt{1+2|\alpha|^2}}\left(e_{00}+\alpha(e_{01}+e_{10})\right),
\]
one has, for all but finitely many $\alpha$,
\[
\chi(\psi_\alpha)=2,\qquad \chi(\psi_\alpha\otimes \psi_\alpha)=4,
\]
and in particular this holds for transcendental $\alpha$ [2110.07781].

That paper also introduced the generic stabilizer rank
\[
\chi_n=\max_{\psi\in S(\mathbb{C}^2)} \chi(\psi^{\otimes n}),
\]
together with its real analogue. Using algebraic geometry, it showed that for each $n$, all but finitely many qubit states modulo phase saturate $\chi_n$, because the set of $n$th tensor powers forms a 1-dimensional irreducible projective variety whose intersection with low-degree secant varieties is finite unless it is the whole variety [2110.07781]. The same framework improved an upper bound on $\chi_n$ from $O((n+1)2^{n/2})$ to $O(2^{n/2})$.

These results changed the interpretation of low-rank stabilizer decomposition. They imply that rank behavior is not exhausted by standard magic states; product states can already display maximal exact rank, constant approximate rank, tensor-product multiplicativity, and generic maximality.

## 6. Orbit-dependent decompositions and higher-dimensional extensions

Recent work generalized the subject from qubit magic states to Clifford orbits of magic states in qubit and qutrit systems. For qubits, the asymptotic per-copy exponent is
\[
\gamma_M=\limsup_{m\to\infty}\frac{\log_p \chi(M^{\otimes m})}{m},
\]
with $p=2$ for qubits and $p=3$ for qutrits [2605.28586]. Lower $\gamma_M$ corresponds to lower asymptotic stabilizer rank growth.

In the qubit case, an explicit decomposition of the qubit $T$-type orbit at four copies gives
\[
\chi(T^{\otimes 4})=3,
\qquad
\gamma_T\le \frac{\log_2(3)}{4}\approx 0.396,
\]
matching the existing best-known exponent but through a direct algebraic identity rather than an entangled cat-state construction [2605.28586]. The same paper emphasized orbit dependence at small tensor powers: at $m=4$, $\chi(H^{\otimes 4})=4$ whereas $\chi(T^{\otimes 4})=3$.

For qutrits, four inequivalent single-qutrit Clifford orbits were analyzed: Strange, Norrell, Hadamard-eigenstate, and qutrit $T$-state. Explicit decompositions yielded
\[
\gamma_S\le \log_3(2)/2\approx 0.316,
\qquad
\gamma_{H_3},\gamma_N\le \log_3(4)/3\approx 0.421,
\]
all strictly below the prior $\gamma_{T_3}\le 1/2$ baseline [2605.28586]. The same work also proved the first nontrivial $\Omega(m/\log m)$ asymptotic lower bounds for the Hadamard-eigenstate and Norrell orbits, and described two-qutrit Clifford circuits that convert two copies of these states into an injectable phase state with constant success probability. By contrast, for the Strange orbit, exhaustive search found no two-copy or three-copy Clifford protocol generating a non-Clifford diagonal gate [2605.28586].

This orbit-by-orbit analysis shows that low-rank stabilizer decomposition is sensitive not only to the amount of magic resource but also to its Clifford orbit. The work formalized this through explicit decompositions, asymptotic exponents, and machine-checked verification in the open-source library *stabrank* with Lean 4 proof formalizations [2605.28586].

## 7. Testing low stabilizer complexity and relation to fidelity

The most recent direction connects low-rank stabilizer decomposition to quantum property testing and pseudorandomness. A 2024 work improved tolerant testing of stabilizer states by showing that one can efficiently distinguish whether the maximum fidelity of a quantum state with a stabilizer state is at least $\epsilon_1$ or at most $\epsilon_2$, provided $\epsilon_2\le \epsilon_1^{O(1)}$ [2410.24202]. The proof combines a random Clifford map with tools from higher-order Fourier analysis.

The same paper studied low stabilizer rank states directly. It showed that if, for an infinite family of quantum states, the stabilizer rank is lower than a constant independent of system size, then stabilizer fidelity is lower bounded by an absolute constant [2410.24202]. Using a result of [GIKL22], the paper drew the implication that low approximate stabilizer rank states are not pseudo-random. In a more quantitative formulation, for every constant $k$ there exists $d_k>0$, independent of system size, such that $\chi(|\psi\rangle)\le k$ implies $F(|\psi\rangle)>d_k$ [2410.24202].

This line of work places low-rank stabilizer decomposition within a broader complexity landscape. Low stabilizer rank is not only a simulation resource measure; it also constrains how “random-looking” a state can be, links to Gowers-norm methods, and supports efficient testing primitives. A plausible implication is that future progress on low-rank stabilizer decomposition will continue to move simultaneously in three directions: sharper constructive decompositions, sharper lower bounds, and sharper certification procedures for identifying states of low stabilizer complexity.

Source: https://www.emergentmind.com/topics/low-rank-stabilizer-decomposition