CI techniques, based on nearly complete expansions over determinants for subsets of molecular orbitals and electrons, have been extended to evaluate, by ® nite di erences, nonadiabatic interactions arising from the nuclear kinetic energy operator. A functional of the non-adiabatic interactions is de® ned and, from the minimum condition of the functional, a simple expression is obtained for the transformation matrix to the quasi-diabatic basis. Two di erent schemes for the minimization of the functional, with and without constraints to the form of the transformation matrix, have been compared. Simple tests are presented to verify the numerical accuracy of the non-adiabatic interactions and to compare the two schemes for the minimization of the functional.