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

Flexible scheme to truncate the hierarchy of pure states

Zhang, P., C. D. Bentley, Alexander Eisfeld

Published Apr 7, 2018• 15 citations
OpenAlex

Abstract

The hierarchy of pure states (HOPS) is a wavefunction-based method that can be used for numerically modeling open quantum systems. Formally, HOPS recovers the exact system dynamics for an infinite depth of the hierarchy. However, truncation of the hierarchy is required to numerically implement HOPS. We want to choose a "good" truncation method, where by "good" we mean that it is numerically feasible to check convergence of the results. For the truncation approximation used in previous applications of HOPS, convergence checks are numerically challenging. In this work, we demonstrate the application of the "n-particle approximation" to HOPS. We also introduce a new approximation, which we call the "n-mode approximation." We then explore the convergence of these truncation approximations with respect to the number of equations required in the hierarchy in two exemplary problems: absorption and energy transfer of molecular aggregates. Published by AIP Publishing.

Related papers

Hierarchy of Stochastic Pure States for Open Quantum System Dynamics

Dieter Suess, Alexander Eisfeld, Walter T. Strunz

We derive a hierarchy of stochastic evolution equations for pure states (quantum trajectories) for open quantum system dynamics with non-Markovian structured environments. This hierarchy of pure states (HOPS) is generally applicable and provides the exact reduced density operator as an ensemble average over normalized states. The corresponding nonlinear equations are presented. We demonstrate that HOPS provides an efficient theoretical tool and apply it to the spin-boson model, the calculation of absorption spectra of molecular aggregates, and energy transfer in a photosynthetic pigment-protein complex.

arXiv:1402.4647•Physical Review Letters•Oct 8, 2014•222 citations

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

A low-rank hierarchical framework for the non-Markovian stochastic Schrödinger equation with convergence analysis

Zhuohan Zhang, Zhenning Cai

We propose and analyze a novel numerical framework for the non-Markovian stochastic Schrödinger equation (NMSSE) based on a low-rank approximation of the bath correlation functions. By decomposing the memory kernel into a finite-dimensional representation, we derive a truncated system of hierarchical equations that effectively balances computational tractability with physical fidelity. A rigorous convergence analysis is established for the hierarchical framework under mild assumptions. We demonstrate that our formulation serves as a mathematical generalization of the Hierarchy of Pure States (HOPS), encompassing it as a special case while offering a more flexible representation of non-Markovian effects. Numerical experiments across several benchmark models are presented to illustrate the validity and efficacy of the proposed method.

arXiv:2607.14689•Jul 16, 2026•0 citations

Non-Markovian stochastic Schrödinger equation: Matrix-product-state approach to the hierarchy of pure states

Xing Bao Gao, Jiajun Ren, Alexander Eisfeld +1 more

We derive a stochastic hierarchy of matrix product states (HOMPS) for non-Markovian dynamics in an open quantum system at finite temperature, which is numerically exact and efficient. HOMPS is obtained from the recently developed stochastic hierarchy of pure states (HOPS) by expressing HOPS in terms of formal creation and annihilation operators. The resulting stochastic first-order differential equation is then formulated in terms of matrix product states and matrix product operators. In this way the exponential complexity of HOPS can be reduced to scale polynomial with the number of particles. The validity and efficiency of HOMPS is demonstrated for the spin-boson model and long chains where each site is coupled to a structured, strongly non-Markovian environment.

arXiv:2109.06393•Physical Review A•Mar 31, 2022•35 citations

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

Richard Hartmann, Michael Werther, Frank Großmann +1 more

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.

The Journal of Chemical Physics•Jun 17, 2019•29 citations

Closures of the functional expansion hierarchy in the non-Markovian quantum state diffusion approach

Gerhard Ritschel, Walter T. Strunz, Alexander Eisfeld

