# Thermal Properties Calculated with Phonopy

In Phonopy the thermal properties, free energy, heat capacity, and entropy, are calculated from their statistical thermodynamic expressions using phonon frequencies on a sampling mesh in the reciprocal space. According to the Phonopy documentation,

The unit systems of free energy, heat capacity, and entropy are kJ/mol, J/K/mol, and J/K/mol, respectively, where 1 mol means $$N_A \times$$ your input unit cell (not formula unit), i.e. you have to divide the value by number of formula unit in your unit cell by yourself. For example, in MgO (conventional) unit cell, if you want to compare with experimental results in kJ/mol, you have to divide the phonopy output by four.

First, what does $$N_A$$ stand for? I am confused with their statement in the documenation, in the example of MgO (conventional) unit cell, does '4' stand for the formula unit?

Second, suppose that I have the thermal properties of the same material in Bulk and monolayer form, both having one unit cell but with different number of atoms or ions in their unit cells, then should I normalize the values and divide by the number of atoms in the unit cell in order to compare the results?

$$N_A$$ stands for Avogadro's number. Energies are reported in kJ/mol, so they are technically for a mole of unit cells. Conversion factors to other units (kcal/mol, eV, Hartree, etc.) can be found or derived easily.