---
title: Optimal Stabilizer Ground State
url: https://www.emergentmind.com/topics/optimal-stabilizer-ground-state
type: topic
---

# Optimal Stabilizer Ground State

An optimal stabilizer ground state is a stabilizer state selected relative to a Hamiltonian expressed in Pauli operators. In the basic variational sense, it is the stabilizer state minimizing $\langle \psi|H|\psi\rangle$ over all stabilizer states; in the refined convention introduced for degenerate stabilizer-ground-state manifolds, it is the member of that manifold with the highest fidelity with the true ground state. The notion turns the stabilizer formalism from a language for Clifford circuits and error correction into a ground-state ansatz with exact energy certificates for Pauli Hamiltonians, an exact linear-scaled algorithm for 1D local sparse systems, periodic extensions, and several roles in quantum simulation, state preparation, and imaginary-time refinement [2403.08441][2603.06286].

## 1. Definition and formal setting

For $n$ qubits, let the set of Hermitian Pauli operators be $\mathcal{P}_n=\pm\{I,X,Y,Z\}^{\otimes n}$. A stabilizer group $S\subset \mathcal{P}_n$ is a commuting subgroup with $-I\notin S$. If $g_1,\ldots,g_m$ are independent commuting generators, one writes $S=\langle g_1,\ldots,g_m\rangle$. A full stabilizer group has $n$ independent generators, and a stabilizer state $|\psi_S\rangle$ is the unique pure state satisfying $g_i|\psi_S\rangle=|\psi_S\rangle$ for all generators. Equivalently,
$$
\mathrm{Stab}(|\psi\rangle)=\{P\in \mathcal{P}_n\mid P|\psi\rangle=|\psi\rangle\},
$$
and a state is a stabilizer state iff $\mathrm{Stab}(|\psi\rangle)$ is a full stabilizer group [2403.08441].

For a Pauli Hamiltonian
$$
H=\sum_{P\in\mathfrak{P}} w_P P,
$$
the stabilizer ground state is defined as the stabilizer state minimizing the energy expectation. In this usage, “optimal” means the energy-minimizing stabilizer state for $H$. If the true ground state of $H$ happens to be a stabilizer state, then the optimal stabilizer ground state coincides with the true ground state; otherwise it is the best stabilizer approximation within the stabilizer ansatz [2403.08441].

A second convention arises when the minimizer is not unique. Over the set $\mathcal{S}_n$ of all $n$-qubit stabilizer states, one defines
$$
E_{\mathrm{SGS}}:=\min_{|\phi\rangle\in\mathcal{S}_n}\langle \phi|H|\phi\rangle,
$$
and the stabilizer-ground-state manifold
$$
\mathcal{G}_{\mathrm{SGS}}:=\{|\phi\rangle\in \mathcal{S}_n:\langle \phi|H|\phi\rangle=E_{\mathrm{SGS}}\}.
$$
The optimal stabilizer ground state is then
$$
|\phi^\ast\rangle:=\arg\max_{|\phi\rangle\in\mathcal{G}_{\mathrm{SGS}}} |\langle \psi_0|\phi\rangle|^2,
$$
where $|\psi_0\rangle$ is the true ground state. In this convention, “optimal” resolves degeneracy by fidelity with the exact ground state rather than by energy alone [2603.06286].

This terminological split is substantial rather than merely stylistic. It separates exact minimization within the stabilizer ansatz from the additional problem of choosing, among energy-degenerate stabilizer minimizers, the one most useful as a proxy for the physical ground state.

## 2. Energy evaluation and algebraic characterization

The basic energy rule is unusually simple. For any stabilizer state $|\psi\rangle$ and Pauli $P\in\mathcal{P}_n$,
- if $P\notin \pm \mathrm{Stab}(|\psi\rangle)$, then $\langle \psi|P|\psi\rangle=0$;
- if $P\in \mathrm{Stab}(|\psi\rangle)$, then $\langle \psi|P|\psi\rangle=+1$;
- if $-P\in \mathrm{Stab}(|\psi\rangle)$, then $\langle \psi|P|\psi\rangle=-1$.

Hence, for $H=\sum_{P\in\mathfrak{P}} w_P P$, only Hamiltonian terms contained in $\pm\mathrm{Stab}(|\psi\rangle)$ contribute to the energy; all others contribute $0$ [2403.08441].

