Exponentially Accelerated Sampling of Pauli Strings for Nonstabilizerness
Pith reviewed 2026-05-16 17:48 UTC · model grok-4.3
The pith
A classical method samples Pauli strings for nonstabilizerness at linear cost in qubit number rather than exponential.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By combining the fast Walsh-Hadamard transform with an exact partition of the Pauli operators, the average computational cost to sample a Pauli string drops from O(2^N) to O(N), and a Clifford-preconditioned Monte Carlo estimator converges with a number of samples that shows no visible growth with N on T-doped Clifford circuits.
What carries the argument
Fast Walsh-Hadamard transform applied to an exact partition of Pauli operators, which generates the probability distribution over Pauli strings needed for the entropies at linear cost per sample.
If this is right
- Stabilizer Rényi entropies become computable for qubit numbers previously out of reach for exact methods.
- Magic growth in hybrid Clifford-T circuits can be tracked quantitatively through the single parameter η.
- Each T gate contributes close to its maximum nonstabilizerness once modest Clifford scrambling is present.
- Long-time nonequilibrium dynamics of nonstabilizerness can be followed in large systems.
Where Pith is reading between the lines
- The linear scaling may make quantitative magic studies feasible for systems with hundreds of qubits if sample count remains bounded.
- Partitioning techniques of this kind could accelerate other computations that sum over Pauli operators.
- If the N-independent sample count holds for typical entangled states, similar preconditioning may help sampling problems in other quantum measures.
Load-bearing premise
The Clifford-preconditioned Monte Carlo estimator stays unbiased and requires a number of samples independent of N for generic highly entangled states beyond the T-doped Clifford circuits examined.
What would settle it
For a random highly entangled state at N=20, run the estimator and observe that the variance of the entropy estimate or the number of samples required for convergence grows exponentially with N.
Figures
read the original abstract
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\'enyi 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 $\eta$ (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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces an efficient classical algorithm for computing stabilizer Rényi entropies and stabilizer nullity of N-qubit states. It combines the fast Walsh-Hadamard transform with an exact partition of Pauli operators to reduce the average cost per sampled Pauli string from O(2^N) to O(N), and develops a Clifford-preconditioned Monte Carlo estimator. Benchmarks on T-doped random Clifford circuits show that the number of required samples exhibits no visible growth with N; the scrambling ratio η is identified as the parameter controlling magic growth, with each T gate approaching its dilute-limit nonstabilizerness contribution under modest Clifford scrambling.
Significance. If the Monte Carlo estimator's variance remains controlled beyond the tested family, the O(N) per-sample cost would enable quantitative studies of nonstabilizerness in large-scale entangled states and long-time dynamics that are inaccessible to exact methods. The empirical N-independence of sample count for T-doped Clifford circuits is a practical strength that could accelerate research on quantum magic.
major comments (2)
- [Benchmarks and Monte Carlo estimator description] The central claim that the Clifford-preconditioned Monte Carlo estimator requires an N-independent number of samples rests on benchmarks for T-doped random Clifford circuits only. No general variance bound or convergence proof is supplied for arbitrary states (e.g., Haar-random or Hamiltonian ground states), which is load-bearing for the assertion of applicability to generic highly entangled wavefunctions.
- [Algorithm and numerical results sections] Validation of the estimator against exact small-N results and explicit error bounds or variance estimates for the Monte Carlo sampling are not provided, leaving the reliability of the reported N-independence insufficiently quantified.
minor comments (2)
- [Abstract] The abstract states that 'the required number of samples shows no visible growth with N' without specifying the range of N tested or the number of circuit realizations used.
- [Introduction and method] Notation for the scrambling ratio η and the precise definition of the Pauli partition should be introduced earlier and used consistently throughout.
Simulated Author's Rebuttal
We thank the referee for the constructive report and positive assessment of the method's potential. We address each major comment below, clarifying the scope of our claims and proposing targeted revisions.
read point-by-point responses
-
Referee: The central claim that the Clifford-preconditioned Monte Carlo estimator requires an N-independent number of samples rests on benchmarks for T-doped random Clifford circuits only. No general variance bound or convergence proof is supplied for arbitrary states (e.g., Haar-random or Hamiltonian ground states), which is load-bearing for the assertion of applicability to generic highly entangled wavefunctions.
Authors: We agree that the observed N-independence of the sample count is an empirical finding specific to the T-doped random Clifford circuit family studied in the benchmarks. The manuscript does not claim or provide a general variance bound or convergence proof that would hold for arbitrary states such as Haar-random states or Hamiltonian ground states. The algorithm and Clifford-preconditioned estimator are formulated for generic wavefunctions, but the sample complexity necessarily depends on the state's nonstabilizerness structure; we will revise the text in the abstract, introduction, and numerical results sections to explicitly state that the N-independence is demonstrated for the tested family and to qualify the broader applicability accordingly. revision: partial
-
Referee: Validation of the estimator against exact small-N results and explicit error bounds or variance estimates for the Monte Carlo sampling are not provided, leaving the reliability of the reported N-independence insufficiently quantified.
Authors: We will add a new subsection in the numerical results that validates the Monte Carlo estimator against exact computations for small N (up to N=8) across several circuit instances, including direct comparisons of the estimated stabilizer Rényi entropies and nullity. We will also include explicit variance estimates derived from the Monte Carlo samples and report standard errors for the N-independence plots to better quantify reliability. revision: yes
Circularity Check
No circularity: acceleration via standard transforms and empirical benchmarks
full rationale
The derivation relies on the fast Walsh-Hadamard transform and an exact Pauli partition to achieve the stated O(N) per-sample cost reduction; this follows directly from linear-algebraic properties of the Hadamard matrix and Pauli-string enumeration without any self-referential fitting or redefinition. The Monte Carlo estimator with Clifford preconditioning is introduced as a standard variance-reduction technique, with N-independent sample counts reported only as an empirical observation on the specific T-doped Clifford ensemble. No load-bearing self-citations, ansatz smuggling, or uniqueness theorems imported from prior author work appear in the chain; the central claims remain independent of the target quantities by construction.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Pauli operators form a complete basis for N-qubit operators and admit exact partitioning compatible with the Walsh-Hadamard transform.
- domain assumption Clifford circuits can be used as preconditioners that do not bias the Monte Carlo estimator for nonstabilizerness.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
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).
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
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.
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 5 Pith papers
-
Nonstabilizerness Mpemba Effects
In U(1)-symmetric random circuits, initial states with lower stabilizer Rényi entropy generate nonstabilizerness faster than those with higher entropy, with the effect also depending on spatial charge structure and ex...
-
Non-stabilizerness and U(1) symmetry in chaotic many-body quantum systems
Exact results show U(1) symmetry substantially suppresses non-stabilizerness in random states, with different leading scaling from entanglement near zero charge density.
-
Computing quantum magic of state vectors
Efficient algorithms compute stabilizer Rényi entropy and mana for quantum states from vectors at O(N d^{2N}) cost using fast Hadamard transform, with open-source implementation.
-
Quantum Magic in early FTQC: From Diagonal Clifford Hierarchy No-Go Theorems to Architecture Design Blueprints
No-go theorems prove hierarchy level and state-independent sequences cannot maximize operational magic in early FTQC, requiring state-aware differentiable optimization and nonlinear phases for scalable magic generation.
-
Non-Local Magic Resources for Fermionic Gaussian States
Closed-form formula computes non-local magic for fermionic Gaussian states from two-point correlations in polynomial time.
Reference graph
Works this paper leans on
-
[1]
P. W. Shor, Polynomial-Time Algorithms for Prime Fac- torization and Discrete Logarithms on a Quantum Com- puter, SIAM J. Comput.26, 1484 (1997)
work page 1997
-
[2]
A. M. Childs and W. Van Dam, Quantum algorithms for algebraic problems, Rev. Mod. Phys.82, 1 (2010)
work page 2010
- [3]
- [4]
-
[5]
J. I. Cirac and P. Zoller, Goals and opportunities in quan- tum simulation, Nature Phys8, 264 (2012)
work page 2012
-
[6]
M. A. Nielsen and I. L. Chuang,Quantum Computa- tion and Quantum Information: 10th Anniversary Edi- tion(Cambridge University Press, 2010)
work page 2010
-
[7]
S. Aaronson and D. Gottesman, Improved simulation of stabilizer circuits, Phys. Rev. A70, 052328 (2004)
work page 2004
-
[8]
The Heisenberg Representation of Quantum Computers
D. Gottesman, The Heisenberg Representation of Quan- tum Computers (1998), arXiv:quant-ph/9807006
work page internal anchor Pith review Pith/arXiv arXiv 1998
-
[9]
S. Bravyi and A. Kitaev, Universal quantum computa- tion with ideal Clifford gates and noisy ancillas, Phys. Rev. A71, 022316 (2005)
work page 2005
- [10]
- [11]
-
[12]
E. Chitambar and G. Gour, Quantum resource theories, Rev. Mod. Phys.91, 025001 (2019)
work page 2019
-
[13]
X. Wang, M. M. Wilde, and Y. Su, Efficiently Com- putable Bounds for Magic State Distillation, Phys. Rev. Lett.124, 090505 (2020)
work page 2020
- [14]
-
[15]
G. E. Fux, E. Tirrito, M. Dalmonte, and R. Fazio, Entan- glement – nonstabilizerness separation in hybrid quan- tum circuits, Phys. Rev. Research6, L042030 (2024)
work page 2024
- [16]
-
[17]
Y. Zhang and Y. Gu, Quantum magic dynamics in ran- dom circuits (2024), arXiv:2410.21128 [quant-ph]
-
[18]
X. Turkeshi, E. Tirrito, and P. Sierant, Magic spread- ing in random quantum circuits, Nat Commun16, 2575 (2025)
work page 2025
- [19]
-
[20]
J. Odavi´ c, T. Haug, G. Torre, A. Hamma, F. Franchini, and S. M. Giampaolo, Complexity of frustration: A new source of non-local non-stabilizerness, SciPost Phys.15, 131 (2023)
work page 2023
-
[21]
G. Lami, T. Haug, and J. De Nardis, Quantum State Designs with Clifford-Enhanced Matrix Product States, PRX Quantum6, 010345 (2025)
work page 2025
-
[22]
P. R. N. Falc˜ ao, P. Sierant, J. Zakrzewski, and E. Tir- rito, Nonstabilizerness Dynamics in Many-Body Local- ized Systems, Phys. Rev. Lett.135, 240404 (2025)
work page 2025
-
[23]
Growth and spreading of quantum resources under random circuit dynamics,
S. Aditya, X. Turkeshi, and P. Sierant, Growth and spreading of quantum resources under random circuit dy- namics (2025), arXiv:2512.14827 [quant-ph]
-
[24]
Mpemba effects in quantum complexity,
S. Aditya, A. Summer, P. Sierant, and X. Turkeshi, Mpemba Effects in Quantum Complexity (2025), arXiv:2509.22176 [quant-ph]
-
[25]
N. Dowling, P. Kos, and X. Turkeshi, Magic Resources of the Heisenberg Picture, Phys. Rev. Lett.135, 050401 (2025)
work page 2025
- [26]
-
[27]
S. Maity and R. Hamazaki, Local spreading of stabi- lizer R´ enyi entropy in a brickwork random Clifford circuit (2025), arXiv:2511.07769 [quant-ph]
-
[28]
D. Szombathy, A. Valli, C. P. Moca, L. Farkas, and G. Zar´ and, Independent stabilizer R´ enyi entropy and entanglement fluctuations in random unitary circuits (2025), arXiv:2501.11489 [quant-ph]
-
[29]
D. Sticlet, B. D´ ora, D. Szombathy, G. Zar´ and, and C. P. Moca, Nonstabilizerness in open XXZ spin chains: Uni- versal scaling and dynamics, Phys. Rev. Research7, 043130 (2025)
work page 2025
-
[30]
G. Passarelli, A. Russomanno, and P. Lucignano, Non- stabilizerness of a Boundary Time Crystal, Phys. Rev. A 111, 062417 (2025), arXiv:2503.05243 [quant-ph]
- [31]
-
[32]
K. Aziz, H. Pan, M. J. Gullans, and J. H. Pixley, Classical Simulations of Low Magic Quantum Dynamics (2025), arXiv:2508.20252 [quant-ph]
work page internal anchor Pith review arXiv 2025
-
[33]
P. S. Tarabunga and E. Tirrito, Magic transition in measurement-only circuits, npj Quantum Inf11, 166 (2025)
work page 2025
-
[34]
D. A. Korbany, M. J. Gullans, and L. Piroli, Long-Range Nonstabilizerness and Phases of Matter, Phys. Rev. Lett. 135, 160404 (2025)
work page 2025
- [35]
- [36]
-
[37]
S. F. E. Oliviero, L. Leone, and A. Hamma, Magic-state resource theory for the ground state of the transverse- field Ising model, Phys. Rev. A106, 042426 (2022)
work page 2022
- [38]
-
[39]
P. R. N. Falc˜ ao, P. S. Tarabunga, M. Frau, E. Tirrito, J. Zakrzewski, and M. Dalmonte, Nonstabilizerness in U(1) lattice gauge theory, Phys. Rev. B111, L081102 (2025)
work page 2025
-
[40]
C. D. White, C. Cao, and B. Swingle, Conformal field theories are magical, Phys. Rev. B103, 075145 (2021)
work page 2021
-
[41]
D. Qian and J. Wang, Quantum nonlocal nonstabilizer- ness, Phys. Rev. A111, 052443 (2025)
work page 2025
-
[42]
M. Frau, P. S. Tarabunga, M. Collura, E. Tirrito, and M. Dalmonte, Stabilizer disentangling of conformal field theories, SciPost Phys.18, 165 (2025)
work page 2025
-
[43]
M. Hoshino, M. Oshikawa, and Y. Ashida, Stabilizer R´ enyi Entropy and Conformal Field Theory (2025), arXiv:2503.13599 [quant-ph]
- [44]
-
[45]
X. Turkeshi, A. Dymarsky, and P. Sierant, Pauli spec- trum and nonstabilizerness of typical quantum many- body states, Phys. Rev. B111, 054301 (2025)
work page 2025
- [46]
- [47]
-
[48]
Non-stabilizerness of sachdev-ye-kitaev model (2025)
S. Bera and M. Schir` o, Non-Stabilizerness of Sachdev- Ye-Kitaev Model, SciPost Phys.19, 159 (2025), arXiv:2502.01582 [quant-ph]
-
[49]
E. Tirrito, X. Turkeshi, and P. Sierant, Anticoncentration and Nonstabilizerness Spreading under Ergodic Quan- tum Dynamics, Phys. Rev. Lett.135, 220401 (2025)
work page 2025
-
[50]
J. Odavi´ c, M. Viscardi, and A. Hamma, Stabilizer en- tropy in nonintegrable quantum evolutions, Phys. Rev. B112, 104301 (2025)
work page 2025
-
[51]
R. Nandkishore and D. A. Huse, Many-Body Localization and Thermalization in Quantum Statistical Mechanics, Annu. Rev. Condens. Matter Phys.6, 15 (2015)
work page 2015
- [53]
-
[54]
M. Beverland, E. Campbell, M. Howard, and V. Kli- uchnikov, Lower bounds on the non-Clifford resources for quantum computations, Quantum Sci. Technol.5, 035009 (2020)
work page 2020
-
[55]
J. Jiang and X. Wang, Lower Bound for the$T$Count Via Unitary Stabilizer Nullity, Phys. Rev. Appl.19, 034052 (2023)
work page 2023
-
[56]
T. Haug and M. Kim, Scalable Measures of Magic Resource for Quantum Computers, PRX Quantum4, 010301 (2023)
work page 2023
-
[57]
G. Lami and M. Collura, Unveiling the Stabilizer Group of a Matrix Product State, Phys. Rev. Lett.133, 010602 (2024)
work page 2024
-
[58]
L. Chen, R. J. Garcia, K. Bu, and A. Jaffe, Magic of random matrix product states, Phys. Rev. B109, 174207 (2024)
work page 2024
-
[59]
T. Haug and L. Piroli, Quantifying nonstabilizerness of matrix product states, Phys. Rev. B107, 035148 (2023)
work page 2023
-
[60]
P. S. Tarabunga, E. Tirrito, M. C. Ba˜ nuls, and M. Dal- monte, Nonstabilizerness via Matrix Product States in the Pauli Basis, Phys. Rev. Lett.133, 010601 (2024)
work page 2024
-
[61]
G. Lami and M. Collura, Nonstabilizerness via Perfect Pauli Sampling of Matrix Product States, Phys. Rev. Lett.131, 180401 (2023)
work page 2023
-
[62]
The quantum magic of fermionic gaussian states
M. Collura, J. D. Nardis, V. Alba, and G. Lami, The non-stabilizerness of fermionic Gaussian states (2025), arXiv:2412.05367 [quant-ph]
-
[63]
C. Wang, Z.-C. Yang, T. Zhou, and X. Chen, Magic transition in monitored free fermion dynamics (2025), arXiv:2507.10688 [quant-ph]
-
[64]
Y.-M. Ding, Z. Wang, and Z. Yan, Evaluating Many- Body Stabilizer R´ enyi Entropy by Sampling Reduced Pauli Strings: Singularities, Volume Law, and Nonlocal Magic, PRX Quantum6, 030328 (2025)
work page 2025
- [65]
-
[66]
Shanks, Computation of the Fast Walsh-Fourier Trans- form, IEEE Trans
J. Shanks, Computation of the Fast Walsh-Fourier Trans- form, IEEE Trans. Comput.C-18, 457 (1969)
work page 1969
-
[67]
T. Haug and L. Piroli, Stabilizer entropies and nonstabi- lizerness monotones, Quantum7, 1092 (2023)
work page 2023
-
[68]
L. Leone and L. Bittel, Stabilizer entropies are mono- tones for magic-state resource theory, Phys. Rev. A110, L040403 (2024)
work page 2024
-
[69]
J. L. Walsh, A Closed Set of Normal Orthogonal Func- tions, Am. J. Math.45, 5 (1923), 2387224
work page 1923
-
[70]
See the Supplemental Material for: (i) the analytic large- NC behavior ofM 2 forT ⊗N |ψ(NC)⟩; (ii) the special treatment of the identity Pauli string in Monte Carlo sampling; and (iii) the extension of the FWHT frame- work to subsystem (mixed-state) stabilizer R´ enyi en- tropies
-
[71]
A. D. C´ orcoles, J. M. Gambetta, J. M. Chow, J. A. Smolin, M. Ware, J. Strand, B. L. T. Plourde, and M. Steffen, Process verification of two-qubit quantum gates by randomized benchmarking, Phys. Rev. A87, 030301 (2013)
work page 2013
- [72]
- [73]
-
[74]
N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, Equation of State Calcula- tions by Fast Computing Machines, J. Chem. Phys.21, 1087 (1953)
work page 1953
-
[75]
W. K. Hastings, Monte Carlo sampling methods using Markov chains and their applications, Biometrika57, 97 (1970)
work page 1970
-
[76]
S. Zhou, Z. Yang, A. Hamma, and C. Chamon, Single T gate in a Clifford circuit drives transition to universal entanglement spectrum statistics, SciPost Phys.9, 087 (2020)
work page 2020
-
[77]
D. Szombathy, A. Valli, C. P. Moca, J. Asb´ oth, L. Farkas, T. Rakovszky, and G. Zar´ and, Spectral Properties Versus Magic Generation in$T$-doped Random Clifford Cir- cuits (2025), arXiv:2412.15912 [quant-ph]
- [78]
- [79]
-
[80]
T. Rakovszky, F. Pollmann, and C. W. Von Keyserlingk, Diffusive Hydrodynamics of Out-of-Time-Ordered Corre- lators with Charge Conservation, Phys. Rev. X8, 031058 (2018)
work page 2018
-
[81]
Sub-ballistic growth of R\'enyi entropies due to diffusion
T. Rakovszky, F. Pollmann, and C. W. von Keyserlingk, Sub-ballistic growth of R´ enyi entropies due to diffusion, Phys. Rev. Lett.122, 250602 (2019), arXiv:1901.10502 [cond-mat]
work page internal anchor Pith review Pith/arXiv arXiv 2019
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.