Ab initio chemisorption studies of H on Fe(110)

Ab initio configuration interaction calculations are performed to study the chemisorption of atomic H on a Fe(110) surface. The lattice is modeled as an embedded three-layer, 40-atom cluster with the Fe atoms fixed at the bulk position. Fe 3d orbitals are explicitly included on five Fe atoms on the surface. Hydrogen strongly binds to the Fe(110) surface at the long-bridge, short-bridge, and quasi three-fold sites. The calculated adsorption energies are 2.76, 2.73, and 2.71 eV, respectively. H-surface bonding at the on-top Fe site is more than 0.4 eV weaker. The calculated H-surface distances are 0.89, 1.03, and 0.87 Å for H at the long-bridge, short-bridge, and quasi three-fold sites, respectively, which agrees well with the LEED value of 0.9 ± 0.1 Å. The H-surface stretching vibrational frequencies are calculated to be 1070, 1066, and 1073 cm⁻¹, at the long-bridge, short-bridge, and quasi three-fold sites, respectively. The work function of Fe(110) decreases on H adsorption. The present calculations indicate that H diffusion into the bulk through the short-bridge site will have a much higher activation barrier than via the long-bridge and quasi three-fold sites.     

Investigation of hydrogen interaction with the Si(1 1 1)-7 × 7 surface by STM with Bi/W tips

The STM technique with a special Bi/W tip was used to study the interaction of hydrogen atoms with the Si(1 1 1)-7 × 7 surface. The reactivity of different room temperature (RT) adsorption sites, such as adatoms (A), rest atoms (R), and corner holes (CH) was investigated. The reactivity of CH sites was found to be ∼2 times less than that of R and A sites. At temperatures higher than RT, hydrogen atoms rearrange among A, R, and CH sites, with increased occupation of R sites (T <  300 °C). Further temperature increase leads to hydrogen desorption, where its surface diffusion plays an active role. We discuss one of the possible desorption mechanisms, with the corner holes surrounded by a high potential barrier. Hydrogen atoms have a higher probability to overcome the desorption barrier rather than diffuse either into or out of the corner hole. The desorption temperature of hydrogen from CH, R, and A sites is about the same, equal to ∼500 °C. Also it is shown that hydrogen adsorption on the CH site causes slight electric charge redistribution over neighbouring adatoms, namely, increases the occupation of electronic states on A sites in the unfaulted halves of the Si(1 1 1)-7 × 7 unit cell. Based on these findings, the indirect method of investigation with conventional W tips was suggested for adsorbate interaction with CH sites.     

First-principles studies of the oxygen adsorption on unreconstructed and reconstructed Ni(1 1 0) surfaces

First-principles calculations have been performed to investigate the adsorption of oxygen on unreconstructed and reconstructed Ni(1 1 0) surfaces. The energetics, structural, electronic and magnetic properties are given in detail. For oxygen adsorption on unreconstructed surface, (n×1)(n=2,3) substrate with oxygen atom on short-bridge site is found to be the most stable adsorption configuration. Whereas energetically most favorable adsorption phase of reconstructed surface is p(n×1) substrate with oxygen atom located at long-bridge site. Our calculations suggest that the surface reconstruction is induced by the oxygen adsorption. We also find there are redistributions of electronic structure and electron transfer from the substrate to adsorbate. Our calculations also indicate surface magnetic moment is enhanced on clean surfaces and oxygen atoms are magnetized weakly after oxygen adsorption. Interestingly, adsorption on unreconstructed surface does not change surface magnetic moment. However, adsorbate leads to reduction of surface magnetic moment in reconstructed system remarkably.     

Sticking of atomic hydrogen on the tungsten (0 0 1) surface

The sticking of hydrogen atoms with kinetic energies in the range 0.003-10 eV on a clean (0 0 1) tungsten surface has been investigated using molecular dynamics simulations. The atoms are found to stick to the surface at 0 and 300 K, with a sticking coefficient smaller than 0.6 for kinetic energies higher than 3 meV. The adsorption sites for H on the W(0 0 1) surface are also presented. The dominant site is in perfect agreement with the experimentally found bridge site.     

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.     

Hydrogen desorption kinetics and band bending for 6H-SiC(0 0 0 1) surfaces

