Full- and low-rank exponential Euler integrators for the Lindblad equation

Abstract

The Lindblad equation is a widely used quantum master equation to model the dynamical evolution of open quantum systems whose states are described by density matrices. These solution matrices are characterized by semi-positiveness and trace preserving properties, which must be guaranteed in any physically meaningful numerical simulation. In this paper, novel full- and low-rank exponential Euler integrators are developed for approximating the Lindblad equation that preserve positivity and trace unconditionally. Theoretical results are presented that provide sharp error estimates for the two classes of exponential integration methods. Results of numerical experiments are discussed that illustrate the effectiveness of the proposed schemes, beyond present state-of-the-art capabilities.

References in corpus (14)

• The density-matrix renormalization group
• Efficient classical simulation of slightly entangled quantum computations
• QuTiP 2: A Python framework for the dynamics of open quantum systems
• A Practical Introduction to Tensor Networks: Matrix Product States and Projected Entangled Pair States
• Efficient simulation of one-dimensional quantum many-body systems
• Neural-Network Approach to Dissipative Quantum Many-Body Dynamics
• Simulation methods for open quantum many-body systems
• Variational Quantum Monte Carlo Method with a Neural-Network Ansatz for Open Quantum Systems
• A positive tensor network approach for simulating open quantum many-body systems
• Low rank approximation for the numerical simulation of high dimensional Lindblad and Riccati equations
• Analyzing the Positivity Preservation of Numerical Methods for the Liouville-von Neumann Equation
• Entanglement of formation of mixed many-body quantum states via Tree Tensor Operators
• Structure-preserving numerical schemes for Lindblad equations
• Adaptive low-rank approximation and denoised Monte-Carlo approach for high-dimensional Lindblad equations