arxiv: 2604.15043 · v1 · submitted 2026-04-16 · 🪐 quant-ph

Recognition: unknown

SyQMA: A memory-efficient, symbolic and exact universal simulator for quantum error correction

George Umbrarescu , David Amaro

Authors on Pith no claims yet

Pith reviewed 2026-05-10 11:35 UTC · model grok-4.3

classification 🪐 quant-ph

keywords quantum simulationquantum error correctionstabilizer simulationsymbolic computationfault tolerancePauli noisemaximum likelihood decodingmagic state preparation

0 comments

The pith

A simulator computes exact symbolic expectation values and measurement probabilities for universal quantum circuits with Pauli noise using only polynomial memory.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper presents a method to simulate universal quantum circuits that include incoherent Pauli noise by extending the stabilizer formalism. Non-Clifford rotations and noise channels are represented compactly using auxiliary qubits along with a modified trace operation. This representation keeps the memory requirement polynomial in the number of qubits. Expectation values and probabilities are obtained as exact symbolic functions of parameters such as rotation angles, measurement outcomes, and noise rates. The approach supports sampling, maximum-likelihood decoding, and exact verification of fault-tolerant protocols without approximations.

Core claim

By representing each non-Clifford Pauli rotation and incoherent Pauli channel via auxiliary qubits and employing a modified trace, the simulator maintains a polynomial-memory description of the quantum state. Expectation values and measurement probabilities are then computed exactly as symbolic expressions in the circuit parameters, with the computation time being exponential only in the number of non-Clifford elements and deterministic measurements.

What carries the argument

Auxiliary-qubit encoding of non-Clifford Pauli rotations and incoherent Pauli channels together with a modified trace.

If this is right

Exact symbolic logical error rates can be derived for quantum error correction circuits.
Maximum-likelihood decoding at the circuit level becomes feasible for noisy quantum programs.
Fault distances of fault-tolerant stabilizer and magic state preparation protocols can be verified exactly.
Disjoint error probabilities of multi-qubit Pauli channels can be converted to independent probabilities for detector error models.
Measurement outcomes can be sampled to simulate adaptive quantum circuits.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

This method could facilitate the design of better error correction codes by allowing parameter optimization through symbolic expressions.
It may connect to tensor network methods for quantum simulation by providing an alternative exact approach for certain noisy circuits.
Further development might reduce the exponential time factor by incorporating approximation techniques for larger systems.
Comparison with numerical simulators on small circuits could validate the symbolic outputs against known results.

Load-bearing premise

The auxiliary-qubit representation of non-Clifford operations and noise, along with the modified trace, produces exact results without leading to exponential complexity in the symbolic expressions themselves.

What would settle it

Running the simulator on a small circuit with one non-Clifford rotation and known Pauli noise, then comparing the symbolic output for a logical error rate against an independent analytical calculation; any discrepancy would falsify the exactness claim.

Figures

Figures reproduced from arXiv: 2604.15043 by David Amaro, George Umbrarescu.

**Figure 2.** Figure 2: Fault-tolerant preparation of |¯0⟩ in the Steane code. (a) Circuit for the FT state preparation of the logical state |¯0⟩ in the Steane code. (b) Discard rate and logical error rates under the three decoding strategies: uncorrected (Eq. (35)), corrected (Eq. (38)), and postselected (Eq. (41)). The uncorrected LER is 0.8 p 1 + 11.8 p 2 + O(p 3 ) and the postselected LER is 1.63 p 3 + 42.98 p 4 + O(p 5 ). Th… view at source ↗

**Figure 3.** Figure 3: Fault-tolerant (FT) magic state preparation of [PITH_FULL_IMAGE:figures/full_fig_p022_3.png] view at source ↗

**Figure 4.** Figure 4: Expectation values (EVs) of the logical operator [PITH_FULL_IMAGE:figures/full_fig_p023_4.png] view at source ↗

**Figure 5.** Figure 5: Fault-tolerant preparation of Clifford states in the [PITH_FULL_IMAGE:figures/full_fig_p025_5.png] view at source ↗

**Figure 6.** Figure 6: Fault-tolerant magic-state preparation of [PITH_FULL_IMAGE:figures/full_fig_p026_6.png] view at source ↗

**Figure 7.** Figure 7: Comparison of memory and runtime between the computation mode with the largest batch size and 1 [PITH_FULL_IMAGE:figures/full_fig_p026_7.png] view at source ↗

