Investigation of the dissociative adsorption for cyclopropane on the copper surface by density functional theory and quantum chemical molecular dynamics method

The dissociative adsorption of cyclopropane on the copper surface was studied using quantum chemical molecular dynamics method with “Colors-Excite” code and density functional theory by Amsterdam Density Functional program (ADF2000). The excited state of cyclopropane was used as adsorbate to simulate the dissociated adsorption under an irradiation energy of ca. 10 eV. One of the C-C bonds in cyclopropane was broken and the two new bonds between cyclopropane and copper surface were formed. The electrons transferred from the copper atoms to cyclopropane with a value of about 0.2e. The shorter distances between the carbons and surface copper atoms showed the existence of strong interaction. Consistently, the results indicated metallacyclopentane was the most possible intermediate species in dissociative adsorption by ADF2000 and “Colors-Excite” method.     

CF interaction with Si(1 0 0)-(2 × 1): Molecular dynamics simulation

In this study, the interaction of CF with the clean Si(1 0 0)-(2 × 1) surface at normal incidence and room temperature was investigated using molecular dynamics simulation. Incident energies of 2, 12 and 50 eV were simulated. C atoms, arising from dissociation, preferentially react with Si to form Si-C bonds. A Si_xC_yF_z interfacial layer is formed, but no etching is observed. The interfacial layer thickness increases with increasing incident energy, mainly through enhanced penetration of the silicon lattice. Silicon carbide and fluorosilyl species are formed at 50 eV, which is in good agreement with available experimental data. The level of agreement between the simulated and experimental results is discussed.     

Adsorption and dissociation of molecular hydrogen on Pt/CeO2 catalyst in the hydrogen spillover process: A quantum chemical molecular dynamics study

Ultra accelerated quantum chemical molecular dynamics method (UA-QCMD) was used to study the dynamics of the hydrogen spillover process on Pt/CeO₂ catalyst surface for the first time. The direct observation of dissociative adsorption of hydrogen on Pt/CeO₂ catalyst surface as well as the diffusion of dissociative hydrogen from the Pt/CeO₂ catalyst surface was simulated. The diffusion of the hydrogen atom in the gas phase explains the high reactivity observed in the hydrogen spillover process. Chemical changes, change of adsorption states and structural changes were investigated. It was observed that parallel adsorption of hydrogen facilitates the dissociative adsorption leading to hydrogen desorption. Impact with perpendicular adsorption of hydrogen causes the molecular adsorption on the surface, which decelerates the hydrogen spillover. The present study also indicates that the CeO₂ support has strong interaction with Pt catalyst, which may cause an increase in Pt activity as well as enhancement of the metal catalyst dispersions and hence increasing the rate of hydrogen spillover reaction.     

DFT study on dissociative adsorption of SiH4 and GeH4 on SiGe(1 0 0)-2 × 1 surface

SiH₄ and GeH₄ dissociative adsorptions on a buckled SiGe(1 0 0)-2 × 1 surface have been analyzed using density functional theory (DFT) at the B3LYP level. The Ge alloying in the Si(1 0 0)-2 × 1 surface affects the dimer buckling and its surface reactivity. Systematic Ge influences on the reaction energetics are found in SiH₄ and GeH₄ reactions with four dimers of Si^∗-Si, Ge^∗-Si, Ge^∗-Ge, and Si^∗-Ge (∗ denotes the protruded atom). On a half H-covered surface, the energy barriers for silane and germane adsorption are higher than those on the pristine surface. The energy barrier for silane adsorption is higher than the corresponding barrier for germane adsorption. Rate constants are also calculated using the transition-state theory. We conclude that the SiGe surface reactivity in adsorption reaction depends on the Ge presence in dimer form. If the surface Ge is present in form of Ge^∗-Ge, the surface reactivity decreases as the Ge^∗-Ge content increases. If the surface Ge prefers to be in form of Ge^∗-Si at low Ge contents, the surface reactivity increases first, then decreases at high Ge surface contents when Ge^∗-Ge prevails. The calculated rate constant ratio of GeH₄ adsorption on Si^∗-Si over Ge^∗-Ge at 650 °C is 2.1, which agrees with the experimental ratio of GeH₄ adsorption probability on Si(1 0 0) over Si(1 0 0) covered by one monolayer Ge. The experimental ratio is 1.7 measured through supersonic molecular beam techniques. This consistency between calculation and experimental results supports that one monolayer of Ge on Si(1 0 0) exists in form of Ge^∗-Ge dimer.     

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.     

