Abstract: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 timesteps. 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 (κCSDM) and curvature-driven trajectory surface hopping (κTSH). We show how powerful these new methods are in terms of computer time and good agreement with methods employing nonadiabatic coupling vectors computed in conventional ways. 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.

Optimization-based Shrinking Dimer Method for Finding Transition States

Global optimization-based dimer method for finding saddle points

Constrained Broyden Minimization Combined with the Dimer Method for Locating Transition State of Complex Reactions

A combined steepest descent and genetic algorithm (SD/GA) approach for the optimization of solvation parameters

Finite-Temperature Dimer Method for Finding Saddle Points on Free Energy Surfaces

Efficient Chain-of-States Approach for Locating Transition State via Spherical Optimization.

Variable-cell Double-Ended Surface Walking Method for Fast Transition State Location of Solid Phase Transitions.

Diagonalization without Diagonalization: A Direct Optimization Approach for Solid-State Density Functional Theory

Adaptive nudged elastic band approach for transition state calculation

Double-Ended Surface Walking Method for Pathway Building and Transition State Location of Complex Reactions

Optimization of Convergence Criteria for Fragmentation Methods

Understanding Solid-State Solvation-Enhanced Thermally Activated Delayed Fluorescence Using a Descriptor-Tuned Screened Range-Separated Functional

A Simple Density Functional Fractional Occupation Number Procedure to Determine the Low Energy Transition Region of Spin-Flip Reactions

Force reversed method for locating transition states

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

Structure Optimization with Stochastic Density Functional Theory

Efficient Softest Mode Finding in Transition States Calculations

High-Index Optimization-Based Shrinking Dimer Method for Finding High-Index Saddle Points

Enhancing initial state overlap through orbital optimization for faster molecular electronic ground-state energy estimation

Efficient Geometry Minimization and Transition Structure Optimization Using Interpolated Potential Energy Surfaces and Iteratively Updated Hessians.

Constrained Broyden Dimer Method with Bias Potential for Exploring Potential Energy Surface of Multistep Reaction Process.