Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The Journal of Chemical Physics

Exact open quantum system dynamics: Optimal frequency vs time representation of bath correlations

Richard Hartmann, Michael Werther, Frank Großmann, Walter T. Strunz

Published Jun 17, 2019• 29 citations
Publisher

Abstract

Two different numerically exact methods for open quantum system dynamics, the hierarchy of pure states (HOPS) method, and the multi-Davydov-Ansatz are discussed. We focus on the suitability of the underlying representations of bath correlations. While in the HOPS case the correct description of the bath correlation function (BCF) in the time domain is decisive, it turns out that a windowed Fourier transform of the BCF is an appropriate indicator of the quality of the discretization in the multi-Davydov-Ansatz. For the spin-boson model with sub-Ohmic spectral density considered here, a discretization of the spectral density based on an exponential distribution, used previously, turns out to be most favorable.

DOI
10.1063/1.5097158

Related papers

Exact Open Quantum System Dynamics Using the Hierarchy of Pure States (HOPS)

Richard Hartmann, Walter T. Strunz

We show that the general and numerically exact Hierarchy of Pure States method (HOPS) is very well applicable to calculate the reduced dynamics of an open quantum system. In particular, we focus on environments with a sub-Ohmic spectral density (SD) resulting in an algebraic decay of the bath correlation function (BCF). The universal applicability of HOPS, reaching from weak to strong coupling for zero and nonzero temperature, is demonstrated by solving the spin-boson model for which we find perfect agreement with other methods, each one suitable for a special regime of parameters. The challenges arising in the strong coupling regime are not only reflected in the computational effort needed for the HOPS method to converge but also in the necessity for an importance sampling mechanism, accounted for by the nonlinear variant of HOPS. In order to include nonzero-temperature effects in the strong coupling regime we found that it is highly favorable for the HOPS method to use the zero-temperature BCF and include temperature via a stochastic Hermitian contribution to the system Hamiltonian.

arXiv:1710.08268•Journal of Chemical Theory and Computation•Oct 10, 2017•73 citations

Testing bath correlation functions for open quantum dynamics simulations

Masaaki Tokieda

Accurate simulations of thermalization in open quantum systems require a reliable representation of the bath correlation function (BCF). Numerical approaches, such as the hierarchical equations of motion and the pseudomode method, inherently approximate the BCF using a finite set of functions, which can impact simulation accuracy. In this work, we propose a practical and rigorous testing framework to assess the validity of approximate BCFs in open quantum dynamics simulations. Our approach employs a harmonic oscillator system, where the computed dynamics can be benchmarked against known exact solutions. To enable practical testing, we make two key methodological advancements. First, we develop numerical techniques to efficiently evaluate these exact solutions across a wide range of BCFs, ensuring broad applicability. Second, we introduce a moment-based state representation that significantly simplifies computations by exploiting the Gaussian nature of the system. Applications to a two-spin system and a transmon-resonator system demonstrate that our testing procedure provides error estimates that capture the qualitative trends observed in thermalization simulations. Using this methodology, we assess the performance of recently proposed BCF construction methods, highlighting both their strengths and a notable challenge posed by sub-Ohmic spectral densities at finite temperatures.

arXiv:2504.08068•Physical Review Research•Nov 18, 2025•1 citation

High accuracy exponential decomposition of bath correlation functions for arbitrary and structured spectral densities: Emerging methodologies and new approaches

H. Takahashi, Samuel L. Rudge, Christoph Kaspar +2 more

This study investigates the decomposition of bath correlation functions (BCFs) in terms of complex exponential functions, with an eye on the realistic modeling of open quantum systems based on the hierarchical equations of motion. We introduce the theoretical background of various decomposition schemes in both time and frequency domains and assess their efficiency and accuracy by demonstrating the decomposition of various BCFs. We further develop a new procedure for the decomposition of BCFs originating from highly structured spectral densities with a high accuracy and compare it with existing fitting techniques. Advantages and disadvantages of each methodology are discussed in detail with special attention to their application to the corresponding quantum dynamical problem. This work provides fundamental tools for choosing and using a variety of decomposition techniques of BCFs for the study of open quantum systems in structured environments.

The Journal of Chemical Physics•May 22, 2024•45 citations

Dynamics of the sub-Ohmic spin-boson model: A comparison of three numerical approaches

Yao Yao, Liwei Duan, Zhiguo Lü +2 more

