Density functional slab model studies of water adsorption on flat and stepped Cu surfaces

Density functional and periodic slab model calculations are performed to study adsorption of water on various Cu surfaces, focusing on monomers and dimers at the planar Cu surfaces and monomers at stepped ones. Single water molecules tend to weakly bind to atop positions with the molecular plane basically parallel to the substrate surface on the planar surfaces or the step plane on the stepped surfaces with negligible structural deformation of water. The experimental adsorption energies of water on the (1 1 1) and (1 0 0) surfaces are about twice as large as the theoretical values of monomerically adsorbed water. This phenomenon is demonstrated to be due to formation of water clusters and/or existence of surface defects. It is revealed that the most favorable hexagonal ring superstructure on Cu(1 1 0) is a four-layer structure, not the commonly accepted bi-layer configuration. We found that the adsorption energy of monomeric water correlates linearly with following quantities, respectively: the bond length and the stretching frequency of the Cu-O bond, the coordination number of the surface Cu atom, the surface work function of the clean surface and the 1b₁ MO energy shift with respect to the value in the gas phase.     

Effects of water and electric field on atomic oxygen adsorption on Pt–Co alloys

Density functional theory is used to study the effect of atomic oxygen adsorption at various coverages with and without the presence of water on ordered and Pt-segregated PtCo surfaces. The strength of O adsorption, as well as surface reconstruction effects due to the adsorbate are strongly influenced by the presence of the oxygen-philic transition metal on the surface or subsurface. At high O coverage, buckling of the Co atom on PtCo surfaces is much smaller than that of Pt on Pt(1 1 1) surfaces, and buckling of Pt atoms on Pt-skin surfaces is negligible. Also, the effect of an electric field perpendicular to the surface on adsorbed water and atomic oxygen is investigated. Spontaneous water dissociation is not found on the ordered and segregated alloy surfaces within the entire applied electric field range (−0.51 to 0.51 V/Å). Water changes orientation under strong negative fields, switching from a metal–O to a metal–H interaction, and the effect is much more pronounced in the low-coordination sites of cluster models.     

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.     

First-principles studies of H2S adsorption and dissociation on metal surfaces

Density functional theory calculations were employed to investigate the molecular and dissociative adsorption of H₂S on the closed packed surfaces of a number of important noble metals (Ag(1 1 1), Au(1 1 1) and Cu(1 1 1)) and transition metals (Ir(1 1 1), Ni(1 1 1), Pd(1 1 1) and Pt(1 1 1)). Energy minima corresponding to adsorbed states were identified with H₂S binding preferentially at the top sites. The adsorption of other S moieties (SH and S) was also examined. SH and S were found to prefer bridge sites and hollow sites, respectively. The binding of H₂S and its S-containing dissociated species is stronger on the transition metals. The elementary reactions of abstraction of H from H₂S to form a surface SH intermediate and abstraction of H from SH to form a surface S intermediate as model pathways for the dissociation of H₂S were examined. Our results suggest that H₂S decomposition on the aforementioned transition metal surfaces is more facile, both thermodynamically and kinetically.     

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.     

Surface relaxation and surface stress of 4d transition metals

Using the density functional theory formulated within the framework of the exact muffin-tin orbitals method, we present a systematic study of the top layer relaxation and surface stress of 4d transition metals. Our calculations predict layer contractions for most surfaces. We also find that the relaxations of the close packed surfaces decrease with increasing atomic number through the 4d series. We propose that the relaxation is mainly due to the reduction of the number of sp electrons in the surface layer relative to bulk. The surface stress is found to be very sensitive to the relaxation and, therefore, an accurate determination of the layer relaxation is necessary for obtaining reliable values for the surface stress. Comparing the top layer relaxations for the close packed surfaces, we see essential deviations between data derived in different ab initio calculations. At the same time, the overall trend for the present surface stress of 4d metals is in reasonable agreement with recent full-potential data.     

Influence of surface stress in the missing row reconstruction of fcc transition metals

Surface stress relief is often invoked as the driving force for surface reconstruction. The present ab initio total energy calculations show that it is not the case, at least for the missing row reconstruction undergone by the (1 1 0) surface of the face centered cubic transition metals. Instead, they show that the tendency to reconstruct or not is closely related to the sign and magnitude of the strain derivative of surface energy. Finally, they unambiguously confirm that relativistic effects are at the origin of this reconstruction.     

