Monte Carlo (MC) simulations were performed with the Metropolis algorithm 30 using the model Hamiltonians described above, i.e. both the one based only on internal energy calculated with DFT, and the one based on free energies. In the second case, the temperature of the Monte Carlo simulation was matched with the free energy from phonon calculations.