A system of vibrations in a crystalline solid lattice can be modeled as an Einstein solid, i.e. by considering N quantum harmonic oscillator potentials along each degree of freedom. Then, the free energy of the system can be written as