To find a practical scheme to numerically solve the non-Markovian Quantum State Diffusion equation (NMQSD), one often uses a functional expansion of the functional derivative that appears in the general NMQSD equation. This expansion leads to a hierarchy of coupled operators. It turned out that if one takes only the zeroth order term into account, one has a very efficient method that agrees remarkably well with the exact results for many cases of interest. We denote this approach as zeroth order functional expansion (ZOFE). In the present work, we investigate two extensions of ZOFE. Firstly, we investigate how the hierarchy converges when taking higher orders into account (which, however, leads to a fast increase in numerical size). Secondly, we demonstrate that by using a terminator that approximates the higher order contributions, one can obtain significant improvement, at hardly any additional computational cost. We carry out our investigations for the case of absorption spectra of molecular aggregates.

The Journal of Chemical Physics•Aug 14, 2017•9 citations
The hierarchy of pure states (HOPS) is a wavefunction-based method that can be used for numerically modeling open quantum systems. Formally, HOPS recovers the exact system dynamics for an infinite depth of the hierarchy. However, truncation of the hierarchy is required to numerically implement HOPS. We want to choose a "good" truncation method, where by "good" we mean that it is numerically feasible to check convergence of the results. For the truncation approximation used in previous applications of HOPS, convergence checks are numerically challenging. In this work, we demonstrate the application of the "n-particle approximation" to HOPS. We also introduce a new approximation, which we call the "n-mode approximation." We then explore the convergence of these truncation approximations with respect to the number of equations required in the hierarchy in two exemplary problems: absorption and energy transfer of molecular aggregates. Published by AIP Publishing.
We derive a hierarchy of stochastic evolution equations for pure states (quantum trajectories) for open quantum system dynamics with non-Markovian structured environments. This hierarchy of pure states (HOPS) is generally applicable and provides the exact reduced density operator as an ensemble average over normalized states. The corresponding nonlinear equations are presented. We demonstrate that HOPS provides an efficient theoretical tool and apply it to the spin-boson model, the calculation of absorption spectra of molecular aggregates, and energy transfer in a photosynthetic pigment-protein complex.
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.
We propose and analyze a novel numerical framework for the non-Markovian stochastic Schrödinger equation (NMSSE) based on a low-rank approximation of the bath correlation functions. By decomposing the memory kernel into a finite-dimensional representation, we derive a truncated system of hierarchical equations that effectively balances computational tractability with physical fidelity. A rigorous convergence analysis is established for the hierarchical framework under mild assumptions. We demonstrate that our formulation serves as a mathematical generalization of the Hierarchy of Pure States (HOPS), encompassing it as a special case while offering a more flexible representation of non-Markovian effects. Numerical experiments across several benchmark models are presented to illustrate the validity and efficacy of the proposed method.
We derive a stochastic hierarchy of matrix product states (HOMPS) for non-Markovian dynamics in an open quantum system at finite temperature, which is numerically exact and efficient. HOMPS is obtained from the recently developed stochastic hierarchy of pure states (HOPS) by expressing HOPS in terms of formal creation and annihilation operators. The resulting stochastic first-order differential equation is then formulated in terms of matrix product states and matrix product operators. In this way the exponential complexity of HOPS can be reduced to scale polynomial with the number of particles. The validity and efficiency of HOMPS is demonstrated for the spin-boson model and long chains where each site is coupled to a structured, strongly non-Markovian environment.
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.
To find a practical scheme to numerically solve the non-Markovian Quantum State Diffusion equation (NMQSD), one often uses a functional expansion of the functional derivative that appears in the general NMQSD equation. This expansion leads to a hierarchy of coupled operators. It turned out that if one takes only the zeroth order term into account, one has a very efficient method that agrees remarkably well with the exact results for many cases of interest. We denote this approach as zeroth order functional expansion (ZOFE). In the present work, we investigate two extensions of ZOFE. Firstly, we investigate how the hierarchy converges when taking higher orders into account (which, however, leads to a fast increase in numerical size). Secondly, we demonstrate that by using a terminator that approximates the higher order contributions, one can obtain significant improvement, at hardly any additional computational cost. We carry out our investigations for the case of absorption spectra of molecular aggregates.