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:1604.00203 · Physical Review A

Digital quantum simulation of many-body non-Markovian dynamics

Ryan Sweke, Mikel Sanz, Ilya Sinayskiy, Francesco Petruccione, E. Solano

Published Aug 15, 2016• 55 citations
arXivHTMLPDF

Abstract

We present an algorithmic method for the digital quantum simulation of many-body locally indivisible non-Markovian open quantum systems. It consists of two parts: first, a Suzuki-Lie-Trotter decomposition of the global system propagator into the product of subsystem propagators, which may not be quantum channels, and second, an algorithmic procedure for the implementation of the subsystem propagators through unitary operations and measurements on a dilated space. By providing rigorous error bounds for the relevant Suzuki-Lie-Trotter decomposition, we are able to analyze the efficiency of the method, and connect it with an appropriate measure of the local indivisibility of the system. In light of our analysis, the proposed method is expected to be experimentally achievable for a variety of interesting cases.

DOI
10.1103/physreva.94.022317

Related papers

Faster Quantum Simulation Of Markovian Open Quantum Systems Via Randomisation

I. J. David, Ilya Sinayskiy, Francesco Petruccione

When simulating the dynamics of open quantum systems with quantum computers, it is essential to accurately approximate the system's behaviour while preserving the physicality of its evolution. Traditionally, for Markovian open quantum systems, this has been achieved using first and second-order Trotter-Suzuki product formulas or probabilistic algorithms. In this work, we introduce novel non-probabilistic algorithms for simulating Markovian open quantum systems using randomisation. Our methods, including first and second-order randomised Trotter-Suzuki formulas and the QDRIFT channel, not only maintain the physicality of the system's evolution but also enhance the scalability and precision of quantum simulations. We derive error bounds and step count limits for these techniques, bypassing the need for the mixing lemma typically employed in Hamiltonian simulation proofs. Furthermore, we implement these randomised algorithms using Classical Sampling (CS), demonstrating their gate complexity advantages over deterministic TS product formulas. This work systematically extends powerful randomisation techniques from Hamiltonian simulation to the general setting of Markovian open quantum systems, highlighting their potential to enable faster and more accurate simulations.

arXiv:2408.11683•Quantum•Sep 3, 2026•1 citation

Qudit-native simulation of the Potts model

Maksim A. Gavreev, Evgeniy O. Kiktenko, Aleksey K. Fedorov +1 more

Simulating entangled, many-body quantum systems is notoriously hard, especially in the case of high-dimensional nature of physical underlying objects. In this work, we propose an approach for simulating the Potts model based on the Suzuki-Trotter decomposition that we construct for qudit systems. Specifically, we introduce two qudit-native decomposition schemes: (i) the first utilizes Molmer-Sorensen gate and additional local levels to encode the Potts interactions, while (ii) the second employs an light-shift gate that naturally fits qudit architectures. These decompositions enable a direct and efficient mapping of the Potts model dynamics into hardware-efficient qudit gate sequences for trapped-ion platform. Furthermore, we demonstrate the use of a Suzuki-Trotter approximation with our evolution-into-gates framework, for detecting the dynamical quantum phase transition. Our results establish a pathway toward qudit-based digital quantum simulation of many-body models and provide a new perspective on probing nonanalytic behavior in high-dimensional quantum many-body models.

quant-pharXiv:2511.13572•Entropy•Nov 18, 2025•0 citations

Efficient Trotter-Suzuki Schemes for Long-time Quantum Dynamics

Marko Maležič, Johann Ostmeyer

Accurately simulating long-time dynamics of many-body systems is a challenge in both classical and quantum computing due to the accumulation of Trotter errors. While low-order Trotter-Suzuki decompositions are straightforward to implement, their rapidly growing error limits access to long-time observables. We present a framework for constructing efficient high-order Trotter-Suzuki schemes by identifying their structure and directly optimizing their parameters over a high-dimensional space. This method enables the discovery of new schemes with significantly improved efficiency compared to traditional constructions, such as those by Suzuki and Yoshida. Based on the theoretical efficiency and practical performance, we recommend two novel highly efficient schemes at $4^{\textrm{th}}$ and $6^{\textrm{th}}$ order. We also demonstrate the effectiveness of these decompositions on the Heisenberg model and the quantum harmonic oscillator, and find that for a fixed final time they perform better across the computational cost. Even when using large time steps, they surpass established low-order schemes like the Leapfrog. Finally, we investigate the in-practice performance of different Trotter schemes and find the decompositions with more uniform coefficients tend to feature improved error accumulation over long times. We have included this observation into our choice of recommended schemes.

