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:2603.23641

QuickQudits: A Framework for Efficient Simulation of Noisy Qudit Clifford Circuits via an Extended Stabilizer Tableau Formalism

Nina Brandl, Mykyta Cherniak, Johannes Kofler, Richard Kueng

quant-phSubmitted Mar 24, 2026• Announced Mar 26, 2026• 0 citations
arXivHTMLPDF

Abstract

We present a comprehensive and self-contained framework for the efficient classical simulation of Clifford circuits acting on $d$-dimensional qudits, including realistic Pauli/Weyl noise via stochastic simulation. Our approach uses the stabilizer tableau formalism for qudits of arbitrary dimension and tracks both stabilizer and destabilizer generators under Clifford updates. The classical simulation remains efficient with simple algebraic Clifford update rules over $\mathbb{Z}_d$. Computational basis measurements in prime dimensions are handled by a generalized Aaronson-Gottesman (CHP) procedure. In composite dimensions, $\mathbb{Z}_d$ is not a field and the standard tableau reduction fails, so we employ an exact Smith normal form decomposition to enable efficient sampling. Noise is modelled as probabilistic mixtures of Weyl operators that act only on the tableau's phase column. For fast simulation of noisy circuits, we support Pauli frames, respectively generalized Weyl frames, and introduce a noise-pushing technique that allows all noise processes to be consolidated into a single phase update at the end of the circuit. Using this representation, circuit fidelity can be determined entirely by the single accumulated phase-shift parameter $Δτ$, reducing fidelity estimation to a simple phase check per shot. Our codebase supports tableau simulation and conventional state-vector and density-matrix backends for qudits, and also includes circuit and tableau visualisations. Additionally, we provide tests and Jupyter notebooks for validation and illustration. This framework forms the basis for a scalable, open-source strong+weak stabilizer simulator including noise and can be found publicly available at https://github.com/QUICK-JKU/QuickQudits.

Comments
26 pages, 5 figures

Related papers

Monte Carlo Qudit Stabilizer Sampling with Pauli Frames

Adeeb Kabir, Yipeng Huang

We demonstrate that Pauli frames, a mathematical object used to efficiently sample qubit stabilizer circuits, can be used to efficiently sample qudit circuits of any prime dimension. Quantum error correcting codes are a central object of study in fault-tolerance quantum computing, which is necessary to realize useful universal quantum computing. Stabilizer codes belong to a large class of well-understood quantum codes that are crucially efficiently simulable on classical hardware using the so-called tableau representation. Popular simulation software like Stim exclusively studies and characterizes qubit codes, which entails repeatedly sampling a noisy quantum circuit with Pauli frames. We highlight the key differences (and fundamental similarities) between qubit and qudit stabilizer tableaus. It then teaches how Pauli frames generate batches of samples given a qudit tableau via two guided examples: one with a random measurement outcome but no noise, and one with a single noise gate.

IEEE International Conference on Quantum Computing and Engineering (QCE)•Aug 30, 2025•0 citations

Scalable Exact Non-Pauli Noise Simulation for Fault-Tolerant Quantum Circuits

Yuewen Hou, Gokul Subramanian Ravi

Fault-tolerant quantum design relies on simulating noisy quantum error correction (QEC) circuits with hundreds to thousands of physical qubits and large circuit depths. High-throughput stabilizer simulators, such as Stim, are the de facto tools for this regime. However, these simulators typically support only Clifford gates and stochastic Pauli noise. As a result, more realistic device noise models-such as coherent control errors and thermal relaxation-are often approximated using Pauli Twirling. This replacement offers speed but can introduce inaccuracies. To address this, we propose QuasiStim for scalable exact QEC simulation of realistic noise, with QuasiFold as a compiler frontend and DeltaFrame as a sampling backend. QuasiFold lowers sampling overhead by merging equivalent non-Pauli sites in the quantum circuit. DeltaFrame samples the resulting stabilizer trajectories by reusing precomputed predicted references instead of resimulating the full tableau for each shot. QuasiFold can reduce the number of shots by over $100\times$ on bivariate-bicycle (BB) codes. Compared to a per-shot tableau simulator, DeltaFrame is $166\times$--$1331\times$ faster on rotated surface-code memory under coherent $R_Z$ noise as $d$ grows from 5 to 35. QuasiStim also makes the PTA bias measurable: the gap is non-negligible and varies with noise strength, decoder, and circuit schedule. QuasiStim will be open-sourced upon publication at https://github.com/isaachyw/QuasiStim

