Quantum Simulation of Lindbladian Dynamics via Repeated Interactions

The Lindblad equation generalizes the Schrödinger equation to quantum systems that undergo dissipative dynamics. The quantum simulation of Lindbladian dynamics is therefore non-unitary, preventing a naive application of state-of-the-art quantum algorithms. Here, we make use of an approximate correspondence between Lindbladian dynamics and evolution based on Repeated Interaction (RI) CPTP maps to write down a Hamiltonian formulation of the Lindblad dynamics and derive a rigorous error bound on the master equation. Specifically, we show that the number of interactions needed to simulate the Liouvillian e^{tmathcal{L}} within error ε scales in a weak coupling limit as νin O\(t²\|mathcal{L}\|_{1→ 1}²/ε\). This is significant because the error in the Lindbladian approximation to the dynamics is not explicitly bounded in existing quantum algorithms for open system simulations. We then provide quantum algorithms to simulate RI maps using an iterative Qubitization approach and Trotter-Suzuki formulas and specifically show that for iterative Qubitization the number of operations needed to simulate the dynamics (for a fixed value of ν) scales in a weak coupling limit as O\(α₀ t + νlog(1/ε\)/loglog(1/ε)) where α₀ is the coefficient 1-norm for the system and bath Hamiltonians. This scaling would appear to be optimal if the complexity of ν is not considered, which underscores the importance of considering the error in the Liouvillian that we reveal in this work.