Non-adiabatic vibrational calculations performed with the accuracy of 0.2 cm−1 spanning the whole energy spectrum up to the dissociation limit for 7LiH are reported. A so far unknown energy level is predicted. The key feature of the approach used in the calculations is a valence-bond (VB) based procedure for determining the effective masses of the two vibrating atoms, which depend on the internuclear distance, R. It is found that all LiH electrons participate in the vibrational motion. The R-dependent masses are obtained from the analysis of the simple VB two-configuration ionic-covalent representation of the electronic wave function. These findings are consistent with an interpretation of the chemical bond in LiH as a quantum mechanical superposition of one-electron ionic and covalent states.