quant-pharXiv:2610.15975•Oct 13, 2026•0 citations

Qudit noisy stabilizer formalism

Paul Aigner, Maria Flors Mor-Ruiz, Wolfgang Dür

We introduce the qudit noisy stabilizer formalism, a framework for efficiently describing the evolution of stabilizer states in prime-power dimensions subject to generalized Pauli-diagonal noise under Clifford operations and generalized Pauli measurements. For arbitrary dimensions, the formalism remains applicable, though restricted to a subset of stabilizer states and operations. The computational complexity scales linearly with the number of qudits in the initial state and exponentially with the number of qudits in the final state. This ensures that when noisy qudit stabilizer states evolve via generalized Pauli measurements and Clifford operations to generate multipartite entangled states of a few qudits, their description remains efficient. We demonstrate this by analyzing the generation of a generalized Bell pair from a noisy linear cluster state subject to two distinct noise sources acting on each of the qudits.

arXiv:2505.03889•Physical Review A•Jul 14, 2025•1 citation

Analytical and Compressed Simulation of Noisy Stabilizer Circuits

Paul Aigner, Jasmin Matti, Maria Flors Mor-Ruiz +2 more

We develop analytical and algorithmic techniques that enable efficient simulation of a broad class of noisy stabilizer circuits. We derive closed-form expressions of expectation values for tensor product of Paulis in circuits with non-deterministic Pauli measurements, yielding an efficient strong simulation method that avoids explicit density matrix construction and enables direct noise parameter sweeps. We introduce a circuit compression framework that reduces the per-sample cost of weak simulation in general noisy stabilizer circuits, including deterministic measurements, by separating parameter-independent preprocessing from sampling. Finally, we extend the analytical framework beyond its standard domain to include a small number of deterministic measurements, general rotations, and non-diagonal noise channels. Our results provide a unified framework for both strong and weak simulation of noisy stabilizer circuits and corresponds to an extension of the noisy stabilizer formalism introduced in \cite{PhysRevA.107.032424}. They offer applications ranging from calculation of the expectation values of entanglement witnesses, determination of reduced states, to energy evaluation.

quant-pharXiv:2604.22588•Apr 27, 2026•0 citations

Clifft: Fast Exact Simulation of Near-Clifford Quantum Circuits

Bradley A. Chase, Farrokh Labib

Exact classical simulation of fault-tolerant quantum circuits remains limited by a tradeoff between exponential state vector scaling, exponential $T$-count scaling in stabilizer-rank approaches, and per-shot tracking overhead in sparse generalized stabilizer simulators. In this work, we introduce Clifft, an open-source simulator that shifts the dominant exponential cost from the total qubit count to a dynamic active subspace by factoring the quantum state into an offline Clifford frame, an online Pauli frame, and a dynamically sized active state vector. This architecture resolves deterministic Clifford coordinate transformations ahead of time, generalizing Stim's compile-once, sample-many execution model to circuits with non-Clifford operations. Consequently, exponential simulation costs are determined by the peak active virtual dimension, which expands during non-Clifford operations and contracts during measurements. Clifft remains within a constant factor of standard tools in the pure-Clifford and non-Clifford limits, while delivering up to orders-of-magnitude throughput gains over GPU-accelerated near-Clifford simulators on low-magic fault-tolerant benchmarks. Executing on commodity CPUs and exposing a Stim-like API, Clifft enables, to our knowledge, the first exact end-to-end simulation of magic state cultivation including the escape stage, over hundreds of billions of shots. These simulations show that escape-stage failures suppress the discrepancy between the true $T$-gate circuit and its $S$-proxy at low decoder-gap thresholds, while at high thresholds the full-protocol behavior approaches the larger discrepancy observed in the cultivation stages alone.

quant-pharXiv:2604.27058•May 1, 2026•0 citations

Exact efficient simulation of noisy logical magic states using Clifford stabilizers

Yugo Takada, Stephen D. Bartlett, Dominic J. Williamson

