Broad background in electron diffraction of 2D materials as a signature of their superior quality.

An unusually broad bell-shaped component (BSC) has been previously observed in surface electron diffraction on different types of 2D systems. It was suggested to be an indicator of uniformity of epitaxial graphene (Gr) and hexagonal boron nitride (hBN). In the current study we use low-energy electron microscopy and micro-diffraction to directly relate the BSC to the crystal quality of the diffracting 2D material. Specially designed lateral heterostructures were used to map the spatial evolution of the diffraction profile across different 2D materials, namely pure hBN, BCN alloy and pure Gr, where the alloy region exhibits deteriorated structural coherency. The presented results show that the BSC intensity has a minimum in the alloyed region, consequently showing that BSC is sensitive to the lateral domain size and homogeneity of the material under examination. This is further confirmed by the presence of a larger number of sharp moiré spots when the BSC is most pronounced in the pure hBN and Gr regions. Consequently, it is proposed that the BSC can be used as a diagnostic tool for determining the quality of the 2D materials.

Encapsulation of metal nanoparticles at the surface of a prototypical layered material.

Encapsulation of metal nanoparticles just below the surface of a prototypical layered material, graphite, is a recently discovered phenomenon. These encapsulation architectures have potential for tuning the properties of two-dimensional or layered materials, and additional applications might exploit the properties of the encapsulated metal nanoclusters themselves. The encapsulation process produces novel surface nanostructures and can be achieved for a variety of metals. Given that these studies of near-surface intercalation are in their infancy, these systems provide a rich area for future studies. This Review presents the current progress on the encapsulation, including experimental strategies and characterization, as well as theoretical understanding which leads to the development of predictive capability. The Review closes with future opportunities where further understanding of the encapsulation is desired to exploit its applications.

Enhanced Nanostructure Dynamics on Au(111) with Adsorbed Sulfur due to Au-S Complex Formation.

Chemisorbed species can enhance the fluxional dynamics of nanostructured metal surfaces which has implications for applications such as catalysis. Scanning tunneling microscopy studies at room temperature reveal that the presence of adsorbed sulfur (S) greatly enhances the decay rate of 2D Au islands in the vicinity of extended step edges on Au(111). This enhancement is already significant at S coverages, θS , of a few hundredths of a monolayer (ML), and is most pronounced for 0.1-0.3 ML where the decay rate is increased by a factor of around 30. For θS close to saturation at about 0.6 ML, sulfur induces pitting and reconstruction of the entire surface, and Au islands are stabilized. Enhanced coarsening at lower θS is attributed to the formation and diffusion across terraces of Au-S complexes, particularly AuS2 and Au4 S4 , with some lesser contribution from Au3 S4 . This picture is supported by density functional theory analysis of complex formation energies and diffusion barriers.

Structure of chalcogen overlayers on Au(111): Density functional theory and lattice-gas modeling.

Ordering of different chalcogens, S, Se, and Te, on Au(111) exhibit broad similarities but also some distinct features, which must reflect subtle differences in relative values of the long-range pair and many-body lateral interactions between adatoms. We develop lattice-gas (LG) models within a cluster expansion framework, which includes about 50 interaction parameters. These LG models are developed based on density functional theory (DFT) analysis of the energetics of key adlayer configurations in combination with the Monte Carlo (MC) simulation of the LG models to identify statistically relevant adlayer motifs, i.e., model development is based entirely on theoretical considerations. The MC simulation guides additional DFT analysis and iterative model refinement. Given their complexity, development of optimal models is also aided by strategies from supervised machine learning. The model for S successfully captures ordering motifs over a broader range of coverage than achieved by previous models, and models for Se and Te capture the features of ordering, which are distinct from those for S. More specifically, the modeling for all three chalcogens successfully explains the linear adatom rows (also subtle differences between them) observed at low coverages of ∼0.1 monolayer. The model for S also leads to a new possible explanation for the experimentally observed phase with a (5 × 5)-type low energy electron diffraction (LEED) pattern at 0.28 ML and to predictions for LEED patterns that would be observed with Se and Te at this coverage.

Identification of an AgS2 Complex on Ag(110).

Adsorbed sulfur has been investigated on the Ag(110) surface at two different coverages, 0.02 and 0.25 monolayers. At the lower coverage, only sulfur adatoms are present. At the higher coverage, there are additional bright features which we identify as linear, independent AgS2 complexes. This identification is based upon density functional theory (DFT) and its comparison with experimental observations including bias dependence and separation between complexes. DFT also predicts the absence of AgS2 complexes at low coverage, and the development of AgS2 complexes around a coverage of 0.25 monolayers of sulfur, as is experimentally observed. To our knowledge, this is the first example of an isolated linear sulfur-metal-sulfur complex.

