Track synthesis, materials and catalysis as it is published.
Search papers, labs, and topics across Lattice.
The authors resolve the long-standing numerical instability in the hierarchical equations of motion (HEOM) caused by finite truncation in bosonic environments by mapping the hierarchy to an auxiliary Fock-space representation. This framework reveals that strong transient amplification stems from the non-normal algebraic structure of unbalanced hierarchy-raising terms in the generator. Applying a non-unitary similarity transformation rebalances these raising and lowering terms, successfully suppressing numerical blowups and enabling stable open quantum system simulations at unprecedentedly strong system-bath couplings with structured spectral densities.
Open quantum system simulations routinely suffer numerical blowups under strong coupling, but rebalancing the generator's non-normal algebraic structure via a similarity transformation finally unlocks these intractable regimes.
The hierarchical equations of motion (HEOM) approach is one of the most powerful methods avail- able to simulate the dynamics of open quantum systems. Although the HEOM approach has been shown to be stable and efficient in many parameter regimes and physical scenarios, it is also well known that the unavoidable finite truncation of the hierarchy can lead to fundamental numeri- cal instabilities for bosonic environments. In this work, we apply recently developed techniques to analyze these issues. We use an auxiliary Fock-space picture of the hierarchy to first identify the unbalanced hierarchy-raising terms responsible for the strong transient amplification associated with the non-normal algebraic structure of the HEOM generator. We then apply a non-unitary similarity transformation that converts these terms into balanced combinations of hierarchy-raising and hierarchy-lowering terms, thereby suppressing their numerical impact. We first demonstrate the efficacy of this approach via a spin-boson model with a simple Brownian oscillator spectral den- sity and subsequently extend the analysis to a structured spectral density, showing for both cases that the transformed HEOM approach allows simulations of much larger system-bath coupling than would otherwise be possible.