Well-conditioned iterative methods for large open quantum systems
Gaspard Beugnot, Paul Gregory, Rémi Robin, Antoine Tilloy
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
Related papers
Continuous variable distributed quantum sensing in integrated photonics
Bethany Puzio, Oliver M. Green, Joel F. Tasker et al.
Securing quantum error correction against misleading advice from AI agents
A. Barış Özgüler
Exact logical error rates for magic state cultivation
Kwok Ho Wan, Ainhoa Zapirain
Hamiltonian engineering via pulses: beyond group averaging
Ivan Beschastnyi, Lucah Patel, David Tinoco
Logarithmic-depth quantum simulation of boson sampling
Changhun Oh
Entanglement swapping across a five-node relay in a multiplexed quantum-classical network
Andrew R. Cameron, Jordan M. Thomas, Alexandru Macridin et al.