Island nucleation in thin-film epitaxy: A first-principles investigation

Real-Time KMC Simulation of Vacancy-Mediated Intermixing in Au@Ag Octahedral Core-Cubic Shell Nanocrystals with Ab Initio-Guided Kinetics

Utilization of core-shell rather than monometallic nanocrystals (NCs) facilitates fine-tuning of NC properties for applications. However, compositional evolution via intermixing can degrade these properties prompting recent experimental studies. We develop an atomistic-level stochastic model for vacancy-mediated intermixing exploiting a formalism which allows incorporation at an ab initio density functional theory level of not just the thermodynamics of vacancy formation, but also relevant diffusion barriers for a vast number of possible local environments (in the core and in the shell, at the interface, and in the intermixed phase). This facilitates a predictive treatment and comprehensive understanding of intermixing on the relevant time scale (e.g., 101-103 s). In contrast, previous modeling at the atomistic level utilized only unrealistic generic prescriptions of barriers or employed simplified continuum treatments. For Au@Ag octahedral core-cubic shell NCs, our modeling not only captures the experimentally observed rate or time scale for intermixing of ∼100 s at 450 °C for 60 nm NCs, but also elucidates the underlying rate controlling processes and the effective intermixing barrier.

Versatile stochastic model for predictive KMC simulation of fcc metal nanostructure evolution with realistic kinetics

Stochastic lattice-gas models provide the natural framework for analysis of the surface diffusion-mediated evolution of crystalline metal nanostructures on the appropriate time scale (often 101-104 s) and length scale. Model behavior can be precisely assessed by kinetic Monte Carlo simulation, typically incorporating a rejection-free algorithm to efficiently handle the broad range of Arrhenius rates for hopping of surface atoms. The model should realistically prescribe these rates, or the associated barriers, for a diversity of local surface environments. However, commonly used generic choices for barriers fail, even qualitatively, to simultaneously describe diffusion for different low-index facets, for terrace vs step edge diffusion, etc. We introduce an alternative Unconventional Interaction-Conventional Interaction formalism to prescribe these barriers, which, even with few parameters, can realistically capture most aspects of behavior. The model is illustrated for single-component fcc metal systems, mainly for the case of Ag. It is quite versatile and can be applied to describe both the post-deposition evolution of 2D nanostructures in homoepitaxial thin films (e.g., reshaping and coalescence of 2D islands) and the post-synthesis evolution of 3D nanocrystals (e.g., reshaping of nanocrystals synthesized with various faceted non-equilibrium shapes back to 3D equilibrium Wulff shapes).

Theory of Anisotropic Metal Nanostructures

A significant challenge in the development of functional materials is understanding the growth and transformations of anisotropic colloidal metal nanocrystals. Theory and simulations can aid in the development and understanding of anisotropic nanocrystal syntheses. The focus of this review is on how results from first-principles calculations and classical techniques, such as Monte Carlo and molecular dynamics simulations, have been integrated into multiscale theoretical predictions useful in understanding shape-selective nanocrystal syntheses. Also, examples are discussed in which machine learning has been useful in this field. There are many areas at the frontier in condensed matter theory and simulation that are or could be beneficial in this area and these prospects for future progress are discussed.

Nucleation-mediated reshaping of facetted metallic nanocrystals: Breakdown of the classical free energy picture

Shape stability is key to avoiding degradation of performance for metallic nanocrystals synthesized with facetted non-equilibrium shapes to optimize properties for catalysis, plasmonics, and so on. Reshaping of facetted nanocrystals is controlled by the surface diffusion-mediated nucleation and growth of new outer layers of atoms. Kinetic Monte Carlo (KMC) simulation of a realistic stochastic atomistic-level model is applied to precisely track the reshaping of Pd octahedra and nanocubes. Unexpectedly, separate constrained equilibrium Monte Carlo analysis of the free energy profile during reshaping reveals a fundamental failure of the classical nucleation theory (CNT) prediction for the reshaping barrier and rate. Why? Nucleation barriers can be relatively low for these processes, so the system is not locally equilibrated before crossing the barrier, as assumed in CNT. This claim is supported by an analysis of a first-passage problem for reshaping within a master equation framework for the model that reasonably captures the behavior in KMC simulations.

