pith. sign in

arxiv: 2408.13601 · v1 · submitted 2024-08-24 · 🧮 math.NA · cs.NA· physics.comp-ph· quant-ph

Full- and low-rank exponential Euler integrators for the Lindblad equation

Hao Chen , Alfio Borz\`i , Denis Jankovi\'c , Jean-Gabriel Hartmann , Paul-Antoine Hervieux This is my paper

Pith reviewed 2026-05-23 21:19 UTC · model grok-4.3

classification 🧮 math.NA cs.NAphysics.comp-phquant-ph
keywords Lindblad equationexponential integratorspositivity preservationtrace preservationlow-rank methodsopen quantum systemsnumerical ODE methods
0
0 comments X

The pith

Exponential Euler integrators for the Lindblad equation preserve positivity and trace for any step size.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The Lindblad equation governs the evolution of open quantum systems via density matrices that must remain positive semi-definite with unit trace. This paper constructs full-rank and low-rank exponential Euler integrators that enforce these properties unconditionally, independent of chosen time step. The schemes rely on exponential approximations of the Lindblad superoperator whose flows stay inside the cone of valid states. Sharp a priori error estimates are derived for both integrator classes. Numerical tests show the methods maintain physical fidelity where earlier approaches fail.

Core claim

The authors develop novel full- and low-rank exponential Euler integrators for approximating the Lindblad equation. These integrators are built so that they map the set of positive semi-definite trace-one matrices into itself for arbitrary positive step sizes. Theoretical results supply sharp error estimates in appropriate matrix norms for the full-rank and low-rank families.

What carries the argument

Exponential Euler integrators obtained by splitting or approximating the Lindblad superoperator so the resulting flow automatically respects the positive semi-definite trace-one cone.

If this is right

  • Simulations of open quantum systems can use arbitrarily large time steps without violating positivity or trace.
  • Low-rank variants reduce storage and arithmetic cost while retaining the preservation property.
  • Sharp error bounds allow direct comparison of accuracy across different step sizes and ranks.
  • The schemes extend reliable long-time integration beyond the range of existing methods.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same splitting idea could be tested on other quantum master equations that share the Lindblad structure.
  • Low-rank preservation might be combined with adaptive rank control to handle systems whose effective dimension changes over time.
  • The unconditional stability suggests these integrators could serve as building blocks inside higher-order or implicit-explicit composite schemes.

Load-bearing premise

The Lindblad operator admits an exponential splitting or approximation whose flow maps valid density matrices to valid density matrices for any step size.

What would settle it

A concrete Lindblad operator and initial density matrix for which the integrator produces a matrix with a negative eigenvalue or trace different from one after one step.

Figures

Figures reproduced from arXiv: 2408.13601 by Alfio Borz\`i, Denis Jankovi\'c, Hao Chen, Jean-Gabriel Hartmann, Paul-Antoine Hervieux.

Figure 6.1
Figure 6.1. Figure 6.1: Numerical results of the FREE integrator for the Lindblad equation with Hamiltonian (6.1) [PITH_FULL_IMAGE:figures/full_fig_p018_6_1.png] view at source ↗
Figure 6.2
Figure 6.2. Figure 6.2: Numerical results of the FREE integrator for the Lindblad equation with time-dependent Hamil [PITH_FULL_IMAGE:figures/full_fig_p019_6_2.png] view at source ↗
Figure 6.3
Figure 6.3. Figure 6.3: Numerical results of the LREE scheme with [PITH_FULL_IMAGE:figures/full_fig_p020_6_3.png] view at source ↗
Figure 6.4
Figure 6.4. Figure 6.4: Numerical results of the LREE scheme with [PITH_FULL_IMAGE:figures/full_fig_p021_6_4.png] view at source ↗
Figure 6.5
Figure 6.5. Figure 6.5: Numerical results of the LREE scheme with [PITH_FULL_IMAGE:figures/full_fig_p022_6_5.png] view at source ↗
Figure 6.6
Figure 6.6. Figure 6.6: Ce vs τ . Left: time-independent case with d = 4, K = 4, a = 1.5, b = 0.5, γk = γ = 0.01, gkl(t) ≡ g = 1. Right: time-dependent case with d = 6, K = 3, a = 1, b = 1, γk = γ = 0.05, gkl(t) = δk,l−1 · (1 + t) 1 4 . 22 [PITH_FULL_IMAGE:figures/full_fig_p022_6_6.png] view at source ↗
Figure 6.7
Figure 6.7. Figure 6.7: Numerical comparison between the proposed exponential schemes and the solvers in QuTip for [PITH_FULL_IMAGE:figures/full_fig_p024_6_7.png] view at source ↗
Figure 6.8
Figure 6.8. Figure 6.8: Numerical comparison between the proposed exponential schemes and the solvers in QuTip for [PITH_FULL_IMAGE:figures/full_fig_p025_6_8.png] view at source ↗
Figure 6.9
Figure 6.9. Figure 6.9: Numerical comparison between the proposed exponential schemes and the solvers in QuTip for [PITH_FULL_IMAGE:figures/full_fig_p025_6_9.png] view at source ↗
Figure 6.10
Figure 6.10. Figure 6.10: Numerical comparison between the proposed exponential schemes and [PITH_FULL_IMAGE:figures/full_fig_p026_6_10.png] view at source ↗
Figure 6.11
Figure 6.11. Figure 6.11: Numerical comparison between the proposed exponential schemes and [PITH_FULL_IMAGE:figures/full_fig_p026_6_11.png] view at source ↗
read the original abstract

The Lindblad equation is a widely used quantum master equation to model the dynamical evolution of open quantum systems whose states are described by density matrices. These solution matrices are characterized by semi-positiveness and trace preserving properties, which must be guaranteed in any physically meaningful numerical simulation. In this paper, novel full- and low-rank exponential Euler integrators are developed for approximating the Lindblad equation that preserve positivity and trace unconditionally. Theoretical results are presented that provide sharp error estimates for the two classes of exponential integration methods. Results of numerical experiments are discussed that illustrate the effectiveness of the proposed schemes, beyond present state-of-the-art capabilities.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript develops novel full-rank and low-rank exponential Euler integrators for the Lindblad master equation. These schemes are asserted to preserve positivity and trace of the density matrix unconditionally for arbitrary step sizes. Sharp a priori error estimates are derived for both integrator classes, and numerical experiments are presented to demonstrate effectiveness beyond current state-of-the-art methods.

Significance. If the unconditional cone preservation and the sharpness of the error bounds can be rigorously established from the Lindblad generator, the work would supply practical, structure-preserving integrators for open quantum systems, particularly valuable for long-time simulations where step-size restrictions are prohibitive. The low-rank variant could additionally reduce computational cost in high-dimensional settings.

major comments (2)
  1. [§3] §3 (exponential Euler step definition): the unconditional positivity and trace preservation is asserted via an exponential splitting, yet no derivation is supplied showing that the effective map sends the PSD trace-1 cone into itself for every h>0; this property is load-bearing for the central claim and must be proved from the complete-positivity of the Lindblad super-operator rather than invoked.
  2. [§4] §4 (error analysis): the sharpness of the error estimates is stated, but the constants in the bounds appear to depend on post-hoc choices in the splitting that are not shown to be independent of the step size or the particular Lindblad coefficients; this undermines the 'sharp' qualifier.
minor comments (2)
  1. [§2] Notation for the low-rank factorization is introduced without a clear distinction from the full-rank case; a dedicated paragraph or table comparing the two would improve readability.
  2. [§5] Figure captions lack explicit mention of the matrix dimensions or the specific Lindblad operators used in the experiments.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We address the two major comments point by point below.

read point-by-point responses

  1. Referee: [§3] §3 (exponential Euler step definition): the unconditional positivity and trace preservation is asserted via an exponential splitting, yet no derivation is supplied showing that the effective map sends the PSD trace-1 cone into itself for every h>0; this property is load-bearing for the central claim and must be proved from the complete-positivity of the Lindblad super-operator rather than invoked.

    Authors: We agree that an explicit derivation is required. The current manuscript invokes the preservation via the splitting but does not supply the full step-by-step argument from complete positivity. In the revised manuscript we will insert a self-contained proof in §3 showing that the effective map sends the PSD trace-1 cone into itself for arbitrary h>0, using the fact that the Lindblad generator generates a completely positive trace-preserving semigroup. This directly addresses the load-bearing claim. revision: yes

  2. Referee: [§4] §4 (error analysis): the sharpness of the error estimates is stated, but the constants in the bounds appear to depend on post-hoc choices in the splitting that are not shown to be independent of the step size or the particular Lindblad coefficients; this undermines the 'sharp' qualifier.

    Authors: The constants arise from the local truncation error of the exponential Euler method and are independent of h by construction (they depend only on the fixed Lipschitz constant of the generator). Dependence on specific Lindblad coefficients is likewise controlled by the generator norm under the standard Lindblad assumptions. To remove any ambiguity about post-hoc choices we will add an explicit remark in the revised §4 stating these independence properties and confirming that the order is optimal. We therefore view the estimates as sharp once this clarification is supplied. revision: partial

Circularity Check

0 steps flagged

No significant circularity; derivations appear self-contained

full rationale

The paper presents exponential Euler integrators derived from the Lindblad equation structure, with claims of unconditional positivity and trace preservation tied to the exponential splitting property. No quoted steps reduce by construction to fitted inputs, self-definitions, or self-citation chains that bear the central load. Error estimates are presented as theoretical results from the method analysis, and the abstract and available context show no renaming of known results or ansatz smuggling. The cone-preservation property is invoked from the equation's form rather than derived tautologically within the paper's own equations.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Based on abstract only; no explicit free parameters, axioms, or invented entities are stated. The central claim rests on the mathematical structure of the Lindblad equation allowing unconditional preservation via exponential Euler steps, but the precise assumptions are not enumerated.

pith-pipeline@v0.9.0 · 5653 in / 1119 out tokens · 17941 ms · 2026-05-23T21:19:30.565842+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages

  1. [1]

    Breuer and F

    H.-P. Breuer and F. Petruccione. The Theory of Open Quantum Systems . Oxford University Press, 2007

    work page 2007

  2. [2]

    E.B. Davies. Quantum Theory of Open Systems . Academic Press, 1976

    work page 1976

  3. [3]

    Gorini, A

    V. Gorini, A. Kossakowski, and E.C.G. Sudarshan. Completely positive dynamical semigroups of N- level systems. Journal of Mathematical Physics , 17(5):821–825, 05 1976

    work page 1976

  4. [4]

    Lindblad

    G. Lindblad. On the generators of quantum dynamical semigroups. Communications in Mathematical Physics, 48(2):119–130, 1976

    work page 1976

  5. [5]

    Riesch and C

    M. Riesch and C. Jirauschek. Analyzing the positivity preservation of numerical methods for the Liouville-von Neumann equation. Journal of Computational Physics , 390:290–296, 2019

    work page 2019

  6. [6]

    Ziolkowski, J.M

    R.W. Ziolkowski, J.M. Arnold, and D.M. Gogny. Ultrafast pulse interactions with two-level atoms. Phys. Rev. A , 52:3082–3094, Oct 1995

    work page 1995

  7. [7]

    Bid´ egaray, A

    B. Bid´ egaray, A. Bourgeade, and D. Reignier. Introducing physical relaxation terms in Bloch equations. Journal of Computational Physics , 170(2):603–613, 2001

    work page 2001

  8. [8]

    Blanes, A

    S. Blanes, A. Iserles, and S. Macnamara. Positivity-preserving methods for ordinary differential equa- tions. ESAIM: M2AN, 56(6):1843–1870, 2022

    work page 2022

  9. [9]

    Colonna, and F

    A Martiradonna, G. Colonna, and F. Diele. Geco: Geometric conservative nonstandard schemes for biochemical systems. Applied Numerical Mathematics, 155:38–57, 2020. Structural Dynamical Systems: Computational Aspects held in Monopoli (Italy) on June 12-15, 2018

    work page 2020

  10. [10]

    Riesch, A

    M. Riesch, A. Pikl, and C. Jirauschek. Completely positive trace preserving methods for the Lind- blad equation. In 2020 International Conference on Numerical Simulation of Optoelectronic Devices (NUSOD), pages 109–110, 2020

    work page 2020

  11. [11]

    Bourgeade and O

    A. Bourgeade and O. Saut. Numerical methods for the bidimensional Maxwell-Bloch equations in nonlinear crystals. Journal of Computational Physics , 213(2):823–843, 2006

    work page 2006

  12. [12]

    Songolo and B

    M.E. Songolo and B. Bid´ egaray-Fesquet. Strang splitting schemes for N-level Bloch models. Interna- tional Journal of Modeling, Simulation, and Scientific Computing , 14(03):2350044, 2023. 27

    work page 2023

  13. [13]

    Songolo and B

    M.E. Songolo and B. Bid´ egaray-Fesquet. Nonstandard finite-difference schemes for the two-level Bloch model. International Journal of Modeling, Simulation, and Scientific Computing , 09(04):1850033, 2018

    work page 2018

  14. [14]

    Le Bris and P

    C. Le Bris and P. Rouchon. Low-rank numerical approximations for high-dimensional Lindblad equa- tions. Phys. Rev. A , 87:022125, Feb 2013

    work page 2013

  15. [15]

    Le Bris, P

    C. Le Bris, P. Rouchon, and J. Roussel. Adaptive low-rank approximation and denoised Monte Carlo approach for high-dimensional Lindblad equations. Phys. Rev. A , 92:062126, Dec 2015

    work page 2015

  16. [16]

    Cao and J

    Y. Cao and J. Lu. Structure-preserving numerical schemes for Lindblad equations, 2024

    work page 2024

  17. [17]

    Schlimgen, K

    A.W. Schlimgen, K. Head-Marsden, L.M. Sager, P. Narang, and D.A. Mazziotti. Quantum simulation of the Lindblad equation using a unitary decomposition of operators. Phys. Rev. Res. , 4:023216, Jun 2022

    work page 2022

  18. [18]

    Weimer, A

    H. Weimer, A. Kshetrimayum, and R. Or´ us. Simulation methods for open quantum many-body systems. Rev. Mod. Phys., 93:015008, Mar 2021

    work page 2021

  19. [19]

    Werner, D

    A.H. Werner, D. Jaschke, P. Silvi, M. Kliesch, T. Calarco, J. Eisert, and S. Montangero. Positive tensor network approach for simulating open quantum many-body systems. Phys. Rev. Lett., 116:237201, Jun 2016

    work page 2016

  20. [20]

    R. Or´ us. A practical introduction to tensor networks: Matrix product states and projected entangled pair states. Annals of Physics , 349:117–158, 2014

    work page 2014

  21. [21]

    Schollw¨ ock

    U. Schollw¨ ock. The density-matrix renormalization group. Rev. Mod. Phys., 77:259–315, Apr 2005

    work page 2005

  22. [22]

    G. Vidal. Efficient classical simulation of slightly entangled quantum computations. Phys. Rev. Lett., 91:147902, Oct 2003

    work page 2003

  23. [23]

    Efficient simulation of one-dimensional quantum many-body systems

    Guifr´ e Vidal. Efficient simulation of one-dimensional quantum many-body systems. Phys. Rev. Lett. , 93:040502, Jul 2004

    work page 2004

  24. [24]

    Hartmann and G

    M.J. Hartmann and G. Carleo. Neural-network approach to dissipative quantum many-body dynamics. Phys. Rev. Lett., 122:250502, Jun 2019

    work page 2019

  25. [25]

    Nagy and V

    A. Nagy and V. Savona. Variational quantum Monte Carlo method with a neural-network ansatz for open quantum systems. Phys. Rev. Lett., 122:250501, Jun 2019

    work page 2019

  26. [26]

    M. Reh, M. Schmitt, and M. G¨ arttner. Time-dependent variational principle for open quantum systems with artificial neural networks. Phys. Rev. Lett., 127:230501, Dec 2021

    work page 2021

  27. [27]

    Arceci, P

    L. Arceci, P. Silvi, and S. Montangero. Entanglement of formation of mixed many-body quantum states via tree tensor operators. Phys. Rev. Lett., 128:040501, Jan 2022

    work page 2022

  28. [28]

    D. Sulz, C. Lubich, G. Ceruti, I. Lesanovsky, and F. Carollo. Numerical simulation of long-range open quantum many-body dynamics with tree tensor networks. Phys. Rev. A , 109:022420, Feb 2024

    work page 2024

  29. [29]

    Chen and A

    H. Chen and A. Borz` ı. Positivity preserving exponential integrators for differential Riccati equations. Journal of Scientific Computing , 96(2):50, 2023

    work page 2023

  30. [30]

    Chen and A

    H. Chen and A. Borz` ı. Low-rank exponential integrators for stiff differential Riccati equations. Sub- mitted, 2024

    work page 2024

  31. [31]

    Nation, and F

    J.R Johansson, P.D. Nation, and F. Nori. Qutip 2: A Python framework for the dynamics of open quantum systems. Computer Physics Communications , 184(4):1234–1240, 2013

    work page 2013

  32. [32]

    Hochbruck and A

    M. Hochbruck and A. Ostermann. Exponential integrators. Acta Numerica, 19:209–286, 2010

    work page 2010

  33. [33]

    N. Lang, H. Mena, and J. Saak. On the benefits of the ldlt factorization for large-scale differential matrix equation solvers. Linear Algebra and its Applications , 480:44–71, 2015

    work page 2015

  34. [34]

    Al-Mohy and N.J

    A.H. Al-Mohy and N.J. Higham. A new scaling and squaring algorithm for the matrix exponential. SIAM Journal on Matrix Analysis and Applications , 31(3):970–989, 2010

    work page 2010

  35. [35]

    Al-Mohy and N.J

    A.H. Al-Mohy and N.J. Higham. Computing the action of the matrix exponential, with an application to exponential integrators. SIAM Journal on Scientific Computing , 33(2):488–511, 2011

    work page 2011

  36. [36]

    Gajic and M.T.J

    Z. Gajic and M.T.J. Qureshi. Lyapunov Matrix Equation in System Stability and Control . Dover Books on Engineering. Dover Publications, 2008

    work page 2008

  37. [37]

    Hamilton

    J.D. Hamilton. Time Series Analysis . Princeton University Press, 2020

    work page 2020

  38. [38]

    Hairer, S.P

    E. Hairer, S.P. Nørsett, and G. Wanner. Solving Ordinary Differential Equations I: Nonstiff Problems . Springer Series in Computational Mathematics. Springer Berlin Heidelberg, 2008

    work page 2008

  39. [39]

    J. H. Verner. Explicit Runge-Kutta methods with estimates of the local truncation error. SIAM Journal 28 on Numerical Analysis, 15(4):772–790, 1978. 29

    work page 1978