Seeking a quantum advantage with trapped-ion quantum simulations of condensed-phase chemical dynamics.

Simulating the quantum dynamics of molecules in the condensed phase represents a longstanding challenge in chemistry. Trapped-ion quantum systems may serve as a platform for the analog-quantum simulation of chemical dynamics that is beyond the reach of current classical-digital simulation. To identify a 'quantum advantage' for these simulations, performance analysis of both analog-quantum simulation on noisy hardware and classical-digital algorithms is needed. In this Review, we make a comparison between a noisy analog trapped-ion simulator and a few choice classical-digital methods on simulating the dynamics of a model molecular Hamiltonian with linear vibronic coupling. We describe several simple Hamiltonians that are commonly used to model molecular systems, which can be simulated with existing or emerging trapped-ion hardware. These Hamiltonians may serve as stepping stones towards the use of trapped-ion simulators for systems beyond the reach of classical-digital methods. Finally, we identify dynamical regimes in which classical-digital simulations seem to have the weakest performance with respect to analog-quantum simulations. These regimes may provide the lowest hanging fruit to make the most of potential quantum advantages.

Non-adiabatic ring polymer molecular dynamics in the phase space of the SU(N) Lie group.

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.

Interference between Molecular and Photon Field-Mediated Electron Transfer Coupling Pathways in Cavities.

Cavity polaritonics creates novel opportunities to direct chemical reactions. Electron transfer (ET) reactions are among the simplest reactions, and they underpin energy conversion. New strategies to manipulate and direct electron flow at the nanoscale are of particular interest in biochemistry, energy science, bioinspired materials science, and chemistry. We show that optical cavities can modulate electron transfer pathway interferences and ET rates in donor-bridge-acceptor (DBA) systems. We derive the rate for DBA electron transfer when the molecules are coupled to cavity modes, emphasizing novel cavity-induced pathway interferences with the molecular electronic coupling pathways, as these interferences allow a new kind of ET rate tuning. The interference between the cavity-induced coupling pathways and the intrinsic molecular coupling pathway is dependent on the cavity properties. Thus, manipulating the interference between the cavity-induced DA coupling and the bridge-mediated coupling offers an approach to direct and manipulate charge flow at the nanoscale.

Non-adiabatic mapping dynamics in the phase space of the SU(N) Lie group.

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 su(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.

Non-adiabatic ring polymer molecular dynamics with spin mapping variables.

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.

Non-adiabatic Matsubara dynamics and non-adiabatic ring-polymer molecular dynamics.

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.

Ring polymer quantization of the photon field in polariton chemistry.

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.

State dependent ring polymer molecular dynamics for investigating excited nonadiabatic dynamics.

A recently proposed nonadiabatic ring polymer molecular dynamics (NRPMD) approach has shown to provide accurate quantum dynamics by incorporating explicit state descriptions and nuclear quantizations. Here, we present a rigorous derivation of the NRPMD Hamiltonian and investigate its performance on simulating excited state nonadiabatic dynamics. Our derivation is based on the Meyer-Miller-Stock-Thoss 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 the photoinduced nonadiabatic 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 state nonadiabatic dynamics. This work provides a solid theoretical foundation of the promising NRPMD Hamiltonian and demonstrates the possibility of using the state-dependent RPMD approach to accurately simulate electronic nonadiabatic dynamics while explicitly quantizing nuclei.

Coherent state mapping ring polymer molecular dynamics for non-adiabatic quantum propagations.

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.