This leads to the stabilizer-group energy functional
$$
E_{\mathrm{stab}}(P,S)=
\begin{cases}
+1,& P\in S,\\
-1,& -P\in S,\\
0,& \text{otherwise},
\end{cases}
$$
and
$$
E_{\mathrm{stab}}(H,S)=\sum_{P\in\mathfrak{P}} w_P E_{\mathrm{stab}}(P,S).
$$
If $S=\langle \mathrm{Stab}(|\psi\rangle)\cap (\pm\mathfrak{P})\rangle$, then $\langle \psi|H|\psi\rangle=E_{\mathrm{stab}}(H,S)$ [2403.08441].

The central structural result is formulated in terms of closed commuting subsets. Let $\pm\mathfrak{P}=\{\pm P\mid P\in\mathfrak{P}\}$ and define
$$
\mathcal{S}(\mathfrak{P})=\{Q\subseteq \pm \mathfrak{P}\mid Q=\langle Q\rangle\cap \pm\mathfrak{P},\ -I\notin \langle Q\rangle\}.
$$
Each $Q\in\mathcal{S}(\mathfrak{P})$ generates a valid stabilizer group built only from Hamiltonian-supported Pauli terms. The theorem states
$$
E_{\mathrm{gs}}^{\mathrm{stab}}=\min_{Q\in \mathcal{S}(\mathfrak{P})} E_{\mathrm{stab}}(H,\langle Q\rangle),
$$
and the minimizing $Q$ is the closed maximally-commuting subset (CMCS). Any stabilizer state stabilized by $S=\langle Q\rangle$ is a stabilizer ground state [2403.08441].

The CMCS formalism replaces brute-force search over all stabilizer states, whose count is of order $2^{(n+1)(n+2)/2}$, by a search over subsets induced by the sparse Pauli support of the Hamiltonian. If $\mathfrak{P}$ acts on at most $n$ qubits, then
$$
|\mathcal{S}(\mathfrak{P})|\le (n+1)|\pm\mathfrak{P}|^n.
$$
For sparse Pauli Hamiltonians with $|\mathfrak{P}|=O(\mathrm{poly}(n))$, this scales as $\exp(C n\log n)$ rather than with the full stabilizer-state count [2403.08441].

A noteworthy methodological point is that this formalism does not rely on the binary symplectic representation. By contrast, other large-scale OSGS procedures do use tableau and symplectic linear algebra to test membership, commutation, and generator independence in polynomial time [2403.08441][2603.06286].

## 3. Exact algorithms in one dimension and periodic extensions

For 1D $k$-local sparse Pauli Hamiltonians, there is an exact algorithm that computes the stabilizer ground state in linear time in the chain length $n$, with a sub-exponential dependence on locality $k$. The starting point is a decomposition of $\pm\mathfrak{P}$ by the last non-identity site:
$$
\pm\mathfrak{P}_m=\{P\in \pm\mathfrak{P}\mid \mathrm{last}(P)=m\}.
$$
For any closed commuting subset $Q$, one writes $Q=\bigcup_m Q_m$ with $Q_m=Q\cap \pm\mathfrak{P}_m$, and the stabilizer energy becomes
$$
E_{\mathrm{stab}}(H,\langle Q\rangle)=\sum_m h_m(Q_m),
$$
where
$$
h_m(Q_m):=\sum_{Q\in Q_m\cap \mathfrak{P}_m} w_Q-\sum_{Q\in Q_m\cap (-\mathfrak{P}_m)} w_{-Q}.
$$
This localizes energy evaluation by last-site index [2403.08441].

The dynamic-programming construction uses a state machine
$$
A_m(Q)=\big(S_{\mathrm{proj}}^m(Q_{\le m}),\ \pm\mathfrak{P}_{\mathrm{invalid}}^m(Q_{\le m}),\ S_{\mathrm{right}}^m(Q_{\ge m})\big),
$$
with
$$
\pm\mathfrak{P}_{\mathrm{invalid}}^m(Q_{\le m})=\{P\in \pm\mathfrak{P}_{>m}\mid [P,Q_{\le m}]\neq 0\},
$$
$$
S_{\mathrm{proj}}^m(Q_{\le m})=\mathfrak{P}_{>m-k}(\langle Q_{\le m}\rangle),
\qquad
S_{\mathrm{right}}^m(Q_{\ge m})=\mathfrak{P}_{\le m}(\langle Q_{\ge m}\rangle),
$$
and
$$
G(A_m)=S_{\mathrm{right}}^m\cap \pm\mathfrak{P}_m=Q_m.
$$
The crucial decoupling property is
$$
\{Q\mid A_m\}=\{Q_{\le m}\mid A_m\}\otimes \{Q_{>m}\mid A_m\},
$$
which turns the search into a Markov-like dynamic program [2403.08441].

