Nonadiabatic surface hopping Herman-Kluk semiclassical initial value representation method revisited: applications to Tully's three model systems

The nonadiabatic surface hopping Herman-Kluk (HK) semiclassical initial value representation (SC-IVR) method for nonadiabatic problems is reformulated. The method has the same spirit as Tully's surface hopping technique [J. Chem. Phys. 93, 1061 (1990)] and almost keeps the same structure as the original single-surface HK SC-IVR method except that trajectories can hop to other surfaces according to the hopping probabilities and phases, which can be easily integrated along the paths. The method is based on a rather general nonadiabatic semiclassical surface hopping theory developed by Herman [J. Chem. Phys. 103, 8081 (1995)], which has been shown to be accurate to the first order in h and through all the orders of the nonadiabatic coupling amplitude. Our simulation studies on the three model systems suggested by Tully demonstrate that this method is practical and capable of describing nonadiabatic quantum dynamics for various coupling situations in very good agreement with benchmark calculations.     

Numerical study of the accuracy and efficiency of various approaches for Monte Carlo surface hopping calculations

A one-dimensional, two-state model problem with two well-separated avoided crossing points is employed to test the efficiency and accuracy of a semiclassical surface hopping technique. The use of a one-dimensional model allows for the accurate numerical evaluation of both fully quantum-mechanical and semiclassical transition probabilities. The calculations demonstrate that the surface hopping procedure employed accounts for the interference between different hopping trajectories very well and provides highly accurate transition probabilities. It is, in general, not computationally feasible to completely sum over all hopping trajectories in the semiclassical calculations for multidimensional problems. In this case, a Monte Carlo procedure for selecting important trajectories can be employed. However, the cancellation due to the different phases associated with different trajectories limits the accuracy and efficiency of the Monte Carlo procedure. Various approaches for improving the accuracy and efficiency of Monte Carlo surface hopping procedures are investigated. These methods are found to significantly reduce the statistical sampling errors in the calculations, thereby increasing the accuracy of the transition probabilities obtained with a fixed number of trajectories sampled.     

A justification for a nonadiabatic surface hopping Herman-Kluk semiclassical initial value representation of the time evolution operator

