First-principles study of benzene on noble metal surfaces: Adsorption states and vacuum level shifts

We have studied the interaction of benzene with Cu(1 1 1), Ag(1 1 1) and Au(1 1 1) surfaces using density functional theory (DFT) within a generalized gradient approximation (GGA) and the van der Waals density functional [vdW-DF; M. Dion, H. Rydberg, E. Schröder, D.C. Langreth, B.I. Lundqvist, Phys. Rev. Lett. 92 (2004) 246401]. The adsorption energies using vdW-DF are significantly more accurate than those using GGA, while the equilibrium adsorption distances between benzene and metal substrates () calculated by both GGA and vdW-DF are almost identical. The work function changes induced by the adsorption of benzene are significantly underestimated compared with the experimental values, as a result of the overestimation of by both GGA and vdW-DF. Instead of determining the values from first-principles calculations, we deduced the most probable adsorption distances in such a way as to reproduce the experimentally-observed work function changes. The deduced adsorption distance () is shortest on Cu(1 1 1) while it is longest on Ag(1 1 1), reflecting the strength of the interactions between benzene and the metal surfaces. It turns out that the substrate dependence of the work function change is mainly ascribed to the difference in the benzene–metal distance (Z_C). Charge transfer and work function changes by the adsorption of benzene were analyzed by means of the induced density of interface states (IDIS) model [H. Vázquez, R. Qszwaldowski, P. Pou, J. Ortega, R. Pérez, F. Flores, A. Kahn, Europhys. Lett. 65 (2004) 802], and compared with the self-consistent GGA calculations. The vacuum level shifts estimated by the IDIS model agree with the GGA results for . On the other hand, the discrepancy between the two methods becomes larger for , where the back donation from the metal substrates to the adsorbate becomes significant. We show that the IDIS model reasonably works well for benzene on Cu(1 1 1), Ag(1 1 1) and Au(1 1 1) surfaces because on all surfaces. However, our analysis reveals that the actual charge density redistribution induced by the adsorption of benzene is more complicated than that assumed in the IDIS model.     

Core-level shifts of InP(1 0 0)(2 × 4) surface: Theory and experiment

Surface core-level shifts (SCLSs) of the (2 × 4)-reconstructed InP(1 0 0) surface with the established mixed In–P dimer structure have been investigated by first-principles calculations and photoelectron spectroscopy. Theoretical values were calculated using both the local density approximation (LDA) and the generalized gradient approximation (GGA) for the exchange-correlation energy functional. The obtained theoretical values are quite similar within both approximations. The found differences originate in the tiny structural differences. It is concluded that the expansion or contraction of the crystal lattice has smaller effect on the SCLSs than the geometrical details of the reconstruction, which suggests that the Madelung potential has the dominant effect on the SCLSs. The results support the presence of a P 2p peak at higher binding energy (BE) compared to bulk peak, as proposed with recent measurements [P. Laukkanen, J. Pakarinen, M. Ahola-Tuomi, M. Kuzmin, R. E. Perälä, I. J. Väyrynen, A. Tukiainen, V. Rimpiläinen, M. Pessa, M. Adell, J. Sadowski, Surf. Sci. 600 (2006) 3022], and reveal several hitherto not reported SCLSs. The calculated SCLSs reproduce the measured spectra within reasonable accuracy. Furthermore, the atomic origins of the InP(1 0 0)(2 × 4) SCLSs are solved. In particular, it is shown that the lowest SCLS of P 2p of the InP(1 0 0)(2 × 4) arises from the topmost In–P dimers.     

Adsorption of methanol and methoxy on NiAl(1 1 0) and Ni3Al(1 1 1): A DFT study

The adsorption of methanol and methoxy on NiAl(1 1 0) and Ni₃Al(1 1 1) has been investigated using density functional theory (DFT). Optimised adsorption geometries and core level shifts are presented. On both surfaces we find that methanol binds to the Al on-top site via its oxygen atom and with the C–O axis tilted away from the surface normal. Methoxy also shows a preference for Al-dominated sites. On NiAl(1 1 0), we predict that methoxy adsorbs with its oxygen atom in the Al–Al bridge site, while it is determined to be adsorbed with its oxygen atom in a 2Ni + Al hollow site on Ni₃Al(1 1 1), closer to Al than Ni. Surface and adsorbate induced binding energy shifts in the Al 2p states are calculated and found to be in good agreement with experimental high resolution photoelectron spectroscopy results.     