Direct Measurements of Surface Strain-Mediated Lateral Interactions between Adsorbates in Colloidal Heteroepitaxy

Surface strain can alter the dynamics of adsorbates, and often, the adsorbates themselves induce and interact via their surface strain fields. In epitaxy, such strain-mediated effects get further compounded when a misfit strain exists due to lattice mismatch between the growing film and substrate. Here, we carry out particle-resolved imaging of heteroepitaxial growth of multilayer colloidal films where the particles interact via a short-range attraction. Surprisingly, although the misfit strain relaxed systematically with film thickness, the adcolloid diffusivity was nonmonotonic. We show that this nonmonotonicity stems from the competition between the spatial extent of self-induced in-layer strain and the short interaction range. Importantly, we provide direct evidence for long-ranged strain-mediated interactions between adsorbates and show that it alters the growing film's morphology.

Shape Stability of Truncated Octahedral fcc Metal Nanocrystals

Metallic nanocrystals (NCs) can be synthesized with tailored nonequilibrium shapes to enhance desired properties, e.g., octahedral fcc metal NCs optimize catalytic activity associated with {111} facets. However, maintenance of optimized properties requires stability against thermal reshaping. Thus, we analyze the reshaping of truncated fcc metal octahedra mediated by surface diffusion using a stochastic atomistic-level model with energetic input parameters for Pd. The model describes NC thermodynamics by an effective nearest-neighbor interaction and includes a realistic treatment of diffusive hopping for undercoordinated surface atoms. Kinetic Monte Carlo simulation reveals that the effective barrier, Eeff, for the initial stage of reshaping is strongly tied to the degree of truncation of the vertices in the synthesized initial octahedral shapes. This feature is elucidated via exact analytic determination of the energy variation along the optimal reshaping pathway at low-temperature (T), which involves transfer of atoms from truncated {100} vertex facets to form new layers on {111} side facets. Deviations from predictions of the low-T analysis due to entropic effects are more prominent for higher T and larger NC sizes.

Shape transition and coalescence of Au islands on Ag (110) by molecular dynamics simulation

We have used molecular dynamics simulations based on many body semi-empirical potentials described by the embedded atom method, to analyze and understand the diffusion and coalescence phenomena of Au-clusters during the heteroepitaxial growth on Ag (110) surface. Temperature ranging from 300 to 700 K were considered. In this study, we examined the heterogeneous system Au_n/Ag(110), where n is the number of atoms in each cluster/island (with n = 15, ….35). Our results show that the clusters diffuse on the Ag (110) surface via different diffusion processes, namely, the exchange mechanism and the simple jump, which generate a 2D to 3D transition. Formation and adsorption energies of clusters with different sizes have been computed using static simulations. The dynamic study of coalescence for two islands of system Au₁₅; Au_0 - 9/Ag(110) at different temperatures makes it possible to deduce the detail of cluster shape and the influence of its temperature on the stability of the system and its growth during this evolution.

Vacancy assisted diffusion on single-atom surface alloys

Bimetallic surfaces can exhibit an improved catalytic activity through tailoring the concentration and/or the arrangement of the two metallic components. However, in order to be catalytically active, the active bimetallic surface structure has to be stable under operating conditions. Typically, structural changes in metals occur via vacancy diffusion. Based on the first-principles determination of formation energies and diffusion barriers we have performed kinetic Monte-Carlo (kMC) simulations to analyse the (meta-)stability of PtRu/Ru(0001), AgPd/Pd(111), PtAu/Au(111) and InCu/Cu(100) surface alloys. In a first step, here we consider single-atom alloys together with one vacancy per simulation cell. We will present results of the time evolution of these structures and analyse them in terms of the interaction between the constituents of the bimetallic surface.

Reshaping of Truncated Pd Nanocubes: Energetic and Kinetic Analysis Integrating Transmission Electron Microscopy with Atomistic-Level and Coarse-Grained Modeling

