We introduce the coherent-state mapping ring polymer molecular dynamics (CS-RPMD), a new method that accurately describes electronic non-adiabatic dynamics with explicit nuclear quantization. This new approach is derived by using coherent-state mapping representation for the electronic degrees of freedom (DOF) and the ring-polymer path-integral representation for the nuclear DOF. The CS-RPMD Hamiltonian does not contain any inter-bead coupling term in the state-dependent potential and correctly describes electronic Rabi oscillations. A classical equation of motion is used to sample initial configurations and propagate the trajectories from the CS-RPMD Hamiltonian. At the time equivalent to zero, the quantum Boltzmann distribution (QBD) is recovered by reweighting the sampled distribution with an additional phase factor. In a special limit that there is one bead for mapping variables and multiple beads for nuclei, CS-RPMD satisfies detailed balance and preserves an approximate QBD. Numerical tests of this method with a two-state model system show very good agreement with exact quantum results over a broad range of electronic couplings.

We use the ring polymer (RP) representation to quantize the radiation field inside an optical cavity to investigate polariton quantum dynamics. Using a charge transfer model coupled to an optical cavity, we demonstrate that the RP quantization of the photon field provides accurate rate constants of the polariton mediated electron transfer reaction compared to Fermi’s golden rule. Because RP quantization uses extended phase space to describe the photon field, it significantly reduces the computational costs compared to the commonly used Fock state description of the radiation field. Compared to the other quasi-classical descriptions of the photon field, such as the classical Wigner based mean-field Ehrenfest model, the RP representation provides a much more accurate description of the polaritonic quantum dynamics because it alleviates the potential quantum distribution leakage problem associated with the photonic degrees of freedom (DOF). This work demonstrates the possibility of using the ring polymer description to treat the quantized radiation field in polariton chemistry, offering an accurate and efficient approach for future investigations in cavity quantum electrodynamics.

Recently proposed non-adiabatic ring polymer molecular dynamics (NRPMD) approach has shown to provide accurate quantum dynamics by incorporating explicit electronic state descriptions and nuclear quantizations. Here, we present a rigorous derivation of the NRPMD Hamiltonian and investigate its performance on simulating excited state non-adiabatic dynamics. Our derivation is based on the Meyer-Miller-Stock-Thoss (MMST) mapping representation for electronic states and the ring-polymer path-integral description for nuclei, resulting in the same Hamiltonian proposed in the original NRPMD approach. In addition, we investigate the accuracy of using NRPMD to simulate photoinduced non-adiabatic dynamics in simple model systems. These model calculations suggest that NRPMD can alleviate the zero-point energy leakage problem that is commonly encountered in the classical Wigner dynamics, and provide accurate excited states non-adiabatic dynamics. This work provides a solid theoretical foundation of the promising NRPMD Hamiltonian and demonstrates the possibility of using state-dependent RPMD approach to accurately simulate electronic nonadiabatic dynamics while explicitly quantize nuclei.

We present a new non-adiabatic ring polymer molecular dynamics (NRPMD) method based on the spin mapping formalism, which we refer to as the spin mapping NRPMD (SM-NRPMD) approach. We derive the path-integral partition function expression using the spin coherent state basis for the electronic states and the ring polymer formalism for the nuclear degrees of freedom. This partition function provides an efficient sampling of the quantum statistics. Using the basic properties of the Stratonovich–Weyl transformation, we further justify a Hamiltonian that we propose for the dynamical propagation of the coupled spin mapping variables and the nuclear ring polymer. The accuracy of the SM-NRPMD method is numerically demonstrated by computing the nuclear position and population auto-correlation functions of non-adiabatic model systems. The results obtained using the SM-NRPMD method agree very well with the numerically exact results. The main advantage of using the spin mapping variables over the harmonic oscillator mapping variables is numerically demonstrated, where the former provides nearly time-independent expectation values of physical observables for systems under thermal equilibrium. We also explicitly demonstrate that SM-NRPMD provides invariant dynamics upon various ways of partitioning the state-dependent and state-independent potentials.

We present the non-adiabatic Matsubara dynamics, a general framework for computing the time-correlation function (TCF) of electronically non-adiabatic systems. This new formalism is derived based on the generalized Kubo-transformed TCF using the Wigner representation for both the nuclear degrees of freedom and the electronic mapping variables. By dropping the non-Matsubara nuclear normal modes in the quantum Liouvillian and explicitly integrating these modes out from the expression of the TCF, we derived the non-adiabatic Matsubara dynamics approach. Further making the approximation to drop the imaginary part of the Matsubara Liouvillian and enforce the nuclear momentum integral to be real, we arrived at the non-adiabatic ring-polymer molecular dynamics (NRPMD) approach. We have further justified the capability of NRPMD for simulating the non-equilibrium TCF. This work provides the rigorous theoretical foundation for several recently proposed state-dependent RPMD approaches and offers a general framework for developing new non-adiabatic quantum dynamics methods in the future.