Comparison of sulfur interaction with hydrogen on Pt(1 1 1), Ni(1 1 1) and Pt3Ni(1 1 1) surfaces: The effect of intermetallic bonding

Adsorption strengths of hydrogen and sulfur both individually and together as co-adsorbates were investigated on Pt(1 1 1), Ni(1 1 1) and Pt₃Ni(1 1 1) surfaces using density functional theory in order to determine the effect of metal alloying on sulfur tolerance. The adsorption strengths of H and S singly follow the same trend: Ni(1 1 1) > Pt(1 1 1) > Pt₃Ni(1 1 1), which correlates well with the respective d-band center positions of each surface. We find that the main effect of alloying is to distort both the sub-layer structure and the Pt overlayer resulting in a lowered d-band. For all three surfaces, the d-band shifts downward non-linearly as a function of S coverage. Nearly identical decreases in d-band position were calculated for each surface, leading to an expectation that subsequent adsorption of H would scale with surface type similarly to single species adsorption. In contradiction to this expectation, there was no clearly discernable difference between the energies of coadsorbed H on Pt(1 1 1) and Ni(1 1 1) and only a slightly lowered energy on Pt₃Ni(1 1 1). This provides evidence that coadsorbed species in close proximity interact directly through itinerant mobile electrons and through electrostatic repulsion rather than solely through the electronic structure of the surface. The combination of the lowered d-band position (arising from distorted geometry) and direct co-adsorbate interactions on Pt₃Ni(1 1 1) leads to a lower energy barrier for H₂S formation on the surface compared to pure Pt(1 1 1). Thus, alloying Pt with Ni both decreases the likelihood of S adsorption and favors S removal through H₂S formation.     

Theoretical study of hydrogenation process of formate on clean and Zn deposited Cu(1 1 1) surfaces

We have studied the effect of Zn on hydrogenation of formate to dioxomethylene on the Cu(1 1 1) surface by using a density functional theory–generalized gradient approximation (DFT–GGA)-pseudopotential method. We show that substitutionally adsorbed Zn changes the stability of intermediate states and the activation barrier of the hydrogenation process only slightly. On the other hand, the Zn atom adsorbed on the Cu surface stabilizes all formate, transition state, and dioxomethylene relative to the gas-phase molecules. Our results support a previously proposed reaction scheme that the adsorption state of Zn changes from substitutional to on-surface adsorption during the methanol synthesis.     

Adsorption structure and spin polarization of pentacene on a magnetized Fe(1 0 0) substrate: SPMDS and ERDA study

The adsorption structure and spin-resolved electronic structure of pentacene on Fe(1 0 0) surfaces were investigated using elastic recoil detection analysis (ERDA) and spin-polarized metastable deexcitation spectroscopy (SPMDS), respectively. It was found that the pentacene molecule adsorbs with its molecular plane parallel to the surface in the initial stage of growth, while the molecular plane tilted from the surface with the increase of the coverage. In the pentacene-Fe surface interaction, a donation and backdonation interaction was indicated. It was found that the polarity of the spin polarization induced in the pentacene molecular orbitals was opposite to each other between the orbitals contributing to the donation and those contributing to the backdonation.     

Entropy-driven adsorption of carbon nanotubes on (0 0 1) and (1 1 1) surfaces of CeO2 islands grown on sapphire substrate

According to the aim to compose combinatorial material by adsorption of carbon nanotubes onto the structured CeO₂ surface the interaction of the armchair (5,5) and zigzag (8,0) nanotubes with the (0 0 1) and (1 1 1) surfaces of CeO₂ islands have been investigated by theoretical methods. The thermodynamics of the adsorption were studied at the low surface coverage region. The interaction energy between the nanotube and the different CeO₂ surfaces shows significant increase when the size of the interface reaches 7–8 unit cells of CeO₂ and it remains unchanged in the larger interface region. However, the entropy term of the adsorption is significantly high when the distances of CeO₂ islands are equal to 27 nm (adsorption of armchair (5,5) nanotube) or 32 nm (adsorption of zigzag (8,0) nanotube). This property supports adsorption of nanotubes onto CeO₂ surfaces which possesses a very specific surface morphology. A long-wave vibration of nanotubes was identified as background of this unexpected phenomenon. This observation could be applicable in the development of such procedures where the nanotube adsorption parallel to the surface is aimed to perform.     

