Pith. sign in

REVIEW 2 major objections 5 minor 83 references

Complex-pole HEOM hides a continuous gauge that can be tuned to stabilize the Liouvillian without changing the exact dynamics.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-11 12:47 UTC pith:7TGCUHTO

load-bearing objection Clean, usable gauge freedom in complex-pole HEOM that actually removes two distinct numerical instabilities and reaches previously inaccessible regimes. the 2 major comments →

arxiv 2607.04834 v1 pith:7TGCUHTO submitted 2026-07-06 physics.chem-ph quant-ph

Hidden Gauge Freedom in Complex-Pole Hierarchical Equations of Motion

Tianchu Li , Andr\'es Montoya-Castillo This is my paper
classification physics.chem-ph quant-ph
keywords hierarchical equations of motionopen quantum systemscomplex-pole decompositiongauge freedomnon-normalityspin-boson modelsub-Ohmic bathBrownian oscillator
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Complex-pole hierarchical equations of motion (HEOM) give numerically exact open-system dynamics for arbitrary baths, but their non-Hermitian Liouvillians often produce unstable eigenvalues or amplify numerical errors. This paper shows that those Liouvillians are not unique: a continuous family of analytically equivalent generators all encode the same bath correlation function yet have radically different spectra and degrees of non-normality. The family is parameterized by a gauge angle that can be optimized once, yielding GO-HEOM. The optimized generator removes spectral divergences for strongly coupled Brownian oscillators and suppresses long-time error growth for sub-Ohmic baths, allowing stable simulations through the delocalized-to-localized quantum phase transition at coupling strengths previously inaccessible to HEOM. Because the gauge is independent of how the bath correlation is decomposed, the same fix applies to existing Prony, AAA and Padé implementations.

Core claim

Complex-pole HEOM possesses a previously unknown continuous gauge freedom: a one-parameter family of Liouvillians, all generating identical exact reduced dynamics for a given bath correlation function, whose eigenspectra and non-normality can be continuously reshaped. Optimizing a global scaling of that gauge produces GO-HEOM, which eliminates positive-real eigenvalues for strong Brownian-oscillator coupling and damps non-normal transient error growth for sub-Ohmic environments.

What carries the argument

The ϕ-dependent gauge transformation of the auxiliary ladder operators (Eq. 7–8 of the main text), which rewrites the Fokker–Planck generator while leaving the bath correlation function invariant and produces the family of HEOM equations (Eq. 10).

Load-bearing premise

That a single global scaling of the gauge angle applied uniformly to every complex pole is already near-optimal for multi-mode baths.

What would settle it

For a fixed multi-pole decomposition of a strongly coupled Brownian or sub-Ohmic bath, show that independent per-mode optimization of each ϕ_k (or the full four-parameter linear map) yields a generator whose largest real eigenvalue or long-time error growth is substantially smaller than the single-parameter GO-HEOM result reported in the paper.

Watch this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 5 minor

Summary. The manuscript reports a continuous gauge freedom in complex-pole hierarchical equations of motion (HEOM). Starting from a phase-space reformulation of the hierarchy, the authors introduce a one-parameter family of linear transformations of the auxiliary ladder operators (Eqs. 7–8) that leaves the bath correlation function C(t) and the exact reduced dynamics invariant while continuously reshaping the HEOM Liouvillian. The resulting ϕ-dependent hierarchy (Eq. 10) recovers conventional complex-pole HEOM at ϕ_k = θ_k. Optimizing a global scaling ϕ_r of the gauge angles yields GO-HEOM, which simultaneously controls the real parts of the Liouvillian spectrum and the Henrici departure from normality. Numerical demonstrations on the spin-boson model show that GO-HEOM eliminates spectral divergences for strongly coupled Brownian-oscillator baths (Figs. 2–3) and suppresses non-normal error growth for sub-Ohmic baths through the delocalized-to-localized quantum phase transition (Fig. 4), extending the accessible coupling range relative to prior complex-pole HEOM.

