Ab-Initio Simulation Laboratory

A simple oscillatory intensity variation in medium-energy electron diffraction found recently [Abukawa ei al., Phys. Rev. Lett. 82, 335 (1999)] was termed correlated thermal diffuse scattering (CTDS). The potential of CTDS as a direct surface structural tool has been fully explored for the Si(001)2 X 1 surface at 300 K in a very-grazing-incidence condition. Nearly 2 pi solid-angle, three-dimensional (3D) CTDS patterns were measured for an energy range of 500-2000 eV. The 3D Patterson functions obtained by Fourier inversion of the measured CTDS patterns clearly revealed the building blocks of the Si(001)2 X 1 surface, i.e., the bond orientations and lengths of the buckled Si dimers, within an accuracy of 1 degrees and 0.1 Angstrom, respectively.

We have observed simple oscillations in three-dimensional (3D) patterns of electron thermal diffuse scattering (separated from electron-electron energy loss) measured on a Si(001) surface. We interpret these oscillations as coherent interference within a small cluster of atoms in which vibrational correlation within the nearest neighbors (NN) is dominant. A 3D Patterson function analysis of the oscillation reveals the atomic structure of the Si(001) surface consisting of NN pairs including dimers. This finding provides a promising new clue to determine the structures of bulk and the surface of solids.

A comprehensive atomic model for the reconstructed surface of Si3N4 thin layer grown on Si(lll) is presented. Kikuchi electron holography images clearly show the existence of adatoms on the Si3N4(0001)/Si(111)-(8 x 8) surface. Compared with the nb initio calculations, more than 30 symmetry-inequivalent atomic pairs in the outmost layers are successfully identified. Scanning tunneling microscopy (STM) images show diamond-shaped unit cells and nine adatoms in each cell. High-resolution STM images reveal extra features and are in good agreement with the partial charge density distribution obtained from total-energy calculations.

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

Atomic structures of the reconstructed Si(113) surfaces were studied by using Kikuchi electron holography (KEH). Three-dimensional images show clearly the characteristics of the puckering model for both Si(113)(3x2) and (3x1) surfaces. The KEH results support the puckering model. Based on bur studies, the tetramers are puckering alternatively in the (3x2) surface. Whereas in (3X1) structures, there are two domains, within each of them, tetramers buckled uniformly, but the overall directions are opposite. When doped with H atoms on a (3x2) surface, the asymmetric tetramers change into symmetric form. [S0163-1829(99)51116-8].

