Density functional study of hypophosphite adsorption on Ni (1 1 1) and Cu (1 1 1) surfaces

Surface structures and electronic properties of hypophosphite, H₂PO₂⁻, molecularly adsorbed on Ni(1 1 1) and Cu(1 1 1) surfaces are investigated in this work by density functional theory at B3LYP/6-31++g(d, p) level. We employ a four-metal-atom cluster as the simplified model for the surface and have fully optimized the geometry and orientation of H₂PO₂⁻ on the metal cluster. Six stable orientations have been discovered on both Ni (1 1 1) and Cu (1 1 1) surfaces. The most stable orientation of H₂PO₂⁻ was found to have its two oxygen atoms interact the surface with two PO bonds pointing downward. Results of the Mulliken population analysis showed that the back donation from 3d orbitals of the transition metal substrate to the unfilled 3d orbital of the phosphorus atom in H₂PO₂⁻ and 4s orbital's acceptance of electron donation from one lone pair of the oxygen atom in H₂PO₂⁻ play very important roles in the H₂PO₂⁻ adsorption on the transition metals. The averaged electron configuration of Ni in Ni₄ cluster is 4s^0.634p^0.023d^9.35 and that of Cu in Cu₄ cluster is 4s^1.004p^0.033d^9.97. Because of this subtle difference of electron configuration, the adsorption energy is larger on the Ni surface than on the Cu surface. The amount of charge transfers due to above two donations is larger from H₂PO₂⁻ to the Ni surface than to the Cu surface, leading to a more positively charged P atom in Ni_nH₂PO₂⁻ than in Cu_nH₂PO₂⁻. These results indicate that the phosphorus atom in Ni_nH₂PO₂⁻ complex is easier to be attacked by a nucleophile such as OH⁻ and subsequent oxidation of H₂PO₂⁻ can take place more favorably on Ni substrate than on Cu substrate.     

Density functional and dynamics study of the dissociative adsorption of hydrogen on Mg (0 0 0 1) surface

A first principles study is performed to investigate the adsorption characteristics of hydrogen on magnesium surface. Substitutional and on-surface adsorption energies are calculated for Mg (0 0 0 1) surface alloyed with the selected elements. To further analyze the hydrogen-magnesium interaction, first principles molecular dynamics method is used which simulates the behavior of H₂ at the surface. Also, charge density differences of substitutionally doped surface configurations were illustrated. Accordingly, Mo and Ni are among the elements yielding lower adsorption energies, which are found to be −9.2626 and −5.2995 eV for substitutionally alloyed surfaces, respectively. In light of the dynamic calculations, Co as an alloying element is found to have a splitting effect on H₂ in 50 fs, where the first hydrogen atom is taken inside the Mg substrate right after the decomposition and the other after 1300 fs. An interesting remark is that, elements which acquire higher chances of adsorption are also seen to be competent at dissociating the hydrogen molecule. Furthermore, charge density distributions support the results of molecular dynamics simulations, by verifying the distinguished effects of most of the 3d and 4d transition metals.     

A density-functional study on the atomic geometry and adsorption of the Cu(1 0 0) c(2 × 2)/N Surface

In this work we have performed total-energy calculations on the geometric structure and adsorption properties of Cu(1 0 0) c(2 × 2)/N surface by using the density-functional theory and the projector-augmented wave method. It is concluded that nitrogen atom was adsorbed on a FFH site with a vertical distance of 0.2 Å towards from surface Cu layer. The bond length of the shortest Cu-N bonding is calculated to be 1.83 Å. Geometry optimization calculations exclude out the possibilities of adsorbate induced reconstruction mode suggested by Driver and Woodruff and the atop structural model. The calculated workfunction for this absorbate-adsorbent system is 4.63 eV which is quite close to that of a clean Cu(1 0 0) surface. The total-energy calculations showed that the average adsorption energy per nitrogen in the case of Cu(1 0 0) c(2 × 2)-N is about 4.88 eV with respect to an isolated N atom. The absorption of nitrogen on Cu(1 0 0) surface yields the hybridization between surface Cu atoms and N, and generates the localized surface states at −1.0 eV relative to Fermi energy E_F. The stretch mode of the adsorbed nitrogen at FFH site is about 30.8 meV. The present study provides a strong criterion to account for the local surface geometry in Cu(1 0 0) c(2 × 2)/N surface.     

