Analyzing the Positivity Preservation of Numerical Methods for the Liouville-von Neumann Equation

Abstract

The density matrix is a widely used tool in quantum mechanics. In order to determine its evolution with respect to time, the Liouville-von Neumann equation must be solved. However, analytic solutions of this differential equation exist only for simple cases. Additionally, if the equation is coupled to Maxwell's equations to model light-matter interaction, the resulting equation set -- the Maxwell-Bloch or Maxwell-Liouville-von Neumann (MLN) equations -- becomes nonlinear. In these advanced cases, numerical methods are required. Since the density matrix has certain mathematical properties, the numerical methods applied should be designed to preserve those properties. We establish the criterion that only methods that have a completely positive trace preserving (CPTP) update map can be used in long-term simulations. Subsequently, we assess the three most widely used methods -- the matrix exponential (ME) method, the Runge-Kutta (RK) method, and the predictor-corrector (PC) approach -- whether they provide this feature, and demonstrate that only the update step of the matrix exponential method is a CPTP map.

References in corpus (4)

• Modeling techniques for quantum cascade lasers
• Procedures for Converting among Lindblad, Kraus and Matrix Representations of Quantum Dynamical Semigroups
• Numerical studies of the interaction of an atomic sample with the electromagnetic field in two dimensions
• Dissipative Quantum Dynamics and Optimal Control using Iterative Time Ordering: An Application to Superconducting Qubits

Cited by in corpus (9)

• A short introduction to the Lindblad Master Equation
• Optoelectronic device simulations based on macroscopic Maxwell-Bloch equations
• Mbsolve: An open-source solver tool for the Maxwell-Bloch equations
• Optical Bloch modeling of few-cycle laser-induced electron dynamics in dielectrics
• Structure-preserving numerical schemes for Lindblad equations
• Dissipative protection of a GKP qubit in a high-impedance superconducting circuit driven by a microwave frequency comb
• Full- and low-rank exponential Euler integrators for the Lindblad equation
• Absence of a dissipative quantum phase transition in Josephson junctions: Theory
• Numerically efficient density-matrix technique for modeling electronic transport in midinfrared quantum cascade lasers