quant-pharXiv:2601.18756•Journal of Physics A Mathematical and Theoretical•Jan 27, 2026•1 citation

Quantum simulation of many-body dynamics with noise-robust Trotter decomposition based on symmetric structures

Bo Yang, Negishi, Naoki

The Suzuki-Trotter decomposition, which digitalizes quantum time evolution, provides a promising framework for simulating quantum dynamics on quantum hardware and exploring quantum advantage over classical computation. However, conventional Trotter circuits require a large number of non-local gates, lowering their faithfulness to the ideal dynamics when implemented on current noisy quantum hardware. While most previous studies have focused on circuit optimization, we instead propose a new Trotter decomposition that is intrinsically circuit-efficient for simulating quantum dynamics on near-term devices. Our method substantially reduces both the residual error by Trotter decomposition and the number of CNOT operations compared to conventional Trotter decompositions by exploiting the symmetry of the target model to construct an effective Hamiltonian with fewer two-qubit gates. We demonstrate the noise robustness of the proposed approach through numerical simulations of a nine-site Heisenberg model under realistic noise, and further validate its experimental practicality on the IBM superconducting device, achieving a state fidelity exceeding $0.98$ when combined with quantum error mitigation in the three-site case. The proposed circuit design is also compatible with existing circuit optimization techniques. Our results establish a practical route toward noise-resilient quantum simulation in many-body dynamics.

arXiv:2505.04552•May 7, 2025•0 citations

Suzuki-Trotter Decompositions and other Methods for Quantum Time Evolution

Johann Ostmeyer

(Suzuki-)Trotter decompositions, splitting methods, (Lie) product formulae... The most common numerical methods for the time evolution of quantum systems come with many names. And they are used practically everywhere with applications ranging from the solution of classical equations of motion and various Monte Carlo simulations to the real and imaginary time evolution on classical as well as quantum computers. Here we review the state of the art of said methods, focussing especially on the progress made over the last few years. We highlight recently discovered efficient time evolution algorithms and explain how best to use them in practice. A central part of this work is the estimation of error bounds that has improved greatly within the past decade. The relevance of time evolution methods for quantum computing is discussed with a focus on noisy hardware. Finally, a comprehensive overview of generalisations, related methods and alternatives to Trotterization is provided. This includes time-dependent Hamiltonian dynamics, processed methods, multi-product formulae, symplectic integrators, TDVP for tensor networks, quantum signal processing, Crouch-Grossman methods and more. The overall perspective in this work is that of a theoretical physicist. All mathematical proofs as well as some technical details are omitted for easier readability. Instead, this review serves as a hands-on guide and, of course, as a starting point for references that provide further details.

quant-pharXiv:2609.19184•Sep 16, 2026•0 citations

Optimizing digital quantum simulation of open quantum lattice models

Xie-Hang Yu, Hongchao Li, J. I. Cirac +1 more

Many-body systems arising in condensed matter physics and quantum optics inevitably couple to the environment and need to be modelled as open quantum systems. While near-optimal algorithms have been developed for simulating many-body quantum dynamics, algorithms for their open system counterparts remain less well investigated. We address the problem of simulating geometrically local many-body open quantum systems interacting with a stationary Gaussian environment. Under a smoothness assumption on the system-environment interaction, we develop near-optimal algorithms that, for a model with $N$ spins and evolution time $t$, attain a simulation error $δ$ in the system-state with $\mathcal{O}(Nt(Nt/δ)^{o(1)})$ gates, $\mathcal{O}(t(Nt/δ)^{o(1)})$ parallelized circuit depth and $\tilde{\mathcal{O}}(N(Nt/δ)^{o(1)})$ ancillas. We additionally show that, if only simulating local observables is of interest, then the circuit depth of the digital algorithm can be chosen to be independent of the system size $N$. This provides theoretical evidence for the utility of these algorithms for simulating physically relevant models, where typically local observables are of interest, on pre-fault tolerant devices. Finally, for the limiting case of Markovian dynamics with commuting jump operators, we propose two algorithms based on sampling a Wiener process and on a locally dilated Hamiltonian construction, respectively. These algorithms reduce the asymptotic gate complexity on $N$ compared to currently available algorithms in terms of the required number of geometrically local gates.

