Electronically nonadiabatic trajectories: Continuous surface switching II

Exciton energy transfer inside cavity-A benchmark study of polaritonic dynamics using the surface hopping method

Strong coupling between the molecular system and photon inside the cavity generates polaritons, which can alter reaction rates by orders of magnitude. In this work, we benchmark the surface hopping method to simulate non-adiabatic dynamics in a cavity. The comparison is made against a numerically exact method (the hierarchical equations of motion) for a model system investigating excitonic energy transfer for a broad range of parameters. Surface hopping captures the effects of the radiation mode well, both at resonance and off-resonance. We have further investigated parameters that can increase or decrease the rate of population transfer, and we find that surface hopping in general can capture both effects well. Finally, we show that the dipole self-energy term within our parameter regime does not significantly affect the system's dynamics.

Generalized Semiclassical Ehrenfest Method: A Route to Wave Function-Free Photochemistry and Nonadiabatic Dynamics with Only Potential Energies and Gradients

We reconsider recent methods by which direct dynamics calculations of electronically nonadiabatic processes can be carried out while requiring only adiabatic potential energies and their gradients. We show that these methods can be understood in terms of a new generalization of the well-known semiclassical Ehrenfest method. This is convenient because it eliminates the need to evaluate electronic wave functions and their matrix elements along the mixed quantum-classical trajectories. The new approximations and procedures enabling this advance are the curvature-driven approximation to the time-derivative coupling, the generalized semiclassical Ehrenfest method, and a new gradient correction scheme called the time-derivative matrix (TDM) scheme. When spin-orbit coupling is present, one can carry out dynamics calculations in the fully adiabatic basis using potential energies and gradients calculated without spin-orbit coupling plus the spin-orbit coupling matrix elements. Even when spin-orbit coupling is neglected, the method is useful because it allows calculations by electronic structure methods for which nonadiabatic coupling vectors are unavailable. In order to place the new considerations in context, the article starts out with a review of background material on trajectory surface hopping, the semiclassical Ehrenfest scheme, and methods for incorporating decoherence. We consider both internal conversion and intersystem crossing. We also review several examples from our group of successful applications of the curvature-driven approximation.

Nonadiabatic Dynamics in a Continuous Circularly Polarized Laser Field with Floquet Phase-Space Surface Hopping

Nonadiabatic chemical reactions involving continuous circularly polarized light (cw CPL) have not attracted as much attention as dynamics in unpolarized/linearly polarized light. However, including circularly (in contrast to linearly) polarized light allows one to effectively introduce a complex-valued time-dependent Hamiltonian, which offers a new path for control or exploration through the introduction of Berry forces. Here, we investigate several inexpensive semiclassical approaches for modeling such nonadiabatic dynamics in the presence of a time-dependent complex-valued Hamiltonian, beginning with a straightforward instantaneous adiabatic fewest-switches surface hopping (IA-FSSH) approach (where the electronic states depend on position and time), continuing to a standard Floquet fewest switches surface hopping (F-FSSH) approach (where the electronic states depend on position and frequency), and ending with an exotic Floquet phase-space surface hopping (F-PSSH) approach (where the electronic states depend on position, frequency, and momentum). Using a set of model systems with time-dependent complex-valued Hamiltonians, we show that the Floquet phase-space adiabats are the optimal choice of basis as far as accounting for Berry phase effects and delivering accuracy. Thus, the F-PSSH algorithm sets the stage for future modeling of nonadiabatic dynamics under strong externally pumped circular polarization.

Decoherence and Its Role in Electronically Nonadiabatic Dynamics

Decoherence is the tendency of a time-evolved reduced density matrix for a subsystem to assume a form corresponding to a statistical ensemble of states rather than a coherent combination of pure-state wave functions. When a molecular process involves changes in the electronic state and the coordinates of the nuclei, as in ultraviolet or visible light photochemistry or electronically inelastic collisions, the reduced density matrix of the electronic subsystem suffers decoherence, due to its interaction with the nuclear subsystem. We present the background necessary to conceptualize this decoherence; in particular, we discuss the density matrix description of pure states and mixed states, and we discuss pointer states and decoherence time. We then discuss how decoherence is treated in the coherent switching with decay of mixing algorithm and the trajectory surface hopping method for semiclassical calculations of electronically nonadiabatic processes.

