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

Simulation of IBM's kicked Ising experiment with Projected Entangled Pair Operator

Haijun Liao, Kang Wang, Zong-Sheng Zhou, Pan Zhang, Tao Xiang

Published Aug 6, 2023• 4 citations
arXivHTMLPDF

Abstract

We perform classical simulations of the 127-qubit kicked Ising model, which was recently emulated using a quantum circuit with error mitigation [Nature 618, 500 (2023)]. Our approach is based on the projected entangled pair operator (PEPO) in the Heisenberg picture. Its main feature is the ability to automatically identify the underlying low-rank and low-entanglement structures in the quantum circuit involving Clifford and near-Clifford gates. We assess our approach using the quantum circuit with 5+1 trotter steps which was previously considered beyond classical verification. We develop a Clifford expansion theory to compute exact expectation values and use them to evaluate algorithms. The results indicate that PEPO significantly outperforms existing methods, including the tensor network with belief propagation, the matrix product operator, and the Clifford perturbation theory, in both efficiency and accuracy. In particular, PEPO with bond dimension $χ=2$ already gives similar accuracy to the CPT with $K=10$ and MPO with bond dimension $χ=1024$. And PEPO with $χ=184$ provides exact results in $3$ seconds using a single CPU. Furthermore, we apply our method to the circuit with 20 Trotter steps. We observe the monotonic and consistent convergence of the results with $χ$, allowing us to estimate the outcome with $χ\to\infty$ through extrapolations. We then compare the extrapolated results to those achieved in quantum hardware and with existing tensor network methods. Additionally, we discuss the potential usefulness of our approach in simulating quantum circuits, especially in scenarios involving near-Clifford circuits and quantum approximate optimization algorithms. Our approach is the first use of PEPO in solving the time evolution problem, and our results suggest it could be a powerful tool for exploring the dynamical properties of quantum many-body systems.

DOI
10.48550/arxiv.2308.03082

Related papers

Complexity of tensor network simulation for noisy quantum circuits

Yuguo Shao, Zishuo Zhao, Song Cheng +1 more

We aim to rigorously address how local noise affects classical simulability of quantum dynamics benchmarked by tensor-network methods. Using operator entanglement entropy (OEE), we prove the following: (1) For single-qubit depolarizing noise on arbitrary circuits, tensor networks with $\mathrm{poly}(n)$ bond dimension suffice for fixed absolute Hilbert-Schmidt error after $\order{1}$ depth, while relative error demands $\order{\log n}$ depth; and this bound is optimal. (2) For single-qubit depolarizing noise on 1D local circuits, the existence of whole-trajectory error-bounded matrix product operator (MPO) of $\mathrm{poly}(n)$ bond dimension at all depths. (3) For general single-qubit noise in 1D brickwall circuits, random two-design gates with contraction coefficient $c<1/3$ yield an $\order{1}$ OEE plateau with probability $1-Te^{-Ω(n)}$, while arbitrary gates with $c<1/48$ give $\order{\log n}$ OEE in the worst case. (4) In higher dimensions, these bounds yield uniform-in-depth $\mathrm{poly}(n)$ average boundary-bond dimensions for projected entangled pair operators~(PEPO) across every cut -- under depolarizing noise at either absolute or relative accuracy, and under general noise with strong contraction at absolute accuracy. Our results establish a rigorous connection between certain noise models, circuit types, and their classical simulability.

cond-mat.str-elCross-listedarXiv:2606.00474•Jun 2, 2026•0 citations

Fast and converged classical simulations of evidence for the utility of quantum computing before fault tolerance

Tomislav Begušić, Johnnie Gray, Garnet Kin‐Lic Chan

A recent quantum simulation of observables of the kicked Ising model on 127 qubits implemented circuits that exceed the capabilities of exact classical simulation. We show that several approximate classical methods, based on sparse Pauli dynamics and tensor network algorithms, can simulate these observables orders of magnitude faster than the quantum experiment and can also be systematically converged beyond the experimental accuracy. Our most accurate technique combines a mixed Schrödinger and Heisenberg tensor network representation with the Bethe free entropy relation of belief propagation to compute expectation values with an effective wave function-operator sandwich bond dimension >16,000,000, achieving an absolute accuracy, without extrapolation, in the observables of <0.01, which is converged for many practical purposes. We thereby identify inaccuracies in the experimental extrapolations and suggest how future experiments can be implemented to increase the classical hardness.

arXiv:2308.05077•Science Advances•Jan 17, 2024•89 citations

Efficient tensor network simulation of IBM's largest quantum processors

Siddhartha Patra, Saeed S. Jahromi, Sukhbinder Singh +1 more

