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.     

Nonadiabatic coupling vectors for excited states within time-dependent density functional theory in the Tamm-Dancoff approximation and beyond

Recently, we have proposed a scheme for the calculation of nonadiabatic couplings and nonadiabatic coupling vectors within linear response time-dependent density functional theory using a set of auxiliary many-electron wavefunctions [I. Tavernelli, E. Tapavicza, and U. Rothlisberger, J. Chem. Phys. 130, 124107 (2009)]. As demonstrated in a later work [I. Tavernelli, B. F. E. Curchod, and U. Rothlisberger, J. Chem. Phys. 131, 196101 (2009)], this approach is rigorous in the case of the calculation of nonadiabatic couplings between the ground state and any excited state. In this work, we extend this formalism to the case of coupling between pairs of singly excited states with the same spin multiplicity. After proving the correctness of our formalism using the electronic oscillator approach by Mukamel and co-workers [S. Tretiak and S. Mukamel, Chem. Rev. (Washington, D.C.) 102, 3171 (2002)], we tested the method on a model system, namely, protonated formaldimine, for which we computed S(1)/S(2) nonadiabatic coupling vectors and compared them with results from high level (MR-CISD) electronic structure calculations.     

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.     

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.     

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.     

Higher order (A+E) multiply sign in circle e pseudo-Jahn-Teller coupling

The pseudo-Jahn-Teller (PJT) coupling of a nondegenerate state A with a twofold degenerate state E by a degenerate vibrational mode e is studied for a general system with a C(3) main rotational axis. The PJT coupling terms up to sixth order are derived by symmetry considerations for this general (A+E) multiply sign in circle e case. The obtained expression for the 3 x 3 diabatic potential energy matrix is found to be closely related to the expression recently developed for the higher order Jahn-Teller case [A. Viel and W. Eisfeld, J. Chem. Phys. 120, 4603 (2004)]. The dynamical PJT coupling, which can arise for states of appropriate symmetry if one of the vibrational modes induces a change of the nuclear point group between D(3h), C(3v), C(3h), and C(3), is discussed. The effect of the higher order PJT coupling is tested by a two-dimensional model study based on the e bending mode of NH(3)(+). The models are analyzed by fitting the two-dimensional potential energy surfaces. The significance of the higher order terms on the nonadiabatic dynamics is demonstrated by quantum wave packet propagations.     

The role of the excited electronic states in the C+ + H2O reaction

The electronic excited states of the [COH2]+ system have been studied in order to establish their role in the dynamics of the C+ + H2O-->[COH]+ +H reaction, which is a prototypical ion-molecule reaction. The most relevant minima and saddle points of the lowest excited state have been determined and energy profiles for the lowest excited doublet and quartet electronic states have been computed along the fragmentation and isomerization coordinates. Also, nonadiabatic coupling strengths between the ground and the first excited state have been computed where they can be large. Our analysis suggests that the first excited state could play an important role in the generation of the formyl isomer, which has been detected in crossed beam experiments [D. M. Sonnenfroh et al., J. Chem. Phys. 83, 3985 (1985)], but could not be explained in quasiclassical trajectory computations [Y. Ishikawa et al., Chem. Phys. Lett. 370, 490 (2003); J. R. Flores, J. Chem. Phys. 125, 164309 (2006)].     

Torsion-vibration coupling in methanol: the adiabatic approximation and intramolecular vibrational redistribution scaling

The four-dimensional model Hamiltonian of Wang and Perry [J. Chem. Phys. 109, 10795 (1998)] is used to compare the approximate adiabatic separation of the torsion and CH stretches in methanol to an exact solution of the same Hamiltonian. The adiabatic approximation accounts for the pattern of the energy levels in the lowest torsional states, including the inverted tunneling splittings, but does not account for the pattern of systematic two- and four-fold near degeneracies at high torsional excitation. In the adiabatic basis, the nonadiabatic couplings mix the torsional and vibrational degrees of freedom and hence are a source for intramolecular vibrational redistribution (IVR). These IVR matrix elements are found to decrease by only a factor of 2 or 3 with each higher coupling order, in agreement with the results of Pearman and Gruebele [Z. Phys. Chem. Munich 214, 1439 (2000)]. This gentle scaling behavior, which contrasts with a steeper falloff with coupling order in more rigid molecules, points to a more important role for direct high-order couplings in torsional molecules. In this model, the scaling behavior derives from a single coupling term that is low order in the torsional angular momentum in combination with one-dimensional torsional functions that include contributions from many torsional angular momenta.     

Combined electronic and nuclear dynamics in a simple model system. II. Spectroscopic transitions

We investigate the correlated electronic and nuclear motion in a model system as proposed by Shin and Metiu [J. Chem. Phys. 102, 9285 (1995)]. The quantum dynamics is studied during laser induced electronic transitions. Here, the influence of nonadiabatic coupling on the absorption spectrum is investigated and the Franck-Condon principle is illustrated in terms of the temporal changes of electronic and nuclear densities. In the case of intense field excitation, multiphoton processes become important, and electronic as well as vibrational wave packets are prepared.     

