Monoclinic phased zirconia (m-ZrO2) thermodynamic property calculations were also carried out using the
Vienna ab initio Simulation Package VASP [17, 19] coupled with PHONOPY [24]. In order to reach the energy
convergence of 10-6 eV, the cutoff energy was set to 600 eV. The Brillouin zone integration was defined as
Monkhorst-Pack mesh with k points of 5