Semiclassical description of nuclear dynamics moving through complex-valued single avoided crossings of two electronic states

The standard fewest-switches surface hopping (FSSH) approach fails to model nonadiabatic dynamics when the electronic Hamiltonian is complex-valued and there are multiple nuclear dimensions; FSSH does not include geometric magnetic effects and does not have access to a gauge independent direction for momentum rescaling. In this paper, for the case of a Hamiltonian with two electronic states, we propose an extension of Tully's FSSH algorithm, which includes geometric magnetic forces and, through diabatization, establishes a well-defined rescaling direction. When combined with a decoherence correction, our new algorithm shows satisfying results for a model set of two-dimensional single avoided crossings.

A "backtracking" correction for the fewest switches surface hopping algorithm

We propose a "backtracking" mechanism within Tully's fewest switches surface hopping (FSSH) algorithm, whereby whenever one detects consecutive (double) hops during a short period of time, one simply rewinds the dynamics backward in time. In doing so, one reduces the number of hopping events and comes closer to a truly fewest switches surface hopping approach with independent trajectories. With this algorithmic change, we demonstrate that surface hopping can be reasonably accurate for nuclear dynamics in a multidimensional configuration space with a complex-valued (i.e., not real-valued) electronic Hamiltonian; without this adjustment, surface hopping often fails. The added computational cost is marginal. Future research will be needed to assess whether or not this backtracking correction can improve the accuracy of a typical FSSH calculation with a real-valued electronic Hamiltonian (that ignores spin).

Nonadiabatic Dynamics in a Laser Field: Using Floquet Fewest Switches Surface Hopping To Calculate Electronic Populations for Slow Nuclear Velocities

We investigate two well-known approaches for extending the fewest switches surface hopping (FSSH) algorithm to periodic time-dependent couplings. The first formalism acts as if the instantaneous adiabatic electronic states were standard adiabatic states, which just happen to evolve in time. The second formalism replaces the role of the usual adiabatic states by the time-independent adiabatic Floquet states. For a set of modified Tully model problems, the Floquet FSSH (F-FSSH) formalism gives a better estimate for both transmission and reflection probabilities than the instantaneous adiabatic FSSH (IA-FSSH) formalism, especially for slow nuclear velocities. More importantly, only F-FSSH predicts the correct final scattering momentum. Finally, in order to use Floquet theory accurately, we find that it is crucial to account for the interference between wavepackets on different Floquet states. Our results should be of interest to all those interested in laser-induced molecular dynamics.

Nonadiabatic Molecular Dynamics at Metal Surfaces

Dynamics at molecule-metal interfaces are a subject of intense current interest and come in many different flavors of experiments: gas-phase scattering, chemisorption, electrochemistry, nanojunction transport, and heterogeneous catalysis, to name a few. These dynamics involve nuclear degrees of freedom entangled with many electronic degrees of freedom (in the metal), and as such there is always the possibility for nonadiabatic phenomena to appear: the nuclei do not necessarily need to move slower than the electrons to break the Born-Oppenheimer (BO) approximation. In this Feature Article, we review a set of dynamical methods developed recently to deal with such nonadiabatic phenomena at a metal surface, methods that serve as alternatives to Tully's independent electron surface hopping (IESH) model. In the weak molecule-metal coupling regime, a classical master equation (CME) can be derived and a simple surface hopping approach is proposed to propagate nuclear and electronic dynamics stochastically. In the strong molecule-metal interaction regime, a Fokker-Planck equation can be derived for the nuclear dynamics, with electronic DoFs incorporated into the overall friction and random force. Lastly, a broadened classical master equation (BCME) can interpolate between the weak and strong molecule-metal interactions. Here, we briefly review these methods and the relevant benchmarking data, showing in particular how the methods can be used to calculate nonequilibrium transport properties. We highlight several open questions and pose several avenues for future study.

