Search

Using scanning tunneling microscopy, we observe an adlayer structure that is dominated by short rows of S atoms, on unreconstructed regions of a Au(111) surface. This structure forms upon adsorption of low S coverage (less than 0.1 monolayer) on a fully reconstructed clean surface at 300 K, then cooling to 5 K for observation. The rows adopt one of three orientations that are rotated by 30 from the close-packed directions of the Au(111) substrate, and adjacent S atoms in the rows are separated by √ 3 times the surface lattice constant, a. Monte Carlo simulations are performed on lattice-gas models, derived using a limited cluster expansion based on density functional theory energetics. Models which include long-range pairwise interactions (extending to 5a), plus selected trio interactions, successfully reproduce the linear rows of S atoms at reasonable temperatures. [ABSTRACT FROM AUTHOR]

A rich menagerie of structures is identified at 5 K following adsorption of low coverages (≤0.05 monolayers) of S on Cu(111) at room temperature. This paper emphasizes the reconstructions at the steps. The A-type close-packed step has 1 row of S atoms along its lower edge, where S atoms occupy alternating pseudo-fourfold-hollow (p4fh) sites. Additionally, there are 2 rows of S atoms of equal density on the upper edge, bridging a row of extra Cu atoms, together creating an extended chain. The B-type close-packed step exhibits an even more complex reconstruction, in which triangle-shaped groups of Cu atoms shift out of their original sites and form a base for S adsorption at (mostly) 4fh sites. We propose a mechanism by which these triangles could generate Cu-S complexes and short chains like those observed on the terraces. [ABSTRACT FROM AUTHOR]

Threshold versions of Schloegl's model on a lattice, which involve autocatalytic creation and spontaneous annihilation of particles, can provide a simple prototype for discontinuous non-equilibrium phase transitions. These models are equivalent to so-called threshold contact processes. A discontinuous transition between populated and vacuum states can occur selecting a threshold of N ≥ 2 for the minimum number, N, of neighboring particles enabling autocatalytic creation at an empty site. Fundamental open questions remain given the lack of a thermodynamic framework for analysis. For a square lattice with N = 2, we show that phase coexistence occurs not at a unique value but for a finite range of particle annihilation rate (the natural control parameter). This generic two-phase coexistence also persists when perturbing the model to allow spontaneous particle creation. Such behavior contrasts both the Gibbs phase rule for thermodynamic systems and also previous analysis for this model. We find metastability near the transition corresponding to a non-zero effective line tension, also contrasting previously suggested critical behavior. Mean-field type analysis, extended to treat spatially heterogeneous states, further elucidates model behavior. [ABSTRACT FROM AUTHOR]

We explore simple lattice-gas reaction models for CO-oxidation on 1D and 2D periodic arrays of surface adsorption sites with CO adsorption and desorption, dissociative O2 adsorption and recombinative desorption (at low rate), and CO + O reaction to form CO2. Adspecies interactions are neglected, and adspecies diffusion is effectively absent. The models are motivated by studies of CO-oxidation on RuO2(110) at high-pressures. Despite the lack of adspecies interactions, negligible adspecies diffusion results in kinetically induced spatial correlations. A transition occurs from a random primarily CO-populated steady-state at high CO-partial pressure, pCO, to a strongly correlated near-O-covered steady-state for low pCO as noted by Matera et al. [J. Chem. Phys. 134, 064713 (2011)]. In addition, we identify a second transition to a random near-O-covered steady-state at very low pCO. Furthermore, we identify and analyze the slow "diffusive dynamics" for very low pCO and provide a detailed characterization of the crossover to the strongly correlated O-covered steady-state as well as of the spatial correlations in that state. [ABSTRACT FROM AUTHOR]

We develop statistical mechanical models amenable to analytic treatment for the dissociative adsorption of O2 at hollow sites on fcc(100) metal surfaces. The models incorporate exclusion of nearest-neighbor pairs of adsorbed O. However, corresponding simple site-blocking models, where adsorption requires a large ensemble of available sites, exhibit an anomalously fast initial decrease in sticking. Thus, in addition to blocking, our models also incorporate more facile adsorption via orientational steering and funneling dynamics (features supported by ab initio Molecular Dynamics studies). Behavior for equilibrated adlayers is distinct from those with finite adspecies mobility. We focus on the low-temperature limited-mobility regime where analysis of the associated master equations readily produces exact results for both short- and long-time behavior. Kinetic Monte Carlo simulation is also utilized to provide a more complete picture of behavior. These models capture both the initial decrease and the saturation of the experimentally observed sticking versus coverage, as well as features of non-equilibrium adlayer ordering as assessed by surface-sensitive diffraction. [ABSTRACT FROM AUTHOR]

