Atomic and electronic structure of unreduced and reduced CeO2 surfaces: a first-principles study

The atomic and electronic structure of (111), (110), and (100) surfaces of ceria (CeO2) were studied using density-functional theory within the generalized gradient approximation. Both stoichiometric surfaces and surfaces with oxygen vacancies (unreduced and reduced surfaces, respectively) have been examined. It is found that the (111) surface is the most stable among the considered surfaces, followed by (110) and (100) surfaces, in agreement with experimental observations and previous theoretical results. Different features of relaxation are found for the three surfaces. While the (111) surface undergoes very small relaxation, considerably larger relaxations are found for the (110) and (100) surfaces. The formation of an oxygen vacancy is closely related to the surface structure and occurs more easily for the (110) surface than for (111). The preferred vacancy location is in the surface layer for CeO2(110) and in the subsurface layer (the second O-atomic layer) for CeO2(111). For both surfaces, the O vacancy forms more readily than in the bulk. An interesting oscillatory behavior is found for the vacancy formation energy in the upper three layers of CeO2(111). Analysis of the reduced surfaces suggests that the additional charge resulting from the formation of the oxygen vacancies is localized in the first three layers of the surface. Furthermore, they are not only trapped in the 4f states of cerium.     

Characterizing TiO2(110) surface states by their work function

The unreconstructed TiO(2)(110) surface is prepared in well-defined states having different characteristic stoichiometries, namely reduced (r-TiO(2), 6 to 9% surface vacancies), hydroxylated (h-TiO(2), vacancies filled with OH), oxygen covered (ox-TiO(2), oxygen adatoms on a stoichiometric surface) and quasi-stoichiometric (qs-TiO(2), a stoichiometric surface with very few defects). The electronic structure and work function of these surfaces and transition states between them are investigated by ultraviolet photoelectron spectroscopy (UPS) and metastable impact electron spectroscopy (MIES). The character of the surface is associated with a specific value of the work function that varies from 4.9 eV for h-TiO(2), 5.2 eV for r-TiO(2), 5.35 eV for ox-TiO(2) to 5.5 eV for qs-TiO(2). We establish the method for an unambiguous characterization of TiO(2)(110) surface states solely based on the secondary electron emission characteristics. This is facilitated by analysing a weak electron emission below the nominal work function energy. The emission in the low energy cut-off region appears correlated with band gap emission found in UPS spectra and is attributed to localised electron emission through Ti(3+)(3d) states.     

First-principles local density approximation (generalized gradient approximation) +U study of catalytic CenOm clusters: U value differs from bulk

Ceria possesses strong catalytic properties for CONO(x) removal and H(2) production. Clusters often show more intriguing functionalities than their bulk counterparts. Here, the geometric and electronic structures of Ce(n)O(m) (n=1-4,m=2n-1,2n) clusters are studied for the first time using the projected augmented wave method in density functional theory with detailed assessment of the exchange-correlation functional and the Hubbard parameter U. We note that the U value strongly affects the electronic structures of the oxygen-deficient Ce(n)O(2n-1) clusters, though less so on the stoichiometric Ce(n)O(2n). Furthermore, the local density approximation (LDA)+U method is more accurate than the generalized gradient approximation+U in describing the localization of the 4f electrons of the Ce(n)O(m) clusters. The calculated vibration frequency of the CeO molecule with the LDA+U (U=4 eV) is 818.4 cm(-1), in close agreement with experimental values of 820-825 cm(-1) for the low lying states. Different optimal U values were noted for the ceria cluster (4 eV) and its bulk (6 eV), due to quantum-size and geometric effects. The largely reduced formation energy of an oxygen vacancy indicates that the catalytic effect of the Ce(n)O(m) clusters are far greater than bulk CeO(2).     

Electronic and atomistic structures of clean and reduced ceria surfaces

