Search

We derive the non-adiabatic ring polymer molecular dynamics (RPMD) approach in the phase space of the SU(N) Lie Group. This method, which we refer to as the spin mapping non-adiabatic RPMD (SM-NRPMD), is based on the spin-mapping formalism for the electronic degrees of freedom (DOFs) and ring polymer path-integral description for the nuclear DOFs. Using the Stratonovich–Weyl transform for the electronic DOFs and the Wigner transform for the nuclear DOFs, we derived an exact expression of the Kubo-transformed time-correlation function (TCF). We further derive the spin mapping non-adiabatic Matsubara dynamics using the Matsubara approximation that removes the high frequency nuclear normal modes in the TCF and derive the SM-NRPMD approach from the non-adiabatic Matsubara dynamics by discarding the imaginary part of the Liouvillian. The SM-NRPMD method has numerical advantages compared to the original NRPMD method based on the Meyer–Miller–Stock–Thoss (MMST) mapping formalism due to a more natural mapping using the SU(N) Lie Group that preserves the symmetry of the original system. We numerically compute the Kubo-transformed position auto-correlation function and electronic population correlation function for three-state model systems. The numerical results demonstrate the accuracy of the SM-NRPMD method, which outperforms the original MMST-based NRPMD. We envision that the SM-NRPMD method will be a powerful approach to simulate electronic non-adiabatic dynamics and nuclear quantum effects accurately. [ABSTRACT FROM AUTHOR]

This is an I easier approach to implement into computer code i , because these equations [Eq. (95) of the paper] are simpler than the corresponding EOMs for { I i SB I n i sb , I i SB I n i sb }. CLARIFICATION ON THE NUMERICAL ALGORITHM USED TO PROPAGATE DYNAMICS We want to clarify the numerical algorithm we used to propagate the EOMs and generate all numerical results presented in the paper. In the above expressions, to compute HT ht , we use Eq. (C2) of the paper. [Extracted from the article]

We present the rigorous theoretical framework of the generalized spin mapping representation for non-adiabatic dynamics. Our work is based upon a new mapping formalism recently introduced by Runeson and Richardson [J. Chem. Phys. 152, 084110 (2020)], which uses the generators of the s u (N) Lie algebra to represent N discrete electronic states, thus preserving the size of the original Hilbert space. Following this interesting idea, the Stratonovich–Weyl transform is used to map an operator in the Hilbert space to a continuous function on the SU(N) Lie group, i.e., a smooth manifold which is a phase space of continuous variables. We further use the Wigner representation to describe the nuclear degrees of freedom and derive an exact expression of the time-correlation function as well as the exact quantum Liouvillian for the non-adiabatic system. Making the linearization approximation, this exact Liouvillian is reduced to the Liouvillian of several recently proposed methods, and the performance of this linearized method is tested using non-adiabatic models. We envision that the theoretical work presented here provides a rigorous and unified framework to formally derive non-adiabatic quantum dynamics approaches with continuous variables and connects the previous methods in a clear and concise manner. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]