Reuse & Permissions

It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 4.0 International license. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

Export citation

Export citation

Choose format for download:

Download Citation
  • Open Access

Long-time exponential decompositions of non-Markovian Gaussian baths

Zhen Huang1,2,3, Zhiyan Ding4, Ke Wang4, Jason Kaye2,3, Xiantao Li5, and Lin Lin6,7,8,*

  • *Contact author: linlin@math.berkeley.edu

Phys. Rev. Research 8, 043053 – Published 15 October, 2026

DOI: https://doi.org/10.1103/7765-s8lk

Abstract

Gaussian baths are widely used to model non-Markovian environments, yet the cost of accurate simulation at long times remains poorly understood, especially when spectral densities exhibit nonanalytic behavior as in a range of realistic models. We rigorously bound the complexity of representing bath correlation functions on a time interval [0,T] by sums of complex exponentials, as employed in recent variants of pseudomode and hierarchical equations of motion methods. These bounds make explicit the dependence on the maximal simulation time T, inverse temperature β, and the type and strength of singularities in an effective spectral density. For a broad class of spectral densities, the required number of exponentials is bounded independently of T, achieving time-uniform complexity. The T dependence emerges only as polylogarithmic factors for spectral densities with strong singularities, such as step discontinuities and inverse power-law divergences. The temperature dependence is mild for bosonic environments and disappears entirely for fermionic environments. Thus, the main bottleneck for long-time simulation is the presence of sharp nonanalytic features in the bath spectrum. Our results are instructive both for long-time simulation of non-Markovian open quantum systems and for Markovian embeddings of classical generalized Langevin equations with memory kernels.

View figure in article

Physics Subject Headings (PhySH)

Article Text

