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.     

CO2 adsorption on the bimetallic Zn-on-Cu(1 1 0) system

Adsorption probability measurements (molecular beam scattering) have been conducted to examine the adsorption dynamics (i.e. the gas-surface energy transfer processes) of CO₂ adsorption on the Zn-on-Cu(1 1 0) bimetallic system. The results indicate surface alloy formation, which is in agreement with prior studies. Depositing Zn at 300 K on Cu(1 1 0), above the condensation temperature of CO₂, leads to a “blocking” of CO₂ adsorption sites by Zn which is incorporated in the Cu(1 1 0) surface. This apparent site blocking effect indicates a lowering of the CO₂ binding energy on the alloyed surface as compared with the clean Cu(1 1 0) support. The Zn coverage has been calibrated by Auger electron spectroscopy and thermal desorption spectroscopy.     

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

Molecularly intact and dissociative adsorption of water on clean Cu(1 1 0): A comparison with the water/Ru(0 0 1) system

An X-ray photoelectron spectroscopy (XPS) study was undertaken of the water/Cu(1 1 0)-system finding non-dissociative adsorption on clean Cu(1 1 0) at temperatures below 150 K. Thermally induced dissociation of D₂O is observed to occur above 150 K, similar to the H₂O/Ru(0 0 1) system, with an experimentally derived activation barrier of 0.53-0.56 eV which is very close in magnitude to the derived activation barrier for desorption of 0.50-0.53 eV. X-ray and electron induced damage to the water overlayer was quantified and used to rationalize the results of a recent XPS study of the water/Cu(1 1 0)-system where partial dissociation was observed already at 90 K.     

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.     

Comparative study of H2 adsorption on W(1 0 0)-c(2 × 2)Cu and W(1 0 0): Surface alloying effects

The interactions of H and H₂ with W(1 0 0)-c(2 × 2)Cu and W(1 0 0) have been investigated through density functional theory (DFT) calculations to elucidate the effect of Cu atoms on the reactivity of the alloy. Cu atoms do not alter the attraction towards top-W sites felt by H₂ molecules approaching the W(1 0 0) surface but make dissociation more difficult due to the rise of late activation barriers. This is mainly due to the strong decrease in the stability of the atomic adsorbed state on bridge sites, the most favourable ones for H adsorption on W(1 0 0). Still, our results show unambiguously that H₂ dissociative adsorption on perfect terraces of the W(1 0 0)-c(2 × 2)Cu surface is a non-activated process which is consistent with the high sticking probability found in molecular beam experiments at low energies.     

Cobaltocene adsorption and dissociation on Cu(1 1 1)

Photoemission results indicate that the initial adsorption of cobaltocene on Cu(1 1 1) at 150 K leads to molecular fragmentation, but with subsequent cobaltocene exposures, molecular absorption occurs. The molecularly adsorbed species is either adsorbed with only a fraction of molecules adopting a preferential orientation along the surface normal or adsorbed with the molecular axis away from the surface normal. This adsorption behavior is compared to nickelocene and ferrocene adsorption.     

The local structure of SO2 and SO3 on Ni(1 1 1)

The normal incidence X-ray standing wave (NIXSW) technique, supported by X-ray photoelectron spectroscopy and near-edge X-ray absorption fine structure (NEXAFS), has been used to determine the local adsorption geometry of SO₂ and SO₃ on Ni(1 1 1). Chemical-state specific NIXSW data for coadsorbed SO₃ and S, formed by the disproportionation of adsorbed SO₂ after heating from 140 K to 270 K, were obtained using S 1s photoemission detection. For adsorbed SO₂ at 140 K the new results confirm those of an earlier study [Jackson et al., Surf. Sci. 389 (1997) 223] that the molecule is located above hollow sites with its molecular plane parallel to the surface and the S and O atoms in off-atop sites; corrections to account for the non-dipole effects in the interpretation of the NIXSW monitored by S 1s and O 1s photoemission, not included in the earlier work, remove the need for any significant adsorption-induced distortion of the SO₂ in this structure. SO₃, not previously investigated, is found to occupy an off-bridge site with the C_3v axis slightly tilted relative to the surface normal and with one O atom in an off-atop site and the other two O atoms roughly between bridge and hollow sites. The O atoms are approximately 0.87 Å closer to the surface than the S atom. This general bonding orientation for SO₃ is similar to that found on Cu(1 1 1) and Cu(1 0 0) both experimentally and theoretically, although the detailed adsorption sites differ.     

