Nonadiabatic Dynamics Algorithms with Only Potential Energies and Gradients: Curvature-Driven Coherent Switching with Decay of Mixing and Curvature-Driven Trajectory Surface Hopping

Direct dynamics by mixed quantum-classical nonadiabatic methods is an important tool for understanding processes involving multiple electronic states. Very often, the computational bottleneck of such direct simulation comes from electronic structure theory. For example, at every time step of a trajectory, nonadiabatic dynamics requires potential energy surfaces, their gradients, and the matrix elements coupling the surfaces. The need for the couplings can be alleviated by employing the time derivatives of the wave functions, which can be evaluated from overlaps of electronic wave functions at successive time steps. However, evaluation of overlap integrals is still expensive for large systems. In addition, for electronic structure methods for which the wave functions or the coupling matrix elements are not available, nonadiabatic dynamics algorithms become inapplicable. In this work, building on recent work by Baeck and An, we propose new nonadiabatic dynamics algorithms that only require adiabatic potential energies and their gradients. The new methods are named curvature-driven coherent switching with decay of mixing (kappa CSDM) and curvature-driven trajectory surface hopping (kappa TSH). We show how powerful these new methods are in terms of computation time and accuracy as compared to previous mixed quantum-classical nonadiabatic dynamics algorithms. The lowering of the computational cost will allow longer nonadiabatic trajectories and greater ensemble averaging to be affordable, and the ability to calculate the dynamics without electronic structure coupling matrix elements extends the dynamics capability to new classes of electronic structure methods.

Automated summary

Direct dynamics by mixed quantum-classical nonadiabatic methods is crucial for understanding processes with multiple electronic states. However, the computational bottleneck lies in the electronic structure theory. In this study, we propose new algorithms that only require adiabatic potential energies and their gradients, resulting in improved computation time and accuracy compared to previous methods. This reduction in computational cost allows for longer nonadiabatic trajectories and extends the dynamics capability to new electronic structure methods.