Driving force for the WO3(0 0 1) surface relaxation

The optimized structure of the WO₃(0 0 1) surface with various types of termination ((1 × 1)O, (1 × 1)WO₂, and c(2 × 2)O) has been simulated using density functional theory with the Perdew-Wang 91 gradient corrected exchange-correlation functional. While the energy of bulk WO₃ depends weakly on the distortions and tilting of the WO₆ octahedra, relaxation of the (0 0 1) surface results in a significant decrease of surface energy (from 10.2 × 10⁻² eV/Ų for the cubic ReO₃-like, c(2 × 2)O-terminated surface to 2.2 × 10⁻² eV/Ų for the relaxed surface). This feature illustrates a potential role of surface relaxation in formation of crystalline nano-size clusters of WO₃. The surface relaxation is accompanied by a dramatic redistribution of the density of states near the Fermi level, in particular a transformation of surface electronic states. This redistribution is responsible for the decrease of electronic energy and therefore is suggested to be the driving force for surface relaxation of the WO₃(0 0 1) surface and, presumably, similar surfaces of other transition metal oxides.     

First principles study of H2S adsorption and dissociation on Fe(1 1 0)

We report first principles density functional theory (DFT) results of H₂S and HS adsorption and dissociation on the Fe(1 1 0) surface. We investigate the site preference of H₂S, HS, and S on Fe(1 1 0). H₂S is found to weakly adsorb on either the short bridge (SB) or long bridge (LB) site of Fe(1 1 0), with a binding energy of no more than 0.50 eV. The diffusion barrier from the LB site to the SB site is found to be small (∼0.10 eV). By contrast to H₂S, HS is predicted to be strongly chemisorbed on Fe(1 1 0), with the S atom in the LB site and the HS bond oriented perpendicular to the surface. Isolated S atoms also are predicted to bind strongly to the LB sites of Fe(1 1 0), where the SB is found to be a transition state for S surface hopping between neighboring LB sites. The minimum energy paths for H₂S and HS dehydrogenation involve rotating an H atom towards a nearby surface Fe atom, with the S-H bonds breaking on the top of one Fe atom. The barrier to break the first S-H bond in H₂S is low at 0.10 eV, and breaking the second S-H bond is barrierless, suggesting deposition of S on Fe(1 1 0) via H₂S is kinetically and thermodynamically facile.     

Phosphine and tertiarybutylphosphine adsorption on the indium-rich InP (0 0 1)-(2 × 4) surface

Phosphine and tertiarybutylphosphine adsorption on the indium-rich InP (0 0 1)-(2 × 4) surface at 25 °C have been studied by internal reflection infrared spectroscopy, X-ray photoelectron spectroscopy, and low energy electron diffraction. Both molecules form a dative bond to the empty dangling bonds on the In-P heterodimers and the second-layer In-In dimers and vibrate symmetrically at 2319 (2315) and 2285 (2281) cm⁻¹ and asymmetrically at 2339 (2339) and 2327 (2323) cm⁻¹. A fraction of these species dissociate into adsorbed PH₂ with the hydrogen and tertiarybutyl ligands transferring to nearby phosphorus sites. The calculated energy barriers for desorption (<11 kcal/mol) of these molecules is less than that for dissociation (>17 kcal/mol) and explains their low sticking probabilities at elevated temperatures under InP growth conditions.     

Formation and characterization of the Cu2O overlayer on Zn-terminated ZnO(0 0 0 1)