Scanning tunneling microscopy studies reveal that trace amounts of adsorbed S below a critical coverage on the order of 10 mML have little effect on the coarsening and decay of monolayer Ag adatom islands on Ag(111) at 300 K. In contrast, above this critical coverage, decay is greatly accelerated. This critical value appears to be determined by whether all S can be accommodated at step edges. Accelerated coarsening derives from the feature that the excess S (above that incorporated at steps) produces significant populations on the terraces of metal-sulfur complexes, which are stabilized by strong Ag–S bonding. These include AgS2, Ag2S2, Ag2S3, and Ag3S3. Such complexes are sufficiently populous and mobile that they can potentially lead to greatly enhanced metal mass transport across the surface. This picture is supported by density functional theory analysis of the relevant energetics, as well as by reaction-diffusion equation modeling to assess the mechanism and degree of enhanced coarsening. [ABSTRACT FROM AUTHOR]

We analyze a discontinuous nonequilibrium phase transition between an active (or reactive) state and a poisoned (or extinguished) state occurring in a stochastic lattice-gas realization of Schloegl’s second model for autocatalysis. This realization, also known as the quadratic contact process, involves spontaneous annihilation, autocatalytic creation, and diffusion of particles on a square lattice, where creation at empty sites requires a suitable nearby pair of particles. The poisoned state exists for all annihilation rates p>0 and is an absorbing particle-free “vacuum” state. The populated active steady state exists only for p below a critical value, pe. If pf denotes the critical value below which a finite population can survive, then we show that pf

Within the framework of a realistic atomistic lattice-gas model, we present the theoretical formulation and simulation procedures for precise analysis of the chemical diffusion flux of highly mobile CO within a nonuniform interacting mixed CO+O adlayer on a Pd(100) surface. The approach applies in both regimes of relatively immobile unequilibrated and fairly mobile near-equilibrated O adlayer distributions. Spatiotemporal behavior in surface reactions is controlled by chemical diffusion in mixed adlayers. Thus, we naturally integrate the above analysis with a previously developed multiscale modeling strategy to describe mesoscale reaction front propagation in CO oxidation on Pd(100). This treatment avoids using a simplified prescription of chemical diffusion and reaction kinetics as in traditional mean-field reaction-diffusion equation approaches. [ABSTRACT FROM AUTHOR]

We have developed an atomistic lattice-gas model for the catalytic oxidation of CO on single-crystal Pd(100) surfaces under ultrahigh vacuum conditions. This model necessarily incorporates an detailed description of adlayer ordering and adsorption-desorption kinetics both for CO on Pd(100), and for oxygen on Pd(100). Relevant energetic parameters are determined by comparing model predictions with experiment, together with some guidance from density functional theory calculations. The latter also facilitates description of the interaction and reaction of adsorbed CO and oxygen. Kinetic Monte Carlo simulations of this reaction model are performed to predict temperature-programed reaction spectra, as well as steady-state bifurcation behavior. [ABSTRACT FROM AUTHOR]

Using a lattice-gas model with pairwise interactions, we study the ordered structures, coverage dependence of the heat of adsorption, and other experimentally observable behavior of adsorbed CO overlayers on Pd(100) single crystal surfaces. Transfer matrix and Monte Carlo methods give accurate information regarding the lattice-gas model that often contradicts simple mean-field-like analysis. We demonstrate the usefulness of the model by reproducing experimental results over a large range of pressures and temperatures. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Using scanning tunneling microscopy, we characterize the size and bias-dependent shape of sulfur atoms on Cu(100) at low coverage (below 0.1 monolayers) and low temperature (quenched from 300 to 5 K). Sulfur atoms populate the Cu(100) terraces more heavily than steps at low coverage, but as coverage approaches 0.1 monolayers, close-packed step edges become fully populated, with sulfur atoms occupying sites on top of the step. Density functional theory (DFT) corroborates the preferential population of terraces at low coverage as well as the step adsorption site. In experiment, small regions with p(2 x 2)-like atomic arrangements emerge on the terraces as sulfur coverage approaches 0.1 monolayer. Using DFT, a lattice gas model has been developed, and Monte Carlo simulations based on this model have been compared with the observed terrace configurations. A model containing eight pairwise interaction energies, all repulsive, gives qualitative agreement. Experiment shows that atomic adsorbed sulfur is the only species on Cu(100) up to a coverage of 0.09 monolayers. There are no Cu-S complexes. In contrast, prior work has shown that a Cu2S3 complex forms on Cu(111) under comparable conditions. On the basis of DFT, this difference can be attributed mainly to stronger adsorption of sulfur on Cu(100) as compared with Cu(111). [ABSTRACT FROM AUTHOR]

