Oxygen vacancies versus fluorine at CeO2(111): a case of mistaken identity?

We propose a resolution to the puzzle presented by the surface defects observed with STM at the (111) surface facet of CeO 2 single crystals. In the seminal paper of Esch et al. [Science 309, 752 (2005)] they were identified with oxygen vacancies, but the observed behavior of these defects is inconsistent with the results of density functional theory (DFT) studies of oxygen vacancies in the literature. We resolve these inconsistencies via DFT calculations of the properties of both oxygen vacancies and fluorine impurities at CeO2(111), the latter having recently been shown to exist in high concentrations in single crystals from a widely used commercial source of such samples. We find that the simulated filled-state STM images of surface-layer oxygen vacancies and fluorine impurities are essentially identical, which would render problematic their experimental distinction by such images alone. However, we find that our theoretical results for the most stable location, mobility, and tendency to cluster, of fluorine impurities are consistent with experimental observations, in contrast to those for oxygen vacancies. Based on these results, we propose that the surface defects observed in STM experiments on CeO2 single crystals reported heretofore were not oxygen vacancies, but fluorine impurities. Since the similarity of the simulated STM images of the two defects is due primarily to the relative energies of the 2p states of oxygen and fluorine ions, this confusion might also occur for other oxides which have been either doped or contaminated with fluorine.

Tuning LDA+U for electron localization and structure at oxygen vacancies in ceria.

We examine the real space structure and the electronic structure (particularly Ce4f electron localization) of oxygen vacancies in CeO(2) (ceria) as a function of U in density functional theory studies with the rotationally invariant forms of the LDA+U and GGA+U functionals. The four nearest neighbor Ce ions always relax outwards, with those not carrying localized Ce4f charge moving furthest. Several quantification schemes show that the charge starts to become localized at U approximately 3 eV and that the degree of localization reaches a maximum at approximately 6 eV for LDA+U or at approximately 5.5 eV for GGA+U. For higher U it decreases rapidly as charge is transferred onto second neighbor O ions and beyond. The localization is never into atomic corelike states; at maximum localization about 80-90% of the Ce4f charge is located on the two nearest neighboring Ce ions. However, if we look at the total atomic charge we find that the two ions only make a net gain of (0.2-0.4)e each, so localization is actually very incomplete, with localization of Ce4f electrons coming at the expense of moving other electrons off the Ce ions. We have also revisited some properties of defect-free ceria and find that with LDA+U the crystal structure is actually best described with U=3-4 eV, while the experimental band structure is obtained with U=7-8 eV. (For GGA+U the lattice parameters worsen for U>0 eV, but the band structure is similar to LDA+U.) The best overall choice is U approximately 6 eV with LDA+U and approximately 5.5 eV for GGA+U, since the localization is most important, but a consistent choice for both CeO(2) and Ce(2)O(3), with and without vacancies, is hard to find.