The preparation of high-fidelity logical magic states is a crucial subroutine for universal fault-tolerant quantum computation (FTQC). Predicting the performance of FTQC and developing improved protocols rely on numerical methods to classically simulate logical magic state preparation in the presence of noise. Clifford logic on Pauli-stabilizer codes with circuit-level Pauli errors can be efficiently simulated using Pauli-stabilizer formalism, but the non-Clifford operations required to prepare logical magic states render generic simulation inefficient. We introduce Clifford-stabilizer simulation, an exact and efficient algorithm based on updating a Clifford-stabilizer group to simulate noisy preparation protocols for a broad class of logical magic states used to implement non-Clifford gates in the third level of the Clifford hierarchy under circuit-level Pauli errors. Clifford-stabilizer simulation applies to a range of operations that commonly appear in preparation protocols for such logical magic states, including Pauli-stabilizer measurements, logical Clifford measurements, and transversal non-Clifford gates. Our algorithm for Clifford-stabilizer simulation maps a non-Clifford circuit with sampled circuit-level Pauli errors to a Clifford circuit that exactly reproduces its measurement outcome distribution, achieving time and space complexities polynomial in relevant protocol parameters. We perform exact simulation of magic state cultivation up to fault distance 7 by Clifford-stabilizer simulation. Our method provides a route to perform exact benchmarking of large-scale logical magic state preparation protocols required for useful FTQC.

