One-to-one correspondence between Hierarchical Equations of Motion and Pseudomodes for Open Quantum System Dynamics
Pith reviewed 2026-05-10 18:35 UTC · model grok-4.3
The pith
Every exponential bath correlation function arises from a Lindblad-damped pseudomode model and maps exactly onto the HEOM hierarchy via a linear transformation.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We prove that every physical bath correlation function that can be written as a sum of N exponential terms can be obtained from a physical model with N interacting pseudomodes which are damped in Lindblad form. For every such bath correlation function there exists a non-unitary, linear transformation which mirrors the evolution of the system-pseudomode state onto the HEOM hierarchy, and vice versa. The proofs are constructive and we give explicit expressions for the mirror transformation as well as for the pseudomode Lindbladian corresponding to a given exponential bath correlation function.
What carries the argument
The non-unitary linear mirror transformation that equates the density-matrix evolution of the system plus Lindblad-damped pseudomodes with the HEOM hierarchy, constructed directly from the exponential coefficients of the bath correlation function.
Load-bearing premise
Any bath correlation function of physical interest can be represented exactly as a finite sum of exponential terms.
What would settle it
A concrete calculation for an exponential bath correlation function in which the transformed pseudomode state fails to satisfy the corresponding HEOM equations at any time after the initial condition.
Figures
read the original abstract
We unite two of the most widely used approaches for strongly damped, non-Markovian open quantum dynamics, the Hierarchical Equations of Motion (HEOM) and the pseudomode method by proving two statements: First, every physical bath correlation function (BCF) that can be written as a sum of $N$ exponential terms can be obtained from a physical model with $N$ interacting pseudomodes which are damped in Lindblad form. Second, for every such BCF there exists a non-unitary, linear transformation which mirrors the evolution of the system-pseudomode state onto the HEOM hierarchy, and vice versa. Our proofs are constructive and we give explicit expressions for the mirror transformation as well as for the pseudomode Lindbladian corresponding to a given exponential BCF. This approach also gives insight and provides elegant derivations of the corresponding Hierarchy of stochastic Pure States (HOPS) method and its nearly-unitary version, nuHOPS. Our result opens several avenues for further optimization of non-Markovian open quantum system dynamics methods.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proves a one-to-one correspondence between the Hierarchical Equations of Motion (HEOM) and the pseudomode method for open quantum system dynamics. For any physical bath correlation function (BCF) expressible as a finite sum of N exponential terms, it constructs an equivalent physical model of N interacting pseudomodes damped via Lindblad operators. It further supplies an explicit non-unitary linear transformation that maps the evolution of the system-pseudomode state onto the HEOM hierarchy (and vice versa). The proofs are constructive, providing closed-form expressions for the Lindbladian and the mirror map; the work also yields derivations of the Hierarchy of stochastic Pure States (HOPS) and its nearly-unitary variant (nuHOPS).
Significance. If the claimed equivalence holds, the result unifies two standard numerical frameworks for strongly non-Markovian open-system dynamics and supplies a rigorous route to hybrid or optimized implementations. The constructive, parameter-free character of the proofs (derived from standard open-system axioms without fitted quantities or ad-hoc entities) is a clear strength, as is the explicit invertibility of the linear map. This should facilitate cross-validation between HEOM and pseudomode codes and clarify the status of stochastic unravelings such as HOPS.
minor comments (2)
- [Abstract / §3] The abstract states that the pseudomode model remains physical for all parameter regimes, but the main text should include a short explicit check (e.g., after Eq. (X) defining the Lindblad operators) that the resulting rates and interaction matrices preserve complete positivity for arbitrary exponential coefficients.
- [§4] Notation for the auxiliary density operators in the HEOM hierarchy and the pseudomode state vector should be aligned more clearly (perhaps in a dedicated comparison table) to make the linear map immediately usable by readers implementing the correspondence.
Simulated Author's Rebuttal
We thank the referee for their positive assessment of our manuscript, accurate summary of the main results, and recommendation for minor revision. The report correctly identifies the constructive nature of the proofs and the potential for unifying HEOM and pseudomode approaches. Below we respond to the referee summary as the primary point raised.
read point-by-point responses
-
Referee: The manuscript proves a one-to-one correspondence between the Hierarchical Equations of Motion (HEOM) and the pseudomode method for open quantum system dynamics. For any physical bath correlation function (BCF) expressible as a finite sum of N exponential terms, it constructs an equivalent physical model of N interacting pseudomodes damped via Lindblad operators. It further supplies an explicit non-unitary linear transformation that maps the evolution of the system-pseudomode state onto the HEOM hierarchy (and vice versa). The proofs are constructive, providing closed-form expressions for the Lindbladian and the mirror map; the work also yields derivations of the Hierarchy of stochastic Pure States (HOPS) and its nearly-unitary variant (nuHOPS).
Authors: We appreciate the referee's concise and accurate encapsulation of our contributions. The one-to-one correspondence is indeed established via the explicit construction of the Lindbladian for the interacting pseudomodes and the invertible non-unitary linear mirror map, both derived directly from the exponential form of the BCF without additional assumptions. The derivations of HOPS and nuHOPS follow naturally from the same framework as stochastic unravelings of the pseudomode dynamics. revision: no
Circularity Check
No significant circularity in the equivalence proof
full rationale
The paper presents constructive proofs deriving explicit Lindblad operators for interacting pseudomodes and a non-unitary linear map between the pseudomode state and HEOM auxiliary operators, directly from the assumption that the bath correlation function is a finite sum of exponentials and standard open-system Lindblad dynamics. No derivations reduce to self-definitions, fitted inputs renamed as predictions, or load-bearing self-citations; the equivalence follows from first-principles mappings without circular reduction. The result is self-contained against external benchmarks such as standard quantum master equations.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Bath correlation functions of physical interest can be represented as finite sums of exponential terms
- domain assumption Pseudomodes damped by Lindblad operators constitute a valid physical model for the environment
Forward citations
Cited by 1 Pith paper
-
Kramers-Kronig Relations and Causality in Non-Markovian Open Quantum Dynamics: Kernel, State, and Effective Kernel
The Nakajima-Zwanzig memory kernel belongs to the operator-valued Hardy space and obeys Kramers-Kronig relations under a real-axis spectral hypothesis, while effective kernels can show upper-half-plane poles from unca...
Reference graph
Works this paper leans on
-
[1]
like uniTEMPO (uniform time evolving matrix product operator) [14] can be understood in this manner. Here we focus specifically on the relation between the pseudomode approach, and the hierarchical methods (HEOM, HOPS, nuHOPS), and we touch chain mapping approaches later. While they all embed the system into an effective environment, the construction and ...
work page internal anchor Pith review Pith/arXiv arXiv 2026
- [2]
-
[3]
I. de Vega and D. Alonso, Dynamics of non-Markovian open quantum systems, Rev. Mod. Phys.89, 015001 (2017)
work page 2017
-
[4]
M. Xu, V . Vadimov, J. T. Stockburger, and J. Ankerhold, Collo- quium: Simulating non-Markovian dynamics in open quantum systems, Rev. Mod. Phys. 10.1103/w3nw-hbjc (2026)
-
[5]
U. Weiss,Quantum Dissipative Systems, Series In Modern Con- densed Matter Physics (World Scientific Publishing Company, Singapore, 2011)
work page 2011
-
[6]
S. Gröblacher, A. Trubarov, N. Prigge, G. D. Cole, M. As- pelmeyer, and J. Eisert, Observation of non-Markovian mi- cromechanical Brownian motion, Nat. Commun.6, 7606 (2015)
work page 2015
-
[7]
X. Dai, R. Trappen, H. Chen, D. Melanson, M. A. Yurtalan, D. M. Tennant, A. J. Martinez, Y . Tang, E. Mozgunov, J. Gib- son, J. A. Grover, S. M. Disseler, J. I. Basham, S. Novikov, R. Das, A. J. Melville, B. M. Niedzielski, C. F. Hirjibehedin, K. Serniak, S. J. Weber, J. L. Yoder, W. D. Oliver, K. M. Zick, D. A. Lidar, and A. Lupascu, Dissipative Landau-Z...
work page 2025
-
[8]
A. Chin, J. Keeling, D. Segal, and H. Wang, Algorithms and software for open quantum system dynamics, J. Chem. Phys. 163, 10.1063/5.0289390 (2025)
-
[9]
J. Cerrillo and J. Cao, Non-markovian dynamical maps: Nu- merical processing of open quantum trajectories, Phys. Rev. Lett.112, 110401 (2014)
work page 2014
- [10]
-
[11]
K. Müller and W. T. Strunz, Quantum trajectory method for highly excited environments in non-Markovian open quantum dynamics, Phys. Rev. A112, 033719 (2025)
work page 2025
-
[12]
V . Link, K. Luoma, and W. T. Strunz, Non-markovian quantum state diffusion for spin environments, New Journal of Physics 25, 093006 (2023)
work page 2023
-
[13]
A. Strathearn, P. Kirton, D. Kilda, J. Keeling, and B. W. Lovett, Efficient non-markovian quantum dynamics using time- evolving matrix product operators, Nat. Commun.9, 1 (2018)
work page 2018
-
[14]
P. Fowler-Wright, B. W. Lovett, and J. Keeling, Efficient many- body non-markovian dynamics of organic polaritons, Phys. Rev. Lett.129, 173001 (2022)
work page 2022
-
[15]
V . Link, H.-H. Tu, and W. T. Strunz, Open quantum system dynamics from infinite tensor network contraction, Phys. Rev. Lett.132, 200403 (2024)
work page 2024
-
[16]
M. Cygorek, M. Cosacchi, A. Vagov, V . M. Axt, B. W. Lovett, J. Keeling, and E. M. Gauger, Simulation of open quantum sys- tems by automated compression of arbitrary environments, Nat. Phys.18, 662 (2022)
work page 2022
- [17]
-
[18]
T. Lacroix, B. Le Dé, A. Riva, A. J. Dunnett, and A. W. Chin, MPSDynamics.jl: Tensor network simulations for finite- temperature (non-Markovian) open quantum system dynamics, J. Chem. Phys.161, 084116 (2024)
work page 2024
-
[19]
B. M. Garraway, Nonperturbative decay of an atomic system in a cavity, Phys. Rev. A55, 2290 (1997). 6
work page 1997
-
[20]
G. Pleasance, B. M. Garraway, and F. Petruccione, Generalized theory of pseudomodes for exact descriptions of non-markovian quantum processes, Phys. Rev. Res.2, 043058 (2020)
work page 2020
-
[21]
F. Mascherpa, A. Smirne, A. D. Somoza, P. Fernández-Acebal, S. Donadi, D. Tamascelli, S. F. Huelga, and M. B. Plenio, Op- timized auxiliary oscillators for the simulation of general open quantum systems, Phys. Rev. A101, 052108 (2020)
work page 2020
-
[22]
Z. Huang, G. Park, G. K.-L. Chan, and L. Lin, Coupled lind- blad pseudomode theory for simulating open quantum systems, arXiv 10.48550/arXiv.2506.10308 (2025), 2506.10308
- [23]
-
[24]
H. Wang and M. Thoss, Multilayer formulation of the multicon- figuration time-dependent Hartree theory, J. Chem. Phys.119, 1289 (2003)
work page 2003
-
[25]
H. Carmichael,Statistical Methods in Quantum Optics 1: Mas- ter Equations and Fokker-Planck Equations, Physics and as- tronomy online library (Springer, 1999)
work page 1999
-
[26]
Y . Zhao, The hierarchy of Davydov’s ansätze: From guesswork to numerically “exact” many-body wave functions, J. Chem. Phys.158, 10.1063/5.0140002 (2023)
-
[27]
M. Werther and F. Großmann, Apoptosis of moving nonorthog- onal basis functions in many-particle quantum dynamics, Phys. Rev. B101, 174315 (2020)
work page 2020
-
[28]
Y . Tanimura, Nonperturbative expansion method for a quantum system coupled to a harmonic-oscillator bath, Phys. Rev. A41, 6676 (1990)
work page 1990
-
[29]
Y . Tanimura, Numerically “exact” approach to open quantum dynamics: The hierarchical equations of motion (HEOM), The Journal of Chemical Physics153, 020901 (2020)
work page 2020
-
[30]
R. P. Feynman and F. L. Vernon, The theory of a general quan- tum system interacting with a linear dissipative system, Ann. Phys.24, 118 (1963)
work page 1963
-
[31]
Imamog¯lu, Stochastic wave-function approach to non- markovian systems, Phys
A. Imamog¯lu, Stochastic wave-function approach to non- markovian systems, Phys. Rev. A50, 3650 (1994)
work page 1994
-
[32]
L. K. Zhou, G. R. Jin, and W. Yang, Systematic and efficient pseudomode method to simulate open quantum systems under a bosonic environment, Phys. Rev. A110, 022221 (2024)
work page 2024
- [33]
-
[34]
D. Lentrodt and J. Evers, Ab initio few-mode theory for quan- tum potential scattering problems, Phys. Rev. X10, 011008 (2020)
work page 2020
-
[35]
R. Bulla, H.-J. Lee, N.-H. Tong, and M. V ojta, Numerical renor- malization group for quantum impurities in a bosonic bath, Phys. Rev. B71, 045122 (2005)
work page 2005
-
[36]
K. H. Hughes, C. D. Christ, and I. Burghardt, Effective-mode representation of non-Markovian dynamics: A hierarchical ap- proximation of the spectral density. I. Application to single sur- face dynamics, J. Chem. Phys.131, 024109 (2009)
work page 2009
-
[37]
A. W. Chin, Á. Rivas, S. F. Huelga, and M. B. Plenio, Exact mapping between system-reservoir quantum models and semi- infinite discrete chains using orthogonal polynomials, J. Math. Phys.51, 092109 (2010)
work page 2010
-
[38]
M. Sánchez-Barquilla and J. Feist, Accurate truncations of chain mapping models for open quantum systems, Nanomateri- als11, 2104 (2021)
work page 2021
-
[39]
C.-Y . Hsieh and J. Cao, A unified stochastic formulation of dis- sipative quantum dynamics. I. Generalized hierarchical equa- tions, J. Chem. Phys.148, 10.1063/1.5018725 (2018)
-
[40]
J. Jin, X. Zheng, and Y . Yan, Exact dynamics of dissipative elec- tronic systems and quantum transport: Hierarchical equations of motion approach, J. Chem. Phys.128, 234703 (2008)
work page 2008
- [41]
-
[42]
B. Debecker, J. Martin, and F. Damanet, Spectral theory of non-Markovian dissipative phase transitions, Phys. Rev. A110, 042201 (2024)
work page 2024
-
[43]
M. Xu, Y . Yan, Q. Shi, J. Ankerhold, and J. T. Stockburger, Taming quantum noise for efficient low temperature simula- tions of open quantum systems, Phys. Rev. Lett.129, 230601 (2022)
work page 2022
-
[44]
R. Hartmann and W. T. Strunz, Exact open quantum system dy- namics using the hierarchy of pure states (HOPS), Journal of Chemical Theory and Computation13, 5834 (2017)
work page 2017
-
[45]
B. Citty, J. K. Lynd, T. Gera, L. Varvelo, and D. I. G. B. Raccah, MesoHOPS: Size-invariant scaling calculations of multi-excitation open quantum systems, J. Chem. Phys.160, 10.1063/5.0197825 (2024)
-
[46]
H. Takahashi, S. Rudge, C. Kaspar, M. Thoss, and R. Borrelli, High accuracy exponential decomposition of bath correlation functions for arbitrary and structured spectral densities: Emerg- ing methodologies and new approaches, J. Chem. Phys.160, 204105 (2024)
work page 2024
-
[47]
I. S. Dunn, R. Tempelaar, and D. R. Reichman, Removing instabilities in the hierarchical equations of motion: Exact and approximate projection approaches, J. Chem. Phys.150, 10.1063/1.5092616 (2019)
-
[48]
Y . Yan, T. Xing, and Q. Shi, A new method to improve the nu- merical stability of the hierarchical equations of motion for dis- crete harmonic oscillator modes, J. Chem. Phys.153, 204109 (2020)
work page 2020
-
[49]
M. Krug and J. Stockburger, On stability issues of the heom method, Eur. Phys. J. Spec. Top.232, 3219 (2023)
work page 2023
-
[50]
W. Alford, L. P. Bettmann, and G. T. Landi, Subtleties in the pseudomodes formalism, arXiv 10.48550/arXiv.2509.16377 (2025), 2509.16377
-
[51]
J. Thoenniss, I. Vilkoviskiy, and D. A. Abanin, Efficient pseu- domode representation and complexity of quantum impurity models, Phys. Rev. B112, 155114 (2025)
work page 2025
-
[52]
M. Yu, W. T. Strunz, and S. Nimmrichter, Non-Markovian dy- namics of the giant atom beyond the rotating-wave approxima- tion, arXiv 10.48550/arXiv.2601.03383 (2026), 2601.03383
-
[53]
X. Gao, J. Ren, A. Eisfeld, and Z. Shuai, Non-Markovian stochastic Schrödinger equation: Matrix-product-state ap- proach to the hierarchy of pure states, Phys. Rev. A105, L030202 (2022)
work page 2022
-
[54]
S. Flannigan, F. Damanet, and A. J. Daley, Many-body quan- tum state diffusion for non-Markovian dynamics in strongly in- teracting systems, Phys. Rev. Lett.128, 063601 (2022)
work page 2022
-
[55]
K. R. Parthasarathy,An Introduction to Quantum Stochastic Calculus(Springer, Basel, Switzerland, 1992)
work page 1992
-
[56]
C. Gardiner and P. Zoller,Quantum Noise(Springer, Berlin, Germany, 2004)
work page 2004
- [57]
-
[58]
A. H. Sayed and T. Kailath, A survey of spectral factorization methods, Numer. Linear Algebra Appl.8, 467 (2001)
work page 2001
-
[59]
J. F. Mahoney and B. D. Sivazlian, Partial fractions expansion: a review of computational methodology and efficiency, J. Com- put. Appl. Math.9, 247 (1983)
work page 1983
-
[60]
G. Park, Z. Huang, Y . Zhu, C. Yang, G. K.-L. Chan, and L. Lin, Quasi-lindblad pseudomode theory for open quantum systems, Phys. Rev. B110, 195148 (2024)
work page 2024
-
[61]
H. M. Wiseman and G. J. Milburn,Quantum Measurement and Control(Cambridge University Press, Cambridge, Eng- 7 land, UK, 2010)
work page 2010
-
[62]
J. Dalibard, Y . Castin, and K. Mølmer, Wave-function approach to dissipative processes in quantum optics, Phys. Rev. Lett.68, 580 (1992)
work page 1992
-
[63]
Carmichael,An Open Systems Approach to Quantum Optics (Springer, Berlin, Germany, 1993)
H. Carmichael,An Open Systems Approach to Quantum Optics (Springer, Berlin, Germany, 1993)
work page 1993
-
[64]
V . Sukharnikov, S. Chuchurka, and F. Schlawin, Non- Markovian dynamics in nonstationary Gaussian baths: A hi- erarchy of pure states approach, Phys. Rev. Res.8, 013123 (2026)
work page 2026
-
[65]
P. C. Parks, A. M. Lyapunov’s stability theory—100 years on∗, IMA J. Math. Control Inf.9, 275 (1992)
work page 1992
-
[66]
Heuser,Lehrbuch der Analysis 1, 16th ed
H. Heuser,Lehrbuch der Analysis 1, 16th ed. (Teubner, Wies- baden, 2006)
work page 2006
-
[67]
Schollwöck, The density-matrix renormalization group in the age of matrix product states, Ann
U. Schollwöck, The density-matrix renormalization group in the age of matrix product states, Ann. Phys.326, 96 (2011)
work page 2011
- [68]
-
[69]
R. Hartmann,Exact Open Quantum System Dynamics - Inves- tigating Environmentally Induced Entanglement, Ph.D. thesis, Technische Universität Dresden (2021)
work page 2021
-
[70]
V . Boettcher, R. Hartmann, K. Beyer, and W. T. Strunz, Dy- namics of a strongly coupled quantum heat engine—computing bath observables from the hierarchy of pure states, J. Chem. Phys.160, 10.1063/5.0192075 (2024). 8 Supplemental Material In the following we present complete proofs of the results stated in the main text. After discussing the interaction ...
-
[71]
The system and pseudomode dynamics is given by the Lindbladian in Eq
PSEUDOMODE INTERACTION PICTURE As mentioned in the main text the pseudomode method uses an Ansatz-Lindbladian describing damped, coupled harmonic os- cillators to approximate the bath correlation function (BCF) of the original Gaussian environment. The system and pseudomode dynamics is given by the Lindbladian in Eq. (6) as ˙ρ=−i[H pm,ρ] + ∑ k LkρL† k − 1...
-
[72]
,zN(t))T ∈C N for anN-dimensional, complex OU-vector
GENERAL RESULTS ON MULTIV ARIATE ORNSTEIN-UHLENBECK PROCESSES In the following we present an analysis of (c-number) multivariate Ornstein-Uhlenbeck (OU) processesz(t):= (z1(t),z 2(t), . . . ,zN(t))T ∈C N for anN-dimensional, complex OU-vector. An insertion of the identity in terms of coherent states in Eq. (S6) makes it clear that the pseudomode BCF is id...
-
[73]
PROOF OF PSEUDOMODE REPRESENTABILITY We start from a bath correlation function (BCF) that can be expressed as a sum of exponential terms αexp(τ) = N ∑ j=1 G je−λ jτ ,G j,λ j ∈C,forτ≥0,(S30) withλ j =γ j +iω j and positive real partsγ j >0, and we setα exp(τ) =α exp(−τ) ∗ forτ≤0, as required from the fundamental relation Eq. (2). Note that this impliesα ex...
-
[74]
DISSIPON TRANSFORMATION Now assume that we have constructed the pseudomode model – as explained above – corresponding to a physical bath correlation function of the usual exponential form (S30). Thus, the parameters of that pseudomode model Eq. (S43), including V,λ i are known. In the following we explicitly derive the HEOM and HOPS equations from the dis...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.