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.     

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.     

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.     

Wavepacket propagation using time-sliced semiclassical initial value methods

A new semiclassical initial value representation (SC-IVR) propagator and a SC-IVR propagator originally introduced by Kay [J. Chem. Phys. 100, 4432 (1994)], are investigated for use in the split-operator method for solving the time-dependent Schrodinger equation. It is shown that the SC-IVR propagators can be derived from a procedure involving modified Filinov filtering of the Van Vleck expression for the semiclassical propagator. The two SC-IVR propagators have been selected for investigation because they avoid the need to perform a coherent state basis set expansion that is necessary in other time-slicing propagation schemes. An efficient scheme for solving the propagators is introduced and can be considered to be a semiclassical form of the effective propagators of Makri [Chem. Phys. Lett. 159, 489 (1989)]. Results from applications to a one-dimensional, two-dimensional, and three-dimensional Hamiltonian for a double-well potential are presented.     

Semiclassical calculation of nonadiabatic thermal rate constants: application to condensed phase reactions

The nonadiabatic transition state theory proposed recently by Zhao et al. [J. Chem. Phys. 121, 8854 (2004)] is extended to calculate rate constants of complex systems by using the Monte Carlo and umbrella sampling methods. Surface hopping molecular dynamics technique is incorporated to take into account the dynamic recrossing effect. A nontrivial benchmark model of the nonadiabatic reaction in the condensed phase is used for the numerical test. It is found that our semiclassical results agree well with those produced by the rigorous quantum mechanical method. Comparing with available analytical approaches, we find that the simple statistical theory proposed by Straub and Berne [J. Chem. Phys. 87, 6111 (1987)] is applicable for a wide friction region although their formula is obtained using Landau-Zener [Phys. Z. Sowjetunion 2, 46 (1932); Proc. R. Soc. London, Ser. A 137, 696 (1932)] nonadiabatic transition probability along a one-dimensional diffusive coordinate. We also investigate how the nuclear tunneling events affect the dependence of the rate constant on the friction.     

Reconciling semiclassical and Bohmian mechanics. V. Wavepacket dynamics

In previous articles [B. Poirier J. Chem. Phys. 121, 4501 (2004); C. Trahan and B. Poirier, ibid. 124, 034115 (2006); 124, 034116 (2006); B. Poirier and G. Parlant, J. Phys. Chem. A 111, 10400 (2007)] a bipolar counterpropagating wave decomposition, psi = psi(+) + psi(-), was presented for stationary states psi of the one-dimensional Schrodinger equation, such that the components psi(+/-) approach their semiclassical Wentzel-Kramers-Brillouin analogs in the large action limit. The corresponding bipolar quantum trajectories are classical-like and well behaved, even when psi has many nodes, or is wildly oscillatory. In this paper, the method is generalized for time-dependent wavepacket dynamics applications and applied to several benchmark problems, including multisurface systems with nonadiabatic coupling.     

Unified semiclassical theory for the two-state system: an analytical solution for general nonadiabatic tunneling

Unified semiclasical solution for general nonadiabatic tunneling between two adiabatic potential energy surfaces is established by employing unified semiclassical solution for pure nonadiabatic transition [C. Zhu, J. Chem. Phys. 105, 4159 (1996)] with the certain symmetry transformation. This symmetry comes from a detailed analysis of the reduced scattering matrix for Landau-Zener type of crossing as a special case of nonadiabatic transition and nonadiabatic tunneling. Traditional classification of crossing and noncrossing types of nonadiabatic transition can be quantitatively defined by the rotation angle of adiabatic-to-diabatic transformation, and this rotational angle enters the analytical solution for general nonadiabatic tunneling. The certain two-state exponential potential models are employed for numerical tests, and the calculations from the present general nonadiabatic tunneling formula are demonstrated in very good agreement with the results from exact quantum mechanical calculations. The present general nonadiabatic tunneling formula can be incorporated with various mixed quantum-classical methods for modeling electronically nonadiabatic processes in photochemistry.     

A new approximation for atom-diatom rotational-relaxation cross sections

