Quantum-classical Liouville dynamics of nonadiabatic proton transfer

Real-Time Dynamics and Detailed Balance in Ring Polymer Surface Hopping: The Impact of Frustrated Hops

Ring polymer surface hopping (RPSH) has been recently introduced as a well-tailored method for incorporating nuclear quantum effects, such as zero-point energy and tunneling, into nonadiabatic molecular dynamics simulations. The practical widespread usage of RPSH demands a comprehensive benchmarking of different reaction regimes and conditions with equal emphasis on demonstrating both the cons and the pros of the method. Here, we investigate the fundamental questions related to the conservation of energy and detailed balance in the context of RPSH. Using Tully's avoided crossing model as well as a 2-state quantum system coupled to a classical bath undergoing Langevin dynamics, we probe the critical problem of the proper treatment of the classically forbidden transitions stemming from the surface hopping algorithm. We show that proper treatment of these frustrated hops is key to the accurate description of real-time dynamics as well as reproducing the correct quantum Boltzmann populations.

DECIDE: A Deterministic Mixed Quantum-Classical Dynamics Approach

Mixed quantum-classical dynamics provides an efficient way of simulating the dynamics of quantum subsystems coupled to many-body environments. Many processes, including proton-transfer reactions, electron-transfer reactions, and vibrational energy transport, for example, take place in such open systems. The most accurate algorithms for performing mixed quantum-classical simulations require very large ensembles of trajectories to obtain converged expectation values, which is computationally prohibitive for quantum subsystems containing even a few degrees of freedom. The recently developed “Deterministic evolution of coordinates with initial decoupled equations” (DECIDE) method has demonstrated high accuracy and low computational cost for a host of model systems; however, these applications relied on representing the equations of motion in subsystem and adiabatic energy bases. While these representations are convenient for certain systems, the position representation is convenient for many other systems, including systems undergoing proton- and electron-transfer reactions. Thus, in this review, we provide a step-by-step derivation of the DECIDE approach and demonstrate how to cast the DECIDE equations in a quantum harmonic oscillator position basis for a simple one-dimensional (1D) hydrogen bond model. After integrating the DECIDE equations of motion on this basis, we show that the total energy of the system is conserved for this model and calculate various quantities of interest. Limitations of casting the equations in an incomplete basis are also discussed.

Theory of vibrational polariton chemistry in the collective coupling regime

We theoretically demonstrate that the chemical reaction rate constant can be significantly suppressed by coupling molecular vibrations with an optical cavity, exhibiting both the collective coupling effect and the cavity frequency modification of the rate constant. When a reaction coordinate is strongly coupled to the solvent molecules, the reaction rate constant is reduced due to the dynamical caging effect. We demonstrate that collectively coupling the solvent to the cavity can further enhance this dynamical caging effect, leading to additional suppression of the chemical kinetics. This effect is further amplified when cavity loss is considered.

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.

Quasi-Lie Brackets and the Breaking of Time-Translation Symmetry for Quantum Systems Embedded in Classical Baths

Many open quantum systems encountered in both natural and synthetic situations are embedded in classical-like baths. Often, the bath degrees of freedom may be represented in terms of canonically conjugate coordinates, but in some cases they may require a non-canonical or non-Hamiltonian representation. Herein, we review an approach to the dynamics and statistical mechanics of quantum subsystems embedded in either non-canonical or non-Hamiltonian classical-like baths which is based on operator-valued quasi-probability functions. These functions typically evolve through the action of quasi-Lie brackets and their associated Quantum-Classical Liouville Equations, or through quasi-Lie brackets augmented by dissipative terms. Quasi-Lie brackets possess the unique feature that, while conserving the energy (which the Noether theorem links to time-translation symmetry), they violate the time-translation symmetry of their algebra. This fact can be heuristically understood in terms of the dynamics of the open quantum subsystem. We then describe an example in which a quantum subsystem is embedded in a bath of classical spins, which are described by non-canonical coordinates. In this case, it has been shown that an off-diagonal open-bath geometric phase enters into the propagation of the quantum-classical dynamics. Next, we discuss how non-Hamiltonian dynamics may be employed to generate the constant-temperature evolution of phase space degrees of freedom coupled to the quantum subsystem. Constant-temperature dynamics may be generated by either a classical Langevin stochastic process or a Nosé–Hoover deterministic thermostat. These two approaches are not equivalent but have different advantages and drawbacks. In all cases, the calculation of the operator-valued quasi-probability function allows one to compute time-dependent statistical averages of observables. This may be accomplished in practice using a hybrid Molecular Dynamics/Monte Carlo algorithms, which we outline herein.

