Modeling Strain Distribution at the Atomic Level in Doped Ceria Films with Extended X-ray Absorption Fine Structure Spectroscopy

Ceria doped with trivalent dopants exhibits nonclassical electrostriction, strong anelasticity, and room-temperature (RT) mechanical creep. These phenomena, unexpected for a ceramic material with a large Young's modulus, have been attributed to the generation of local strain in the vicinity of the host Ce cations due to symmetry-breaking point defects, including oxygen vacancies. However, understanding why strain is generated at the host rather than at the dopant site, as well as predicting these effects as a function of dopant size and concentration, remains a challenge. We have used the evolutionary-algorithm-based reverse Monte Carlo modeling to reconcile the experimental data of extended X-ray absorption fine structure and X-ray diffraction in a combined model structure. By extracting the details of the radial distribution function (RDF) around the host (Ce) and trivalent dopants (Sm or Y), we find that RDF of the first-nearest neighbor (1NN) of host and dopant cations as well as the second-nearest neighbor (2NN) of the dopant are each best modeled with two separate populations corresponding to short and long interatomic distances. This heterogeneity indicates that fluorite symmetry is not preserved locally, especially for the dopant first-and second-NN sites, appearing at surprisingly low doping fractions (5 mol % Sm and 10 mol % Y). Given that Ce rather than dopant sites act as the source of local strain for electrostriction and RT creep, we conclude that the environment around the dopant does not respond to electrical and mechanical excitations, likely because of its similarity to the double fluorite structure which has poor electrostrictive and anelastic properties. The trends we observe in the RDFs around the Ce sites as a function of dopant size and concentration suggest that the response of these sites can be controlled by the extent of doping: Increasing dopant size to increase strain magnitude at the 1NN shell of Ce and decreasing dopant fraction to decrease strain propagation to the 2NN shell of Ce should produce stronger electrostrictive response and RT creep.