Abstract:Imaginary-time path integral (PI) is a rigorous tool to compute static properties at finite temperatures. However, the stiff PI internal modes poses a sampling challenge. This is commonly tackled using staging coordinates, in which the free particle (FP) term of the PI action is diagonalized. We introduce novel and simple staging coordinates that diagonalize the entire action of the harmonic oscillator (HO) model, rendering it efficiently applicable to systems with harmonic character, such as quantum oscillators and crystals. The method is not applicable to fluids or systems with imaginary modes. Unlike FP staging, the HO staging provides a unique treatment of the centroid mode. We provide implementation schemes for PIMC and PIMD in NVT ensemble. Sampling efficiency is assessed in terms of precision and accuracy of estimating the energy and heat capacity of a one-dimensional HO and an asymmetric anharmonic oscillator (AO). In PIMC, the HO coordinates propose collective moves that perfectly sample the HO contribution, then (for AO) the anharmonic term is sampled using standard Metropolis method. This results in a high acceptance rate and, hence, high precision, in comparison to the FP staging. In PIMD, the HO coordinates prescribe definitions for the fictitious masses, yielding equal frequencies when applied to HO model. This allows for larger time step sizes relative to standard staging, without affecting accuracy or integrator stability. We also present results using normal mode (NM) coordinates, based on both HO and FP models. While staging and NM coordinates show similar performance (for FP or HO), staging is computationally preferable due to its cheaper scaling with Trotter number. The enhanced sampling of HO coordinates open avenues for efficient estimation of nuclear quantum effects in more complex systems with harmonic character, such as real molecular bonds and quantum crystals.

Two-Tier Modular Anharmonic Small Matrix Path Integral with Composite Spin-Boson Baths

Essential vs Removable Bath Anharmonicity: Path Integral Results with Model Electronic-Vibrational Baths

Tree-based Implementation of the Small Matrix Path Integral for System-Bath Dynamics

Kink Sum for Long-Memory Small Matrix Path Integral Dynamics

Nonadiabatic Dynamics Via the Symmetrical Quasi-Classical Method in the Presence of Anharmonicity.

A Unified Stochastic Formulation of Dissipative Quantum Dynamics. II. Beyond Linear Response of Spin Baths.

Role of Bath-Induced Many-Body Interactions in the Dissipative Phases of the Su-Schrieffer-Heeger Model

Ab-Initio Approach to Many-Body Quantum Spin Dynamics

Path-integral methodology and simulations of quantum thermal transport: Full counting statistics approach

Accurate Calculation of Equilibrium Reduced Density Matrix for the System-Bath Model: A Multilayer Multiconfiguration Time-Dependent Hartree Approach and Its Comparison to a Multi-Electronic-state Path Integral Molecular Dynamics Approach

Hierarchical equations of motion for multiple baths (HEOM-MB) and their application to Carnot cycle

Generalized Path Integral Energy and Heat Capacity Estimators of Quantum Oscillators and Crystals using Harmonic Mapping

Open quantum dynamics theory for a complex subenvironment system with a quantum thermostat: Application to a spin heat bath

Harmonic Oscillator Staging Coordinates for Efficient Path Integral Simulations of Quantum Oscillators and Crystals

Quadratic Scaling Bosonic Path Integral Molecular Dynamics

Path integral spin dynamics for quantum paramagnets

Path integral Lindblad master equation through transfer tensor method & the generalized quantum master equation

An effective Hamiltonian for the simulation of open quantum molecular systems

A unified stochastic formulation of dissipative quantum dynamics. I. Generalized hierarchical equations

About the performance of perturbative treatments of the spin-boson dynamics within the hierarchical equations of motion approach

Spin-mapping approach for nonadiabatic molecular dynamics