Low-energy electron diffraction, X-ray photoelectron spectroscopy and synchrotron-radiation-excited angle-resolved photoelectron spectroscopy have been used to characterize Cu-oxide overlayers on the Zn-terminated ZnO(0 0 0 1) surface. Deposition of Cu on the ZnO(0 0 0 1)-Zn surface results in the formation of Cu clusters with (1 1 1) top terraces. Oxidation of these clusters by annealing at 650 K in O₂ atmosphere (1.3 × 10⁻⁴ Pa) leads to an ordered Cu₂O overlayer with (1 1 1) orientation. Good crystallinity of the Cu₂O(1 1 1) overlayer is proved by energy dispersion of one of Cu₂O valence bands. The Cu₂O(1 1 1) film exhibits a strong p-type semiconducting nature with the valence band maximum (VBM) of 0.1 eV below the Fermi level. The VBM of ZnO at the Cu₂O(1 1 1)/ZnO(0 0 0 1)-Zn interface is estimated to be 2.4 eV, yielding the valence-band offset of 2.3 eV.     

Theoretical studies on the adsorption and decomposition of H2O on Pd(1 1 1) surface

To provide information about the chemistry of water on Pd surfaces, we performed density functional slab model studies on water adsorption and decomposition at Pd(1 1 1) surface. We located transition states of a series of elementary steps and calculated activation energies and rate constants with and without quantum tunneling effect included. Water was found to weakly bind to the Pd surface. Co-adsorbed species OH and O that are derivable from H₂O stabilize the adsorbed water molecules via formation of hydrogen bonds. On the clean surface, the favorable sites are top and bridge for H₂O and OH, respectively. Calculated kinetic parameters indicate that dehydrogenation of water is unlikely on the clean regular Pd(1 1 1) surface. The barrier for the hydrogen abstraction of H₂O at the OH covered surface is approximately 0.2-0.3 eV higher than the value at the clean surface. Similar trend is computed for the hydroxyl group dissociation at H₂O or O covered surfaces. In contrast, the O-H bond breaking of water on oxygen covered Pd surfaces, H₂O_ad + O_ad → 2OH_ad, is predicted to be likely with a barrier of ∼0.3 eV. The reverse reaction, 2OH_ad → H₂O_ad + O_ad, is also found to be very feasible with a barrier of ∼0.1 eV. These results show that on oxygen-covered surfaces production of hydroxyl species is highly likely, supporting previous experimental findings.     

Cu adsorption on pyrite (1 0 0): Ab initio and spectroscopic studies

Surface optimised S 2p photoelectron spectra show that both surface S²⁻ monomers and (S-S)²⁻ dimers are present at pyrite (1 0 0) fracture surfaces. In order to determine which sulfur species are involved in Cu adsorption, fresh pyrite surfaces were exposed to Cu²⁺ in solution. The S 2p spectra suggest that both types of S surface species are involved in the mechanism of Cu adsorption (activation). Ab initio density functional theory was used to model Cu adsorbed onto pyrite (1 0 0) to support the interpretation of the spectroscopy. Mulliken population analysis confirms the charge distribution suggested by the core line shifts as observed in the photoelectron spectra. The ab initio calculations were consistent with a two-coordinate bond between Cu(I), a surface S monomer and a surface S dimer.     

Theoretical study on adsorption of Au and hydrated Au cations on clean Si(1 1 1) surface

To model the adsorption of Au⁺ cation in aqueous solution on the semiconductor surface, the interactions of Au⁺ and hydrated Au⁺ cations with clean Si(1 1 1) surface were investigated by using hybrid density functional theory (B3LYP) and Møller-Plesset second-order perturbation (MP2) methods. Si(1 1 1) surface was described with Si₇H₁₁, Si₁₁H₁₇ and Si₂₂H₂₁ clusters. The effect of the basis set superposition error (BSSE) was taken into account by applying the counterpoise correction. The calculated results indicated that the binding energies between hydrated Au⁺ cations and clean Si(1 1 1) surface are large, suggesting a strong interaction between hydrated Au⁺ cations and the semiconductor surface. The bonding nature of the chemical adsorption of Au⁺ to Si surface can be classified as partial covalent as well as ionic bonding. As the number of water molecules increases, the water molecules form hydrogen bond network with one another and only one water molecule binds directly to the Au⁺ cation. The Au⁺ cation in aqueous solution will safely attach to the clean Si(1 1 1) surface.     

Probing the properties of the (1 1 1) and (1 0 0) surfaces of LaB6 through infrared spectroscopy of adsorbed CO

