Atomistic calculation of association energy in doped ceria
Abstract Atomistic simulation based on energy minimization techniques was carried out to study the association energies of more than 300 defect structures of ceria doped with Lu 2 O 3 , Yb 2 O 3 , Er 2 O 3 , Y 2 O 3 , Gd 2 O 3 , Eu 2 O 3 , Sm 2 O 3 , Nd 2 O 3 and La 2 O 3 . The calculation ensemble includes all the possible defect structures in 6.25 mol% doped ceria. It is revealed that a clear preference for vacancy–dopant association occurs energetically for small dopant cations and vacancy–Ce 4+ association for large dopant cations, with the crossover at Gd 3+ . The defect structures with the maximum association effect for different trivalent cations were determined. For Gd 3+ doped ceria which presents the lowest association effect in the trivalent doped ceria, the association energies of four different defect structures are extremely close, which agrees well with a recent EXAFS investigation, that is, Gd 3+ is distributed randomly in the ceria lattice. This work provides insight for the understanding of ionic conductivity in doped ceria.