On the inclusion of the diagonal Born-Oppenheimer correction in surface hopping methods

Abstract

The diagonal Born-Oppenheimer correction (DBOC) stems from the diagonal second derivative coupling term in the adiabatic representation, and it can have an arbitrary large magnitude when a gap between neighbouring Born-Oppenheimer (BO) potential energy surfaces (PESs) is closing. Nevertheless, DBOC is typically neglected in mixed quantum-classical methods of simulating nonadiabatic dynamics (e.g., fewest-switch surface hopping (FSSH) method). A straightforward addition of DBOC to BO PESs in the FSSH method, FSSH+D, has been shown to lead to numerically much inferior results for models containing conical intersections. More sophisticated variation of the DBOC inclusion, phase-space surface-hopping (PSSH) was more successful than FSSH+D but on model problems without conical intersections. This work comprehensively assesses the role of DBOC in nonadiabatic dynamics of two electronic state problems and the performance of FSSH, FSSH+D, and PSSH methods in variety of one- and two-dimensional models. Our results show that the inclusion of DBOC can enhance the accuracy of surface hopping simulations when two conditions are simultaneously satisfied: 1) nuclei have kinetic energy lower than DBOC and 2) PESs are not strongly nonadiabatically coupled. The inclusion of DBOC is detrimental in situations where its energy scale becomes very high or even diverges, because in these regions PESs are also very strongly coupled. In this case, the true quantum formalism heavily relies on an interplay between diagonal and off-diagonal nonadiabatic couplings while surface hopping approaches treat diagonal terms as PESs and off-diagonal ones stochastically.

10 pages, 16 figures, the version is aligned with the J. Chem. Phys. publication

References in corpus (5)

• Surface hopping from the perspective of quantum-classical Liouville dynamics
• When do we need to account for the geometric phase in excited state dynamics ?
• Why do mixed quantum-classical methods describe short-time dynamics through conical intersections so well? Analysis of geometric phase effects
• On the inclusion of the diagonal Born-Oppenheimer correction in surface hopping methods
• Analysis of geometric phase effects in the quantum-classical Liouville formalism

Cited by in corpus (12)

• DFTBaby: A software package for non-adiabatic molecular dynamics simulations based on long-range corrected tight-binding TD-DFT(B)
• A Walk Through the Approximations of Ab Initio Multiple Spawning
• On the inclusion of the diagonal Born-Oppenheimer correction in surface hopping methods
• Exploring the many-body dynamics near a conical intersection with trapped Rydberg ions
• A Phase-Space Semiclassical Approach for Modeling Nonadiabatic Nuclear Dynamics with Electronic Spin
• Calculating nonadiabatic couplings and Berry's phase by variational quantum eigensolvers
• Geometric phase effects in excited state dynamics through a conical intersection in large molecules: N-dimensional linear vibronic coupling model study
• Shortcuts to Adiabaticity Assisted by Counterdiabatic Born-Oppenheimer Dynamics
• Exploring vibrational ladder climbing in vibronic coupling models: Toward experimental observation of a geometric phase signature of a conical intersection
• A Quantum-Classical Liouville Formalism in a Preconditioned Basis and Its Connection with Phase-Space Surface Hopping
• Collisional effects in the blue wing of Lyman-alpha
• Mixed quantum-classical approach to model non-adiabatic electron-nuclear dynamics: Detailed balance and improved surface hopping method