© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2501.09607 · Quantum

A posteriori error estimates for the Lindblad master equation

Paul-Louis Etienney, Rémi Robin, Pierre Rouchon

Published Mar 16, 2026• 0 citations

Abstract

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.

DOI: 10.22331/q-2026-03-16-2031

Related papers

Unconditionally stable time discretization of Lindblad master equations in infinite dimension using quantum channels

Rémi Robin, Pierre Rouchon, Lev-Arcady Sellem

We examine the time discretization of Lindblad master equations in infinite-dimensional Hilbert spaces. Our study is motivated by the fact that, with unbounded Lindbladian, projecting the evolution onto a finite-dimensional subspace using a Galerkin approximation inherently introduces stiffness, leading to a Courant--Friedrichs--Lewy type condition for explicit integration schemes. We propose and establish the convergence of a family of explicit numerical schemes for time discretization adapted to infinite dimension. These schemes correspond to quantum channels and thus preserve the physical properties of quantum evolutions on the set of density operators: linearity, complete positivity and trace. Numerical experiments inspired by bosonic quantum codes illustrate the practical interest of this approach when approximating the solution of infinite dimensional problems by that of finite dimensional problems of increasing dimension.

arXiv:2503.01712•Mar 3, 2025•0 citations

Error Bounds for Finite-Dimensional Approximations of Input-Output Open Quantum Systems by Subspace Truncation and Adiabatic Elimination

Onvaree Techakesari, Hendra I. Nurdin

An important class of physical systems that are of interest in practice are input-output open quantum systems that can be described by quantum stochastic differential equations and defined on an infinite-dimensional underlying Hilbert space. Most commonly, these systems involve coupling to a quantum harmonic oscillator as a system component. This paper is concerned with error bounds in the finite-dimensional approximations of input-output open quantum systems defined on an infinite-dimensional Hilbert space. We develop a framework for developing error bounds between the time evolution of the state of a class of infinite-dimensional quantum systems and its approximation on a finite-dimensional subspace of the original, when both are initialized in the latter subspace. This framework is then applied to two approaches for obtaining finite-dimensional approximations: subspace truncation and adiabatic elimination. Applications of the bounds to some physical examples drawn from the literature are provided to illustrate our results.

arXiv:1509.02629•Sep 9, 2015•0 citations

Convergence Analysis of Galerkin Approximations for the Lindblad Master Equation

Rémi Robin, Pierre Rouchon

This paper analyzes the numerical approximation of the Lindblad master equation on infinite-dimensional Hilbert spaces. We employ a classical Galerkin approach for spatial discretization and investigate the convergence of the discretized solution to the exact solution. Using \textit{a priori} estimates, we derive explicit convergence rates and demonstrate the effectiveness of our method through examples motivated by autonomous quantum error correction.

arXiv:2510.11416•Oct 13, 2025•0 citations

Error bounds on finite-dimensional approximations of input-output open quantum systems

Onvaree Techakesari, Hendra I. Nurdin

Many physical systems of interest that are encountered in practice are input-output open quantum systems described by quantum stochastic differential equations and defined on an infinite-dimensional underlying Hilbert space. Most commonly, these systems involve coupling to a quantum harmonic oscillator as a system component. This paper is concerned with the error in the finite-dimensional approximation of input-output open quantum systems defined on an infinite-dimensional underlying Hilbert space. We present explicit error bounds between the time evolution of the state of a class of infinite-dimensional quantum systems and its approximation on a finite-dimensional subspace of the original, when both are initialized in the latter subspace. Application to a physical example drawn from the literature is provided to illustrate our results.

IEEE Conference on Decision and Control (CDC)•Dec 1, 2015•3 citations

Generic method for integrating Lindblad master equations

Jiayin Gu, Fan Zhang

The time evolution of Markovian open quantum systems is governed by Lindblad master equations, whose solution can be formally written as the Lindbladian exponential acting on the initial density matrix. By expanding this Lindbladian exponential into the Taylor series, we propose a generic method for integrating Lindblad master equations. In this method, the series is truncated, retaining a finite number of terms, and the iterative actions of Lindbladian on the density matrix follow the corresponding master equation. Our method offers significant improvements in numerical efficiencies both in memory cost and computation time, especially for systems with many degrees of freedom. Moreover, our proposed method can be integrated seamlessly with tensor networks. Two illustrative examples, a two-level system exhibiting damped Rabi oscillations and a driven dissipative Heisenberg chain, are used to demonstrate the validity of our method. The superiority of our method is benchmarked with detailed performance tests.

