Self-diffusion of small clusters on fcc metal (111) surfaces

We use ab initio density-functional theory supplemented with the embedded-atom method to study the self-diffusion of small clusters on the (111) surface of eight fcc metals. A zigzag motion is found to be important in the dimer and tetramer diffusions. The dimer diffuses by a zigzag and concerted motion. The trimer diffuses by a concerted three-atom motion. The tetramer diffuses through a zigzag motion where only two atoms move simultaneously in each step. Thus, instead of increasing, the migration energy is lowered (or stays constant) for the tetramer as compared to that for the trimer. This novel break of the upwards trend in migration energy is predicted to be a general phenomenon.     

Adsorption sites of111In and111mCd on Pd(111) surfaces characterized by the electric field gradient

PAC measurements of the electric field gradient at the nucleus of isolated probe atoms on Pd(111) surfaces lead to the identification of five adsorption sites successively occupied by the parents¹¹¹In and^111mCd during annealing between 80K and 600K. The data are consistent with an estimated activation energy for Pd surface self diffusion of E_a=0.76(8)eV.     

Surface interactions of Au(I) cyclo-trimer with Au(111) and Al(111) surfaces: A computational study

A plane-wave density functional theory (DFT) study on surface interactions of a cyclo-[Au(μ-Pz)]₃ monolayer (denoted as T), Pz = pyrazolate, with Au(111) and Al(111) surfaces (denoted as M′) has been performed. Structural and electronic properties at the M′–T interfaces are determined from individually optimized structures of M′, T and M′–T. Results show that the gold pyrazolate trimer (T) binds more strongly on the Au(111) surface than on Al(111). Charge redistribution has been observed at both M′–T interfaces, where charge is “pushed” back towards the Au(111) surface from the trimer monolayer in Au(111)–T system, while the opposite happens in the Al(111)–T system where the charge is being pushed toward the trimer monolayer from the Al(111) surface. Considerable changes to the work function of Au(111) and Al(111) surfaces upon the trimer adsorption which arise from monolayer vacuum level shifts and dipole formation at the interfaces are calculated. The interaction between cyclo-[Au(μ-Pz)]₃ with metal surfaces causes band broadening of the gold pyrazolate trimer in M′–T systems. The present study aids better understanding of the role of intermolecular interactions, bond dipoles, energy-level alignment and electronic coupling at the interface of metal electrodes and organometallic semiconductor to help design metal–organic field effect transistors (MOFETs) and other organometallic electronic devices.     

Diffusion of atoms on an fcc(111) surface

The motion of adatoms on an fcc(111) surface is modeled. In light of recently revealed specific features of the potential energy surface, a new empirical potential taking into account the nonequivalence of fcc and hcp adsorption sites was proposed. The modeling was performed by the method of stochastic cycles, a new variety of MDM. Calculations were performed for the diffusion of oxygen atoms over a platinum surface. The potential energy surface was constructed based on experimental data and DFT calculations. An analysis demonstrated that the diffusion activation energy coincides with the static barrier while the preexponential factor is weakly temperature-dependent. Note, however, that the degree of asymmetry of the potential produces a significant effect on the preexponential factor: the changeover from a symmetric two-well potential to a strongly asymmetric single-well potential causes an increase in the preexponential factor by an order of magnitude.     

室温下Ge(111)表面原子的协调运动

用扫描隧道显微镜观察了Ge(111)表面的由快速降温至室温而保留下来的中温相.发现除了单个Ge增原子的扩散运动外,主要的运动是发生在畴界上的Ge增原子列的开环和闭环协调运动.从实验上证实了这种最近引起很大实验和理论兴趣的运动方式.通过观察又提出了有待进一步实验和理论研究的问题.     

Diffusion and pair interactions of CO molecules on Pd(111)