The atomistic and electronic structures of oxygen vacancies on the (111) and (110) surfaces of ceria are studied by means of periodic density functional calculations. The removal of a neutral surface oxygen atom leaves back two excess electrons that are shown to localize on two cerium ions neighboring the defect. The resulting change of valency of these Ce ions (Ce4+ --> Ce3+) originates from populating tightly bound Ce 4f states and is modeled by adding a Hubbard U term to the traditional energy functionals. The calculated atomistic and electronic structures of the defect-free and reduced surfaces are shown to agree with spectroscopic and microscopic measurements. The preferential defect segregation and the different chemical reactivity of the (111) and (110) surfaces are discussed in terms of energetics and features in the electronic structure.     

NbC固体表面能带和电子结构的理论研究

采用DFT/BLYP方法对NbC(001)和(111)面的电子结构进行研究。计算结果表明,对于NbC(001)表面,其表面态主要集中于费米能级(EF)下方约4.5eV附近区域,并以表面Nb原子和C原子为主要成分。O2分子在该表面吸附时,趋向于吸附在表面Nb原子上。对于NbC(111)表面,其表面态集中在EF下方0.02.0eV区域,靠近EF的态具有较高的表面活性,其主要成分为表面Nb原子的4dxz/dyz成分。上述结论与光电子能谱实验结果基本一致;但由于金属原子d电子数的差异导致NbC(111)表面态成分与类似的TiC化合物并不相同。     

High-pressure STM of the interaction of oxygen with Ag(111)

To identify surface phases that could play a role for the epoxidation of ethylene on Ag catalysts we have studied the interaction of Ag(111) with O(2) at elevated pressures. Experiments were performed using high-pressure scanning tunneling microscopy (STM) at temperatures between 450 and 480 K and O(2) pressures in the mbar range. Below p(O(2)) approximately 1 mbar the surface largely showed the structure of bare Ag(111). At p(O(2)) above approximately 1 mbar the (4 x 4)O structure and the closely related (4 x 5 radical 3)rect structure were observed. The findings confirm theoretical predictions that the (4 x 4)O structure is thermodynamically stable at the oxygen partial pressure of the industrial ethylene oxide synthesis. However, in other experiments only a rough, disordered structure was observed. The difference is caused by the chemical state of the STM cell that depends on the pretreatment and on previous experiments. The surface was further analyzed by X-ray photoelectron spectroscopy (XPS). Although these measurements were performed after sample transfer to ultra-high vacuum (UHV), so that the surface composition was modified, the two surface states could still be identified by the presence of carbonate or a carbonaceous species, and by the absence or presence of a high-binding energy oxygen species, respectively. It turns out that the (4 x 4)O structure only forms under extremely clean conditions, indicating that the (4 x 4)O phase and similar oxygen-induced reconstructions of the Ag(111) surface are chemically unstable. Chemical reactions at the inner surfaces of the STM cell also complicate the detection of the catalytic formation of ethylene oxide.     

甲烷在清洁 Pd(111) 及氧改性的 Pd(111) 表面解离的密度泛函理论研究

 采用广义梯度近似 (GGA) 的密度泛函理论 (DFT) 并结合平板模型, 研究了 CH4 在清洁 Pd(111) 及 O 改性的 Pd(111) 表面发生 C朒 键断裂的反应历程. 优化了裂解过程中反应物、过渡态和产物的几何构型, 获得了反应路径上各物种的吸附能及反应的活化能. 结果表明, CH4 采用一个 H 原子指向表面的构型在 Pd(111) 表面的顶位吸附, CH3 的最稳定的吸附位置为顶位, OH, O 和 H 的最稳定吸附位置均为面心立方. CH4 在清洁 Pd(111) 表面裂解的活化能为 0.97 eV, 低于它在 O 原子改性 (O 没有参与反应) 的 Pd(111) 表面的活化能 1.42 eV, 说明表面氧原子抑制了 CH4 中 C朒 键的断裂. 当亚表面 O 原子和表面 O 原子 (O 参与反应) 共同存在时, C朒 键断裂的活化能为 0.72 eV, 低于只有表层氧存在时的活化能 (1.43 eV), 说明亚表面的 O 原子对 CH4 分子的活化具有促进作用. CH4 在 O 原子改性的 Pd(111) 表面裂解生成 CH3 和 H, 以及生成 CH3 和 OH 的反应活化能分别为 1.42 和 1.43 eV, 说明 CH4 在 O 原子改性的 Pd(111) 表面发生这两种反应的难易程度相当.     

