Phase corrected higher-order expression for surface hopping transition amplitudes in nonadiabatic scattering problems

Improving the efficiency of Monte Carlo surface hopping calculations

A surface hopping method with a Monte Carlo procedure for deciding whether to hop at each step along the classical trajectories used in the semiclassical calculation is discussed. It is shown for a simple one-dimensional model problem that the numerical efficiency of the method can be improved by averaging over several copies of the sections of each trajectory that span the interaction regions. The use of Sobol sequences in the selection of the initial momentum for the trajectories is also explored. It is found that accurate results can be obtained with relatively small trajectory samples.

A semiclassical model for the calculation of nonadiabatic transition probabilities for classically forbidden transitions

A semiclassical surface hopping model is presented for the calculation of nonadiabatic transition probabilities for the case in which the avoided crossing point is in the classically forbidden regions. The exact potentials and coupling are replaced with simple functional forms that are fitted to the values, evaluated at the turning point in the classical motion, of the Born-Oppenheimer potentials, the nonadiabatic coupling, and their first few derivatives. For the one-dimensional model considered, reasonably accurate results for transition probabilities are obtained down to around 10(-10). The possible extension of this model to many dimensional problems is discussed. The fact that the model requires only information at the turning point, a point that the trajectories encounter would be a significant advantage in many dimensional problems over Landau-Zener type models, which require information at the avoided crossing seam, which is in the forbidden region where the trajectories do not go.

Semiclassical nonadiabatic surface-hopping wave function expansion at low energies: hops in the forbidden region

The accuracy of a semiclassical surface-hopping expansion of the time-independent wave function for problems in which the nonadiabatic coupling is peaked in the classically forbidden regions is studied numerically for a one-dimensional curve-crossing problem. This surface-hopping expansion has recently been shown to satisfy the Schrodinger equation to all orders in h and all orders in the nonadiabatic coupling. It has also been found to provide very accurate transition probabilities for problems in which the crossing points of the diabatic energy surfaces are classically allowed. In the numerical study reported here, transition probabilities are evaluated for energies well below the crossing point energy. It is found that the expansion provides accurate results for transition probabilities as small as 10(-11).

Toward an accurate and efficient semiclassical surface hopping procedure for nonadiabatic problems

The derivation of a semiclassical surface hopping procedure from a formally exact solution of the Schrodinger equation is discussed. The fact that the derivation proceeds from an exact solution guarantees that all phase terms are completely and accurately included. Numerical evidence shows the method to be highly accurate. A Monte Carlo implementation of this method is considered, and recent work to significantly improve the statistical accuracy of the Monte Carlo approach is discussed.

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.

An analysis of the accuracy of an initial value representation surface hopping wave function in the interaction and asymptotic regions

The behavior of an initial value representation surface hopping wave function is examined. Since this method is an initial value representation for the semiclassical solution of the time independent Schrodinger equation for nonadiabatic problems, it has computational advantages over the primitive surface hopping wave function. The primitive wave function has been shown to provide transition probabilities that accurately compare with quantum results for model problems. The analysis presented in this work shows that the multistate initial value representation surface hopping wave function should approach the primitive result in asymptotic regions and provide transition probabilities with the same level of accuracy for scattering problems as the primitive method.

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.

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.