Efficient and Deterministic Propagation of Mixed Quantum-Classical Liouville Dynamics

We propose a highly efficient mixed quantum-classical molecular dynamics scheme based on a solution of the quantum-classical Liouville equation (QCLE). By casting the equations of motion for the quantum subsystem and classical bath degrees of freedom onto an approximate set of coupled first-order differential equations for c-numbers, this scheme propagates the composite system in time deterministically in terms of independent classical-like trajectories. To demonstrate its performance, we apply the method to the spin-boson model, a photoinduced electron transfer model, and a Fenna-Matthews-Olsen complex model, and find excellent agreement out to long times with the numerically exact results, using several orders of magnitude fewer trajectories than surface-hopping solutions of the QCLE. Owing to its accuracy and efficiency, this method promises to be very useful for studying the dynamics of mixed quantum-classical systems.

Ultrafast Reactive Quantum Dynamics Coupled to Classical Solvent Dynamics Using an Ehrenfest Approach

The inclusion of solvent effects in the theoretical analysis of molecular processes becomes increasingly important. Currently, it is not feasible to directly include the solvent on the quantum level. We use an Ehrenfest approach to study the coupled time evolution of quantum dynamically treated solutes and classical solvents system. The classical dynamics of the solvent is coupled to the wavepacket dynamics of the solute and rotational and translational degrees of freedom of the solute are included classically. This allows quantum dynamics simulations for ultrafast processes that are decided by environment interactions without explicit separation of time scales. We show the application to the dissociation of ICN in liquid Ar as a proof of principal system and to the more applied example of uracil in water.

Accurate Long-Time Mixed Quantum-Classical Liouville Dynamics via the Transfer Tensor Method

In this Letter, we combine the recently introduced transfer tensor method with the mixed quantum-classical Liouville method. The resulting protocol provides an accurate, general, flexible and robust new route for simulating the reduced dynamics of the quantum subsystem for arbitrarily long times, starting with computationally feasible short-time mixed quantum-classical Liouville dynamical maps. The accuracy and feasibility of the methodology are demonstrated on a spin-boson benchmark model.

Quantum dynamics in open quantum-classical systems

Often quantum systems are not isolated and interactions with their environments must be taken into account. In such open quantum systems these environmental interactions can lead to decoherence and dissipation, which have a marked influence on the properties of the quantum system. In many instances the environment is well-approximated by classical mechanics, so that one is led to consider the dynamics of open quantum-classical systems. Since a full quantum dynamical description of large many-body systems is not currently feasible, mixed quantum-classical methods can provide accurate and computationally tractable ways to follow the dynamics of both the system and its environment. This review focuses on quantum-classical Liouville dynamics, one of several quantum-classical descriptions, and discusses the problems that arise when one attempts to combine quantum and classical mechanics, coherence and decoherence in quantum-classical systems, nonadiabatic dynamics, surface-hopping and mean-field theories and their relation to quantum-classical Liouville dynamics, as well as methods for simulating the dynamics.

Can we derive Tully's surface-hopping algorithm from the semiclassical quantum Liouville equation? Almost, but only with decoherence

In this article, we demonstrate that Tully's fewest-switches surface hopping (FSSH) algorithm approximately obeys the mixed quantum-classical Liouville equation (QCLE), provided that several conditions are satisfied--some major conditions, and some minor. The major conditions are: (1) nuclei must be moving quickly with large momenta; (2) there cannot be explicit recoherences or interference effects between nuclear wave packets; (3) force-based decoherence must be added to the FSSH algorithm, and the trajectories can no longer rigorously be independent (though approximations for independent trajectories are possible). We furthermore expect that FSSH (with decoherence) will be most robust when nonadiabatic transitions in an adiabatic basis are dictated primarily by derivative couplings that are presumably localized to crossing regions, rather than by small but pervasive off-diagonal force matrix elements. In the end, our results emphasize the strengths of and possibilities for the FSSH algorithm when decoherence is included, while also demonstrating the limitations of the FSSH algorithm and its inherent inability to follow the QCLE exactly.