The chemisorption of pentacene on Si(0 0 1)-2 × 1

Adsorption structures of the pentacene (C₂₂H₁₄) molecule on the clean Si(0 0 1)-2 × 1 surface were investigated by scanning tunneling microscopy (STM) in conjunction with density functional theory calculations and STM image simulations. The pentacene molecules were found to adsorb on four major sites and four minor sites. The adsorption structures of the pentacene molecules at the four major sites were determined by comparison between the experimental and the simulated STM images. Three out of the four theoretically identified adsorption structures are different from the previously proposed adsorption structures. They involve six to eight Si-C covalent chemical bonds. The adsorption energies of the major four structures are calculated to be in the range 67-128 kcal/mol. It was also found that the pentacene molecule hardly hopped on the surface when applying pulse bias voltages on the molecule, but was mostly decomposed.     

Density functional calculations for Ni adsorption on Al(1 1 0)

The adsorption of 0.25, 0.5 and 1 monolayer (ML) of the transition metal Ni on the metal substrate Al(1 1 0) was studied using first-principles calculations at the level of density functional theory. The metal–metal system was analyzed with the generalized gradient approximation. Four stable atomic configurations were considered, and the optimized geometries and adsorption energies of different Ni adsorption sites on the Al(1 1 0) surface at selected levels of coverage were calculated and compared. The four-fold hollow site was determined to be the most stable adsorption site with adsorption energy of 5.101 eV at 0.25 ML, 3.874 eV at 0.5 ML and 3.665 eV at 1 ML. The adsorption energies of the four sites slightly decreased as the Ni coverage increased. Work function analysis showed that when Ni is adsorbed on the Al(1 1 0) surface, the work function decreased as the coverage increased due to depolarization. The Mulliken population and density of states were calculated to determine the charge distribution of the adsorption site, confirming that a chemisorption interaction exists between the adsorbed Ni atom and Al(1 1 0) surface atoms.     

Theoretical studies of Sm/Si(1 1 1)–n × 1 (n = 5, 7) and Sm/Si(1 1 1)–(8 × 2) structures

In this study we employ a state-of-the-art pseudopotential method to perform local density of states (LDOS) calculations of n × 1 (n = 5, 7) and (8 × 2) reconstructions induced by the adsorption of rare-earth samarium (RE) in the submonolayer range. We conducted a full comparison between images from scanning tunneling microscopy (STM) and theoretical LDOS. Images taken of both filled and empty states show the effects induced by honeycomb chains and Seiwatz chains. We conclude that LDOS calculations are consistent with the assignment of features observed experimentally by STM.     

Comparative dynamical study of the HD and isotopomers interacting with Pd(1 1 1) and Cu(1 1 0) surfaces

We study the dynamics of HD and H₂ molecules interacting with Pd(1 1 1) and Cu(1 1 0) using the classical trajectory method based on potential energy surfaces obtained from Density Functional Theory calculations. Our results predict a negligible isotopic effect on the dissociative adsorption probability on Pd(1 1 1) whereas on Cu(1 1 0), the adsorption probability for HD(ν_i=0) is slightly lower than for H₂(ν_i=0), mainly due to its lower initial vibrational zero point energy. The final rotational energy distribution of scattered HD and H₂ molecules are very similar. This shows that the asymmetric mass distribution of HD, barely affects the fraction of initial translational energy transferred to rotation during the scattering process. Our calculations point to the larger number of open rotational excitation channels for HD, as the main cause of rotational excitation probabilities larger than for H₂. The theoretical apparent rotational temperature, T_rot, of HD molecules scattered from Pd(1 1 1) at impact energy , is in good agreement with the experimental value. In contrast, for Cu(1 1 0) the theoretical T_rot is much lower than the value measured for Cu(1 0 0). Possible reasons for such a discrepancy between theory and experiments are discussed.     

