where ω qv are the phonon frequencies at wave vector q and band index v, are then calculated with PHONOPY. Due to the computational cost of calculating phonon dispersions for supercells, we opt to calculate them only for a subset of the structures. We then construct a surrogate model based only on the internal energy, E 0 , and equilibrium volume, V 0 relative to the ground state structure, to reproduce the relative free energy as