Signatures of nanoconfinement on the linear and nonlinear vibrational spectroscopy of a model hydrogen-bonded complex dissolved in a polar solvent

The one-dimensional IR (1D-IR) absorption and IR pump-probe spectra of a hydrogen stretch in a model hydrogen-bonded complex dissolved in a polar solvent confined in spherical hydrophobic cavities of different sizes were simulated using ground-state mixed quantum-classical dynamics. Due to a thorough analysis of key properties of the complex and solvent from equilibrium trajectory data, we were able to gain insight into the microscopic details underlying the spectra. Both the 1D-IR and IR pump-probe spectra manifested the effects of confinement on the relative stabilities of the covalent and ionic forms of the complex through pronounced changes in their peak intensities and numbers. However, in contrast to the 1D-IR spectra, the time-resolved pump-probe spectra were found to be uniquely sensitive to the changes in the molecular dynamics as the cavity size is varied. In particular, it was found that the variations in the time evolutions of the peak intensities in the pump-probe spectra reflect the differences in the solvation dynamics associated with the various forms of the complex in different locations within the cavities. The ability to detect these differences underscores the advantage of using pump-probe spectroscopy for studying nanoconfined systems.

Computational study of the one- and two-dimensional infrared spectra of a proton-transfer mode in a hydrogen-bonded complex dissolved in a polar nanocluster

The signatures of nanosolvation on the one- and two-dimensional (1D and 2D) IR spectra of a proton-transfer mode in a hydrogen-bonded complex dissolved in polar solvent molecule nanoclusters of varying size are elucidated by using mixed quantum-classical molecular dynamics simulations. For this particular system, increasing the number of solvent molecules successively from N=7 to N=9 initiates the transition of the system from a cluster state to a bulk-like state. Both the 1D and 2D IR spectra reflect this transition through pronounced changes in their peak intensities and numbers, but the time-resolved 2D IR spectra also manifest spectral features that uniquely identify the onset of the cluster-to-bulk transition. In particular, it is observed that in the 1D IR spectra, the relative intensities of the peaks change such that the number of peaks decreases from three to two as the size of the cluster increases from N=7 to N=9. In the 2D IR spectra, off-diagonal peaks are observed in the N=7 and N=8 cases at zero waiting time, but not in the N=9 case. It is known that there are no off-diagonal peaks in the 2D IR spectrum of the bulk version of this system at zero waiting time, so the disappearance of these peaks is a unique signature of the onset of bulk-like behavior. Through an examination of the trajectories of various properties of the complex and solvent, it is possible to relate the emergence of these off-diagonal peaks to an interplay between the vibrations of the complex and the solvent polarization dynamics.

Filtering schemes in the quantum-classical Liouville approach to nonadiabatic dynamics

We study a number of filtering schemes for the reduction of the statistical error in nonadiabatic calculations by means of the quantum-classical Liouville equation. In particular, we focus on a scheme based on setting a threshold value on the sampling weights, so that when the threshold is overcome the value of the weight is reset, and on another approach which prunes the ensemble of the allowed nonadiabatic transitions according to a generalized sampling probability. Both methods have advantages and drawbacks, however, their combination drastically improves the performance of an algorithm known as the sequential short-time step propagation [MacKernan et al., J. Phys: Condens. Matter 14, 9069 (2002)], which is derived from a simple first order expansion of the quantum-classical propagator. Such an algorithm together with the combined filtering procedures produce results that compare very well with those obtained by means of numerically accurate path integral quantum calculations for the spin-boson model, even for intermediate and strong coupling regimes.

Dynamical steps that bridge piecewise adiabatic shapes in the exact time-dependent potential energy surface

