Density functional theory investigation of CN on Cu(1 1 1), Ni(1 1 1) and Ni(1 0 0)

The adsorption of CN on Cu(1 1 1), Ni(1 1 1) and Ni(1 0 0) has been investigated using density functional theory (DFT). While experimental studies of CN on Cu(1 1 1) show the molecular axis to be essentially parallel to the surface, the normally-preferred DFT approach using the generalised gradient approximation (GGA) yields a lowest energy configuration with the C-N axis perpendicular to the surface, although calculations using the local density approximation (LDA) do indicate that the experimental geometry is energetically favoured. The same conclusions are found for CN on Ni(1 1 1); on both surfaces bonding through the N atom is always unfavourable, in contrast to some earlier published results of ab initio calculations for Ni(1 1 1)/CN and Ni(1 0 0)/CN. The different predictions of the GGA and LDA approaches may lie in subtly different relative energies of the CN 5σ and 1π orbitals, a situation somewhat similar to that for CO adsorbed on Pt(1 1 1) which has proved challenging for DFT calculations. On Ni(1 0 0) GGA calculations favour a lying-down species in a hollow site in a geometry rather similar to that found experimentally and in GGA calculations for CN on Ni(1 1 0).     

Adsorption of thiophene on Ni(1 0 0), Cu(1 0 0), and Pd(1 0 0) surfaces: ab initio periodic density functional study

Adsorption of thiophene on the (1 0 0) surfaces of Ni, Cu, and Pd has been investigated by the ab initio density functional theory method (periodic DMol³). Several parallel and perpendicular adsorption geometries are examined in detail. For Ni(1 0 0), both dissociative and molecular adsorption structures are found with small difference in energy. Thiophene adsorbs only molecularly on Cu(1 0 0) and Pd(1 0 0). The most stable molecular adsorption structures on all the surfaces are quite similar, where thiophene adsorbs on top of a 4-fold hollow with the symmetry axis rotated 45° from the metal rows. These stable structures arise from a good matching of the thiophene molecule to the metal surfaces. The calculated adsorption geometries are in reasonable agreement with XAFS experiments.     

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.     

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.     

O2 dissociation on Pd(2 1 1) and Cu(2 1 1) surfaces

We have performed ab initio Density Functional Theory (DFT) based calculations to observe the reactivity of the Pd(2 1 1) and Cu(2 1 1) surfaces towards O₂. In order to properly address the adsorption dynamics, the static potential energy surface calculations have been complemented with first principles molecular dynamics calculations, which reveal interesting steering effects that complicate the dissociation dynamics. We have found that on both surfaces the step microfacets are very reactive and the dissociation of the O₂ molecule at room temperature occurs mostly on those sites.     

Density functional theory investigation of the structure of SO2 and SO3 on Cu(1 1 1) and Ni(1 1 1)

Density functional theory (DFT) slab calculations, mainly using the generalised gradient approximation, have been used to investigate the minimum energy structures of molecular SO₂ and SO₃ on Cu(1 1 1) and Ni(1 1 1) surfaces. On Ni(1 1 1) the optimal local adsorption structures are in close agreement with experimental results for both molecular species obtained using the X-ray standing wavefield technique, although for adsorbed SO₂ the energetic difference between two alternative lateral positions of the lying-down molecule on the surface is marginally significant. On Cu(1 1 1) the results for adsorbed SO₂, in particular, were sensitive to the DFT functional used in the calculations, but in all cases failed to reproduce the experimentally-established preference for adsorption with the molecular plane perpendicular to the surface. This result is discussed in the context of previously published DFT results for these species adsorbed on Cu(1 0 0). The optimal geometry found for SO₃ on Cu(1 1 1) is similar to that on Ni(1 1 1), providing agreement with experiment regarding the molecular orientation but not the adsorption site.     

Structure of enantiopure and racemic alanine adlayers on Cu(1 1 0)

Plane wave density functional theory has been employed to analyze the structure of alanine adlayers on the Cu(1 1 0) surface. Alanine forms (3 × 2) adlayers on Cu(1 1 0) that are closely related to the structures of glycine on the same surface. There is essentially no energy difference between the most stable racemic and enantiopure alanine adlayers. This observation implies that adsorption of racemic alanine on Cu(1 1 0) will result in a pseudoracemate adlayer.     

The adsorption site of methoxy and ethoxy on Cu(1 0 0) as determined by vibrational spectroscopy and first principle calculations

In the present work we determine the adsorption site of two polyatomic molecules, methoxy and ethoxy, on Cu(1 0 0). This is accomplished by comparing experimental intramolecular vibrational modes to the corresponding modes calculated by first principle methods. We explore the three different high symmetry adsorption sites on top, bridge and hollow using several different metal clusters to represent the Cu(1 0 0) surface. The experimental results for both methoxy and ethoxy are best reproduced by the most realistic representation of the hollow position, a Cu₅ cluster.     

Structure and dynamics of the missing-row reconstruction on O/Cu(0 0 1) and O/Ag(0 0 1)

The geometry and the vibrational properties of missing row reconstructed O/Cu(0 0 1) and O/Ag(0 0 1) surfaces are investigated by means of density functional theory and density functional perturbation theory, using the local density and the generalized-gradient approximations. Our results predict very similar structural and vibrational properties for the two reconstructed surfaces. In the case of copper our calculations reproduce quite accurately the experimental results, while for the missing row reconstructed O/Ag(0 0 1) surface the agreement between theory and experiment is less satisfactory.     

