“…Since the lattice parameter that minimizes the free energy of the system depends on temperature (thermal expansion), it would be necessary, in order to refine the thermodynamic description of the system at zero pressure, to calculate phonon frequencies for different lattice parameters and take those values that minimize the vibration free energy, which constitutes the quasiharmonic approximation. However, this kind of study has shown [44] that even at relatively elevated temperature (800 K) the free energy difference between the harmonic (adopted in this work and that consists in calculating F vib at fixed volume) and quasiharmonic approximation is only a few meV/atom, a degree of refinement of second 7 order with respect to other uncertainties in the ab initio input data. We therefore opted to work at constant lattice parameter (the 0 K equilibrium value) to preserve computational resources.…”