Stability against reshaping of metallic fcc nanocrystals synthesized with tailored far-from-equilibrium shapes is key to maintaining optimal properties for applications such as catalysis. Yet Arrhenius analysis of experimental reshaping kinetics, and appropriate theory and simulation, is lacking. Thus, we use TEM to monitor the reshaping of Pd nanocubes of ∼25 nm side length between 410 °C (over ∼4.5 h) and 440 °C (over ∼0.25 h), extracting a high effective energy barrier of E_eff ≈ 4.6 eV. We also provide an analytic determination of the energy variation along the optimal pathway for reshaping that involves transfer of atoms across the nanocube surface from edges or corners to form new layers on side {100} facets. The effective barrier from this analysis is shown to increase strongly with the degree of truncation of edges and corners in the synthesized nanocube. Theory matches experiment for the appropriate degree of truncation. In addition, we perform simulations of a stochastic atomistic-level model incorporating a realistic description of diffusive hopping for undercoordinated surface atoms, thereby providing a visualization of the initial reshaping process.

Engineering of entanglement and spin state transfer via quantum chains of atomic spins at large separations

Several recent experiments have shown that long-range exchange interactions can determine collective magnetic ground states of nanostructures in bulk and on surfaces. The ability to generate and control entanglement in a system with long-range interaction will be of great importance for future quantum technology. An important step forward to reach this goal is the creation of entangled states for spins of distant magnetic atoms. Herein, the generation of long-distance entanglement between remote spins at large separations in bulk and on surface is studied theoretically, based on a quantum spin Hamiltonian and time-dependent Schrödinger equation for experimentally realized conditions. We demonstrate that long-distance entanglement can be generated between remote spins by using an appropriate quantum spin chain (a quantum mediator), composed by sets of antiferromagnetically coupled spin dimers. Ground state properties and quantum spin dynamics of entangled atoms are studied. We demonstrate that one can increase or suppress entanglement by adding a single spin in the mediator. The obtained result is explained by monogamy property of entanglement distribution inside a quantum spin system. We present a novel approach for non-local sensing of remote magnetic adatoms via spin entanglement.

The diffusion of a Ga atom on GaAs(001)β2(2 × 4): Local superbasin kinetic Monte Carlo

We use first-principles density-functional theory to characterize the binding sites and diffusion mechanisms for a Ga adatom on the GaAs(001)β2(2 × 4) surface. Diffusion in this system is a complex process involving eleven unique binding sites and sixteen different hops between neighboring binding sites. Among the binding sites, we can identify four different superbasins such that the motion between binding sites within a superbasin is much faster than hops exiting the superbasin. To describe diffusion, we use a recently developed local superbasin kinetic Monte Carlo (LSKMC) method, which accelerates a conventional kinetic Monte Carlo (KMC) simulation by describing the superbasins as absorbing Markov chains. We find that LSKMC is up to 4300 times faster than KMC for the conditions probed in this study. We characterize the distribution of exit times from the superbasins and find that these are sometimes, but not always, exponential and we characterize the conditions under which the superbasin exit-time distribution should be exponential. We demonstrate that LSKMC simulations assuming an exponential superbasin exit-time distribution yield the same diffusion coefficients as conventional KMC.

Charge Transport in Molecular Materials: An Assessment of Computational Methods

The booming field of molecular electronics has fostered a surge of computational research on electronic properties of organic molecular solids. In particular, with respect to a microscopic understanding of transport and loss mechanisms, theoretical studies assume an ever-increasing role. Owing to the tremendous diversity of organic molecular materials, a great number of computational methods have been put forward to suit every possible charge transport regime, material, and need for accuracy. With this review article we aim at providing a compendium of the available methods, their theoretical foundations, and their ranges of validity. We illustrate these through applications found in the literature. The focus is on methods available for organic molecular crystals, but mention is made wherever techniques are suitable for use in other related materials such as disordered or polymeric systems.

Adsorbate interactions on the GaN(0001) surface and their effect on diffusion barriers and growth morphology

