where R and P are, respectively, coordinates and momenta of the nuclear DoFs. x i and p i are coordinates and momenta of the mapping electronic variables. i labels the electronic states, F is the total number of these states, V i (R) is the potential energy of i th electronic state, γ denotes the ZPE correction and d i j represents the first-order non-adiabatic coupling (NAC) vector.