pith. sign in

arxiv: 2604.22568 · v1 · submitted 2026-04-24 · 🪐 quant-ph · math-ph· math.MP

On truncations of hierarchical equations of motion for finite-dimensional systems

Vasilii Vadimov This is my paper

Pith reviewed 2026-05-08 12:08 UTC · model grok-4.3

classification 🪐 quant-ph math-phmath.MP
keywords hierarchical equations of motiontruncationsspectral convergenceopen quantum systemsSchur complementspectral pollutionfinite-dimensional systems
0
0 comments X

The pith

Finite-dimensional truncations of hierarchical equations of motion using a Schur-complement terminator converge in spectrum to the full equations and introduce no spurious unstable modes when the original system is stable.

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

The paper proves that approximations to the hierarchical equations of motion for finite-dimensional open quantum systems, obtained by truncating the hierarchy at finite depth with a Schur-complement terminator, have eigenvalues that approach those of the untruncated equations as depth increases. It further shows that deep truncations produce no artificial eigenvalues with positive real parts if the exact equations are stable. This matters for numerical work because these equations model dissipative quantum dynamics, and truncations must be accurate enough for reliable simulation of systems like spins interacting with environments without introducing numerical artifacts. The proofs apply specifically when the system Hilbert space is finite-dimensional and the terminator closes the hierarchy via a Schur-complement relation.

Core claim

For finite-dimensional open quantum systems, finite truncations of the hierarchical equations of motion constructed with a Schur-complement terminator have spectra that converge to the spectrum of the full HEOM as the truncation depth increases. These approximations are free of spectral pollution: if the exact HEOM is stable, then sufficiently deep truncations do not generate spurious unstable modes. The result is illustrated explicitly for the spin-boson model.

What carries the argument

The Schur-complement terminator, a closure relation that expresses higher-tier auxiliary operators in terms of lower-tier ones to produce a finite matrix representation of the HEOM.

If this is right

  • Eigenvalues of the dynamics can be obtained to arbitrary accuracy by solving finite matrix problems at large enough truncation depth.
  • Numerical integration of the truncated equations will not exhibit artificial exponential growth if the original system is stable.
  • Truncation depth can be increased systematically until the computed spectrum stabilizes.
  • The method supplies a rigorous justification for using these truncations in simulations of finite-dimensional open quantum systems.

Where Pith is reading between the lines

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

  • Similar convergence arguments might be developed for other closure relations beyond the Schur-complement type.
  • The result suggests monitoring the spectrum during numerical runs to confirm that truncation depth is sufficient to avoid any residual pollution.
  • Extensions could address how the convergence rate depends on the specific bath parameters or system dimension.

Load-bearing premise

The exact hierarchical equations of motion are stable, the quantum system is finite-dimensional, and the truncation uses a Schur-complement terminator.

What would settle it

Finding an eigenvalue with positive real part in a sufficiently deep Schur-complement truncation of a stable HEOM would falsify the no-spectral-pollution result.

Figures

Figures reproduced from arXiv: 2604.22568 by Vasilii Vadimov.

Figure 1
Figure 1. Figure 1: Eigenvalues of truncated Liouvillians LT(γ∗ ) . Here, α = 2, ω0 = 2, and η = 0.5. Poles and residues of coth approximation are given in Tab. 1. where σˆ x and σˆz are Pauli matrices. We choose the spectral density in the form J(ω) = αωω4 0 view at source ↗
read the original abstract

We study truncations of hierarchical equations of motion (HEOM) for finite-dimensional open quantum systems. We prove that for finite-dimensional approximations constructed with a Schur-complement type of terminator, the spectrum converges to that of the full HEOM as the truncation depth increases. We also prove that this approximation is free of spectral pollution: sufficiently deep truncations do not produce spurious unstable modes, provided the exact HEOM is stable. We illustrate the results for the spin-boson model.

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 / 3 minor