A first principles study of oxygen reduction reaction on a Pt(111) surface modified by a subsurface transition metal M (M = Ni, Co, or Fe)

We have performed first-principle density functional theory calculations to investigate how a subsurface transition metal M (M = Ni, Co, or Fe) affects the energetics and mechanisms of oxygen reduction reaction (ORR) on the outermost Pt mono-surface layer of Pt/M(111) surfaces. In this work, we found that the subsurface Ni, Co, and Fe could down-shift the d-band center of the Pt surface layer and thus weaken the binding of chemical species to the Pt/M(111) surface. Moreover, the subsurface Ni, Co, and Fe could modify the heat of reaction and activation energy of various elementary reactions of ORR on these Pt/M(111) surfaces. Our DFT results revealed that, due to the influence of the subsurface Ni, Co, and Fe, ORR would adopt a hydrogen peroxide dissociation mechanism with an activation energy of 0.15 eV on Pt/Ni(111), 0.17 eV on Pt/Co(111), and 0.16 eV on Pt/Fe(111) surface, respectively, for their rate-determining O2 protonation reaction. In contrast, ORR would follow a peroxyl dissociation mechanism on a pure Pt(111) surface with an activation energy of 0.79 eV for its rate-determining O protonation reaction. Thus, our theoretical study explained why the subsurface Ni, Co, and Fe could lead to multi-fold enhancement in catalytic activity for ORR on the Pt mono-surface layer of Pt/M(111) surfaces.     

Size-dependent surface reactions of Ag nanoparticles supported on highly oriented pyrolytic graphite

Various sizes of Ag particles were grown on highly oriented pyrolytic graphite (HOPG) surfaces, which had previously been modified with nanopits to act as anchoring sites. Surface reactions of O2, CHCl3, and CCl4 on the Ag particles and bulk Ag(111) surfaces were studied by X-ray photoelectron spectroscopy (XPS), and it has been shown that size dependence of O2 and CHCl3 reactions on Ag differs from that of CCl4. Weak reactions of O2 and CHCl3 were observed on the bulk Ag(111) surfaces, while strong reactions occur on Ag particles with medium Ag coverage, suggesting that the reactions are controlled by the number of surface defect sites. On the contrary, the dissociation of CCl4 is mainly determined by the exposed Ag facet area, mainly Ag(111) facet, and strong dissociation reaction happens on the bulk Ag(111) surface. The results suggest that the size effects, which are often discussed in heterogeneous catalysis, are strongly dependent on the reaction mechanism.     

Probing the 4f states of ceria by tunneling spectroscopy

Low-temperature scanning tunneling microscopy and spectroscopy have been employed to analyze the local electronic structure of the (111) surface of a ceria thin film grown on Ru(0001). On pristine, defect-free oxide terraces, the empty 4f states of Ce(4+) ions appear as the only spectral feature inside the 6 eV oxide band gap. In contrast, occupied states are detected between -1.0 and -1.5 eV below E(Fermi) in conductance spectra of different point and line defects, such as surface oxygen vacancies, grain boundaries and step edges. They are assigned to partially filled 4f states localized at the Ce(3+) ions. The presence of excess electrons indicates the oxygen-deficient nature of the direct oxide environment. The f state spectroscopy with the STM allows us to probe the spatial distribution of Ce(3+) ions in the ceria surface, providing unique insight into the local reduction state of this chemically important material system.     

DFT+U study of defects in bulk rutile TiO(2)

We present a systematic study of electronic gap states in defected titania using our implementation of the Hubbard-U approximation in the grid-based projector-augmented wave density functional theory code, GPAW. The defects considered are Ti interstitials, O vacancies, and H dopants in the rutile phase of bulk titanium dioxide. We find that by applying a sufficiently large value for the Hubbard-U parameter of the Ti 3d states, the excess electrons localize spatially at the Ti sites and appear as states in the band gap. At U=2.5?eV, the position in energy of these gap states are in fair agreement with the experimental observations. In calculations with several excess electrons and U=2.5?eV, all of these end up in gap states that are spatially localized around specific Ti atoms, thus effectively creating one Ti(3+) ion per excess electron. An important result of this investigation is that regardless of which structural defect is the origin of the gap states, at U=2.5?eV, these states are found to have their mean energies within a few hundredths of an eV from 0.94 eV below the conduction band minimum.     