The conditioned recurrence is
$$
E_{\mathrm{gs}}^{\le m}(A_m)=h_m(G(A_m))+\min_{A_{m-1}\in F_{\mathrm{back}}(A_m)} E_{\mathrm{gs}}^{\le m-1}(A_{m-1}),
$$
with initial condition $E_{\mathrm{gs}}^{\le 1}(A_1)=h_1(G(A_1))$. After iterating from $m=1$ to $n$,
$$
E_{\mathrm{gs}}^{\mathrm{stab}}=\min_{A_n\in \tilde{\mathcal{A}}_n} E_{\mathrm{gs}}^{\le n}(A_n),
$$
and the minimizing CMCS is reconstructed by backtracking [2403.08441].

The candidate-state bound is
$$
|\tilde{\mathcal{A}}_m'|<N_A=(4kM)^{3k}=\exp\!\big(3k\log(4kM)\big),
$$
where $|\mathfrak{P}_m|\le M$ for all $m$. For local sparse Hamiltonians with $|\mathfrak{P}_m|=O(\mathrm{poly}(k))$, the overall runtime is
$$
T=O\!\big(n\exp(Ck\log k)\big),
$$
linear in system size and sub-exponential in locality [2403.08441].

The same framework extends to infinite 1D periodic local sparse Hamiltonians of period $l$ by imposing periodic boundary conditions on a supercell of size $c\,l$. The associated state machines satisfy $A_{m+cl}\simeq A_m$ for some $c\le N_A$, and the energy density is obtained from a periodic dynamic program with $A_{cl}\simeq A_0$. In higher-dimensional periodic settings, the analogous search is restricted to periodic closed commuting subsets $\mathcal{S}_c(\mathfrak{P})$, whose minimizer is the closed maximally-commuting periodic subset (CMCPS) [2403.08441].

Because the algorithm returns the explicit generator set $S^\star=\langle Q^\star\rangle$, the output is not only an energy value but also a certificate of optimality within the stabilizer ansatz. If $S^\star$ has $n-s$ independent generators, the degeneracy is $2^s$ [2403.08441].

## 4. Degeneracy resolution, overlap criteria, and state selection

When many stabilizer states attain the same minimum energy, a purely energetic definition is insufficient for simulation tasks that need large overlap with the true ground state. The refined OSGS program therefore adds an overlap criterion on top of the stabilizer minimum. Energy is still evaluated through stabilizer expectations,
$$
\langle \phi_G|H|\phi_G\rangle=\sum_{P\in S\cap \langle G\rangle}\mathrm{sgn}(P)\,h_P,
$$
which is computable in time polynomial in $n$ using the stabilizer tableau, because membership $P\in \langle G\rangle$ is a symplectic linear-algebra query and the eigenvalue $\mathrm{sgn}(P)$ is read off from the tableau phase bit [2603.06286].

For small to medium scale problems, one route uses a robustness-of-magic-inspired linear algebra. Writing the Hamiltonian coefficients as a vector $b$ under the normalized Pauli trace inner product and using the stabilizer-incidence matrix $A$, one obtains
$$
E_{\mathrm{SGS}}=\min_m x_m/2^n,\qquad x=A^\dagger b,
$$
which yields the SGS energy and the minimizing maximally commuting subgroups. Degeneracy is then resolved by a two-step filter: first form the common subgroup $G_{\mathrm{com}}$ across the minimizers and a reduced Hamiltonian $H_{\mathrm{sub}}$; then retain generator sets satisfying
$$
\xi(\mathcal{G},H)=
\begin{cases}
1,& \exists g\in \mathcal{G}\ \text{with}\ [g,H]=0,\\
0,& \text{else}.
\end{cases}
$$
The selection is rerun on $H_{\mathrm{sub}}$ to fix phases and complete the full generator set $\mathcal{G}_{\mathrm{op}}$ [2603.06286].

