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

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.     

Molecular precursor-mediated methanol dissociation on Si(1 1 1)7 × 7: ab initio study

We present an ab initio study of methanol interaction with the Si(1 1 1)7 × 7 surface using a Si(1 1 1)4 × 2 model. The study of the methanol dissociation on Si(1 1 1)4 × 2 shows that pair dissociation on adatom-restatom dangling bonds is largely favoured, in agreement with the experimental observations. The “center” type adatom is slightly more reactive than the “corner” type one, although the difference is weak. Similar behaviour is observed in both adatom types. Our results for a direct CH₃OH dissociation favouring a basic cleavage (adsorption of OH and CH₃ fragments) rather than an acidic one (adsorption of H and OCH₃ fragments), we are finally led to take a kinetic effect into consideration to reconcile theory with experiment. We show that the presence of molecular precursor states is possible. Different orientations with respect to the silicon dangling bonds of these molecular precursors are investigated. However, the corresponding energies are very close and, considering their relative energies, it is finally difficult to discriminate between acidic and basic cleavages.     

Comparative DFT study of the adsorption of 1,3-butadiene, 1-butene and 2-cis/trans-butenes on the Pt(1 1 1) and Pd(1 1 1) surfaces

The interaction of 1,3-butadiene, 1-butene and 2-cis/trans-butenes on the Pt(1 1 1) and Pd(1 1 1) surfaces has been studied with density functional theory methods (DFT). The same most stable adsorption modes have been found on both metal surfaces with similar adsorption energies. For 1,3-butadiene the 1,2,3,4-tetra-σ adsorption structure is shown to be the most stable one, in competition with a 1,4-metallacycle-type mode, which is only less stable by 10-12 kJ mol⁻¹. On Pt(1 1 1) these total energy calculations were combined with simulations of the vibrational spectra. This confirms that the 1,2,3,4-tetra-σ adsorption is the most probable adsorption structure, but cannot exclude the 1,4-metallacycle as a minority species. Although similar in type and energy, the adsorption on the Pd(1 1 1) surface shows a markedly different geometry, with a smaller molecular distortion upon adsorption. The most stable adsorption structure for the butene isomers is the di-σ-mode. Similarly to the case of the 1,3-butadiene, the adsorption geometry is closer to the gas phase one on Pd than on Pt, hence explaining the different spectroscopic results, without the previously assumed requirement of a different binding mode. Moreover the present study has shown that the different selectivity observed on Pt(1 1 1) and Pd(1 1 1) for the hydrogenation reaction of butadiene cannot be satisfactory explained by the single comparison of the relative stabilities of 1,3-butadiene and 1-butene on these metals.     

First principle calculations of benzotriazole adsorption onto clean Cu(1 1 1)

The adsorption of benzotriazole (BTAH or C₆N₃H₅) on a Cu(1 1 1) surface is investigated by using first principle density functional theory calculations (VASP). It is found that BTAH can be physisorbed (<0.1 eV) or weakly chemisorbed (∼0.43 eV) onto Cu(1 1 1), and the chemical bond is formed through nitrogen sp² lone pairs. The weak chemisorption can be stabilized by reaction with neighboring protonphilic radicals, like OH⁻. Furthermore, the geometries and associated energies of intermolecular hydrogen bonds between adsorbates on Cu(1 1 1) are also calculated. A model of the first layer of BTAH/BTA⁻ on Cu(1 1 1) surface is developed based on a hydrogen bond network structure.     

On-surface and sub-surface oxygen on ideal and reconstructed Cu(1 0 0)

In order to understand the first steps of the Cu(1 0 0) oxidation we performed first principles calculations for on-surface and sub-surface oxygen on this surface. According to our calculations, the adsorption energies for all on-surface site oxygen atoms increase, whereas the energies of the sub-surface atoms decrease with the increasing oxygen coverage. At coverage 1 ML and higher on the reconstructed surface, structures including both on- and sub-surface atoms are energetically more favourable than structures consisting only of on-surface adsorbates. On the ideal (1 0 0) surface this change can be perceived at coverage 0.75 ML.     

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.     

Atomic and molecular adsorption on Pt(1 1 1)

The adsorption of several atomic (H, O, N, S, and C) and molecular (N₂, HCN, CO, NO, and NH₃) species and molecular fragments (CN, CNH₂, NH_2, NH, CH₃, CH₂, CH, HNO, NOH, and OH) on the (1 1 1) facet of platinum, an important industrial and fuel cell catalyst, was studied using self-consistent periodic density functional theory (DFT-GGA) calculations at a coverage of 1/4 ML. The best binding site, energy, and position, as well as an estimated diffusion barrier, of each species were determined. The binding strength for all the species can be ordered as follows: N₂ < NH₃ < HCN < NO < CO < CH₃ < OH < NH₂ < H < CN < NH < O < HNO < CH₂ < NOH < CNH₂ < N < S < CH < C. Although the atomic species generally preferred fcc sites, there was no clear trend in site preference by the molecular species or molecular fragments. The vibrational frequencies of all the stable adsorbates in their best and second best adsorption sites were calculated and found to be in good agreement with experimental values reported in the literature. Finally, the decomposition thermochemistry of NOH, HNO, NO, NH₃, N₂, CO, and CH₃ was analyzed.     

CO and hydrogen adsorption on Pd(2 1 0)