The diffusion and interactions of CO molecules on Pd(111) were studied by scanning tunneling microscopy. By following the random walk motion of individual molecules as a function of temperature, an activation energy barrier for diffusion of 118 +/- 5 meV was determined. The interaction between CO molecules was found to be repulsive for pairs separated by one or two Pd(111) lattice distances, and weakly attractive at a separation of sqrt[3].     

Surface mass transport of Ag and In on vicinal Si(111)

Mass transport of Ag and In on vicinal Si(111) has been investigated by scanning Auger microscopy (SAM). Highly anisotropic surface diffusion and surface electromigration due to direct current were observed for Ag and In adatoms on 0°−, 0.5°−, 3°− and 6°−off vicinal Si(111) surfaces. The diffusion on the intermediate layer is strongly enhanced in the direction parallel to the step edge for Ag adatoms, while it is remarkably suppressed in the direction perpendicular to the step edge for In adatoms. The activation energy of the diffusion for the Ag adatoms ranged between 0.81 and 1.3 eV, while that for In adatoms increased from 0.31 to 0.66 eV with increasing the vicinal angle. The anisotropic diffusion transport is explained in terms of the step structure and the difference in the binding energy at the step site and the terrace site.     

Initial stages in the Sm-Si(111) interface formation

The initial stages in the formation of the Sm-Si(111) interface have been studied by thermal desorption, atomic beam modulation, and low-energy-electron-diffraction spectroscopy. The structure of adsorbed films and samarium silicide films, as well as the Sm atom desorption kinetics have been investigated within a broad range of surface coverages and temperatures. The activation energy of desorption from the thermally most stable 3×2 submonolayer structure, as well as the binding energy of a single samarium atom with the substrate, have been measured. The temperature of the onset of silicide decomposition and the activation energy of this process have been determined. It is shown that the Sm-Si(111) interface forms by a mechanism close to that of Stransky-Krastanov. Fiz. Tverd. Tela (St. Petersburg) 40, 371–378 (February 1998)     

Adsorption of water monomer and clusters on platinum(111) terrace and related steps and kinks II. Surface diffusion

Surface diffusion of water monomer, dimer, and trimer on the (111) terrace, (221) and (322) stepped, and (763) and (854) kinked surfaces of platinum was studied by density functional theory using the PW91 approximation to the energy functional. Monomer diffusion on the terrace is facile, with an activation barrier of 0.20 eV, while dimer and trimer diffusions are restricted due to their high activation barriers of 0.43 and 0.48 eV, respectively. During monomer diffusion on the terrace the O–Pt distance increases by 0.54 Å, about 23% of the initial distance of 2.34 Å. The calculated rate of monomer diffusion hops is in good agreement with the onset temperature of diffusion measurements of Daschbach et al., J. Chem. Phys., 120 (2004) 1516. Alternative monomer diffusion pathways, in which the molecule rolls or flips, were also found. These pathways have diffusion barriers of 0.22 eV. During dimer diffusion on the terrace, the donor molecule rises 0.4 Å at the saddle point, while the acceptor rises by only 0.03 Å. Monomer diffusion up to steps and kinks, with activation barriers of 0.11–0.13 eV, facilitate chain formation on top of step edges. The energy landscape of monomer diffusion from terrace to step to kink sites is downhill with a maximum activation barrier of 0.26 eV. A model for water adsorption is presented in which monomer diffusion leads to concurrent formation of terrace clusters and population of steps/kinks, the latter consistent with the STM measurements of Morgenstern et al., Phys. Rev. Lett., 77 (1996) 703.     

STM-mediated atom motion: a Co atom and mixed CoCu<Subscript>n</Subscript> chains on a Cu(111) surface

