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.

Orbital Surface Hopping from Orbital Quantum-Classical Liouville Equation for Nonadiabatic Dynamics of Many-electron Systems

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

Efficient and accurate surface hopping for long time nonadiabatic quantum dynamics

Nonadiabatic Dynamics Near Metal Surface With Periodic Drivings: A Floquet Surface Hopping Algorithm

Surface hopping methods for nonadiabatic dynamics in extended systems

Multilayer Subsystem Surface Hopping Method for Large-Scale Nonadiabatic Dynamics Simulation with Hundreds of Thousands of States

Recent Progress in Surface Hopping: 2011–2015

Nonadiabatic Dynamics of Molecules Interacting with Metal Surfaces: A Quantum-Classical Approach Based on Langevin Dynamics and the Hierarchical Equations of Motion

Subspace Surface Hopping with Size-Independent Dynamics.

Nonadiabatic dynamics near metal surfaces with periodic drivings: A generalized surface hopping in Floquet representation

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

Branching and Phase Corrected Surface Hopping: A Benchmark of Nonadiabatic Dynamics in Multilevel Systems

Nonadiabatic dynamics with quantum nuclei: simulating charge transfer with ring polymer surface hopping

Surface Hopping without Momentum Jumps: A Quantum-Trajectory-Based Approach to Nonadiabatic Dynamics

A semiclassical non-adiabatic phase-space approach to molecular translations and rotations: A new picture of surface hopping and electronic inertial effects

An efficient implementation of spin–orbit coupling within the framework of semiempirical orthogonalization-corrected methods for ultrafast intersystem crossing dynamics

First-principles Nonadiabatic Dynamics of Molecules at Metal Surfaces with Vibrationally Coupled Electron Transfer

A generalized surface hopping method

An Efficient, Augmented Surface Hopping Algorithm That Includes Decoherence for Use in Large-Scale Simulations

Quantum Quality with Classical Cost: Ab Initio Nonadiabatic Dynamics Simulations using the Mapping Approach to Surface Hopping