A posteriori error estimates for the Lindblad master equation

We are interested in the simulation of open quantum systems governed by the Lindblad master equation in an infinite-dimensional Hilbert space. To simulate the solution of this equation, the standard approach involves two sequential approximations: first, we truncate the Hilbert space to derive a differential equation in a finite-dimensional subspace. Then, we use discrete time-step to obtain a numerical solution to the finite-dimensional evolution. In this paper, we establish bounds for these two approximations that can be explicitly computed to guarantee the accuracy of the numerical results. Through numerical examples, we demonstrate the efficiency of our method, empirically highlighting the tightness of the upper bound. While adaptive time-stepping is already a common practice in the time discretization of the Lindblad equation, we extend this approach by showing how to dynamically adjust the truncation of the Hilbert space. This enables fully adaptive simulations of the density matrix. For large-scale simulations, this approach can significantly reduce computational time and relieves users of the challenge of selecting an appropriate truncation. Furthermore, as a special case, our method naturally applies to Hamiltonian (unitary) dynamics.

💡 Research Summary

This paper addresses the numerical simulation of open quantum systems described by the Lindblad master equation in infinite‑dimensional Hilbert spaces. The authors focus on the two standard approximation steps: (i) truncating the infinite‑dimensional Hilbert space to a finite‑dimensional subspace, and (ii) discretising time with a numerical integrator. For each step they derive explicit a‑posteriori error bounds that can be evaluated using only the computed trajectory in the truncated space, thereby avoiding any dependence on the exact (inaccessible) solution.

The analysis begins by defining the truncated operators (H_N=P_NHP_N) and (\Gamma_{i,N}=P_N\Gamma_iP_N) and the corresponding finite‑dimensional Lindbladian (L_N). Using Duhamel’s formula and the contractivity of completely positive trace‑preserving (CPTP) maps in the trace norm, they prove the key inequality
\