Saved in:
| Main Authors: | , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | https://arxiv.org/abs/2501.09607 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866918392393367552 |
|---|---|
| author | Etienney, Paul-Louis Robin, Rémi Rouchon, Pierre |
| author_facet | Etienney, Paul-Louis Robin, Rémi Rouchon, Pierre |
| contents | 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. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2501_09607 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | A posteriori error estimates for the Lindblad master equation Etienney, Paul-Louis Robin, Rémi Rouchon, Pierre Numerical Analysis Quantum Physics 65, 81 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. |
| title | A posteriori error estimates for the Lindblad master equation |
| topic | Numerical Analysis Quantum Physics 65, 81 |
| url | https://arxiv.org/abs/2501.09607 |