The desorption kinetics of hydrogen from polished 6H-SiC(0 0 0 1) surfaces exposed to various sources of hydrogen have been determined using temperature programmed desorption (TPD). For (3 × 3) 6H-SiC(0 0 0 1) surfaces prepared via annealing and cooling in SiH₄, desorption of 0.2 ± 0.05 monolayer of molecular hydrogen was observed to occur at ≈590 °C. This β₁ H₂ desorption peak exhibited second order kinetics with an activation energy of 2.4 ± 0.2 eV. For (3 × 3) 6H-SiC surfaces exposed to atomic hydrogen generated via either a hot rhenium filament or remote hydrogen plasma, low energy electron diffraction patterns showed an eventual conversion back to (1 × 1) symmetry. Spectra acquired using Auger electron and X-ray photoelectron spectroscopies revealed that the atomic hydrogen exposure removed the excess Si. Photoelectron spectroscopy results also showed a 0.5 eV increase in binding energy for the Si2p and C1s core levels after removal of the Si-Si bilayer that is indicative of a decrease in band bending at the SiC surface. TPD from the (3 × 3) 6H-SiC(0 0 0 1) surfaces exposed to atomic hydrogen showed substantially more molecular hydrogen desorption (1-2 ML) through the appearance of a new desorption peak (β_2,3) that started at ≈200 °C. The β_2,3 peak exhibited second order desorption kinetics and a much lower activation energy of 0.6 ± 0.2 eV. A third smaller hydrogen desorption state was also detected in the 650-850 °C range. This last feature could be resolved into two separate desorption peaks (α₁ and α₂) both of which exhibited second order kinetics with activation energies of 4.15 ± 0.15 and 4.3 ± 0.15 eV, respectively. Based on comparisons to hydrogen desorption from Si and diamond surfaces, the β and α desorption peaks were assigned to hydrogen desorption from Si and C sites, respectively.     

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.     

Study on the adsorption of fluorescein on Ag(1 1 0) substrate

The adsorption of fluorescein on the Ag(1 1 0) surface has been investigated by the first-principles pseudopotential method. Various adsorption geometries have been calculated and the energetically most favorable structure of fluorescein/Ag(1 1 0) was identified. The fluorescein molecule, in most favorable structure, is on hollow site, and the adsorption energy is 2.34 eV. Here the adsorption sites refer to the positions at the first layer of the substrate where the middle carbon atom of the fluorescein molecule is located. The bonding strength of the fluorescein molecule to the Ag substrate is site selective, being determined by electron transfer to the oxygen atoms of the molecule and local electrostatic attraction between the oxygen atoms and the silver atoms.     

Density functional theory study of selenium adsorption on Fe(1 1 0)

The adsorption of atomic Se on a Fe(1 1 0) surface is examined using the density functional theory (DFT). Selenium is adsorbed in high-symmetry adsorption sites: the -short and long-bridge, and atop sites at 1/2, 1/4, and 1 monolayer (ML) coverages. The long bridge (LB) site is found to be the most stable, followed by the short bridge (SB) and top sites (T). The following overlayer structures were examined, p(2 × 2), c(2 × 2), and p(1 × 1), which correspond to 1/4 ML, 1/2 ML, and 1 ML respectively. Adsorption energy is −5.23 eV at 1/4 ML. Se adsorption results in surface reconstruction, being more extensive for adsorption in the long bridge site at 1/2 ML, with vertical displacements between +8.63 and −6.69% -with regard to the original Fe position-, affecting the 1st and 2nd neighbours. The largest displacement in x or y-directions was determined to be 0.011, 0.030, and 0.021 Å for atop and bridge sites. Comparisons between Se-adsorbed and pure Fe surfaces revealed reductions in the magnetic moments of surface-layer Fe atoms in the vicinity of the Se. At the long bridge site, the presence of Se causes a decrease in the surface Fe d-orbital density of states between 4 and 5 eV below Fermi level. The density of states present a contribution of Se states at −3.1 eV and −12.9 eV. stabilized after adsorption. The Fe-Fe overlap population decrease and a Fe-Se bond are formed at the expense of the metallic bond.     

Hydrogen adsorption and diffusion on Pd(1 1 1)

The adsorption, diffusion and ordering of hydrogen on Pd(1 1 1) was studied by scanning tunneling microscopy in the temperature range of 37-90 K. At low coverage isolated hydrogen atoms were observed. They formed √3×√3-1H islands as the coverage increased. Above 1/3 monolayer (ML) coverage areas of a new phase with √3×√3-2H structure were formed, with both structures coexisting between 1/3 and 2/3 ML. Finally a 1 × 1 structure was formed after high exposures of hydrogen above 50 K, with a coverage close to 1 ML. Atomically resolved images reveal that H binds to fcc hollow sites.     

A density functional study of atomic oxygen and water molecule adsorption on Ni(1 1 1) and chromium-substituted Ni(1 1 1) surfaces

