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.     

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.     

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.     

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

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.     

Adsorption and dissociation of O2 on CuCl(1 1 1) surface: A density functional theory study

The adsorption and dissociation of O₂ on CuCl(1 1 1) surface have been systematically studied by the density functional theory (DFT) slab calculations. Different kinds of possible modes of atomic O and molecular O₂ adsorbed on CuCl(1 1 1) surface and possible dissociation pathways are identified, and the optimized geometry, adsorption energy, vibrational frequency and Mulliken charge are obtained. The calculated results show that the favorable adsorption occurs at hollow site for O atom, and molecular O₂ lying flatly on the surface with one O atom binding with top Cu atom is the most stable adsorption configuration. The O-O stretching vibrational frequencies are significantly red-shifted, and the charges transferred from CuCl to oxygen. Upon O₂ adsorption, the oxygen species adsorbed on CuCl(1 1 1) surface mainly shows the characteristic of the superoxo (O₂⁻), which primarily contributes to improving the catalytic activity of CuCl, meanwhile, a small quantity of O₂ dissociation into atomic O also occur, which need to overcome very large activation barrier. Our results can provide some microscopic information for the catalytic mechanism of DMC synthesis over CuCl catalyst from oxidative carbonylation of methanol.     

The local structure of SO2 and SO3 on Ni(1 1 1): A scanned-energy mode photoelectron diffraction study

O 1s and S 2p scanned-energy mode photoelectron diffraction (PhD) data, combined with multiple-scattering simulations, have been used to determine the local adsorption geometry of the SO₂ and SO₃ species on a Ni(1 1 1) surface. For SO₂, the application of reasonable constraints on the molecular conformation used in the simulations leads to the conclusion that the molecule is centred over hollow sites on the surface, with the molecular plane essentially parallel to the surface, and with both S and O atoms offset from atop sites by almost the same distance of 0.65 Å. For SO₃, the results are consistent with earlier work which concluded that surface bonding is through the O atoms, with the S atom higher above the surface and the molecular symmetry axis almost perpendicular to the surface. Based on the O 1s PhD data alone, three local adsorption geometries are comparably acceptable, but only one of these is consistent with the results of an earlier normal-incidence X-ray standing wave (NIXSW) study. This optimised structural model differs somewhat from that originally proposed in the NIXSW investigation.     

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.     

Adsorption dynamics of O2 on Cu(1 0 0)

We have studied the adsorption of O₂ on the Cu(1 0 0) surface using both static potential energy surface (PES) calculations and ab initio molecular dynamics. The dynamical calculations complement the PES results, revealing steering effects which could not be predicted based on the static calculations only. We study the effect of oxidation and Ag doping on O₂ adsorption dynamics. The results are discussed in the light of recent molecular beam experiments.     

Theoretical study of CO adsorption on the illuminated TiO2 (0 0 1) surface

A quantum modeling of the CO adsorption on illuminated anatase TiO₂ (0 0 1) is presented. The calculated adsorption energy and geometries of illuminated case are compared with the ground state case. The calculations were achieved by using DFT formalism and the BH and HLYP. Upon photoexcitation, an electron-hole pair is generated. Comparing of natural population in the ground state and the exited state, shows that an electron is trapped in a Ti⁴⁺ ion and a hole is localized in an oxygen ion. The photoelectron helps generation of a CO₂ molecule on the TiO₂ surface. As shown by optimization of these systems, the CO molecule adsorbed vertically on the TiO₂ (0 0 1) surface in the ground state case while the CO molecule made an angle of 134.3° to this surface at the excited state case. Based on the here used model the obtained adsorption energy was 0.36 eV which is in excellent agreement with the reported experimental value. In the present work the C-O stretch IR frequencies are calculated which are 1366.53 and 1423.16 cm⁻¹. These results are in good agreement with the earlier reported works for the surface carbonaceous compounds, and oxygenated carbon species.     

Reaction of SO2 with Cu/TiO2(1 1 0): Effects of size and metal-oxide interactions in the chemical activation of Cu nanoparticles