A comparison of surface hopping approaches for capturing metal-molecule electron transfer: A broadened classical master equation versus independent electron surface hopping

Within a generalized Anderson-Holstein model, we investigate electron transfer rates using two different surface hopping algorithms: a broadened classical master equation (BCME) and independent electron surface hopping (IESH). We find that for large enough bandwidth and density of one electron states, and in the presence of external friction, the IESH results converge to the BCME results for impurity-bath model systems, recovering both relaxation rates and equilibrium populations. Without external friction, however, the BCME and IESH results can strongly disagree, and preliminary evidence suggests that IESH does not always recover the correct equilibrium state. Finally, we also demonstrate that adding an electronic thermostat to IESH does help drive the metallic substrate to the correct equilibrium state, but this improvement can sometimes come at the cost of worse short time dynamics. Overall, our results should be of use for all computational chemists looking to model either gas phase scattering or electrochemical dynamics at a metal interface.

An Efficient, Augmented Surface Hopping Algorithm That Includes Decoherence for Use in Large-Scale Simulations

We propose and implement a highly efficient augmented surface hopping algorithm that (i) can be used for large simulations (with many nuclei and many electronic states) and (ii) includes the effects of decoherence without parametrization. Our protocol is based on three key modifications of the surface hopping methodology: (a) a novel separation of classical and quantum degrees of freedom that treats avoided and trivial crossings efficiently, (b) a multidimensional approximation of the time derivative matrix that avoids explicit construction of the derivative coupling at most time steps, and (c) an efficient approximation for the augmented fewest-switches surface hopping decoherence rate. We will show that this protocol can be several orders of magnitude more efficient than the traditional protocol for large multidimensional problems. Furthermore, the marginal cost for including decoherence effects is now negligible.

Development of semiclassical molecular dynamics simulation method

Various quantum mechanical effects such as nonadiabatic transitions, quantum mechanical tunneling and coherence play crucial roles in a variety of chemical and biological systems.

Surface hopping with a manifold of electronic states. I. Incorporating surface-leaking to capture lifetimes

We investigate the incorporation of the surface-leaking (SL) algorithm into Tully's fewest-switches surface hopping (FSSH) algorithm to simulate some electronic relaxation induced by an electronic bath in conjunction with some electronic transitions between discrete states. The resulting SL-FSSH algorithm is benchmarked against exact quantum scattering calculations for three one-dimensional model problems. The results show excellent agreement between SL-FSSH and exact quantum dynamics in the wide band limit, suggesting the potential for a SL-FSSH algorithm. Discrepancies and failures are investigated in detail to understand the factors that will limit the reliability of SL-FSSH, especially the wide band approximation. Considering the easiness of implementation and the low computational cost, we expect this method to be useful in studying processes involving both a continuum of electronic states (where electronic dynamics are probabilistic) and processes involving only a few electronic states (where non-adiabatic processes cannot ignore short-time coherence).

Surface hopping outperforms secular Redfield theory when reorganization energies range from small to moderate (and nuclei are classical)

We evaluate the accuracy of Tully's surface hopping algorithm for the spin-boson model in the limit of small to moderate reorganization energy. We calculate transition rates between diabatic surfaces in the exciton basis and compare against exact results from the hierarchical equations of motion; we also compare against approximate rates from the secular Redfield equation and Ehrenfest dynamics. We show that decoherence-corrected surface hopping performs very well in this regime, agreeing with secular Redfield theory for very weak system-bath coupling and outperforming secular Redfield theory for moderate system-bath coupling. Surface hopping can also be extended beyond the Markovian limits of standard Redfield theory. Given previous work [B. R. Landry and J. E. Subotnik, J. Chem. Phys. 137, 22A513 (2012)] that establishes the accuracy of decoherence-corrected surface-hopping in the Marcus regime, this work suggests that surface hopping may well have a very wide range of applicability.

