The matrix in right hand side can be seen as the nuclear kinetic operator in the adiabatic-Fock states basis except a constant -2 /2M . Since we already derived the formula of V in the same basis in Eq. 36, we can directly plug these expressions into the TISE Ĥ|Ψ = ε|Ψ , and obtain the following eigenvalue problem