Simple Flux-Side Formulation of State-Resolved Thermal Reaction Rates for Ring-Polymer Surface Hopping

Coherence in Chemistry: Foundations and Frontiers

Coherence refers to correlations in waves. Because matter has a wave-particle nature, it is unsurprising that coherence has deep connections with the most contemporary issues in chemistry research (e.g., energy harvesting, femtosecond spectroscopy, molecular qubits and more). But what does the word "coherence" really mean in the context of molecules and other quantum systems? We provide a review of key concepts, definitions, and methodologies, surrounding coherence phenomena in chemistry, and we describe how the terms "coherence" and "quantum coherence" refer to many different phenomena in chemistry. Moreover, we show how these notions are related to the concept of an interference pattern. Coherence phenomena are indeed complex, and ambiguous definitions may spawn confusion. By describing the many definitions and contexts for coherence in the molecular sciences, we aim to enhance understanding and communication in this broad and active area of chemistry.

Mean-Field Ring Polymer Rates Using a Population Dividing Surface

Mean-field ring polymer molecular dynamics offers a computationally efficient method for the simulation of reaction rates in multilevel systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, it is necessary to ensure reactive trajectories form kinked ring polymer configurations at the dividing surface. Building on this idea, we introduce a population difference coordinate and a reactive flux expression modified to only include contributions from kinked configurations. We test the accuracy of the resulting mean-field rate theory on a series of linear vibronic coupling model systems. We demonstrate that this new kMF-RP rate approach is efficient to implement and quantitatively accurate for models over a wide range of driving forces, coupling strengths, and temperatures.

A size-consistent multi-state mapping approach to surface hopping

We develop a multi-state generalization of the recently proposed mapping approach to surface hopping (MASH) for the simulation of electronically nonadiabatic dynamics. This new approach extends the original MASH method to be able to treat systems with more than two electronic states. It differs from previous approaches in that it is size consistent and rigorously recovers the original two-state MASH in the appropriate limits. We demonstrate the accuracy of the method by applying it to a series of model systems for which exact benchmark results are available, and we find that the method is well suited to the simulation of photochemical relaxation processes.

Recovering Marcus Theory Rates and Beyond without the Need for Decoherence Corrections: The Mapping Approach to Surface Hopping

It is well-known that fewest-switches surface hopping (FSSH) fails to correctly capture the quadratic scaling of rate constants with diabatic coupling in the weak-coupling limit, as expected from Fermi's golden rule and Marcus theory. To address this deficiency, the most widely used approach is to introduce a "decoherence correction", which removes the inconsistency between the wave function coefficients and the active state. Here we investigate the behavior of a new nonadiabatic trajectory method, called the mapping approach to surface hopping (MASH), on systems that exhibit an incoherent rate behavior. Unlike FSSH, MASH hops between active surfaces deterministically and can never have an inconsistency between the wave function coefficients and the active state. We show that MASH not only can describe rates for intermediate and strong diabatic coupling but also can accurately reproduce the results of Marcus theory in the golden-rule limit, without the need for a decoherence correction. MASH is therefore a significant improvement over FSSH in the simulation of nonadiabatic reactions.

On detailed balance in nonadiabatic dynamics: From spin spheres to equilibrium ellipsoids

Trajectory-based methods that propagate classical nuclei on multiple quantum electronic states are often used to simulate nonadiabatic processes in the condensed phase. A long-standing problem of these methods is their lack of detailed balance, meaning that they do not conserve the equilibrium distribution. In this article, we investigate ideas for restoring detailed balance in mixed quantum-classical systems by tailoring the previously proposed spin-mapping approach to thermal equilibrium. We find that adapting the spin magnitude can recover the correct long-time populations but is insufficient to conserve the full equilibrium distribution. The latter can however be achieved by a more flexible mapping of the spin onto an ellipsoid, which is constructed to fulfill detailed balance for arbitrary potentials. This ellipsoid approach solves the problem of negative populations that has plagued previous mapping approaches and can therefore be applied also to strongly asymmetric and anharmonic systems. Because it conserves the thermal distribution, the method can also exploit efficient sampling schemes used in standard molecular dynamics, which drastically reduces the number of trajectories needed for convergence. The dynamics does however still have mean-field character, as is observed most clearly by evaluating reaction rates in the golden-rule limit. This implies that although the ellipsoid mapping provides a rigorous framework, further work is required to find an accurate classical-trajectory approximation that captures more properties of the true quantum dynamics.

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.

Nonadiabatic transition paths from quantum jump trajectories

We present a means of studying rare reactive pathways in open quantum systems using transition path theory and ensembles of quantum jump trajectories. This approach allows for the elucidation of reactive paths for dissipative, nonadiabatic dynamics when the system is embedded in a Markovian environment. We detail the dominant pathways and rates of thermally activated processes and the relaxation pathways and photoyields following vertical excitation in a minimal model of a conical intersection. We find that the geometry of the conical intersection affects the electronic character of the transition state as defined through a generalization of a committor function for a thermal barrier crossing event. Similarly, the geometry changes the mechanism of relaxation following a vertical excitation. Relaxation in models resulting from small diabatic coupling proceeds through pathways dominated by pure dephasing, while those with large diabatic coupling proceed through pathways limited by dissipation. The perspective introduced here for the nonadiabatic dynamics of open quantum systems generalizes classical notions of reactive paths to fundamentally quantum mechanical processes.