Summary. The paper proves that, for finite-dimensional open quantum systems, HEOM truncations constructed via Schur-complement terminators have spectra that converge to the spectrum of the untruncated HEOM as truncation depth increases. It further proves that these truncations introduce no spectral pollution (i.e., no spurious unstable eigenvalues appear for sufficiently deep truncations) provided the exact infinite HEOM is stable. The claims are illustrated numerically on the spin-boson model.

Significance. If the proofs are correct, the work supplies rigorous justification for a common class of HEOM truncations, directly addressing convergence and stability questions that arise in numerical applications to open quantum systems. The finite-dimensional setting permits an algebraic treatment via matrices and Schur complements, and the conditional result (no pollution when the exact hierarchy is stable) is a useful, falsifiable statement. The spin-boson illustration, while not the central contribution, provides a concrete check.

major comments (2)
  1. [§3] §3 (or the main theorem statement): the proof of spectral convergence relies on the finite-dimensional assumption to treat the hierarchy operators as matrices; it is not immediately clear whether the argument extends to the infinite-dimensional case or requires additional compactness arguments. This is load-bearing for the scope of the claim.
  2. [Theorem 2] Theorem on absence of spectral pollution: the stability hypothesis on the exact HEOM is assumed rather than derived; the manuscript should explicitly state whether this hypothesis can be verified a priori for typical bath correlation functions or must be checked numerically for each application.
minor comments (3)
  1. [Abstract] The abstract and introduction should clarify that the results are conditional on the exact HEOM being stable; this is already stated in the body but is easy to miss on first reading.
  2. [Figure 1] Figure 1 (spin-boson spectra): the caption should report the truncation depths used and the norm or distance metric employed to quantify convergence to the reference spectrum.
  3. [Eq. (X)] Notation: the definition of the Schur-complement terminator (Eq. (X)) should be cross-referenced in the statement of the main theorems to avoid ambiguity about which closure is being analyzed.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment and the recommendation for minor revision. We address the two major comments point by point below. Both points are constructive, and we will incorporate clarifications into the revised manuscript.

read point-by-point responses

  1. Referee: [§3] §3 (or the main theorem statement): the proof of spectral convergence relies on the finite-dimensional assumption to treat the hierarchy operators as matrices; it is not immediately clear whether the argument extends to the infinite-dimensional case or requires additional compactness arguments. This is load-bearing for the scope of the claim.

    Authors: We thank the referee for highlighting this point. The manuscript is deliberately restricted to finite-dimensional open quantum systems, as indicated by the title, abstract, and the opening paragraph of Section 1. This restriction permits an entirely algebraic treatment in which the (truncated) HEOM operators are finite matrices, allowing direct application of Schur-complement identities and standard matrix perturbation theory for the spectral convergence result (Theorem 1). We agree that the same argument does not immediately carry over to infinite-dimensional baths or system operators; additional functional-analytic ingredients such as compactness or sectorial properties would be needed. Because the paper makes no claim about the infinite-dimensional setting, we view the finite-dimensional hypothesis as part of the stated scope rather than a hidden limitation. In the revised manuscript we will add a short clarifying paragraph at the end of Section 3 that explicitly notes the reliance on finite dimensionality and states that extension to the infinite-dimensional case is left for future work. revision: partial

  2. Referee: [Theorem 2] Theorem on absence of spectral pollution: the stability hypothesis on the exact HEOM is assumed rather than derived; the manuscript should explicitly state whether this hypothesis can be verified a priori for typical bath correlation functions or must be checked numerically for each application.

    Authors: The referee correctly observes that Theorem 2 takes stability of the untruncated HEOM as an assumption rather than proving it. This is intentional: stability is a property of the infinite hierarchy and therefore depends on the concrete form of the bath correlation functions, the system-bath coupling strengths, and the temperature. For many standard spectral densities (Drude-Lorentz, underdamped Brownian, etc.) that produce exponentially decaying or finite-sum correlation functions, stability can be verified a priori by inspecting the spectrum of the auxiliary-mode operators or by appealing to existing results in the HEOM literature on the absence of unstable modes. For more general or numerically generated correlation functions, however, one must in practice compute or bound the eigenvalues of a sufficiently deep truncation and check that no eigenvalues with positive real part appear. In the revised manuscript we will expand the paragraph immediately following Theorem 2 to make this distinction explicit, giving brief guidance on when analytic verification is feasible and when numerical confirmation for the specific parameters is advisable. revision: yes