Significance. If the gauge equivalence and the reported numerical gains hold, the work supplies a practical, decomposition-independent stabilization strategy for complex-pole HEOM that is immediately compatible with existing AAA/Prony decompositions, MPS/TDVP propagators, and filtering techniques. The explicit separation of spectral divergence from non-normal error amplification clarifies two distinct instability mechanisms that have limited HEOM in strong-coupling and long-memory regimes. The phase-space derivation and the alternative generalized-decomposition route in the SM give a transparent algebraic foundation rather than an ad-hoc fix. The demonstrated access to sub-Ohmic dynamics through the quantum phase transition and to previously divergent Brownian-oscillator reorganization energies is of direct interest for electron-transfer and quantum-critical open-system problems.

major comments (2)
  1. The paper deliberately restricts the gauge to the one-parameter family ϕ_k = ϕ_r θ_k (text after Eq. 10 and SM Table I). While the SM correctly notes the larger four-parameter linear transformation (SM Eqs. S35–S43), no numerical evidence is given that the restricted family is near-optimal, nor that the chosen ϕ_r remains stable under changes of hierarchy truncation N_b or MPS bond dimension. A short scan of independent per-mode ϕ_k (or a statement that such a scan yields no further gain) would strengthen the claim that GO-HEOM is a robust, general strategy rather than a convenient but possibly suboptimal slice of the gauge space.
  2. For the sub-Ohmic case the improvement is attributed to reduced non-normality rather than spectral relocation (Fig. 4b and accompanying text). The only quantitative non-normality diagnostic shown is the Henrici measure for the Brownian-oscillator Liouvillian (Fig. 2b). A corresponding plot (or table) of η_H versus ϕ_r for the sub-Ohmic hierarchy at α = 0.2 would make the mechanistic distinction between the two instability channels fully quantitative and would allow readers to verify that the chosen ϕ_r indeed minimizes non-normality.
minor comments (5)
  1. Eq. (10) is dense; a brief parenthetical reminder that the conventional HEOM is recovered at ϕ_k = θ_k (already stated in the text) would help readers navigate the multi-line expression.
  2. Figure 2(d) caption and panel labels use both φ_r and ϕ_r; consistent notation throughout the figures and SM would avoid minor confusion.
  3. SM Table I lists optimized ϕ_r for each α; a one-sentence statement in the main text that these values were obtained by the same short-trajectory / spectrum scan described after Eq. 10 would make the protocol fully self-contained.
  4. The abstract and introduction claim that GO-HEOM is “independent of the bath-correlation decomposition scheme.” A short explicit check (or citation) that the same ϕ_r optimization works for a Prony-fitted decomposition, not only AAA, would make this claim concrete.
  5. Typographical: “Andr´ es” and similar accented characters appear inconsistently in the author line and SM; a uniform encoding would improve readability.

Circularity Check

0 steps flagged

No significant circularity: the gauge family is an exact algebraic equivalence preserving C(t) by construction, and φ_r optimization is independent numerical post-processing that does not redefine the dynamics or observables.

full rationale

The derivation begins from the standard complex-pole HEOM (SM Eq. S6), maps it to an equivalent phase-space/Fokker–Planck form, then applies a continuous linear transformation of the auxiliary bosons (main-text Eqs. 5–9 and the one-parameter family M(φ) of Eq. 8) that leaves the bath correlation function and the exact reduced dynamics invariant. The resulting φ-dependent hierarchy (Eq. 10) recovers the conventional complex-pole HEOM exactly when φ_k= heta_k; an independent derivation via a generalized exponential decomposition of C(t) is also given in the SM. All subsequent claims (control of spectrum and Henrici non-normality, elimination of positive-Re eigenvalues for strong BO coupling, suppression of long-time error growth for sub-Ohmic baths) are demonstrated by direct numerical comparison of the transformed versus untransformed generators (Figs. 2–4 and S1–S3). The single global scaling φ_r is a free numerical hyper-parameter scanned for stability; it is never fitted to physical data nor used to define any observable. No load-bearing uniqueness theorem, ansatz, or self-citation is required for the existence of the gauge or for the reported improvements. The construction is therefore self-contained and non-circular.

