Ordering of oxygen vacancies and excess charge localization in bulk ceria: A DFT+U study

G. E. MURGIDA; V. FERRARI; M. VERÓNICA GANDUGLIA PIROVANO; A. M. LLOIS

PHYSICAL REVIEW B

AMER PHYSICAL SOC

Lugar: New York; Año: 2014 vol. 90 p. 115120 - 115120

1098-0121

The importance of ceria (CeO$_2$) in many applications originates from the ease of oxygen vacancy formation and healing. The ordering of vacancies and the whereabouts of the excess charge in bulk CeO$_2$ are of no less significance than at ceria surfaces, but they have not received the same attention.  In this work, the formation of neutral oxygen vacancies in bulk CeO$_2$ are investigated using density-functional theory (DFT) in the DFT$+U$ ($U$ is an effective onsite Coulomb interaction parameter) approach for a broad range of vacancy concentrations, $Theta$ $(rac{1}{64}leqThetaleqrac{1}{4})$. We find   that the excess charge prefers to be localized in cation sites such that the mean Ce$^{3+}$ coordination number is maximized, and if nearest neighbor cation sites are reduced, they rather be non-uniformly distributed. Furthermore, we show that a vacancy repels other vacancies from its nearest neighbor shell and that the [110] and  [111] directions are possible directions for clustering of second and third neighbor vacancies,  respectively. Vacancies prefer not to share cations. The results are discussed in a simple physical picture which enables the separation of the different contributions to the averaged vacancy formation energy.  We also consider cells with fluorite structure and same stoichiometries as in existing bulk phases, i.e.,  Ce$_{11}$O$_{20}$ ($Theta=rac{1}{11}$), Ce$_{7}$O$_{12}$ ($Theta=rac{1}{7}$),  and Ce$_{2}$O$_{3}$ ($Theta=rac{1}{4}$), as well as the corresponding real structures. We find that  the vacancy ordering and the location of the excess electrons are consistent with the results for single-phase reduced CeO$_2$, but  the Ce$_{11}$O$_{20}$, Ce$_{7}$O$_{12}$,  and  Ce$_{2}$O$_{3}$ structures are substantially more stable. The stability of these phases as a function of pressure and temperature is discussed. Vacancy-induced lattice relaxations effects are crucial for the interpretation of the results.