Density functional study of oxygen adsorption on the Mg3Nd (0 0 1) surface

The adsorption of oxygen atoms on Mg₃Nd (0 0 1) surface was studied based on density function theory (DFT), in which the exchange-correlation potential was chosen as the generalized gradient approximation (GGA) in the Perdew and Wang (PW91). The most preferred adsorption position was at the top-hollow site. Upon the optimization on top-hollow site with different coverage, it was found that the adsorption energy decreased with oxygen coverage. The density of states analysis showed that obvious charge transfer took place between O atom and the nearest Nd atom and chemical bond formed between O atom and the nearest Nd atom after O adsorption. The result of surface energy as a function of chemical potential change of oxygen indicated the clean Mg₃Nd (0 0 1) surface was easy to adsorb oxygen and form 1.00 ML surface.     

Gd adsorption on the Mo(2 1 1) surface

The adsorption of Gd thin layers on the Mo(2 1 1) face was investigated by using Auger electron spectroscopy (AES), low electron energy diffraction (LEED), scanning tunneling microscopy (STM), X-ray photoemission spectroscopy (XPS) and measurements of the work function changes (Δφ). It was found that at 300 K Gd does not form any dilute chain structures and from the very beginning of the adsorption process Gd forms a densely packed layer. The dilute p(4 × 1) chain structure was observed by LEED after annealing thin layers (θ < 1 ML) to temperatures above 770 K. STM images confirm the existence of the p(4 × 1) structure islands. The intermixing of the substrate and adsorbate atoms takes place.     

Low reactivity of molecular oxygen with Si(1 1 1)-c(2 × 8)

The adsorption of molecular oxygen on the c(2 × 8) reconstruction of quenched Si(1 1 1) surfaces has been studied at the atomic scale using scanning tunneling microscopy (STM) at room temperature (RT). It has been found that clean well reconstructed c(2 × 8) adatoms do not react with O₂ molecules but that a limited oxidation can start where adatom sites arranged in reconstructed structures are present. Comparison between O₂ adsorption on Si(1 1 1)-c(2 × 8) and Si(1 1 1)-7 × 7 reconstructions coexisting on the same quenched silicon surface has been carried out in detail. For each atomic site present on the surface the variation of reacted sites with exposure has been measured. For low O₂ exposures, bright and dark oxygen induced sites appear on the Si(1 1 1)-7 × 7, while Si(1 1 1)-c(2 × 8) does not oxidized at all. At high O₂ exposures, large oxidation areas have spread on the 7 × 7 reconstruction, preferentially on the faulted halves of the unit cell, and smaller oxidation areas induced by topological defects have grown all around clean un-reacted c(2 × 8) regions.     

Structural stabilities and diffusion of small Fe clusters on Fe (1 1 0) surface: A molecular dynamics study

Using Embedded-atom-method (EAM) potential of iron, structural stabilities of small Fe clusters on a Fe (1 1 0) surface have been investigated by molecular dynamics studies. It is presented that a tetramer and heptamer clusters are more stable than other sizes. These two clusters have high transition energies. They can be a critical nucleus at low and high temperature, respectively. A dimer diffuses more easily with lower energy barrier than single adatom. The trimer's rotation and dimer shearing mechanisms have been investigated in this paper.     

The surface stress of the (1 1 0) and (1 0 0) surfaces of rutile and the effect of water adsorbents

The surface stress on clean TiO₂ (1 1 0) and (1 0 0) surfaces, and those with four types of adsorbent - (i) molecularly adsorbed water, (ii) dissociatively adsorbed water, (iii) dissociatively adsorbed water at an oxygen vacancy, and (iv) adsorbed hydrogen - was investigated in the framework of density functional theory using a slab model. The calculations were intended to rationalize the effect of the artificially introduced stress that occurs in experimentally photoinduced hydrophilicity. Tensile stress was observed for a clean (1 1 0) surface, and a mixture of tensile and compressive stress for a clean (1 0 0) surface. The adsorbate-induced surface stresses were analyzed in terms of the sixfold coordinated character of the surface titanium atoms, hydrogen bonds between the adsorbents and the bridging oxygen atoms, and the change in electron density in the vicinity of the surface.     

Water adsorption and dissociation on BeO(0 0 1) and (1 0 0) surfaces