The reaction of SO₂ with stoichiometric TiO₂(1 1 0), partially reduced TiO_2 − x(1 1 0) and Cu/TiO₂(1 1 0) was investigated using synchrotron based X-ray photoemission spectroscopy. SO₂ adsorbs on perfect TiO₂(1 1 0) forming SO₄ species at room temperature, while SO₂ dissociatively adsorbs on partially reduced TiO_2 − x(1 1 0) forming SO₄, SO₃ as well as two sulfide species. SO₂ exposure to Cu particles supported on perfect TiO₂(1 1 0) can lead to the formation of SO₄, SO₃ and sulfide species. When depositing Cu on SO₄/TiO₂(1 1 0) at room temperature, the dissociation efficiency of Cu atoms is much higher than that of Cu deposited on TiO₂(1 1 0) prior to SO₂ dosing. The post-deposited Cu atoms can efficiently contact and react with SO₄ species before they form Cu-Cu bonds and big clusters. Small Cu nanoparticles supported on TiO₂(1 1 0) are more reactive towards SO₂ than surfaces of bulk copper. The chemical reactivity of the Cu/TiO₂(1 1 0) system increases with Cu coverage until reaching a maximum at θ_Cu = 0.5-0.8 ML. After this point, an increase in Cu coverage leads to the formation of big Cu particles and the reactivity of the system decrease to that typical of bulk Cu. A comparison with results for SO₂/Cu/MgO(0 0 1) indicates that the effects of size and metal ↔ oxide interactions are important for the chemical activation of Cu nanoparticles on titania.     

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.     

Adsorption of SO2 on Li atoms deposited on MgO (1 0 0) surface: DFT calculations

The adsorption of sulfur dioxide molecule (SO₂) on Li atom deposited on the surfaces of metal oxide MgO (1 0 0) on both anionic and defect (F_s-center) sites located on various geometrical defects (terrace, edge and corner) has been studied using density functional theory (DFT) in combination with embedded cluster model. The adsorption energy (E_ads) of SO₂ molecule (S-atom down as well as O-atom down) in different positions on both of O⁻² and F_s sites is considered. The spin density (SD) distribution due to the presence of Li atom is discussed. The geometrical optimizations have been done for the additive materials and MgO substrate surfaces (terrace, edge and corner). The oxygen vacancy formation energies have been evaluated for MgO substrate surfaces. The ionization potential (IP) for defect free and defect containing of the MgO surfaces has been calculated. The adsorption properties of SO₂ are analyzed in terms of the E_ads, the electron donation (basicity), the elongation of S-O bond length and the atomic charges on adsorbed materials. The presence of the Li atom increases the catalytic effect of the anionic O⁻² site of MgO substrate surfaces (converted from physisorption to chemisorption). On the other hand, the presence of the Li atom decreases the catalytic effect of the F_s-site of MgO substrate surfaces. Generally, the SO₂ molecule is strongly adsorbed (chemisorption) on the MgO substrate surfaces containing F_s-center.     

CO2 dissociation on Ni(2 1 1)

The direct and H-mediated dissociation of CO₂ on Ni(2 1 1) were investigated at the level of density functional theory. Although formate (HCOO) formation via CO₂ hydrogenation was widely reported for CO₂ adsorption on metal surfaces, it is found that on Ni(2 1 1) HCOO dissociation into CHO and O is much difficult, while direct dissociation of adsorbed CO₂ into CO and O is more favorable. It is also found that the degree of electron transfer from surface to adsorbed CO₂ correlates with the elongation of C-O bond lengths and the reduction of the CO₂ dissociation barrier.     

Adsorption of O2 on the CaO and SrO (1 0 0) surfaces: A first-principles study

First-principles calculations based on density functional theory (DFT) have been performed to investigate the adsorption of O₂ on the CaO and SrO (1 0 0) surfaces. The present results indicate that the bridge-top site for both the CaO and SrO (1 0 0) surfaces is the most favorable site for O₂ adsorption, with predicted adsorption energies of 1.437 eV and 1.236 eV, respectively. Detailed analysis of density of states, Mullikan population and vibration frequency are performed. The calculated results imply the possible formation of a peroxo (O₂²⁻) when O₂ is adsorbed with the bridge-top mode on both CaO and SrO (1 0 0) surfaces.     

First principle study of SrTiO3 (0 0 1) surface and adsorption of NO on SrTiO3 (0 0 1)

The results of first-principles calculations about the two possible terminations of (0 0 1) surface of SrTiO₃ perovskite and adsorption of NO on SrTiO₃ (0 0 1) surface were presented. Both surface parameters (atomic structures and electronic configurations) and adsorption parameters (bond, energy and charge) of NO on SrTiO₃ (0 0 1) surface, which have never been investigated before as far as we know were investigated using density functional theory calculations with the local-density approximation (DFT-LDA). It was found that the two possible terminations of SrTiO₃ (0 0 1) surface have large surface relaxation, which leads to surface polarization and exhibits different reactivity toward the dissociative adsorption of NO. The electronic states of TiO₂-terminated surface have larger difference than that of bulk, so it is more favorable for adsorption of NO than SrO-terminated 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.     