Studying the adsorbate interactions on a surface helps in understanding the growing surface morphologies and calculating the effective surface diffusion barriers. We study the interaction between Ga-Ga, N-N and Ga-N adatom pairs on the polar GaN(0001) surface using ab initio calculations based on density functional theory. The interaction energy between two adatoms on the surface does not seem to follow definite trends with increasing distance between the adatoms. The presence of a number of possible reconstructions on clean GaN(0001) and periodic effects due to the finite size complicate the analysis of the interactions. Various components of the total interaction energy are separated. We find that there is a large substrate lattice distortion caused due to Ga and N adatoms. The resulting elastic interaction is a major component of the interactions between the adatoms on the GaN(0001) surface. The dipolar interaction is much smaller in magnitude. We also evaluate the component of the interaction energy due to the substrate-mediated electronic interactions. The barriers for surface hopping of adatoms are significantly modified in the presence of other adatoms. We identify several possible surface hopping processes for Ga and N adatoms and calculate their barriers. In particular, we find that the N adatom has a lower barrier to move to an adjoining site on the other side of a neighboring Ga adatom. Kinetic Monte Carlo simulations are performed to see the effect of adatom interactions on the growing surface morphologies of GaN(0001). At the submonolayer growth stage, the fast diffusion of N adatoms located near Ga adatoms leads to more regular island features. In this way, we illustrate the role of adatom interactions in the initial surface nucleation and the morphologies of the growing GaN(0001) film.

Point island models for nucleation and growth of supported nanoclusters during surface deposition

Point island models (PIMs) are presented for the formation of supported nanoclusters (or islands) during deposition on flat crystalline substrates at lower submonolayer coverages. These models treat islands as occupying a single adsorption site, although carrying a label to track their size (i.e., they suppress island structure). However, they are particularly effective in describing the island size and spatial distributions. In fact, these PIMs provide fundamental insight into the key features for homogeneous nucleation and growth processes on surfaces. PIMs are also versatile being readily adapted to treat both diffusion-limited and attachment-limited growth and also a variety of other nucleation processes with modified mechanisms. Their behavior is readily and precisely assessed by kinetic Monte Carlo simulation.

The role of ion transport phenomena in memristive double barrier devices

In this work we report on the role of ion transport for the dynamic behavior of a double barrier quantum mechanical Al/Al₂O₃/Nb_xO_y/Au memristive device based on numerical simulations in conjunction with experimental measurements. The device consists of an ultra-thin Nb_xO_y solid state electrolyte between an Al₂O₃ tunnel barrier and a semiconductor metal interface at an Au electrode. It is shown that the device provides a number of interesting features such as an intrinsic current compliance, a relatively long retention time, and no need for an initialization step. Therefore, it is particularly attractive for applications in highly dense random access memories or neuromorphic mixed signal circuits. However, the underlying physical mechanisms of the resistive switching are still not completely understood yet. To investigate the interplay between the current transport mechanisms and the inner atomistic device structure a lumped element circuit model is consistently coupled with 3D kinetic Monte Carlo model for the ion transport. The simulation results indicate that the drift of charged point defects within the Nb_xO_y is the key factor for the resistive switching behavior. It is shown in detail that the diffusion of oxygen modifies the local electronic interface states resulting in a change of the interface properties.

Imaging the Bonds of Dehalogenated Benzene Radicals on Cu(111) and Au(111)

Dissociative adsorption of doubly substituted benzene molecules leads to formation of benzyne radicals. In this study, co-adsorbed hydrogen molecules are used in scanning tunneling hydrogen microscopy to enhance the contrast of the meta- and the para-isomers of these radicals on Cu(111) and Au(111). Up to three hydrogen molecules are attached to one radical. One hydrogen molecule reveals the orientation of the carbon ring and its adsorption site, allowing discrimination between the two radicals. Two hydrogen molecules reflect the bond picture of the carbon skeleton and reveals that adsorption on Cu(111) distorts the meta- isomer differently from its gas-phase distortion. Three hydrogen molecules allow us to determine the bond picture of a minor species.

Interfacial free energy adjustable phase field crystal model for homogeneous nucleation

To describe the homogeneous nucleation process, an interfacial free energy adjustable phase-field crystal model (IPFC) was proposed by reconstructing the energy functional of the original phase field crystal (PFC) methodology. Compared with the original PFC model, the additional interface term in the IPFC model effectively can adjust the magnitude of the interfacial free energy, but does not affect the equilibrium phase diagram and the interfacial energy anisotropy. The IPFC model overcame the limitation that the interfacial free energy of the original PFC model is much less than the theoretical results. Using the IPFC model, we investigated some basic issues in homogeneous nucleation. From the viewpoint of simulation, we proceeded with an in situ observation of the process of cluster fluctuation and obtained quite similar snapshots to colloidal crystallization experiments. We also counted the size distribution of crystal-like clusters and the nucleation rate. Our simulations show that the size distribution is independent of the evolution time, and the nucleation rate remains constant after a period of relaxation, which are consistent with experimental observations. The linear relation between logarithmic nucleation rate and reciprocal driving force also conforms to the steady state nucleation theory.

