How important are the residual nonadiabatic couplings for an accurate simulation of nonadiabatic quantum dynamics in a quasidiabatic representation?

Abstract

Diabatization of the molecular Hamiltonian is a standard approach to removing the singularities of nonadiabatic couplings at conical intersections of adiabatic potential energy surfaces. In general, it is impossible to eliminate the nonadiabatic couplings entirely -- the resulting "quasidiabatic" states are still coupled by smaller but nonvanishing residual nonadiabatic couplings, which are typically neglected. Here, we propose a general method for assessing the validity of this potentially drastic approximation by comparing quantum dynamics simulated either with or without the residual couplings. To make the numerical errors negligible to the errors due to neglecting the residual couplings, we use the highly accurate and general eighth-order composition of the implicit midpoint method. The usefulness of the proposed method is demonstrated on nonadiabatic simulations in the cubic Jahn--Teller model of nitrogen trioxide and in the induced Renner--Teller model of hydrogen cyanide. We find that, depending on the system, initial state, and employed quasidiabatization scheme, neglecting the residual couplings can result in wrong dynamics. In contrast, simulations with the exact quasidiabatic Hamiltonian, which contains the residual couplings, always yield accurate results.

Improved the main text. Added Eqs. (10) and (31)

Cited by in corpus (7)

• Roadmap for Molecular Benchmarks in Nonadiabatic Dynamics
• Representation and Conservation of Angular Momentum in the Born-Oppenheimer Theory of Polyatomic Molecules
• Local diabatic representation of conical intersection quantum dynamics
• Nonadiabatic conical intersection dynamics in the local diabatic representation with Strang splitting and Fourier basis
• The Parallel-Transported (Quasi)-Diabatic Basis
• Ehrenfest dynamics accelerated with SPEED
• On the single-Hessian Gaussian wavepacket dynamics