Many problems in quantumchemistry deal with the computation of fundamental or excited states of molecules and lead to the resolution of eigenvalue problems. One of the major difficulties in these computations lies in the very largedimension of the systems to be solved. Indeed these eigenfunctions dependon 3n variables where n stands for the number of particles (electrons and/or nucleari) in the molecule. In order to diminish the size of the systems to be solved, the chemists have proposed manyinteresting ideas. Among those stands the adiabatic variable method; we present in this paper a mathematical analysis of thisapproximation and propose, in particular, an a posteriori estimate thatmight allow for verifying the adiabaticity hypothesis that is done on some variables; numerical simulations that support thea posteriori estimators obtained theoretically are also presented.