Axiom & Free-Parameter Ledger

3 free parameters · 3 axioms · 1 invented entities

The central claim rests on the standard complex-pole HEOM construction, the existence of a linear canonical transformation that preserves the Fokker–Planck operator, and the empirical observation that a one-parameter slice of the full gauge space is already sufficient for practical stabilization. No new physical entities are postulated; the free parameters are purely numerical optimization knobs.

free parameters (3)
  • ϕ_r (global gauge scaling) = 1.1–2.1 (case-dependent)
    Chosen by scanning a one-dimensional interval and selecting the value that minimizes max(Re λ) and/or Henrici non-normality for each spectral density and coupling strength; values range from 1.1 to 2.1 (main text and SM Table I).
  • AAA fitting tolerance ε = 10^{-6}
    Fixed at 10^{-6} to obtain the exponential decomposition of C(t); controls the number of poles (3 for BO, 25 for sub-Ohmic).
  • Hierarchy truncation N_b and MPS bond dimension N_r = N_b=10–25, N_r=200–300
    Increased with coupling strength (N_b = 10–25, N_r up to 300) until dynamics converge; standard numerical cut-offs, not physical parameters.
axioms (3)
  • domain assumption The bath correlation function admits an exponential (complex-pole) decomposition C(t)=∑ d_k e^{-z_k t} that can be obtained by AAA or Prony methods.
    Standard starting point of all complex-pole HEOM formulations (cited as Refs. 36–38, 48, 49).
  • standard math A linear invertible transformation of the auxiliary bosonic operators that leaves the Fokker–Planck operator form-invariant generates an analytically equivalent hierarchy.
    Follows from the phase-space representation already used in earlier HEOM literature (Eqs. 5–9 and SM).
  • domain assumption The reduced system dynamics obtained from any member of the gauge family converge to the same physical trajectory once the hierarchy is sufficiently deep.
    Verified numerically for intermediate coupling (Fig. 2a) and assumed to hold in the strong-coupling regime where conventional HEOM diverges.
invented entities (1)
  • GO-HEOM (gauge-optimized hierarchical equations of motion) independent evidence
    purpose: Name for the practical algorithm that selects the numerically most stable member of the gauge family before propagation.
    A methodological label rather than a new physical object; independent evidence is the set of stable trajectories shown in Figs. 3–4.

pith-pipeline@v1.1.0-grok45 · 23198 in / 3091 out tokens · 25496 ms · 2026-07-11T12:47:02.356264+00:00 · methodology

0 comments
read the original abstract

While complex-pole hierarchical equations of motion (HEOM) have dramatically expanded the reach of numerically exact quantum dynamics simulations of open quantum systems, they suffer from numerical instabilities rooted in the non-Hermitian structure of their Liouvillian. Yet, the origin of this structure remains obscure. Here, we report a previously unknown gauge freedom in complex-pole HEOM: a continuous family of analytically equivalent Liouvillians, all encoding the same bath correlation function, whose numerical properties vary dramatically. This gauge controls both the eigenspectrum and non-normality of the hierarchy generator, revealing spectral divergence and non-normal error amplification as two distinct instability mechanisms. By optimizing this gauge, we introduce GO--HEOM, which eliminates divergences in strongly coupled Brownian oscillator environments and extends numerically exact simulations of sub-Ohmic dynamics -- including through the delocalized-to-localized quantum phase transition -- to previously inaccessible coupling strengths. Because this gauge transformation is independent of the bath-correlation decomposition scheme, our GO--HEOM becomes a general, broadly compatible strategy for accessing numerically exact quantum dynamics of open quantum systems over arbitrary coupling and highly non-Markovian regimes.

Figures

Figures reproduced from arXiv: 2607.04834 by Andr\'es Montoya-Castillo, Tianchu Li.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic derivation of the GO–HEOM and the two [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (a) [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Spin-boson dynamics with a sub-Ohmic spectral den [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗

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