Tailored multisite lattice-gas (msLG) models are developed for CO oxidation on Pd(100) at low-pressures. These models include multiple adsorption site types and superlattice adlayer ordering due to short-range exclusion for highly mobile reactant adspecies. However, they are simplified to neglect longer-range weaker adspecies interactions, so that the key energetic parameters are the CO desorption barrier and the reaction barrier. We discuss existing density functional theory results for these energies and present additional analysis for CO adsorption. After also including an appropriate nontrivial specification of the dynamics of adsorption onto mixed reactant adlayers, we develop rate equations for the reaction kinetics. Our formulation goes beyond traditional mean-field (MF) Langmuirian treatments by accounting for multiple adsorption sites and for the strong spatial correlations associated with superlattice ordering. Specifically, we utilize factorization approximations based on appropriate site motifs, and also Padé resummation of exact low-coverage expansions for sticking coefficients. Our beyond-MF rate equations are successful in accurately predicting key aspects of reactive steady-state behavior, and thus expand the utility of rate equation formulations in surface chemistry. This is confirmed by comparison with precise kinetic Monte Carlo simulation results. Specifically, we not only assess bistability and criticality observed for CO oxidation but also find more complex multistability associated with symmetry-breaking transitions in high-coverage CO adlayers. [ABSTRACT FROM AUTHOR]

The article discusses the stochastic modeling of catalytic reaction-diffusion processes in nonporous materials and 2D surface. Topics discussed include analysis of basic irreversible first-order conversion reaction in a 1D linear nanopore, multi-site lattice gas modeling of oxidation reactions on metal catalyst surfaces at low pressure and analysis of catalytic reactions on metal oxide surfaces.

Sulfur-metal complexes, containing only a few atoms, can open new, highly efficient pathways for transport of metal atoms on surfaces. For example, they can accelerate changes in the shape and size of morphological features, such as two-dimensional nanoclusters, over time. In this study we perform STM under conditions that are designed to specifically isolate such complexes. We find a new, unexpected S-Cu complex on the Cu(111) surface, which we identify as Cu2S3. We propose that Cu2S3 enhances mass transport in this system, which contradicts a previous proposal based on Cu3S3. We analyze bonding within these Cu-S complexes, identifying a principle for stabilization of sulfur complexes on coinage metal surfaces. [ABSTRACT FROM AUTHOR]

An accurate description of oxygen dissociation pathways and kinetics for various local adlayer environments is key for an understanding not just of the coverage dependence of oxygen sticking, but also of reactive steady states in oxidation reactions. Density functional theory analysis for M (100) surfaces with M = Pd, Rh, and Ni, where O prefers the fourfold hollow adsorption site, does not support the traditional Brundle-Behm-Barker picture of dissociative adsorption onto second-nearest-neighbor hollow sites with an additional blocking constraint. Rather adsorption via neighboring vicinal bridge sites dominates, although other pathways can be active. The same conclusion also applies for M = Pt and Ir, where oxygen prefers the bridge adsorption site. Statistical mechanical analysis is performed based on kinetic Monte Carlo simulation of a multisite lattice-gas model consistent with our revised picture of adsorption. This analysis determines the coverage and temperature dependence of sticking for a realistic treatment of the oxygen adlayer structure. [ABSTRACT FROM AUTHOR]