Growth of pentacene on Ag(1 1 1) surface: A NEXAFS study

Thin films of pentacene (C₂₂H₁₄) have become widely used in the field of organic electronics. Here films of C₂₂H₁₄ of thickness ranging from submonolayer to multilayer were thermally deposited on Ag(1 1 1) surface. The determination of molecular geometry in pentacene films on Ag(1 1 1) studied by X-ray absorption at different stages of growth up to one monolayer is presented.XAS spectra at the C K-edge were collected as a function of the direction of the electric field at the surface. The different features of the spectra were assigned to resonances related to the various molecular unoccupied states by the comparison with the absorption coefficient of the pentacene gas phase. The transitions involving antibonding π states show a pronounced angular dependence for all the measured coverages, from submonolayer to multilayer. The spectra analysis indicates a nearly planar chemisorption of the first pentacene layer with a tilt angle of 10°.     

Interaction of atomic H and CO with Cu(1 1 0) and bimetallic Ni/Cu(1 1 0)

Co-adsorption of hydrogen and CO on Cu(1 1 0) and on a bimetallic Ni/Cu(1 1 0) surface was studied by thermal desorption spectroscopy. Hydrogen was exposed in atomic form as generated in a hot tungsten tube. The Ni/Cu surface alloy was prepared by physical vapor deposition of nickel. It turned out that extended exposure of atomic hydrogen leads not only to adsorption at surface and sub-surface sites, but also to a roughening of the Cu(1 1 0) surface, which results in a decrease of the desorption temperature for surface hydrogen. Exposure of a CO saturated Cu(1 1 0) surface to atomic H leads to a removal of the more strongly bonded on-top CO (α₁ peak) only, whereas the more weakly adsorbed CO molecules in the pseudo threefold hollow sites (α₂ peak) are hardly influenced. No reaction between CO and H could be observed. The modification of the Cu(1 1 0) surface with Ni has a strong influence on CO adsorption, leading to three new, distinct desorption peaks, but has little influence on hydrogen desorption. Co-adsorption of H and CO on the Ni/Cu(1 1 0) bimetallic surface leads to desorption of CO and H₂ in the same temperature regime, but again no reaction between the two species is observed.     

Observation of hydrogen adsorption on 6H-SiC(0 0 0 1) surface

We have used coaxial impact-collision ion scattering spectroscopy (CAICISS) and time-of-flight elastic recoil detection analysis (TOF-ERDA) to investigate the adsorption of atomic hydrogen on the 6H-SiC(0 0 0 1)√3×√3 surface. It has been found that the saturation coverage of hydrogen on the 6H-SiC(0 0 0 1)√3×√3 surface is about 1.7 ML. Upon saturated adsorption of atomic hydrogen, the √3×√3 surface structure changes to the 1×1 structure. The data of the CAICISS measurements have indicated that as a result of the hydrogen adsorption, Si adatoms on the √3×√3 surface move from T₄ to on-top sites.     

Adsorption of methanol on Ni3Al(1 1 1) and NiAl(1 1 0): A high resolution PES study

The adsorption of methanol on Ni₃Al(1 1 1) and NiAl(1 1 0) has been studied using high resolution photoemission spectroscopy (HR-PES) and density functional theory (DFT). Both methanol and methoxy are formed on these surfaces after the initial methanol exposure at low temperatures. Heating to 200 K leads to further formation of methoxy. On NiAl(1 1 0) two different methoxy species are observed where the first is formed upon methanol adsorption, and the other results from methanol decomposition during heating. The DFT calculations show that methanol and methoxy interacts with the Al atoms on both surfaces. Methanol is found to bond through the oxygen atom to the Al on-top site on Ni₃Al(1 1 1) and NiAl(1 1 0) with the C–O axis tilted with respect to the surface normal. On Ni₃Al(1 1 1) methoxy is situated in a 2Ni+Al hollow site, whereas on NiAl(1 1 0) the Al–Al bridge site is preferred.     

Chlorine adsorption on Si(1 1 1) studied by optical methods