Sulfur adsorption on coinage metal(100) surfaces: propensity for metal-sulfur complex formation relative to (111) surfaces.

Experimental data from low-temperature Scanning Tunneling Microscopy (LTSTM) studies on coinage metal surfaces with very low coverages of S is providing new insights into metal-S interactions. A previous LTSTM study for Cu(100), and a new analysis reported here for Ag(100), both indicate no metal-sulfur complex formation, but an Au4S5 complex was observed previously on Au(100). In marked contrast, various complexes have been proposed and/or observed on Ag(111) and Cu(111), but not on Au(111). Also, exposure to trace amounts of S appears to enhance mass transport far more dramatically on (111) than on (100) surfaces for Cu and Ag, a feature tied to the propensity for complex formation. Motivated by these observations, we present a comprehensive assessment at the level of DFT to assess the existence and stability of complexes on (100) surfaces, and compare results with previous analyses for (111) surfaces. Consistent with experiment, our DFT analysis finds no stable complexes on Ag(100) and Cu(100), but several exist for Au(100). In addition, we systematically relate stability for adsorbed and gas-phase species within the framework of Hess's law. We thereby provide key insight into the various energetic contributions to stability which in turn elucidates the difference in behavior between (100) and (111) surfaces.

Reshaping, Intermixing, and Coarsening for Metallic Nanocrystals: Nonequilibrium Statistical Mechanical and Coarse-Grained Modeling.

Self-assembly of supported 2D or 3D nanocrystals (NCs) by vacuum deposition and of 3D NCs by solution-phase synthesis (with possible subsequent transfer to a support) produces intrinsically nonequilibrium systems. Individual NCs can have far-from-equilibrium shapes and composition profiles. The free energy of NC ensembles is lowered by coarsening which can involve Ostwald ripening or Smoluchowski ripening (NC diffusion and coalescence). Preservation of individual NC structure and inhibition of coarsening are key, e.g., for avoiding catalyst degradation. This review focuses on postsynthesis evolution of metallic NCs. Atomistic-level modeling typically utilizes stochastic lattice-gas models to access appropriate time and length scales. However, predictive modeling requires incorporation of realistic rates for relaxation mechanisms, e.g., periphery diffusion and intermixing, in numerous local environments (rather than the use of generic prescriptions). Alternative coarse-grained modeling must also incorporate appropriate mechanisms and kinetics. At the level of individual NCs, we present analyses of reshaping, including sintering and pinch-off, and of compositional evolution in a vacuum environment. We also discuss modeling of coarsening including diffusion and decay of individual NCs and unconventional coarsening processes. We describe high-level modeling integrated with scanning tunneling microscopy (STM) studies for supported 2D epitaxial nanoclusters and developments in modeling for 3D NCs motivated by in situ transmission electron microscopy (TEM) studies.

Characteristics of sulfur atoms adsorbed on Ag(100), Ag(110), and Ag(111) as probed with scanning tunneling microscopy: experiment and theory.

In this paper, we report that S atoms on Ag(100) and Ag(110) exhibit a distinctive range of appearances in scanning tunneling microscopy (STM) images, depending on the sample bias voltage, VS. Progressing from negative to positive VS, the atomic shape can be described as a round protrusion surrounded by a dark halo (sombrero) in which the central protrusion shrinks, leaving only a round depression. This progression resembles that reported previously for S atoms on Cu(100). We test whether DFT can reproduce these shapes and the transition between them, using a modified version of the Lang-Tersoff-Hamann method to simulate STM images. The sombrero shape is easily reproduced, but the sombrero-depression transition appears only for relatively low tunneling current and correspondingly realistic tip-sample separation, dT, of 0.5-0.8 nm. Achieving these conditions in the calculations requires sufficiently large separation (vacuum) between slabs, together with high energy cutoff, to ensure appropriate exponential decay of electron density into vacuum. From DFT, we also predict that an analogous transition is not expected for S atoms on Ag(111) surfaces. The results are explained in terms of the through-surface conductance, which defines the background level in STM, and through-adsorbate conductance, which defines the apparent height at the point directly above the adsorbate. With increasing VS, for Ag(100) and Ag(110), we show that through-surface conductance increases much more rapidly than through-adsorbate conductance, so the apparent adsorbate height drops below background. In contrast, for Ag(111) the two contributions increase at more comparable rates, so the adsorbate level always remains above background and no transition is seen.

Squeezed nanocrystals: equilibrium configuration of metal clusters embedded beneath the surface of a layered material.

