Correlated Displacement of Dynamic Elastic Dipoles Produces Nonclassical Electrostriction in Zr-Doped Ceria

By combining experimental data with density functional theory-based ab initio molecular dynamics modeling, this work provides evidence that nonclassical electrostriction in isovalent Zr-doped ceria is due to the correlated anharmonic motion of dynamic elastic dipoles associated with multiple [ZrO₈]-local bonding units with a high Zr concentration (Zr_0.1Ce_0.9O₂). Introduction of 0.5 mol % trivalent or divalent codopants (Sc, Yb, La, or Ca) reduces the longitudinal electrostriction strain coefficient by more than a factor of 10, produces a 3-fold decrease in the relative dielectric permittivity, and increases the elastic modulus. Since these changes depend neither on the radius nor on the valency of the codopant, we conclude that the responsible species are charge-compensating oxygen vacancies (V_O). For trivalent dopants (Do_0.005Zr_0.1Ce_0.895O_1.9975), oxygen vacancies are present at a concentration ratio 1:40 with respect to Zr, giving, for random distribution, a characteristic interaction distance of ≤2.3 unit cells (1.2 nm). Oxygen vacancies participate in [ZrO₇-V_O] local bonding units, disrupting the correlated dynamic displacements of the connected [ZrO₈]-local bonding units. Such correlated motion of dynamic elastic dipoles may also explain the exponential increase in the longitudinal electrostriction strain coefficient with an increase in Zr concentration to <0.2 mole fraction and must be taken into account for further development of nonclassical electrostrictors based on Zr-doped ceria.