Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2608.14798

Stabilizer Statistical Mechanics: A Framework for Efficient Quantification and Classification of Magic States

William E. Salazar, Gaurav Saxena, Jack S. Baker, Leong Chuan Kwek, Thi Ha Kyaw

quant-phSubmitted Aug 14, 2026• Announced Aug 18, 2026• 0 citations
arXivHTMLPDF

Abstract

The partition function is statistical mechanics' answer to an exponentially large spectrum, distilling it into a single analytic object whose temperature dependence resolves the full structure of the underlying ensemble. We show that magic, the resource separating universal quantum computation from classically simulable stabilizer dynamics, admits precisely such a description. Mapping the Pauli spectrum of a quantum state onto the energy levels of a fictitious many-body system, the Pauli gas, we construct its canonical partition function, the stabilizer partition function, and from its associated free energy a magic monotone that we call the stabilizer work. Both these objects are analytic functions of an inverse-temperature-like parameter and are efficiently estimable via Bell sampling. The framework is analytically tractable. We derive exact ensemble-averaged partition functions for Haar-random, $ν$-compressible, and pseudomagic states, together with concentration guarantees. We show that for every value of its parameter, the stabilizer work is a faithful, Subadditive magic monotone, while remaining efficiently accessible on quantum hardware and admitting an operational interpretation. Unlike measures that probe a single moment of the Pauli distribution, the stabilizer work is intrinsically moment-generating. As the temperature is tuned from high to low, it interpolates continuously between the stabilizer 2-Rényi entropy and the stabilizer nullity, revealing two previously disconnected monotones as limiting cases of a single object. We demonstrate the framework on low-rank stabilizer simulation, resource interconversion, molecular ground states, and quantum many-body systems. A thermodynamics of magic is therefore not merely an analogy but a working toolkit, opening a statistical-mechanical route to magic properties of quantum systems that no single measure can access.

Comments
42 pages, 11 figures. Comments are welcome. William E. Salazar and Gaurav Saxena contributed equally to this work

Related papers

Universal equilibrium magic in quantum many-body systems

Soumyadeep Sarma, Tobias Haug, John Preskill +1 more

Thermalization conventionally describes local properties of isolated many-body systems at equilibrium. Magic, or nonstabilizerness $\unicode{x2013}$ the resource enabling universal quantum computation $\unicode{x2013}$ is by contrast encoded in the global structure of the many-body wavefunction. We show that, despite its global nature, the magic of equilibrium pure states of chaotic many-body systems, including late-time evolved states and energy eigenstates, is universal: it is captured by the thermal Scrooge ensemble, the minimally informative ensemble of pure states consistent with the Gibbs state at the same effective temperature. Therefore, for systems with no conserved quantities other than the total energy, equilibrium magic is a function of temperature alone, independent of the initial state and other microscopic features of the equilibrium state. This yields concrete universal predictions for the stabilizer Rényi entropies (SREs). At infinite temperature, the SRE is set by Haar-like fluctuations of the Pauli spectrum, while at finite temperature energy conservation induces a volume-law thermodynamic correction controlled by the thermal Pauli spectrum. We support these predictions with analytical arguments and extensive numerical simulations. We further show that chaotic many-body systems at high temperatures possess long-range magic and entanglement that cannot be removed by finite-depth local quantum circuits. Our results establish magic as a thermodynamic property of chaotic many-body systems and suggest that Scrooge ensembles may provide a unified framework for quantum many-body resources.

quant-pharXiv:2608.22939•Aug 25, 2026•0 citations

Pauli spectrum and nonstabilizerness of random fermionic Gaussian states

Xhek Turkeshi, Piotr Sierant, Poetri Sonya Tarabunga

We characterize the magic of random fermionic Gaussian states through their Pauli spectrum and their stabilizer entropies. For the Haar ensemble of Majorana Gaussian states we derive closed finite-$N$ expressions for the average stabilizer purities and reconstruct the Pauli spectrum exactly, as a mixture of products of beta-distributed variables resolved by Majorana weight. The stabilizer entropies, and their filtered versions, freeze at Rényi index $q_c=2$: for $q>2$ the magic density is $1/(q-1)$, controlled by rare weight-two Majorana strings. Resolving these contributions requires exponentially many characteristic samples, making direct sampling estimates of the stabilizer entropies inefficient in this frozen regime. For Haar-random Gaussian states at fixed filling, we express the averaged stabilizer purities at positive integer $q$ as coefficient integrals over $\mathrm{SU}(2q)$, with dimension independent of $N$, and evaluate the $q=2$ case exactly. We establish that the frozen density describes typical states in the Majorana and half-filled number-conserving ensembles and relate these ensembles to the corresponding $\mathrm{SYK}_2$ ground states. Finally, we show that typical fermionic Gaussian states retain certified magic after discarding any fixed fraction $κ<κ_\star^{\rm G}\simeq0.7613$ of their modes in the thermodynamic limit, remaining magical even beyond the $2/3$ threshold for Haar-random states in the full Hilbert space.