We study the exact time-dependent potential energy surface (TDPES) in the presence of strong nonadiabatic coupling between the electronic and nuclear motion. The concept of the TDPES emerges from the exact factorization of the full electron-nuclear wave function [A. Abedi, N. T. Maitra, and E. K. U. Gross, Phys. Rev. Lett. 105, 123002 (2010)]. Employing a one-dimensional model system, we show that the TDPES exhibits a dynamical step that bridges between piecewise adiabatic shapes. We analytically investigate the position of the steps and the nature of the switching between the adiabatic pieces of the TDPES.

Mixed quantum-classical molecular dynamics study of the hydroxyl stretch in methanol/carbon-tetrachloride mixtures II: excited state hydrogen bonding structure and dynamics, infrared emission spectrum, and excited state lifetime

We present a mixed quantum-classical molecular dynamics study of the hydrogen-bonding structure and dynamics of a vibrationally excited hydroxyl stretch in methanol/carbon-tetrachloride mixtures. The adiabatic Hamiltonian of the quantum-mechanical hydroxyl is diagonalized on-the-fly to obtain the ground and first-excited adiabatic energy levels and wave functions which depend parametrically on the instantaneous configuration of the classical degrees of freedom. The dynamics of the classical degrees of freedom are determined by Hellmann-Feynman forces obtained by taking the expectation value of the force with respect to the ground or excited vibrational wave functions. Polarizable force fields are used which were previously shown to reproduce the experimental infrared absorption spectrum rather well, for different isotopomers and over a wide composition range [Kwac, K.; Geva, E. J. Phys. Chem. B 2011, 115, 9184]. We show that the agreement of the absorption spectra with experiment can be further improved by accounting for the dependence of the dipole moment derivatives on the configuration of the classical degrees of freedom. We find that the propensity of a methanol molecule to form hydrogen bonds increases upon photoexcitation of its hydroxyl stretch, thereby leading to a sizable red-shift of the corresponding emission spectrum relative to the absorption spectrum. Treating the relaxation from the first excited to the ground state as a nonadiabatic process, and calculating its rate within the framework of Fermi's golden rule and the harmonic-Schofield quantum correction factor, we were able to predict a lifetime which is of the same order of magnitude as the experimental value. The experimental dependence of the lifetime on the transition frequency is also reproduced. Nonlinear mapping relations between the hydroxyl transition frequency and bond length in the excited state and the electric field along the hydroxyl bond axis are established. These mapping relations make it possible to reduce the computational cost of the mixed quantum-classical treatment to that of a fully classical treatment.

TOPOLOGICAL STUDY OF REACTIVE POTENTIAL-ENERGY SURFACES FROM MODEL DIABATIC STATES COUPLED TO SURROUNDING MEDIA

Reaction paths are studied with the help of diabatic potential-energy surfaces coupled in a generic external field. We show that all putative geometrical and topological features of two-dimensional (2D) potential-energy surfaces for an isomerization can be generated with a model consisting of strictly diabatic harmonic 2D wells coupled in a static (yet uniform) external field. For a large class of three-state models, we provide a phase diagram of possible topologies in the relevant parameter space for three-state models. By following the shift in diabatic electronic amplitudes (denoted by {|ck|2}) in coherent quantum states produced in the field, we assess whether the models within each phase can produce reaction paths visiting regions of configurational space that are structurally similar to the diabatic transition state.

Finite temperature application of the corrected propagator method to reactive dynamics in a condensed-phase environment

The recently proposed mixed quantum-classical method is extended to applications at finite temperatures. The method is designed to treat complex systems consisting of a low-dimensional quantum part (the primary system) coupled to a dissipative bath described classically. The method is based on a formalism showing how to systematically correct the approximate zeroth-order evolution rule. The corrections are defined in terms of the total quantum Hamiltonian and are taken to the classical limit by introducing the frozen Gaussian approximation for the bath degrees of freedom. The evolution of the primary system is governed by the corrected propagator yielding the exact quantum dynamics. The method has been tested on a standard model system describing proton transfer in a condensed-phase environment: a symmetric double-well potential bilinearly coupled to a bath of harmonic oscillators. Flux correlation functions and thermal rate constants have been calculated at two different temperatures for a range of coupling strengths. The results have been compared to the fully quantum simulations of Topaler and Makri [J. Chem. Phys. 101, 7500 (1994)] with the real path integral method.