Skip to content

Well-conditioned iterative methods for large open quantum systems

Gaspard Beugnot, Paul Gregory, Rémi Robin, Antoine Tilloy

quant-pharXiv:2608.30860

Abstract

Markovian open quantum systems are well modeled by the Lindblad Master Equation (ME) dd t ρt = L ρt, where L is a linear (super-)operator and ρt is the system state, a positive matrix. When designing or characterizing a quantum system, one is usually interested in the steady state ρ∞ (such that L ρ∞ = 0), the first few excited states, and trajectories t ρt. In finite dimension, ρt is an n× n matrix, L thus typically costs n4 to store explicitly as a dense matrix, and O(n6) to diagonalize or invert exactly, making standard linear algebraic techniques expensive for large systems. However, L usually costs only O(n3) to apply. This makes iterative methods appealing, but they do not work without a good preconditioner. In this article, our main observation is that a part of the Lindblad equation, corresponding to the so-called no-jump evolution S, can be inverted efficiently. Using this inverse map, we introduce an auxiliary completely positive trace-preserving (CPTP) map Φ whose fixed point is directly related to ρ∞, all the other eigenvalues having smaller magnitude. The map Φ is thus well suited to iterative methods, and ρ∞ can be found in a few Arnoldi iterations. Using the same inverse map S-1 as preconditioner, we compute the low-lying spectrum efficiently via shift-invert Arnoldi, and, as a proof of concept, build an implicit time integrator that is competitive on stiff systems in the low-precision regime. For the steady-state and low excited states problems, our methods scale like O(n3) per iteration and offer state-of-the-art performance on CPU and GPU.

Create a lesson