Circularity Check

0 steps flagged

No significant circularity in mathematical proofs

full rationale

The paper establishes mathematical proofs of spectral convergence and absence of spurious unstable modes for Schur-complement truncations of the HEOM as truncation depth increases, conditional on stability of the exact (infinite) HEOM and finite-dimensionality of the system. These results follow from algebraic properties of the Schur complement applied to the hierarchy operators, without any reduction to fitted parameters, self-definitional closures, or load-bearing self-citations. The derivation chain is self-contained as a conditional proof on matrix spectra and does not invoke any of the enumerated circularity patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work consists of mathematical proofs on operator spectra and truncations; it relies on standard results from linear algebra and functional analysis rather than new postulates or fitted quantities.

axioms (2)
  • standard math Schur complement preserves relevant spectral properties for the terminator construction
    Invoked in the definition of the finite-dimensional approximation.
  • domain assumption Stability of the exact infinite HEOM implies stability of sufficiently deep truncations
    Explicit condition stated for the no-spectral-pollution result.

pith-pipeline@v0.9.0 · 5366 in / 1272 out tokens · 52717 ms · 2026-05-08T12:08:26.046077+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

31 extracted references · 31 canonical work pages

  1. [1]

    Tanimura and R

    Y. Tanimura and R. Kubo, Time evolution of a quantum system in contact with a nearly gaussian-markoffian noise bath, Journal of the Physical Society of Japan 58(1), 101 (1989), doi:10.1143/jpsj.58.101

    work page doi:10.1143/jpsj.58.101 1989

  2. [2]

    Tanimura, Numerically “exact” approach to open quantum dynamics: The hierarchical equations of mo- tion (heom), The Journal of Chemical Physics153, 10.1063/5.0011599 (2020)

    Y. Tanimura, Numerically “exact” approach to open quantum dynamics: The hierarchical equations of motion (heom), The Journal of Chemical Physics 153(2) (2020), doi:10.1063/5.0011599

    work page doi:10.1063/5.0011599 2020

  3. [3]

    S. Bai, S. Zhang, C. Huang and Q. Shi, Hierarchical equations of motion for quantum chemical dynamics: Recent methodology developments and applications, Accounts of Chemical Research 57(21), 3151 (2024), doi:10.1021/acs.accounts.4c00492

    work page doi:10.1021/acs.accounts.4c00492 2024

  4. [4]

    H. Yang, Z. Sheng, L. Zhou and Z. Shuai, Hierarchical equations of motion in matrix product states: Formalism and applications for charge transport, Journal of Chemical Theory and Computation 21(20), 10108 (2025), doi:10.1021/acs.jctc.5c01385

    work page doi:10.1021/acs.jctc.5c01385 2025

  5. [5]

    L. Ye, X. Wang, D. Hou, R. Xu, X. Zheng and Y. Yan, Heom‐quick: a program for accurate, efficient, and universal characterization of strongly correlated quantum impurity systems, WIREs Computational Molecular Science 6(6), 608 (2016), doi:10.1002/wcms.1269

    work page doi:10.1002/wcms.1269 2016

  6. [6]

    X. Dan, M. Xu, J. T. Stockburger, J. Ankerhold and Q. Shi, Efficient low-temperature simulations for fermionic reservoirs with the hierarchical equations of motion method: Application to the anderson impurity model, Physical Review B 107(19), 195429 (2023), doi:10.1103/physrevb.107.195429

    work page doi:10.1103/physrevb.107.195429 2023

  7. [7]

    Strümpfer and K

    J. Strümpfer and K. Schulten, Open quantum dynamics calculations with the hierarchy equations of motion on parallel computers, Journal of Chemical Theory and Computation 8(8), 2808 (2012), doi:10.1021/ct3003833

    work page doi:10.1021/ct3003833 2012

  8. [8]

    QuTiP- BoFiN: A bosonic and fermionic numerical hierarchical-equations-of-motion library with applications in light-harvesting, quantum control, and single-molecule electronics,

    N. Lambert, T. Raheja, S. Cross, P. Menczel, S. Ahmed, A. Pitchford, D. Burgarth and F. Nori, Qutip-bofin: A bosonic and fermionic numerical hierarchical-equations-of-motion library with applications in light-harvesting, quantum control, and single-molecule electronics, Physical Review Research 5(1), 013181 (2023), doi:10.1103/physrevresearch.5.013181

    work page doi:10.1103/physrevresearch.5.013181 2023

  9. [9]

    Huang, P.-C

    Y.-T. Huang, P.-C. Kuo, N. Lambert, M. Cirio, S. Cross, S.-L. Yang, F. Nori and Y.-N. Chen, An efficient julia framework for hierarchical equations of motion in open quantum systems, Communications Physics 6(1) (2023), doi:10.1038/s42005-023-01427-2

    work page doi:10.1038/s42005-023-01427-2 2023

  10. [10]

    Zhang, L

    D. Zhang, L. Ye, J. Cao, Y. Wang, R. Xu, X. Zheng and Y. Yan, Heom‐quick2: A general‐purpose simulator for fermionic many‐body open quantum systems—an update, WIREs Computational Molecular Science 14(4) (2024), doi:10.1002/wcms.1727

    work page doi:10.1002/wcms.1727 2024

  11. [11]

    M. Xu, V. Vadimov, M. Krug, J. T. Stockburger and J. Ankerhold, A universal framework for quantum dissipation:minimally extended state space and exact time-local dynamics, doi:10.48550/ARXIV.2307.16790 (2023), 2307.16790

    work page doi:10.48550/arxiv.2307.16790 2023

  12. [12]

    M. Xu, V. Vadimov, J. T. Stockburger and J. Ankerhold, Colloquium: Simulating non-markovian dynamics in open quantum systems, Reviews of Modern Physics (2026), doi:10.1103/w3nw-hbjc

    work page doi:10.1103/w3nw-hbjc 2026

  13. [13]

    B. Witt, L. Rudnicki, Y. Tanimura and F. Mintert, Exploring complete positivity in hierarchy equations of motion, New Journal of Physics 19(1), 013007 (2017), doi:10.1088/1367-2630/19/1/013007

    work page doi:10.1088/1367-2630/19/1/013007 2017

  14. [14]

    M. Xu, Y. Yan, Q. Shi, J. Ankerhold and J. Stockburger, Taming quantum noise for efficient low temperature simulations of open quantum systems, Physical Review Letters 129(23), 230601 (2022), doi:10.1103/physrevlett.129.230601

    work page doi:10.1103/physrevlett.129.230601 2022

  15. [15]

    Uni- fied analysis of non-Markovian open quantum systems in Gaus- sian environment using superoperator formalism,

    Z. Huang, Y. Zhu, G. Park and L. Lin, Unified error bounds for perturbations of non-markovian open quantum systems in gaussian environments, doi:10.48550/ARXIV.2411.08741 (2024)

    work page doi:10.48550/arxiv.2411.08741 2024

  16. [16]

    Chen and I

    X. Chen and I. Franco, Bexcitonics: Quasiparticle approach to open quantum dynamics, The Journal of Chemical Physics 160(20) (2024), doi:10.1063/5.0198567

    work page doi:10.1063/5.0198567 2024

  17. [17]

    Q. Shi, Y. Xu, Y. Yan and M. Xu, Efficient propagation of the hierarchical equations of motion using the matrix product state method, The Journal of Chemical Physics 148(17) (2018), doi:10.1063/1.5026753

    work page doi:10.1063/1.5026753 2018

  18. [18]

    I. S. Dunn, R. Tempelaar and D. R. Reichman, Removing instabilities in the hierarchical equations of motion: Exact and approximate projection approaches, The Journal of Chemical Physics 150(18) (2019), doi:10.1063/1.5092616

    work page doi:10.1063/1.5092616 2019

  19. [19]

    Y. Yan, T. Xing and Q. Shi, A new method to improve the numerical stability of the hierarchical equations of motion for discrete harmonic oscillator modes, The Journal of Chemical Physics 153(20) (2020), doi:10.1063/5.0027962

    work page doi:10.1063/5.0027962 2020

  20. [20]

    T. Li, Y. Yan and Q. Shi, A low-temperature quantum fokker–planck equation that improves the numerical stability of the hierarchical equations of motion for the brownian oscillator spectral density, The Journal of Chemical Physics 156(6) (2022), doi:10.1063/5.0082108

    work page doi:10.1063/5.0082108 2022

  21. [21]

    Krug and J

    M. Krug and J. Stockburger, On stability issues of the heom method, The European Physical Journal Special Topics 232(20–22), 3219 (2023), doi:10.1140/epjs/s11734-023-00972-9

    work page doi:10.1140/epjs/s11734-023-00972-9 2023

  22. [22]

    S. A. Gershgorin, Uber die abgrenzung der eigenwerte einer matrix, Izvestiya Rossiyskoy akademii nauk. Seriya matematicheskaya (6), 749 (1931)

    work page 1931

  23. [23]

    R. S. Varga, Geršgorin and his circles, No. 36 in Springer series in computational mathematics. Springer, Berlin [u.a.], corr. 2. print. edn., ISBN 9783540211006 (2011)

    work page 2011

  24. [24]

    and Hasler, D

    M. Griesemer and D. Hasler, On the smooth feshbach–schur map, Journal of Functional Analysis 254(9), 2329 (2008), doi:10.1016/j.jfa.2008.01.015

    work page doi:10.1016/j.jfa.2008.01.015 2008

  25. [25]

    Li, S.-X

    X. Li, S.-X. Lyu, Y. Wang, R.-X. Xu, X. Zheng and Y. Yan, Towards quantum simulation of non-markovian open quantum dynamics: A universal and compact theory, doi:10.48550/ARXIV.2401.17255 (2024)

    work page doi:10.48550/arxiv.2401.17255 2024

  26. [26]

    Vadimov, M

    V. Vadimov, M. Xu, J. T. Stockburger, J. Ankerhold and M. Möttönen, Nonlinear-response theory for lossy superconducting quantum circuits, Physical Review Research 7(1), 013317 (2025), doi:10.1103/physrevresearch.7.013317

    work page doi:10.1103/physrevresearch.7.013317 2025

  27. [27]

    D. G. Feingold and R. S. Varga, Block diagonally dominant matrices and generalizations of the gerschgorin circle theorem., Pacific Journal of Mathematics 12(4), 1241 (1962)

    work page 1962

  28. [28]

    Varah, A lower bound for the smallest singular value of a matrix, Linear Algebra and its Applications 11(1), 3 (1975), doi:10.1016/0024-3795(75)90112-3

    J. Varah, A lower bound for the smallest singular value of a matrix, Linear Algebra and its Applications 11(1), 3 (1975), doi:10.1016/0024-3795(75)90112-3

    work page doi:10.1016/0024-3795(75)90112-3 1975

  29. [29]

    H. N. Salas, Gershgorin’s theorem for matrices of operators, Linear Algebra and its Applications 291(1–3), 15 (1999), doi:10.1016/s0024-3795(98)10219-7

    work page doi:10.1016/s0024-3795(98)10219-7 1999

  30. [30]

    M. J. Ablowitz and A. S. Fokas, Complex variables: introduction and applications, Cambridge university press (2003)

    work page 2003

  31. [31]

    Trefethen

    Y. Nakatsukasa, O. S \` e te and L. N. Trefethen, The AAA algorithm for rational approximation , SIAM Journal on Scientific Computing 40(3), A1494 (2018), doi:10.1137/16m1106122

    work page doi:10.1137/16m1106122 2018