Nonadiabatic quantum dynamics in O(3P)+H2→OH+H: a revisited study

To investigate the extent of nonadiabatic effects in the title reaction, quasi-classical trajectory and nonadiabatic quantum scattering as well as the nonadiabatic quantum-classical trajectory calculations were performed on the accurate ab initio benchmark potential energy surfaces of the lowest (3)A' and (3)A" electronic states [Rogers et al., J Phys Chem A 2000, 104, 2308], together with the spin-orbit coupling matrix [Maiti and Schatz, J Chem Phys 2003, 119, 12360] and the lowest singlet (1) A' potential energy surface [Dobby and Knowles, Faraday Discuss 1998, 110, 247]. Comparison of the calculated total cross sections from both adiabatic and nonadiabatic calculations has demonstrated that for adiabatic channels including (3)A'→(3)A' and (3)A"→(3)A", difference does exist between the two kinds of adiabatic and nonadiabatic calculations, showing nonadiabatic effects to some extent. Such nonadiabatic effects tend to become more conspicuous at high collision energies and are found to be more pronounced with trajectories/quantum wave packet initiated on (3)A' than on (3)A". Furthermore, the present study also showed that nonadiabatic effects can bring the component of forward-scattering in the product angular distributions.     

The electronic nonadiabatic coupling term: can it be ignored in dynamic calculations?

Whereas the search for the degeneracy points which are better known as conical intersections (or ci-points) is usually carried out with a lot of devotion, the nonadiabatic coupling terms (NACTs) which together with the adiabatic potential energy surfaces appear in the nuclear Born-Oppenheimer-Schrodinger equation are ignored in most dynamical calculations. In the present article we consider two well known frameworks, namely, the semiclassical surface hopping method and the vibrational coupling model Hamiltonian that avoid the NACTs and examine to what extent, this procedure is justified.     

Single switch surface hopping for a model of pyrazine

The single switch trajectory surface hopping algorithm is tested for numerical simulations of a two-state three-mode model for the internal conversion of pyrazine through a conical intersection of potential energy surfaces. The algorithm is compared to two other surface hopping approaches, namely, Tully's method of the fewest switches [J. Chem. Phys. 93, 1061 (1990)] and the method by Voronin et al. [J. Phys. Chem. A 102, 6057 (1998)]. The single switch algorithm achieves the most accurate results. Replacing its deterministic nonadiabatic branching condition by a probabilistic accept-reject criterion, one obtains the method of Voronin et al. without momentum adjustment. This probabilistic version of the single switch approach outperforms the considered algorithms in terms of accuracy, memory requirement, and runtime.     

Space-time contours to treat the interaction between an intense electric field and a molecular system

A few years ago, we developed an approach to treat molecular systems exposed to an external, intense, time-dependent field (J. Phys. Chem. A 2003, 107, 4724; J. Chem. Phys. 2003, 119, 6998). Within this study, we encountered two novel concepts: the dressed (namely, field affected) time-dependent nonadiabatic coupling term and the space-time contours. In the present article, we analyze the newly introduced nonadiabatic coupling term and discuss its importance for dynamical studies. We also refer to the just mentioned space-time contour and present the more efficient contour for realistic situations. The scope of the above-mentioned articles is extended with the aim of defining quasi-adiabatic states for such situations. Strictly speaking, molecular systems exposed to intense, fast oscillating fields are not expected to form adiabatic states. Still we consider such a situation and end up with three possibilities for quasi-adiabatical time-dependent states eventually to be used within the Born-Oppenheimer approximation.     

Semiclassical treatment of thermally activated electron transfer in the inverted region under the fast dielectric relaxation

The previously formulated semiclassical theory (Zhao, Liang, and Nakamura, J. Phys. Chem. A 2006, 110, 8204) is used to study electron transfer in the Marcus inverted case by considering multidimensional potential energy surfaces of donor and acceptor. The Zhu-Nakamura formulas of nonadiabatic transition in the case of Landau-Zener type are incorporated into the approach. The theory properly takes into account the nonadiabatic transition coupled with the nuclear tunneling and can cover the whole range from weak to strong coupling regime uniformly under the assumption of fast solvent relaxation. The numerical calculations are performed for the 12-dimensional model of shifted harmonic oscillators and demonstrate that the reaction rate with respect to the electronic coupling shows a maximum, confirming the adiabatic suppression in the strong coupling limit. The adiabatic suppression is dramatically reduced by the effect of nuclear tunneling compared to the case that the Landau-Zener formula is used. The possible extension and applications to the case of the slow solvent dynamics are discussed.     

Ab initio potential energy surfaces and nonadiabatic interactions in the H+ +NO collision system