Dynamics of the sub-Ohmic spin-boson model is examined using three numerical approaches, namely the Dirac-Frenkel time-dependent variation with the Davydov D(1) ansatz, the adaptive time-dependent density matrix renormalization group method within the representation of orthogonal polynomials, and a perturbative approach based on a unitary transformation. In order to probe the validity regimes of the three approaches, we study the dynamics of a qubit coupled to a bosonic bath with and without a local field. Comparison of the up-state population evolution shows that the three approaches are in agreement in the weak-coupling regime but exhibit marked differences when the coupling strength is large. The Davydov D(1) ansatz and the time-dependent density matrix renormalization group can both be reliably employed in the weak-coupling regime, while the former is also valid in the strong-coupling regime as judged by how faithfully the trial state follows the Schrödinger equation. We further explore the bipartite entanglement dynamics between two qubits coupled with individual bosonic baths which reveals entanglement sudden death and revival.

Physical Review E•Aug 16, 2013•44 citations

Non-Markovian dynamics in nonstationary Gaussian baths: a hierarchy of pure states approach

Vladislav Sukharnikov, Stasis Chuchurka, Frank Schlawin

Building on the standard hierarchy of pure states (HOPS) approach, we construct a generalized formulation suitable for open quantum systems interacting with nonstationary Gaussian baths, potentially extending its applicability to nonequilibrium baths. This is achieved by extending the conventional exponential decomposition of a bath correlation function (BCF) for nonstationary cases. Using our formulation of HOPS, we derive the corresponding hierarchy of master equations and, when each term in the BCF expansion can be associated with an independent physical bath, we show how the formalism connects to the well-known pseudomode representation. We demonstrate the method's performance on two examples of nonstationary squeezed reservoirs generated via uniform squeezing and degenerate parametric amplification in a one-sided cavity. Benchmarking against the hierarchy of master equations shows that HOPS is more efficient under hierarchy truncation. The pseudomode representation is shown to be more efficient in the strongly non-Markovian regime. Our results highlight HOPS as a versatile and powerful tool for simulating open quantum systems in nonstationary baths, with potential applications ranging from squeezed light-matter interactions to driven quantum materials and dissipative phase transitions.

quant-pharXiv:2511.01358•Physical Review Research•Nov 4, 2025•0 citations

Discretization of Structured Bosonic Environments at Finite Temperature by Interpolative Decomposition: Theory and Application

H. Takahashi, Raffaele Borrelli

2024, 161, 151101]. The approach leverages a low-rank decomposition of the Fourier-transform relation connecting the bath correlation function to its spectral density. By capturing the time, frequency, and temperature dependencies encoded in the spectral density-autocorrelation function relation, our method significantly reduces the degrees of freedom required for simulating open quantum system dynamics. We benchmark our approach against existing methods and demonstrate its efficacy through applications to both simple models and a realistic electron transfer process in biological systems. Additionally, we show that this new approach can be effectively combined with the tensor-train formalism to investigate the quantum dynamics of systems interacting with complex non-Markovian environments. Finally, we provide a perspective on the selection and application of various spectral density discretization techniques.

