In what has been described as the fourth age of quantum chemistry, variational nuclear motion programs are now routinely being used to obtain the vibration-rotation levels and corresponding wavefunctions of small molecules to the sort of high accuracy demanded by comparison with spectroscopy. In this perspective, I will discuss the current state-of-the-art which, for example, shows that these calculations are increasingly competitive with measurements or, indeed, replacing them and thus becoming the primary source of data on key processes. To achieve this accuracy ab initio requires consideration of small effects, routinely ignored in standard calculations, such as those due to quantum electrodynamics. Variational calculations are being used to generate huge lists of transitions which provide the input for models of radiative transport through hot atmospheres and to fill in or even replace measured transition intensities. Future prospects such as the study of molecular states near dissociation, which can provide a link with low-energy chemical reactions, are discussed. Published by AIP Publishing. [http://dx