Performing atomic scale simulations, we study the effect of the scanning tunneling microscopy tip on atom motion on a metal surface at zero bias voltage. We concentrate on a Co atom and mixed CoCu _n (n ? 68) chains on a Cu(111) surface. It is revealed that the atom motion can be tuned by adjusting the tip-substrate distance. The change in the potential landscape induced by the tip is found to depend on the tip height. In the presence of the tip, the Co atom can freely jump from the fcc site to the hcp site or vice versa when putting the tip above the adatom at a certain height. For the mixed CoCu _n chains on the Cu(111) surface, the diffusion barrier of the end Co atom from the fcc site to the nearby hcp site increases with the increasing chain length and reaches the limit when the chain length is beyond CoCu₇ without the tip. Especially, the short chains can perform a collective motion with the help of the tip. The importance of the relaxation induced by the tip-adatom interaction is demonstrated.     

The mechanism of H2 dissociation and adsorption on Mn-modified Ni(111) surface: A density functional theory-based investigation

The mechanism of H₂ dissociative adsorption on Mn-modified Ni(111) surface is investigated and explained using spin-polarized density functional theory (DFT). Potential energy surface (PES) is used to determine the efficient reaction pathway of H₂ on the surface. The dissociative adsorption of H₂ in the hollow sites with its center-of-mass (CM) positioned on top of Ni atom has low activation barrier. This is lower compared if its CM is on top of the Mn atom. The difference in the reactivity of H₂ with Ni and Mn as the CM is corroborated by the positions of the bonding and antibonding orbitals of H₂ as it approaches the surface which is verified from local density of states (LDOS). The greater density of states in the region around the Fermi level of the d_zz, d_xz, and d_yz orbitals of the Ni atom explains the low activation barrier obtained for the dissociation of H₂ on top of the Ni atom in the Mn-modified Ni(111) surface.     

A study of the (00) LEED beam intensity at normal incidence from CdS(0001), Cu(001), Cu(111), and Ni(111)

