Lattice dynamics calculations are performed using the finite displacement method as implemented in PHONOPY . For a subset of structures for each system, the phonon dispersions are calculated with DFT for validation using the same settings as mentioned above in supercells containing around 100 atoms. Lattice dynamics calculations are then performed for all structures used for training the cluster expansions using the trained MLFFs, and the vibrational free energy is calculated.