On surface hopping and time-reversal

While individual Tully-style fewest switches surface hopping (FSSH) trajectories are stochastic and cannot be inverted in time, it is possible to reverse in time the dynamics of a swarm of FSSH trajectories. Here we show exactly how to invert such dynamics, and we investigate the stability of such time-reversed surface hopping dynamics. We demonstrate that FSSH trajectories can be inverted successfully for short time periods, but the time-reversed dynamics become unstable for long times with multiple hopping events. We argue that this instability of FSSH going backward in time can be correlated with the stability of the FSSH algorithm going forward in time.

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.

FIRST PRINCIPLES MOLECULAR DYNAMICS INVOLVING EXCITED STATES AND NONADIABATIC TRANSITIONS

Extensions of traditional molecular dynamics to excited electronic states and non-Born–Oppenheimer dynamics are reviewed focusing on applicability to chemical reactions of large molecules, possibly in condensed phases. The latter imposes restrictions on both the level of accuracy of the underlying electronic structure theory and the treatment of nonadiabaticity. This review, therefore, exclusively deals with ab initio "on the fly" molecular dynamics methods. For the same reason, mainly mixed quantum-classical approaches to nonadiabatic dynamics are considered.

Quantum and semiclassical theories for nonadiabatic transitions based on overlap integrals related to fast degrees of freedom

Alternative treatments of quantum and semiclassical theories for nonadiabatic dynamics are presented. These treatments require no derivative couplings and instead are based on overlap integrals between eigenstates corresponding to fast degrees of freedom, such as electronic states. Derived from mathematical transformations of the Schrödinger equation, the theories describe nonlocal characteristics of nonadiabatic transitions. The idea that overlap integrals can be used for nonadiabatic transitions stems from an article by Johnson and Levine [Chem. Phys. Lett. 13, 168 (1972)]. Furthermore, overlap integrals in path-integral form have been recently made available by Schmidt and Tully [J. Chem. Phys. 127, 094103 (2007)] to analyze nonadiabatic effects in thermal equilibrium systems. The present paper expands this idea to dynamic problems presented in path-integral form that involve nonadiabatic semiclassical propagators. Applications to one-dimensional nonadiabatic transitions have provided excellent results, thereby verifying the procedure. In principle these theories that are presented can be applied to multidimensional systems, although numerical costs could be quite expensive.

A new approach to decoherence and momentum rescaling in the surface hopping algorithm

As originally proposed, the fewest switches surface hopping (FSSH) algorithm does not allow for decoherence between wavefunction amplitudes on different adiabatic surfaces. In this paper, we propose an inexpensive correction to standard FSSH dynamics wherein we explicitly model the decoherence of nuclear wave packets on distinct electronic surfaces. Our augmented fewest switches surface hopping approach is conceptually simple and, thus far, it has allowed us to capture several key features of the exact quantum results. Two points in particular merit attention. First, we obtain the correct branching ratios when a quantum particle passes through more than one region of nonadiabatic coupling. Second, our formalism provides a new and natural approach for rescaling nuclear momenta after a surface hop. Both of these features should become increasingly important as surface hopping schemes are applied to higher-dimensional problems.

Semiclassical description of electronically nonadiabatic dynamics via the initial value representation

The initial value representation (IVR) of semiclassical (SC) theory is used in conjunction with the Meyer-Miller/Stock-Thoss description of electronic degrees of freedom in order to treat electronically nonadiabatic processes. It is emphasized that the classical equations of motion for the nuclear and electronic degrees of freedom that emerge in this description are precisely the Ehrenfest equations of motion (the force on the nuclei is the force averaged over the electronic wave function) but that the trajectories given by these equations of motion do not have the usual shortcomings of the traditional Ehrenfest model when they are used within the SC-IVR framework. For example, in the traditional Ehrenfest model (a mixed quantum-classical approach) the nuclear motion emerges from a nonadiabatic encounter on an average potential energy surface (a weighted average according to the population in the various electronic states), while the SC-IVR describes the correct correlation between electronic and nuclear dynamics, i.e., the nuclear motion is on one potential energy surface or the other depending on the electronic state. Calculations using forward-backward versions of SC-IVR theory are presented to illustrate this behavior. An even more approximate version of the SC-IVR, the linearized approximation (LSC-IVR), is slightly better than the traditional Ehrenfest model, but since it cannot describe quantum coherence effects, the LSC-IVR is also not able to describe the correct correlation between nuclear and electronic dynamics.

