ABSTRACT:We propose an accurate and efficient method to compute vibrational spectra of molecules, based on exact diagonalization of an algebraically calculated matrix based on powers of Morse coordinate. The present work focuses on the 1D potential of diatomic molecules: as typical examples, we apply this method to the standard Lennard-Jones oscillator, and to the ab initio potential of the H 2 molecule. Global cm −1 accuracy is exhibited through the H 2 spectrum, obtained through the diagonalization of a 30 × 30 matrix. This theory is at the root of a new method to obtain globally accurate vibrational spectral data in the context of the multi-dimensional potential of polyatomic molecules, at an affordable computational cost.