First principles based mean field model for oxygen reduction reaction

A first principles-based mean field model was developed for the oxygen reduction reaction (ORR) taking account of the coverage- and material-dependent reversible potentials of the elementary steps. This model was applied to the simulation of single crystal surfaces of Pt, Pt alloy and Pt core-shell catalysts under Ar and O(2) atmospheres. The results are consistent with those shown by past experimental and theoretical studies on surface coverages under Ar atmosphere, the shape of the current-voltage curve for the ORR on Pt(111) and the material-dependence of the ORR activity. This model suggests that the oxygen associative pathway including HO(2)(ads) formation is the main pathway on Pt(111), and that the rate determining step (RDS) is the removal step of O(ads) on Pt(111). This RDS is accelerated on several highly active Pt alloys and core-shell surfaces, and this acceleration decreases the reaction intermediate O(ads). The increase in the partial pressure of O(2)(g) increases the surface coverage with O(ads) and OH(ads), and this coverage increase reduces the apparent reaction order with respect to the partial pressure to less than unity. This model shows details on how the reaction pathway, RDS, surface coverages, Tafel slope, reaction order and material-dependent activity are interrelated.     

Comparative study of water dissociation on Rh(111) and Ni(111) studied with first principles calculations

The dissociation and formation of water on the Rh(111) and Ni(111) surfaces have been studied using density functional theory with generalized gradient approximation and ultrasoft pseudopotentials. Calculations have been performed on 2x2 surface unit cells, corresponding to coverages of 0.25 ML, with spot checks on 3x3 surface unit cells (0.11 ML). On both surfaces, the authors find that water adsorbs flat on top of a surface atom, with binding energies of 0.35 and 0.25 eV, respectively, on Rh(111) and Ni(111), and is free to rotate in the surface plane. Barriers of 0.92 and 0.89 eV have to be overcome to dissociate the molecule into OH and H on the Rh(111) and Ni(111) surfaces, respectively. Further barriers of 1.03 and 0.97 eV need to be overcome to dissociate OH into O and H. The barriers for the formation of the OH molecule from isolated adsorbed O and H are found to be 1.1 and 1.3 eV, and the barriers for the formation of the water molecule from isolated adsorbed OH and H are 0.82 and 1.05 eV on the two surfaces. These barriers are found to vary very little as coverage is changed from 0.25 to 0.11 ML. The authors have also studied the dissociation of OH in the presence of coadsorbed H or O. The presence of a coadsorbed H atom only weakly affects the energy barriers, but the effect of O is significant, changing the dissociation barrier from 1.03 to 1.37 and 1.15 eV at 0.25 or 0.11 ML coverage on the Rh(111) surface. Finally, the authors have studied the dissociation of water in the presence of one O atom on Rh(111), at 0.11 ML coverage, and the authors find a barrier of 0.56 eV to dissociate the molecule into OH+OH.     

Probing water interactions and vacancy production on gadolinia-doped ceria surfaces using electron stimulated desorption

Polycrystalline gadolinia-doped ceria (GDC) surfaces were studied using low-energy (5-400 eV) electron stimulated desorption (ESD). H(+), O(+), and H(3)O(+) were the primary cationic desorption products with H(+) as the dominant channel. H(+), H(3)O(+), and O(+) have a 22 eV threshold followed by a yield change around 40 eV. H(+) also has an additional yield change approximately 75 eV and O(+) has an additional change approximately 150 eV. The O(+) ESD yield change approximately 150 eV may indicate bond breaking of Gd-O and the involvement of oxygen vacancies. The H(+) and H(3)O(+) threshold data collectively indicate the presence of hydroxyl groups and chemisorbed water molecules on the GDC surfaces. ESD temperature dependence measurements show that the interaction of water with GDC surface defect sites, mainly oxygen vacancies, influences the desorption of H(+), O(+), and H(3)O(+). The temperature dependence of the O(+) ESD at 400 eV incident electron energy yields a 0.21 eV activation energy. This is close to the energy needed for oxygen vacancy production next to a pair of Ce(3+) on a CeO(2) surface. These results may indicate a correlation between the O(+) ESD yield and oxygen vacancy density on GDC surfaces and a potential correlation of O(+) ESD and GDC ionic conductivity.     

