(11) These are the equations implemented in Phonopy, and can be used to describe harmonic properties of any real phonon spectrum, including volume-dependent phonon spectrum representing the effects of thermal expansion. However, when phonon frequencies themselves become explicitly temperature-dependent, ω → ω(T ), Eq. 10 is no longer an accurate description of vibrational free energy, and must be corrected.