quant-pharXiv:2609.16929•Sep 16, 2026•0 citations
We present a comprehensive and self-contained framework for the efficient classical simulation of Clifford circuits acting on ddd-dimensional qudits, including realistic Pauli/Weyl noise via stochastic simulation. Our approach uses the stabilizer tableau formalism for qudits of arbitrary dimension and tracks both stabilizer and destabilizer generators under Clifford updates. The classical simulation remains efficient with simple algebraic Clifford update rules over Zd\mathbb{Z}_dZd​. Computational basis measurements in prime dimensions are handled by a generalized Aaronson-Gottesman (CHP) procedure. In composite dimensions, Zd\mathbb{Z}_dZd​ is not a field and the standard tableau reduction fails, so we employ an exact Smith normal form decomposition to enable efficient sampling. Noise is modelled as probabilistic mixtures of Weyl operators that act only on the tableau's phase column. For fast simulation of noisy circuits, we support Pauli frames, respectively generalized Weyl frames, and introduce a noise-pushing technique that allows all noise processes to be consolidated into a single phase update at the end of the circuit. Using this representation, circuit fidelity can be determined entirely by the single accumulated phase-shift parameter ΔτΔτΔτ, reducing fidelity estimation to a simple phase check per shot. Our codebase supports tableau simulation and conventional state-vector and density-matrix backends for qudits, and also includes circuit and tableau visualisations. Additionally, we provide tests and Jupyter notebooks for validation and illustration. This framework forms the basis for a scalable, open-source strong+weak stabilizer simulator including noise and can be found publicly available at https://github.com/QUICK-JKU/QuickQudits.
We demonstrate that Pauli frames, a mathematical object used to efficiently sample qubit stabilizer circuits, can be used to efficiently sample qudit circuits of any prime dimension. Quantum error correcting codes are a central object of study in fault-tolerance quantum computing, which is necessary to realize useful universal quantum computing. Stabilizer codes belong to a large class of well-understood quantum codes that are crucially efficiently simulable on classical hardware using the so-called tableau representation. Popular simulation software like Stim exclusively studies and characterizes qubit codes, which entails repeatedly sampling a noisy quantum circuit with Pauli frames. We highlight the key differences (and fundamental similarities) between qubit and qudit stabilizer tableaus. It then teaches how Pauli frames generate batches of samples given a qudit tableau via two guided examples: one with a random measurement outcome but no noise, and one with a single noise gate.
Fault-tolerant quantum design relies on simulating noisy quantum error correction (QEC) circuits with hundreds to thousands of physical qubits and large circuit depths. High-throughput stabilizer simulators, such as Stim, are the de facto tools for this regime. However, these simulators typically support only Clifford gates and stochastic Pauli noise. As a result, more realistic device noise models-such as coherent control errors and thermal relaxation-are often approximated using Pauli Twirling. This replacement offers speed but can introduce inaccuracies. To address this, we propose QuasiStim for scalable exact QEC simulation of realistic noise, with QuasiFold as a compiler frontend and DeltaFrame as a sampling backend. QuasiFold lowers sampling overhead by merging equivalent non-Pauli sites in the quantum circuit. DeltaFrame samples the resulting stabilizer trajectories by reusing precomputed predicted references instead of resimulating the full tableau for each shot. QuasiFold can reduce the number of shots by over 100×100\times100× on bivariate-bicycle (BB) codes. Compared to a per-shot tableau simulator, DeltaFrame is 166×166\times166×--1331×1331\times1331× faster on rotated surface-code memory under coherent RZR_ZRZ​ noise as ddd grows from 5 to 35. QuasiStim also makes the PTA bias measurable: the gap is non-negligible and varies with noise strength, decoder, and circuit schedule. QuasiStim will be open-sourced upon publication at https://github.com/isaachyw/QuasiStim
We introduce the qudit noisy stabilizer formalism, a framework for efficiently describing the evolution of stabilizer states in prime-power dimensions subject to generalized Pauli-diagonal noise under Clifford operations and generalized Pauli measurements. For arbitrary dimensions, the formalism remains applicable, though restricted to a subset of stabilizer states and operations. The computational complexity scales linearly with the number of qudits in the initial state and exponentially with the number of qudits in the final state. This ensures that when noisy qudit stabilizer states evolve via generalized Pauli measurements and Clifford operations to generate multipartite entangled states of a few qudits, their description remains efficient. We demonstrate this by analyzing the generation of a generalized Bell pair from a noisy linear cluster state subject to two distinct noise sources acting on each of the qudits.
We develop analytical and algorithmic techniques that enable efficient simulation of a broad class of noisy stabilizer circuits. We derive closed-form expressions of expectation values for tensor product of Paulis in circuits with non-deterministic Pauli measurements, yielding an efficient strong simulation method that avoids explicit density matrix construction and enables direct noise parameter sweeps. We introduce a circuit compression framework that reduces the per-sample cost of weak simulation in general noisy stabilizer circuits, including deterministic measurements, by separating parameter-independent preprocessing from sampling. Finally, we extend the analytical framework beyond its standard domain to include a small number of deterministic measurements, general rotations, and non-diagonal noise channels. Our results provide a unified framework for both strong and weak simulation of noisy stabilizer circuits and corresponds to an extension of the noisy stabilizer formalism introduced in \cite{PhysRevA.107.032424}. They offer applications ranging from calculation of the expectation values of entanglement witnesses, determination of reduced states, to energy evaluation.
Exact classical simulation of fault-tolerant quantum circuits remains limited by a tradeoff between exponential state vector scaling, exponential TTT-count scaling in stabilizer-rank approaches, and per-shot tracking overhead in sparse generalized stabilizer simulators. In this work, we introduce Clifft, an open-source simulator that shifts the dominant exponential cost from the total qubit count to a dynamic active subspace by factoring the quantum state into an offline Clifford frame, an online Pauli frame, and a dynamically sized active state vector. This architecture resolves deterministic Clifford coordinate transformations ahead of time, generalizing Stim's compile-once, sample-many execution model to circuits with non-Clifford operations. Consequently, exponential simulation costs are determined by the peak active virtual dimension, which expands during non-Clifford operations and contracts during measurements. Clifft remains within a constant factor of standard tools in the pure-Clifford and non-Clifford limits, while delivering up to orders-of-magnitude throughput gains over GPU-accelerated near-Clifford simulators on low-magic fault-tolerant benchmarks. Executing on commodity CPUs and exposing a Stim-like API, Clifft enables, to our knowledge, the first exact end-to-end simulation of magic state cultivation including the escape stage, over hundreds of billions of shots. These simulations show that escape-stage failures suppress the discrepancy between the true TTT-gate circuit and its SSS-proxy at low decoder-gap thresholds, while at high thresholds the full-protocol behavior approaches the larger discrepancy observed in the cultivation stages alone.
The preparation of high-fidelity logical magic states is a crucial subroutine for universal fault-tolerant quantum computation (FTQC). Predicting the performance of FTQC and developing improved protocols rely on numerical methods to classically simulate logical magic state preparation in the presence of noise. Clifford logic on Pauli-stabilizer codes with circuit-level Pauli errors can be efficiently simulated using Pauli-stabilizer formalism, but the non-Clifford operations required to prepare logical magic states render generic simulation inefficient. We introduce Clifford-stabilizer simulation, an exact and efficient algorithm based on updating a Clifford-stabilizer group to simulate noisy preparation protocols for a broad class of logical magic states used to implement non-Clifford gates in the third level of the Clifford hierarchy under circuit-level Pauli errors. Clifford-stabilizer simulation applies to a range of operations that commonly appear in preparation protocols for such logical magic states, including Pauli-stabilizer measurements, logical Clifford measurements, and transversal non-Clifford gates. Our algorithm for Clifford-stabilizer simulation maps a non-Clifford circuit with sampled circuit-level Pauli errors to a Clifford circuit that exactly reproduces its measurement outcome distribution, achieving time and space complexities polynomial in relevant protocol parameters. We perform exact simulation of magic state cultivation up to fault distance 7 by Clifford-stabilizer simulation. Our method provides a route to perform exact benchmarking of large-scale logical magic state preparation protocols required for useful FTQC.