A theoretic insight into the catalytic activity promotion of CeO2 surfaces by Mn doping

In this paper, we investigated the primary reduction and oxygen replenishing processes over Mn substitutionally doped CeO(2)(111) surfaces by density functional theory with the on-site Coulomb correction (DFT + U). The results indicated that Mn doping could make the surface much more reducible and the adsorbed O(2) could be effectively activated to form superoxo (O(2)(-)) and/or peroxo species (O(2)(2-)). The Mn doping induced the Mn 3d-O 2p gap state instead of Ce 4f acting as an electrons acceptor and donor during the first oxygen vacancy formation and O(2) replenishing, which helped to lower the formation energy of the first and second oxygen vacancies to -0.46 eV and 1.40 eV, respectively. In contrast, the formation energy of a single oxygen vacancy in the pure ceria surface was 2.08 eV and only peroxo species were identified as the O(2) molecule adsorbed. Our work provides a theoretical and electronic insight into the catalytic redox processes of Mn doped ceria surfaces, which may help to understand the enhanced catalytic performances of MnO(x)-CeO(2) oxides, as reported in previous experimental works.     

Orientation dependent oxygen exchange kinetics on single crystal SrTiO(3) surfaces

The perovskite SrTiO(3) is arguably one of the most important oxide systems in condensed matter research. In this study, we report measurement of the orientation dependence of oxygen exchange on SrTiO(3) single crystal surfaces by dynamic conductivity measurements under electrochemical perturbations. Activation energy for electrical conduction in the 923-1223 K range at an oxygen partial pressure of ～10(-11) Pa of (100), (111), and (110) single crystals was found to be 2.6 eV, 2.7 eV, and 3.1 eV, respectively. The equilibration kinetics show profound dependence on the surface orientation and are modelled using a heterogeneous relaxation process. All surfaces show similar cationic sub-lattice limited rate behavior with (111), (100), and (110) having the fastest, intermediate, and slowest rates, respectively. We discuss the orientation dependence and its relation to local atomic structure in light of previous experimental and theoretical studies.     

Shape evolution of Cu2O nanostructures via kinetic and thermodynamic controlled growth

We report the shape evolution process of Cu(2)O nanocrystals upon slow oxidation of Cu under ambient conditions, yielding novel hexagonal and triangular platelike morphologies. The shape of the obtained nanocrystals evolves from hexagonal to triangular to octahedral; the growth patterns are governed by kinetically and thermodynamically controlled growth. Preferential adsorption of I(-) on {111} planes of Cu(2)O nanoparticles induced the selective crystal growth of metastable platelike structures with {111} faces as the basal planes. On aging, the growth process appeared to shift into the thermodynamic regime and the thermodynamically stable octahedral shape is obtained. The possible growth mechanisms were investigated by varying the synthetic conditions. The band gap of Cu(2)O nanooctahedrons was determined by the classical Tauc approach to be 2.24 eV, which is blue shifted with respect to the bulk Cu(2)O value (2.17 eV). Results suggest that the slow oxidation process and use of crystallographic selective surfactants are essential for the appearance of anisotropic metastable shapes. In general, surface energy control by surfactant molecules might provide a convenient channel for tailoring nanocrystal shapes of metal oxides.     

The surface dependence of CO adsorption on Ceria