Theoretical study of the chemisorption of CO on bimetallic RhCu surfaces and nanoparticles

The CO interaction with bimetallic RhCu surface models representing several compositions has been studied by first principles density functional theory calculations. The analysis of the bare bimetallic clusters Rh(4s) and Cu(3s) core-level binding energies indicates that is not possible to extract information about the oxidation state of the alloy components. The present calculations predict that CO does always sit on top sites, the influence of the alloy composition on the equilibrium geometry and vibrational frequency of CO chemisorbed at a given Rh or Cu site being very small. However, there is a large difference in the structural properties corresponding to CO chemisorption above either Rh or Cu. Therefore, the absolute value of the vibrational frequency of chemisorbed CO does not permit to extract any information about the alloy composition but afford to assign the chemisorption site. Finally, the CO adsorption energy does not follow a monotonic trend with composition. The use of the Constrained Space Orbital Variation analysis permits one to firmly establish that the difference in adsorption energy for different compositions cannot be explained through differences in the σ-donation and π-backdonation mechanisms.     

Dissociative adsorption and adsorbate-induced reconstruction on Fe{2 1 1}

The products of CO, NO, O₂ and N₂ dissociation on Fe{2 1 1} have been studied by means of first-principles density functional theory. Preferred adsorption sites for adatoms C, N and O are identified, and trends in charge transfer and surface magnetism described. An experimentally observed (2 × 1) reconstruction induced by O is confirmed to be energetically stable, and a similar reconstruction induced by N is tentatively predicted. It is argued that these reconstructions may be important not only in the context of the catalytic reactivity of the Fe{2 1 1} surface, but also for the initial stages of surface nitridation and oxidation.     

CO adsorption on the Pt(1 1 1) surface: a comparison of a gradient corrected functional and a hybrid functional

The adsorption of CO on the Pt(1 1 1) surface in a pattern has been studied with the gradient corrected functional of Perdew and Wang and the B3LYP hybrid functional. A slab which is periodic in two dimensions is used to model the system. The Perdew-Wang functional incorrectly gives the fcc site as the most favorable adsorption site, in accord with a set of previous studies. The B3LYP functional gives the top site as the preferred site. This confirms results from cluster studies where it was suggested that the different splitting, dependent on the functional, between highest occupied and lowest unoccupied molecular orbital, could be the reason for this change of the adsorption site. This is supported by an analysis based on the projected density of states and the Mulliken population.     

Cluster and periodic DFT calculations of adsorption and activation of CO2 on the Cu(hkl) surfaces

The adsorption behavior and thermal activation of carbon dioxide on the Cu(1 1 1), Cu(1 0 0), and Cu(1 1 0) surfaces have been investigated by means of density functional theory calculations and cluster models and periodic slabs. According to the cluster models, the optimized results indicate that the basis set of C and O atoms has a distinct effect on the adsorption energy, but an indistinct one on the equilibrium geometry. For the CO₂/Cu(hkl) adsorption systems studied here, the final structure of adsorbed CO₂ is near linear and the preferred modes for the adsorption of CO₂ onto the Cu(1 1 1), Cu(1 0 0), and Cu(1 1 0) surfaces are the side-on adsorption at the cross bridge site with an adsorption energy of 13.06 kJ/mol, the side-on adsorption at the short bridge site (13.54 kJ/mol), and the end-on adsorption on the on-top site with C–O bonds located along the short bridge site (26.01 kJ/mol), respectively. However, the calculated adsorption energies from periodic slabs are lower as compared to the experimental data as well as the cluster model data, indicating that the periodic slab approach of generalized gradient approximation in the density function theory may be not suitable to obtain quantitative information on the interaction of CO₂ with Cu(hkl) surfaces.     

Diffusion of Pb adatom and ad-dimer on Si(1 0 0) from ab initio studies

The diffusion pathways of Pb adatoms and ad-dimers on Si(1 0 0) are investigated by first-principles calculations. Pb adatoms are found to diffuse on top of the Si(1 0 0) dimer row with an energy barrier of 0.31 eV. However, Pb dimers are energetically more stable. Pb dimers on top of the dimer row have a high energy barrier (0.95 eV) to rotate from the lowest energy configuration to the orientation parallel to the underlying Si(1 0 0) dimer row. Once the ad-dimer is oriented parallel to Si(1 0 0) dimer row, they can diffuse along the dimer row with an energy barrier of only 0.32 eV.     

