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 (16)
- The density-matrix renormalization group in the age of matrix product states
- QuTiP 2: A Python framework for the dynamics of open quantum systems
- Real time evolution using the density matrix renormalization group
- Matrix Product Density Operators: Simulation of finite-T and dissipative systems
- Quantum trajectories and open many-body quantum systems
- Matrix product states represent ground states faithfully
- Open XXZ spin chain: Nonequilibrium steady state and strict bound on ballistic transport
- Exact relaxation in a class of non-equilibrium quantum lattice systems
- Time-evolving a matrix product state with long-ranged interactions
- Variational Matrix Product Operators for the Steady State of Dissipative Quantum 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
- On entropy growth and the hardness of simulating time evolution
- Matrix product operators and states: NP-hardness and undecidability
- Unifying variational methods for simulating quantum many-body systems
- A Riemannian Approach to the Lindbladian Dynamics of a Locally Purified Tensor Network