; the sums over Latin indices go over the electronic states. The eigenvalues of ĤBO , i.e. the adiabatic or BO potential energy surfaces (PES), are denoted (α) l and ∆ lk is the energy difference between states l and k at the nuclear configuration R (α) (t). The second terms in Eqs (13) and (14) are the corrections arising from EF