Comparison of the different distribution functions in Gd-doped ceria system by molecular dynamics simulations


 Cerium dioxide-based materials are among the most studied for applications in the energy and environmental fields and are also of interest in biology and medicine. The fluorite structure of CeO2 is locally distorted by the concomitant presence of doping cations, such as Gd3+ and oxygen vacancies. The cation-anion bond length distribution then becomes increasingly asymmetric with the doping ratio and temperature. In these cases, the Molecular Dynamics simulation results indicate that the commonly used maximum of the pair correlation function g(r) first peak can no longer be adopted to estimate the mean bond length. To determine the true cation-anion bond length, the analysis of the Radial Distribution Function R(r)first peak is necessary. Furthermore, the asymmetry of this peak must be accounted for when extracting the mean value of the distribution. The gap between the g(r)maximum and the R(r)mean position derived from the fit using a skewed Gaussian function clearly increases with the doping ratio and temperature, leading to different conclusions concerning bond length evolution. The present study also suggests that care must be taken when the bond length is deduced from the pair distribution functions Gpdf(r) as is the case in total scattering experiments (X-ray and neutrons). Finally, relations between the cumulants of the effective distribution of distances as determined in EXAFS experiments and the cumulants describing the real distance distribution are proposed considering that both these distributions are modelled by a skewed Gaussian function.