Nonadiabatic instanton rate theory beyond the golden-rule limit

Fermi's golden rule describes the leading-order behaviour of the reaction rate as a function of the diabatic coupling. Its asymptotic (ℏ →0) limit is the semiclassical golden-rule instanton rate theory, which rigorously approximates nuclear quantum effects, lends itself to efficient numerical computation and gives physical insight into reaction mechanisms. However the golden rule by itself becomes insufficient as the strength of the diabatic coupling increases, so higher-order terms must be additionally considered. In this work we give a first-principles derivation of the next-order term beyond the golden rule, represented as a sum of three components. Two of them lead to new instanton pathways that extend the golden-rule case and, among other factors, account for the effects of recrossing on the full rate. The remaining component derives from the equilibrium partition function and accounts for changes in potential energy around the reactant and product wells due to diabatic coupling. The new semiclassical theory demands little computational effort beyond a golden-rule instanton calculation. It makes it possible to rigorously assess the accuracy of the golden-rule approximation and sets the stage for future work on general semiclassical nonadiabatic rate theories.

Equilibrium-nonequilibrium ring-polymer molecular dynamics for nonlinear spectroscopy

Two-dimensional Raman and hybrid terahertz-Raman spectroscopic techniques provide invaluable insight into molecular structures and dynamics of condensed-phase systems. However, corroborating experimental results with theory is difficult due to the high computational cost of incorporating quantum-mechanical effects in the simulations. Here, we present the equilibrium-nonequilibrium ring-polymer molecular dynamics (RPMD), a practical computational method that can account for nuclear quantum effects on the two-time response function of nonlinear optical spectroscopy. Unlike a recently developed approach based on the double Kubo transformed (DKT) correlation function, our method is exact in the classical limit, where it reduces to the established equilibrium-nonequilibrium classical molecular dynamics method. Using benchmark model calculations, we demonstrate the advantages of the equilibrium-nonequilibrium RPMD over classical and DKT-based approaches. Importantly, its derivation, which is based on the nonequilibrium RPMD, obviates the need for identifying an appropriate Kubo transformed correlation function and paves the way for applying real-time path-integral techniques to multidimensional spectroscopy.

Path Integrals for Nonadiabatic Dynamics: Multistate Ring Polymer Molecular Dynamics

This review focuses on a recent class of path-integral-based methods that simulate nonadiabatic dynamics in the condensed phase using only classical molecular dynamics trajectories in an extended phase space. Specifically, a semiclassical mapping protocol is used to derive an exact, continuous, Cartesian variable path-integral representation for the canonical partition function of a system in which multiple electronic states are coupled to nuclear degrees of freedom. Building on this exact statistical foundation, multistate ring polymer molecular dynamics methods are developed for the approximate calculation of real-time thermal correlation functions. The remarkable promise of these multistate ring polymer methods, their successful applications, and their limitations are discussed in detail.Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

Recent progress in approximate quantum dynamics methods for the study of proton-coupled electron transfer reactions

Proton-coupled electron transfer (PCET) reactions are ubiquitous natural processes at the heart of energy conversion reactions in photosynthesis and respiration, DNA repair, and diverse enzymatic reactions. Theoretical formulation and computational method developments have eyed modeling of thermal and photoinduced PCET for the last three decades. The accumulation of these studies, collected in dozens of reviews, accounts, and perspectives, has firmly established the influence of quantum effects, including non-adiabatic electronic transitions, vibrational relaxation, zero-point energy, and proton tunneling, on the rate and mechanism of PCET reactions. Here, we focus on some recently-developed methods, spanning the last eight years, that can quantitatively capture these effects in the PCET context and provide efficient means for their qualitative description in complex systems. The theoretical background of each method and their accuracy with respect to exact results are discussed and the results of relevant PCET simulations based on each method are presented.

An improved path-integral method for golden-rule rates

We present a simple method for the calculation of reaction rates in the Fermi golden-rule limit, which accurately captures the effects of tunneling and zero-point energy. The method is based on a modification of the recently proposed golden-rule quantum transition state theory (GR-QTST) of Thapa, Fang, and Richardson [J. Chem. Phys. 150, 104107 (2019)]. While GR-QTST is not size consistent, leading to the possibility of unbounded errors in the rate, our modified method has no such issue and so can be reliably applied to condensed phase systems. Both methods involve path-integral sampling in a constrained ensemble; the two methods differ, however, in the choice of constraint functional. We demonstrate numerically that our modified method is as accurate as GR-QTST for the one-dimensional model considered by Thapa and co-workers. We then study a multidimensional spin-boson model, for which our method accurately predicts the true quantum rate, while GR-QTST breaks down with an increasing number of boson modes in the discretization of the spectral density. Our method is able to accurately predict reaction rates in the Marcus inverted regime without the need for the analytic continuation required by Wolynes theory.

Nonadiabatic dynamics with quantum nuclei: simulating charge transfer with ring polymer surface hopping

Exploring effects of quantizing nuclei in non-adiabatic dynamics for simulating charge transfer in a dimer of “ethylene-like-molecules” at different temperatures.

Revisiting nuclear tunnelling in the aqueous ferrous–ferric electron transfer

We find that golden-rule quantum transition-state theory predicts nearly an order of magnitude less tunnelling than some of the previous estimates. This may indicate that the spin-boson model of electron transfer is not valid in the quantum regime.