Abstract:We have recently introduced an efficient semi-empirical non-adiabatic molecular dynamics method for the simulation of charge transfer/transport in molecules and molecular materials, denoted fragment orbital-based surface hopping (FOB-SH) [J. Spencer et al., J. Chem. Phys. 145, 064102 (2016)]. In this method, the charge carrier wavefunction is expanded in a set of charge localized, diabatic electronic states and propagated in the time-dependent potential due to classical nuclear motion. Here we derive and implement an exact expression for the non-adiabatic coupling vectors between the adiabatic electronic states in terms of nuclear gradients of the diabatic electronic states. With the non-adiabatic coupling vectors (NACVs) available, we investigate how different flavours of fewest switches surface hopping affect detailed balance, internal consistency, and total energy conservation for electron hole transfer in a molecular dimer with two electronic states. We find that FOB-SH satisfies detailed balance across a wide range of diabatic electronic coupling strengths provided that the velocities are adjusted along the direction of the NACV to satisfy total energy conservation upon a surface hop. This criterion produces the right fraction of energy-forbidden (frustrated) hops, which is essential for correct population of excited states, especially when diabatic couplings are on the order of the thermal energy or larger, as in organic semiconductors and DNA. Furthermore, we find that FOB-SH is internally consistent, that is, the electronic surface population matches the average quantum amplitudes, but only in the limit of small diabatic couplings. For large diabatic couplings, inconsistencies are observed as the decrease in excited state population due to frustrated hops is not matched by a corresponding decrease in quantum amplitudes. The derivation provided here for the NACV should be generally applicable to any electronic structure approach where the electronic Hamiltonian is constructed in a diabatic electronic state basis.

Adiabatic electronic flux density: a Born-Oppenheimer Broken Symmetry ansatz

Nonadiabatic Dynamics Via the Symmetrical Quasi-Classical Method in the Presence of Anharmonicity.

Can The Mystery of The Born-Oppenheimer Electronic Current Density Be Explained With A Simple Phase Space Electronic Hamiltonian? Yes (And A Lot More Too)

Density functional theory beyond the Born-Oppenheimer approximation: Exact mapping onto an electronically non-interacting Kohn-Sham molecule

Mixed quantum–classical approach to model non-adiabatic electron–nuclear dynamics: Detailed balance and improved surface hopping method

Different flavors of nonadiabatic molecular dynamics

Born-Oppenheimer Dynamics, Electronic Friction, and the Inclusion of Electron-Electron Interactions

Gaussians for Electronic and Rovibrational Quantum Dynamics

Fully quantum non-adiabatic dynamics in electronic-nuclear coherent state basis

Diagonalizing the Born-Oppenheimer Hamiltonian via Moyal Perturbation Theory, Nonadiabatic Corrections and Translational Degrees of Freedom

Detailed balance, internal consistency, and energy conservation in fragment orbital-based surface hopping

Application of the Adiabatic Connection Random Phase Approximation to Electron-Nucleus Hyperfine Coupling Constants

Surface Hopping, Electron Translation Factors, Electron Rotation Factors, Momentum Conservation, and Size Consistency

Density functional theory of electron transfer beyond the Born-Oppenheimer approximation: Case study of LiF

On the representation of coupled adiabatic potential energy surfaces using quasi-diabatic Hamiltonians: a distributed origins expansion approach

Simulating Attochemistry: Which Dynamics Method to Use?

Total angular momentum conservation in Ehrenfest dynamics with a truncated basis of adiabatic states

The Numerical Equivalence of Diabatic and Adiabatic Representations in Diatomic Molecules

Adiabatic Dynamics of Coupled Spins and Phonons in Magnetic Insulators

Interstate Vibronic Coupling Constants Between Electronic Excited States for Complex Molecules

Diagonalizing the Born–Oppenheimer Hamiltonian via Moyal perturbation theory, nonadiabatic corrections, and translational degrees of freedom