The atomistic structure of yttria stabilised zirconia at 6.7 mol%: an ab initio study.
File(s)
Accepted version
Author(s)
Type
Journal Article
Abstract
Yttria stabilized zirconia (YSZ) is an important oxide ion conductor used in solid oxide fuel cells, oxygen sensing devices, and for oxygen separation. Doping pure zirconia (ZrO2) with yttria (Y2O3) stabilizes the cubic structure against phonon induced distortions and this facilitates high oxide ion conductivity. The local atomic structure of the dopant is, however, not fully understood. X-ray and neutron diffraction experiments have established that, for dopant concentrations below 40 mol% Y2O3, no long range order is established. A variety of local structures have been suggested on the basis of theoretical and computational models of dopant energetics. These studies have been restricted by the difficulty of establishing force field models with predictive accuracy or exploring the large space of dopant configurations with first principles theory. In the current study a comprehensive search for all symmetry independent configurations (2857 candidates) is performed for 6.7 mol% YSZ modelled in a 2 × 2 × 2 periodic supercell using gradient corrected density functional theory. The lowest energy dopant structures are found to have oxygen vacancy pairs preferentially aligned along the 〈210〉 crystallographic direction in contrast to previous results which have suggested that orientation along the 〈111〉 orientation is favourable. Analysis of the defect structures suggests that the Y(3+)-Ovac interatomic separation is an important parameter for determining the relative configurational energies. Current force field models are found to be poor predictors of the lowest energy structures. It is suggested that the energies from a simple point charge model evaluated at unrelaxed geometries is actually a better descriptor of the energy ordering of dopant structures. Using these observations a pragmatic procedure for identifying low energy structures in more complicated material models is suggested. Calculation of the oxygen vacancy migration activation energies within the lowest energy 〈210〉 oriented structures gives results consistent with experimental observations.
Date Issued
2016-12-07
Date Acceptance
2016-11-02
Citation
Physical Chemistry Chemical Physics, 2016, 18, pp.31277-31285
ISSN
1463-9084
Publisher
Royal Society of Chemistry
Start Page
31277
End Page
31285
Journal / Book Title
Physical Chemistry Chemical Physics
Volume
18
Copyright Statement
© 2016 Royal Society of Chemistry
Sponsor
Engineering & Physical Science Research Council (EPSRC)
Identifier
http://www.ncbi.nlm.nih.gov/pubmed/27819074
Grant Number
EP/I00422X/1
Subjects
Chemical Physics
02 Physical Sciences
03 Chemical Sciences
Publication Status
Published
Coverage Spatial
England
Date Publish Online
2016-11-02