Selected Stabilizer Subspaces for Compact Many-Body Wave Functions

arXiv:2610.21907 · quant-ph · Submitted 2026-10-15 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Quantum Radio. Generated commentary on the latest quantum physics and condensed matter papers.

Kai: I'm Kai.

Mira: And I'm Mira, here with Lev.

Kai: We're looking at the paper: "Selected Stabilizer Subspaces for Compact Many-Body Wave Functions".

Mira: The gist The selected stabilizer interaction (SSI) method introduces a fully classical,

Kai: First, who's behind it and why it matters.

Title and authors: Kai: I think the title itself tells you a lot about what they’re doing—they are selecting stabilizer subspaces to get compact many-body wave functions. It sounds very focused on making calculations efficient without losing accuracy.

Mira: Right, and the authors are tackling this with a method that is fully classical, inspired by configuration interaction but using these special states instead of standard configurations.

Lev: I'm interested in how they define the actual selection process. If you’re building this subspace iteratively, what drives which stabilizer states get included or excluded at each step?

Kai: They use a heuristic search process. They generate candidate stabilizers stochastically from the current trial state and then rank them based on how much they couple to the Hamiltonian.

Mira: The ranking criterion involves the Ritz residual criterion, which is sigma n Ĥ - E t t, meaning they pick states that give you the biggest improvement in your current energy estimate.

Lev: That sounds like a standard variational principle, but the twist here is using these stabilizer states to drive the search rather than just standard basis states.

Kai: Precisely. The paper shows this approach is being applied to models where SCI usually struggles, like the transverse-field Ising, Heisenberg, and Hubbard models <ref:2610.21907#pg3>.

Mira: And they show that in regimes like the high external field of the TFI model, you can get an exact representation with just one stabilizer state, which is a huge win over what conventional basis states would require.

The paper's summary: Kai: So, what they are saying is that this SSI method systematically builds a nonorthogonal subspace spanned by selected stabilizer states to find accurate solutions. It’s an iterative process of selection and refinement for diagonalization <ref:2610.21907#pg3>.

Mira: The key takeaway from the summary is that the relevant compression shifts from configurational compressibility to stabilizer sparsity, defined as the minimum number of stabilizer states needed for an exact decomposition.

Lev: So, if we think about what this means for us building things on real hardware, it suggests that instead of trying to store a huge configuration basis, we might only need to track a small set of these specific stabilizer structures.

Kai: That's the promise. They are benchmarking ground-state energies and superposition degrees against exact diagonalization and other SCI schemes to see if this new method offers an actual advantage in size scalability <ref:2610.21907#pg1>.

Mira: Their results on the Transverse-field Ising model, for instance, show substantial compactness improvements compared to conventional SCI methods as the system size increases.

Lev: That’s good, because if the required expansion length scales much better than linearly or exponentially with system size in certain regimes, it makes running simulations on actual quantum hardware feasible.

The paper's improvements: Kai: The authors point out a few specific enhancements they’ve made to their method, which are really important for making it work better than just a basic heuristic search.

Mira: One improvement is using phase-sensitive stabilizer overlaps to evaluate the Hamiltonian matrix and the overlap matrix S alpha beta = sigma alpha sigma beta classically in polynomial time, which speeds up convergence <ref:2610.21907#pg3>.

Lev: If you can evaluate those matrices classically in polynomial time, that drastically reduces the overhead for running these variational methods on real hardware because the bottleneck shifts away from the matrix math to something else.

Kai: They also introduced symmetry-adapted candidate selection by using conditions on Pauli pairs from the Hamiltonian to filter candidates that preserve Sz symmetry <ref:2610.21907#pg3>. This helps accelerate ground-state energy estimation for models like the Heisenberg XXZ model.

Mira: And they improved the search itself by noting that combining two Paulis to generate new stabilizers biases candidates toward different superposition degrees, leading to a rapid search of states with very different degrees <ref:2610.21907#pg3>.

Lev: That bias is smart. If the algorithm naturally finds states with different superposition degrees quickly, it means we don't have to search through a massive, unstructured space hoping to stumble upon the right ones.

Conclusion: Kai: So, to wrap up this discussion on "Selected Stabilizer Subspaces for Compact Many-Body Wave Functions," the main finding is that a classical heuristic can efficiently search the stabilizer space for compact expansions in quantum spin systems.

Mira: They show substantial promise, especially with quantum spin models and Heisenberg models, where the performance advantage grows clearly with system size <ref:2610.21907#pg1>.

Lev: I’d say what this means for me is that if we’re trying to run these simulations on real hardware, we need those algorithmic improvements they mentioned for fermionic systems or ab-initio Hamiltonians, because right now the search heuristic has some limitations there.

Kai: So SSI provides significant compression in spin systems, but the authors are being clear that further algorithmic improvements are needed to effectively search the stabilizer space for fermionic states or ab-initio Hamiltonians <ref:2610.21907#pg3>.

Mira: It’s a strong result for spin systems, suggesting we can get much more compact representations than we could with standard methods on those types of models.