Quantum states of hydrogen (H, D, T) atoms on Cu(1 0 0) and (1 1 0) surfaces

We investigate the quantum mechanical behavior of adsorbed hydrogen (H, D, T) on Cu(1 0 0) and (1 1 0) surfaces. We construct potential energy surfaces (PESs) for the motion of the hydrogen H atom on Cu(1 0 0) and (1 1 0) surfaces within the framework of density functional theory. The potential energy takes a minimum value on the hollow site of Cu(1 0 0) and on the short bridge site of Cu(1 1 0). Moreover, we calculate the quantum states of hydrogen atom motion on these calculated PESs. The ground state wave function of the hydrogen atom motion is strongly localized around the hollow site on the Cu(1 0 0) surface. On the other hand, the ground state wave function of the hydrogen atom motion on Cu(1 1 0) is distributed from the short bridge site to two neighboring pseudo-threefold sites. We finally show isotope effects on the quantum states of the motion of hydrogen on both surfaces.     

Density functional theory studies on stress stabilization of the Cu(1 1 0) striped phase

Under oxygen exposure, the Cu(1 1 0) surface shows a striped phase consisting of alternating bare (1 1 0) surface areas and added-row (2 × 1)O reconstructions. Density functional theory is used to show that the major origin for the formation of the striped phase is the elastic interaction between these areas. The difference between the surface stress of the bare surface and the (2 × 1)O covered surface is predicted to be 1.3 N/m in reasonable agreement with values derived from grazing incidence X-ray diffraction. Supercell calculations for periods of up to 62 Å confirm that the formation of the striped phase is favorable compared to an added-row (4 × 1)O reconstruction with the same coverage. But the predicted equilibrium period of roughly 30 Å is significantly smaller than in experiment. The calculations are impeded by the surface energy alternating with the number of layers in the slab. This behavior is related to a quantum well behavior of the Cu 4s-electrons. A simple model for this behavior is discussed and compared to ab initio results.     

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.     

CO methanation on Ni(1 1 1) and modified Ni3Al(1 1 1) surfaces: A first-principle study

The adsorption energies of intermediates in CO methanation on the modified Ni₃Al(1 1 1) surface and the Ni(1 1 1) surface are calculated using density functional theory. A microkinetic analysis based on the calculated adsorption energies is performed to explain the different kinetics of CO methanation catalyzed by Ni₃Al and Ni powders. The electronic structures of different atoms on the modified Ni₃Al surface are also presented, and correlate well with the adsorption energies and geometries.     

Adsorption of chlorophenol on the Cu(1 1 1) surface: A first-principles density functional theory study

The interaction between a 2-chlorophenol (C₆H₄OHCl) molecule and the Cu(1 1 1) surface has been investigated using density functional theory as an initial step in gaining a better understanding of the catalyzed formation of dioxin compounds on a clean copper surface. The 2-chlorophenol molecule is found to form several weakly bonded, horizontally and vertically oriented configurations. Dissociative modes have also been investigated. For the latter, the formation of phenyl and benzyne fragments is found to be more energetically favourable than the formation of 2-chlorophenoxy radicals.     

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.     

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.     

SiO adsorption on a p(2 × 2) reconstructed Si(1 0 0) surface

We have investigated the adsorption mechanism of SiO molecule incident on a clean Si(1 0 0) p(2 × 2) reconstructed surface using density functional theory based methods. Stable adsorption geometries of SiO on Si surface, as well as their corresponding activation and adsorption energies are identified. We found that the SiO molecule is adsorbed on the Si(1 0 0) surface with almost no activation energy. An adsorption configuration where the SiO binds on the channel separating the dimer rows, forming a Si-O-Si bridge on the surface, is the energetically most favourable geometry found. A substantial red-shift in the calculated vibrational frequencies of the adsorbed SiO molecule in the bridging configurations is observed. Comparison of adsorption energies shows that SiO adsorption on a Si(1 0 0) surface is energetically less favourable than the comparable O₂ adsorption. However, the role of SiO in the growth of silicon sub-oxides during reactive magnetron plasma deposition is expected to be significant due to the relatively large amount of SiO molecules incident on the deposition surface and its considerable sticking probability. The stable adsorption geometries found here exhibit structural properties similar to the Si/SiO₂ interface and may be used for studying SiO_x growth.     

Methanol adsorption on silicon (0 0 1)

Using a first-principles pseudopotential technique, we have investigated the adsorption of CH₃OH on the Si(0 0 1) surface. We have found that, in agreement with the overall experimental picture, the most probable chemisorption path for methanol adsorption on silicon (0 0 1) is as follows: the gas phase CH₃OH adsorbs molecularly to the electrophilic surface Si atom via the oxygen atom and then dissociates into Si-OCH₃ and H, bonded to the electrophilic and nucleophilic surface silicon dimer atoms, respectively. Other possible adsorption models and dissociation paths are also discussed. Our calculations also suggest that the most probable methanol coverage is 0.5 ML, i.e., one molecule per Si-Si dimer, in agreement with experimental evidences. The surface atomic and electronic structures are discussed and compared to available theoretical and experimental data. In addition, we propose that a comparison of our theoretical STM images and calculated vibrational modes for the adsorbed systems with detailed experimental investigations could possibly confirm the presented adsorption picture.     

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 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.