A semiclassical approximation to the S matrix of the infinite-order-sudden approximation is introduced. This is employed to yield for the energy-transfer effective cross section a purely classical approximation, analogous to the Mason-Monchick approximation [J. Chem. Phys. 36, 1622 (1962)] for traditional collision integrals. Constraints on energy and on angular momentum transfer are included. Numerical evaluation of this new approximation can readily be performed alongside that for traditional collision integrals. The new result is tested against full classical trajectory calculations for six potential energy surfaces for the collision systems H-N(2), He-N(2), He-CO, and Ar-CO(2). Differences of no more than 15% from the classical trajectory calculations have been obtained.     

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.     

A semiclassical hybrid approach to many particle quantum dynamics

We analytically derive a correlated approach for a mixed semiclassical many particle dynamics, treating a fraction of the degrees of freedom by the multitrajectory semiclassical initial value method of Herman and Kluk [Chem. Phys. 91, 27 (1984)] while approximately treating the dynamics of the remaining degrees of freedom with fixed initial phase space variables, analogously to the thawed Gaussian wave packet dynamics of Heller [J. Chem. Phys. 62, 1544 (1975)]. A first application of this hybrid approach to the well studied Secrest-Johnson [J. Chem. Phys. 45, 4556 (1966)] model of atom-diatomic collisions is promising. Results close to the quantum ones for correlation functions as well as scattering probabilities could be gained with considerably reduced numerical effort as compared to the full semiclassical Herman-Kluk approach. Furthermore, the harmonic nature of the different degrees of freedom can be determined a posteriori by comparing results with and without the additional approximation.     

Semiclassical initial value representation treatment of a hydrogen bonded complex of rigid water molecules from a single trajectory in Cartesian coordinates

A semiclassical initial value representation approach for molecular systems in Cartesian coordinates is combined with a recently proposed time averaging technique [J. Chem. Phys. 118, 7174 (2003)]. It is shown that a single trajectory can yield the zero-point energy of the water dimer with good accuracy for the model chosen when compared to fully constrained Cartesian semiclassical calculations. The convergence with respect to the number of averaging time origins is discussed.     

Molecular symmetry and ab initio calculations: IV. Symmetry-matrix and symmetry-supermatrix in calculations of two-electron repulsion integrals

The product of two Gaussians having different centers is itself a one-center Gaussian, thus multicenter integrals with a Cartesian Gaussian basis can be reduced to one-center integrals. Recurrence relations for overlap integrals and electron repulsion integrals (ERIs) are derived at these centers. The calculations of overlap integrals and ERIs are carried out step by step from the highest symmetry case (one center) to required cases (different centers) by using the translation of Cartesian Gaussians. Full exploitation of symmetry in calculation processes can result in optimal use of these recurrence relations. Compared with the recently published algorithms, based on the recurrence relations derived by Obara and Saika [J. Chem. Phys., 84 , 3963 (1986)], the floating point operations (FLOPs) for ERI calculations (having four different centers) can be reduced by a factor of ca. 2. A significant extra saving in calculations and storage can be obtained if atoms, linear, or planar molecules are discussed. © 1997 John Wiley & Sons, Inc.     

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.     

Generalization of classical mechanics for nuclear motions on nonadiabatically coupled potential energy surfaces in chemical reactions

Classical trajectory study of nuclear motion on the Born-Oppenheimer potential energy surfaces is now one of the standard methods of chemical dynamics. In particular, this approach is inevitable in the studies of large molecular systems. However, as soon as more than a single potential energy surface is involved due to nonadiabatic coupling, such a naive application of classical mechanics loses its theoretical foundation. This is a classic and fundamental issue in the foundation of chemistry. To cope with this problem, we propose a generalization of classical mechanics that provides a path even in cases where multiple potential energy surfaces are involved in a single event and the Born-Oppenheimer approximation breaks down. This generalization is made by diagonalization of the matrix representation of nuclear forces in nonadiabatic dynamics, which is derived from a mixed quantum-classical representation of the electron-nucleus entangled Hamiltonian [Takatsuka, K. J. Chem. Phys. 2006, 124, 064111]. A manifestation of quantum fluctuation on a classical subsystem that directly contacts with a quantum subsystem is discussed. We also show that the Hamiltonian thus represented gives a theoretical foundation to examine the validity of the so-called semiclassical Ehrenfest theory (or mean-field theory) for electron quantum wavepacket dynamics, and indeed, it is pointed out that the electronic Hamiltonian to be used in this theory should be slightly modified.     