83 extracted references · 2 linked inside Pith

  1. [1]

    Tanaka and Y

    M. Tanaka and Y. Tanimura, Quantum dissipative dy- namics of electron transfer reaction system: Nonpertur- bative hierarchy equation approach, J. Phys. Soc. Jpn. 78, 073802 (2009)

    2009

  2. [2]

    Tanaka and Y

    M. Tanaka and Y. Tanimura, Multistate electron trans- fer dynamics in the condensed phase: Exact calcula- tions from the reduced hierarchy equations of motion ap- proach, J. Chem. Phys.132, 214502 (2010)

    2010

  3. [3]

    Y. Liu, Y. Yan, T. Xing, and Q. Shi, Understanding the large kinetic isotope effect of hydrogen tunneling in condensed phases by using double-well model systems, J. Phys. Chem. B125, 5959 (2021)

    2021

  4. [4]

    Huang, S

    C. Huang, S. Bai, and Q. Shi, Simulation of the pump–probe spectra and excitation energy relaxation of the b850 band of the lh2 complex in purple bacteria, J. Phys. Chem. B128, 7467 (2024)

    2024

  5. [5]

    T. Li, P. Venkatesh, Q. Shi, and A. Montoya-Castillo, For molecular polaritons, disorder and phonon timescales control the activation of dark states in the thermody- namic limit, arXiv preprint arXiv:2603.06868 (2026)

    arXiv 2026

  6. [6]

    Strathearn, P

    A. Strathearn, P. Kirton, D. Kilda, and B. W. Lovett, Efficient non-markovian quantum dynamics using time- evolving matrix product operators, Nat. Commun.9, 3322 (2018)

    2018

  7. [7]

    Leng, Y.-M

    X. Leng, Y.-M. Yan, R.-D. Zhu, K. Song, Y.-X. Weng, and Q. Shi, Simulation of the two-dimensional elec- tronic spectroscopy and energy transfer dynamics of light-harvesting complex ii at ambient temperature, J. Phys. Chem. B122, 4642 (2018)

    2018

  8. [8]

    H¨ artle, G

    R. H¨ artle, G. Cohen, D. R. Reichman, and A. J. Mil- lis, Decoherence and lead-induced interdot coupling in nonequilibrium electron transport through interacting quantum dots: A hierarchical quantum master equation approach, Phys. Rev. B88, 235426 (2013)

    2013

  9. [9]

    Song and Q

    L.-Z. Song and Q. Shi, A new approach to calculate charge carrier transport mobility in organic molecular crystals from imaginary time path integral simulations, J. Chem. Phys.142, 174103 (2015)

    2015

  10. [10]

    Y.-M. Yan, M. Xu, Y.-Y. Liu, and Q. Shi, Theoreti- cal study of charge carrier transport in organic molecu- lar crystals using the nakajima-zwanzig-mori generalized master equation, J. Chem. Phys.150, 234101 (2019)

    2019

  11. [11]

    T. Xing, T. Li, Y. Yan, S. Bai, and Q. Shi, Applica- tion of the imaginary time hierarchical equations of mo- tion method to calculate real time correlation functions, J. Chem. Phys.156, 244102 (2022)

    2022

  12. [12]

    Jankovi´ c and N

    V. Jankovi´ c and N. Vukmirovi´ c, Spectral and thermo- dynamic properties of the holstein polaron: Hierarchical equations of motion approach, Phys. Rev. B105, 054311 (2022)

    2022

  13. [13]

    Kloss, D

    B. Kloss, D. R. Reichman, and R. Tempelaar, Multiset matrix product state calculations reveal mobile franck- condon excitations under strong holstein-type coupling, Phys. Rev. Lett.123, 126601 (2019)

    2019

  14. [14]

    T. Li, Y. Yan, and Q. Shi, Is there a finite mobility for the one vibrational mode Holstein model? Implications from real time simulations, J. Chem. Phys.160, 111102 (2024)

    2024

  15. [15]

    Jankovi´ c, Charge transport limited by nonlocal electron-phonon interaction

    V. Jankovi´ c, Charge transport limited by nonlocal electron-phonon interaction. i. hierarchical equations of motion approach, Phys. Rev. B112, 035111 (2025)

    2025

  16. [16]

    Bhattacharyya, T

    S. Bhattacharyya, T. Sayer, and A. Montoya-Castillo, Nonequilibrium relaxation exponentially delays the onset of quantum diffusion, Proc. Natl. Acad. Sci. USA122, e2424582122 (2025)

    2025

  17. [17]

    W. Li, J. Ren, and Z. Shuai, A general charge trans- port picture for organic semiconductors with nonlo- cal electron-phonon couplings, Nat. Commun.12, 4260 (2021)

    2021

  18. [18]

    Zhang, Y

    S. Zhang, Y. Chen, and Q. Shi, Simulating the opera- tion of a quantum computer in a dissipative environment, J. Chem. Phys.160, 054101 (2024)

    2024

  19. [19]

    Nakamura and J

    K. Nakamura and J. Ankerhold, Impact of time-retarded noise on dynamical decoupling schemes for qubits, Phys. Rev. B111, 064503 (2025)

    2025

  20. [20]

    Nakamura and J

    K. Nakamura and J. Ankerhold, Entanglement dynamics and performance of two-qubit gates for superconducting qubits under non-markovian effects, Phys. Rev. Res.8, 013337 (2026)

    2026

  21. [21]

    Y. Chen, S. Zhang, and Q. Shi, Simulation of one and two qubit superconducting quantum gates in the presence of non-markovian 1/f noise, Phys. Rev. B113, 094302 (2026)

    2026

  22. [22]

    T. Li, P. Venkatesh, N. Shitara, and A. Montoya-Castillo, Numerically exact quantum dynamics with tensor net- works: Predicting the decoherence of interacting spin systems, J. Chem. Phys.164, 091103 (2026)

    2026

  23. [24]

    Tanimura, Stochastic liouville, langevin, fokker–planck, and master equation approaches to quantum dissipative systems, J

    Y. Tanimura, Stochastic liouville, langevin, fokker–planck, and master equation approaches to quantum dissipative systems, J. Phys. Soc. Jpn.75, 082001 (2006)

    2006

  24. [25]

    Y.-A. Yan, F. Yang, Y. Liu, and J.-S. Shao, Hierarchical approach based on stochastic decoupling to dissipative systems, Chem. Phys. Lett.395, 216 (2004)

    2004

  25. [26]

    Tanimura, Numerically “exact” approach to open quantum dynamics: The hierarchical equations of mo- tion (HEOM), J

    Y. Tanimura, Numerically “exact” approach to open quantum dynamics: The hierarchical equations of mo- tion (HEOM), J. Chem. Phys.153, 020901 (2020)

    2020

  26. [27]

    S. Bai, S. Zhang, C. Huang, and Q. Shi, Hierarchi- cal equations of motion for quantum chemical dynam- ics: Recent methodology developments and applications, Acc. Chem. Res.57, 3151 (2024)

    2024

  27. [28]

    Q. Shi, L. Chen, G. Nan, R.-X. Xu, and Y. Yan, Efficient hierarchical liouville space propagator to quantum dissi- pative dynamics, J. Chem. Phys.130, 084105 (2009)

    2009

  28. [29]

    Q. Shi, Y. Xu, Y. Yan, and M. Xu, Efficient propaga- tion of the hierarchical equations of motion using the ma- trix product state method, J. Chem. Phys.148, 174102 (2018)

    2018

  29. [30]

    Borrelli and S

    R. Borrelli and S. Dolgov, Expanding the range of hier- archical equations of motion by tensor-train implemen- tation, J. Phys. Chem. B125, 5397 (2021)

    2021

  30. [31]

    Y. Yan, Y. Liu, T. Xing, and Q. Shi, Theoretical study of excitation energy transfer and nonlinear spectroscopy of photosynthetic light-harvesting complexes using the nonperturbative reduced dynamics method, Wiley Inter- discip. Rev.: Comput. Mol. Sci.11, e1498 (2021)

    2021

  31. [32]

    Y. Ke, R. Borrelli, and M. Thoss, Hierarchical equations of motion approach to hybrid fermionic and bosonic en- vironments: Matrix product state formulation in twin space, J. Chem. Phys.156(2022). 7

    2022

  32. [33]

    Ke, Tree tensor network state approach for solving hierarchical equations of motion, J

    Y. Ke, Tree tensor network state approach for solving hierarchical equations of motion, J. Chem. Phys.158 (2023)

    2023

  33. [34]

    Chen and I

    X. Chen and I. Franco, Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in struc- tured thermal environments, J. Chem. Phys.163(2025)

    2025

  34. [35]

    W. Guan, P. Bao, J. Peng, Z. Lan, and Q. Shi, mp- sqd: A matrix product state based python package to simulate closed and open system quantum dynamics, J. Chem. Phys.161, 122501 (2024)

    2024

  35. [36]

    Z.-H. Chen, Y. Wang, X. Zheng, R.-X. Xu, and Y. Yan, Universal time-domain Prony fitting decompo- sition for optimized hierarchical quantum master equa- tions, J. Chem. Phys.156, 221102 (2022)

    2022

  36. [37]

    Zhang, A

    L. Zhang, A. Erpenbeck, Y. Yu, and E. Gull, Minimal pole representation for spectral functions, J. Chem. Phys. 162(2025)

    2025

  37. [38]

    M. Xu, Y. Yan, Q. Shi, J. Ankerhold, and J. Stockburger, Taming quantum noise for efficient low temperature sim- ulations of open quantum systems, Phys. Rev. Lett.129, 230601 (2022)

    2022

  38. [39]

    A. C. Hunt and S. C. Althorpe, Exploiting the path- integral radius of gyration in open quantum dynamics, J. Chem. Phys.164(2026)

    2026

  39. [40]

    T. Li, Y. Yan, and Q. Shi, A low-temperature quan- tum fokker–planck equation that improves the numerical stability of the hierarchical equations of motion for the brownian oscillator spectral density, J. Chem. Phys.156, 064107 (2022)

    2022

  40. [41]

    Y. Yan, M. Xu, T. Li, and Q. Shi, Efficient propagation of the hierarchical equations of motion using the tucker and hierarchical tucker tensors, J. Chem. Phys.154, 194104 (2021)

    2021

  41. [42]

    I. S. Dunn, R. Tempelaar, and D. R. Reichman, Re- moving instabilities in the hierarchical equations of mo- tion: Exact and approximate projection approaches, J. Chem. Phys.150, 184109 (2019)

    2019

  42. [43]

    Firmino, E

    T. Firmino, E. Mangaud, F. Cailliez, A. Devolder, D. Mendive-Tapia, F. Gatti, C. Meier, M. Desouter- Lecomte, and A. De La Lande, Quantum effects in ul- trafast electron transfers within cryptochromes, Phys. Chem. Chem. Phys.18, 21442 (2016)

    2016

  43. [44]

    Winter, H

    A. Winter, H. Rieger, M. Vojta, and R. Bulla, Quan- tum phase transition in the sub-ohmic spin-boson model: Quantum monte carlo study with a continuous imaginary time cluster algorithm, Phys. Rev. Lett.102, 030601 (2009)

    2009

  44. [45]

    Goulko, H.-T

    O. Goulko, H.-T. Chen, M. Goldstein, and G. Cohen, Transient dynamical phase diagram of the spin-boson model, Phys. Rev. Lett.134, 056502 (2025)

    2025

  45. [46]

    Wang and M

    H. Wang and M. Thoss, From coherent motion to local- ization: Ii. dynamics of the spin-boson model with sub- ohmic spectral density at zero temperature, Chem. Phys. 370, 78 (2010)

    2010

  46. [47]

    Bulla, N

    R. Bulla, N. Tong, and M. Vojta, Numerical renormal- ization group for bosonic systems and application to the sub-ohmic spin-boson model, Phys. Rev. Lett.91, 170601 (2003)

    2003

  47. [48]

    Nakatsukasa, O

    Y. Nakatsukasa, O. S` ete, and L. N. Trefethen, The aaa algorithm for rational approximation, SIAM J. Sci. Com- put.40, A1494 (2018)

    2018

  48. [49]

    X. Dan, M. Xu, J. 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, Phys. Rev. B107, 195429 (2023)

    2023

  49. [50]

    This is structurally analogous to the Lindblad jump-operator gauge [?], and distinct from the trivial rescaling of ADOs discussed in Ref

    We use the term gauge in the sense standard in physics: a transformation acting on unphysical (here, auxiliary) degrees of freedom that leaves all physical observables invariant while reshaping the mathematical representa- tion of the generator. This is structurally analogous to the Lindblad jump-operator gauge [?], and distinct from the trivial rescaling...

  50. [51]

    A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. Fisher, A. Garg, and W. Zwerger, Dynamics of the dissi- pative two-level system, Rev. Mod. Phys.59, 1 (1987)

    1987

  51. [52]

    May and O

    V. May and O. K¨ uhn,Charge and Energy Transfer Dynamics in Molecular Systems, 2nd ed. (Wiley-VCH, Weinheim, 2004)

    2004

  52. [53]

    H. Liu, L. Zhu, S. Bai, and Q. Shi, Reduced Quantum Dynamics with Arbitrary Bath Spectral Densities: Hi- erarchical Equations of Motion Based on Several Differ- ent Bath Decomposition Schemes., J. Chem. Phys.140, 134106 (2014)

    2014

  53. [54]

    Shi, L.-P

    Q. Shi, L.-P. Chen, G.-J. Nan, R.-X. Xu, and Y.-J. Yan, Electron transfer dynamics: Zusman equation versus ex- act theory, J. Chem. Phys.130, 164518 (2009)

    2009

  54. [55]

    Y. Yan, T. Xing, and Q. Shi, A new method to im- prove the numerical stability of the hierarchical equa- tions of motion for discrete harmonic oscillator modes, J. Chem. Phys.153, 214109 (2020)

    2020

  55. [56]

    T. Li, C. Huang, S. Bai, and Q. Shi, Theoretical meth- ods based on linear response theory to simulate dy- namics and absorption spectra of molecular polaritons, J. Chem. Phys.162(2025)

    2025

  56. [58]

    A. Garg, J. N. Onuchic, and V. Ambegaokar, Ef- fect of friction on electron transfer in biomolecules, J. Chem. Phys.83, 4491 (1985)

    1985

  57. [59]

    Ikeda and G

    T. Ikeda and G. D. Scholes, Generalization of the hierar- chical equations of motion theory for efficient calculations with arbitrary correlation functions, J. Chem. Phys.152, 204101 (2020)

    2020

  58. [60]

    Ito and Y

    H. Ito and Y. Tanimura, Simulating two-dimensional infrared-raman and raman spectroscopies for inter- molecular and intramolecular modes of liquid water, J. Chem. Phys.144, 074201 (2016)

    2016

  59. [61]

    Meier and D

    C. Meier and D. Tannor, Non-markovian evolution of the density operator in presence of strong laser fields, J. Chem. Phys.111, 3365 (1999)

    1999

  60. [62]

    L. N. Trefethen, A. E. Trefethen, S. C. Reddy, and T. A. Driscoll, Hydrodynamic stability without eigenval- ues, Science261, 578 (1993)

    1993

  61. [63]

    P. J. Schmid, Nonmodal stability theory, Annu. Rev. Fluid Mech.39, 129 (2007)

    2007

  62. [64]

    C. Duan, Z. Tang, J. Cao, and J. Wu, Zero-temperature localization in a sub-ohmic spin-boson model investi- gated by an extended hierarchy equation of motion, Phys. Rev. B95, 214308 (2017)

    2017

  63. [65]

    C. Gong, Z. Sheng, and W. Li, Entanglement structure of the dynamical phases in the sub-ohmic spin-boson model, arXiv preprint arXiv:2606.20313 (2026)

    Pith/arXiv arXiv 2026

  64. [66]

    Hartmann and W

    R. Hartmann and W. T. Strunz, Exact open quan- tum system dynamics using the hierarchy of pure states 8 (hops), J. Chem. Theory Comput.13, 5834 (2017)

    2017

  65. [67]

    D. E. Makarov and N. Makri, Path integrals for dis- sipative systems by tensor multiplication. condensed phase quantum dynamics for arbitrarily long time, Chem. Phys. Lett.221, 482 (1994)

    1994

  66. [68]

    D. E. Makarov and N. Makri, Stochastic resonances and nonlinear response in double-quantum-well structures, Phys. Rev. B52, R2257 (1995)

    1995

  67. [69]

    G. E. Fux, D. Kilda, B. W. Lovett, and J. Keeling, Ten- sor network simulation of chains of non-markovian open quantum systems, Phys. Rev. Res.5, 033078 (2023)

    2023

  68. [70]

    Suess, A

    D. Suess, A. Eisfeld, and W. T. Strunz, Hierarchy of stochastic pure states for open quantum system dynam- ics, Phys. Rev. Lett.113, 150403 (2014)

    2014

  69. [71]

    Hidden Gauge Freedom in Complex-Pole Hierarchical Equations of Motion

    R. Hartmann and W. T. Strunz, Exact open quan- tum system dynamics using the hierarchy of pure states (hops), J. Chem. Theory Comput.13, 5834 (2017). Supplemental Material for “Hidden Gauge Freedom in Complex-Pole Hierarchical Equations of Motion” Tianchu Li 1 and Andr´ es Montoya-Castillo1,∗ 1Department of Chemistry, University of Colorado Boulder, Bould...

    2017

  70. [72]

    Q. Shi, Y. Xu, Y. Yan, and M. Xu, Efficient propagation of the hierarchical equations of motion using the matrix product state method, J. Chem. Phys.148, 174102 (2018)

    2018

  71. [73]

    Y. Ke, R. Borrelli, and M. Thoss, Hierarchical equations of motion approach to hybrid fermionic and bosonic environments: Matrix product state formulation in twin space, J. Chem. Phys.156(2022)

    2022

  72. [74]

    Borrelli and S

    R. Borrelli and S. Dolgov, Expanding the range of hierarchical equations of motion by tensor- train implementation, J. Phys. Chem. B125, 5397 (2021)

    2021

  73. [75]

    Lubich, I

    C. Lubich, I. Oseledets, and B. Vandereycken, Time integration of tensor trains, SIAM J. Numer. Anal.53, 917 (2015)

    2015

  74. [76]

    M. Xu, Y. Yan, Q. Shi, J. Ankerhold, and J. Stockburger, Taming quantum noise for efficient low temperature simulations of open quantum systems, Phys. Rev. Lett.129, 230601 (2022)

    2022

  75. [77]

    Shi, L.-P

    Q. Shi, L.-P. Chen, G.-J. Nan, R.-X. Xu, and Y.-J. Yan, Electron transfer dynamics: Zusman equation versus exact theory, J. Chem. Phys.130, 164518 (2009)

    2009

  76. [78]

    H. Liu, L. Zhu, S. Bai, and Q. Shi, Reduced Quantum Dynamics with Arbitrary Bath Spectral Densities: Hierarchical Equations of Motion Based on Several Different Bath Decomposition Schemes., J. Chem. Phys.140, 134106 (2014)

    2014

  77. [79]

    Y. Yan, T. Xing, and Q. Shi, A new method to improve the numerical stability of the hierar- chical equations of motion for discrete harmonic oscillator modes, J. Chem. Phys.153, 214109 (2020)

    2020

  78. [80]

    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, J. Chem. Phys.156, 064107 (2022)

    2022

  79. [81]

    A. Garg, J. N. Onuchic, and V. Ambegaokar, Effect of friction on electron transfer in biomolecules, J. Chem. Phys.83, 4491 (1985)

    1985

  80. [82]

    Tanimura and R

    Y. Tanimura and R. Kubo, Time evolution of a quantum system in contact with a nearly gaussian-markoffian noise bath, J. Phys. Soc. Jpn.58, 101 (1989)

    1989

Showing first 80 references.