Keywords: solid oxide fuel cell, yttria stabilised zirconia, cubic zirconia, bixbyite yttria, ab initio, density functional theory, empirical potential, Born-Mayer-Huggins, point charge model, phonons, harmonic approximation. been observed in either X-ray or neutron diffraction experiments. The prediction of local defect structure and the interaction between defects is therefore of great interest. This has not been possible to date as the number of possible defect topologies is very large and to perform reliable total energy calculations for all of them would be prohibitively expensive. Previous theoretical studies have only considered a selection of representative structures. In this study, a comprehensive search for low energy defect structures using a combined classical modelling and density functional theory approach is used to identify the low energy isolated defect structures at the dilute limit, 3.2mol%. Through analysis of energetics computed using the best available Born-Mayer-Huggins empirical potential model, a point charge model, DFT, and a local strain energy estimated in the harmonic approximation, the main chemical and physical descriptors that correlate to the low energy DFT structures are discussed. It is found that the empirical potential model reproduces a general trend of increasing DFT energetics across a series of locally strain relaxed structures, but is unreliable both in predicting some incorrect low energy structures, and 3 in finding some meta-stable structures to be unstable. A better predictor of low energy defect structures is found to be the total electrostatic energy of a simple point charge model calculated at the unrelaxed geometries of the defects. In addition, the strain relaxation energy is estimated effectively in the harmonic approximation to the imaginary phonon modes of undoped c-ZrO 2 , but is found to be unimportant in determining the low energy defect structures. These results allow us to propose a set of easily computed descriptors that can be used to identify the low energy YSZ defect structures, negating the combinatorial complexity and number of defect structures that need to be considered.
Abstract