Approach to nonadiabatic transitions by density matrix evolution and molecular dynamics simulations

Many biological processes are characterized by an essentially quantum dynamical event, such as electron or proton transfer, in a complex classical environment. To treat such processes properly by computer simulation, allowing nonadiabatic transitions involving excited states, we recently developed a density matrix evolution (DME ) method [H. J. C. Berendsen and J. Mavri, J. Phys. Chem, 97 , 13464 (1993)] which simulates the dynamics of quantum systems embedded in a classical environment. The formalism of the method is presented and an overview of the applications ranging from collisions of a quantum harmonic oscillator with noble gas atoms to proton tunneling in a double-well hydrogen bond is given. The methodology for treatment of proton-transfer processes with inclusion of excited states is presented. Future application of the method on biologically interesting processes, such as proton transfer in enzymatic reactions, is discussed. © 1996 John Wiley & Sons, Inc.     

Vibrational energy relaxation of polyatomic molecules in liquid solution via the linearized semiclassical method

Vibrational energy relaxation (VER) of polyatomic, as opposed to diatomic, molecules can occur via different, often solvent assisted, intramolecular and/or intermolecular pathways. In this paper, we apply the linearized semiclassical (LSC) method for calculating VER rates in the prototypical case of a rigid, symmetrical and linear triatomic molecule (A-B-A) in a monatomic liquid. Starting at the first excited state of either the symmetric or asymmetric stretches, VER can occur either directly to the ground state or indirectly via intramolecular vibrational relaxation (IVR). The VER rate constants for the various pathways are calculated within the framework of the Landau-Teller formalism, where they are expressed in terms of two-time quantum-mechanical correlation functions. The latter are calculated by the LHA-LSC method, which puts them in a "Wignerized" form, and employs a local harmonic approximation (LHA) in order to compute the necessary multidimensional Wigner integrals. Results are reported for the LHL/Ar model of Deng and Stratt [J. Chem. Phys. 2002, 117, 1735], as well as for CO(2) in liquid argon and in liquid neon. The LHA-LSC method is shown to give rise to significantly faster VER and IVR rates in comparison to the classical treatment, particularly at lower temperatures. We also find that the type and extent of the quantum rate enhancement is strongly dependent on the particular VER pathway. Finally, we find that the classical and semiclassical treatments can give rise to opposite trends when it comes to the dependence of the VER rates on the solvent.     

A coherent state approach to semiclassical nonadiabatic dynamics

A semiclassical (SC) approximation to the quantum mechanical propagator for nonadiabatic systems is derived. Our derivation starts with an exact path integral expression that uses canonical coherent states for the nuclear degrees of freedom and spin coherent states for the electronic degrees of freedom. A stationary path approximation (SPA) is then applied to the path integral to obtain the SC approximation. The SPA results in complex classical trajectories of both nuclear and electronic degrees of freedom and a double ended boundary condition. The root search problem is solved using the previously proposed "real trajectory local search" algorithm. The SC approximation is tested on three simple one dimensional two-state systems proposed by Tully [J. Chem. Phys. 93, 1061 (1990)], and the SC results are compared to Ehrenfest and surface hopping predictions. Excellent agreement with quantum results is reached when the SC trajectory is far away from caustics. We discuss the origin of caustics in this SC formalism and the strengths and weaknesses of this approach.     

Statistical mechanical theory for steady state systems. VI. Variational principles

Several variational principles that have been proposed for nonequilibrium systems are analyzed. These include the principle of minimum rate of entropy production due to Prigogine [Introduction to Thermodynamics of Irreversible Processes (Interscience, New York, 1967)], the principle of maximum rate of entropy production, which is common on the internet and in the natural sciences, two principles of minimum dissipation due to Onsager [Phys. Rev. 37, 405 (1931)] and to Onsager and Machlup [Phys. Rev. 91, 1505 (1953)], and the principle of maximum second entropy due to Attard [J. Chem.. Phys. 122, 154101 (2005); Phys. Chem. Chem. Phys. 8, 3585 (2006)]. The approaches of Onsager and Attard are argued to be the only viable theories. These two are related, although their physical interpretation and mathematical approximations differ. A numerical comparison with computer simulation results indicates that Attard's expression is the only accurate theory. The implications for the Langevin and other stochastic differential equations are discussed.     

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)].