The adsorption process of chlorine on Si(1 1 1) has been studied by means of real time surface differential reflectance (SDR) spectroscopy and second harmonic generation (SHG). The structure observed at 3.6 eV in SDR spectra is attributed to transitions including Si–Cl antibonding states. However, the overall feature is due to the removal of the electronic states of the clean surface. Developments of adsorption on Si adatom dangling bonds and breaking of adatom back bonds are obtained from SDR spectra and second harmonic (SH) intensity. They are well fit by the solutions of the rate equations under the assumption of adsorption of atoms without migration, and the initial sticking probability on the dangling bonds and the initial breaking probability of the back bonds are determined. Dependence of the adsorption kinetics on the carrier concentration is briefly reported.     

First-principles study of small palladium clusters on NiAl(1 1 0) alloy surface

Theoretical calculations focused on the geometry, stability, electronic and magnetic properties of small palladium clusters Pd_n (n=1–5) adsorbed on the NiAl(1 1 0) alloy surface were carried out within the framework of density functional theory (DFT). In agreement with the experimental observations, both Ni-bridge and Al-bridge sites are preferential for the adsorption of single palladium atom, with an adsorption energy difference of 0.04 eV. Among the possible structures considered for Pd_n (n=1–5) clusters adsorbed on NiAl(1 1 0) surface, Pd atoms tend to form one-dimensional (1D) chain structure at low coverage (from Pd₁ to Pd₃) and two-dimensional (2D) structures are more stable than three-dimensional (3D) structures for Pd₄ and Pd₅. Furthermore, metal-substrate bonding prevails over metal–metal bonding for Pd cluster adsorbed on NiAl(1 1 0) surface. The density of states for Pd atoms of Pd/NiAl(1 1 0) system are strongly affected by their chemical environment. The magnetic feature emerged upon the adsorption of Pd clusters on NiAl(1 1 0) surface was due to the charge transfer between Pd atoms and the substrate. These findings may shade light on the understanding of the growth of Pd metal clusters on alloy surface and the construction of nanoscale devices.     

Adsorption and thermal decomposition of N-methylaniline on Pt(1 1 1)

The adsorption and thermal decomposition of N-methylaniline (NMA) on the Pt(1 1 1) surface has been studied with reflection absorption infrared spectroscopy (RAIRS), temperature programmed desorption (TPD), and X-ray photoelectron spectroscopy (XPS). NMA adsorbs molecularly at 85 K through the nitrogen lone pair and is stable up to 300 K. At temperatures of 300–350 K it converts to two or more surface intermediates including the N-methyleneaniline (NMEA) species. This NMEA intermediate dissociates upon annealing to 450 K, and further annealing leads to the desorption of HCN and H₂, leaving only C on the surface at 800 K.     

Observation of charge transfer states of F4-TCNQ on the 2-methylpropene chemisorbed Si(1 0 0)(2 × 1) surface

We have investigated the valence electronic states of 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F₄-TCNQ) on the 2-methylpropene chemisorbed Si(1 0 0)(2 × 1) surface using valence photoelectron spectroscopy. Since the electron affinity of condensed F₄-TCNQ is 5.24 eV and the energy from the valence band maximum of the 2-methylpropene saturated Si(1 0 0)(2 × 1) surface to the vacuum level is 4.1 eV, spontaneous charge transfer would be expected in the present system. At sub-monolayer coverage of F₄-TCNQ, characteristic peaks are observed at 1.1 and 2.5 eV below Fermi energy. The former peak is assigned to a singly occupied affinity level, and the latter is ascribed to a relaxed highest occupied molecular orbital of adsorbed F₄-TCNQ. The work function change is increased up to +1.3 eV as a function of F₄-TCNQ coverage. These results support the occurrence of charge transfer into F₄-TCNQ on the 2-methylpropene saturated Si(1 0 0)(2 × 1) surface.     

Temperature dependent LEED study of Pb{1 1 1}

The thermal behavior of Pb{1 1 1} was studied using low-energy electron diffraction (LEED) in the temperature range 11–323 K. The surface interlayer spacings increase with temperature at about the same rate as the bulk up to 0.5 T_m, and then increase faster. The relaxation of the surface, which is larger than for other fcc {1 1 1} surfaces, is maintained in the temperature range studied. Although Pb has a larger expansion coefficient than other metals, the surface thermal expansion behavior is in line with other surfaces, and is consistent with harmonic interplanar potentials.