We have studied the adsorption of CO on Pd(2 1 0) by performing density functional theory (DFT) calculations within the generalized gradient approximation. We find a relatively small corrugation in the CO adsorption energies with the two bridge sites being energetically almost degenerate. CO is furthermore known as a strong poison in heterogeneous catalysis. We have therefore also addressed the coadsorption of CO with atomic hydrogen. There is a significant inhibition of the hydrogen adsorption due to the presence of CO which is analysed in terms of the electronic structure of the adsorbate system.     

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.     

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.     

First-principles study on field evaporation of surface atoms from W(0 1 1) and Mo(0 1 1) surfaces

The simulations of field-evaporation processes for surface atoms on W(0 1 1) and Mo(0 1 1) surfaces are implemented using first-principles calculations based on the real-space finite-difference method. The threshold values of the external electric field for evaporation of the surface atoms, which are ∼6 V/Å for tungsten and ∼5 V/Å for molybdenum, are in agreement with the experimental results. While the threshold value of the electric field and the local-field enhancement around the evaporating atoms agree with those expected from the conclusion of the previous study using structureless jellium, the induced charge around the surface atom has a significant difference from that obtained by the jellium model.     

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

Adsorption of cyclohexadiene, cyclohexene and cyclohexane on Pt(1 1 1)

The adsorption of 1,3-cyclohexadiene, 1,4-cyclohexadiene, cyclohexene and cyclohexane on Pt(1 1 1) was studied using ab initio density functional theory. For 1,3-cyclohexadiene three adsorption modes were distinguished: bridge 1,2-di-σ/3,4-π, hollow 1,4-di-σ/2,3-π and bridge 1,4-di-σ/2,3-π with adsorption energies of −155, −147 and −75 kJ/mol, respectively. Three stable adsorption modes were also identified for 1,4-cyclohexadiene: bridge quadra-σ, hollow di-σ/π and bridge di-π with adsorption energies of −146 kJ/mol, −142 kJ/mol and −88 kJ/mol, respectively. Cyclohexene was found to adsorb in six modes: 4 di-σ and 2 π-adsorption modes. The preferred configuration was found to be boat di-σ with an adsorption energy of −81 kJ/mol. The three other di-σ adsorption modes have comparable adsorption energies, ranging from −64 to −69 kJ/mol. Molecular strain and CPt bonding energies are used to elucidate stability trends. Cyclohexane is found to adsorb only at the hollow site whereby the axial hydrogen atoms are positioned over surface Pt-atoms with an adsorption energy of −37 kJ/mol. The calculations correctly predict the weakening of the axial CH bonds and provide a possible explanation for the large shift in the vibrational frequencies.     

Cu(1 1 1) and Cu(0 0 1) surface electronic states. Comparison between theory and experiment

Recent advances in both the experimental resolution and in the computational capabilities motivate new studies of surface properties. In particular, a detailed comparison between theoretical and experimental data is expected to provide a better insight into surface and image states. In this work we present a joint effort analyzing such features of the Cu(1 1 1) and Cu(0 0 1) surfaces. The experiments are performed by using both linear and non-linear angle-resolved photoemission. From the theoretical point of view, we make use of the Green function embedding technique within density functional theory. We are able to account for the image states by suitably modifying the effective potential in the Kohn-Sham equation and the generalized boundary conditions on the Green function. Comprehensive theoretical and experimental results on the effective mass and the binding energy of the Shockley state and the first image states are reported.     

Coadsorption of methyl radicals and oxygen on Rh(1 1 1)

The chemisorption of CH₃ on Rh(1 1 1) is studied to understand the origin of the weakened symmetric stretch mode. A few different explanations for this weakened mode have been suggested in previous studies. These include C-H bond depletion and donation into C-H anti-bond orbitals either in an upright or tilted geometry. We investigate these possibilities by performing first-principles density functional calculations. Our results show strong adsorption at all high-symmetry sites with methyl in two possible orientations. A thorough analysis of the adsorption geometry shows that C_3v symmetry is preferred over a tilted species, ruling out tilting as a mechanism for C-H mode softening. Evidence of a multi-center bond between methyl and the surface rhodium atoms (similar to the kind shown recently by Michaelides and Hu for methyl on Ni(1 1 1)) is presented, showing that C-H bond depletion is the cause of mode-softening for methyl on Rh(1 1 1). Experimental results have shown that mode-softening diminishes when an electronegative species is coadsorbed, suggesting that donation into C-H anti-bonding orbitals is the mechanism for mode-softening. We therefore examine the coadsorption of oxygen and methyl on Rh(1 1 1). Our results suggest a new model for the effect of O on CH₃. Analysis of charge density differences shows that the dominant initial effects of O coadsorption are the removal of charge from the C-surface bond and the transfer of charge to the C-H bond. Subsequent increase of the H-Rh distance further reduces mode softening.     

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.     

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.     

Scanning tunneling microscopy study on initial stage of atomic hydrogen adsorption on the Si(1 1 1) 7 × 7 surface

Initial hydrogen adsorption on the Si(1 1 1) 7 × 7 surface was studied by scanning tunneling microscopy (STM) in an ultrahigh vacuum. Room temperature adsorbed hydrogen on the adatom in the 7 × 7 reconstruction led to depression of adatoms in the STM images. The hydrogen uptake curve at the adatom site as a function of hydrogen exposure time was well represented by Langmuir adsorption. No preferential adsorption was seen among four inequivalent adatoms in the 7 × 7 reconstruction. Adsorption of the adjacent center and corner adatoms respectively showed ∼10% higher adsorption. Even though the number of reacted adatoms in the half unit of the 7 × 7 reconstruction was statistically random, the number of reacted adatoms in the nearest neighbor half unit was enhanced as the number of reacted sites increased in the half unit.