We show how quantum-inspired 2D tensor networks can be used to efficiently and accurately simulate the largest quantum processors from IBM, namely Eagle (127 qubits), Osprey (433 qubits), and Condor (1121 qubits). We simulate the dynamics of a complex quantum many-body system—specifically, the kicked Ising experiment considered recently by IBM in Y. Kim , —using graph-based projected entangled pair states (gPEPS), which was proposed by some of us in . Our results show that simple tensor updates are already sufficient to achieve very large unprecedented accuracy with remarkably low computational resources for this model. Apart from simulating the original experiment for 127 qubits, we also extend our results to 433 and 1121 qubits, and for evolution times around eight times longer, thus setting a benchmark for the newest IBM quantum machines. We also report accurate simulations for infinitely many qubits. Our results show that gPEPS are a natural tool to efficiently simulate quantum computers with an underlying lattice-based qubit connectivity, such as all quantum processors based on superconducting qubits. Published by the American Physical Society 2024

arXiv:2309.15642•Physical Review Research•Mar 26, 2024•63 citations

Clifford Perturbation Approximation for Quantum Error Mitigation

Zhang, Ruiqi, Yuguo Shao, Fuchuan Wei +3 more

Quantum error mitigation (QEM) is critical for harnessing the potential of near-term quantum devices. Particularly, QEM protocols can be designed based on machine learning, where the mapping between noisy computational outputs and ideal ones can be learned on a training set consisting of Clifford circuits or near-Clifford circuits that contain only a limited number of non-Clifford gates. This learned mapping is then applied to noisy target circuits to estimate the ideal computational output. In this work, we propose a learning-based error mitigation framework called Clifford Perturbation Data Regression (CPDR), which constructs training sets by Clifford circuits with small perturbations. Specifically, these circuits are parameterized quantum circuits, where the rotation angles of the gates are restricted to a narrow range, ensuring that the gates remain close to Clifford gates. This design enables the efficient simulation of the training circuits using the Sparse Pauli Dynamics method. As a result, CPDR is able to utilize training sets with a better diversity to train the model, compared with previous learning-based QEM models that construct training sets with only Clifford or near-Clifford circuits. Numerical simulations on small-scale Ising model circuits demonstrate that the performance of CPDR dramatically outperforms that of existing methods such as Zero-Noise Extrapolation and learning-based Probabilistic Error Cancellation. Furthermore, using the experimental data from IBM's 127-qubit Eagle processor, our findings suggest that CPDR demonstrates improved accuracy compared to the original mitigation results reported in [Nature 618, 500 (2023)].

arXiv:2412.09518•Dec 12, 2024•0 citations

Simulating quantum circuit expectation values by Clifford perturbation theory

Tomislav Begušić, Kasra Hejazi, Garnet Kin‐Lic Chan

The classical simulation of quantum circuits is of central importance for benchmarking near-term quantum devices. The fact that gates belonging to the Clifford group can be simulated efficiently on classical computers has motivated a range of methods that scale exponentially only in the number of non-Clifford gates. Here, we consider the expectation value problem for circuits composed of Clifford gates and non-Clifford Pauli rotations and introduce a heuristic perturbative approach based on the truncation of the exponentially growing sum of Pauli terms in the Heisenberg picture. Numerical results are shown on a quantum approximate optimization algorithm benchmark for the E3LIN2 problem, and we also demonstrate how this method can be used to quantify coherent and incoherent errors of local observables in Clifford circuits. Our results indicate that this systematically improvable perturbative method offers a viable alternative to exact methods for approximating expectation values of large near-Clifford circuits.

arXiv:2306.04797•The Journal of Chemical Physics•Apr 15, 2025•16 citations

Combining tensor networks and Pauli propagation to simulate quantum dynamics

Alessandro Cerati

The study of classical simulability is a key part of quantum information theory, providing new tools for fundamental research and the development of quantum computing. An increasingly successful framework for the study of simulability is resource theory. The best-known and best-understood quantum resource is entanglement, but in recent years nonstabilizerness, or magic, has also been gathering attention from the community. It is not yet well understood to what extent different resources play off each other in making systems hard to classically simulate. It has recently been shown that Clifford-augmented matrix product states (CAMPS) allow for efficient circuit simulation in some cases where both entanglement and magic are present. Another emerging method which has shown the ability to simulate systems with both resources is Pauli propagation (PP). Furthermore, the interplay between entanglement and magic at the level of physical phenomena is also poorly understood, in particular in the context of real-time dynamics. In this work, I implement an analytical algorithm proposed by Fux et al. (2025) which allows CAMPS to efficiently simulate essentially the first N non-Clifford Pauli exponential gates in an N-qubit quantum circuit, reproducing results obtained by greedy-search entanglement cooling in the paper. I then develop CAMPS-PP, a method for expectation value dynamics which evolves the state as a CAMPS and the observable in the Heisenberg picture using PP. I show that this method is significantly faster than pure CAMPS for slightly deep random Pauli rotation circuits at the tested system sizes, up to N=100, with improvements of up to two orders of magnitude for N=20. Finally, in order to investigate the physical relevance of magic I apply CAMPS-PP to an XXZ Heisenberg Floquet circuit with Clifford and non-Clifford angles, finding that CAMPS-PP slightly extends the simulable non-Clifford gate density and significantly reduces memory requirements compared to PP.