A Faraday cage apparatus is used for the measurement of the (00) LEED beam intensity, I₍₀₀₎, and the total secondary emission coefficient, δ(E_k), for angles of incidence from 0° ± 2° to 8° ± 2°, with an energy resolution of ± 0.037 of the incident beam energy, in the energy range 1 to 200 eV. The data are normalized and expressed as a fraction of the incident beam intensity. The basic principle of operation is the separation of the incident and specularly diffracted beams in a uniform magnetic field. Monolayer, or in-plane, resonances associated with the emergence of nonspecular beams, as well as beam threshold minima, are observed in I₍₀₀₎ at normal incidence from clean CdS(0001), Cu(111), and Ni(111). Some major differences are observed in the I₍₀₀₎ profiles for the clean (111) surfaces of nickel and copper. All secondary Bragg peaks, except the 2$\text{2}\text{3}$ order, have greater intensities for Ni(111) in the energy range 50–150 eV, thus indicating that the atomic scattering cross-section for electrons in this energy range is larger for nickel than for copper. For the (111) surface of nickel, the (11) resonance is missing, but the (10) resonance and all $\text{1}\text{3}$ order secondary Bragg peaks between the second and fifth orders are observed. For Cu(111) both the (10) and (11) resonances are observed, but the $\text{1}\text{3}$, $\text{2}\text{3}$, 1$\text{2}\text{3}$, and 3$\text{1}\text{3}$ order secondary Bragg peaks are missing in this energy range. These data indicate that multiple scattering with evanescent intermediate waves, or “shadowing”, is predominate on the (111) surfaces on nickel and copper for energies above 30 eV, and that below 30 eV multiple scattering with propagating intermediate waves is predominate on Cu(111). Correlation of the (00) beam intensity profiles from clean Ni(111) at 0°, 2°, and 6° with the intensity profiles of the (10). (1&#x0304;0), and (11) non-specular beams is nearly one-to-one from 30 eV to 100 eV, thus supporting the dynamical theories of LEED in which peaks in the (00) beam are expected to occur at nearly the same energies as peaks in the non-specular beams.     

Atomic and molecular adsorption on Pd(111)

The adsorption properties of a variety of atomic species (H, O, N, S, and C), molecular species (N₂, HCN, CO, NO, and NH₃) and molecular fragments (CN, NH₂, NH, CH₃, CH₂, CH, HNO, NOH, and OH) are calculated on the (111) facet of palladium using periodic self-consistent density functional theory (DFT–GGA) calculations at ¼ ML coverage. For each species, we determine the optimal binding geometry and corresponding binding energy. The vibrational frequencies of these adsorbed species are calculated and are found to be in good agreement with experimental values that have been reported in literature. From the binding energies, we calculate potential energy surfaces for the decomposition of NO, CO, N₂, NH₃, and CH₄ on Pd(111), showing that only the decomposition of NO is thermochemically preferred to its molecular desorption.     

Nickel carbonyl adsorption on the Ni(111) surface

Overlayers formed by the adsorption of Ni(CO)₄ in CO on the Ni(111) surface at 100 K were characterized using high resolution electron energy loss spectroscopy and thermal desorption spectroscopy. At temperatures below 135 K, molecular nickel carbonyl adsorbs on the CO saturated Ni(111) surface as suggested by several observations. Vibrational transitions characteristic of molecular Ni(CO)₄ are dominant. The energy dependence of both the elastic and inelastic electron scattering cross sections are dramatically altered by Ni(CO)₄ adsorption. All of the mass spectrometer ionization fragments typical of molecular Ni(CO)₄ are observed in the narrow thermal desorption peak at 150 K. The inelastic scattering cross sections for both adsorbed nickel carbonyl and adsorbed CO on the Ni(111) surface suggest that a nonresonant dipole scattering mechanism is dominant.     

Diffusion processes of trimers on missing row surfaces: \mathrm{Cu }_{3}/\mathrm{Ag }(110) and \mathrm{Ag }_{3}/\mathrm{Cu }(110)

A semi-empirical potential according to the embedded atom, has been applied to investigate the diffusion of trimers by computing the energy barriers for different mechanisms. Our attention was more focused on the leapfrog process which is likely to occur on missing row surfaces. The activation barriers of this mechanism are calculated using drag method at 0K. These barriers are found to be 0.64 and 0.68 eV for hopping out the channel for $\mathrm{Cu }_{3}/\mathrm{Ag }(110) \mathrm{and Ag }_{3}/\mathrm{Cu }$ (110) respectively. While for hopping down at the other side they are about 0.42 and 0.32 eV. Moreover, a deep metastable position is observed during leapfrog diffusion leading to some spectacular trimer motion. At high temperature and essentially for $\mathrm{Cu }_{3}/\mathrm{Ag }$ (110), we also observed a competition between leapfrog process and concerted jump mechanism with a deformation of trimer geometry. Implications of these findings are briefly discussed.     

Reactive molecular dynamic simulations of hydrocarbon dissociations on Ni(111) surfaces

Empirical potential parameters for H, C and Ni elements have been developed for the ReaxFF force field in order to study the decomposition of small hydrocarbon molecules on nickel using molecular dynamics simulations. These parameters were optimized using the geometrical and energetic information obtained from density functional (DFT) calculations on a subset of hydrogen and methane reactions with nickel (111) surfaces. The resulting force field was then used to obtain a molecular perspective of the dynamics of the methane dissociative adsorption on Ni(111) as well as two other small alkane molecules, ethane and n-butane. NVT simulations of dissociative adsorption of methane over a range of temperatures enabled the estimation of the sticking coefficient for the adsorption as well as the activation energy of the first C–H bond breaking. The rate constants of each elementary step (both forward and reverse) of CH_x dissociation on Ni(111) were obtained by monitoring the surface species and a microkinetic model was constructed as a result. Qualitative analyses of the simulations of ethane and n-butane decompositions on Ni(111) demonstrate that such reactive MD technique can also be used to obtain useful information on complex reaction networks.     

Initial stages of formation of a Yb-Si(111) interface

The initial stages of formation of an Yb-Si(111) interface are investigated by several methods: thermal desorption spectroscopy, atomic beam modulation, and low-energy electron diffraction. The structure of the adsorbed films and ytterbium silicide films is analyzed over a wide range of surface coverage ratios, along with the desorption kinetics of Yb atoms. The desorption activation energies of Yb atoms are measured for 3×2, 5×1, and 2×1 submonolayer structures. The temperature interval in which ytterbium silicide decomposes and the activation energy of this process are determined. It is shown that the Yb-Si(111) phase interface evolves by a mechanism similar to the Stransky-Krastanov mechanism. Fiz. Tverd. Tela (St. Petersburg) 39, 256–263 (February 1997)     

Molecular orbital study of co chemisorption and oxidation on a Pt(111) surface

A molecular orbital study was made, using an atom superposition and electron delocalization (ASED) technique, of the structures and energy levels of CO on Pt(111) surface. CO is predicted to be preferentially adsorbed at a height of 2.05 Å from the surface at a 1-fold position with the carbon end down. The calculated binding energy (1.7 eV) is in good agreement with the recent experimental result (1.5 eV) of Campbell et al. Calculated binding energies for bridging (1.3 eV) and high coordinate (1.1 eV) sites are predicted to be smaller in magnitude. Calculated results are used to discuss the ordering of energy levels of adsorbed CO. The interaction between CO (adsorbed) and O (adsorbed) has been studied to estimate the energy of activation for the oxidation of CO on Pt(111) surface. The calculated activation energy (1.6 eV) is in reasonable agreement with the recent experimental result (1 eV) of Campbell et al. The Langmuir-Hinshelwood mechanism is found to be favored. We predict CO₂ bonds vertically.     

Mechanistic study of the CO oxidation reaction on the CuO (111) surface during chemical looping combustion

Cu-based oxides oxygen carriers and catalysts are found to exhibit attractive activity for CO oxidation, but the dispute with respect to the reaction mechanism of CO and O₂ on the CuO surface still remains. This work reports the kinetic study of CO oxidation on the CuO (111) surface by considering the adsorption, reaction and desorption processes based on density functional theory calculations with dispersion correction (DFT-D). The Eley–Rideal (ER) CO oxidation mechanism was found to be more feasible than the Mars-van-Krevelen (MvK) and Langmuir–Hinshelwood (LH) mechanisms, which is quite different from previous knowledge. The energy barrier of ER, LH, and MvK mechanisms are 0.557, 0.965, and 0.999 eV respectively at 0 K. The energy barrier of CO reaction with the adsorbed O species on the surface is as low as 0.106 eV, which is much more active in reacting with CO molecules than the lattice O of CuO (111) surface (0.999 eV). A comparison with the catalytic activity of the perfect Cu₂O (111) surface shows that the ER mechanism dictates both the perfect Cu₂O (111) and the CuO (111) surface activity for CO oxidation. The activity of the perfect Cu₂O (111) surface is higher than that of the perfect CuO (111) surface at elevated temperatures. A micro-kinetic model of CO oxidation on the perfect CuO (111) surface is established by providing the rate constants of elementary reaction steps in the Arrhenius form, which could be helpful for the modeling work of CO catalytic oxidation.     

Magnetic properties of (Co Pd ) superstructures on Pd(100) and Pd(111)

The magnetic properties of (Co_nPd_m)_r superstructures on Pd(100) and Pd(111) are evaluated using the fully-relativistic spin-polarized screened Korringa-Kohn-Rostoker method. It is found that only in the case of a Pd(111) substrate such superstructures exhibit perpendicular magnetism, while on a Pd(100) substrate the magnetization is oriented in-plane. Also investigated is the effect of interdiffusion in repeated superstructures. By using the inhomogeneous coherent potential approximation (CPA) for layered systems the effect of ordering into (repeated) superstructures can be described in an ab-initio-like manner. It is found that already small amounts of interdiffusion can be decisive for the actual value of the magnetic anisotropy energy. Received 3 November 1999 and Received in final form 18 January 2000