Symmetry Adapted Hierarchical Equations of Motion for Exact Simulations of Large Polariton Systems
Abstract
Hierarchical equations of motion(HEOM) provide exact dynamics of open quantum systems coupled to harmonic baths, but their computational cost becomes prohibitive for systems with many independent local environments. In this work, we develop a symmetry-adapted HEOM formalism to significantly reduce the computational cost for the permutationally invariant Holstein-Tavis-Cummings (HTC) model. The method removes redundant information in two stages. First, all auxiliary density operators (ADOs) related only by relabeling identical molecules and their bath channels are replaced by a single canonical occupation-pattern representative. Second, molecules with the same local hierarchy occupation produce repeated matrix elements within each representative, allowing only the distinct complex variables to be propagated instead of the full (N + 1) × (N + 1) ADO matrices. The resulting matrix-free equations are evaluated using precomputed connections and molecular multiplicities. At fixed hierarchy depth L and number of bath correlation exponentials m, the number of canonical representatives becomes independent of the ensemble size for N ≥ L and the number of unique variables saturates for N ≥ L + 2. The formulation easily extends to multiple-exponential bath decompositions, arbitrary initial density operators, static disorders, and cavity loss. Our benchmarks reproduce conventional HEOM dynamics while requiring far fewer propagated variables and substantially less memory.
Related
- Nonlinear Freezing of Vibrational Polariton Transport via Mesoscale Simulations
- Comparing Symmetry-Adapted Perturbation Theory to Intermolecular Energy Decomposition Analysis Based on Absolutely Localized Molecular Orbitals
- Direct calculation of dipole moment surfaces through second-order symmetry-adapted perturbation theory for properties
- First-order Degenerate Symmetry-Adapted Perturbation Theory
- SALTED: a symmetry-adapted machine-learning program for predicting electron-densities in molecules and materials
BibTeX
@misc{mondal2026symmetry,
title = {{Symmetry Adapted Hierarchical Equations of Motion for Exact Simulations of Large Polariton Systems}},
author = {M. Elious Mondal and Pengfei Huo},
howpublished = {arXiv},
year = {2026},
month = aug,
doi = {10.48550/arxiv.2608.20657},
eprint = {2608.20657},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2608.20657},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Symmetry Adapted Hierarchical Equations of Motion for Exact Simulations of Large Polariton Systems
AU - M. Elious Mondal
AU - Pengfei Huo
PY - 2026
DA - 2026/08/30
PB - arXiv
DO - 10.48550/arxiv.2608.20657
UR - https://arxiv.org/abs/2608.20657
AB - Hierarchical equations of motion(HEOM) provide exact dynamics of open quantum systems coupled to harmonic baths, but their computational cost becomes prohibitive for systems with many independent local environments. In this work, we develop a symmetry-adapted HEOM formalism to significantly reduce the computational cost for the permutationally invariant Holstein-Tavis-Cummings (HTC) model. The method removes redundant information in two stages. First, all auxiliary density operators (ADOs) related only by relabeling identical molecules and their bath channels are replaced by a single canonical occupation-pattern representative. Second, molecules with the same local hierarchy occupation produce repeated matrix elements within each representative, allowing only the distinct complex variables to be propagated instead of the full (N + 1) × (N + 1) ADO matrices. The resulting matrix-free equations are evaluated using precomputed connections and molecular multiplicities. At fixed hierarchy depth L and number of bath correlation exponentials m, the number of canonical representatives becomes independent of the ensemble size for N ≥ L and the number of unique variables saturates for N ≥ L + 2. The formulation easily extends to multiple-exponential bath decompositions, arbitrary initial density operators, static disorders, and cavity loss. Our benchmarks reproduce conventional HEOM dynamics while requiring far fewer propagated variables and substantially less memory.
ER -