“…A finely tuned ZrO 2 and Y 2 O 3 interatomic potential that describes the entire Zr 1−x Y x O 2−x/2 system, with x/2 being the Y 2 O 3 dopant concentration of YSZ has to be established. According to the reported literature (Bush et al, 1994;Devanathan et al, 2006;Dwivedi et al, 1990;Fisher et al, 1998Fisher et al, , 1999Khan et al, 1998;Kilo et al, 2003;Lau et al, 2011;Lewis et al, 1985;Li et al, 1995;Minervini et al, 2000;Okazaki et al, 1994;Sawaguchi et al, 2000;Schelling et al, 2001;Shimojo et al, 1992;;van Duin et al, 2008;Yamamura et al, 1999;Zacate et al, 2000), the typical "semi-empirically" fitted properties of ZrO 2 and Y 2 O 3 crystals that are chosen for the fitting dataset can be the lattice parameters, lattice elastic properties, dielectric constants, defect formation energies (e.g. vacancies and interstitials), and phonon frequencies of cubic (c-ZrO 2 , space group Fm3m), tetragonal (t-ZrO 2 , space group P42/nmc) (Ackermann et al, 1975;Aldebert et al, 1985;Boysen et al, 1991;Dash et al, 2004;Howard et al, 1988;Smith et al, 1965;Zhao et al, 2002;), monoclinic (m-ZrO2, space group P2 1 /c ), and yttria (Baller et al, 2000; (Bogicevic et al, 2001).…”