Ab initio Kinetic Monte Carlo simulations of dissolution at the NaCl-water interface

We have used ab initio molecular dynamics (AIMD) simulations to study the interaction of water with the NaCl surface. As expected, we find that water forms several ordered hydration layers, with the first hydration layer having water molecules aligned so that oxygen atoms are on average situated above Na sites. In an attempt to understand the dissolution of NaCl in water, we have then combined AIMD with constrained barrier searches, to calculate the dissolution energetics of Na(+) and Cl(-) ions from terraces, steps, corners and kinks of the (100) surface. We find that the barrier heights show a systematic reduction from the most stable flat terrace sites, through steps to the smallest barriers for corner and kink sites. Generally, the barriers for removal of Na(+) ions are slightly lower than for Cl(-) ions. Finally, we use our calculated barriers in a Kinetic Monte Carlo as a first order model of the dissolution process.

Kinetic Monte Carlo simulation of CO adsorption on sulfur-covered Pd(100)

The use of atomistic Kinetic Monte Carlo method was explored to examine the influence of sulfur poisoning on CO adsorption on Pd(100) surface. The model explicitly incorporates key elementary processes such as CO adsorption and CO desorption including diffusion of surface CO and S species. Relevant energetic and kinetic parameters were derived using information calculated from density functional theory as a starting point. Kinetic Monte Carlo simulation was performed to determine relevant observables such as CO saturation coverage as a function of amount of preadsorbed sulfur and to predict temperature programmed desorption spectra.

Real-time ab initio KMC simulation of the self-assembly and sintering of bimetallic epitaxial nanoclusters: Au + Ag on Ag(100)

Far-from-equilibrium shape and structure evolution during formation and post-assembly sintering of bimetallic nanoclusters is extremely sensitive to the periphery diffusion and intermixing kinetics. Precise characterization of the many distinct local-environment-dependent diffusion barriers is achieved for epitaxial nanoclusters using density functional theory to assess interaction energies both with atoms at adsorption sites and at transition states. Kinetic Monte Carlo simulation incorporating these barriers then captures structure evolution on the appropriate time scale for two-dimensional core-ring and intermixed Au-Ag nanoclusters on Ag(100).

Enhanced atom mobility on the surface of a metastable film

A remarkable enhancement of atomic diffusion is highlighted by scanning tunneling microscopy performed on ultrathin metastable body-centered tetragonal Co films grown on Fe(001). The films follow a nearly perfect layer-by-layer growth mode with a saturation island density strongly dependent on the layer on which the nucleation occurs, indicating a lowering of the diffusion barrier. Density functional theory calculations reveal that this phenomenon is driven by the increasing capability of the film to accommodate large deformations as the thickness approaches the limit at which a structural transition occurs. These results disclose the possibility of tuning surface diffusion dynamics and controlling cluster nucleation and self-organization.

Thermal string excitations in artificial spin-ice square dipolar array

We report on a theoretical investigation of artificial spin-ice dipolar arrays using a nanoisland shape adopted from recent experiments (Farhan et al 2013 Nature Phys. 9 375). The number of thermal magnetic string excitations in the square lattice is drastically increased by a vertical displacement of rows and columns. We find large increments especially for low temperatures and for string excitations with quasi-monopoles of charges ± 4. By kinetic Monte Carlo simulations we address the thermal stability of such excitations, thereby providing time scales for their experimental observation.

Acetylene on Cu(111): imaging a molecular surface arrangement with a constantly rearranging tip

Acetylene on Cu(111) is investigated by scanning tunnelling microscopy (STM); a surface pattern previously derived from diffraction measurements can be validated, if the variation of the STM image transfer function through absorption of an acetylene molecule onto the tip apex is taken into account. Density functional theory simulations point to a balance between short-range repulsive interactions of acetylene/Cu(111) associated with surface stress and longer range attractive interactions as the origin of the ordering.

