arXiv: Liouvillian Gauge Impacts HEOM Numerical Stability, Researchers Find

Researchers Tianchu Li and Andrés Montoya-Castillo report the discovery of a “continuous family of analytically equivalent Liouvillians” within complex-pole hierarchical equations of motion, or HEOM, a numerical method used to simulate open quantum systems. Despite representing the same underlying physics, these mathematically distinct arrangements exhibit different numerical stability, revealing a hidden influence on simulation results. The team identified “spectral divergence and non-normal error amplification” as two pathways to inaccuracy within HEOM, offering a more nuanced understanding of simulation failures. By optimizing this “gauge,” the researchers introduce “GO–HEOM,” which they say eliminates divergences in strongly coupled environments and extends simulations to previously inaccessible coupling strengths, enabling the study of a “delocalized-to-localized quantum phase transition.”

This means the specific mathematical formulation chosen can influence simulation outcomes, even if the underlying physics remains constant, a finding that challenges conventional assumptions about numerical precision in quantum dynamics. The team identified that these instabilities stem from two separate, yet interconnected, mechanisms; this nuanced understanding is crucial for improving the accuracy of modeling open quantum systems, where interactions with the environment are complex and often difficult to predict.

This understanding led to the development of GO–HEOM, a gauge optimization strategy designed to address these issues directly. The team reports that GO–HEOM eliminates divergences and broadens the range of accessible coupling strengths for simulations.