Adsorption and dissociation of O2 on the Cu2O(1 1 1) surface: Thermochemistry, reaction barrier

The adsorption and dissociation of O₂ on the perfect and oxygen-deficient Cu₂O(1 1 1) surface have been systematically studied using periodic density functional calculations. Different kinds of possible modes of atomic O and molecular O₂ adsorbed on the Cu₂O(1 1 1) surface are identified: atomic O is found to prefer threefold 3Cu site on the perfect surface and O_vacancy site on the deficient surface, respectively. Cu_CUS is the most advantageous site with molecularly adsorbed O₂ lying flatly over singly coordinate Cu_CUS-Cu_CSA bridge on the perfect surface. O₂ adsorbed dissociatively on the deficient surface, which is the main dissociation pathway of O₂, and a small quantity of molecularly adsorbed O₂ has been obtained. Further, possible dissociation pathways of molecularly adsorbed O₂ on the Cu₂O(1 1 1) surface are explored, the reaction energies and relevant barriers show that a small quantity of molecularly adsorbed O₂ dissociation into two O atoms on the deficient surface is favorable both thermodynamically and kinetically in comparison with the dissociation of O₂ on the perfect surface. The calculated results suggest that the presence of oxygen vacancy exhibits a strong chemical reactivity towards the dissociation of O₂ and can obviously improve the catalytic activity of Cu₂O, which is in agreement with the experimental observation.     

Reaction mechanism for CO oxidation on Cu(3 1 1): A density functional theory study

The microscopic reaction mechanism for CO oxidation on Cu(3 1 1) surface has been investigated by means of comprehensive density functional theory (DFT) calculations. The elementary steps studied include O₂ adsorption and dissociation, dissociated O atom adsorption and diffusion, as well as CO adsorption and oxidation on the metal. Our results reveal that O₂ is considerably reactive on the Cu(3 1 1) surface and will spontaneously dissociate at several adsorption states, which process are highly dependent on the orientation and site of the adsorbed oxygen molecule. The dissociated O atom may likely diffuse via inner terrace sites or from a terrace site to a step site due to the low barriers. Furthermore, we find that the energetically most favorable site for CO molecule on Cu(3 1 1) is the step edge site. According to our calculations, the reaction barrier of CO + O → CO₂ is about 0.3 eV lower in energy than that of CO + O₂ → CO₂ + O, suggesting the former mechanism play a main role in CO oxidation on the Cu(3 1 1) surface.     

Adsorption and reactions of HN3 on Si(1 0 0)-2 × 1: A computational study

We have studied the adsorption and decomposition of HN₃ on Si(1 0 0)-2 × 1 surface using the hybrid density functional B3LYP method and effective core potential basis, LanL2DZ, with Si₁₅H₁₆ as a double dimer surface model for cluster calculations. The result shows that the barriers for the dissociative adsorption of HN₃ forming HN(a) + N₂(g) are quite low by stepwise dissociation processes occurring either on a dimer or across the dimers. The low activation energies are consistent with previous experimental observations that the molecularly adsorbed HN₃(a) can undergo decomposition producing HN(a) at low surface temperatures. On the other hand, the predicted activation energies for the N₃(a) + H(a) formation processes are all relatively higher. These results also explain the absence of the N₃(a) species in HREELS measurements following each annealing experiment. Several selected reaction paths were also confirmed with slab model calculations using an optimization approach coupling the energy and gradient calculations by the slab model with the geometrical optimization using Berny algorithm.In addition, the adsorbate effect was examined for the end-on and side-on molecular configurations. For the side-on adsorption configuration, all possible combinations with 1-4 adsorbates can exist on the four surface Si sites of the double dimers, with adsorption energies lying closely to the multiples of that of a single side-on adsorbate (LM2); i.e., adsorption energies are nearly additive. Interestingly, for the end-on adsorption, only 1 and 2 HN₃ molecules can adsorb on a dimer due to the presence of the negative charges on the remaining Si sites in the neighboring dimer. For the two end-on adsorbates on the same dimer, the total adsorption energy is close to two times that of HN₃(a) or LM1. For the mixed end-on/side-on configurations, only one of each type can co-exist on a single dimer pair (Si₁-Si₂ or Si₃-Si₄) sites with adsorption energy close to the sum of those of one end-on and one side-on adsorbate. Finally, the predicted reaction routes and vibrational frequencies showed good agreement with previous experimental results. The stabilities of many ad-species involved in these reactions with end-on and/or side-on configurations have been predicted together with the transition states connecting those species.