For large systems, the search is recast as a combinatorial optimization over an $n\times 2n$ binary symplectic tableau with a phase vector. A commutation graph is built on the Hamiltonian-supported Pauli set, and a genetic algorithm optimizes the fitness
$$
\mathrm{Fitness}(\mathcal{G})=-E(\mathcal{G};H)-\lambda_{\mathrm{pen}}N_{\mathrm{nc}}-\mu_{\mathrm{pen}}N_{\mathrm{dep}}+\rho\,\xi(\mathcal{G},H)+\rho_{\mathrm{sub}}\xi(\mathcal{G},H_{\mathrm{sub}}),
$$
with an additional fidelity term in benchmarks where the true ground state is known. The resulting classical preprocessing scales as
$$
O\!\big(M\,G_{\max}(p+n^2)\big),
$$
for population size $M$, generations $G_{\max}$, and $p=|S|$ nonzero Pauli terms. Once the generator set is fixed, standard tableau synthesis yields a Clifford circuit $U_G$ with worst-case $O(n^2/\log n)$ gates and $O(n^2)$ classical synthesis time [2603.06286].

This overlap-sensitive perspective aligns naturally with older stabilizer-manipulation results. Efficient inner-product algorithms compute overlaps between stabilizer states in $O(n^3)$ time in general and show quadratic behavior for many practical instances, which is directly relevant whenever one wishes to compare multiple stabilizer minimizers against a target state or against each other [1210.6646].

## 5. Applications in quantum simulation and many-body physics

The stabilizer-ground-state framework was proposed not only as a variational endpoint but also as a building block for broader simulation workflows. In 1D many-body problems, stabilizer ground states can track phases and topological features. A generalized 1D cluster model,
$$
H=\sum_n\big(-X_{n-1}Z_nX_{n+1}-J_y Y_nY_{n+1}+h_y Y_n\big),
$$
shows cluster, ferromagnetic, and polarized phases consistent with DMRG, with simple stabilizer generators in each phase. For the 2D toric code in fields,
$$
H=-(\sum_v A_v+\sum_p B_p)-h_x\sum_j X_j-h_z\sum_j Z_j,
$$
one can scan single-qubit rotations $U(\alpha,\beta)$ and compute the stabilizer ground state of $H'(\alpha,\beta)=U(\alpha,\beta)^\dagger H U(\alpha,\beta)$, producing an extended stabilizer ground state phase diagram consistent with Monte Carlo and capturing transitions along the $h_x=h_z$ line qualitatively [2403.08441].

These states also function as algorithmic reference states. Stabilizer initializations improve VQE energy optimization versus zero-state or product-state starts; stabilizer CMCS can guide a splitting $H=H_0+H_1$ in which $H_0$ aligns with stabilizers and $H_1$ is treated perturbatively; and the framework has stated synergies with stabilizer tensor networks, Clifford-augmented DMRG, and stabilizer-based Monte Carlo [2403.08441].

In measurement-based deterministic imaginary time evolution, the OSGS is used as the anchor state. Imaginary-time evolution
$$
|\psi(\beta)\rangle=\frac{e^{-\beta H}|\psi_{\mathrm{init}}\rangle}{\|e^{-\beta H}|\psi_{\mathrm{init}}\rangle\|}
$$
is implemented through weak measurements and deterministic feedback. The threshold is set to the stabilizer-ground-state energy $E_{\mathrm{th}}=E_{\mathrm{SGS}}$, and when the running energy estimate exceeds $E_{\mathrm{th}}$, a Clifford reset $U_c$ reprepares $|\phi^\ast\rangle$. Because $|\phi^\ast\rangle$ maximizes the initial ground-state overlap inside the SGS manifold, fewer measurements are needed to make the $|E_0\rangle$ component dominant [2603.06286].