Ab initio calculations on the H(+)+NO system have been carried out in Jacobi coordinates at the multireference configuration interaction level employing Dunning's correlation-consistent polarized valence triple zeta basis set to analyze the role of low-lying electronic excited states in influencing the collision dynamics relevant to the experimental collision energy range of 9.5-30 eV. The lowest two adiabatic potential energy surfaces, asymptotically correlating to H(+)+NO(X (2)Pi) and H((2)S)+NO(+)(X (1)Sigma(+)), have been obtained. Using ab initio procedures, the (radial) nonadiabatic couplings and the mixing angle between the lowest two electronic states (1 (2)A' and 2 (2)A') have been obtained to yield the corresponding quasidiabatic potential energy matrix. The strengths of the computed vibrational coupling matrix elements reflect a similar trend, as has been observed experimentally in the magnitudes of the state-to-state transition probability for the inelastic vibrational excitations [J. Krutein and F. Linder, J. Chem. Phys. 71, 559 (1979); F. A. Gianturco et al., J. Phys. B 14, 667 (1981)].     

A field theoretical approach to calculate electronic Born-Oppenheimer coupling terms

In this paper we suggest to consider the spatial distribution of the Born-Oppenheimer nonadiabatic coupling terms as fields which are created by sources, located at degeneracy points, and which can be derived using the ordinary mathematical tools of field theory. It is shown that the curl-divergence equations as formed within a given Hilbert space [M. Baer, Chem. Phys. Lett. 35, 112 (1975)] can be converted into a set of inhomogeneous coupled Poisson equations which are solved for a given set of boundary conditions. The method is applied to the three-state Hilbert subspace of the H(3) system. The numerical results are compared with ab initio calculations for which a very encouraging fit is found.     

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.     

Renner-Teller nonadiabatic coupling terms: An ab-initio study of the HNH molecule

In this Communication we present the first theoretical/numerical treatment of nonadiabatic coupling terms (NACT) that originate from the Renner-Teller (RT) model, namely, those that follow from the splitting of an electronic level of a linear molecule when it becomes bent. These two newly formed states are characterized by different symmetries and are designated as A and B. Our main findings: (1) The RT NACTs are quantized as long as they are calculated close enough to collinear configuration of the molecule (in this case HNH). Their value is tau = 1 (the Jahn-Teller values in similar situations, are tau = (1/2)). (2) Calculation of RT NACTs at bent configurations (i.e., at a distance from the linear axis) yield decreased values, sometimes by more than 50%. This last finding implies that in strongly bent configurations the two-state Hilbert subspace (formed by the above mentioned A and B states) is affected by upper states, most likely via Jahn-Teller conical intersections. (3) This study has also important practical implications. The fact that the RT NACTs decrease in (strongly) bent situations implies that analyzing spectroscopic data employing only the two Pi-states may not be sufficient in order to achieve the required accuracy.     

The E x e dynamic Jahn-Teller problem: a new insight from the strong coupling limit

Correct boundary conditions for the E x e dynamic Jahn-Teller problem are considered explicitly for the first time to obtain approximate analytical solutions in the strong coupling limit. Numerical solutions for the decoupled equations using the finite difference method are also presented. The numerical solutions for the decoupled equations exhibit avoided crossings in the weak coupling region, which explains the oscillating behavior of the solutions obtained by Longuet-Higgins et al. for the coupled equations. The obtained analytical energy expressions show improved agreement with the numerical calculations as compared with the previous treatment in which the potentials were assumed to be harmonic. We demonstrate that the pseudorotational energy j(2)/(2g(2)), where g is the dimensionless vibronic coupling constant, and j total angular momentum: j=+/-1/2,+/-3/2,..., in the conventional strong coupling expression for the vibronic levels of the lower sheet is exact. Non-Hermitian first-order perturbation theory gives the energy which is correct up to 1/g(4). The asymptotic behavior of the wave function at the origin does not influence the corrected energy up to order of 1/g(4). At the same time the treatment of the upper sheet with correct boundary conditions gives solutions which are entirely different from the corresponding Slonczewski's solutions. Besides, the correct boundary conditions enable us to evaluate the nonadiabatic coupling between the lower and upper potential sheets. The energy correction due to the nonadiabatic coupling is estimated to be of order 1/g(6).     

Recovering the Crooks equation for dynamical systems in the isothermal-isobaric ensemble: a strategy based on the equations of motion

The Crooks equation [Eq. (10) in J. Stat. Phys. 90, 1481 (1998)] relates the work done on a system during a nonequilibrium transformation to the free energy difference between the final and the initial state of the transformation. Recently, the authors have derived the Crooks equation for systems in the canonical ensemble thermostatted by the Nose-Hoover or Nose-Hoover chain method [P. Procacci et al., J. Chem. Phys. 125, 164101 (2006)]. That proof is essentially based on the fluctuation theorem by Evans and Searles [Adv. Phys. 51, 1529 (2002)] and on the equations of motion. Following an analogous approach, the authors derive here the Crooks equation in the context of molecular dynamics simulations of systems in the isothermal-isobaric (NPT) ensemble, whose dynamics is regulated by the Martyna-Tobias-Klein algorithm [J. Chem. Phys. 101, 4177 (1994)]. Their present derivation of the Crooks equation correlates to the demonstration of the Jarzynski identity for NPT systems recently proposed by Cuendet [J. Chem. Phys. 125, 144109 (2006)].