Schloegl's second model on a (d ≥ 2)-dimensional hypercubic lattice involves: (i) spontaneous annihilation of particles with rate p and (ii) autocatalytic creation of particles at vacant sites at a rate proportional to the number of suitable pairs of neighboring particles. This model provides a prototype for nonequilibrium discontinuous phase transitions. However, it also exhibits nontrivial generic two-phase coexistence: Stable populated and vacuum states coexist for a finite range, pf(d) < p < pe(d), spanned by the orientation-dependent stationary points for planar interfaces separating these states. Analysis of interface dynamics from kinetic Monte Carlo simulation and from discrete reaction-diffusion equations (dRDEs) obtained from truncation of the exact master equation, reveals that pe(f) ∼ 0.211 3765 + ce(f),/d as d →∞, where Δc = ce -- Cf &thkap; 0.014. A metastable populated state persists above pe(d) up to a spinodal p = ps(d), which has a well-defined limit ps(d →∞) = ¼. The dRDEs display artificial propagation failure, absent in the stochastic model due to fluctuations. This feature is amplified for increasing d, thus complicating our analysis. [ABSTRACT FROM AUTHOR]

A two-dimensional atomistic realization of Schlögl’s second model for autocatalysis is implemented and studied on a square lattice as a prototypical nonequilibrium model with first-order transition. The model has no explicit symmetry and its phase transition can be viewed as the nonequilibrium counterpart of liquid-vapor phase separations. We show some familiar concepts from study of equilibrium systems need to be modified. Most importantly, phase coexistence can be a generic feature of the model, occurring over a finite region of the parameter space. The first-order transition becomes continuous as a temperature-like variable increases. The associated critical behavior is studied through Monte Carlo simulations and shown to be in the two-dimensional Ising universality class. However, some common expectations regarding finite-size corrections and fractal properties of geometric clusters for equilibrium systems seems to be inapplicable. [ABSTRACT FROM AUTHOR]

The dynamics of nanoscale clusters can be distinct from macroscale behavior described by continuum formalisms. For diffusion of 2D clusters of N atoms in homoepitaxial systems mediated by edge atom hopping, macroscale theory predicts simple monotonic size scaling of the diffusion coefficient, DN - Nβ , with β = 3/2. However, modeling for nanoclusters on metal(100) surfaces reveals that slow nucleation-mediated diffusion displaying weak size scaling β < 1 occurs for "perfect" sizes Np = L2 and L(L+1) for integer L = 3,4,. . . (with unique square or near-square ground state shapes) and also for Np+3, Np+4,. . . . In contrast, fast facile nucleation-free diffusion displaying strong size scaling β≈2.5 occurs for sizesNp+1 andNp+2.DN versusNoscillates strongly between the slowest branch (for Np+3) and the fastest branch (for Np+1). All branches merge for N = O(102), but macroscale behavior is only achieved for much larger N = O(103). This analysis reveals the unprecedented diversity of behavior on the nanoscale. [ABSTRACT FROM AUTHOR]

Lattice-gas models for CO oxidation can exhibit a discontinuous nonequilibrium transition between reactive and inactive states, which disappears above a critical CO-desorption rate. Using finite-size-scaling analysis, we demonstrate a crossover from Ising to mean-field behavior at the critical point, with increasing surface mobility of adsorbed CO or with decreasing system size. This behavior is elucidated by analogy with that of equilibrium Ising-type systems with long-range interactions. [ABSTRACT FROM AUTHOR]

Studies the reaction mechanisms for oxidation of the silicon (100) surface by atomic oxygen with high-level quantum mechanical methods. Use of quantum mechanics and molecular mechanics method; Initial surface products for oxygen adsorption; Ground state of the separated reactions, oxygen and silicon; Two channels for silicon oxide desorption.

To elucidate surface diffusion in the presence of a coadsorbate with superlattice ordering, we consider particle hopping on a square lattice with some fraction, θ[sub B], of quenched blocking sites arranged with checkerboard or c(2x2) ordering. Behavior for low θ[sub B] corresponds to diffusion around isolated obstacles, and can be described by exact density expansions. Behavior for high θ[sub B] corresponds to percolative diffusion along (or sometimes away from) domain boundaries. The connectivity of these domain boundaries is closely related to the existence of symmetry breaking [i.e., long-range c(2x2) order] in the distribution of blocking sites. In some cases, symmetry breaking induces critical behavior for diffusive transport which is fundamentally different from that for the conventional "ant in the labyrinth" problem. Our results apply to recently developed models for CO oxidation, where CO(ads) diffuses rapidly through coadsorbed relatively immobile c(2x2)-O(ads). The characterization of CO diffusion in these systems is key to describing spatial pattern formation. © 2000 American Institute of Physics. [ABSTRACT FROM AUTHOR]

