Materials with high oxygen ion conductivity and low electronic conductivity are required for electrolytes in solid oxide fuel cells (SOFC) and high-temperature electrolysis (SOEC). A potential candidate for the electrolytes, which separate oxidation and reduction processes, is rare-earth doped ceria. The prediction of the ionic conductivity of the electrolytes and a better understanding of the underlying atomistic mechanisms provide an important contribution to the future of sustainable and efficient energy conversion and storage. The central aim of this paper is the detailed investigation of the relationship between defect interactions at the microscopic level and the macroscopic oxygen ion conductivity in the bulk of doped ceria. By combining ab initio density functional theory (DFT) with Kinetic Monte Carlo (KMC) simulations, the oxygen ion conductivity is predicted as a function of the doping concentration. Migration barriers are analyzed for energy contributions, which are caused by the interactions of dopants and vacancies with the migrating oxygen vacancy. We clearly distinguish between energy contributions that are either uniform for forward and backward jumps or favor one migration direction over the reverse direction. If the presence of a dopant changes the migration energy identically for forward and backward jumps, the resulting energy contribution is referred to as blocking. If the change in migration energy due to doping is different for forward and backward jumps of a specific ionic configuration, the resulting energy contributions are referred to as trapping. The influence of both effects on the ionic conductivity is analyzed: blocking determines the dopant fraction where the ionic conductivity exhibits the maximum. Trapping limits the maximum ionic conductivity value. In this way, a deeper understanding of the underlying mechanisms determining the influence of dopants on the ionic conductivity is obtained and the ionic conductivity is predicted more accurately. The detailed results and insights obtained here for doped ceria can be generalized and applied to other ion conductors that are important for SOFCs and SOECs as well as solid state batteries.

The rate of oxygen ion jumps in a solid oxide depends not only on the activation energy but also on the pre-exponential factor of diffusion. In order to allow a fully ab initio prediction of the oxygen ion conductivity in pure and samarium doped ceria, we calculated the attempt frequency for an oxygen ion jump from first principles combining DFT+U, the NEB method, phonon calculations and the transition state theory. Different definitions of the jump attempt frequency are presented. The equivalence of the Eyring and the Vineyard method is shown without restriction to the Gamma point. Convergence checks of the phonon mesh reveal that the common reduction to the Gamma point is not sufficient to calculate the attempt frequency. Calculations of Sm doped ceria revealed an increase of the prefactor. The attempt frequency for the constant pressure case in quasi-harmonic approximation is larger than the attempt frequency at constant volume in harmonic approximation. The calculated electronic energies, enthalpies and entropies of migration are in agreement with the experimental diffusion coefficients and activation energies.

The oxygen ion conductivity of polycrystalline samples of Sm‐doped ceria and of Gd‐doped ceria is studied as a function of doping fraction and temperature using impedance spectroscopy allowing the separation of bulk and grain boundary conductivity. The introduction of a fine spacing for the Sm dopant fraction allows the clear identification of the dopant fraction leading to the largest bulk conductivity. At 267°C, the largest bulk conductivity is shown for Ce _0.93 Sm _0.07 O _1.965 . With increasing temperature, indications of an increase in the dopant fraction, which leads to the maximum in conductivity, are found. For the grain boundary conductivity, the maximum appears at larger dopant fractions compared to the bulk conductivity. The largest total conductivity for both dopants is again found for Sm‐doped ceria. In literature, different syntheses and sample preparation methods led to larger total conductivities for Gd‐doped ceria. In this work, we demonstrate that the variation of sintering conditions leads to scattering in the conductivity over one order of magnitude. Finally, we demonstrate that, in nominally pure cerium oxide, impurities dominate the ionic conductivity.

We have studied the free energy migration barriers $\mathrm{\Delta }{F}^{‡}$ for oxygen diffusion in pure ceria and Sm-doped ceria for the temperatures 300, 700, and 1000 K. We used the density functional theory in the generalized gradient approximation and an additional Hubbard $U$ parameter for the Ce $4f$ electronic states. We compare the results for the free energy deduced from three different methods. First, a static harmonic approach is applied in which the temperature dependent vibrational contributions to energy and entropy are deduced from the phonon frequencies of supercells with a fixed volume. Second, a static quasiharmonic approach is used in which a part of the anharmonicity effect is introduced via an implicit dependence of the harmonic frequencies on the thermally expanding cell volume. Third, the free energy barriers are calculated using metadynamics and molecular dynamics in which anharmonicity effects are naturally taken into account. The three methods examined in this study lead to distinctly different results. According to the harmonic approximation, the migration free energy difference $\mathrm{\Delta }{F}^{‡}$ increases with increasing temperature due to an increasing entropic contribution. According to the quasiharmonic approximation, the migration free energy is independent of temperature. Finally, molecular dynamics predicts a thermally induced increase in the migration free energy. We conclude that temperature dependent experimental lattice constants cancel out the increasing entropic contribution with increasing temperature in the static quasiharmonic approach. The full consideration of anharmonicity effects in the metadynamics method again leads to a temperature dependent migration free energy.

Calcium dodecahydro-closo-dodecaborate, CaB₁₂H₁₂, has a percolating Ca migration path with low activation barrier (650 meV) and can be doped with Al, Bi, or a number of trivalent rare-earth cations.