arXiv:2412.13661•Physical review. E•Sep 8, 2025•1 citation

Simulating Open Quantum Systems Using Hamiltonian Simulations

Zhiyan Ding, Xiantao Li, Lin Lin

We present a novel method to simulate the Lindblad equation, drawing on the relationship between Lindblad dynamics, stochastic differential equations, and Hamiltonian simulations. We derive a sequence of unitary dynamics in an enlarged Hilbert space that can approximate the Lindblad dynamics up to an arbitrarily high order. This unitary representation can then be simulated using a quantum circuit that involves only Hamiltonian simulation and tracing out the ancilla qubits. There is no need for additional postselection in measurement outcomes, ensuring a success probability of one at each stage. Our method can be directly generalized to the time-dependent setting. We provide numerical examples that simulate both time-independent and time-dependent Lindbladian dynamics with accuracy up to the third order. Published by the American Physical Society 2024

arXiv:2311.15533•PRX Quantum•May 10, 2024•66 citations

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.

We examine the time discretization of Lindblad master equations in infinite-dimensional Hilbert spaces. Our study is motivated by the fact that, with unbounded Lindbladian, projecting the evolution onto a finite-dimensional subspace using a Galerkin approximation inherently introduces stiffness, leading to a Courant--Friedrichs--Lewy type condition for explicit integration schemes. We propose and establish the convergence of a family of explicit numerical schemes for time discretization adapted to infinite dimension. These schemes correspond to quantum channels and thus preserve the physical properties of quantum evolutions on the set of density operators: linearity, complete positivity and trace. Numerical experiments inspired by bosonic quantum codes illustrate the practical interest of this approach when approximating the solution of infinite dimensional problems by that of finite dimensional problems of increasing dimension.

An important class of physical systems that are of interest in practice are input-output open quantum systems that can be described by quantum stochastic differential equations and defined on an infinite-dimensional underlying Hilbert space. Most commonly, these systems involve coupling to a quantum harmonic oscillator as a system component. This paper is concerned with error bounds in the finite-dimensional approximations of input-output open quantum systems defined on an infinite-dimensional Hilbert space. We develop a framework for developing error bounds between the time evolution of the state of a class of infinite-dimensional quantum systems and its approximation on a finite-dimensional subspace of the original, when both are initialized in the latter subspace. This framework is then applied to two approaches for obtaining finite-dimensional approximations: subspace truncation and adiabatic elimination. Applications of the bounds to some physical examples drawn from the literature are provided to illustrate our results.

Many physical systems of interest that are encountered in practice are input-output open quantum systems described by quantum stochastic differential equations and defined on an infinite-dimensional underlying Hilbert space. Most commonly, these systems involve coupling to a quantum harmonic oscillator as a system component. This paper is concerned with the error in the finite-dimensional approximation of input-output open quantum systems defined on an infinite-dimensional underlying Hilbert space. We present explicit error bounds between the time evolution of the state of a class of infinite-dimensional quantum systems and its approximation on a finite-dimensional subspace of the original, when both are initialized in the latter subspace. Application to a physical example drawn from the literature is provided to illustrate our results.

The time evolution of Markovian open quantum systems is governed by Lindblad master equations, whose solution can be formally written as the Lindbladian exponential acting on the initial density matrix. By expanding this Lindbladian exponential into the Taylor series, we propose a generic method for integrating Lindblad master equations. In this method, the series is truncated, retaining a finite number of terms, and the iterative actions of Lindbladian on the density matrix follow the corresponding master equation. Our method offers significant improvements in numerical efficiencies both in memory cost and computation time, especially for systems with many degrees of freedom. Moreover, our proposed method can be integrated seamlessly with tensor networks. Two illustrative examples, a two-level system exhibiting damped Rabi oscillations and a driven dissipative Heisenberg chain, are used to demonstrate the validity of our method. The superiority of our method is benchmarked with detailed performance tests.

We present a novel method to simulate the Lindblad equation, drawing on the relationship between Lindblad dynamics, stochastic differential equations, and Hamiltonian simulations. We derive a sequence of unitary dynamics in an enlarged Hilbert space that can approximate the Lindblad dynamics up to an arbitrarily high order. This unitary representation can then be simulated using a quantum circuit that involves only Hamiltonian simulation and tracing out the ancilla qubits. There is no need for additional postselection in measurement outcomes, ensuring a success probability of one at each stage. Our method can be directly generalized to the time-dependent setting. We provide numerical examples that simulate both time-independent and time-dependent Lindbladian dynamics with accuracy up to the third order. Published by the American Physical Society 2024