**Figure 8.** Figure 8: Fault-tolerant preparation of |¯0⟩ in the J17, 1, 5K 2D colour code. (a) Flagged preparation and syndrome extraction circuit adapted from [FA25], where the accepted branch postselects on all measured ancilla and syndrome outcomes being +1. (b) Exact discard rate and Z¯ logical error rates under the uncorrected (Eq. (35)), corrected (Eq. (38)), and postselected (Eq. (41)) decoding strategies. 27 [PITH_FULL… view at source ↗

read the original abstract

The classical simulation of universal quantum circuits is crucial both fundamentally and practically for quantum computation. We propose SyQMA, a simulator with several convenient features, particularly suited for quantum error correction (QEC). SyQMA simulates universal quantum circuits with incoherent Pauli noise and computes exact expectation values and measurement probabilities as symbolic functions of circuit parameters: rotation angles, measurement outcomes, and noise rates. This simulator can sample measurement outcomes, enabling the simulation of dynamic quantum programs where circuit composition depends on prior measurement outputs. For QEC, it performs circuit-level maximum-likelihood decoding, provides exact symbolic expressions for logical error rates, and verifies the fault distance of fault-tolerant (FT) stabiliser and magic state preparation protocols. These features are enabled by an intuitive extension of stabiliser simulators, where each non-Clifford Pauli rotation and incoherent Pauli channel is compactly represented via auxiliary qubits and a modified trace. Representing the state requires only polynomial memory and time, while computing expectation values and measurement probabilities takes exponential time in the number of non-Clifford rotations and deterministic measurements, but only polynomial memory. The FT preparation of stabiliser and magic states, including the first stage of magic state cultivation, is analysed without approximations. We also exactly convert the disjoint error probabilities of a general multi-qubit Pauli channel to independent ones, a key step for creating and sampling from detector error models. The code is publicly available and open-source.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

SyQMA adds auxiliary qubits to stabilizer tableaux to get exact symbolic probabilities for noisy universal circuits, which is a concrete extension worth checking, though the size of the resulting expressions may limit its reach.

read the letter

SyQMA extends the stabilizer formalism by encoding each non-Clifford Pauli rotation and each incoherent Pauli channel with an auxiliary qubit, then pulls out symbolic expectation values and probabilities through a modified trace. The state stays a tableau on n plus O(k) qubits, so memory stays polynomial while the final symbolic outputs take exponential time in the number of such operations. They also add sampling for adaptive circuits and a conversion from general multi-qubit Pauli channels to independent ones for detector models. The open-source code is a clear plus for anyone who wants to test the claims directly. They apply the method to circuit-level maximum-likelihood decoding and to exact fault-distance checks on stabilizer and magic-state preparation protocols, including the first stage of cultivation, without falling back on approximations. That part is useful for QEC work where closed-form logical error rates matter. The soft spot is the practical size of the symbolic expressions themselves. Even if the tableau stays small, a probability with k non-Cliffords can in principle expand to a sum with 2^k terms unless the implementation keeps everything factored or applies aggressive simplification. The paper acknowledges the exponential runtime but does not spell out the output representation or simplification steps in enough detail to judge how far the method scales beyond a few parameters. For circuits with more than a handful of non-Cliffords the exact expressions could become unwieldy even if memory does not. This is aimed at QEC researchers who need precise symbolic results on modest-depth circuits rather than Monte Carlo sampling. A reader already working with stabilizer simulators or fault-tolerant protocol verification will find the concrete examples and the public code worth looking at. I would send it to referees. The algorithmic idea is new enough and the QEC applications are specific enough that the claims can be checked properly in review.

Referee Report

1 major / 1 minor

Summary. The manuscript introduces SyQMA, an extension of stabilizer simulation that represents each non-Clifford Pauli rotation and incoherent Pauli channel using one auxiliary qubit together with a modified trace. This construction is claimed to simulate universal circuits with Pauli noise while returning exact symbolic expressions for expectation values and measurement probabilities (as functions of rotation angles, noise rates, and measurement outcomes) using only polynomial memory for the state tableau; computation time is exponential in the number of non-Clifford operations. The simulator supports sampling for dynamic circuits and is applied to quantum error correction tasks, including circuit-level maximum-likelihood decoding, exact symbolic logical error rates, and fault-distance verification for stabilizer and magic-state preparation protocols. An exact conversion from disjoint to independent error probabilities for multi-qubit Pauli channels is also presented. The code is stated to be publicly available.

Significance. If the auxiliary-qubit representation and modified trace preserve exactness for arbitrary numbers of non-Clifford operations without hidden exponential growth in the size or complexity of the final symbolic expressions, the work would supply a useful exact tool for symbolic analysis of QEC protocols where only numerical or approximate methods are currently available. The open-source implementation and the exact channel-conversion step are concrete strengths that support reproducibility and downstream use in detector-error-model construction.

major comments (1)

[Abstract] Abstract: The central claim that exact symbolic functions are obtained with only polynomial memory is load-bearing for all QEC applications. The construction adds one auxiliary qubit per non-Clifford rotation or channel, after which the modified trace encodes a sum/product over those auxiliaries. The manuscript does not specify the representation of the resulting symbolic expression (e.g., whether it is kept in a compact factored, circuit-encoded, or sparse-polynomial form, or whether automatic simplification is performed). Without such a representation, the output size grows exponentially in the number of auxiliaries, rendering the 'exact symbolic' results impractical for circuits with more than a few non-Clifford gates even though the intermediate tableau remains polynomial-sized. This must be clarified with an explicit statement of the output data structure and any complexity bounds.

minor comments (1)

[Abstract] The abstract states that 'representing the state requires only polynomial memory and time' while 'computing expectation values... takes exponential time... but only polynomial memory.' The distinction between state representation and output expression size should be made explicit in the main text to avoid reader confusion.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for identifying this important point of clarification. We address the comment below and will revise the manuscript to incorporate the requested details.

read point-by-point responses

Referee: [Abstract] Abstract: The central claim that exact symbolic functions are obtained with only polynomial memory is load-bearing for all QEC applications. The construction adds one auxiliary qubit per non-Clifford rotation or channel, after which the modified trace encodes a sum/product over those auxiliaries. The manuscript does not specify the representation of the resulting symbolic expression (e.g., whether it is kept in a compact factored, circuit-encoded, or sparse-polynomial form, or whether automatic simplification is performed). Without such a representation, the output size grows exponentially in the number of auxiliaries, rendering the 'exact symbolic' results impractical for circuits with more than a few non-Clifford gates even though the intermediate tableau remains polynomial-sized. This must be clarified with an explicit statement of the output data structure and any complex

Authors: We agree that an explicit description of the output representation is needed. In SyQMA the modified trace evaluates the symbolic expectation or probability by enumerating all 2^k assignments to the auxiliary qubits (k = number of non-Clifford rotations/channels). Each assignment contributes a single monomial whose coefficient is the stabilizer expectation value (±1 or 0) and whose variables are the circuit parameters (sin(θ/2), cos(θ/2), p, 1-p, etc.). The resulting expression is therefore stored as a sparse list of at most 2^k monomials; no dense polynomial or automatic simplification is performed. This representation uses only polynomial memory during tableau evolution and per-configuration evaluation, while the final list size and evaluation time are O(2^k poly(n)). We will add a dedicated paragraph (and pseudocode) in the Methods section describing this data structure, its worst-case size, and how it is emitted for downstream use (e.g., logical-error-rate polynomials). The abstract will also be updated to reference the exponential dependence on k for the symbolic output. revision: yes

Circularity Check

0 steps flagged

No significant circularity; direct algorithmic extension of stabilizer simulation

full rationale

The paper presents SyQMA as an explicit construction: each non-Clifford Pauli rotation or incoherent Pauli channel is represented by adding one auxiliary qubit, the state remains a stabilizer tableau on n + O(k) qubits, and expectation values are obtained via a modified trace. This is a direct, self-contained algorithmic extension with explicitly stated resource bounds (polynomial memory, exponential time in k). No equations reduce a claimed prediction to a fitted parameter or to a self-referential definition. No load-bearing self-citations or uniqueness theorems from prior author work are invoked to justify the core method. The derivation chain is therefore independent of its own outputs and does not collapse by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The method rests on the standard stabilizer formalism plus a new auxiliary-qubit encoding whose correctness is asserted but not derived in the abstract.

axioms (1)

standard math Stabilizer formalism and Pauli channel representation are sufficient to capture the relevant quantum operations
The simulator is described as an intuitive extension of stabilizer simulators.

invented entities (1)

Auxiliary qubits with modified trace for non-Clifford rotations and Pauli channels no independent evidence
purpose: Compact representation that enables polynomial memory and exact symbolic computation
Introduced as the core technical device allowing the claimed features.

pith-pipeline@v0.9.0 · 5552 in / 1176 out tokens · 51908 ms · 2026-05-10T11:35:23.525341+00:00 · methodology

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Analytical and Compressed Simulation of Noisy Stabilizer Circuits
quant-ph 2026-04 unverdicted novelty 6.0

Closed-form expressions and circuit compression enable efficient strong and weak simulation of noisy stabilizer circuits with non-deterministic measurements.

Reference graph

Works this paper leans on

3 extracted references · cited by 1 Pith paper

[1]

as Pauli operators they satisfyF 2 f =I f,C 2 r =O 2 r =S 2 r =I r, andC rOr =iS r,
[2]

their generalised trace isTr f(Ff) = 2ϵf andTrf(If) = 2for all flip qubits,
[3]

We choose the definition of the generalised trace so that the trace of the identity isTr(I) = 2for every qubit

andTr r(Cr) = 2 cos(θr),Trr(Or) = 0,Trr(Sr) = 2 sin(θr),Trr(Ir) = 2for all rotation qubits. We choose the definition of the generalised trace so that the trace of the identity isTr(I) = 2for every qubit. The tracetris then replaced by the generalised traceTr, and we simulate the evolution of the joint stateρon all qubits inNwith SyQMA’s framework. The (ge...