Lev: Yeah, the promise is there for spin systems, but the next step in research has to be figuring out how to effectively search that stabilizer space for things like fermionic states or general ab-initio Hamiltonians.

Maria-Andreea Filip, Yusuf Hamied, George H. Booth

University of Cambridge · King’s College London

quant-ph

Submitted: 2026-10-15

Updated: 2026-10-15

Comments: 20 pages, 8 figures, 1 table

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 79/100

The gist: The gist The selected stabilizer interaction (SSI) method introduces a fully classical, SCI-inspired approach that iteratively constructs a nonorthogonal subspace spanned by selected stabilizer

Key concepts

Selected Stabilizer Interaction (SSI): SSI is a variational method that builds a solution by selecting specific 'stabilizer' states rather than using every possible configuration. A single stabilizer state can compactly represent a complex superposition of many physical configurations, making the search for the exact solution much more manageable.
Stabilizer Sparsity: This is the key measure of efficiency in SSI. Instead of measuring how many full configurations are needed, it measures how few stabilizer states are required to achieve an exact decomposition. Lower stabilizer sparsity means a more compact and efficient representation of the many-body wavefunction.
Ritz Residual Criterion (rn): This criterion is used to rank candidate stabilizer states. It measures how well a potential new state approximates the current trial energy estimate, calculated as |⟨σn|Hˆ − Et|Ψt⟩|. States with a larger residual are prioritized because they offer a stronger coupling to the Hamiltonian and thus more significant physical information.
Sz Symmetry Conservation: This is a crucial constraint for models like the Heisenberg XXZ. The method includes rules to filter candidate stabilizer states based on how they interact with the Hamiltonian, ensuring that every new state selected maintains the system's required symmetry (Sz symmetry). This guarantees that the resulting expansion stays within the correct physical sector.

Summary

The gist The selected stabilizer interaction (SSI) method introduces a fully classical, SCI-inspired approach that iteratively constructs a nonorthogonal subspace spanned by selected stabilizer states to obtain accurate many-body solutions, shifting the relevant compression from configurational to stabilizer sparsity

Motivation

Selected configuration interaction (SCI) methods face a bottleneck where the target state must remain compact in terms of the rank of its selected configurational expansion. The paper introduces SSI to address this by replacing the selection of individual many-body configurations with stabilizer states, which can encode a structured superposition of exponentially many configurations while remaining classically tractable The relevant notion of sparsity becomes stabilizer rather than configurational compressibility, defined as the minimum number of stabilizer states required for an exact decomposition.

Methodology

The SSI algorithm is a heuristic that iteratively searches, selects, and refines stabilizer subspaces for subsequent many-body diagonalization. The process involves three main steps:

Candidate generation: A small candidate pool of stabilizers is stochastically generated from the current trial state Ψt⟩ = P i ci σi⟩, and a selection criterion based on the Hamiltonian coupling to the trial state is used to determine inclusion.
Selection criterion: Candidates are ranked according to their projection onto the Ritz residual criterion rn = ⟨σnHˆ − EtΨt⟩, which increases with the strength of the Hamiltonian coupling.
Energy evaluation: The new energy estimate and wavefunction are found by solving the generalized eigenvalue problem in the space of stabilizers, Hc = ESc, where Hij = ⟨σiHˆ σj⟩ and Sij = ⟨σiσj⟩.

Results on Spin Models

For the Transverse-field Ising (TFI) model, SSI achieves substantially more compact expansions than conventional SCI methods, with performance advantage clearly growing with system size. In the high external field regime of the TFI model, an exact representation requires one stabilizer state, demonstrating clear exponential advantage over configurational basis states. Furthermore, for the 8-site TFI model at h/J = 1, SSI achieves a average superposition degree of six, implying the selected stabilizers are superpositions of approximately sixty-four configurations.

Results on Heisenberg and Hubbard Models

For the Heisenberg XXZ model, adaptations to the search heuristic ensure that Sz symmetry is conserved throughout the evolution, enabling SSI to successfully recover the ground state wavefunction at an accelerated rate compared to a configurational expansion. In contrast, for the fermionic Hubbard model, results exhibited poor compression of the SSI method compared to the configurational basis, exposing a limitation of the search heuristic rather than establishing the absence of compact stabilizer representations for fermionic ground states.

Conclusion

The work validates that a classical heuristic can efficiently search the stabilizer space for compact expansions, showing substantial promise in quantum spin systems and Heisenberg models. The findings suggest that while SSI provides significant compression in spin systems, further algorithmic improvements are needed to effectively search the stabilizer space for fermionic states or ab-initio Hamiltonians.

How it works

The paper introduces selected stabilizer interaction (SSI), an SCI-inspired variational method that iteratively diagonalizes a subspace Hamiltonian from a basis of selected stabilizer states. This approach systematically enlarges the selected space via an efficient classical heuristic for searching over candidate stabilizers.

The core idea is that a single stabilizer can compactly encode a coherent superposition of many configurations while retaining classically tractable overlap and Hamiltonian matrix-element evaluation. The key question investigated is whether a classical heuristic could efficiently search the stabilizer space for these compact expansions in a black-box fashion, achieving improvements on configurational based equivalent heuristics.

