Mechanical motions of atoms near their equilibrium geometry are approximated by a universal mapping on a collection of independent harmonic oscillators. The respective frequencies are characteristic to the material, for example, lattice phonons or molecular vibrations. The quantum description of atomic motions involves the solution of the Schrödinger equation for the harmonic oscillator in the space of proper wave functions. The energy levels obtained are shown to be equally spaced where the level spacing is proportional to the oscillator frequency. The stationary solutions are identified in terms of Hermite polynomials, demonstrating remarkable differences from the classical harmonic oscillator. The importance of the zero-point energy is emphasized in the context of the stability of chemical bonds, where the Harmonic approximation is shown to be reasonably valid for typical interatomic bonds in standard thermal conditions. This explains the relevance of the harmonic approximation for analyzing the absorption spectrum of infrared radiation by molecules.