The adsorption of carbon monoxide on the LaB₆(1 0 0) and LaB₆(1 1 1) surfaces was studied experimentally with the techniques of reflection absorption infrared spectroscopy and X-ray photoelectron spectroscopy. The interaction of CO with the two surfaces was also studied with density functional theory. Both surfaces adsorb CO molecularly at low temperatures but in markedly different forms. On the LaB₆(1 1 1) surface CO initially adsorbs at 90 K in a form that yields a CO stretching mode at 1502-1512 cm⁻¹. With gentle annealing to 120 K, the CO switches to a bonding environment characterized by multiple CO stretch values from 1980 to 2080 cm⁻¹, assigned to one, two, or three CO molecules terminally bonded to the B atoms of a triangular B₃ unit at the (1 1 1) surface. In contrast, on the LaB₆(1 0 0) surface only a single CO stretch is observed at 2094 cm⁻¹, which is assigned to an atop CO molecule bonded to a La atom. The maximum intensity of the CO stretch vibration on the (1 0 0) surface is higher than on the (1 1 1) surface by a factor of 5. This difference is related to the different orientations of the CO molecules on the two surfaces and to reduced screening of the CO dynamic dipole moment on the (1 0 0) surface, where the bonding occurs further from the surface plane. On LaB₆(1 0 0), XPS measurements indicate that CO dissociates on the surface at temperatures above 400 K.     

Initial stages of oxidation of Cu(1 1 1)

Several surface analysis techniques were combined to study the initial stages of oxidation of Cu(1 1 1) surfaces exposed to O₂ at low pressure (<5 × 10⁻⁶ mbar) and room temperature. Scanning tunneling microscopy (STM) results show that the reactivity is governed by the restructuring of the Cu(1 1 1) surface. On the terraces, oxygen dissociative adsorption leads to the formation of isolated O adatoms and clusters weakly bound to the surface. The O adatoms are located in the fcc threefold hollow sites of the unrestructured terraces. Friedel oscillations with an amplitude lower than 5 pm have been measured around the adatoms. At step edges, surface restructuring is initiated and leads to the nucleation and growth of a two-dimensional disordered layer of oxide precursor. The electronic structure of this oxide layer is characterised by a band gap measured by scanning tunneling spectroscopy to be ∼1.5 eV wide. The growth of the oxide islands progresses by consumption of the upper metal terraces to form triangular indents. The extraction of the Cu atoms at this interface generates a preferential orientation of the interface along the close-packed directions of the metal. A second growth front corresponds to the step edges of the oxide islands and progresses above the lower metal terraces. This is where the excess Cu atoms extracted at the first growth front are incorporated. STM shows that the growing disordered oxide layer consists of units of hexagonal structure with a first nearest neighbour distance characteristic of a relaxed Cu-Cu distance (∼0.3 nm), consistent with local Cu₂O(1 1 1)-like elements. Exposure at 300 °C is necessary to form an ordered two-dimensional layer of oxide precursor. It forms the so-called “29” superstructure assigned to a periodic distorted Cu₂O(1 1 1)-like structure.     

Adsorption of CO, CO2 and NO molecules on a BaTiO3 (0 0 1) surface

Since the development of Scanning Tunnelling Microscopy (STM) technique, considerable attention has been devoted to various molecules adsorbed on various surfaces. Also, a new concept emerged with molecules on surfaces considered as nano machines by themselves. In this context, a thorough knowledge of surfaces and adsorbed molecules at an atomic scale are thus particularly invaluable. The present work describes the first Density Functional Theory (DFT) study of adsorption of CO, CO₂ and NO molecules on a BaTiO₃ surface following a first preliminary calculation of O and O₂ adsorption on the same surface. In the previously considered work, we found that a (0 0 1) surface with BaO termination is more stable than the one with TiO₂-termination. Consequently, we extended our study to CO, CO₂ and NO molecules adsorbed on a (0 0 1) surface with BaO termination. The present calculation was performed on a (1 × 1) cell with one monolayer of adsorbed molecules. Especially, a series of cases implying CO molecules adsorbed in various geometrical configurations has been examined. The corresponding adsorption energy varies in the range of −0.17 to −0.10 eV. The adsorption energy of a CO₂ molecule directly located above an O surface atom (called O_s) is of the order of −0.18 eV. The O-C distance length is then 1.24 Å and the O-C-O and O-C-O_s angles are 134.0° and 113.0°, respectively. For NO adsorption, the most important induced structural changes are the followings: (i) the N-O bond is broken when a NO molecule is absorbed on a Ba-O_s bridge site. In that case, N and O atoms are located above an O and a Ba surface atom, respectively, whereas the O-Ba-O_s and N-O_s-Ba angles are 106.5° and 63.0°, respectively. The N-O distance is as large as 2.58 Å and the adsorption energy is as much as −2.28 eV. (ii) In the second stable position, the NO molecule has its N atom adsorbed above an O_s atom, the N-O axis being tilted toward the Ba atom. The N-O_s-Ba angle is then 41.1° while the adsorption energy is only −0.10 eV. At last, the local densities of states around C, O as well as N atoms of the considered adsorbed molecules have also been discussed.     

