Frozen Gaussian approximation with surface hopping for mixed quantum-classical dynamics: A mathematical justification of fewest switches surface hopping algorithms

Abstract

We develop a surface hopping algorithm based on frozen Gaussian approximation for semiclassical matrix Schrödinger equations, in the spirit of Tully's fewest switches surface hopping method. The algorithm is asymptotically derived from the Schrödinger equation with rigorous approximation error analysis. The resulting algorithm can be viewed as a path integral stochastic representation of the semiclassical matrix Schrödinger equations. Our results provide mathematical understanding to and shed new light on the important class of surface hopping methods in theoretical and computational chemistry.

arXiv comment35 pages

Cited by in corpus (5)

• Surface hopping from the perspective of quantum-classical Liouville dynamics
• Path integral molecular dynamics with surface hopping for thermal equilibrium sampling of nonadiabatic systems
• Improved sampling and validation of frozen Gaussian approximation with surface hopping algorithm for nonadiabatic dynamics
• Accelerated sampling by infinite swapping of path integral molecular dynamics with surface hopping
• A Diabatic Surface Hopping Algorithm based on Time Dependent Perturbation Theory and Semiclassical Analysis