arXiv:2509.02268•Sep 2, 2025•0 citations
We present an algorithmic method for the digital quantum simulation of many-body locally indivisible non-Markovian open quantum systems. It consists of two parts: first, a Suzuki-Lie-Trotter decomposition of the global system propagator into the product of subsystem propagators, which may not be quantum channels, and second, an algorithmic procedure for the implementation of the subsystem propagators through unitary operations and measurements on a dilated space. By providing rigorous error bounds for the relevant Suzuki-Lie-Trotter decomposition, we are able to analyze the efficiency of the method, and connect it with an appropriate measure of the local indivisibility of the system. In light of our analysis, the proposed method is expected to be experimentally achievable for a variety of interesting cases.
When simulating the dynamics of open quantum systems with quantum computers, it is essential to accurately approximate the system's behaviour while preserving the physicality of its evolution. Traditionally, for Markovian open quantum systems, this has been achieved using first and second-order Trotter-Suzuki product formulas or probabilistic algorithms. In this work, we introduce novel non-probabilistic algorithms for simulating Markovian open quantum systems using randomisation. Our methods, including first and second-order randomised Trotter-Suzuki formulas and the QDRIFT channel, not only maintain the physicality of the system's evolution but also enhance the scalability and precision of quantum simulations. We derive error bounds and step count limits for these techniques, bypassing the need for the mixing lemma typically employed in Hamiltonian simulation proofs. Furthermore, we implement these randomised algorithms using Classical Sampling (CS), demonstrating their gate complexity advantages over deterministic TS product formulas. This work systematically extends powerful randomisation techniques from Hamiltonian simulation to the general setting of Markovian open quantum systems, highlighting their potential to enable faster and more accurate simulations.
Simulating entangled, many-body quantum systems is notoriously hard, especially in the case of high-dimensional nature of physical underlying objects. In this work, we propose an approach for simulating the Potts model based on the Suzuki-Trotter decomposition that we construct for qudit systems. Specifically, we introduce two qudit-native decomposition schemes: (i) the first utilizes Molmer-Sorensen gate and additional local levels to encode the Potts interactions, while (ii) the second employs an light-shift gate that naturally fits qudit architectures. These decompositions enable a direct and efficient mapping of the Potts model dynamics into hardware-efficient qudit gate sequences for trapped-ion platform. Furthermore, we demonstrate the use of a Suzuki-Trotter approximation with our evolution-into-gates framework, for detecting the dynamical quantum phase transition. Our results establish a pathway toward qudit-based digital quantum simulation of many-body models and provide a new perspective on probing nonanalytic behavior in high-dimensional quantum many-body models.
Accurately simulating long-time dynamics of many-body systems is a challenge in both classical and quantum computing due to the accumulation of Trotter errors. While low-order Trotter-Suzuki decompositions are straightforward to implement, their rapidly growing error limits access to long-time observables. We present a framework for constructing efficient high-order Trotter-Suzuki schemes by identifying their structure and directly optimizing their parameters over a high-dimensional space. This method enables the discovery of new schemes with significantly improved efficiency compared to traditional constructions, such as those by Suzuki and Yoshida. Based on the theoretical efficiency and practical performance, we recommend two novel highly efficient schemes at 4th4^{\textrm{th}}4th and 6th6^{\textrm{th}}6th order. We also demonstrate the effectiveness of these decompositions on the Heisenberg model and the quantum harmonic oscillator, and find that for a fixed final time they perform better across the computational cost. Even when using large time steps, they surpass established low-order schemes like the Leapfrog. Finally, we investigate the in-practice performance of different Trotter schemes and find the decompositions with more uniform coefficients tend to feature improved error accumulation over long times. We have included this observation into our choice of recommended schemes.
The Suzuki-Trotter decomposition, which digitalizes quantum time evolution, provides a promising framework for simulating quantum dynamics on quantum hardware and exploring quantum advantage over classical computation. However, conventional Trotter circuits require a large number of non-local gates, lowering their faithfulness to the ideal dynamics when implemented on current noisy quantum hardware. While most previous studies have focused on circuit optimization, we instead propose a new Trotter decomposition that is intrinsically circuit-efficient for simulating quantum dynamics on near-term devices. Our method substantially reduces both the residual error by Trotter decomposition and the number of CNOT operations compared to conventional Trotter decompositions by exploiting the symmetry of the target model to construct an effective Hamiltonian with fewer two-qubit gates. We demonstrate the noise robustness of the proposed approach through numerical simulations of a nine-site Heisenberg model under realistic noise, and further validate its experimental practicality on the IBM superconducting device, achieving a state fidelity exceeding 0.980.980.98 when combined with quantum error mitigation in the three-site case. The proposed circuit design is also compatible with existing circuit optimization techniques. Our results establish a practical route toward noise-resilient quantum simulation in many-body dynamics.
(Suzuki-)Trotter decompositions, splitting methods, (Lie) product formulae... The most common numerical methods for the time evolution of quantum systems come with many names. And they are used practically everywhere with applications ranging from the solution of classical equations of motion and various Monte Carlo simulations to the real and imaginary time evolution on classical as well as quantum computers. Here we review the state of the art of said methods, focussing especially on the progress made over the last few years. We highlight recently discovered efficient time evolution algorithms and explain how best to use them in practice. A central part of this work is the estimation of error bounds that has improved greatly within the past decade. The relevance of time evolution methods for quantum computing is discussed with a focus on noisy hardware. Finally, a comprehensive overview of generalisations, related methods and alternatives to Trotterization is provided. This includes time-dependent Hamiltonian dynamics, processed methods, multi-product formulae, symplectic integrators, TDVP for tensor networks, quantum signal processing, Crouch-Grossman methods and more. The overall perspective in this work is that of a theoretical physicist. All mathematical proofs as well as some technical details are omitted for easier readability. Instead, this review serves as a hands-on guide and, of course, as a starting point for references that provide further details.
Many-body systems arising in condensed matter physics and quantum optics inevitably couple to the environment and need to be modelled as open quantum systems. While near-optimal algorithms have been developed for simulating many-body quantum dynamics, algorithms for their open system counterparts remain less well investigated. We address the problem of simulating geometrically local many-body open quantum systems interacting with a stationary Gaussian environment. Under a smoothness assumption on the system-environment interaction, we develop near-optimal algorithms that, for a model with NNN spins and evolution time ttt, attain a simulation error δδδ in the system-state with O(Nt(Nt/δ)o(1))\mathcal{O}(Nt(Nt/δ)^{o(1)})O(Nt(Nt/δ)o(1)) gates, O(t(Nt/δ)o(1))\mathcal{O}(t(Nt/δ)^{o(1)})O(t(Nt/δ)o(1)) parallelized circuit depth and O~(N(Nt/δ)o(1))\tilde{\mathcal{O}}(N(Nt/δ)^{o(1)})O~(N(Nt/δ)o(1)) ancillas. We additionally show that, if only simulating local observables is of interest, then the circuit depth of the digital algorithm can be chosen to be independent of the system size NNN. This provides theoretical evidence for the utility of these algorithms for simulating physically relevant models, where typically local observables are of interest, on pre-fault tolerant devices. Finally, for the limiting case of Markovian dynamics with commuting jump operators, we propose two algorithms based on sampling a Wiener process and on a locally dilated Hamiltonian construction, respectively. These algorithms reduce the asymptotic gate complexity on NNN compared to currently available algorithms in terms of the required number of geometrically local gates.