Adsorption of carbon monoxide on Pt{100} surfaces: dependence of the CO stretching vibrational frequency on surface coverage

We have used the ab initio cluster model approach to study the dependence of the CO stretching frequency on CO surface coverage. We have also investigated the relative importance of the various factors that can affect the position of the CO stretching band as coverage increases. Two effects can change the CO stretching frequency: the adsorbate–adsorbate dipole coupling, which is a purely physical effect, and the changes in the 2π^* CO molecular orbitals, due to the different chemical environment at higher coverages. From our vibrational analysis, we conclude that CO–CO dipole coupling is the main cause of the upward shift of the CO stretching band when the CO coverage is increased. The population of the 2π^* CO molecular orbitals does not change at any coverage within the region considered. We have also estimated the ¹²CO–¹³CO dipole coupling, which previous studies have assumed to be weak. Our results demonstrate that the ¹²CO–¹³CO dipole coupling is indeed weak compared with the ¹²CO–¹²CO dipole coupling. At a CO surface coverage of 0.5 monolayers (ML), we have calculated a band shift of 40 cm⁻¹ to higher frequency. However, we should point out that when one ¹²CO molecule is surrounded by a ¹³CO environment, the ¹²CO stretching band shifts 10 cm⁻¹ upwards. We have also computed the heat of adsorption of CO on Pt{100}-(1×1) as a function of CO coverage. The initial heat of adsorption is calculated to be about 192 kJ mol⁻¹ and then drops to 180 kJ mol⁻¹ at 0.5 ML. These results agree quite well with recent calorimetric measurements. Besides that, we have estimated that the CO–CO interaction energy at 0.5 ML is repulsive and has a value of 5 kJ mol⁻¹.     

Electronic structures of the suboxide films formed on TiC(1 0 0) and ZrC(1 0 0) surfaces: Density functional theory studies

Density functional theory (DFT) calculations have been performed to elucidate the electronic structures of the TiO-like film on TiC(1 0 0) and the ZrO-like film on ZrC(1 0 0), which are assumed to be monolayers of suboxide films with (1 × 1) periodicity with respect to the substrate (1 0 0) surfaces. It was revealed that the electronic structures of both films were characterized by the existence of a band around 6 eV and a band around the Fermi level. The former and latter bands were mostly composed of O 2p and metal d orbitals, respectively, indicating the substantial ionic nature of the film. The calculated DOS well reproduced the previously obtained photoelectron spectra. From the inspection of the optimized structures, it was found that the both suboxide films have rippled structures; the metal and oxygen atoms are displaced vertically downward and upward, respectively, maintaining the (1 × 1) structures.     

Vacancy formation and CO adsorption on gold-doped ceria surfaces

We study the effect of gold doping on oxygen vacancy formation and CO adsorption on the (1 1 0) and (1 0 0) surfaces of ceria by using density functional theory, corrected for on-site Coulomb interactions (DFT + U). The Au dopant substitutes a Ce atom in the surface layer, leading to strong structural distortions. The formation of one oxygen vacancy near a dopant atom is energetically “downhill” while the formation of a second vacancy around the same dopant requires energy. When the surface is in equilibrium with gaseous oxygen at 1 atm and room temperature there is a 0.4 probability that no oxygen atom left the neighborhood of a dopant. This means that the sites where the dopant has not lost oxygen are very active in oxidation reactions. Above 400 K almost all dopants have an oxygen vacancy next to them and an oxidation reaction in such a system takes place by creating a second vacancy. The energy required to form a second vacancy is smaller on (1 1 0) than on (1 0 0). On the (1 1 0) surface, it is much easier to form a second vacancy on the doped surface than the first vacancy on the undoped surface. The energy required to form a second oxygen vacancy on (1 0 0) is comparable to that of forming the first vacancy on the undoped surface. Thus doping makes the (1 1 0) surface a better oxidant but it has a small effect on the oxidative power of the (1 0 0) surface. On the (1 1 0) surface CO adsorption results in formation of a carbonate-like structure, similar to the undoped surface, while on the (1 0 0) surface direct formation of CO₂ is observed, in contrast to the undoped surface. The Au dopant weakens the bond of the surrounding oxygen atoms to the oxide making it a better oxidant, facilitating CO oxidation.     

