where ĤM = T0 + V ( R0 ) is the Hamiltonian of the reactive molecule whose ground state potential energy surface is taken to be a symmetric double well potential along the reaction coordinate R0 , see Methods, Eq. 25 for details. We use a set of localized diabatic vibrational states