Shapes of functional metallic nanocrystals, typically synthesized either free in solution or supported on surfaces, are key for controlling properties. Here, we consider a novel new class of metallic nanocrystals, copper clusters embedded near the surface of graphite, which can be considered a model system for metals embedded beneath surfaces of layered materials, or beneath supported membranes. We develop a continuum elasticity (CE) model for the equilibrium shape of these islands, and compare its predictions with experimental data. The CE model incorporates appropriate surface energy, adhesion energies, and strain energy. The agreement between the CE model and the data is-with one exception-excellent, both qualitatively and quantitatively, and is achieved with a single adjustable parameter. The model predicts that the embedded island shape is invariant with size, manifest both by constant side slope and by constant aspect ratio. This prediction is rationalized by dimensional analysis of the relevant energetic contributions. The aspect ratio (width : height) of an embedded Cu cluster is much larger than that of a supported but non-embedded Cu cluster, due to resistance of the graphene membrane to deformation. Experimental data diverge from the model predictions only in the case of the aspect ratio of small islands, below a critical height of ∼10 nm. The divergence may be due to bending strain, which is treated only approximately in the model. Strong support for the CE model and its interpretation is provided by additional data for embedded Fe clusters. Most of these observations and insights should be generally applicable to systems where a metal cluster is embedded beneath a layered material or supported membrane, provided that shape equilibration is possible.

Reverse-engineering of graphene on metal surfaces: a case study of embedded ruthenium.

Using scanning tunneling microscopy, x-ray photoelectron spectroscopy, and x-ray absorption spectroscopy, we show that Ru forms metallic nanoislands on graphite, covered by a graphene monolayer. These islands are air-stable, contain 2-4 layers of Ru, and have diameters on the order of 10 nm. To produce these nanoislands two conditions must be met during synthesis. The graphite surface must be ion-bombarded, and subsequently held at an elevated temperature (1000-1180 K) during Ru deposition. A coincidence lattice forms between the graphene overlayer and the Ru island top. Its characteristics-coincidence lattice constant, corrugation amplitude, and variation of carbon lattice appearance within the unit cell-closely resemble the well-established characteristics of single-layer graphene on the (0001) surface of bulk Ru. Quantitative analysis of the graphene lattice in relation to the coincidence lattice on the island tops show that the two-dimensional lattice constant of the underlying metal equals that of bulk Ru(0001), within experimental error. The embedded Ru islands are energetically favored over on-top (adsorbed) islands, based on density-functional-theory calculations for Ru films with 1-3 Ru layers. We propose a formation mechanism in which Ru atoms intercalate via defects that act as entry portals to the carbon galleries, followed by nucleation and growth in the galleries. In this model, high deposition temperature is necessary to prevent blockage of entry portals.

Oxygen and sulfur adsorption on vicinal surfaces of copper and silver: Preferred adsorption sites.

We present an extensive density functional theory (DFT) study of adsorption site energetics for oxygen and sulfur adsorbed on two vicinal surfaces of Cu and Ag, with the goal of identifying the most stable adsorption site(s), identifying trends and common themes, and comparing with experimental work in the literature where possible. We also present benchmark calculations for adsorption on the flat (111) and (100) surfaces. The first vicinal surface is the (211), and results are similar for both metals. We find that the step-doubling reconstruction is favored with both adsorbates and is driven by the creation of a special stable fourfold hollow (4fh) site at the reconstructed step. Zig-zag chain structures consisting of X-M-X units (X = chalcogen, M = metal) at the step edge are considered, in which the special 4fh site is partially occupied. The zig-zag configuration is energetically competitive for oxygen but not sulfur. DFT results for oxygen agree with experiment in terms of the stability of the reconstruction, but contradict the original site assignment. The second vicinal surface is the (410), where again results are similar for both metals. For oxygen, DFT predicts that step sites are filled preferentially even at lowest coverage, followed by terrace sites, consistent with the experiment. For sulfur, in contrast, DFT predicts that terrace sites fill first. Oxygen forms O-M-O rows on the top edge of the step, where it occupies incomplete 4fh sites. This resolves an experimental ambiguity in the site assignment. For both the (211) and (410) surfaces, the interaction energy that stabilizes the X-M-X chain or row correlates with the linearity of the X-M-X unit, which may explain key differences between oxygen and sulfur.

Sulfur Atoms Adsorbed on Cu(100) at Low Coverage: Characterization and Stability against Complexation.

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 × 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).

Comparison of S-adsorption on (111) and (100) facets of Cu nanoclusters.

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.

Interplay between surface and surface resonance states on height selective stability of fcc Dy(111) film at nanoscale.

Using first-principles calculations we show that face-centered cubic Dy(111) ultrathin films exhibit height selective stability. The origin of such height selection can be attributed to the interplay between the localized surface states and surface resonance states due to electron confinement effects. Such effect could be utilized to manipulate the film thickness at the atomic level to achieve desirable film properties or to control the growth of nanostructures on the thin film for various applications.