Plateaus in water adsorption isotherms on hydroxylated BeO surfaces suggest significant differences between the hydroxylated (1 0 0) and (0 0 1) surface structures and reactivities. Density functional theory structures and energies clarify these differences. Using relaxed surface energies, a Wulff construction yields a prism crystal shape exposing long (1 0 0) sides and much smaller (0 0 1) faces. This is consistent with the BeO prisms observed when beryllium metal is oxidized. A water oxygen atom binds to a single surface beryllium ion in the preferred adsorption geometry on either surface. The water oxygen/beryllium bonding is stronger on the surface with greater beryllium atom exposure, namely the less-stable (0 0 1) surface. Water/beryllium coordination facilitates water dissociation. On the (0 0 1) surface, the dissociation products are a hydroxide bridging two beryllium ions and a metal-coordinated hydride with some surface charge depletion. On the (1 0 0) surface, water dissociates into a hydroxide ligating a Be atom and a proton coordinated to a surface oxygen but the lowest energy water state on the (1 0 0) surface is the undissociated metal-coordinated water. The (1 0 0) fully hydroxylated surface structure has a hydrogen bonding network which facilitates rapid proton shuffling within the network. The corresponding (0 0 1) hydroxylated surface is fairly open and lacks internal hydrogen bonding. This supports previous experimental interpretations of the step in water adsorption isotherms. Further, when the (1 0 0) surface is heated to 1000 K, hydroxides and protons associate and water desorbs. The more open (0 0 1) hydroxylated surface is stable at 1000 K. This is consistent with the experimental disappearance of the isotherm step when heating to 973 K.     

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.     

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.     

Theoretical study of hydrogen adsorption on the GaN(0 0 0 1) surface

Ab initio density functional theory, using the B3LYP hybrid functional with all-electron basis sets, has been applied to the adsorption of H on the (0 0 0 1) surface of wurtzite GaN. For bulk GaN, good agreement is obtained with photoemission and X-ray emission data for the valence band and for the Ga 3d and N 2s shallow core levels. A band gap of E_g = 4.14 eV is computed vs the experimental value (at 0 K) of 3.50 eV. A simple model, consisting of a (2 × 2) structure with 3/4-monolayer (ML) of adsorbed H, is found to yield a density of states in poor agreement with photoemission data for H adsorbed on surfaces prepared by ion bombardment and annealing. A new model, consisting of co-adsorbed Ga (1/4 ML) and H (1/2 ML), is proposed to account for these data.     

Photoemission study of the (2 × 2) structure formed by H2O adsorption on the Zr(0 0 0 1) surface

We have studied adsorption of water on the Zr(0 0 0 1) surface at sub-monolayer coverage by means of LEED and photoemission spectroscopy. An ordered (2 × 2) structure is formed after adsorption of 0.6-1.4 Langmuirs at 473 K. The sharpest LEED pattern was observed at an exposure of 1.2 L implying a coverage of 0.5 ML of oxygen. The same exposure at 293 K gives only a weak and diffuse (2 × 2) pattern. In addition, the sharp (2 × 2) pattern obtained at 473 K can be reversibly weakened by cooling to 293 K and subsequently sharpened by heating. For the sharp (2 × 2) structure, valence band spectra indicated dissociation of water and showed a peak composed mainly of O 2p derived states with two components at 6.0 eV and 6.6 eV binding energy. On cooling to 293 K, the O 2p peak became narrower and a new state appeared at 7.9 eV. Two components of the O 1s core level were resolved for the (2 × 2) structure, assigned to oxide and hydroxyl groups. The hydrogen on the surface of Zr(0 0 0 1) resulting from the dissociation of water and from bulk segregation strongly influenced the formation of the (2 × 2) structure of oxygen, and caused a temperature instability of the structure.     

CH4 dissociation on NiCo (1 1 1) surface: A first-principles study

A density-functional theory method has been conducted to investigate the adsorption of CH_x (x = 0-4) as well as the dissociation of CH_x (x = 1-4) on (1 1 1) facets of ordered NiCo alloy. The results have been compared with those obtained on pure Ni (1 1 1) surface. It shows that the adsorption energies of C and CH are decreased while it is increased for CH₃ on NiCo (1 1 1) compared to those on pure Ni (1 1 1). Furthermore, on NiCo (1 1 1), dissociation of CH_x prefers not to the top of Ni, but to the top of Co. The rate-determining step for CH₄ dissociation is considered as the first step of dehydrogenation on NiCo (1 1 1), while it is the fourth step of dehydrogenation on Ni (1 1 1). Furthermore, the activation barrier in rate-determining step is slightly higher by 0.07 eV on Ni (1 1 1) than that on NiCo (1 1 1). From above results, it is important to point out that carbon is easy to form on NiCo (1 1 1) although the adsorption energy of C atom is slightly decreased compared to that on Ni (1 1 1).     