Sum frequency generation and density functional studies of CO-H interaction and hydrogen bulk dissolution on Pd(1 1 1)

CO-H interaction and H bulk dissolution on Pd(1 1 1) were studied by sum frequency generation (SFG) vibrational spectroscopy and density functional theory (DFT). The theoretical findings are particularly important to rationalize the experimentally observed mutual site blocking of CO and H and the effect of H dissolution on coadsorbate structures. Dissociative hydrogen adsorption on CO-precovered Pd(1 1 1) is impeded due to an activation barrier of ∼2.5 eV for a CO coverage of 0.75 ML, an effect which is maintained down to 0.33 ML CO. Preadsorbed hydrogen prevented CO adsorption at 100 K, while hydrogen was replaced from the surface by CO above 125 K. The temperature-dependent site blocking of hydrogen originates from the onset of hydrogen diffusion into the Pd bulk around 125 K, as shown by SFG and theoretical calculations using various approaches. When Pd(1 1 1) was exposed to 1:1 CO/H₂ mixtures at 100 K, on-top CO was absent in the SFG spectra although hydrogen occupies only threefold hollow sites on Pd(1 1 1). DFT attributes the absence of on-top CO to H atoms diffusing between hollow sites via bridge sites, thereby destabilizing neighboring on-top CO molecules. According to the calculations, the stretching frequency of bridge-bonded CO with a neighboring bridge-bonded hydrogen atom is redshifted by 16 cm⁻¹ when compared to bridging CO on the clean surface. Implications of the observed effects on hydrogenation reactions are discussed and compared to the C₂H₄-H coadsorption system.     

Adsorption of O2 on the CaO and SrO (1 0 0) surfaces: A first-principles study

First-principles calculations based on density functional theory (DFT) have been performed to investigate the adsorption of O₂ on the CaO and SrO (1 0 0) surfaces. The present results indicate that the bridge-top site for both the CaO and SrO (1 0 0) surfaces is the most favorable site for O₂ adsorption, with predicted adsorption energies of 1.437 eV and 1.236 eV, respectively. Detailed analysis of density of states, Mullikan population and vibration frequency are performed. The calculated results imply the possible formation of a peroxo (O₂²⁻) when O₂ is adsorbed with the bridge-top mode on both CaO and SrO (1 0 0) surfaces.     

Density functional theory study of hydrogenation mechanism in Fe-doped Mg(0 0 0 1) surface

Using density functional theory (DFT) in combination with nudged elastic band (NEB) method, the dissociative chemisorptions and diffusion processes of hydrogen on both pure and Fe-doped Mg(0 0 0 1) surfaces are studied. Firstly, the dissociation pathway of H₂ and the relative barrier were investigated. The calculated dissociation barrier (1.08 eV) of hydrogen molecule on a pure Mg(0 0 0 1) surface is in good agreement with comparable experimental and theoretical studies. For the Fe-doped Mg(0 0 0 1) surface, the activated barrier decreases to 0.101 eV due to the strong interaction between the s orbital of H and the d orbital of Fe. Then, the diffusion processes of atomic hydrogen on pure and Fe-doped Mg(0 0 0 1) are presented. The obtained diffusion barrier to the first subsurface is 0.45 eV and 0.98 eV, respectively. Finally, Chou method was used to investigate the hydrogen sorption kinetic mechanism of pure MgH₂ and Mg mixed with 5 at.% Fe atoms composites. The obtained activation energies are 0.87 ± 0.02 and 0.31 ± 0.01 eV for H₂ dissociation on the pure surface and H atom diffusion in Fe-doped Mg surfaces, respectively. It suggests that the rate-controlling step is dissociation of H₂ on the pure Mg surface while it is diffusion of H atom in the Fe-doped Mg surface. And both of fitting data are matching well with our calculation results.     