A closely related stabilizer-first strategy appears in stabilizer-accelerated many-body estimation. There the Hamiltonian is split as
$$
H=H_{\mathrm{stab}}+H_{\mathrm{res}},
$$
with $H_{\mathrm{stab}}$ determined by an optimal stabilizer group and $H_{\mathrm{res}}$ carrying the non-stabilizer residue. In the Lipkin–Meshkov–Glick model,
$$
\tilde H=\frac12\sum_i Z_i-\frac{v_x}{2(N-1)}\sum_{i<j}(X_iX_j+\chi Y_iY_j),
$$
the optimal stabilizer ground state is $|1\rangle^{\otimes N}$ for $v_x<2$ and an equal-weight parity superposition stabilized by $\langle X_1X_N,\ldots,X_{N-1}X_N,(-1)^N\prod_i Z_i\rangle$ for $v_x>2$. In this setting, stabilizer ground states accelerate imaginary-time convergence because the initial overlap with the exact ground state is larger, and a single non-unitary reweighting $e^{-\theta J_z}$ captures much of the remaining correlation [2505.02923].

Molecular applications further broaden the scope. Stabilizer configuration interaction finds the best stabilizer approximations to molecular ground states up to 36 qubits and constructs generalized stabilizer states that improve the approximation further, while adaptive SCI scales as $O(n^3)$ [2410.21125]. Independent stabilizer-approximation studies report that water in STO-3G requires 14 qubits and benzene 72 qubits, and that the best stabilizer states approximate the true ground states very well, especially when the molecules are strongly distorted [2302.11734].

## 6. Preparation, verification, limitations, and related directions

Once an optimal stabilizer ground state has been identified, it remains to prepare and certify it. For 1D $k$-local Hamiltonians with an $n-s$ generator output, the stabilizer ground states can be prepared with
$$
O\!\big(n k'/\log k'\big)
$$
single- and two-qubit Clifford gates, where $k'=k+s$; in the non-degenerate case this becomes $O(nk/\log k)$ [2403.08441]. More generally, any stabilizer state is locally Clifford equivalent to a graph state, and graph-decimation methods synthesize preparation circuits by reducing graph edges with two-qubit Clifford gates. The AI-guided method QuSynth combines reinforcement learning and Monte Carlo tree search and reports reductions in two-qubit gate count by up to a factor of $2.5$ compared to previous approaches while retaining low circuit depth [2603.17743].

Experimental certification is also unusually structured. For a target stabilizer state $\rho=|\Psi\rangle\langle\Psi|$, a Pauli-measurement verification scheme has acceptance operator
$$
\Omega=\sum_j p_j E_j,
$$
spectral gap $\nu(\Omega)=1-\beta(\Omega)$, and exact sample complexity
$$
N(\epsilon,\delta,\Omega)=\left\lceil \frac{\ln\delta}{\ln[1-\nu(\Omega)\epsilon]} \right\rceil.
$$
For any entangled stabilizer state, any separable-measurement protocol satisfies the universal bound $\nu(\Omega)\le 2/3$, yielding
$$
N\ge \left\lceil \frac{\ln\delta}{\ln(1-2\epsilon/3)} \right\rceil \approx \left\lceil \frac{3}{2\epsilon}\ln(1/\delta)\right\rceil.
$$
Optimal Pauli-measurement protocols saturate this bound for all connected graph states up to seven qubits, and X/Z-only protocols satisfy $\nu\le 1/2$ [2007.09713].

The main conceptual limitation is explicit. The method is exact within the stabilizer ansatz, not for arbitrary quantum states. In noncommuting or frustrated Pauli Hamiltonians, or whenever the true ground state is strongly non-stabilizer, the optimal stabilizer ground state need not coincide with the physical ground state. The formalism still returns the best stabilizer state and its energy, but not a claim of exactness beyond that ansatz [2403.08441][2603.06286].

A second limitation is computational. For 1D local sparse Hamiltonians the state-machine algorithm is exact and efficient, but in higher dimensions exact computation is generally NP-hard, so approximate or heuristic extensions may be needed [2403.08441]. This suggests a practical division of labor: exact CMCS-based methods where locality and dimensionality permit them, and overlap-aware heuristic or hybrid procedures when the search space or non-stabilizer residue becomes dominant.

In that sense, the optimal stabilizer ground state is best understood as a certified Clifford backbone for ground-state estimation: exact for commuting-stabilizer instances, variationally controlled for general Pauli Hamiltonians, and frequently useful as the reference state from which more expressive classical or quantum refinements proceed.

Source: https://www.emergentmind.com/topics/optimal-stabilizer-ground-state