Theoretical study of oxygen adsorption at the Fe(1 1 0) and (1 0 0) surfaces

Electronic, magnetic and structural properties of atomic oxygen adsorbed in on-surface and subsurface sites at the two most densely packed iron surfaces are investigated using density functional theory combined with a thermodynamics formalism. Oxygen coverages varying from a quarter to two monolayers (MLs) are considered. At a 1/4 ML coverage, the most stable on-surface adsorption sites are the twofold long bridge sites on the (1 1 0), and the fourfold-hollow sites on the (1 0 0) surface. The presence of on-surface oxygen atoms enhances the magnetic moments of the atoms of the two topmost Fe layers. Detailed results on the surface magnetic properties, due to O incorporation, are presented as well. Subsurface adsorption is found unfavored. The most stable subsurface O, in tetrahedral positions at the (1 0 0) and octahedral ones at the (1 1 0) surface, are characterized by substantially lower binding than that in the on-surface sites. Subsurface oxygen increases the interplanar distance between the uppermost Fe layers. The preadsorbed oxygen overlayer enhances binding of subsurface O atoms, particularly for tetrahedral sites beneath the (1 1 0) surface.     

Molecular dynamics study on surface structure and surface energy of rutile TiO2 (1 1 0)

The formula for surface energy was modified in accordance with the slab model of molecular dynamics (MDs) simulations, and MD simulations were performed to investigate the relaxed structure and surface energy of perfect and pit rutile TiO₂(1 1 0). Simulation results indicate that the slab with a surface more than four layers away from the fixed layer expresses well the surface characteristics of rutile TiO₂ (1 1 0) surface; and the surface energy of perfect rutile TiO₂ (1 1 0) surface converges to 1.801±0.001 J m⁻². The study on perfect and pit slab models proves the effectiveness of the modified formula for surface energy. Moreover, the surface energy of pit surface is higher than that of perfect surface and exhibits an upper-concave parabolic increase and a step-like increase with increasing the number of units deleted along [0 0 1] and [1 1 0], respectively. Therefore, in order to obtain a higher surface energy, the direction along which atoms are cut out should be chosen in accordance with the pit sizes: [] direction for a small pit size and [0 0 1] direction for a big pit size; or alternatively the odd units of atoms along [1 1 0] direction are removed.     

A theoretical study of hydrogen adsorption and diffusion on a W(1 1 0) surface

We have used density functional theory to investigate hydrogen adsorption and diffusion on a W(1 1 0) surface. Hydrogen adsorption structures were examined from low coverage to one monolayer, and a threefold hollow site was found to be the most stable site at all coverages. In contrast to previous assertions, the work function decrease is not due to electron transfer from the hydrogen atoms to the W surface, but due to electron depletion at the vacuum region above the hydrogen atoms. Hydrogen atoms can diffuse via short-bridge sites and long-bridge sites at a coverage of θ = 1.0. Although the calculated activation energy for hydrogen diffusion via a short-bridge site is as small as 0.05 eV, field emission microscope experiments have shown that the activation energy for hydrogen diffusion is about 0.20 eV, which agrees fairly well with our calculated value of the activation energy via a long-bridge site. This discrepancy can be related to hydrogen delocalization on the W(1 1 0) surface, which has been suggested by electron energy loss spectroscopy experiments.     

Substitutional adsorption geometry for Pb(1 1 1)-(√3 × √3)R30°-K

Intermixed structures for alkalis (larger than Li) on close-packed substrates have previously been observed only on Al(1 1 1). This study shows that K forms an ordered intermixed structure on Pb(1 1 1). The structures of clean Pb(1 1 1) and Pb(1 1 1)-(√3 × √3)R30°-K were studied using dynamical low-energy electron diffraction (LEED). The clean Pb(1 1 1) surface at 47 K was found to be a relaxed version of the bulk structure, in agreement with an earlier study of the same surface [Y.S. Li, F. Jona, P.M. Marcus, Phys. Rev. B 43 (1991) 6337]. At room temperature, adsorption of K on this surface results in a (√3 × √3)R30° structure, which was shown using dynamical LEED to consist of K atoms substituted in surface vacancies. The K-Pb bond length was found to be 3.62 ± 0.3 Å, with no significant change to the Pb interlayer spacings.