{

Trimer is the smallest cluster that can have a one-dimensional or a two-dimensional structure on fee (111) surface. Using first-principles density-functional-theory calculations, the structural and dynamical properties of Al trimer on Al(111) surface have been studied in detail. Al trimer on Al(111) surface has four close-packed (compact) triangular configurations, two linear configurations, and some other noncompact triangular configurations. The close-packed triangular trimers are more stable than the noncompact triangular trimers as well as the linear trimers. For the dynamics of Al trimer on Al(111) surface, the diffusion processes are much more complicated than the adatom and dimer diffusions. There are three different kinds of diffusion mechanisms: concerted translations and rotation of compact triangular trimers (the energy barrier

Single differential and total cross sections of the electron-impact ionization of the hydrogen atom are calculated numerically in the two-potential distorted-wave approximation for excess energies from 0 to 1 eV. Partial-wave contributions to the cross sections are also investigated. The near-threshold cross section is parametrized by the power gamma and the proportionality constant a(0) for models with various asymptotic charges, and the dependence of a(0) on the asymptotic charge is studied. The validity range of the threshold power law is also discussed.

Using ab initio calculations, we fmd that the calculated energy barrier for exchange diffusion of Ni adatoms on Ni(100) surfaces shows a surprisingly large dependence on the size of the surface unit cell. It decreases from 1.39 to 0.78 eV when the cell size changes from (2 x 2) to (6 x 6). This is due to the long-ranged strain field created by the transition state for atomic exchange, which needs a larger cell to relax. The hopping diffusion energy, on the other hand, shows only a very small size effect and remains approximately constant at 0.82-0.86 eV, independently of the cell size. Our results indicate that Ni diffusion on Ni(100) occurs by the exchange mechanism and this is consistent with recent experiments. Previous results obtained using (3 x 3) or (4 x 4) unit cells did not converge sufficiently well to yield correct conclusions.

We propose all the diffraction patterns can be directly transformed to provide three-dimensional atomic structures for the system studied. Depending on the scattering process, either the holography or Patterson transform scheme is used. For diffraction patterns which are generated from a localized emitter source or dominated by an inelastic-scattering feature like core-level photoelectron or low-energy Kikuchi electron, holography transform is needed. On the other hand, for diffraction patterns which were dominated by elastic-scattering, like grazing-incidence X-ray diffraction, electron correlated thermal diffuse scattering or low-energy electron diffraction curves, Patterson transform is needed. To prove our point, high-fidelity and artifact-free three-dimensional atomic structures obtained by transform of low-energy Kikuchi electron patterns and low-energy electron diffraction curves are presented. The future of these direct methods by transforming diffraction patterns will be discussed. (C) 2001 Elsevier Science Ltd. All rights reserved.

We present first-principles calculations that provide a detailed diffusion picture of an adsorbed Si atom on the Si(111)-(7x7) surface. Several diffusion paths for the adsorbed Si atom are established by mapping out the total energy as a function of its positions on the surface. For diffusion between the faulted and unfaulted halves, the energy barriers range from 0.96 to 1.21 eV, while remarkable low-energy barriers from 0.3 to 0.7 eV are discovered within the faulted and unfaulted regions.

More than 50 symmetry-inequivalent emitter-scatterer (E-S) pairs were observed for the Si(111)-(7 x 7) surface using Kikchi electron holography (KEH) with various incident/detection configurations. For different configurations, the E-S pairs in the backscattering direction are preferentially enhanced. Thus, one can obtain detailed structural information by changing the incident/detection direction, which will be helpful in sorting out the correct surface structure model.

Using ab initio density-functional theory, the self-diffusion of adatom and dimer on fcc(100) metal surfaces are studied. For adatom diffusion, we find that the exchange mechanism is favored for Al, Ir, Ni, Pd, Pt and Au, while the hopping mechanism is favored for Rh, Cu, and Ag. Except for Ir/Ir(100), the exchange diffusion energy has a surprising large size-effect and decreases as the surface unit cell increases. This is due to the long-ranged strain-field created at the exchange transition state, which needs a larger cell to relax. The hopping diffusion energy, on the other hand, has a very small size-effect and keeps approximately the same value for various surface unit cells. For self-diffusion on lr(100), the formation of covalent bonds are found at the exchange transition state, and thus the exchange diffusion energy has a little size-effect. Our results also indicate that the exchange mechanism is energetically more favorable for dimer diffusion on fcc(100) surface whenever it is favored for adatom diffusion on fcc(100) surface.

A Patterson-like scheme is proposed for direct inversion of the conventional low-energy electron diffraction (LEED) intensity versus energy I(E) curves, which is in contrast with the previously suggested holographic scheme. Using the Si(111)-(7 X 7) and Si(113)-(3 X 2) surfaces as examples, high quality three-dimensional images, with a resolution better than 0.5 Angstrom, of both surface atoms and bulk atoms are obtained from the direct Patterson inversion of LEED-I(E) curves with the integral-energy phase-summing method.

{

Trimer is the smallest cluster that can have a one-dimensional or two-dimensional structure on surfaces, and it can diffuse and transform between these structures. Using first-principles density-functional theory (DFT) calculations, the structural and dynamical behaviors of Al trimer on Al(111) surface have been studied in detail. Al trimer on Al(111) surface has three different kinds of structure conformations (groups with similar configurations): close-packed (compact) triangular trimers, non-compact triangular trimers, and linear trimers. The close-packed triangular trimers are more stable than the non-compact triangular trimers and the linear trimers, while most of the non-compact triangular trimers are as stable as the linear trimers. For the dynamics of Al trimer on Al(111) surface, there are three different kinds of diffusion mechanisms: (1) concerted translations and rotation of compact triangular trimers (the highest energy barrier by DFT calculation

Using density functional calculations, we demonstrate a catalytic reaction path with activation barriers of less than 0.5 eV for CO oxidation on the neutral and unsupported icosahedral nanoclusters of Au(55), Ag(55), and Au(25)Ag(30). Both CO and O(2) adsorb more strongly on these clusters than on the corresponding bulk surfaces. The reaction path consists of an intermediate involving OOCO complex through which the coadsorption energy of CO and O(2) on these clusters is expected to play an important role in the reaction. Based on the studies for the Au and Ag nanoclusters, a model alloy nanocluster of Au(25)Ag(30) was designed to provide a larger coadsorption energy for CO and O2 and was anticipated to be a better catalyst for CO oxidation from energetic analysis. (C) 2008 American Institute of Physics.

We find that electron diffraction patterns can be directly inverted to provide three-dimensional atomic structures for the system studied. Depending on the scattering process, either holography or a Patterson inversion scheme is used. For diffraction patterns which were generated from a localized emitter source or predominantly by an inelastic-scattering feature like low-energy Kikuchi electrons, holography inversion is needed. The information obtained from Kikuchi electron holography includes the building blocks on the surface and their relative position to the atoms below the surface layer. On the other hand, for diffraction patterns generated predominantly by an inelastic-scattering feature like low-energy electron diffraction (LEED), a Patterson inversion is needed. The information obtained from the Patterson transform of the LEED I(E) curves is the relative positions of surface atoms to the atoms in underlying layers; no intra-layer information can be extracted with this method. High-fidelity and artifact-free three-dimensional atomic structures obtained by inversion of low-energy Kikuchi electron patterns and low-energy electron diffraction curves are presented. The results from the two inversion methods are complementary and can be used to construct or to discriminate the surface atomic structural models. The future of these direct methods by inverting diffraction patterns is discussed.

We used the embedded-atom method potential to study the structures, adsorption energies, binding energies, migration paths, and energy barriers of the Ir adatom and small clusters on fcc Ir (100), (110), and (111) surfaces. We found that the barrier for single-adatom diffusion is lowest on the (111) surface, higher on the (110) surface, and highest on the (100) surface. The exchange mechanisms of adatom diffusion on (100) and (110) surfaces are energetically favored. On all three Ir surfaces, Ir-2 dimers with nearest-neighbor spacing are the most stable. On the (110) surface, the Ir-2 dimer diffuses collectively along the (110) channel, while motion perpendicular to the channel walls is achieved by successive one-atom and correlated jumps. On (111) surface, the Ir-2 dimer diffuses in a zigzag motion on hcp and fee sites without breaking into two single atoms. On the (100) surface, diffusion of the Ir-2 dimer is achieved by successive one-atom exchange with the substrate atom accompanying by a 90 degrees rotation of the Ir-2 dimer. This mechanism has a surprisingly low activation energy of 0.65 eV, which is 0.14 eV lower than the energy for single adatom exchange on the (100) surface. Trimers were found to have a one-dimensional (1D) structure on (100) and (110) surfaces, and a 2D structure on the (111) surface. The observed abrupt drop of the diffusion barrier of tetramer, I-gamma 4 on the Ir (111) surface was confirmed theoretically.

The structure of a type of surface magic cluster is determined by a combination of scanning tunneling microscopy, density-functional calculations, and dynamical low energy electron diffraction. The diffraction method is applicable because these clusters created through hierarchical self-organization of Ga deposited onto a Si(111)-7x7 surface have identical size and structure and form an ordered array with exact translational symmetry. The unprecedented detailed structure information provided by the diffraction measurement is consistent with direct microscopic imaging and theoretical calculations.

Several long-range superstructures have been observed with the scanning tunneling microscopy on the reconstructed Pt(100) surface at room temperature. They are present in strained domains and involve both the Fermi electrons and the concomitant lattice distortions. A first-principles calculation shows that the top layer expanded similar to 18% on average and the Fermi surface for a single hexagon layer displays some nesting portions, which can be related to the wave vectors of the observed superstructures. Thus, these superstructures existing in the local domains of the reconstructed surface have the likely origin of incipent charge density waves.

n/a

Surface-supported planar clusters can sprout active research and create numerous applications in the realm of nanotechnology. Exploitation of these clusters will be more extended if their properties on a supported substrate are thoroughly apprehended, and if they can be fabricated in a controllable way. Here we report finding the magic numbers in two-dimensional Ag clusters grown on Pb quantum islands. We demonstrate, with the images and energy spectra of atomic precision, the transition from electronic origin to a geometric one within the same system. Applying the magic nature, we can also produce a large array of planar clusters with well-defined sizes and shapes.

Experimentally, self-organized Ag planar clusters have been observed on the periodic template found on the Pb quantum islands, which are grown on the Si(111) surface. These planar clusters register a remarkable abundance variation at some specific atomic numbers and possess enhanced stability. They are thus denoted as two-dimensional magic Ag nanoclusters (or nanopucks). In this work, detailed calculations based on ab initio density functional theory are made to illuminate how the size and shape effects related to electronic confinement influence the sequence of these two-dimensional Ag nanostructures. The simulation results demonstrate that the evolution of a sequence of planar magic Ag clusters is strongly correlated with their electronic structures. Meanwhile, the role of substrate in the formation of magic Ag clusters is also examined. The symmetry and size of the periodic pattern on the substrate have helped to build up the distinguishable geometric structures in experiment. Further analysis of the related electronic and geometrical properties of these clusters not only explains the occurrence and sequence of the magic numbers but also helps to elucidate the mechanism of their formation.

In this study, the 13-atom cluster structures of alkaline metals, alkaline-earth metals, boron group metals, carbon group metals, and 3d, 4d, and 5d transition metals in the periodic table are investigated by density functional theory with three kinds of exchange-correlation (XC) functionals: (i) local-density approximation (LDA); (ii) generalized gradient approximation (GGA) with Perdew-Wang 91; and (iii) generalized gradient approximation with Perdew-Burke-Ernzerhof. The dependence on pseudopotentials (PPs) with and without semicore electrons is also examined. The relative energies of five selected high-symmetry three-dimensional and four low-symmetry layer-type isomers for each element of interest are calculated and studied. Among the 44 metallic 13-atom clusters, our results show that the two GGA XC functionals have a great consistency; LDA and GGA results also reveal a great consistency, apart from the Cr, Mn, Fe, Co, Ni, and Rh 13-atom clusters, for which the results show a significant difference. Meanwhile, for most of the elements, the calculations with and without semicore PPs also produce consistent results, except for Cr, Mo, and V, which require a careful treatment of semicore states in the PPs.

The 13-atom metal clusters of fcc elements (Al, Rh, Ir, Ni, Pd, Pt, Cu, Ag, Au) were studied by density functional theory calculations. The global minima were searched for by the ab initio random structure searching method. In addition to some new lowest-energy structures for Pd13 and Au13, we found that the effective coordination numbers of the lowest-energy clusters would increase with the ratio of the dimer-to-bulk bond length. This correlation, together with the electronic structures of the lowest-energy clusters, divides the 13-atom clusters of these fcc elements into two groups (except for Au13, which prefers a two-dimensional structure due to the relativistic effect). Compact-like clusters that are composed exclusively of triangular motifs are preferred for elements without d-electrons (Al) or with (nearly) filled d-band electrons (Ni, Pd, Cu, Ag). Non-compact clusters composed mainly of square motifs connected by some triangular motifs (Rh, Ir, Pt) are favored for elements with unfilled d-band electrons.

n/a

The clean-off reaction of AgO added rows by CO on Ag(110) and Au/Ag(110) bimetallic surfaces was studied by scanning tunneling microscopy (STM) and compared with density functional theory (DFT). This combined study of a model system illustrated the complexity of catalytic enhancement in bimetallic systems. By analyzing in situ time-lapsed STM image series, we found that CO oxidation on a Au-enriched Ag(110) surface leads to an exponential depletion of oxygen with time and a reaction rate that is synergistically enhanced by the presence of Au. First principles calculations indicate that the local atomic configuration around the active reaction sites at the chain ends and the preference of An atom substitution into the subsurface second Ag layer are of critical importance. By calculating CO adsorption energies and reaction barriers for plausible reaction pathways, a detailed description of the CO oxidation reaction emerges, For the optimal reaction pathway, a large (similar to 0.09 eV) barrier reduction and a small barrier of similar to 0.01 eV were found for the Eley-Rideal (ER) mechanism. In contrast, a small (similar to 0.03 eV) barrier reduction and a moderate barrier of similar to 0.23 eV were obtained for the Langmuir-Hinshelwood (LH) mechanism. The ER transitional state was also found to be lower in energy. We conclude that, irrespective of whether the ER mechanism is actually rate dominating, it is definitively enhanced.

n/a

Three-dimensional atomic images of a Pt(111) surface are obtained by direct inversion of multiple low-energy Kikuchi electron-diffraction patterns. The images are in the backscattering direction. and the positions of the images are consistent with those expected from the atomic structure near the Pt(111) surface. The strong electron scattering of the Pt atoms causes no observable problems in the Kikuchi electron holography.

We report calculations of the dissociative adsorption of H-2 at Pd (100) covered with 1/4 monolayer of sulfur using quantum dynamics as well as molecular dynamics and taking all six degrees of freedom of the two H atoms fully into account. The ab-initio potential-energy surface (PES) is found to be very strongly corrugated. In particular, we discuss the influence of tunneling, zero-point vibrations due to the localization of the wave function of the nuclei when narrow valleys of the PES are passed, steering of the approaching H-2 molecules towards low-energy barrier configurations, and the important role of subsurface absorbates for the hydrogen dissociation. It is shown that "established" concepts derived from low-dimensional dynamical studies are not necessarily valid in a high-dimensional treatment. (C) 1998 Elsevier Science B.V. All rights reserved.

n/a

Adsorption of C-60 onto Si(111)root 3 x root 3-In surface presents a fascinating example of interplay between molecular adsorbate and surface structural defects. It has been found that adsorbing C-60 molecules are trapped by the substitutional Si-defects. In turn, the group of a few adsorbed C-60 can act as a trap for the mobile vacancies of the root 3 x root 3-In reconstruction. Namely, adsorbed C-60 induces a strain in the indium layer, and when a mobile vacancy happens to get into the surface area surrounded by fullerenes, the In atoms between the C-60 and the vacancy shift from the T-4 to the H-3 sites, fixing a vacancy in a given location. (C) 2011 Elsevier B.V. All rights reserved.

n/a

Scanning tunneling microscopy (STM) observations of the close-packed C60 fullerene arrays on Si(111)
R3xR3-Ag surface have revealed the presence of dim C60 molecules which constitute 9–12% of all fullerenes. The dim C60 fullerenes reside  1.6 A lower than the bright (“normal”) C60.While the brightC60 are in continuous rotation, the dim C60 are fixed in one of the single orientations, indicating a more tight bonding to the surface. At room temperature (RT), the dynamic switching from bright to dim C60 and vice versa has been detected. Switching slows down with decreasing temperature and becomes completely frozen at 110 K, which implies that the switching is a thermally driven process. RT deposition of  0.1 monolayer of Ag onto C60 array eliminates completely the dim C60 molecules. Experimental results can be understood if one assumes that formation of the dim C60 is associated with disintegration of Ag trimer on Si(111)R3xR3- Ag surface under a given C60 fullerene.

Self-assembly of atoms or molecules on a crystal surface is considered one of the most promising methods to create molecular devices. Here we report a stepwise self-assembly of C60 molecules into islands with unusual shapes and preferred sizes on a gold–indium-covered Si(111) surface. Specifically, 19-mer islands prefer a non-compact boomerang shape, whereas hexagonal 37-mer islands exhibit extraordinarily enhanced stability and abundance. The stepwise self-assembly is mediated by the moiré interference between an island with its underlying lattice, which essentially maps out the adsorption-energy landscape of a C60 on different positions of the surface with a lateral magnification factor and dictates the probability for the subsequent attachment of C60 to an island’s periphery. Our discovery suggests a new method for exploiting the moiré interference to dynamically assist the self-assembly of particles and provides an unexplored tactic of engineering atomic scale moiré magnifiers to facilitate the growth of monodispersed mesoscopic structures.