Benzene adsorption and oxidation on Ir(1 1 1)

Adsorption, decomposition and oxidation of benzene on Ir(1 1 1) was studied by high resolution (synchrotron) XPS, temperature programmed desorption and low energy electron diffraction. Molecular adsorption of benzene on Ir(1 1 1) is observed between 170 K and 350 K. Above this temperature both desorption and decomposition of benzene take place. An ordered adsorbate structure was observed upon adsorption around 335 K. Decomposition involves C-C bond breaking as the formation of CH_ad is observed. The presence of a saturated O_ad layer (0.5 ML) weakens molecular benzene adsorption and suppresses decomposition.     

Photoemission study of glycine adsorption on Cu/Au(1 1 1) interfaces

The adsorption of glycine on Au(1 1 1) pre-deposited with different amounts of Cu was investigated with both conventional X-ray photoelectron spectroscopy (XPS) and synchrotron-based photoemission. In the Cu submonolayer range, glycine physically adsorbs on the Cu/Au(1 1 1) surfaces in its zwitterionic form and completely desorbs at 350 K. The C 1s, O 1s and N 1s core level binding energies monotonically increase with Cu coverage. This indicates that, in the Cu submonolayer range, the admetal is alloyed with Au rather than forming overlayers on the Au(1 1 1) substrate, consistent with our recent experimental and theoretical results [X. Zhao, P. Liu, J. Hrbek, J.A. Rodriguez, M. Pérez, Surf. Sci. 592 (2005) 25]. Upon increasing the amount of deposited Cu over 1 ML, part of the glycine overlayer transforms from the zwitterionic form to the anionic form (NH₂CH₂COO⁻) and adsorbs chemically on the Cu/Au(1 1 1) surface with the N 1s binding energy shifted by −2.3 eV. When the amount of deposited Cu is at 3.0 or 6.0 ML, the intensity of the N 1s chemisorption peak increases with aging time at 300 K. It indicates that glycine adsorption induces Cu segregation from the subsurface region onto the top layer of the substrate. Judging from the initial N 1s peak intensities, it is concluded that 64% and 36% of the top layer are still occupied by Au atoms before glycine adsorption even when the amounts of deposited Cu are 3.0 and 6.0 ML, respectively. On the Au(1 1 1) surface pre-dosed with 6.0 ML of Cu, part of the chemisorbed glycine will desorb and part will decompose upon heating to 450-500 K. In addition, about 20% of the glycine exists in the neutral form when the glycine overlayer was dosed on Cu/Au(1 1 1) held at 100 K.     

Theoretical study of CO and Pb adsorption on the Ni(1 1 1) and Ni3Al(1 1 1) surfaces

Adsorption of CO molecules and Pb atoms on the Ni(1 1 1) and Ni₃Al(1 1 1) substrates is studied theoretically within an ab initio density-functional-theory approach. Stable adsorption sites and the corresponding adsorption energies are first determined for stoichiometric surfaces. The three-fold hollow sites (fcc for Pb and hcp for CO) are found most favourable on both substrates. Next, the effect of surface alloying by a substitution of selected topmost substrate atoms by Pb or Ni atoms on the adsorption characteristics is investigated. When the surface Al atoms of the Ni₃Al(1 1 1) substrate are replaced by Ni atoms, the Pb and CO adsorption energies approach those for a pure Ni(1 1 1) substrate. The Pb alloying has a more substantial effect. On the Ni₃Al(1 1 1) substrate, it reduces considerably adsorption energy of CO. On the Ni(1 1 1) substrate, CO binding strengthens slightly upon the formation of the Ni(1 1 1)p(2×2)-Pb surface alloy, whereas it weakens drastically when the Ni(1 1 1)-Pb surface alloy is formed.     

Determination of a (4 × 4) structure formed on a Cu(0 0 1) surface by adsorption of calcium