Formation of Two-Dimensional Copper Selenide on Cu(111) at Very Low Selenium Coverage.

Using scanning tunneling microscopy (STM), we observed that adsorption of Se on Cu(111) produced islands with a (√3×√3)R30° structure at Se coverages far below the structure's ideal coverage of 1/3 monolayer. On the basis of density functional theory (DFT), these islands cannot form due to attractive interactions between chemisorbed Se atoms. DFT showed that incorporating Cu atoms into the √3-Se lattice stabilizes the structure, which provided a plausible explanation for the experimental observations. STM revealed three types of √3 textures. We assigned two of these as two-dimensional layers of strained CuSe, analogous to dense planes of bulk klockmannite (CuSe). Klockmannite has a bulk lattice constant that is 11 % shorter than √3 times the surface lattice constant of Cu(111). This offers a rationale for the differences observed between these textures, for which strain limits the island size or distorts the √3 lattice. STM showed that existing step edges adsorb Se and facet toward ⟨12‾ 1⟩, which is consistent with DFT.

Adsorption of dysprosium on the graphite (0001) surface: Nucleation and growth at 300 K.

We have studied nucleation and growth of Dy islands on the basal plane of graphite at 300 K using scanning tunneling microscopy, density functional theory (DFT) in a form that includes van der Waals interactions, and analytic theory. The interaction of atomic Dy with graphite is strong, while the diffusion barrier is small. Experiment shows that at 300 K, the density of nucleated islands is close to the value predicted for homogeneous nucleation, using critical nucleus size of 1 and the DFT-derived diffusion barrier. Homogeneous nucleation is also supported by the monomodal shape of the island size distributions. Comparison with the published island density of Dy on graphene shows that the value is about two orders of magnitude smaller on graphite, which can be attributed to more effective charge screening in graphite. The base of each island is 3 atomic layers high and atomically ordered, forming a coincidence lattice with the graphite. Islands resist coalescence, probably due to multiple rotational orientations associated with the coincidence lattice. Upper levels grow as discernible single-atom layers. Analysis of the level populations reveals significant downward interlayer transport, which facilitates growth of the base. This island shape is metastable, since more compact three-dimensional islands form at elevated growth temperature.

Preface: special topic on supramolecular self-assembly at surfaces.

Supramolecular self-assembly at surfaces is one of the most exciting and active fields in Surface Science today. Applications can take advantage of two key properties: (i) versatile pattern formation over a broad length scale and (ii) tunability of electronic structure and transport properties, as well as frontier orbital alignment. It provides a new frontier for Chemical Physics as it uniquely combines the versatility of Organic Synthesis and the Physics of Interfaces. The Journal of Chemical Physics is pleased to publish this Special Topic Issue, showcasing recent advances and new directions.

The (111) surface of NaAu2: structure, composition, and stability.

The (111) surface of single-crystal NaAu(2) is a model for catalytically active, powdered NaAu(2). We prepare and characterize this surface with a broad suite of techniques. Preparation in ultrahigh vacuum consists of the traditional approach of ion bombardment (to remove impurities) and thermal annealing (to restore surface order). Both of these steps, however, cause loss of sodium (Na), and repeated treatments eventually trigger conversion of the surface and near-surface regions to crystalline gold. The bulk has a limited ability to repopulate the surface Na. Under conditions where Na depletion is minimized, electron diffraction patterns are consistent with the bulk-terminated structure, and scanning tunneling microscopy reveals mesa-like features with lateral dimensions of a few tens of nanometers. The tops of the mesas do not possess fine structure characteristic of a periodic lattice, suggesting that the surface layer is disordered under the conditions of these experiments.

Structure of the clean Gd5Ge4(010) surface.

We have characterized the (010) surface of Gd5Ge4 using scanning tunneling microscopy (STM) and x-ray photoelectron spectroscopy. Data from different samples have the following features in common: (1) the surface composition equals the bulk composition to within 5 at.%, both after ion etching and after annealing at temperatures of 400-1200 K; and (2) the surface exhibits terraces of two types. The height of the steps between similar terraces corresponds well to the separation between equivalent layers along the <010> direction in the bulk structure. Density functional theory (DFT) shows that the surface energy of the (0001) plane of hexagonal close-packed Gd is lower than that of the (111) plane of diamond-type Ge, suggesting that surfaces of Gd5Ge4 (for comparable density) should be rich in Gd. Indeed, DFT shows that among the bulk terminations of Gd5Ge4, a pure Ge termination is not favored. Each of the three remaining terminations (two pure Gd and one mixed, Gd-Ge) has its minimum surface energy in a different range of the possible Gd chemical potentials, indicating that different terminations may be stable under different conditions. DFT shows that the heights of the steps between dissimilar terraces, measured in STM, are consistent with the two pure Gd terminations.