Semiclassical theory of electronically nonadiabatic chemical dynamics: incorporation of the Zhu-Nakamura theory into the frozen Gaussian propagation method

The title theory is developed by combining the Herman-Kluk semiclassical theory for adiabatic propagation on single potential-energy surface and the semiclassical Zhu-Nakamura theory for nonadiabatic transition. The formulation with use of natural mathematical principles leads to a quite simple expression for the propagator based on classical trajectories and simple formulas are derived for overall adiabatic and nonadiabatic processes. The theory is applied to electronically nonadiabatic photodissociation processes: a one-dimensional problem of H2+ in a cw (continuous wave) laser field and a two-dimensional model problem of H2O in a cw laser field. The theory is found to work well for the propagation duration of several molecular vibrational periods and wide energy range. Although the formulation is made for the case of laser induced nonadiabatic processes, it is straightforwardly applicable to ordinary electronically nonadiabatic chemical dynamics.

Generalized trajectory surface hopping method based on the Zhu-Nakamura theory

We present a generalized formulation of the trajectory surface hopping method applicable to a general multidimensional system. The method is based on the Zhu-Nakamura theory of a nonadiabatic transition and therefore includes the treatment of classically forbidden hops. The method uses a generalized recipe for the conservation of angular momentum after forbidden hops and an approximation for determining a nonadiabatic transition direction which is crucial when the coupling vector is unavailable. This method also eliminates the need for a rigorous location of the seam surface, thereby ensuring its applicability to a wide class of chemical systems. In a test calculation, we implement the method for the DH(2) (+) system, and it shows a remarkable agreement with the previous results of C. Zhu, H. Kamisaka, and H. Nakamura, [J. Chem. Phys. 116, 3234 (2002)]. We then apply it to a diatomic-in-molecule model system with a conical intersection, and the results compare well with exact quantum calculations. The successful application to the conical intersection system confirms the possibility of directly extending the present method to an arbitrary potential of general topology.

Army ants algorithm for rare event sampling of delocalized nonadiabatic transitions by trajectory surface hopping and the estimation of sampling errors by the bootstrap method

The most widely used algorithm for Monte Carlo sampling of electronic transitions in trajectory surface hopping (TSH) calculations is the so-called anteater algorithm, which is inefficient for sampling low-probability nonadiabatic events. We present a new sampling scheme (called the army ants algorithm) for carrying out TSH calculations that is applicable to systems with any strength of coupling. The army ants algorithm is a form of rare event sampling whose efficiency is controlled by an input parameter. By choosing a suitable value of the input parameter the army ants algorithm can be reduced to the anteater algorithm (which is efficient for strongly coupled cases), and by optimizing the parameter the army ants algorithm may be efficiently applied to systems with low-probability events. To demonstrate the efficiency of the army ants algorithm, we performed atom-diatom scattering calculations on a model system involving weakly coupled electronic states. Fully converged quantum mechanical calculations were performed, and the probabilities for nonadiabatic reaction and nonreactive deexcitation (quenching) were found to be on the order of 10(-8). For such low-probability events the anteater sampling scheme requires a large number of trajectories ( approximately 10(10)) to obtain good statistics and converged semiclassical results. In contrast by using the new army ants algorithm converged results were obtained by running 10(5) trajectories. Furthermore, the results were found to be in excellent agreement with the quantum mechanical results. Sampling errors were estimated using the bootstrap method, which is validated for use with the army ants algorithm.