A second-order Kubo response theory-centroid approach to vibrational energy relaxation for single-mode excitations

Ring-polymer, centroid, and mean-field approximations to multi-time Matsubara dynamics

Based on a recently developed generalization of Matsubara dynamics to the multi-time realm, we present a formal derivation of multi-time generalizations of ring-polymer molecular dynamics, thermostatted ring-polymer molecular dynamics (TRPMD), centroid molecular dynamics (CMD), and mean-field Matsubara dynamics. Additionally, we analyze the short-time accuracy of each methodology. We find that for multi-time correlation functions of linear operators, (T)RPMD is accurate up to order t³, while CMD is only correct up to t, indicating a degradation in the accuracy of these methodologies with respect to the single-time counterparts. The present work provides a firm justification for the use of path-integral-based approximations for the calculation of multi-time correlation functions.

Inclusion of nuclear quantum effects for simulations of nonlinear spectroscopy

The computation and interpretation of nonlinear vibrational spectroscopy is of vital importance for understanding a wide range of dynamical processes in molecular systems. Here, we introduce an approach to evaluate multi-time response functions in terms of multi-time double symmetrized Kubo transformed thermal correlation functions. Furthermore, we introduce a multi-time extension of ring polymer molecular dynamics to evaluate these Kubo transforms. Benchmark calculations show that the approximations are useful for short times even for nonlinear operators, providing a consistent improvement over classical simulations of multi-time correlation functions. The introduced methodology thus provides a practical way of including nuclear quantum effects in multi-time response functions of non-linear optical spectroscopy.

A canonical averaging in the second-order quantized Hamilton dynamics

Quantized Hamilton dynamics (QHD) is a simple and elegant extension of classical Hamilton dynamics that accurately includes zero-point energy, tunneling, dephasing, and other quantum effects. Formulated as a hierarchy of approximations to exact quantum dynamics in the Heisenberg formulation, QHD has been used to study evolution of observables subject to a single initial condition. In present, we develop a practical solution for generating canonical ensembles in the second-order QHD for position and momentum operators, which can be mapped onto classical phase space in doubled dimensionality and which in certain limits is equivalent to thawed Gaussian. We define a thermal distribution in the space of the QHD-2 variables and show that the standard beta=1/kT relationship becomes beta'=2/kT in the high temperature limit due to an overcounting of states in the extended phase space, and a more complicated function at low temperatures. The QHD thermal distribution is used to compute total energy, kinetic energy, heat capacity, and other canonical averages for a series of quartic potentials, showing good agreement with the quantum results.