High island densities and long range repulsive interactions: Fe on epitaxial graphene

The understanding of metal nucleation on graphene is essential for promising future applications, especially of magnetic metals which can be used in spintronics or computer storage media. A common method to study the grown morphology is to measure the nucleated island density n as a function of growth parameters. Surprisingly, the growth of Fe on graphene is found to be unusual because it does not follow classical nucleation: n is unexpectedtly high, it increases continuously with the deposited amount θ and shows no temperature dependence. These unusual results indicate the presence of long range repulsive interactions. Kinetic Monte Carlo simulations and density functional theory calculations support this conclusion. In addition to answering an outstanding question in epitaxial growth, i.e., to find systems where long range interactions are present, the high density of magnetic islands, tunable with θ, is of interest for nanomagnetism applications.

Surface nucleation and growth in the system of interacting particles

Recent experiments on epitaxial growth of metals on graphene have shown a strong dependence of island densities on coverage. These investigations cannot be explained by the standard mean-field nucleation theories. To understand them, we extend to higher coverage the former theory of rate equations developed for the initial state of nucleation, in a system where adsorbate interaction is included. We account for that, in the case of high coverage, the repulsive interaction influences both the attachment of monomers to clusters and the mobility of atoms. In our work we analyze the modification of the dependence of the island density on coverage, temperature and F/D ratio. In some regimes our theory results in the experimentally observed substantial growth of island density with coverage for a high deposited amount and a weak dependence on deposition rate F. We also find out the local maxima in temperature dependence of island density, as a consequence of long-range repulsive interactions.

Do two-dimensional "noble gas atoms" produce molecular honeycombs at a metal surface?

Anthraquinone self-assembles on Cu(111) into a giant honeycomb network with exactly three molecules on each side. Here we propose that the exceptional degree of order achieved in this system can be explained as a consequence of the confinement of substrate electrons in the pores, with the pore size tailored so that the confined electrons can adopt a noble-gas-like two-dimensional quasi-atom configuration with two filled shells. Formation of identical pores in a related adsorption system (at different overall periodicity due to the different molecule size) corroborates this concept. A combination of photoemission spectroscopy with density functional theory computations (including van der Waals interactions) of adsorbate-substrate interactions allows quantum mechanical modeling of the spectra of the resultant quasi-atoms and their energetics.

Adsorbates in a box: titration of substrate electronic states

Nanoscale confinement of adsorbed CO molecules in an anthraquinone network on Cu(111) with a pore size of ≈4  nm arranges the CO molecules in a shell structure that coincides with the distribution of substrate confined electronic states. Molecules occupy the states approximately in the sequence of rising electron energy. Despite the sixfold symmetry of the pore boundary itself, the adsorbate distribution adopts the threefold symmetry of the network-substrate system, highlighting the importance of the substrate even for such quasi-free-electron systems.

Melting of two-dimensional adatom superlattices stabilized by long-range electronic interactions

The melting transition of Ce adatom superlattices stabilized by long-range substrate-mediated electronic interactions on Cu(111) and Ag(111) noble metal surfaces has been investigated by low-temperature scanning tunneling microscopy, density functional theory calculations, and kinetic Monte Carlo simulations. Intriguingly, owing to the interaction between Ce adatoms and substrate, these superlattices undergo two-dimensional melting to a liquid without transition through the hexatic phase. The crucial parameters for this direct solid to liquid transition are identified.

Formation of complex wedding-cake morphologies during homoepitaxial film growth of Ag on Ag(111): atomistic, step-dynamics, and continuum modeling

An atomistic lattice-gas model is developed which successfully describes all key features of the complex mounded morphologies which develop during deposition of Ag films on Ag(111) surfaces. We focus on this homoepitaxial thin film growth process below 200 K. The unstable multilayer growth mode derives from the presence of a large Ehrlich-Schwoebel step-edge barrier, for which we characterize both the step-orientation dependence and the magnitude. Step-dynamics modeling is applied to further characterize and elucidate the evolution of the vertical profiles of these wedding-cake-like mounds. Suitable coarse-graining of these step-dynamics equations leads to instructive continuum formulations for mound evolution.

