On the contrary, the most computationally inexpensive and simple vibrational model is the harmonic approximation. This approach assumes that the potential energy surface (PES) in the vicinity of the equilibrium geometry (local minimum) can be expressed with the second-order Taylor expansion. For a diatomic molecule, in which the only vibrational coordinate is the interatomic distance r, it reads as