Formation and hydrogenation of p(2 × 2)-N on Pt(1 1 1)

The formation of a well-ordered p(2 × 2) overlayer of atomic nitrogen on the Pt(1 1 1) surface and its reaction with hydrogen were characterized with reflection absorption infrared spectroscopy (RAIRS), temperature programmed desorption (TPD), low energy electron diffraction (LEED), Auger electron spectroscopy (AES), and X-ray photoelectron spectroscopy (XPS). The p(2 × 2)-N overlayer is formed by exposure of ammonia to a surface at 85 K that is covered with 0.44 monolayer (ML) of molecular oxygen and then heating to 400 K. The reaction between ammonia and oxygen produces water, which desorbs below 400 K. The only desorption product observed above 400 K is molecular nitrogen, which has a peak desorption temperature of 453 K. The absence of oxygen after the 400 K anneal is confirmed with AES. Although atomic nitrogen can also be produced on the surface through the reaction of ammonia with an atomic, rather than molecular, oxygen overlayer at a saturation coverage of 0.25 ML, the yield of surface nitrogen is significantly less, as indicated by the N₂ TPD peak area. Atomic nitrogen readily reacts with hydrogen to produce the NH species, which is characterized with RAIRS by an intense and narrow (FWHM ∼ 4 cm⁻¹) peak at 3322 cm⁻¹. The areas of the H₂ TPD peak associated with NH dissociation and the XPS N 1s peak associated with the NH species indicate that not all of the surface N atoms can be converted to NH by the methods used here.     

Assessment of competing mechanisms of the abstraction of hydrogen from CH4 on Li/MgO(0 0 1)

First-principles calculations based on the density functional theory (DFT)-pseudopotential method with plane wave basis sets are used to study the energetics of H abstraction from CH₄ on the Li-doped MgO(0 0 1) surface. Experimental work has led to the proposal of two competing mechanisms for the reaction: either direct interaction of CH₄ with an O⁻ hole state bound to the Li dopant, or interaction of CH₄ with a surface F-centre whose charge state is modified by the presence of the dopant. The calculations are performed in periodically repeated slab geometry, with attention given to system-size errors. For direct interaction with the hole state, the calculations indicate that the reaction is weakly endothermic, by about 0.2 eV, but for the F-centre mechanism it is endothermic by over 1 eV. Even allowing for likely DFT errors of a few tenths of an electron volt, this gives strong evidence against the F-centre mechanism.     

Formamide reactions on rutile TiO2(0 1 1) surface

The reaction of formamide over the (0 1 1) faceted TiO₂(0 0 1) surface has been studied by Temperature Programmed Desorption (TPD) and X-ray Photoelectron Spectroscopy (XPS). Two main reactions were observed: dehydration to HCN and H₂O and decomposition to NH₃ and CO. The dehydration reaction was found to be three to four times larger than the decomposition at all coverages. Each of these reactions is found to occur in two temperature domains which are dependent upon surface coverage. The low temperature pathway (at about 400 K) is largely insensitive to surface coverage while the high temperature pathway (at about 500 K) shifts to lower temperatures with increasing surface coverage. These two temperature pathways may indicate two adsorption modes of formamide: molecular (via an η¹(O) mode of adsorption) and dissociative (via an η²(O,N) mode of adsorption). C1s and N1s XPS scans indicated the presence of multiple species after formamide absorption at 300 K. These occurred at ca. 288.5 eV (-CONH-) and 285 eV (sp³/sp² C) for the C1s and 400 eV-(NH₂), 398 eV (-NH) and 396 eV (N) for the N1s and result from further reaction of formamide with the surface.