La energía obtenida para el oscilador armónico depende del desarrollo (ecuación 2) que se truncó por el segundo término. Añadiendo las derivadas tercera y cuarta el desarrollo se obtiene una mejor aproximación para la energía. El término de corrección para la ecuación de la energía vibracional es: $-h\nu_e x_e(v+\frac{1}{2})^2$, donde $\nu_e x_e$ se conoce como constante de anarmonicidad. La energía vibracional corregida es: $ E_{vib}=(v+\frac{1}{2})h\nu_e-h\nu_e x_e(v+\frac{1}{2})^2$
La inclusión de la anarmonicidad nos deja la energía interna de la molécula como sigue:
\begin{equation}\label{13} E_{int}=B_e h J(J+1)+(v+\frac{1}{2})h\nu_e-h\nu_e x_e(v+\frac{1}{2})^2+E_{ele} \end{equation}