We consider Schloegl models (or contact processes) where particles on a square grid annihilate at a rate p and are created at a rate of kn=n(n-1)/[N(N-1)] at empty sites with n particles in a neighborhood ΩN of size N. Simulation reveals a discontinuous transition between populated and vacuum states, but equistable p=peq determined by the stationarity of planar interfaces between these states depends on the interface orientation and on ΩN. The behavior for large ΩN follows from continuum equations. These also depend on the interface orientation and on ΩN shape, but a unique peq=0.211 376 320 4 emerges imposing a Gibbs phase rule. [ABSTRACT FROM AUTHOR]

For diffusion of two-dimensional homoepitaxial clusters of N atoms on metal (100) surfaces mediated by edge atom hopping, macroscale continuum theory suggests that the diffusion coefficient scales like DN ~N -β with β = 3/2. However, we find quite different and diverse behavior in multiple size regimes. These include: (i) facile diffusion for small sizes N < 9 ; (ii) slow nucleation-mediated diffusion with small β < 1 for "perfect" sizes N = Np = L2 or L(L+1), for L = 3, 4, ... having unique ground-state shapes, for moderate sizes 9 ≤ N ≤ O (102); the same also applies for N = Np+3,Np+4, ... (iii) facile diffusion but with large β > 2 for N = Np+1 andNp+ 2 also for moderate sizes 9 ≤ N ≤ O (102); (iv) merging of the above distinct branches and subsequent anomalous scaling with 1 ≲ β < 3/2, reflecting the quasifacetted structure of clusters, for larger N = O (102) to N = O (103) ; (v) classic scaling with β = 3/2 for very large N = O (103) and above. The specified size ranges apply for typical model parameters. We focus on the moderate size regime where we show that diffusivity cycles quasiperiodically from the slowest branch forNp+ 3 (notNp) to the fastest branch forNp+ 1. Behavior is quantified by kinetic Monte Carlo simulation of an appropriate stochastic lattice-gas model. However, precise analysis must account for a strong enhancement of diffusivity for short time increments due to back correlation in the cluster motion. Further understanding of this enhancement, of anomalous size scaling behavior, and of the merging of various branches, is facilitated by combinatorial analysis of the number of the ground-state and low-lying excited state cluster configurations, and also of kink populations. [ABSTRACT FROM AUTHOR]

In order to gain insight into the nature of chemical bonding of sulfur atoms on coinage metal surfaces, we compare the adsorption energy and structural parameters for sulfur at four-fold hollow (4fh) sites on (100) facets and at three-fold hollow (3fh) sites on (111) facets of Cu nanoclusters. Consistent results are obtained from localized atomic orbital and plane-wave based density functional theory using the same functionals. PBE and its hybrid counterpart (PBE0 or HSE06) also give similar results. 4fh sites are preferred over 3fh sites with stronger bonding by ~0.6 eV for nanocluster sizes above ~280 atoms. However, for smaller sizes there are strong variations in the binding strength and the extent of the binding site preference. We show that suitable averaging over clusters of different sizes, or smearing the occupancy of orbitals, provide useful strategies to aid assessment of the behavior in extended surface systems. From site-projected density of states analysis using the smearing technique, we show that S adsorbed on a 4fh site has similar bonding interactions with the substrate as that on a 3fh site, but with much weaker antibonding interactions. [ABSTRACT FROM AUTHOR]

Confinement of electrons can occur in metal islands or in continuous films grown heteroepitaxially upon a substrate of a different metal or on a metallic alloy. Associated quantum size effects (QSE) can produce a significant height-dependence of the surface free energy for nanoscale thicknesses of up to 10-20 layers. This may suffice to induce height selection during film growth. Scanning STM analysis has revealed remarkable flat-topped or mesa-like island and film morphologies in various systems. We discuss in detail observations of QSE and associated film growth behavior for Pb/Cu(111), Ag/Fe(100), and Cu/fcc-Fe/Cu(100) [A/B or A/B/A], and for Ag/NiAl(110) with brief comments offered for Fe/Cu3Au(001) [A/BC binary alloys]. We also describe these issues for Ag/5-fold i-Al-Pd-Mn and Bi/5-fold i-Al-Cu-Fe [A/BCD ternary icosohedral quasicrystals]. Electronic structure theory analysis, either at the level of simple free electron gas models or more sophisticated Density Functional Theory calculations, can provide insight into the QSE-mediated thermodynamic driving force underlying height selection. [ABSTRACT FROM AUTHOR]