Facet-dependent stability of near-surface oxygen vacancies and excess charge localization at CeO2surfaces.

To study the dependence of the relative stability of surface (VA) and subsurface (VB) oxygen vacancies with the crystal facet of CeO2, the reduced (100), (110) and (111) surfaces, with two different concentrations of vacancies, were investigated by means of density functional theory (DFT + U) calculations. The results show that the trend in the near-surface vacancy formation energies for comparable vacancy spacings, i.e. (110) < (100) < (111), does not follow the one in the surface stability of the facets, i.e. (111) < (110) < (100). The results also reveal that the preference of vacancies for surface or subsurface sites, as well as the preferred location of the associated Ce3+polarons, are facet- and concentration-dependent. At the higher vacancy concentration, theVAis more stable than theVBat the (110) facet whereas at the (111), it is the other way around, and at the (100) facet, both theVAand theVBhave similar stability. The stability of theVAvacancies, compared to that of theVB, is accentuated as the concentration decreases. Nearest neighbor polarons to the vacant sites are only observed for the less densely packed (110) and (100) facets. These findings are rationalized in terms of the packing density of the facets, the lattice relaxation effects induced by vacancy formation and the localization of the excess charge, as well as the repulsive Ce3+-Ce3+interactions.