quant-pharXiv:2609.15739•Sep 15, 2026•0 citations

Exponentially Accelerated Sampling of Pauli Strings for Nonstabilizerness

Zhenyu Xiao, Shinsei Ryu

Quantum magic, quantified by nonstabilizerness, measures departures from stabilizer structure and underlies potential quantum speedups. We introduce an efficient classical framework for computing stabilizer Rényi entropies and stabilizer nullity of generic $N$-qubit wavefunctions. The method combines the fast Walsh-Hadamard transform with an exact partition of Pauli operators, reducing the average cost per sampled Pauli string from $\mathcal{O}(2^N)$ to $\mathcal{O}(N)$. We further develop a Monte Carlo estimator with Clifford preconditioning and find that the required number of samples shows no visible growth with $N$ in our benchmarks. Applying the method to $T$-doped random Clifford circuits, we identify the scrambling ratio $η$ (Clifford gates per $T$ gate) as the key parameter governing magic growth. Each $T$ gate approaches its dilute-limit nonstabilizerness power with only modest Clifford scrambling. Our approach enables quantitative studies of magic in highly entangled states and long-time nonequilibrium dynamics.

quant-pharXiv:2601.00761•Jan 5, 2026•0 citations

Operational interpretation of the Stabilizer Entropy

Lennart Bittel, Lorenzo Leone

Magic-state resource theory is a fundamental framework with far-reaching applications in quantum error correction and the classical simulation of quantum systems. Recent advances have significantly deepened our understanding of magic as a resource across diverse domains, including many-body physics, nuclear and particle physics, and quantum chemistry. Central to this progress is the stabilizer Rényi entropy, a computable and experimentally accessible magic monotone. Despite its widespread adoption, a rigorous operational interpretation of the stabilizer entropy has remained an open problem. In this work, we provide such an interpretation in the context of quantum property testing. By showing that the stabilizer entropy is the most robust measurable magic monotone, we demonstrate that the Clifford orbit of a quantum state becomes exponentially indistinguishable from Haar-random states, at a rate governed by the stabilizer entropy M α ( ψ ) and the number of available copies. This implies that the Clifford orbit forms an approximate state k -design, with an approximation error exp ⁡ ( − Θ ( M α ( ψ ) ) ) for α ≥ 2 . Conversely, we establish that the optimal probability of distinguishing a given quantum state from the set of stabilizer states is also governed by its stabilizer entropy. These results reveal that the stabilizer entropy quantitatively characterizes the transition from stabilizer states to universal quantum states, thereby offering a comprehensive operational perspective of the stabilizer entropy as a quantum resource.

arXiv:2507.22883•Quantum•Apr 15, 2026•8 citations

Quantum magic dynamics in random circuits

Yuzhen Zhang, Yingfei Gu

Magic refers to the degree of “quantumness” in a system that cannot be fully described by stabilizer states and Clifford operations alone. In quantum computing, stabilizer states and Clifford operations can be efficiently simulated on a classical computer, even though they may appear complicated from the perspective of entanglement. In this sense, magic is a crucial resource for unlocking the unique computational power of quantum computers to address problems that are classically intractable. Magic can be quantified by measures such as Wigner negativity and mana that satisfy fundamental properties such as monotonicity under Clifford operations. In this paper, we generalize the statistical mechanical mapping methods of large- q random circuits to the calculation of Rényi Wigner negativity and mana. Based on this, we find: (1) a precise formula describing the competition between magic and entanglement in many-body states prepared under Haar random circuits; (2) a formula describing the spreading and scrambling of magic in states evolved under random Clifford circuits; (3) a quantitative description of magic “squeezing” and “teleportation” under measurements. Finally, we comment on the relation between coherent information and magic.

arXiv:2410.21128•npj Quantum Information•May 14, 2026•7 citations

Exponentially Accelerated Sampling of Pauli Strings for Nonstabilizerness

Anonymous, Shinsei Ryu