The adsorption of calcium (Ca) atoms on a Cu(0 0 1) surface has been studied by low-energy electron diffraction (LEED) at 130, 300 and 400 K. It is found that a (4 × 4) was the only LEED pattern appeared at 400 K while a quasi-hexagonal structure was formed in a wide range of submonolayer coverage at 130 K. At 300 K, the (4 × 4) LEED spots were broad and weak. The (4 × 4) structure formed at 400 K was determined by a tensor LEED I-V analysis. It is a new-type of surface alloys consisting of five substitutional Ca atoms, nine surface Cu atoms, and two atomic vacancies in the unit cell. In spite of a quite large size-difference between Ca (3.94 Å) and Cu (2.55 Å) atoms, all Ca atoms are located at the substitutional sites. Among surface alloys so far reported, the atomic size ratio between Cu and Ca in the (4 × 4), 1.54, is the largest. Optimized structural parameters reveal that large lateral displacements of surface Cu atoms, being enabled by the appearance of the vacancies, allow the formation of the (4 × 4) structure.     

Adsorption and reactivity of SO2 on Ir(1 1 1) and Rh(1 1 1)

The adsorption and reactivity of SO₂ on the Ir(1 1 1) and Rh(1 1 1) surfaces were studied by surface science techniques. X-ray photoelectron spectroscopy measurements showed that SO₂ was molecularly adsorbed on both the Ir(1 1 1) surface and the Rh(1 1 1) surface at 200 K. Adsorbed SO₂ on the Ir(1 1 1) surface disproportionated to atomic sulfur and SO₃ at 300 K, whereas adsorbed SO₂ on the Rh(1 1 1) surface dissociated to atomic sulfur and oxygen above 250 K. Only atomic sulfur was present on both surfaces above 500 K, but the formation process and structure of the adsorbed atomic sulfur on Ir(1 1 1) were different from those on Rh(1 1 1). On Ir(1 1 1), atomic sulfur reacted with surface oxygen and was completely removed from the surface, whereas on Rh(1 1 1), sulfur did not react with oxygen.     

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.     

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.     

Surface and subsurface oxide formation on Ni(1 0 0) and Ni(1 1 1)

The oxidation of Ni(1 0 0) and Ni(1 1 1) at elevated temperatures and large oxygen exposures, typical of the methods used in the preparation of NiO(1 0 0) films for surface studies, has been investigated by medium energy ion scattering (MEIS) using 100 keV H⁺ incident ions. Oxide film growth proceeds significantly faster on Ni(1 1 1) than on Ni(1 0 0), but on both surfaces oxide penetration occurs to depths significantly greater than 100 Å with total exposures of 1200 and 6000 L respectively. The metal/oxide interface is extremely rough, with metallic Ni extending to the surface, even for much thicker oxide films on Ni(1 1 1). On Ni(1 1 1), NiO growth occurs with the (1 0 0) face parallel to the Ni(1 1 1) surface and the close-packed 〈1 1 0〉 directions parallel. On Ni(1 0 0) the MEIS blocking curves cannot be reconciled with a single orientation of NiO(1 0 0) (with the 〈1 1 0〉 directions parallel) on the surface, but is consistent with the substantial orientational disorder (including tilt) previously identified by spot-profile analysis LEED.     

Comparison of O adsorption on Ni3Al (0 0 1), (0 1 1), and (1 1 1) surfaces through first-principles calculations

First-principles calculations were performed to study the properties of O adsorption on Ni₃Al (0 0 1), (0 1 1), and (1 1 1) surfaces using the Cambridge serial total package (CASTEP) code. Stable adsorption sites are identified. The atomic and electronic structures and adsorption energies are predicted. The adsorption sites for O on the Ni₃Al (0 0 1) surface are at the 2Ni–2Al fourfold hollow site, whereas O prefers to adsorb at the Ni–Al bridge site on (0 1 1) surface and 2Ni–Al threefold hollow site on (1 1 1) surface. It is found that O shows the strongest affinity for Al and the state of O is the most stabilized when O adsorbs on (0 0 1) surface, while the affinity of O for Al on (0 1 1) surface is weaker than (0 0 1) surface, and (1 1 1) surface is the weakest. The stronger O and Al affinity indicates more stable Al₂O₃ when oxidized. The experiment has shown that the oxidation resistance of single crystal superalloy in different orientations improves in the order of (1 1 1), (0 1 1), and (0 0 1) surface, suggesting that the oxidation in different crystallographic orientations may be related to the affinity of O for Al in the surface.     