Direct evidence for the effect of quantum confinement of surface-state electrons on atomic diffusion

We report on the direct observations of the effect of quantum confinement of surface-state electrons on atomic diffusion. Confined electronic states induced by open nanoscale resonators [consisting of two parallel monatomic Cu chains on Cu(111)] are studied by means of scanning tunneling microscope measurements and first-principles calculations. Strongly anisotropic diffusion of adatoms around and inside resonators is revealed at low temperatures. The formation of diffusion channels and empty zones is demonstrated. We show that it is possible to engineer atomic diffusion by varying the distance between the resonator walls.

Monodisperse microisland formation on Ni/Ru(0001) monolayers

The creation of identical microislands consisting of Ni trimers and multiples thereof on Ru(0001) induced by oxygen adsorption has been observed using scanning tunnelling microscopy. The island formation is caused by an oxygen induced expulsion of Ni atoms or trimers out of the moiré-distorted (densified) Ni monolayer. The exceptional stability of the Ni trimers is attributed to oxygen attachment, forming Ni-oxygen composites, as verified by detailed density functional theory calculations. The high density, identical structure, and notable thermal stability of these islands open up new perspectives for the study of the properties of nanostructured surfaces.

Adsorption-induced desorption of benzene on Si(111)-7 x 7 by substrate-mediated electronic interactions

The process of benzene adsorption on an adjacent adatom-rest atom pair on Si(111)-7 x 7 at room temperature was studied using in-situ scanning tunneling microscopy (STM). Both adsorption and desorption of benzene were observed to take place mostly at adjacent sites during the process. DFT calculation results show that the bond length between the rest atom and the carbon atom in a pre-adsorbed benzene molecule increases due to the charge transfer from a neighboring rest atom in response to an approaching benzene molecule. Such increase in the bond length, when coupled resonantly to the C-Si thermal vibration, could result in bond breakage and desorption of the adsorbate. The studies provide evidence for the desorption of a chemisorbed benzene caused by an adsorbing benzene at a neighboring site through a substrate-mediated electronic interaction.

Restructuring of cobalt nanoparticles induced by formation and diffusion of monodisperse metal-sulfur complexes

Time-resolved scanning tunneling microscopy (STM) is used to investigate a massive sulfur-induced transformation of a homogeneous array of approximately 2 nm Co nanoparticles into a new cobalt sulfide phase. The underlying atomistic mass-transport process is revealed and, surprisingly, found to be mediated exclusively by the formation and detachment of monosized Co3S4 complexes at the perimeter of the Co nanoparticles. The process is followed by fast diffusion, agglomeration of the complexes, and subsequent crystallization into a cobalt sulfide phase.

Monomer adsorption on equilateral triangular lattices with attractive first-neighbor interactions

We have recently studied a model of monomer adsorption on infinitely long equilateral triangular lattices with terraces of finite width M and nonperiodic boundaries. This study was restricted to the case of repulsive adsorbate-adsorbate first-neighbor interactions but included attractive, repulsive, and negligible second-neighbor interactions. The present work extends this study to the case of attractive first-neighbors, and the phases are determined, as before, with a confidence exceeding 10 significant figures. Phase diagrams are included for terrace widths M < or =11. Most of the occupational characteristics of the phases fit exact analytic expressions in M. The infinite-M limit of these expressions, combined with other analyses, provide the complete phase diagram for the infinite two-dimensional lattice. In addition to the empty and full coverage phases, there are three phases exhibiting stripe and cluster features that were not observed in the case of repulsive first-neighbors.

Coverage and nearest-neighbor dependence of adsorbate diffusion

We present data on the coverage and nearest-neighbor dependences of the diffusion of CO on Cu(111) by time-lapsed scanning tunneling microscope (STM) imaging. Most notable is a maximum in diffusivity of CO at a local coverage of one molecule per 20 substrate atoms and a repulsion between CO molecules upon approach closer than three adsites, which in combination with a less pronounced increase in potential energy at the diffusion transition state, leads to rapid diffusion of CO molecules around one another. We propose a new method of evaluating STM-based diffusion data that provides all parameters necessary for the modeling of the dynamics of an adsorbate population.