Quantum magic, quantified by nonstabilizerness, measures departures from stabilizer structure and underlies potential quantum speedups. We introduce an efficient classical framework for computing stabilizer Rényi entropies and stabilizer nullity of generic N-qubit wave functions. The method combines the fast Walsh-Hadamard transform with an exact partition of Pauli operators, reducing the average cost per sampled Pauli string from O(2^{N}) to O(N). We further develop a Monte Carlo estimator with Clifford preconditioning and find that the required number of samples shows no visible growth with N in our benchmarks. Applying the method to T-doped random Clifford circuits, we identify the scrambling ratio η (Clifford gates per T gate) as the key parameter governing magic growth. Each T gate approaches its dilute-limit nonstabilizerness power with only modest Clifford scrambling. Our approach enables quantitative studies of magic in highly entangled states and long-time nonequilibrium dynamics.

Physical Review Letters•Jul 20, 2026•1 citation
The partition function is statistical mechanics' answer to an exponentially large spectrum, distilling it into a single analytic object whose temperature dependence resolves the full structure of the underlying ensemble. We show that magic, the resource separating universal quantum computation from classically simulable stabilizer dynamics, admits precisely such a description. Mapping the Pauli spectrum of a quantum state onto the energy levels of a fictitious many-body system, the Pauli gas, we construct its canonical partition function, the stabilizer partition function, and from its associated free energy a magic monotone that we call the stabilizer work. Both these objects are analytic functions of an inverse-temperature-like parameter and are efficiently estimable via Bell sampling. The framework is analytically tractable. We derive exact ensemble-averaged partition functions for Haar-random, ννν-compressible, and pseudomagic states, together with concentration guarantees. We show that for every value of its parameter, the stabilizer work is a faithful, Subadditive magic monotone, while remaining efficiently accessible on quantum hardware and admitting an operational interpretation. Unlike measures that probe a single moment of the Pauli distribution, the stabilizer work is intrinsically moment-generating. As the temperature is tuned from high to low, it interpolates continuously between the stabilizer 2-Rényi entropy and the stabilizer nullity, revealing two previously disconnected monotones as limiting cases of a single object. We demonstrate the framework on low-rank stabilizer simulation, resource interconversion, molecular ground states, and quantum many-body systems. A thermodynamics of magic is therefore not merely an analogy but a working toolkit, opening a statistical-mechanical route to magic properties of quantum systems that no single measure can access.
Thermalization conventionally describes local properties of isolated many-body systems at equilibrium. Magic, or nonstabilizerness $\unicode{x2013}$ the resource enabling universal quantum computation $\unicode{x2013}$ is by contrast encoded in the global structure of the many-body wavefunction. We show that, despite its global nature, the magic of equilibrium pure states of chaotic many-body systems, including late-time evolved states and energy eigenstates, is universal: it is captured by the thermal Scrooge ensemble, the minimally informative ensemble of pure states consistent with the Gibbs state at the same effective temperature. Therefore, for systems with no conserved quantities other than the total energy, equilibrium magic is a function of temperature alone, independent of the initial state and other microscopic features of the equilibrium state. This yields concrete universal predictions for the stabilizer Rényi entropies (SREs). At infinite temperature, the SRE is set by Haar-like fluctuations of the Pauli spectrum, while at finite temperature energy conservation induces a volume-law thermodynamic correction controlled by the thermal Pauli spectrum. We support these predictions with analytical arguments and extensive numerical simulations. We further show that chaotic many-body systems at high temperatures possess long-range magic and entanglement that cannot be removed by finite-depth local quantum circuits. Our results establish magic as a thermodynamic property of chaotic many-body systems and suggest that Scrooge ensembles may provide a unified framework for quantum many-body resources.
We characterize the magic of random fermionic Gaussian states through their Pauli spectrum and their stabilizer entropies. For the Haar ensemble of Majorana Gaussian states we derive closed finite-NNN expressions for the average stabilizer purities and reconstruct the Pauli spectrum exactly, as a mixture of products of beta-distributed variables resolved by Majorana weight. The stabilizer entropies, and their filtered versions, freeze at Rényi index qc=2q_c=2qc​=2: for q>2q>2q>2 the magic density is 1/(q−1)1/(q-1)1/(q−1), controlled by rare weight-two Majorana strings. Resolving these contributions requires exponentially many characteristic samples, making direct sampling estimates of the stabilizer entropies inefficient in this frozen regime. For Haar-random Gaussian states at fixed filling, we express the averaged stabilizer purities at positive integer qqq as coefficient integrals over SU(2q)\mathrm{SU}(2q)SU(2q), with dimension independent of NNN, and evaluate the q=2q=2q=2 case exactly. We establish that the frozen density describes typical states in the Majorana and half-filled number-conserving ensembles and relate these ensembles to the corresponding SYK2\mathrm{SYK}_2SYK2​ ground states. Finally, we show that typical fermionic Gaussian states retain certified magic after discarding any fixed fraction κ<κ⋆G≃0.7613κ<κ_\star^{\rm G}\simeq0.7613κ<κ⋆G​≃0.7613 of their modes in the thermodynamic limit, remaining magical even beyond the 2/32/32/3 threshold for Haar-random states in the full Hilbert space.
Quantum magic, quantified by nonstabilizerness, measures departures from stabilizer structure and underlies potential quantum speedups. We introduce an efficient classical framework for computing stabilizer Rényi entropies and stabilizer nullity of generic NNN-qubit wavefunctions. The method combines the fast Walsh-Hadamard transform with an exact partition of Pauli operators, reducing the average cost per sampled Pauli string from O(2N)\mathcal{O}(2^N)O(2N) to O(N)\mathcal{O}(N)O(N). We further develop a Monte Carlo estimator with Clifford preconditioning and find that the required number of samples shows no visible growth with NNN in our benchmarks. Applying the method to TTT-doped random Clifford circuits, we identify the scrambling ratio ηηη (Clifford gates per TTT gate) as the key parameter governing magic growth. Each TTT gate approaches its dilute-limit nonstabilizerness power with only modest Clifford scrambling. Our approach enables quantitative studies of magic in highly entangled states and long-time nonequilibrium dynamics.
Magic-state resource theory is a fundamental framework with far-reaching applications in quantum error correction and the classical simulation of quantum systems. Recent advances have significantly deepened our understanding of magic as a resource across diverse domains, including many-body physics, nuclear and particle physics, and quantum chemistry. Central to this progress is the stabilizer Rényi entropy, a computable and experimentally accessible magic monotone. Despite its widespread adoption, a rigorous operational interpretation of the stabilizer entropy has remained an open problem. In this work, we provide such an interpretation in the context of quantum property testing. By showing that the stabilizer entropy is the most robust measurable magic monotone, we demonstrate that the Clifford orbit of a quantum state becomes exponentially indistinguishable from Haar-random states, at a rate governed by the stabilizer entropy M α ( ψ ) and the number of available copies. This implies that the Clifford orbit forms an approximate state k -design, with an approximation error exp ⁡ ( − Θ ( M α ( ψ ) ) ) for α ≥ 2 . Conversely, we establish that the optimal probability of distinguishing a given quantum state from the set of stabilizer states is also governed by its stabilizer entropy. These results reveal that the stabilizer entropy quantitatively characterizes the transition from stabilizer states to universal quantum states, thereby offering a comprehensive operational perspective of the stabilizer entropy as a quantum resource.
Magic refers to the degree of “quantumness” in a system that cannot be fully described by stabilizer states and Clifford operations alone. In quantum computing, stabilizer states and Clifford operations can be efficiently simulated on a classical computer, even though they may appear complicated from the perspective of entanglement. In this sense, magic is a crucial resource for unlocking the unique computational power of quantum computers to address problems that are classically intractable. Magic can be quantified by measures such as Wigner negativity and mana that satisfy fundamental properties such as monotonicity under Clifford operations. In this paper, we generalize the statistical mechanical mapping methods of large- q random circuits to the calculation of Rényi Wigner negativity and mana. Based on this, we find: (1) a precise formula describing the competition between magic and entanglement in many-body states prepared under Haar random circuits; (2) a formula describing the spreading and scrambling of magic in states evolved under random Clifford circuits; (3) a quantitative description of magic “squeezing” and “teleportation” under measurements. Finally, we comment on the relation between coherent information and magic.
Quantum magic, quantified by nonstabilizerness, measures departures from stabilizer structure and underlies potential quantum speedups. We introduce an efficient classical framework for computing stabilizer Rényi entropies and stabilizer nullity of generic N-qubit wave functions. The method combines the fast Walsh-Hadamard transform with an exact partition of Pauli operators, reducing the average cost per sampled Pauli string from O(2^{N}) to O(N). We further develop a Monte Carlo estimator with Clifford preconditioning and find that the required number of samples shows no visible growth with N in our benchmarks. Applying the method to T-doped random Clifford circuits, we identify the scrambling ratio η (Clifford gates per T gate) as the key parameter governing magic growth. Each T gate approaches its dilute-limit nonstabilizerness power with only modest Clifford scrambling. Our approach enables quantitative studies of magic in highly entangled states and long-time nonequilibrium dynamics.