NO adsorption and diffusion on unreconstructed Pt{1 0 0} surface. A density functional theory investigation

Ab initio density functional theory was used to investigate the adsorption and diffusion of a single NO molecule on the unreconstructed Pt{1 0 0}-(1 × 1) surface. To our knowledge this is the first theoretical study of the NO diffusion activation energy on the Pt{1 0 0} surface. The most stable adsorption position for NO corresponds to the bridge site with the axis of the molecule perpendicular to the surface. The bond of the NO molecule to the surface is through the N-atom. We found that there is a low adsorption energy when the NO molecule is bonded through the O-atom and the axis is perpendicular to the surface, for the three high symmetry sites investigated. NO diffusion between bridge-hollow sites, bridge-atop sites, and hollow-atop sites was also investigated. The barrier for NO diffusion is 0.41 eV, which corresponds to the energy difference between the bridge and hollow sites. This value is around 15% of the highest adsorption energy found on this surface. NO stretch frequencies are also calculated for the three high symmetry sites investigated.     

First principles calculations of the adsorption of acetylene on the Si(0 0 1) surface at low and full coverage

The adsorption of acetylene on the Si(0 0 1)-c(2 × 4) surface at low and full coverage is studied by first principles density functional calculations using the generalized gradient approximation. For a single acetylene molecule, the most stable configuration corresponds to the di-σ site, on-top of a silicon dimer. This configuration is 0.36 eV more stable than the end-bridge site between two adjacent Si dimers. However, if there are two acetylene molecules, the paired end bridge configuration becomes the most stable. We have also studied the kinetics of the adsorption of a single acetylene molecule. Our calculations show that the reaction is barrier-free for adsorption in the di-σ configuration, while there is an energy barrier of 0.19 eV for adsorption in the end-bridge site. At monolayer coverage, the most stable configuration corresponds to acetylene molecules in the pair-end bridge configuration, in agreement with previous calculations. We have found a noticeable coverage dependence only for the end-bridge site, but not for the di-σ. Our results show that to have an accurate picture of the adsorption of acetylene on the Si(0 0 1) surface, a large unit cell is needed.     

Initial incorporation of sulfur into the Pd(1 1 1) surface: A theoretical study

Density functional theory is used to investigate the initial inclusion of sulfur into the subsurface interstitial sites of Pd(1 1 1) surface. Pure subsurface adsorption is found to be less energetically favorable than on-surface adsorption. The incorporation of sulfur into the metal becomes more favorable than continuous adsorption on the surface after a critical on-surface sulfur coverage. We find subsurface sulfur occupation to be energetically favorable after adsorption of more than half a monolayer on the surface. Occupation of subsurface sites induces a pronounced structural distortion of the Pd(1 1 1) surface. We find significant expansion of interplanar spacing between the uppermost surface metal layers and rearrangement of the S overlayer. The interplay between the energy cost due to structural distortion of Pd(1 1 1) and the energy gain due to bond formation for different structures is discussed.     

Theoretical studies on the adsorption and decomposition of H2O on Pd(1 1 1) surface

To provide information about the chemistry of water on Pd surfaces, we performed density functional slab model studies on water adsorption and decomposition at Pd(1 1 1) surface. We located transition states of a series of elementary steps and calculated activation energies and rate constants with and without quantum tunneling effect included. Water was found to weakly bind to the Pd surface. Co-adsorbed species OH and O that are derivable from H₂O stabilize the adsorbed water molecules via formation of hydrogen bonds. On the clean surface, the favorable sites are top and bridge for H₂O and OH, respectively. Calculated kinetic parameters indicate that dehydrogenation of water is unlikely on the clean regular Pd(1 1 1) surface. The barrier for the hydrogen abstraction of H₂O at the OH covered surface is approximately 0.2-0.3 eV higher than the value at the clean surface. Similar trend is computed for the hydroxyl group dissociation at H₂O or O covered surfaces. In contrast, the O-H bond breaking of water on oxygen covered Pd surfaces, H₂O_ad + O_ad → 2OH_ad, is predicted to be likely with a barrier of ∼0.3 eV. The reverse reaction, 2OH_ad → H₂O_ad + O_ad, is also found to be very feasible with a barrier of ∼0.1 eV. These results show that on oxygen-covered surfaces production of hydroxyl species is highly likely, supporting previous experimental findings.