Classical Simulations of Low Magic Quantum Dynamics
Pith reviewed 2026-05-22 12:20 UTC · model grok-4.3
The pith
Classical simulation algorithms efficiently simulate adaptive quantum circuits with low magic states by exploiting frequent Pauli measurements.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We develop classical simulation algorithms for adaptive quantum circuits that produce states with low levels of magic. These algorithms are particularly well-suited to circuits with high rates of Pauli measurements, such as those encountered in quantum error correction and monitored quantum circuits. The measurements serve to limit the buildup of magic induced by non-Clifford operations arising from generic noise processes or unitary gates. Our algorithms also allow a systematic truncation procedure to achieve approximate simulation. To benchmark our approach, we study the dynamics of all-to-all monitored quantum circuits with a sub-extensive rate of T-gates per unit of circuit depth, where
What carries the argument
Pauli measurements that limit magic accumulation from non-Clifford operations in adaptive circuits, enabling polynomial-cost classical simulation.
Load-bearing premise
Frequent Pauli measurements are assumed to sufficiently suppress the accumulation of magic from non-Clifford operations or noise, maintaining polynomial simulation cost.
What would settle it
A calculation showing exponential growth in effective magic or simulation runtime with system size despite high rates of Pauli measurements would invalidate the efficiency claims.
Figures
read the original abstract
We develop classical simulation algorithms for adaptive quantum circuits that produce states with low levels of ``magic'' (i.e., non-stabilizerness). These algorithms are particularly well-suited to circuits with high rates of Pauli measurements, such as those encountered in quantum error correction and monitored quantum circuits. The measurements serve to limit the buildup of magic induced by non-Clifford operations arising from generic noise processes or unitary gates, respectively. Our algorithms also allow a systematic truncation procedure to achieve approximate simulation. To benchmark our approach, we study the dynamics of all-to-all monitored quantum circuits with a sub-extensive rate of T-gates per unit of circuit depth, where we can simulate previously inaccessible system sizes and depths. We characterize measurement-induced phase transitions in the output wavefunction, including in the entanglement, purification, and magic. We outline the utility of our algorithms to simulate dynamics with low magic and high entanglement, complementary to the leading matrix-product state approaches.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops classical simulation algorithms for adaptive quantum circuits producing low-magic states, leveraging high rates of Pauli measurements to suppress magic buildup from non-Clifford gates or noise. It introduces a systematic truncation procedure for approximate simulation and benchmarks the approach on all-to-all monitored circuits with sub-extensive T-gates per depth, enabling access to larger system sizes and depths. The work characterizes measurement-induced phase transitions in entanglement, purification, and magic, positioning the method as complementary to matrix-product states for low-magic high-entanglement regimes.
Significance. If the central claims hold, the algorithms offer a practical tool for classically simulating previously inaccessible regimes of monitored quantum circuits and error-corrected dynamics, particularly where entanglement is high but magic is controlled by measurements. The benchmarking results on phase transitions provide concrete data on larger scales, and the complementarity to MPS methods addresses a genuine gap in simulation techniques for low-magic states.
major comments (2)
- [§3] §3 (algorithm description): The truncation procedure is introduced as systematic, but no explicit error bound or convergence analysis is given relating the truncation threshold (a free parameter) to the approximation error in the output state or observables; this is load-bearing for the claim of controlled approximate simulation at polynomial cost.
- [§5] §5 (benchmarking and phase transitions): While sub-extensive T-gates per depth allow simulation of larger sizes, the efficiency claim rests on the assumption that Pauli measurements keep effective magic (e.g., via a monotone such as mana or stabilizer extent) from growing exponentially; no derivation or bound is provided showing magic scaling with measurement rate, depth, and noise strength for generic adaptive circuits, as required to guarantee polynomial cost beyond the specific benchmark.
minor comments (2)
- Figure captions for the phase-transition plots could more explicitly state the system sizes, depths, and number of samples used to generate the data points.
- The abstract mentions 'previously inaccessible system sizes and depths' but the main text would benefit from a direct comparison table to prior simulation limits in the literature.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the constructive comments, which help clarify the scope and rigor of our claims. We respond to each major comment below, indicating where revisions have been made to the manuscript.
read point-by-point responses
-
Referee: [§3] §3 (algorithm description): The truncation procedure is introduced as systematic, but no explicit error bound or convergence analysis is given relating the truncation threshold (a free parameter) to the approximation error in the output state or observables; this is load-bearing for the claim of controlled approximate simulation at polynomial cost.
Authors: We agree that an explicit error bound strengthens the presentation of the truncation procedure. In the revised manuscript we have added a derivation in §3 relating the truncation threshold to the approximation error. Specifically, when truncating the stabilizer decomposition by discarding components below the threshold, the error in the expectation value of any Pauli observable is bounded by twice the total discarded weight (i.e., the sum of the stabilizer extents of the discarded terms). This bound is independent of system size for normalized states and directly controls the observable error. We also include a brief convergence statement: as the threshold tends to zero the approximate state converges to the exact state in the 1-norm for the relevant observables. These additions make the controlled-approximation claim explicit while preserving the polynomial-cost scaling whenever the retained magic remains sub-exponential. revision: yes
-
Referee: [§5] §5 (benchmarking and phase transitions): While sub-extensive T-gates per depth allow simulation of larger sizes, the efficiency claim rests on the assumption that Pauli measurements keep effective magic (e.g., via a monotone such as mana or stabilizer extent) from growing exponentially; no derivation or bound is provided showing magic scaling with measurement rate, depth, and noise strength for generic adaptive circuits, as required to guarantee polynomial cost beyond the specific benchmark.
Authors: Our efficiency claim is made for the concrete regime of all-to-all monitored circuits with a sub-extensive density of T gates per layer, where high-rate Pauli measurements are shown numerically to keep the effective magic (mana and stabilizer extent) from growing exponentially with depth. We have expanded the discussion in §5 to include additional scaling plots of mana versus depth and measurement probability, together with references to existing results on measurement-induced suppression of magic. A fully general analytical bound that holds for arbitrary adaptive circuits, arbitrary measurement rates, and arbitrary noise strengths would require model-specific assumptions on the adaptivity and the distribution of measurement outcomes; such a derivation lies outside the scope of the present work. The manuscript therefore positions the algorithm as complementary to MPS methods precisely in the low-magic, high-entanglement regime that our benchmarks access, rather than claiming universality for all adaptive circuits. revision: partial
- A general derivation or bound on magic scaling with measurement rate, depth, and noise strength that holds for arbitrary adaptive circuits (beyond the specific all-to-all monitored ensemble with sub-extensive T gates) is not provided, as it would require additional model-dependent assumptions not addressed in the current manuscript.
Circularity Check
No significant circularity; derivation is self-contained
full rationale
The paper develops classical simulation algorithms for low-magic adaptive circuits by extending the stabilizer formalism with magic measures and a truncation procedure. The central approach relies on the physical assumption that frequent Pauli measurements limit magic accumulation from non-Clifford gates or noise, which is then benchmarked on all-to-all monitored circuits with sub-extensive T-gates. This assumption is tested empirically rather than defined into the result by construction, and the characterizations of entanglement, purification, and magic phase transitions follow from the low-magic representation without reducing to fitted inputs or self-referential definitions. No load-bearing step equates a prediction to its own input via equation or self-citation chain; the work is complementary to matrix-product states and builds on independently established concepts.
Axiom & Free-Parameter Ledger
free parameters (1)
- truncation threshold
axioms (2)
- domain assumption Pauli measurements limit the buildup of magic induced by non-Clifford operations
- standard math Stabilizer formalism plus magic measures allow efficient classical representation when magic is low
Forward citations
Cited by 1 Pith paper
-
Exponentially Accelerated Sampling of Pauli Strings for Nonstabilizerness
A sampling method combining fast Walsh-Hadamard transform and Clifford-preconditioned Monte Carlo reduces Pauli-string sampling cost from O(2^N) to O(N) with sample count independent of N for stabilizer Rényi entropie...
Reference graph
Works this paper leans on
-
[1]
U. Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Annals of Physics 326, 96 (2011)
work page 2011
-
[2]
G. Vidal, Efficient simulation of one-dimensional quan- tum many-body systems, Physical review letters93, 040502 (2004)
work page 2004
-
[3]
A. G. Fowler, M. Mariantoni, J. M. Martinis, and A. N. Cleland, Surface codes: Towards practical large-scale quantum computation, Physical Review A86, 032324 (2012)
work page 2012
-
[4]
B. M. Terhal, Quantum error correction for quantum memories, Reviews of Modern Physics87, 307 (2015)
work page 2015
-
[5]
Y. Li, X. Chen, and M. P. Fisher, Quantum zeno ef- fect and the many-body entanglement transition, Physi- cal Review B98, 205136 (2018)
work page 2018
-
[6]
B. Skinner, J. Ruhman, and A. Nahum, Measurement- induced phase transitions in the dynamics of entangle- ment, Phys. Rev. X9, 031009 (2019)
work page 2019
-
[7]
S. Choi, Y. Bao, and X.-L. Qi, Quantum information dynamics in multipartite monitoring of quantum circuits, Physical Review Letters125, 030505 (2020)
work page 2020
- [8]
-
[9]
The Heisenberg Representation of Quantum Computers
D. Gottesman, The heisenberg representation of quan- tum computers (1998), arXiv:quant-ph/9807006 [quant- ph]
work page internal anchor Pith review Pith/arXiv arXiv 1998
-
[11]
S. Bravyi and A. Kitaev, Universal quantum computa- tion with ideal clifford gates and noisy ancillas, Phys. Rev. A71, 022316 (2005)
work page 2005
- [12]
- [13]
-
[14]
R. Y. Huang, M. Newman, Z. Cai, S. Gharibian, and J. Yirka, Classical simulation of quantum circuits by low-rank pauli matrix decomposition, Quantum4, 310 (2020)
work page 2020
- [15]
- [16]
- [17]
- [18]
- [19]
-
[20]
Classical simulability of Clifford+T circuits with Clifford-augmented matrix product states
Z. Liu and B. K. Clark, Classical simulability of clif- ford+t circuits with clifford-augmented matrix prod- uct states, arXiv preprint arXiv:2412.17209 (2025), arXiv:2412.17209 [quant-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[21]
X. Qian, J. Huang, and M. Qin, Augmenting density ma- trix renormalization group with clifford circuits, Phys. Rev. Lett.133, 190402 (2024)
work page 2024
-
[22]
S. Masot-Llima and A. Garcia-Saez, Stabilizer tensor net- works: A universal quantum circuit simulator based on a complete basis of stabilizer states, Quantum5, 559 (2021)
work page 2021
-
[23]
Y. Li, X. Chen, and M. P. A. Fisher, Measurement- driven entanglement transition in hybrid quantum cir- cuits, Phys. Rev. B100, 134306 (2019)
work page 2019
- [24]
-
[25]
T. Haug and M. Kim, Scalable measures of magic re- source for quantum computers, PRX Quantum4, 010301 (2023)
work page 2023
-
[26]
L. Leone, S. F. Oliviero, and A. Hamma, Stabilizer r´ enyi entropy, Physical Review Letters128, 10.1103/phys- revlett.128.050402 (2022)
-
[27]
J. Bermejo-Vega, D. Hangleiter, M. Schwarz, R. Raussendorf, and J. Eisert, Architectures for quantum simulation showing a quantum speedup, Phys. Rev. X8, 021010 (2018)
work page 2018
-
[28]
Y. Wang, Y. Wang, Y.-A. Chen, W. Zhang, T. Zhang, J. Hu, W. Chen, Y. Gu, and Z.-W. Liu, Efficient fault- tolerant implementations of non-clifford gates with re- configurable atom arrays, npj Quantum Information10, 10.1038/s41534-024-00945-3 (2024)
-
[29]
M. Beverland, E. Campbell, M. Howard, and V. Kli- uchnikov, Lower bounds on the non-clifford resources for quantum computations, Quantum Science and Technol- ogy5, 035009 (2020)
work page 2020
-
[30]
M. Howard and E. Campbell, Application of a resource theory for magic states to fault-tolerant quantum com- puting, Physical Review Letters118, 090501 (2017)
work page 2017
-
[31]
J. R. Seddon and E. T. Campbell, Quantifying magic for multi-qubit operations, PRX Quantum2, 010345 (2021)
work page 2021
-
[32]
R. K. Huang, J. Barry, and S. T. Flammia, Efficient clas- sical simulation of quantum circuits with magic states, PRX Quantum3, 010345 (2022)
work page 2022
-
[33]
S. Vijay, Measurement-driven phase transition within a volume-law entangled phase (2020), arXiv:2005.03052 [quant-ph]
-
[34]
H. J. Kimble, The quantum internet, Nature453, 1023–1030 (2008)
work page 2008
-
[35]
C. Noel, P. Niroula, D. Zhu, A. Risinger, L. Egan, D. Biswas, M. Cetina, A. V. Gorshkov, M. J. Gullans, D. A. Huse, and C. Monroe, Measurement-induced quan- tum phases realized in a trapped-ion quantum computer, Nature Physics18, 760–764 (2022)
work page 2022
- [36]
-
[37]
This choice ensures that the circuit is in a nontrivial en- tangling regime: choosing the probability of a CZ gate arbitrarily asp u, ifp u is too small it drives the system into a trivial area-law entangled phase for all non-zero measurement rates
-
[38]
C.-M. Jian, Y.-Z. You, R. Vasseur, and A. W. W. Ludwig, Measurement-induced criticality in random quantum cir- cuits, Phys. Rev. B101, 104302 (2020)
work page 2020
-
[39]
I. Sfiligoi, D. C. Bradley, B. Holzman, P. Mhashilkar, S. Padhi, and F. Wurthwein, The pilot way to grid re- sources using glideinWMS, in2009 WRI World Congress on Computer Science and Information Engineering, 2, Vol. 2, pp. 428–432
- [40]
-
[41]
OSG, Open science data federation ()
-
[42]
K. M. R. Audenaert and M. B. Plenio, Entanglement on mixed stabilizer states: normal forms and reduction procedures, New Journal of Physics7, 170–170 (2005)
work page 2005
-
[43]
Learning stabilizer states by Bell sampling
A. Montanaro, Learning stabilizer states by bell sam- pling, arXiv preprint arXiv:1707.04012 (2017)
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[44]
C.-Y. Lai and H.-C. Cheng, Learning quantum circuits of some T gates, IEEE Trans. Inf. Theory68, 3951 (2022), 2106.12524
-
[45]
M. Van den Nest, J. Dehaene, and B. De Moor, An ef- ficient algorithm to recognize local Clifford equivalence of graph states, Phys. Rev. A70, 034302 (2004), quant- ph/0405023
-
[46]
S. Anders and H. J. Briegel, Fast simulation of stabilizer circuits using a graph-state representation, Phys. Rev. A 73, 022334 (2006)
work page 2006
-
[47]
M. J. Gullans and D. A. Huse, Private communication (2023). Appendix A: Algorithms In this section, we provide an algorithm to update the LRSD in Eq. (4) under Clifford operations. These are the Clifford unitaries and Pauli measurements. We also explain how to compute the Born probabilities of mea- surement outcomes, and update the state after the mea- su...
work page 2023
-
[48]
Clifford Evolution In this section, we describe how to update the LRSD of a quantum state as it evolves under Clifford unitaries and Pauli measurements, see Eqs. (14) and (15). Through- out, we write the density matrix in the form of Eq. (4) of the main text, |ψ⟩⟨ψ|= X l∈LLRSD λl σl ρS,(A1) where the coefficientsλ l are real, the operatorsσ l are logical ...
-
[49]
R´ enyi Entropy of LRSD We now present algorithms to compute the R´ enyi en- tanglement entropies using the representation of the state in Eq. (4). We present an algorithm for computing the reduced density matrix for an arbitrary state, from which we can compute the entanglement. We also present an algorithm which is valid only for [LS,S] = 0.These meth- ...
-
[50]
Bell sampling |ۧ𝜓 𝐻 𝐻 𝐻 𝐻 𝐻 |ۧ𝜓 𝑇 𝐻 1 𝑇 x t 1 0 0 FIG
Bell Sampling and Magic Computation Algorithms a. Bell sampling |ۧ𝜓 𝐻 𝐻 𝐻 𝐻 𝐻 |ۧ𝜓 𝑇 𝐻 1 𝑇 x t 1 0 0 FIG. 8. Bell sampling setup for the single-pair all-to-all model. The quantum circuit to measure the magicMof an L-qubit state|ψ⟩=K|+x⟩ ⊗N whereKis the Kraus operator of the circuit gates and measurements. We prepare two copies |ψ⟩ ⊗ |ψ⟩, and perform a meas...
-
[51]
Cluster We begin with the connectivity transition, which is a percolation transition from a connected state (when all qubits are entangled) to a disconnected state (when the state effectively breaks into multiple disentangled clus- ters). The quantum state factorizes as |ψ⟩= O i |ψCi ⟩,(B1) where eachC i denotes a disjoint subset of qubits, or clus- ter, ...
-
[52]
Entanglement Transition To study the entanglement transition, we need to gen- eralize the concept of entanglement entropy and the corresponding volume/area-law phase to the all-to-all model, as there is no clear notion of surface and bulk in this effective 0D system. Therefore, we generalize the entanglement entropy to the following max-min form [48] Smm,...
-
[53]
Purification Transition This section provides additional data supporting the purification transition results reported in Sec. V. We give more details on how we locate the purification transition in Section V. We present the data used to numerically extract the phase boundary of the purification transition in Fig. 5(a). We also explain the fitting procedur...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.