Search papers, labs, and topics across Lattice.
This paper investigates approximations to the path-integral radius of gyration, ${\mathcal R}^2(ω)$, within the Hierarchical Equations of Motion (HEOM) approach to open quantum dynamics, focusing on the treatment of Matsubara-decay terms. It demonstrates that the Ishizaki-Tanimura correction corresponds to separating smooth and Brownian contributions to ${\mathcal R}^2(ω)$ and proposes a modification for faster baths. The authors further introduce an 'A4' adaptation of the AAA algorithm to efficiently fit ${\mathcal R}^2(ω)$ to a sum of poles, leading to improved HEOM performance at low temperatures.
A simple tweak to the Ishizaki-Tanimura correction unlocks faster and more efficient Hierarchical Equations of Motion (HEOM) calculations, especially at low temperatures.
A major challenge in open quantum dynamics is the inclusion of Matsubara-decay terms in the memory kernel, which arise from the quantum-Boltzmann delocalisation of the bath modes. This delocalisation can be quantified by the radius of gyration squared ${\mathcal R}^2(ω)$ of the imaginary-time Feynman paths of the bath modes as a function of the frequency $ω$. In a Hierarchical Equations of Motion (HEOM) calculation with a Debye--Drude spectral density, ${\mathcal R}^2(ω)$ is the only quantity that is treated approximately (assuming convergence with respect to hierarchy depth). Here, we show that the well-known Ishizaki--Tanimura correction is equivalent to separating smooth from `Brownian' contributions to ${\mathcal R}^2(ω)$, and that modifying the correction leads to a more efficient HEOM in the case of fast baths. We also develop a simple `A4' adaptation of the `AAA' (Adaptive Antoulas--Anderson) algorithm in order to fit ${\mathcal R}^2(ω)$ to a sum over poles, which results in an extremely efficient implementation of the standard HEOM method at low temperatures.