Density functional theory (DFT) for generalized gradient approximation calculations has been used to study the adsorption of atomic oxygen and water molecules on Ni(1 1 1) and different kind of Ni-Cr(1 1 1) surfaces. The fcc hollow site is energetically the most favorable for atomic oxygen adsorption and on top site is favorable for water adsorption. The Ni-Cr surface has the highest absorption energy for oxygen at 6.86 eV, followed by the hcp site, whereas the absorption energy is 5.56 eV for the Ni surface. The Ni-O bond distance is 1.85 Å for the Ni surface. On the other hand, the result concerning the Ni-Cr surface implies that the bond distances are 1.93-1.95 Å and 1.75 Å for Ni-O and Cr-O, respectively. The surface adsorption energy for water on top site for two Cr atom substituted Ni-Cr surface is 0.85 eV. Oxygen atoms prefer to bond with Cr rather than Ni atoms. Atomic charge analysis demonstrates that charge transfer increases due to the addition of Cr. Moreover, a local density of states (LDOS) study examines the hybridization occurring between the metal d orbital and the oxygen p orbital; the bonding is mainly ionic, and water bonds weakly in both cases.     

Coverage effects on phenol adsorption on Si(1 0 0)2 × 1 as: A first principle calculation

We investigated the adsorption of a 6-dimers Si(1 0 0)2 × 1 surface as a function of coverage and adsorption type (molecular/dissociative) by first principle calculations. In particular, we performed calculations on models with 2, 3, 4 and 6 phenol molecules, corresponding to coverage Θ = 0.34, 0.5, 0.67 and 1. We found that total adsorption energy, when at least one phenol is in a molecular state is lower than the sum of the corresponding singly adsorbed molecules. The dissociative adsorption of multiple molecules, both in parallel and switched configuration is most favoured for a coverage Θ = 0.34 (2.6 eV per adsorbed molecule). This values decreases to 2.0 eV and remains constant till the coverage 1 is reached.The energy barrier for the molecular-to-dissociated transition of a phenol molecule, in presence of another dissociatively adsorbed molecule is ∼0.008 eV and it is similar to the value in case of single adsorption. Possible hydrogen displacements were also considered.     

The interaction of carbon with Si(1 1 1):As and Si(1 1 1):H surfaces, a theoretical study

Density functional theory calculations have been performed to determine the adsorption site of carbon at the Si(1 1 1):As and Si(1 1 1):H surfaces at different coverages. The As- and H-passivated surfaces were simulated by replacing the topmost Si layer by As or by saturating the Si dangling bonds with hydrogen atoms, respectively. Different high symmetry sites were considered. Carbon was placed successively in the fourfold (T₄) or threefold coordinated (H₃), the ontop (T₁) sites or substituted for a Si atom in the S₅ position located underneath the Si adatom in the T₄ site. We found that the preferred carbon adsorption site depends on the coverage of the passivated surfaces. At low coverages i.e. at 1/16 ML and 1/3 ML, it prefers a distorted T₄ position whereas at 1 ML, it occupies an H₃ site. This contrasts with the clean surface where the most energetically favored site is the S₅ at all coverages. Carbon adsorption induces a significant change in the structural geometry of the surface atoms, leading to a charge re-arrangement in the surface layers.     

First-principles calculations of C diffusion through the surface and subsurface of Ag/Ni(1 0 0) and reconstructed Ag/Ni(1 0 0)

Density functional theory calculations are performed to investigate the C diffusion through the surface and subsurface of Ag/Ni(1 0 0) and reconstructed Ag/Ni(1 0 0). The calculated geometric parameters indicate the center of doped Ag is located above the Ni(1 0 0) surface owing to the size mismatch. The C binding on the alloy surface is substantially weakened, arising from the less attractive interaction between C and Ag atoms, while in the subsurface, the C adsorption is promoted as the Ag coverage is increased. The effect of substitutional Ag on the adsorption property of Ni(1 0 0) is rather short-range, which agrees well with the analysis of the projected density of states. Seven pathways are constructed to explore the C diffusion behavior on the bimetallic surface. Along the most kinetically favorable pathway, a C atom hops between two fourfold hollow sites via an adjacent octahedral site in the subsurface of reconstructed Ag/Ni(1 0 0). The “clock” reconstruction which tends to improve the surface mobility, is more favorable on the alloy surface because the c(2 × 2) symmetry is inherently broken by the Ag impurity. As a consequence, the local lattice strain induced by the C transport is effectively relieved by the Ag-enhanced surface mobility and the C diffusion barrier is lowered from 1.16 to 0.76 eV.     

Aluminium adsorption on Ir(1 1 1) at a quarter monolayer coverage: A first-principles study

We present an ab initio density-functional study for aluminium adsorption on Ir(1 1 1) at high symmetry sites, namely, the fcc-, hcp-hollow, top and bridge sites. In each case, we calculate the atomic geometry, average binding energy, work function, and surface dipole moment at the coverage of 0.25 monolayer. We find the favourable structure to be Al at threefold hcp-hollow site, with a corresponding binding energy of 4.46 eV. We present and compare the electronic properties of the two lowest energy structures, i.e., at the threefold hollow sites and discuss the nature of the Al-Ir bond and binding site preference. In particular, we observe a large hybridization of Al-3s, 3p and Ir-5d states near Fermi level, forming an inter-metallic bonds. This results in a significant electron transfer from the Al atoms to the Ir(1 1 1) substrate, inducing an outward pointing surface dipole moment and a large decrease in the work function of 1.69 eV for Al in the hcp-hollow site. Compared to the fcc-hollow site, adsorption in the hcp-hollow site results in a lower density-of-states at the Fermi level, as well as a greater hybridization in the bonding states.     

Dynamic study of W atoms and clusters on W (1 1 1) surfaces

Using a field ion microscope, the diffusion behaviors and atomic processes of W atoms and clusters on W (1 1 1) surfaces were observed directly. The activation energy of W clusters diffusion on W (1 1 1) as a function of cluster size has an oscillatory and increasing behavior. But, the activation energy of a single W atom is especially high. The compact geometric structures are more stable and have higher activation energies of surface diffusion than structures with extra atoms at the periphery. Besides the terrace diffusion, other atomic processes such as the ascending, descending, detachment motion on W (1 1 1) surfaces were also observed. Unlike the general systems, their occurrence temperatures are quite near. These experimental results were used to discuss the formation mechanism of single atom W tips.     

NO adsorption and diffusion on unreconstructed Pt{1 0 0} surface. A density functional theory investigation

Ab initio density functional theory was used to investigate the adsorption and diffusion of a single NO molecule on the unreconstructed Pt{1 0 0}-(1 × 1) surface. To our knowledge this is the first theoretical study of the NO diffusion activation energy on the Pt{1 0 0} surface. The most stable adsorption position for NO corresponds to the bridge site with the axis of the molecule perpendicular to the surface. The bond of the NO molecule to the surface is through the N-atom. We found that there is a low adsorption energy when the NO molecule is bonded through the O-atom and the axis is perpendicular to the surface, for the three high symmetry sites investigated. NO diffusion between bridge-hollow sites, bridge-atop sites, and hollow-atop sites was also investigated. The barrier for NO diffusion is 0.41 eV, which corresponds to the energy difference between the bridge and hollow sites. This value is around 15% of the highest adsorption energy found on this surface. NO stretch frequencies are also calculated for the three high symmetry sites investigated.     

Diffusion of single adatom Cu on Cu (0 0 1) and (1 1 0) surfaces

With static relaxation, the surface diffusion activation energies of a single Cu adatom migrated by both atomic exchange and hopping mechanisms and the forces acted on the diffusing adatom from other atoms of Cu (0 0 1) or (1 1 0) surface are calculated by using the MAEAM. When adatom migrated on Cu (0 0 1) or (1 1 0) surface, the increment curves of the system energy by hopping mechanism are symmetrical and the saddle points are in the midpoints of the migration path, but the ones by the exchange mechanism are dissymmetrical and the saddle points are always close to the initial hole positions of the adatom and away from the initial equilibrium positions of the exchanged atom. From minimization of both the diffusion activation energy and the force acted on the diffusing adatom from other atoms, we found that, on Cu (0 0 1) surface the favorable diffusion mechanism is hopping mechanism, however, on Cu (1 1 0) surface, hopping via long bridge is easier than the exchange mechanism but the hopping via short bridge is more difficult than the exchange mechanism.     

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.     

Effect of hydrogenation on B/Si(0 0 1)-(1 × 2)

Ab initio calculations, based on pseudopotentials and density functional theory, have been performed to investigate the effect of hydrogenation on the atomic geometries and the energetics of substitutional boron on the generic Si(0 0 1)-(1 × 2) surface. For a single B atom substitution corresponding to 0.5 ML coverage, we have considered two different sites: (i) the mixed Si-B dimer structure and (ii) boron substituting for the second-layer Si to form Si-B back-bond structure, which is energetically more favorable than the mixed Si-B dimer by 0.1 eV/dimer. However, when both of these cases are passivated by hydrogen atoms, the situation is reversed and the Si-B back-bond case becomes 0.1 eV/dimer higher in energy than the mixed Si-B dimer case. For the B incorporation corresponding to 1 ML coverage, among the substitutional cases, 100% interdiffusion into the third layer of Si and 50% interdiffusion into the second layer of Si are energetically similar and more favorable than the other cases that are considered. However, when the surface is passivated with hydrogen, the B atoms energetically prefer to stay at the third layer of the Si substrate.