The approximation of the full nuclear potential of a molecular system in its equilibrium geometry and in the vicinity of the respective minimum on the potential energy surface by the function \[V=\frac{1}{2}\ \sum _{\begin{array}{c} i,j=1 \end{array}}^{3\ N}\frac{\partial ^{2}V}{\partial q_{i}\partial q_{j}}\] where q i are mass-weighted cartesian displacements of nuclei relative to their equilibrium positions. The approximation allows one to describe vibrational motion in terms of independent vibrational modes (normal modes) each of which is governed by a simple one-dimensional harmonic potential.