Derivative Couplings with Built-In Electron-Translation Factors: Application to Benzene

Derivative couplings are the essential quantities at the interface between electronic-structure calculations and nonadiabatic dynamics. Unfortunately, standard approaches for calculating these couplings usually neglect electronic motion, which can lead to spurious electronic transitions. Here we provide a general framework for correcting these anomalies by incorporating perturbative electron-translation factors (ETFs) into the atomic-orbital basis. For a range of representative organic molecules, we find that our ETF correction is often small but can be qualitatively important, especially for few-atom systems or highly symmetric molecules. Our method entails no additional computational cost, such that ETFs are “built-in,” and it is equivalent to a simple rule of thumb: We should set the antisymmetrized version of the nuclear overlap-matrix derivative to zero wherever it appears. Thus, we expect that built-in ETFs will be regularly incorporated into future studies of nonadiabatic dynamics

Derivative couplings and analytic gradients for diabatic states, with an implementation for Boys-localized configuration-interaction singles

We demonstrate that Boys-localized diabatic states do indeed exhibit small derivative couplings, as is required of quasidiabatic states. In doing so, we present a general formalism for calculating derivative couplings and analytic gradients for diabatic states. We then develop additional equations specific to the case of Boys-localized configuration-interaction singles (CIS)—in particular, the analytic gradient of the CIS dipole matrix—and we validate our implementation against finite-difference results. In a forthcoming paper, we will publish additional algorithmic and computational details and apply our method to the Closs energy-transfer systems as a further test of the validity of Boys-localized diabatic states