A kinetic Monte Carlo/UBI-QEP model of O2 adsorption on fcc (1 1 1) metal surfaces

The previously developed kinetic Monte Carlo model of molecular oxygen adsorption on fcc (1 0 0) metal surfaces has been extended to fcc (1 1 1) surfaces. The model treats uniformly all elementary steps of the process—O₂ adsorption, dissociation, recombination, desorption, and atomic oxygen hopping—at various coverages and temperatures. The model employs the unity bond index—quadratic exponential potential (UBI-QEP) formalism to calculate coverage-dependent energetics (atomic and molecular binding energies and activation barriers of elementary steps) and a Metropolis-type algorithm including the Arrhenius-type reaction rates to calculate coverage- and temperature-dependent features, particularly the adsorbate distribution over the surface. Optimal values of non-energetic model parameters (the spatial constraint, a travel distance of “hot” atoms, attempt frequencies of elementary steps) have been chosen. Proper modifications of the fcc (1 0 0) model have been made to reflect structural differences in the fcc (1 1 1) surface, in particular the presence of two different hollow sites (fcc and hcp). Detailed simulations were performed for molecular oxygen adsorption on Ni(1 1 1). We found that at very low coverages, only O₂ adsorption and dissociation were effective, while O₂ desorption and O₂ and O diffusion practically did not occur. At a certain O + O₂ coverage, the O₂ dissociation becomes the fastest process with a rate one-two orders of magnitude higher than adsorption. Dissociation continuously slows down due to an increase in the activation energy of dissociation and due to the exhaustion of free sites. The binding energies of both molecular and atomic oxygen decrease with coverage, and this leads to greater mobility of atomic oxygen and more pronounced desorption of molecular oxygen. Saturation is observed when the number of adsorbed molecules becomes approximately equal to the number of desorbed molecules. Simulated coverage dependences of the sticking probability and of the atomic binding energy are in reasonable agreement with experimental data. From comparison with the results of the previous work, it appears that the binding energy profiles for Ni(1 1 1) and Ni(1 0 0) have similar shapes, although at any coverage the absolute values of the oxygen binding energy are higher for the (1 0 0) surface. For metals other than Ni, particularly Pt, the model projections were found to be too parameter-dependent and therefore less certain. In such cases further model developments are needed, and we briefly comment on this situation.     

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

First-principles pseudo-potential calculations within density-functional theory framework are performed in order to study the structural and electronic properties of nickel adsorption and diffusion on a GaN(0 0 0 1)-2×2 surface. The adsorption energies and potential energy surfaces are investigated for a Ni adatom on the Ga-terminated (0 0 0 1) surface of GaN. This surface is also used to study the effect of the nickel surface coverage. The results show that the most stable positions of a Ni adatom on GaN(0 0 0 1) are at the H₃ sites and T₄ sites, for low and high Ni coverage respectively. In addition, confirming previous experimental results, we have found that the growth of Ni monolayers on the GaN(0 0 0 1) surface is possible.     

Self-assembly of 2,6-dimethylpyridine on Cu(1 1 0) directed by weak hydrogen bonding

Sequential stages of formation of a self-assembled monolayer of flat-lying 2,6-dimethylpyridine molecules on a single crystal Cu(1 1 0) surface have been observed by low-temperature scanning tunneling microscopy (LT-STM). At an adsorption temperature of 10 K, all of the molecules are randomly distributed at low coverage upon adsorption. The isolated molecules align their molecular axes parallel to the 〈0 0 1〉 azimuth of the Cu lattice. The nitrogen atom in the molecule is located at the four-fold hollow site. Upon annealing to 100 K, the molecules associate to form head-to-head dimers. The dimer units involve a pair of weak hydrogen bonds between methyl group-hydrogen atoms and N moieties on adjacent molecules, forming a core structure for further growth. In a later stage of self-assembly, single head-to-tail weak hydrogen bonds between ring C-H bonds and N moieties form in chains on the periphery of the central cores, leading to larger domains with a c(6 × 2) overlayer structure.