arXiv:2412.13793•Journal of Chemical Theory and Computation•Feb 17, 2025•6 citations
Two different numerically exact methods for open quantum system dynamics, the hierarchy of pure states (HOPS) method, and the multi-Davydov-Ansatz are discussed. We focus on the suitability of the underlying representations of bath correlations. While in the HOPS case the correct description of the bath correlation function (BCF) in the time domain is decisive, it turns out that a windowed Fourier transform of the BCF is an appropriate indicator of the quality of the discretization in the multi-Davydov-Ansatz. For the spin-boson model with sub-Ohmic spectral density considered here, a discretization of the spectral density based on an exponential distribution, used previously, turns out to be most favorable.
We show that the general and numerically exact Hierarchy of Pure States method (HOPS) is very well applicable to calculate the reduced dynamics of an open quantum system. In particular, we focus on environments with a sub-Ohmic spectral density (SD) resulting in an algebraic decay of the bath correlation function (BCF). The universal applicability of HOPS, reaching from weak to strong coupling for zero and nonzero temperature, is demonstrated by solving the spin-boson model for which we find perfect agreement with other methods, each one suitable for a special regime of parameters. The challenges arising in the strong coupling regime are not only reflected in the computational effort needed for the HOPS method to converge but also in the necessity for an importance sampling mechanism, accounted for by the nonlinear variant of HOPS. In order to include nonzero-temperature effects in the strong coupling regime we found that it is highly favorable for the HOPS method to use the zero-temperature BCF and include temperature via a stochastic Hermitian contribution to the system Hamiltonian.
Accurate simulations of thermalization in open quantum systems require a reliable representation of the bath correlation function (BCF). Numerical approaches, such as the hierarchical equations of motion and the pseudomode method, inherently approximate the BCF using a finite set of functions, which can impact simulation accuracy. In this work, we propose a practical and rigorous testing framework to assess the validity of approximate BCFs in open quantum dynamics simulations. Our approach employs a harmonic oscillator system, where the computed dynamics can be benchmarked against known exact solutions. To enable practical testing, we make two key methodological advancements. First, we develop numerical techniques to efficiently evaluate these exact solutions across a wide range of BCFs, ensuring broad applicability. Second, we introduce a moment-based state representation that significantly simplifies computations by exploiting the Gaussian nature of the system. Applications to a two-spin system and a transmon-resonator system demonstrate that our testing procedure provides error estimates that capture the qualitative trends observed in thermalization simulations. Using this methodology, we assess the performance of recently proposed BCF construction methods, highlighting both their strengths and a notable challenge posed by sub-Ohmic spectral densities at finite temperatures.
This study investigates the decomposition of bath correlation functions (BCFs) in terms of complex exponential functions, with an eye on the realistic modeling of open quantum systems based on the hierarchical equations of motion. We introduce the theoretical background of various decomposition schemes in both time and frequency domains and assess their efficiency and accuracy by demonstrating the decomposition of various BCFs. We further develop a new procedure for the decomposition of BCFs originating from highly structured spectral densities with a high accuracy and compare it with existing fitting techniques. Advantages and disadvantages of each methodology are discussed in detail with special attention to their application to the corresponding quantum dynamical problem. This work provides fundamental tools for choosing and using a variety of decomposition techniques of BCFs for the study of open quantum systems in structured environments.
Dynamics of the sub-Ohmic spin-boson model is examined using three numerical approaches, namely the Dirac-Frenkel time-dependent variation with the Davydov D(1) ansatz, the adaptive time-dependent density matrix renormalization group method within the representation of orthogonal polynomials, and a perturbative approach based on a unitary transformation. In order to probe the validity regimes of the three approaches, we study the dynamics of a qubit coupled to a bosonic bath with and without a local field. Comparison of the up-state population evolution shows that the three approaches are in agreement in the weak-coupling regime but exhibit marked differences when the coupling strength is large. The Davydov D(1) ansatz and the time-dependent density matrix renormalization group can both be reliably employed in the weak-coupling regime, while the former is also valid in the strong-coupling regime as judged by how faithfully the trial state follows the Schrödinger equation. We further explore the bipartite entanglement dynamics between two qubits coupled with individual bosonic baths which reveals entanglement sudden death and revival.
Building on the standard hierarchy of pure states (HOPS) approach, we construct a generalized formulation suitable for open quantum systems interacting with nonstationary Gaussian baths, potentially extending its applicability to nonequilibrium baths. This is achieved by extending the conventional exponential decomposition of a bath correlation function (BCF) for nonstationary cases. Using our formulation of HOPS, we derive the corresponding hierarchy of master equations and, when each term in the BCF expansion can be associated with an independent physical bath, we show how the formalism connects to the well-known pseudomode representation. We demonstrate the method's performance on two examples of nonstationary squeezed reservoirs generated via uniform squeezing and degenerate parametric amplification in a one-sided cavity. Benchmarking against the hierarchy of master equations shows that HOPS is more efficient under hierarchy truncation. The pseudomode representation is shown to be more efficient in the strongly non-Markovian regime. Our results highlight HOPS as a versatile and powerful tool for simulating open quantum systems in nonstationary baths, with potential applications ranging from squeezed light-matter interactions to driven quantum materials and dissipative phase transitions.
2024, 161, 151101]. The approach leverages a low-rank decomposition of the Fourier-transform relation connecting the bath correlation function to its spectral density. By capturing the time, frequency, and temperature dependencies encoded in the spectral density-autocorrelation function relation, our method significantly reduces the degrees of freedom required for simulating open quantum system dynamics. We benchmark our approach against existing methods and demonstrate its efficacy through applications to both simple models and a realistic electron transfer process in biological systems. Additionally, we show that this new approach can be effectively combined with the tensor-train formalism to investigate the quantum dynamics of systems interacting with complex non-Markovian environments. Finally, we provide a perspective on the selection and application of various spectral density discretization techniques.