References (64)

  1. Á. Rivas and S. F. Huelga, Open Quantum Systems: An Introduction, SpringerBriefs in Physics (Springer, Berlin, 2012).
  2. P. Forn-Díaz, L. Lamata, E. Rico, J. Kono, and E. Solano, Ultrastrong coupling regimes of light-matter interaction, Rev. Mod. Phys. 91, 025005 (2019).
  3. A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, Dynamical mean-field theory of strongly correlated fermion systems and the limit of infinite dimensions, Rev. Mod. Phys. 68, 13 (1996).
  4. R. Härtle, G. Cohen, D. R. Reichman, and A. J. Millis, Decoherence and lead-induced interdot coupling in nonequilibrium electron transport through interacting quantum dots: A hierarchical quantum master equation approach, Phys. Rev. B 88, 235426 (2013).
  5. S. Mukamel and D. Abramavicius, Many-body approaches for simulating coherent nonlinear spectroscopies of electronic and vibrational excitons, Chem. Rev. 104, 2073 (2004).
  6. M. Tanaka and Y. Tanimura, Quantum dissipative dynamics of electron transfer reaction system: Nonperturbative hierarchy equations approach, J. Phys. Soc. Jpn. 78, 073802 (2009).
  7. B. M. Garraway, Nonperturbative decay of an atomic system in a cavity, Phys. Rev. A 55, 2290 (1997).
  8. A. Dorda, M. Nuss, W. von der Linden, and E. Arrigoni, Auxiliary master equation approach to nonequilibrium correlated impurities, Phys. Rev. B 89, 165105 (2014).
  9. A. Dorda, M. Ganahl, H. G. Evertz, W. von der Linden, and E. Arrigoni, Auxiliary master equation approach within matrix product states: Spectral properties of the nonequilibrium Anderson impurity model, Phys. Rev. B 92, 125145 (2015).
  10. M. P. Woods and M. B. Plenio, Dynamical error bounds for continuum discretisation via Gauss quadrature rules—A Lieb-Robinson bound approach, J. Math. Phys. 57, 022105 (2016).
  11. R. Rosenbach, J. Cerrillo, S. F. Huelga, J. Cao, and M. B. Plenio, Efficient simulation of non-Markovian system-environment interaction, New J. Phys. 18, 023035 (2016).
  12. F. Mascherpa, A. Smirne, S. F. Huelga, and M. B. Plenio, Open systems with error bounds: Spin-boson model with spectral density variations, Phys. Rev. Lett. 118, 100401 (2017).
  13. D. Tamascelli, A. Smirne, S. F. Huelga, and M. B. Plenio, Nonperturbative treatment of non-Markovian dynamics of open quantum systems, Phys. Rev. Lett. 120, 030402 (2018).
  14. D. Tamascelli, A. Smirne, J. Lim, S. F. Huelga, and M. B. Plenio, Efficient simulation of finite-temperature open quantum systems, Phys. Rev. Lett. 123, 090402 (2019).
  15. A. Nüßeler, I. Dhand, S. F. Huelga, and M. B. Plenio, Efficient simulation of open quantum systems coupled to a fermionic bath, Phys. Rev. B 101, 155134 (2020).
  16. Y. Tanimura, Numerically “exact” approach to open quantum dynamics: The hierarchical equations of motion (HEOM), J. Chem. Phys. 153, 020901 (2020).
  17. M. Xu, Y. Yan, Q. Shi, J. Ankerhold, and J. T. Stockburger, Taming quantum noise for efficient low temperature simulations of open quantum systems, Phys. Rev. Lett. 129, 230601 (2022).
  18. 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. B 110, 195148 (2024).
  19. J. Thoenniss, I. Vilkoviskiy, and D. A. Abanin, Efficient pseudomode representation and complexity of quantum impurity models, Phys. Rev. B 112, 155114 (2025).
  20. M. Lednev, F. J. García-Vidal, and J. Feist, Lindblad master equation capable of describing hybrid quantum systems in the ultrastrong coupling regime, Phys. Rev. Lett. 132, 106902 (2024).
  21. Z. Huang, G. Park, G. K.-L. Chan, and L. Lin, Coupled Lindblad pseudomode theory for simulating open quantum systems, Phys. Rev. Lett. 136, 090403 (2026).
  22. I. de Vega, U. Schollwöck, and F. Alexander Wolf, How to discretize a quantum bath for real-time evolution, Phys. Rev. B 92, 155126 (2015).
  23. A. Nüßeler, D. Tamascelli, A. Smirne, J. Lim, S. F. Huelga, and M. B. Plenio, Fingerprint and universal Markovian closure of structured bosonic environments, Phys. Rev. Lett. 129, 140604 (2022).
  24. P. Menczel, K. Funo, M. Cirio, N. Lambert, and F. Nori, Non-Hermitian pseudomodes for strongly coupled open quantum systems: Unravelings, correlations, and thermodynamics, Phys. Rev. Res. 6, 033237 (2024).
  25. L. Van Hove, The occurrence of singularities in the elastic frequency distribution of a crystal, Phys. Rev. 89, 1189 (1953).
  26. J. Kaye, K. Chen, and O. Parcollet, Discrete Lehmann representation of imaginary time Green’s functions, Phys. Rev. B 105, 235115 (2022).
  27. Y. Nakatsukasa, O. Sète, and L. N. Trefethen, The AAA algorithm for rational approximation, SIAM J. Sci. Comput. 40, A1494 (2018).
  28. C. Mejuto-Zaera, L. Zepeda-Núñez, M. Lindsey, N. Tubman, B. Whaley, and L. Lin, Efficient hybridization fitting for dynamical mean-field theory via semi-definite relaxation, Phys. Rev. B 101, 035143 (2020).
  29. H. Shinaoka and Y. Nagai, Sparse modeling of large-scale quantum impurity models with low symmetries, Phys. Rev. B 103, 045120 (2021).
  30. Z. Huang, E. Gull, and L. Lin, Robust analytic continuation of Green’s functions via projection, pole estimation, and semidefinite relaxation, Phys. Rev. B 107, 075151 (2023).
  31. L. Ying, Pole recovery from noisy data on imaginary axis, J. Sci. Comput. 92, 107 (2022).
  32. J. Kaye, Z. Huang, H. U. R. Strand, and D. Golež, Decomposing imaginary-time Feynman diagrams using separable basis functions: Anderson impurity model strong-coupling expansion, Phys. Rev. X 14, 031034 (2024).
  33. Z. Huang, D. Golez, H. U. R. Strand, and J. Kaye, Automated evaluation of imaginary time strong coupling diagrams by sum-of-exponentials hybridization fitting, SciPost Phys. 19, 121 (2025).
  34. D. Gazizova, L. Zhang, E. Gull, and J. P. F. LeBlanc, Feynman diagrammatics based on discrete pole representations: A path to renormalized perturbation theories, Phys. Rev. B 110, 075158 (2024).
  35. D. Gazizova, R. Farid, B. D. E. McNiven, I. Assi, E. G. Armstrong, and J. P. F. LeBlanc, Computation kernel for Feynman diagrams, Phys. Rev. B 112, 035172 (2025).
  36. L. Zhang and E. Gull, Minimal pole representation and controlled analytic continuation of Matsubara response functions, Phys. Rev. B 110, 035154 (2024).
  37. L. Zhang, Y. Yu, and E. Gull, Minimal pole representation and analytic continuation of matrix-valued correlation functions, Phys. Rev. B 110, 235131 (2024).
  38. L. Ying, Analytic continuation from limited noisy Matsubara data, J. Comput. Phys. 469, 111549 (2022).
  39. H. LaBollita, J. Kaye, and A. Hampel, Stabilizing the calculation of the self-energy in dynamical mean-field theory using constrained residual minimization, Phys. Rev. B 111, 115155 (2025).
  40. N. Sheng, A. Hampel, S. Beck, O. Parcollet, N. Wentzell, J. Kaye, and K. Chen, Low-rank Green’s function representations applied to dynamical mean-field theory, Phys. Rev. B 107, 245123 (2023).
  41. M. Caffarel and W. Krauth, Exact diagonalization approach to correlated fermions in infinite dimensions: Mott transition and superconductivity, Phys. Rev. Lett. 72, 1545 (1994).
  42. A. Liebsch and H. Ishida, Temperature and bath size in exact diagonalization dynamical mean field theory, J. Phys.: Condens. Matter 24, 053201 (2012).
  43. F. A. Wolf, A. Go, I. P. McCulloch, A. J. Millis, and U. Schollwöck, Imaginary-time matrix product state impurity solver for dynamical mean-field theory, Phys. Rev. X 5, 041032 (2015).
  44. A. Valenti, I. Park, A. Georges, A. J. Millis, and O. Parcollet, Neural-network quantum embedding solvers for correlated materials, arXiv:2603.15741.
  45. L. Zhang, A. Erpenbeck, Y. Yu, and E. Gull, Minimal pole representation for spectral functions, J. Chem. Phys. 162, 214111 (2025).
  46. E. Paprotzki and M. Eckstein, High-order strong-coupling expansion for x-ray absorption on a dynamically screened impurity, Phys. Rev. B 112, 205126 (2025).
  47. Z. Huang, L. Lin, G. Park, and Y. Zhu, Unified error bounds for perturbations of non-Markovian open quantum systems in Gaussian environments, arXiv:2411.08741.
  48. L. N. Trefethen and J. A. C. Weideman, The exponentially convergent trapezoidal rule, SIAM Rev. 56, 385 (2014).
  49. A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. Fisher, A. Garg, and W. Zwerger, Dynamics of the dissipative two-state system, Rev. Mod. Phys. 59, 1 (1987).
  50. M. Thorwart, E. Paladino, and M. Grifoni, Dynamics of the spin-boson model with a structured environment, Chem. Phys. 296, 333 (2004).
  51. F. K. Wilhelm, S. Kleff, and J. von Delft, The spin-boson model with a structured environment: A comparison of approaches, Chem. Phys. 296, 345 (2004).
  52. D. I. G. Bennett, P. Malý, C. Kreisbeck, R. van Grondelle, and A. Aspuru-Guzik, Mechanistic regimes of vibronic transport in a heterodimer and the design principle of incoherent vibronic transport in phycobiliproteins, J. Phys. Chem. Lett. 9, 2665 (2018).
  53. P. Siegle, I. Goychuk, P. Talkner, and P. Hänggi, Markovian embedding of non-Markovian superdiffusion, Phys. Rev. E 81, 011136 (2010).
  54. I. Goychuk, Viscoelastic subdiffusion: Generalized Langevin equation approach, Adv. Chem. Phys. 150, 187 (2012).
  55. R. Kupferman, Fractional kinetics in Kac–Zwanzig heat bath models, J. Stat. Phys. 114, 291 (2004).
  56. B. A. Dalton, C. Ayaz, H. Kiefer, A. Klimek, L. Tepper, and R. R. Netz, Fast protein folding is governed by memory-dependent friction, Proc. Natl. Acad. Sci. USA 120, e2220068120 (2023).
  57. G. Jung, M. Hanke, and F. Schmid, Iterative reconstruction of memory kernels, J. Chem. Theory Comput. 13, 2481 (2017).
  58. J. O. Daldrop, B. G. Kowalik, and R. R. Netz, External potential modifies friction of molecular solutes in water, Phys. Rev. X 7, 041065 (2017).
  59. I. Vilkoviskiy and D. A. Abanin, Bound on approximating non-Markovian dynamics by tensor networks in the time domain, Phys. Rev. B 109, 205126 (2024).
  60. K. Müller and W. T. Strunz, One-to-one correspondence between hierarchical equations of motion and pseudomodes for open quantum system dynamics, arXiv:2604.06466.
  61. Z. Huang, 2026, https://github.com/Hertz4/collections/tree/main/hyb_complexity.
  62. S. A. Adelman and J. D. Doll, Generalized Langevin equation approach for atom/solid-surface scattering: Collinear atom/harmonic chain model, J. Chem. Phys. 61, 4242 (1974).
  63. E. Lutz, Fractional Langevin equation, Phys. Rev. E 64, 051106 (2001).
  64. R. Roy and T. Kailath, ESPRIT-estimation of signal parameters via rotational invariance techniques, IEEE Trans. Acoust. Speech Signal Process. 37, 984 (1989).

Outline

Information

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation