DFT study of NH3 dissociation on Si(1 1 1)-7 × 7. The role of intermolecular interactions

The adsorption of NH₃ molecule on the Si(1 1 1)-7 × 7 surface modelled with a cluster has been studied using density functional theory (DFT). The results indicate the existence of a precursor state for the non-dissociative chemisorption. The active site for the molecular chemisorption is the adatom; while the NH₃ molecule adsorbs on the Si restatom via this preadsorbed state, the adsorption on the Si adatom is produced practically without an energy barrier. The ammonia adsorption on the adatom induces an electron transfer from the dangling bond of this atom to the dangling bond of the adjacent Si restatom, hindering this site for the adsorption of a second NH₃ incoming molecule. However, this second molecule links strongly by means of two H-bonds. The dissociative chemisorption process was studied considering one and two ammonia molecules. For the dissociation of a lonely NH₃ molecule an energy barrier of ∼0.3 eV was calculated, yielding NH₂ on the adatom and H on the restatom. When two molecules are adsorbed, the NH₃-NH₃ interaction yields the weakening of a N-H bond of the ammonia molecule adsorbed closer the Si surface. As a consequence, the dissociation barrier practically disappears. Thus, the presence of a second NH₃ molecule at the adatom-restatom pair of the Si(1 1 1)-7 × 7 surface makes the dissociative reaction self-assisted, the total adsorption process elapsing with a negligible activation barrier (less than 0.01 eV).     

Characterization of diamond (1 0 0) surface with oxygen termination

Ab initio density functional theory (DFT) was employed to study reconstructions of diamond (1 0 0) surfaces in the presence of hydrogen, oxygen and hydroxyl. Clean and (2 × 1):1H surfaces are taken as reference. The properties of oxidization diamond surfaces with several adsorption structures, namely, O-on-top (OT) site, O-bridge (BR) site, hydroxyl (-OH), hydroxyl/hydroxyl, OT/hydroxyl, BR/hydroxyl have been considered. The calculated results indicate that the BR model is much more stable than the OT model, and the most energetically favorable structures of oxygenated surfaces are those with chemisorbed hydroxyl (-OH) group. Furthermore, the stability of the structures is also discussed from the point of HOMO-LUMO gap. Analysis of electronic structures shows that the presence of hydrogen induces surface conductivity whereas oxygen weakens it.     

Olefin bonding mechanism and growth to (1 0 0)(2 × 1):H diamond

The bonding and growth mechanism of photochemically attached olefin molecules to (1 0 0)(2 × 1):H diamond is characterized using atomic force (AFM) and scanning tunneling microscopy (STM) experiments in combination with molecular orbital calculations. To identify growth schemas, diamond surfaces after 10, 40 and 90 min of photo-chemically stimulated growth have been characterized. These data show clearly island formation which is discussed taking into account a growth model from silicon. The island growth shows no directional properties which are attributed to arrangement and geometrical properties of hydrogen terminated carbon bonds at the surface of (1 0 0) oriented (2 × 1) reconstructed diamond.     

Study of vacancy on diamond (1 0 0) (2 × 1) surface from first-principles

Structure and energy related properties of neutral and charged vacancies on relaxed diamond (1 0 0) (2 × 1) surface were investigated by means of density functional theory. Calculations indicate that the diffusion of a single vacancy from the top surface layer to the second layer is not energetically favored. Analysis of energies in charged system shows that neutral state is most stable on diamond (1 0 0) (2 × 1) surface. The multiplicity of possible states can exist on diamond (1 0 0) surface in dependence on the surface Fermi level, which supports that surface diffusion of a vacancy is mediated by the change of vacancy charge states. Analysis of density of states shows surface vacancy can be effectively measured by photoelectricity technology.     

Atomic and electronic structure of methanol on Ge(1 0 0)

We have performed density-functional theory (DFT) calculations to investigate the adsorption structures of methanol on a Ge(1 0 0) surface. Among many possible adsorption configurations, the most favorable configurations at room temperature were found to be those in which the OH-dissociated methanol molecule forms O-Ge bonds, with the methoxy group either parallel or perpendicular to the Ge surface. The spatial arrangement of methoxy group relative to the Ge(1 0 0) surface is not critical. The dissociated H is bonded to an adjacent up-Ge atom, passivating the dangling bond. The possibility of H diffusion to other Ge atoms is also investigated. The corresponding simulated images explain well the adsorption features observed experimentally. The reaction pathways explain the feasibility of OH-dissociative structures at room temperature. The two OH-dissociative configurations where methoxy groups are either parallel or perpendicular to Ge surfaces are similar in thermodynamic and kinetic aspects.     

Characteristic energy band values and electron attenuation length of a chemical-vapor-deposition diamond (0 0 1)2 × 1 surface

The characteristic energy band values such as the Fermi-level position with respect to valence band top for a boron-doped p-type hydrogen-terminated chemical-vapor-deposition (CVD) diamond (0 0 1)2 × 1 surface and for a clean CVD diamond (0 0 1)2 × 1 surface have been determined by a new method with an accuracy of ±0.02 eV. The electron attenuation length for the clean diamond (0 0 1)2 × 1 surface for the electron kinetic energy of C 1s X-ray photoemission peak by Mg Kα excitation is experimentally determined to be 2.1-2.2 nm. These values are compared and discussed with the previously reported experimental and simulation values.     

The dissociative adsorption of borane on the Ge(1 0 0)-2 × 1 surface: A density functional theory study

By means of cluster models coupled with density functional theory, we have studied the hydroboration of the Ge(1 0 0)-2 × 1 surface with BH₃. It was found that the Ge(1 0 0) surface exhibits rather different surface reactivity toward the dissociative adsorption of BH₃ compared to the C(1 0 0) and Si(1 0 0) surfaces. The strong interaction still exists between the as-formed BH₂ and H adspeices although the dissociative adsorption of BH₃ on the Ge(1 0 0) surface occurs readily, which is in distinct contrast to that on the C(1 0 0) and Si(1 0 0) surfaces. This can be understood by the electrophilic nature of the down Ge atom, which makes it unfavourable to form a GeH bond with the dissociating proton-like hydrogen. Alternatively, it can be attributed to the weak proton affinity of the Ge(1 0 0) surface. Nevertheless, the overall dissociative adsorption of BH₃ on group IV semiconductor surfaces is favourable both thermodynamically and kinetically, suggesting the interesting analogy and similar diversity chemistry of solid surface in the same group.     

Ammonia on Ni(1 1 1) surface studied by first principles: Bonding, multilayers structure and comparison with experimental IR and XPS data

DFT calculations have been performed on the adsorption of NH₃ on Ni(1 1 1) to obtain information on the structure of the absorbed species, the nature of the chemical interactions between the adsorbate and the surface and the structure of multilayers formed at high coverage. A cluster model, using localized basis functions as well as an approach based on plane waves and periodic boundary conditions have been considered. The two approaches lead to similar results for the relative stabilities of investigated adsorption sites (atop > fcc > hcp) with an adsorption energy of about 15-24 kcal/mol (depending on the coverage). On the atop site, α-ammonia adsorbs molecularly with an equilibrium distance between the nitrogen atom and the surface of 2.03 ± 0.02 Å and a geometry close to the one of the molecule in the gas phase. The good agreement between the two DFT approaches clearly underlines the local nature of the adsorption reaction. The vibrational frequencies computed for NH₃ adsorbed in this site are in good agreement with experimental values. They show that the interaction with the surface leads to a weakening of the strength of the N-H bond while the angular stretching is stronger. Both orbital and topological analyses were used to investigate chemical interactions between the cluster and the molecule. The results strongly suggest an electrostatic (non-covalent) interaction between the substrate and the molecule. Calculations with NH₃ coverages above 0.25 confirm that saturation occurs at a coverage of 0.25. Above the saturation coverage, ammonia molecules in excess form multiple layers (β and γ ammonia) bonded to the first layer by intermolecular hydrogen bonds. N 1s core level shifts calculations performed for the several investigated coverages are also in good agreement with experimental XPS data. It is shown that the H-bond network more than the bond to the surface allows to understand the N 1s core level variations.     

Layer-resolved photoelectron diffraction from Si(0 0 1) and GaAs(0 0 1)

Photoelectron diffraction in the layer-resolved mode brings more detailed information about local atomic arrangement than is obtained in the standard mode. This is demonstrated in crystals with diamond and zinc-blende structures, both for unpolarized photon excitation as well as for circularly polarized excitation. The full angular distributions of photoemission intensities are evaluated for large atomic clusters representing ideally truncated surfaces of Si(0 0 1) and GaAs(0 0 1). Highly structured layer-resolved patterns enable a more detailed understanding of the standard mode outcomes. Photoelectron intensities from atomic layers placed at different depths under the crystal surface provide direct evidence about electron attenuation and its anisotropy in crystals.     

Vibrationally-assisted dissociative adsorption of oxygen on Ru(0 0 0 1)-p(2 × 1)-O

Supersonic molecular beam technique combined with high resolution X-ray photoelectron spectroscopy using synchrotron radiation was applied to the study of the dynamics of dissociative adsorption of oxygen on Ru(0 0 0 1) surface in high coverage region. The Ru(0 0 0 1) surface pre-covered with oxygen atoms of 0.5 monolayer, which corresponds to the p(2 × 1)-O structure, was dosed to oxygen molecules with translational energy of 0.5 eV. Oxygen uptake was compared between the cases with and without the beam source heated in order to verify the effects of internal energy of oxygen. We found drastic enhancement in initial sticking probability of oxygen when the beam source was heated to 1400 K. We concluded that the enhancement of sticking probability is mainly caused by molecular vibrational excitation, indicating that dissociation barrier is located in the exit channel on potential energy surface.     

CO2 adsorption on Cr(1 1 0) and Cr2O3(0 0 0 1)/Cr(1 1 0)

We attempt to correlate qualitatively the surface structure with the chemical activity for a metal surface, Cr(1 1 0), and one of its surface oxides, Cr₂O₃(0 0 0 1)/Cr(1 1 0). The kinetics and dynamics of CO₂ adsorption have been studied by low energy electron diffraction (LEED), Aug er electron spectroscopy (AES), and thermal desorption spectroscopy (TDS), as well as adsorption probability measurements conducted for impact energies of E_i = 0.1-1.1 eV and adsorption temperatures of T_s = 92-135 K. The Cr(1 1 0) surface is characterized by a square shaped LEED pattern, contamination free Cr AES, and a single dominant TDS peak (binding energy E_d = 33.3 kJ/mol, first order pre-exponential 1 × 10¹³ s⁻¹). The oxide exhibits a hexagonal shaped LEED pattern, Cr AES with an additional O-line, and two TDS peaks (E_d = 39.5 and 30.5 kJ/mol). The initial adsorption probability, S₀, is independent of T_s for both systems and decreases exponentially from 0.69 to 0.22 for Cr(1 1 0) with increasing E_i, with S₀ smaller by ∼0.15 for the surface oxide. The coverage dependence of the adsorption probability, S(Θ), at low E_i is approx. independent of coverage (Kisliuk-shape) and increases initially at large E_i with coverage (adsorbate-assisted adsorption). CO₂ physisorbs on both systems and the adsorption is non-activated and precursor mediated. Monte Carlo simulations (MCS) have been used to parameterize the beam scattering data. The coverage dependence of E_d has been obtained by means of a Redhead analysis of the TDS curves.     

Interaction of 2,2,6,6-tetramethyl-3,5-heptanedione with the Si(1 0 0)-2 × 1 surface: Scanning tunneling microscopy and density functional theory study

Room temperature adsorption and reaction of 2,2,6,6-tetramethyl-3,5-heptanedione (dpmH) on the Si(1 0 0)-2 × 1 surface has been studied with ultra-high vacuum scanning tunneling microscopy (UHV-STM) and temperature programmed desorption (TPD). The molecule is found to chemisorb as a mixture of at least five distinct species. Density functional theory (DFT) was used to calculate the structures and adsorption energies of 12 possible addition products. Unique bonding assignments for each experimental feature are proposed by consideration of a common intermediate reaction network, and a comparison of possible reaction pathways leading to the final products. These assignments are: OH inter-dimer dissociation, OH intra-dimer dissociation, 1,5 intra-dimer addition, 1,5 inter-dimer addition, and intra-dimer [2 + 2]CO addition with OH dissociation on an adjacent dimer. TPD and STM results show that the molecule dissociates completely upon annealing to 700 °C with formation of the c(4 × 4) phase at low exposures, and SiC islands for exposures exceeding 0.15 L.     

Adsorption of carbon monoxide on Pd(3 1 1) and (2 1 1) surfaces

Adsorption of carbon monoxide on Pd(3 1 1) and (2 1 1) stepped surfaces has been investigated by the extended London-Eyring-Polyani-Sato (LEPS) method constructed using a 5-parameter Morse potential. The calculated results show that there exist common characteristics of CO adsorption on the two surfaces. At low coverage, CO occupies threefold hollow site of the (1 1 1) terrace and is tilted with respect to the surface normal. Among the threefold hollow sites on the (1 1 1) terrace, the nearer the site is to the step, the greater is the influence of the step. The twofold bridge site on the (1 0 0) step is also a stable adsorption site at high coverage. Because of the different lengths of the (1 1 1) terraces, the (3 1 1) and (2 1 1) stepped surfaces have different characteristics. A number of new sites are exposed on the boundary regions, including the fourfold hollow site (H₄) of the (3 1 1) surface and the fivefold hollow site (H₅) of the (2 1 1) surface. At high coverage, CO resides in the H₅ site of the (2 1 1) surface, but the H₄ site of the (3 1 1) surface is not a stable adsorption site. This study further shows that the on-top site on the (1 0 0) step of Pd(3 1 1) is a stable adsorption site, but the same type of site on Pd(2 1 1) is not.     

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.     

Coadsorption of CN and O on Cu (1 0 0) surface: A density functional study

The adsorption of cyanide (CN) or oxygen atom, as well as the coadsorption of CN + O on Cu (1 0 0) surface is studied by using density functional theory (DFT) and the cluster model method. Cu₁₄ cluster is used to simulate the surface. Perpendicular and parallel bonding geometries of CN adsorbed on Cu (1 0 0) surface are considered, respectively. The present calculations show that the CN may be absorbed on top and bridge sites by carbon atom of cyanide (C-down), and C-down on top site is the most favorable. The adsorbed C-N stretch frequencies compared with that of the gaseous CN species are all red-shifted, except the C-down on top site. The charge transfer from the surface to the CN species leads to an increase in work function for the Cu surface. The oxygen atom adsorbed on the four-fold hollow site of Cu (1 0 0) is the most favorable, and is consistent with the experimental study. The coadsorption of O at a four-fold hollow site tends to block adsorption of CN at the nearby sites. If O coverage increases, the CN may be adsorbed on the top and bridges sites with the C-down model. The reaction CN + O → OCN on the Cu (1 0 0) is predicted to be exothermic, and formed OCN species may be stably absorbed on the Cu (1 0 0).     

First-principle study of NO adsorption on the LaFeO3 (0 1 0) surface

The adsorption of NO molecule on the LaFeO₃ (0 1 0) surface was studied using first-principle calculations based on density functional theory. The calculated results indicate that the Fe-top site is the most favorable for NO adsorption. The N-O bond length, Mulliken charge, and the N-O vibration frequency of the NO molecule are discussed after adsorption. The analysis results of the density of the states show that when NO is adsorbed with the Fe-NO configuration, the bonding mechanism is mainly from the interaction between the NO and the Fe d orbit.     

Optical second-harmonic generation study of incorporation of nitrogen atoms at Si(1 0 0) surfaces

The kinetics of N incorporation into Si layers has been studied on the Si(1 0 0) surface probing surface dangling bonds with an optical second-harmonic (SH) generation during surface exposure to N atoms generated by a radio frequency N₂ plasma. It is observed that SH intensity decreases with N dose. The rate of decrease in SH intensity apparently decreases with surface temperature, whereas total amount of N atoms taken on the surface remains constant. This fact suggests that N atoms are incorporated at the subsurface layers at higher temperatures. It is shown that the N incorporation at the subsurface layers proceeds by the indiffusion of N atoms either directly or indirectly via intermediate, metastable adsorption at the first surface layer. Applying Kisliuk adsorption model, activation energies of the N indiffusion are evaluated to be 0.30 ± 0.03 and 0.34 ± 0.05 eV for the indirect and direct path, respectively.     

A density functional theory study of ethylene adsorption on Ni10(1 1 1), Ni13(1 0 0) and Ni10(1 1 0) surface cluster models and Ni13 nanocluster

Ethylene adsorption was studied by use of DFT/B3LYP with basis set 6-31G(d,p) in Gaussian’03 software. It was found that ethylene has adsorbed molecularly on all clusters with π adsorption mode. Relative energy values were calculated to be −50.86 kcal/mol, −20.48 kcal/mol, −32.44 kcal/mol and −39.27 kcal/mol for Ni₁₃ nanocluster, Ni₁₀(1 1 1), Ni₁₃(1 0 0) and Ni₁₀(1 1 0) surface cluster models, respectively. Ethylene adsorption energy is inversely proportional to Ni coordination number when Ni₁₀(1 1 1), Ni₁₃(1 0 0) and Ni₁₀(1 1 0) cluster models and Ni₁₃ nanocluster are compared with each other.     

Carbon monoxide adsorption on the metastable fcc-Co/Cu(1 0 0) surface

The adsorption properties of CO on the epitaxial five-monolayer Co/Cu(1 0 0) system, where the Co overlayer has stabilized in the metastable fcc-phase, are reported. This system is known to exhibit metallic quantum well (MQW) states at energies 1 eV or greater above the Fermi level, which may influence CO adsorption. The CO/fcc-Co/Cu(1 0 0) system was explored with low energy electron diffraction (LEED), inverse photoemission (IPE), reflection-absorption infrared spectroscopy (RAIRS) and temperature programmed desorption (TPD). Upon CO adsorption, a new feature is observed in IPE at 4.4 eV above E_F and is interpreted as the CO 2π^∗ level. When adsorbed at room temperature, TPD exhibits a CO desorption peak at ∼355 K, while low temperature adsorption reveals additional binding configurations with TPD features at ∼220 K and ∼265 K. These TPD peak temperatures are correlated with different C-O stretch vibrational frequencies observed in the IR spectra. The adsorption properties of this surface are compared to those of the surfaces of single crystal hcp-Co, as well as other metastable thin film systems.     

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.