On the electronic path integral normal modes of the Meyer-Miller-Stock-Thoss representation of nonadiabatic dynamics

Abstract

Accurate and efficient simulation of nonadiabatic dynamics is highly desirable for understanding charge and energy transfer in complex systems. A key criterion for obtaining an accurate method is conservation of the Quantum Boltzmann Distribution (QBD). For a single surface, Matsubara dynamics is known to conserve the QBD, as a consequence of truncating the dynamics in the higher normal modes of the imaginary-time path integral. Recently, a nonadiabatic Matsubara (NA-Mats) dynamics has been proposed (J. Chem. Phys., 2021, 154, 124124) which truncates in the normal modes of the nuclear variables but not in the electronic variables, which are described with the Meyer-Miller-Stock-Thoss (MMST) representation. Surprisingly, this NA-Mats method does not appear to conserve the QBD for a general system. This poses the question of the effect of truncating the higher path integral normal modes of the electronic variables in the MMST representation. In this article, we present what we believe is the first study of electronic normal modes of the MMST representation. We find that observables are not usually a function of a finite number of normal modes and that the higher normal modes are not constrained by the distribution, unlike in conventional nuclear normal modes. Furthermore, truncating the dynamics in MMST normal modes leads to inaccurate correlation functions and while the QBD appears conserved for an ensemble of trajectories, it is not for a single trajectory. Overall, this suggests that MMST path integral normal modes are not optimal for obtaining an accurate, QBD conserving nonadiabatic dynamics method.

14 pages, 7 figures, submitted to JCP, SI available as ancillary file

References in corpus (11)

• Efficient stochastic thermostatting of path integral molecular dynamics
• How to remove the spurious resonances from ring polymer molecular dynamics
• Boltzmann-conserving classical dynamics in quantum time-correlation functions: Matsubara dynamics
• Generalized spin mapping for quantum-classical dynamics
• Relation of centroid molecular dynamics and ring-polymer molecular dynamics to exact quantum dynamics
• A mapping approach to surface hopping
• Nonadiabatic semiclassical dynamics in the mixed quantum-classical initial value representation
• Exciton dynamics from the mapping approach to surface hopping: Comparison with Förster and Redfield theories
• On detailed balance in nonadiabatic dynamics: From spin spheres to equilibrium ellipsoids
• Which Algorithm Best Propagates the Meyer-Miller-Stock-Thoss Mapping Hamiltonian for Non-Adiabatic Dynamics?
• Time-reversible implementation of MASH for efficient nonadiabatic molecular dynamics