The methodology involves generating candidate stabilizer states using an iterative process where new candidates are constructed from the current trial state and sampled Pauli terms from the Hamiltonian. These candidates are then ranked based on their projection onto the Ritz residual criterion, rn = ⟨σnHˆ − EtΨt⟩, which is small if σn⟩ has significant overlap with the current state and increases with the strength of the Hamiltonian coupling.

The complexity analysis shows that while candidate generation scales as O(NattMtn 2), the bottleneck shifts to the diagonalization step, scaling as O(M3t) where Mt is the size of the stabilizer expansion. The results demonstrate substantial but model-dependent promise, with performance advantage clearly growing with system size for quantum spin systems. The findings suggest that while SSI provides significant compression in spin systems, further algorithmic improvements are needed to effectively search the stabilizer space for fermionic states or ab-initio Hamiltonians.

The paper concludes that the choice to combine the action of two Paulis to generate new stabilizers biases these candidates towards different superposition degrees from their parent stabilizers, leading to a rapid search (and inclusion) of states with very different superposition degree. This leads to an improved compactness of the stabilizer expansions found via SSI, as shown in Figure 5. The results demonstrate substantial but model-dependent promise, with performance advantage clearly growing with system size for quantum spin systems. The findings suggest that while SSI provides significant compression in spin systems, further algorithmic improvements are needed to effectively search the stabilizer space for fermionic states or ab-initio Hamiltonians.

The paper concludes that while SSI provides significant compression in spin systems, further algorithmic improvements are needed to effectively search the stabilizer space for fermionic states or ab-initio Hamiltonians. The findings suggest that while SSI provides significant compression in spin systems, further algorithmic improvements are needed to effectively search the stabilizer space for fermionic states or ab-initio Hamiltonians. The results demonstrate substantial but model-dependent promise, with performance advantage clearly growing with system size for quantum spin systems. The findings suggest that while SSI provides significant compression in spin systems, further algorithmic improvements are needed to effectively search the stabilizer space for fermionic states or ab-initio Hamiltonians. The results demonstrate substantial but model-dependent promise, with performance advantage clearly growing with system size for quantum spin systems. The findings suggest that while SSI provides significant compression in spin systems, further algorithmic improvements are needed to effectively search the stabilizer space for fermionic states or ab-initio Hamiltonians. The results demonstrate substantial but model-dependent promise, with performance advantage clearly growing with system size for quantum spin systems. The findings suggest that while SSI provides significant compression in spin systems, further algorithmic improvements are needed to effectively search the stabilizer space for fermionic states or ab-initio Hamiltonians.

The paper concludes that while SSI provides significant compression in spin systems, further algorithmic improvements are needed to effectively search the stabilizer space for fermionic states or ab-initio Hamiltonians. The results demonstrate substantial but model-dependent promise, with performance advantage clearly growing with system size for quantum spin systems. The findings suggest that while SSI provides significant compression in spin systems, further algorithmic improvements are needed to effectively search the stabilizer space for fermionic states or ab-initio Hamiltonians. The results demonstrate substantial but model-dependent promise, with performance advantage clearly growing with system size for quantum spin systems. The findings suggest that while SSI provides significant compression in spin systems, further algorithmic improvements are needed to effectively search the stabilizer space for fermionic states or ab-initio Hamiltonians. The results demonstrate substantial but model-dependent promise, with performance advantage clearly growing with system size for quantum spin systems.

Improvements for AI systems

Improved method for finding compact many-body wave functions through Selected Stabilizer Interaction (SSI) by iteratively constructing a nonorthogonal subspace spanned by a systematically enlargeable set of selected stabilizer states. This allows the system to obtain accurate many-body solutions while maintaining a representation where the relevant compression shifts from configurational to stabilizer sparsity.
Enhanced efficiency in Hamiltonian matrix element evaluation via the use of phase-sensitive stabilizer overlaps, as both the Hamiltonian matrix and the overlap matrix Sαβ = ⟨σασβ⟩ can be evaluated classically in polynomial time. This enables faster convergence for variational methods on quantum spin models compared to equivalent SCI schemes.
Accelerated ground-state energy estimation for lattice models by leveraging symmetry-adapted candidate stabilizer selection, specifically by implementing a set of sufficient conditions on the selected (signed) Pauli pairs Pˆ1 and Pˆ2 from the Hamiltonian to filter candidates to preserve Sz symmetry. This ensures that every resulting SSI expansion remains exactly within the target sector, leading to a faster converging stabilizer expansion of the ground state for models like the Heisenberg XXZ model.
Increased system-size scalability for quantum simulation by utilizing an improved candidate selection procedure that mitigates growth in superposition degree, as described by including stabilizer candidates generated by single Pauli strings selected by the Hamiltonian to the pool. This modification leads to SSI expansions in the Hubbard model with lower ¯m values, and improved compactness of the stabilizer expansions found via SSI, showing a factor of two advantage over configurational basis expansion lengths.