Large-scale stochastic simulation of open quantum systems
arXiv:2501.17913 · doi:10.1038/s41467-025-66846-x
Abstract
Understanding the precise interaction mechanisms between quantum systems and their environment is crucial for advancing stable quantum technologies, designing reliable experimental frameworks, and building accurate models of real-world phenomena. However, simulating open quantum systems, which feature complex non-unitary dynamics, poses significant computational challenges that require innovative methods to overcome. In this work, we introduce the tensor jump method (TJM), a scalable, embarrassingly parallel algorithm for stochastically simulating large-scale open quantum systems, specifically Markovian dynamics captured by Lindbladians. This method is built on three core principles where, in particular, we extend the Monte Carlo wave function (MCWF) method to matrix product states, use a dynamic time-dependent variational principle (TDVP) to significantly reduce errors during time evolution, and introduce what we call a sampling MPS to drastically reduce the dependence on the simulation's time step size. We demonstrate that this method scales more effectively than previous methods and ensures convergence to the Lindbladian solution independent of system size, which we show both rigorously and numerically. Finally, we provide evidence of its utility by simulating Lindbladian dynamics of XXX Heisenberg models up to a thousand spins using a consumer-grade CPU. This work represents a significant step forward in the simulation of large-scale open quantum systems, with the potential to enable discoveries across various domains of quantum physics, particularly those where the environment plays a fundamental role, and to both dequantize and facilitate the development of more stable quantum hardware.
24 pages, 13 figures, 1 table (includes Methods and Appendix)
References in corpus (45)
- The density-matrix renormalization group in the age of matrix product states
- Area laws for the entanglement entropy - a review
- QuTiP 2: A Python framework for the dynamics of open quantum systems
- QuTiP: An open-source Python framework for the dynamics of open quantum systems
- Efficient simulation of one-dimensional quantum many-body systems
- Real time evolution using the density matrix renormalization group
- Time-dependent density-matrix renormalization-group using adaptive effective Hilbert spaces
- Matrix Product Density Operators: Simulation of finite-T and dissipative systems
- Probing the relaxation towards equilibrium in an isolated strongly correlated 1D Bose gas
- Quantum trajectories and open many-body quantum systems
- Time-dependent variational principle for quantum lattices
- Matrix Product States and Projected Entangled Pair States: Concepts, Symmetries, and Theorems
- Unifying time evolution and optimization with matrix product states
- Matrix product states represent ground states faithfully
- Time-evolution methods for matrix-product states
- Mixed-state dynamics in one-dimensional quantum lattice systems: a time-dependent superoperator renormalization algorithm
- A short introduction to the Lindblad Master Equation
- Tensor networks for complex quantum systems
- Open XXZ spin chain: Nonequilibrium steady state and strict bound on ballistic transport
- Exact relaxation in a class of non-equilibrium quantum lattice systems
- Hand-waving and Interpretive Dance: An Introductory Course on Tensor Networks
- Matrix product operator representations
- Non-Markovian quantum processes: complete framework and efficient characterisation
- Time-evolving a matrix product state with long-ranged interactions
- Exact nonequilibrium steady state of a strongly driven open XXZ chain
- Algebraic versus exponential decoherence in dissipative many-particle systems
- Simulation methods for open quantum many-body systems
- Post-Matrix Product State Methods: To tangent space and beyond
- Variational Matrix Product Operators for the Steady State of Dissipative Quantum Systems
- A positive tensor network approach for simulating open quantum many-body systems
- Non-equilibrium dynamics of bosonic atoms in optical lattices: Decoherence of many-body states due to spontaneous emission
- A Strictly Single-Site DMRG Algorithm with Subspace Expansion
- Efficient tensor network simulation of IBM's Eagle kicked Ising experiment
- On entropy growth and the hardness of simulating time evolution
- Thermofield-based chain mapping approach for open quantum systems
- Matrix product operators and states: NP-hardness and undecidability
- One-dimensional many-body entangled open quantum systems with tensor network methods
- Unifying variational methods for simulating quantum many-body systems
- The MQT Handbook: A Summary of Design Automation Tools and Software for Quantum Computing
- Simulating Noisy Quantum Circuits with Matrix Product Density Operators
- Entanglement-Optimal Trajectories of Many-Body Quantum Markov Processes
- Efficient sampling of noisy shallow circuits via monitored unraveling
- Stable bipolarons in open quantum systems
- Optimized trajectory unraveling for classical simulation of noisy quantum dynamics
- A Riemannian Approach to the Lindbladian Dynamics of a Locally Purified Tensor Network