•Sep 24, 2026•0 citations
We perform classical simulations of the 127-qubit kicked Ising model, which was recently emulated using a quantum circuit with error mitigation [Nature 618, 500 (2023)]. Our approach is based on the projected entangled pair operator (PEPO) in the Heisenberg picture. Its main feature is the ability to automatically identify the underlying low-rank and low-entanglement structures in the quantum circuit involving Clifford and near-Clifford gates. We assess our approach using the quantum circuit with 5+1 trotter steps which was previously considered beyond classical verification. We develop a Clifford expansion theory to compute exact expectation values and use them to evaluate algorithms. The results indicate that PEPO significantly outperforms existing methods, including the tensor network with belief propagation, the matrix product operator, and the Clifford perturbation theory, in both efficiency and accuracy. In particular, PEPO with bond dimension χ=2χ=2χ=2 already gives similar accuracy to the CPT with K=10K=10K=10 and MPO with bond dimension χ=1024χ=1024χ=1024. And PEPO with χ=184χ=184χ=184 provides exact results in 333 seconds using a single CPU. Furthermore, we apply our method to the circuit with 20 Trotter steps. We observe the monotonic and consistent convergence of the results with χχχ, allowing us to estimate the outcome with χ→∞χ\to\inftyχ→∞ through extrapolations. We then compare the extrapolated results to those achieved in quantum hardware and with existing tensor network methods. Additionally, we discuss the potential usefulness of our approach in simulating quantum circuits, especially in scenarios involving near-Clifford circuits and quantum approximate optimization algorithms. Our approach is the first use of PEPO in solving the time evolution problem, and our results suggest it could be a powerful tool for exploring the dynamical properties of quantum many-body systems.
We aim to rigorously address how local noise affects classical simulability of quantum dynamics benchmarked by tensor-network methods. Using operator entanglement entropy (OEE), we prove the following: (1) For single-qubit depolarizing noise on arbitrary circuits, tensor networks with poly(n)\mathrm{poly}(n)poly(n) bond dimension suffice for fixed absolute Hilbert-Schmidt error after $\order{1}$ depth, while relative error demands $\order{\log n}$ depth; and this bound is optimal. (2) For single-qubit depolarizing noise on 1D local circuits, the existence of whole-trajectory error-bounded matrix product operator (MPO) of poly(n)\mathrm{poly}(n)poly(n) bond dimension at all depths. (3) For general single-qubit noise in 1D brickwall circuits, random two-design gates with contraction coefficient c<1/3c<1/3c<1/3 yield an $\order{1}$ OEE plateau with probability 1−Te−Ω(n)1-Te^{-Ω(n)}1−Te−Ω(n), while arbitrary gates with c<1/48c<1/48c<1/48 give $\order{\log n}$ OEE in the worst case. (4) In higher dimensions, these bounds yield uniform-in-depth poly(n)\mathrm{poly}(n)poly(n) average boundary-bond dimensions for projected entangled pair operators~(PEPO) across every cut -- under depolarizing noise at either absolute or relative accuracy, and under general noise with strong contraction at absolute accuracy. Our results establish a rigorous connection between certain noise models, circuit types, and their classical simulability.
A recent quantum simulation of observables of the kicked Ising model on 127 qubits implemented circuits that exceed the capabilities of exact classical simulation. We show that several approximate classical methods, based on sparse Pauli dynamics and tensor network algorithms, can simulate these observables orders of magnitude faster than the quantum experiment and can also be systematically converged beyond the experimental accuracy. Our most accurate technique combines a mixed Schrödinger and Heisenberg tensor network representation with the Bethe free entropy relation of belief propagation to compute expectation values with an effective wave function-operator sandwich bond dimension >16,000,000, achieving an absolute accuracy, without extrapolation, in the observables of <0.01, which is converged for many practical purposes. We thereby identify inaccuracies in the experimental extrapolations and suggest how future experiments can be implemented to increase the classical hardness.
We show how quantum-inspired 2D tensor networks can be used to efficiently and accurately simulate the largest quantum processors from IBM, namely Eagle (127 qubits), Osprey (433 qubits), and Condor (1121 qubits). We simulate the dynamics of a complex quantum many-body system—specifically, the kicked Ising experiment considered recently by IBM in Y. Kim , —using graph-based projected entangled pair states (gPEPS), which was proposed by some of us in . Our results show that simple tensor updates are already sufficient to achieve very large unprecedented accuracy with remarkably low computational resources for this model. Apart from simulating the original experiment for 127 qubits, we also extend our results to 433 and 1121 qubits, and for evolution times around eight times longer, thus setting a benchmark for the newest IBM quantum machines. We also report accurate simulations for infinitely many qubits. Our results show that gPEPS are a natural tool to efficiently simulate quantum computers with an underlying lattice-based qubit connectivity, such as all quantum processors based on superconducting qubits. Published by the American Physical Society 2024
Quantum error mitigation (QEM) is critical for harnessing the potential of near-term quantum devices. Particularly, QEM protocols can be designed based on machine learning, where the mapping between noisy computational outputs and ideal ones can be learned on a training set consisting of Clifford circuits or near-Clifford circuits that contain only a limited number of non-Clifford gates. This learned mapping is then applied to noisy target circuits to estimate the ideal computational output. In this work, we propose a learning-based error mitigation framework called Clifford Perturbation Data Regression (CPDR), which constructs training sets by Clifford circuits with small perturbations. Specifically, these circuits are parameterized quantum circuits, where the rotation angles of the gates are restricted to a narrow range, ensuring that the gates remain close to Clifford gates. This design enables the efficient simulation of the training circuits using the Sparse Pauli Dynamics method. As a result, CPDR is able to utilize training sets with a better diversity to train the model, compared with previous learning-based QEM models that construct training sets with only Clifford or near-Clifford circuits. Numerical simulations on small-scale Ising model circuits demonstrate that the performance of CPDR dramatically outperforms that of existing methods such as Zero-Noise Extrapolation and learning-based Probabilistic Error Cancellation. Furthermore, using the experimental data from IBM's 127-qubit Eagle processor, our findings suggest that CPDR demonstrates improved accuracy compared to the original mitigation results reported in [Nature 618, 500 (2023)].
The classical simulation of quantum circuits is of central importance for benchmarking near-term quantum devices. The fact that gates belonging to the Clifford group can be simulated efficiently on classical computers has motivated a range of methods that scale exponentially only in the number of non-Clifford gates. Here, we consider the expectation value problem for circuits composed of Clifford gates and non-Clifford Pauli rotations and introduce a heuristic perturbative approach based on the truncation of the exponentially growing sum of Pauli terms in the Heisenberg picture. Numerical results are shown on a quantum approximate optimization algorithm benchmark for the E3LIN2 problem, and we also demonstrate how this method can be used to quantify coherent and incoherent errors of local observables in Clifford circuits. Our results indicate that this systematically improvable perturbative method offers a viable alternative to exact methods for approximating expectation values of large near-Clifford circuits.
The study of classical simulability is a key part of quantum information theory, providing new tools for fundamental research and the development of quantum computing. An increasingly successful framework for the study of simulability is resource theory. The best-known and best-understood quantum resource is entanglement, but in recent years nonstabilizerness, or magic, has also been gathering attention from the community. It is not yet well understood to what extent different resources play off each other in making systems hard to classically simulate. It has recently been shown that Clifford-augmented matrix product states (CAMPS) allow for efficient circuit simulation in some cases where both entanglement and magic are present. Another emerging method which has shown the ability to simulate systems with both resources is Pauli propagation (PP). Furthermore, the interplay between entanglement and magic at the level of physical phenomena is also poorly understood, in particular in the context of real-time dynamics. In this work, I implement an analytical algorithm proposed by Fux et al. (2025) which allows CAMPS to efficiently simulate essentially the first N non-Clifford Pauli exponential gates in an N-qubit quantum circuit, reproducing results obtained by greedy-search entanglement cooling in the paper. I then develop CAMPS-PP, a method for expectation value dynamics which evolves the state as a CAMPS and the observable in the Heisenberg picture using PP. I show that this method is significantly faster than pure CAMPS for slightly deep random Pauli rotation circuits at the tested system sizes, up to N=100, with improvements of up to two orders of magnitude for N=20. Finally, in order to investigate the physical relevance of magic I apply CAMPS-PP to an XXZ Heisenberg Floquet circuit with Clifford and non-Clifford angles, finding that CAMPS-PP slightly extends the simulable non-Clifford gate density and significantly reduces memory requirements compared to PP.