Permutationally Restrained Diabatization by Machine Intelligence

Simulations of electronically nonadiabatic processes may employ either the adiabatic or diabatic representation. Direct dynamics calculations are usually carried out in the adiabatic basis because the energy, force, and state coupling can be evaluated directly by many electronic structure methods. However, although its straightforwardness is appealing, direct dynamics is expensive when combined with quantitatively accurate electronic structure theories. This generates interest in analytically fitted surfaces to cut the expense, but the cuspidal ridges of the potentials and the singularities and vector nature of the couplings at high-dimensional, nonsymmetry-determined intersections in the adiabatic representation make accurate fitting almost impossible. This motivates using diabatic representations, where the surfaces are smooth and the couplings are also smooth and-importantly-scalar. In a recent previous work, we have developed a method called diabatization by deep neural network (DDNN) that takes advantage of the smoothness and nonuniqueness of diabatic bases to obtain them by machine learning. The diabatic potential energy matrices (DPEMs) learned by the DDNN method yield not only diabatic potential energy surfaces (PESs) and couplings in an analytic form useful for dynamics calculations, but also adiabatic surfaces and couplings in the adiabatic representation can be calculated inexpensively from the transformation. In the present work, we show how to extend the DDNN method to produce good approximations to global permutationally invariant adiabatic PESs simultaneously with DPEMs. The extended method is called permutationally restrained DDNN.