Long-range electronic interactions at a high temperature: bromine adatom islands on Cu(111)

Long-range electronic interactions between Br adatom islands, which are formed at approximately 600 K, on Cu(111) are mediated by substrate surface-state electrons at that elevated temperature. Using scanning tunneling microscopy at 4 K, we have quantified nearest neighbor island separations and found favored spacings to be half-multiples of the Fermi wavelength of Cu(111). The strong interaction potential and decay length of the interisland interactions are discussed in terms of the interaction of Br with the substrate surface state.

Low-temperature ultrafast mobility in systems with long-range repulsive interactions: Pb/Si(111)

A realization of the numerous phases predicted in systems with long-range repulsive interactions was recently found in Pb/Si(111). Surprisingly, these numerous phases can be grown at low temperatures approximately 40 K over macroscopic distances. This unusual observation can be explained from theoretical calculations of the collective diffusion coefficient D(c) in systems with long-range repulsive interactions. Instead of a gradual dependence of D(c) on coverage, it was found that D(c) has sharp maxima at low temperatures for every stable phase (i.e., for every rational value of the coverage theta=p/q) in agreement with the experiment.

Synchronous relaxation algorithm for parallel kinetic Monte Carlo simulations of thin film growth

We present an optimistic synchronous relaxation algorithm for parallel kinetic Monte Carlo (KMC) simulations of thin film growth. This algorithm is based on spatial decomposition of the KMC lattice and it employs two measures aimed at improving the parallel efficiency: dynamic global updating and domain boundary shifting. We utilize this algorithm to simulate two different growth models, which represent the growth of Ag on Ag(111) and the heteroepitaxial growth of Ag on Pt(111). We show that these simulations can achieve good efficiency-especially for large domain sizes with a moderate number of processors. We find that domain boundary shifting can improve efficiency-especially for simulations of growth in the AgPt(111) system, where the potential-energy surface topology creates areas of rapid, localized motion. We analyze the origins of parallel efficiency in these simulations.

Adatom self-organization induced by quantum confinement of surface electrons

We present a novel mechanism of nanostructure growth based on quantum confinement of surface-state electrons. Ab initio calculations and the kinetic Monte Carlo simulations reveal the phenomenon of confinement-induced adatom self-organization in quantum corrals. Our studies indicate that new atomic-scale nanostructures can be engineered exploiting the quantum confinement of surface electrons.

Trapping-mediated chemisorption of ethylene on Si(001)-c(4 x 2)

Adsorption of ethylene molecules on Si(001)-c(4 x 2) was studied using scanning tunneling microscopy at low temperatures. Ethylene molecules trapped at the surface at 50 K were imaged only after decay to chemisorption, each bonding to a Si dimer. Atomic-scale observations of temperature-dependent kinetics show that the decay exhibited Arrhenius behavior with the reaction barrier of 128 meV in clear evidence of the trapping-mediated chemisorption, however, with an anomalously small preexponential factor of 300 Hz. Such a small prefactor is attributed to the entropic bottleneck at the transition state caused by the free-molecule-like trap state.

Critical epinucleation on reconstructured surfaces and first-principle calculation of homonucleation on Si(100)

We introduce the concept of "critical epinucleation to distinguish nucleation on surfaces with and without reconstruction. On a reconstructed surface, the critical classical nucleus is stable against dissociation, but may not yet break the underlying surface reconstruction. Consequently, there must exist a "critical epinucleus" that is not only stable but also has established the epiconfiguration by unreconstructing the underlying substrate. We illustrate this concept by first-principle calculation of homonucleation on reconstructed Si(001) surface where the critical epinucleus consists of six adatoms.

Self-organized growth of nanopucks on Pb quantum islands

Electronic Moirè patterns found on lead (Pb) quantum islands can serve as a template to grow self-organized cluster (nanopucks) arrays of various materials. These patterns can be divided into fcc- and hcp-stacked areas, which exhibit different binding strengths to the deposited adatoms. For Ag adatoms, the binding energy can differ substantially and the confined nucleation thus occurs in the fcc sites. Both the size distribution and spatial arrangement of the Ag nanopucks are analyzed and found to be commensurate with the characteristics of the template island, which exhibits a bilayer oscillatory behavior.