An understanding of the interaction between ceria and environmentally sensitive molecules is vital for developing its role in catalysis. We present the structure and energetics of CO adsorbed onto stoichiometric (111), (110), and (100) surfaces of ceria from first principles density functional theory corrected for on-site Coulomb interactions, DFT+U. DFT+U is applied because it can describe consistently the properties of both the stoichiometric and reduced surfaces. Our major finding is that the interaction is strongly surface dependent, consistent with experiment. Upon interaction of CO with the (111) surface, weak binding is found, with little perturbation to the surface or the molecule. For the (110) and (100) surfaces, the most stable adsorbate is that in which the CO molecule bridges two oxygen atoms and pulls these atoms out of their lattice sites, with formation of a (CO(3)) species. This results in a strong modification to the surface structure, consistent with that resulting from mild reduction. The electronic structure also demonstrates reduction of the ceria surface and consequent localization of charge on cerium atoms neighboring the vacancy sites. The surface-bound (CO(3)) species is identified as a carbonate, (CO(3))(2-) group, which is formed along with two reduced surface Ce(III) ions, in good agreement with experimental infrared data. These results provide a detailed investigation of the interactions involved in the adsorption of CO on ceria surfaces, allowing a rationalization of experimental findings and demonstrate further the applicability of the DFT+U approach to the study of systems in which reduced ceria surfaces play a role.     

A theoretical study of surface reduction mechanisms of CeO(2)(111) and (110) by H(2).

Reaction mechanisms for the interactions between CeO(2)(111) and (110) surfaces are investigated using periodic density functional theory (DFT) calculations. Both standard DFT and DFT+U calculations to examine the effect of the localization of Ce 4f states on the redox chemistry of H(2)-CeO(2) interactions are described. For mechanistic studies, molecular and dissociative local minima are initially located by placing an H(2) molecule at various active sites of the CeO(2) surfaces. The binding energies of physisorbed species optimized using the DFT and DFT+U methods are very weak. The dissociative adsorption reactions producing hydroxylated surfaces are all exothermic; exothermicities at the DFT level range from 4.1 kcal mol(-1) for the (111) to 26.5 kcal mol(-1) for the (110) surface, while those at the DFT+U level are between 65.0 kcal mol(-1) for the (111) and 81.8 kcal mol(-1) for the (110) surface. Predicted vibrational frequencies of adsorbed OH and H(2)O species on the surfaces are in line with available experimental and theoretical results. Potential energy profiles are constructed by connecting molecularly adsorbed and dissociatively adsorbed intermediates on each CeO(2) surface with tight transition states using the nudged elastic band (NEB) method. It is found that the U correction method plays a significant role in energetics, especially for the intermediates of the exit channels and products that are partially reduced. The surface reduction reaction on CeO(2)(110) is energetically much more favorable. Accordingly, oxygen vacancies are more easily formed on the (110) surface than on the (111) surface.     

氧化铜表面臭氧分解路径及表面氧物种生成机理研究

基于密度泛函理论（DFT）计算研究了O3在完整和具有氧空位的CuO(111)表面吸附的吸附位、吸附结构、吸附能和电子转移情况,比较了O3在完整表面和具有氧空位的表面分解的路径和能垒,分析了氧空位和表面吸附氧的生成机理。结果表明,在完整CuO表面,O3分子通过化学吸附或物理吸附表面结合,吸附能最高为-1.22eV（构型bri(2)）。O3在具有氧空位的CuO表面均为化学吸附,吸附能最高为-2.95eV（构型ovbri(3)）,显著高于完整表面的吸附能。O3吸附后,Cu吸附位的电荷密度减小,O3中的O原子附近的电荷密度显著增强,电荷从CuO表面转移到O3,并形成Cu-O离子键。O3分解后形成了超氧物种,提高了表面的氧化活性。在完整表面,以构型bri(2)为起始构型的路径反应能垒最低,为0.52eV;O2*在完整表面的脱附所需要的最低能量为0.42eV,形成氧空位的O2*脱附能为2.06eV。在具有氧空位的表面,O3分解的反应能垒为0.30eV（构型ovbri(1)）和0.12eV（构型ovbri(3)）,均低于完整表面的反应能垒;分解形成的O2*的最低脱附能也低于完整表面,为0.27eV。可见,氧空位的形成提高了吸附能,降低了反应能垒,使O3分子更容易吸附在CuO表面,并加快了O3的催化分解。