A justification is given for the validity of a nonadiabatic surface hopping Herman-Kluk (HK) semiclassical initial value representation (SC-IVR) method. The method is based on a propagator that combines the single surface HK SC-IVR method [J. Chem. Phys. 84, 326 (1986)] and Herman's nonadiabatic semiclassical surface hopping theory [J. Chem. Phys. 103, 8081 (1995)], which was originally developed using the primitive semiclassical Van Vleck propagator. We show that the nonadiabatic HK SC-IVR propagator satisfies the time-dependent Schrodinger equation to the first order of variant Planck's over 2pi and the error is O(variant Planck's over 2pi(2)). As a required lemma, we show that the stationary phase approximation, under current assumptions, has an error term variant Planck's over 2pi(1) order higher than the leading term. Our derivation suggests some changes to the previous development, and it is shown that the numerical accuracy in applications to Tully's three model systems in low energies is improved.     

Exploring the role of decoherence in condensed-phase nonadiabatic dynamics: a comparison of different mixed quantum/classical simulation algorithms for the excited hydrated electron

Mixed quantum/classical (MQC) molecular dynamics simulation has become the method of choice for simulating the dynamics of quantum mechanical objects that interact with condensed-phase systems. There are many MQC algorithms available, however, and in cases where nonadiabatic coupling is important, different algorithms may lead to different results. Thus, it has been difficult to reach definitive conclusions about relaxation dynamics using nonadiabatic MQC methods because one is never certain whether any given algorithm includes enough of the necessary physics. In this paper, we explore the physics underlying different nonadiabatic MQC algorithms by comparing and contrasting the excited-state relaxation dynamics of the prototypical condensed-phase MQC system, the hydrated electron, calculated using different algorithms, including: fewest-switches surface hopping, stationary-phase surface hopping, and mean-field dynamics with surface hopping. We also describe in detail how a new nonadiabatic algorithm, mean-field dynamics with stochastic decoherence (MF-SD), is to be implemented for condensed-phase problems, and we apply MF-SD to the excited-state relaxation of the hydrated electron. Our discussion emphasizes the different ways quantum decoherence is treated in each algorithm and the resulting implications for hydrated-electron relaxation dynamics. We find that for three MQC methods that use Tully's fewest-switches criterion to determine surface hopping probabilities, the excited-state lifetime of the electron is the same. Moreover, the nonequilibrium solvent response function of the excited hydrated electron is the same with all of the nonadiabatic MQC algorithms discussed here, so that all of the algorithms would produce similar agreement with experiment. Despite the identical solvent response predicted by each MQC algorithm, we find that MF-SD allows much more mixing of multiple basis states into the quantum wave function than do other methods. This leads to an excited-state lifetime that is longer with MF-SD than with any method that incorporates nonadiabatic effects with the fewest-switches surface hopping criterion.     

Using semiclassical surface hopping for coupled partial wave calculations on systems with non-spherically symmetric potentials

It is shown that a semiclassical surface hopping (SH) approach provides a simple and efficient method for scattering calculations with non-spherically symmetric potentials. The calculations are performed by expanding the wave function in an angular momentum state basis. Since the potential is not spherically symmetric, the different angular states are coupled. The semiclassical SH method, which is typically used for problems with coupled electronic states, can, in principle, be employed for any coupled state problem. The particular SH method employed is known to provide highly accurate results for coupled electronic state problems. The method is tested on model two angular state problems using potential surfaces and couplings arising from a non-spherically symmetric scattering problem. The results for these model problems are in excellent agreement with exact quantum calculations. Full calculations, which are converged with regard to the number of angular basis states, are also performed for the non-spherically symmetric problem. It is shown that an approximation to the surface hopping amplitudes that simplifies the numerical implementation of the method provides results in excellent agreement with the full surface hopping calculation.     

Measuring nonadiabaticity of molecular quantum dynamics with quantum fidelity and with its efficient semiclassical approximation

We propose to measure nonadiabaticity of molecular quantum dynamics rigorously with the quantum fidelity between the Born-Oppenheimer and fully nonadiabatic dynamics. It is shown that this measure of nonadiabaticity applies in situations where other criteria, such as the energy gap criterion or the extent of population transfer, fail. We further propose to estimate this quantum fidelity efficiently with a generalization of the dephasing representation to multiple surfaces. Two variants of the multiple-surface dephasing representation (MSDR) are introduced, in which the nuclei are propagated either with the fewest-switches surface hopping or with the locally mean field dynamics (LMFD). The LMFD can be interpreted as the Ehrenfest dynamics of an ensemble of nuclear trajectories, and has been used previously in the nonadiabatic semiclassical initial value representation. In addition to propagating an ensemble of classical trajectories, the MSDR requires evaluating nonadiabatic couplings and solving the Schro?dinger (or more generally, the quantum Liouville-von Neumann) equation for a single discrete degree of freedom. The MSDR can be also used in the diabatic basis to measure the importance of the diabatic couplings. The method is tested on three model problems introduced by Tully and on a two-surface model of dissociation of NaI.     

On the properties of a primitive semiclassical surface hopping propagator for nonadiabatic quantum dynamics

A previously developed nonadiabatic semiclassical surface hopping propagator [M. F. Herman J. Chem. Phys. 103, 8081 (1995)] is further studied. The propagator has been shown to satisfy the time-dependent Schrodinger equation (TDSE) through order h, and the O(h2) terms are treated as small errors, consistent with standard semiclassical analysis. Energy is conserved at each hopping point and the change in momentum accompanying each hop is parallel to the direction of the nonadiabatic coupling vector resulting in both transmission and reflection types of hops. Quantum mechanical analysis and numerical calculations presented in this paper show that the h2 terms involving the interstate coupling functions have significant effects on the quantum transition probabilities. Motivated by these data, the h2 terms are analyzed for the nonadiabatic semiclassical propagator. It is shown that the propagator can satisfy the TDSE for multidimensional systems by including another type of nonclassical trajectories that reflect on the same surfaces. This h2 analysis gives three conditions for these three types of trajectories so that their coefficients are uniquely determined. Besides the nonadiabatic semiclassical propagator, a numerically useful quantum propagator in the adiabatic representation is developed to describe nonadiabatic transitions.     

Decoherence and surface hopping: when can averaging over initial conditions help capture the effects of wave packet separation?

Fewest-switches surface hopping (FSSH) is a popular nonadiabatic dynamics method which treats nuclei with classical mechanics and electrons with quantum mechanics. In order to simulate the motion of a wave packet as accurately as possible, standard FSSH requires a stochastic sampling of the trajectories over a distribution of initial conditions corresponding, e.g., to the Wigner distribution of the initial quantum wave packet. Although it is well-known that FSSH does not properly account for decoherence effects, there is some confusion in the literature about whether or not this averaging over a distribution of initial conditions can approximate some of the effects of decoherence. In this paper, we not only show that averaging over initial conditions does not generally account for decoherence, but also why it fails to do so. We also show how an apparent improvement in accuracy can be obtained for a fortuitous choice of model problems, even though this improvement is not possible, in general. For a basic set of one-dimensional and two-dimensional examples, we find significantly improved results using our recently introduced augmented FSSH algorithm.     

Globally uniform semiclassical surface-hopping wave function for nonadiabatic scattering

A globally uniform time-independent semiclassical wave function for nonadiabatic scattering is presented. This wave function, which takes the form of a surface-hopping expansion, is motivated by the globally uniform semiclassical wave function of Kay and co-workers for the single-surface case. The surface-hopping expansion is similar to a previously presented primitive semiclassical wave function for nonadiabatic problems. This earlier wave function has the important feature that it correctly incorporates all phase terms, allowing for an accurate treatment of quantum interference effects. The globally uniform expression has important numerical advantages over the primitive formulation. The globally uniform wave function does not have caustic singularities, and the globally uniform calculation avoids a root search for trajectories obeying double-ended boundary conditions that is required by the primitive semiclassical calculation.     

Quantum dynamics study of H+NH3-->H2+NH2 reaction

We report in this paper a quantum dynamics study for the reaction H+NH3-->NH2+H2 on the potential energy surface of Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The quantum dynamics calculation employs the semirigid vibrating rotor target model [J. Z. H. Zhang, J. Chem. Phys. 111, 3929 (1999)] and time-dependent wave packet method to propagate the wave function. Initial state-specific reaction probabilities are obtained, and an energy correction scheme is employed to account for zero point energy changes for the neglected degrees of freedom in the dynamics treatment. Tunneling effect is observed in the energy dependency of reaction probability, similar to those found in H+CH4 reaction. The influence of rovibrational excitation on reaction probability and stereodynamical effect are investigated. Reaction rate constants from the initial ground state are calculated and are compared to those from the transition state theory and experimental measurement.     

Semiclassical nonadiabatic dynamics based on quantum trajectories for the O(3P,1D) + H2 system

The O(3P,1D) + H2 --> OH + H reaction is studied using trajectory dynamics within the approximate quantum potential approach. Calculations of the wave-packet reaction probabilities are performed for four coupled electronic states for total angular momentum J = 0 using a mixed coordinate/polar representation of the wave function. Semiclassical dynamics is based on a single set of trajectories evolving on an effective potential-energy surface and in the presence of the approximate quantum potential. Population functions associated with each trajectory are computed for each electronic state. The effective surface is a linear combination of the electronic states with the contributions of individual components defined by their time-dependent average populations. The wave-packet reaction probabilities are in good agreement with the quantum-mechanical results. Intersystem crossing is found to have negligible effect on reaction probabilities summed over final electronic states.     

Surface temperature effect on the scattering of D2(v = 0, j = 0)-Cu(111) system

We perform four-dimensional (4D?2D) as well as six-dimensional (6D) quantum dynamics on a parametrically time- and temperature-dependent effective Hamiltonian for D(2)(v, j)-Cu(111) system, where such effective potential has been derived through a mean-field approach between molecular degrees of freedom and surface modes with Bose-Einstein probability factor for their initial state distribution. We present the convergence of the theoretically calculated sticking probabilities employing 4D?2D quantum dynamics with increasing number of surface atoms as well as layers for rigid surface and the surface at a particular temperature, where the temperature-dependent sticking probabilities appear exclusively dictated by those surface modes directed along the Z-axis. The sticking and state-to-state transition probabilities obtained from 6D quantum dynamics are shown as a function of initial kinetic energy of the diatom at different surface temperature. Theoretically calculated sticking probabilities display the similar trend with the experimentally measured one.     

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.     

Reactive scattering in hyperspherical coordinates

Preliminary results are presented for transition probabilities in the H + H₂ system derived from an adiabatic representation in terms of surface functions on hyperspheres. Special attention is given to the resonance structure for transition probabilities in the first vibrational level.     

Conical intersections and semiclassical trajectories: comparison to accurate quantum dynamics and analyses of the trajectories

Semiclassical trajectory methods are tested for electronically nonadiabatic systems with conical intersections. Five triatomic model systems are presented, and each system features two electronic states that intersect via a seam of conical intersections (CIs). Fully converged, full-dimensional quantum mechanical scattering calculations are carried out for all five systems at energies that allow for electronic de-excitation via the seam of CIs. Several semiclassical trajectory methods are tested against the accurate quantum mechanical results. For four of the five model systems, the diabatic representation is the preferred (most accurate) representation for semiclassical trajectories, as correctly predicted by the Calaveras County criterion. Four surface hopping methods are tested and have overall relative errors of 40%-60%. The semiclassical Ehrenfest method has an overall error of 66%, and the self-consistent decay of mixing (SCDM) and coherent switches with decay of mixing (CSDM) methods are the most accurate methods overall with relative errors of approximately 32%. Furthermore, the CSDM method is less representation dependent than both the SCDM and the surface hopping methods, making it the preferred semiclassical trajectory method. Finally, the behavior of semiclassical trajectories near conical intersections is discussed.     

处理量子反应散射的负虚数势方法

处理量子反应散射的负虚数势方法蔡宇民，胡旭光，李前树（西安石油学院化学工程系．西安，７１００６１）（吉林大学理论化学研究所，长春，１３００２３）关键词反应散射，反应几率，光学势近年来，处理反应散射问题的各种严格的和近似的量子力学方法发展迅速［１～５］...     

A line integral reaction path approximation for large systems via nonlinear constrained optimization: application to alanine dipeptide and the beta hairpin of protein G

A variation of the line integral method of Elber with self-avoiding walk has been implemented using a state of the art nonlinear constrained optimization procedure. The new implementation appears to be robust in finding approximate reaction paths for small and large systems. Exact transition states and intermediates for the resulting paths can easily be pinpointed with subsequent application of the conjugate peak refinement method [S. Fischer and M. Karplus, Chem. Phys. Lett. 194, 252 (1992)] and unconstrained minimization, respectively. Unlike previous implementations utilizing a penalty function approach, the present implementation generates an exact solution of the underlying problem. Most importantly, this formulation does not require an initial guess for the path, which makes it particularly useful for studying complex molecular rearrangements. The method has been applied to conformational rearrangements of the alanine dipeptide in the gas phase and in water, and folding of the beta hairpin of protein G in water. In the latter case a procedure was developed to systematically sample the potential energy surface underlying folding and reconstruct folding pathways within the nearest-neighbor hopping approximation.     

Quantum scattering calculation of the O(1D)+HBr reaction

Three-dimensional time-dependent quantum wave packet calculation for the O((1)D)+HBr reaction has been carried out using an accurate ab initio global potential energy surface [K. A. Peterson, J. Chem. Phys. 113, 4598 (2000)]. The calculations show that the initial state-selected reaction probabilities are dominated by resonance structures, and the lifetime of the resonance is generally in the subpicosecond time scale. The energy dependence of the reaction cross section is computed, which manifests still resonance structures, and is a decreasing function of the translational energy. The thermal rate constants are also computed, which are nearly independent on the temperature. The calculation results are discussed and compared to similar reaction with deep well.