Many technologies require stable or metastable surface morphology. In this paper we study the factors that control the metastability of a common feature of rough surfaces: "hillocks."We use low energy electron microscopy to follow the evolution of the individual atomic steps in hillocks on Pd(111). We show that the uppermost island in the stack often adopts a static, metastable configuration. Modeling this result shows that the degree of the metastability depends on the configuration of steps dozens of atomic layers lower. Our model allows us to link surface metastability to the atomic processes of surface evolution.
Efficient prediction of sampling-intensive thermodynamic properties is needed to evaluate material performance and permit high-throughput materials modeling for a diverse array of technology applications. To alleviate the prohibitive computational expense of high-throughput configurational sampling with density functional theory (DFT), surrogate modeling strategies like cluster expansion are many orders of magnitude more efficient but can be difficult to construct in systems with high compositional complexity. We therefore employ minimal-complexity graph neural network models that accurately predict and can even extrapolate to out-of-train distribution formation energies of DFT-relaxed structures from an ideal (unrelaxed) crystallographic representation. This enables the large-scale sampling necessary for various thermodynamic property predictions that may otherwise be intractable and can be achieved with small training data sets. Two exemplars, optimizing the thermodynamic stability of low-density high-entropy alloys and modulating the plateau pressure of hydrogen in metal alloys, demonstrate the power of this approach, which can be extended to a variety of materials discovery and modeling problems.
The fabrication of long-lived electrical contacts to thermoelectric Bi2Te3-based modules is a challenging problem due to chemical incompatibilities and rapid diffusion rates. Previously, technical guidance from SAND report 2015-7203 selected electroplated Au as the preferred method for fabrication of long-lived contacts because of concerns that the grain structure of sputtered/physical vapor deposited (PVD) Au contacts can evolve during aging. We have re-evaluated PVD Au contacts and show that they are appropriate for long-life service. We measure grain size and morphology at different aging times under accelerated temperature gradient conditions, and we show that the PVD Au contacts are stable and remain relatively unchanged. The PVD Au fabricated here is not subject to the deterioration observed in the previous report.
Understanding hydrogen incorporation into palladium requires detailed knowledge of surface and subsurface structure and atomic interactions as surface hydrogen is being embedded. Using density functional theory (DFT), we examine the energies of hydrogen layers of varying coverage adsorbed on Pd(111). Here we find that H–H and H–Pd interactions promote the formation of the well-known ($\sqrt{3}$ x $\sqrt{3}$) phases but also favor an unreported (3 × 3) phase at high H coverages for which we present experimental evidence. We relate the stability of isolated H vacancies of the (3 × 3) phase to the need of H2 molecules to access bare Pd before they can dissociate. Following higher hydrogen dosage, we observe initial steps of hydride formation, starting with small clusters of subsurface hydrogen. The interaction between H and Pd is complicated by the persistent presence of carbon at the surface. X-ray photoelectron spectroscopy experiments show that trace amounts of carbon, emerging from the Pd bulk despite many surface cleaning cycles, become mobile enough to repopulate the C-depleted surface at temperatures above 200 K. When exposed to hydrogen, these surface carbon atoms react to form benzene, as evidenced by scanning tunneling microscopy observations interpreted with DFT.
We study the structure of the threading edge dislocations, or “elbows,” which are an essential component of the well-known herringbone reconstruction of the (111) surface of Au. Previous work had shown that these dislocations can be stabilized by long-range elastic relaxations into the bulk. However, the validity of the harmonic spring model that had been used to estimate the energies of the dislocations is uncertain. To enable a more refined model of the dislocation energetics, we have imaged the atomic structure of these dislocations using scanning tunneling microscopy. We find that the harmonic spring model does not adequately reproduce the observed structure. We are able to reproduce the structure, however, with a two-dimensional Frenkel-Kontorova (FK) model that uses a pairwise Morse potential to describe the interactions between the top layer Au atoms on a rigid substrate. The parameters of the potential were obtained by fitting the energy of uniaxially compressed phases, or “stripes”, computed with density functional theory, as a function of surface Au density. Within this model, the formation of the threading dislocations remains unfavorable. However, the large forces on the substrate atoms near the threading-dislocation cores, render the assumption of a completely rigid substrate questionable. Indeed, if the FK parameters are modified to account for the relaxation of just one more atomic layer, threading dislocations can, in principle, become favorable, even without bulk elastic relaxations. Additional evidence for a small elbow energy is that our computed change in the Au(111) surface stress tensor caused by the (√ 3 × 22) reconstruction is considerably smaller than previous estimates.
The understanding and control of charge carrier interactions with defects at buried insulator/semiconductor interfaces is essential for achieving optimum performance in modern electronics. Here, we report on the use of scanning ultrafast electron microscopy (SUEM) to remotely probe the dynamics of excited carriers at a Si surface buried below a thick thermal oxide. Our measurements illustrate a previously unidentified SUEM contrast mechanism, whereby optical modulation of the space-charge field in the semiconductor modulates the electric field in the thick oxide, thus affecting its secondary electron yield. By analyzing the SUEM contrast as a function of time and laser fluence we demonstrate the diffusion mediated capture of excited carriers by interfacial traps.
We successfully demonstrated the utility of surface science techniques - namely scanning probe microscopy and thermal desorption spectroscopy - on three different material systems: incipient soot formed during fossil fuel combustion, surface oxides passivating polycrystalline nickel hydrogen uptake, and aluminum hydride cluster formation underpinning solid-state hydrogen fuel storage. For all three material systems, surface science techniques haven proven to probe intricate nanoscale phenomena that are critical to macroscale material behavior. This LDRD has gained insight into early-stage pollution formation, the impacts of common contaminants on tritium flow regulation, and the limitations of solid-state hydrogen fuel storage. Our results support the diversification of national energy technologies.
The predictive understanding of catalytic surface reactions requires accurate microkinetic models, and while decades of work has been devoted to the elucidation of the reaction steps in these models, many open questions remain. One key issue is a lack of approaches enabling the local spatially resolved assessment of catalytic activity over a surface. In this report, we detail efforts to develop a new diagnostic approach to solve this problem. The approach is based upon laser resonance enhanced multiphoton ionization of reaction products emitted into the gas phase followed by spatially resolved imaging of the resultant ions or electrons. Ion imaging is pursued with a velocity-selected spatially resolved ion imaging microscope, while electron imaging was attempted in a low energy electron microscope. Successful demonstration of the ion imaging microscope coupled with the development of transport simulations shows promise for a revolutionary new tool to assess local catalytic activity
This project aimed to identify the performance-limiting mechanisms in mid- to far infrared (IR) sensors by probing photogenerated free carrier dynamics in model detector materials using scanning ultrafast electron microscopy (SUEM). SUEM is a recently developed method based on using ultrafast electron pulses in combination with optical excitations in a pump- probe configuration to examine charge dynamics with high spatial and temporal resolution and without the need for microfabrication. Five material systems were examined using SUEM in this project: polycrystalline lead zirconium titanate (a pyroelectric), polycrystalline vanadium dioxide (a bolometric material), GaAs (near IR), InAs (mid IR), and Si/SiO 2 system as a prototypical system for interface charge dynamics. The report provides detailed results for the Si/SiO 2 and the lead zirconium titanate systems.
Pd readily absorbs hydrogen and its isotopes, and can be used to purify gas mixtures involving tritium. Tritium decays to He, forming He bubbles. Bubbles causes possible PCT effects swelling, He release, all leading to failures. Radioactive decay experiments take many years. Molecular dynamics (MD) studies can be quickly done. No previous MD methods can simulate He bubble nucleation and growth.
Pd readily absorbs hydrogen and its isotopes, and can be used to purify gas mixtures involving tritium. Tritium decays to He, forming He bubbles. Bubbles causes possible PCT effects swelling, He release, all leading to failures. Radioactive decay experiments take many years. Molecular dynamics (MD) studies can be quickly done. No previous MD methods can simulate He bubble nucleation and growth.
We have used scanning tunneling microscopy and density functional theory calculations to study molecular layers of coronene on Cu(111). The structure and stability of these layers is determined by the balance between coronene-substrate and coronene-coronene interactions. Here, we characterize this balance by measuring the maximum coverage before coronene dewets the substrate to form three-dimensional islands. We find that coronene molecules lie parallel to the substrate at the maximum coverage, in contrast to previous observations of tilted coronene on metal surfaces. We attribute this previously reported tilt to a metastability caused by an activation barrier to nucleate three-dimensional islands.
Atmospheric ice affects Earth's radiative properties and initiates most precipitation. Growing ice typically requires a particle, often airborne mineral dust, e.g., to catalyze freezing of supercooled cloud droplets. How chemistry, structure and morphology determine the ice-nucleating ability of minerals remains elusive. Not surprisingly, poor understanding of a erosol-cloud interactions is a major source of uncertainty in climate models. In this project, we combine d optical microscopy with atomic force microscopy to explore the mechanisms of initial ice formation on alkali feldspar, a mineral proposed to dominate ice nucleation in Earth's atmosphere. When cold air becomes supersaturated with respect to water, we discovered that supercooled liquid water condenses at steps without having to overcome a nucleation barrier, and subsequently freezes quickly. Our results imply that steps, common even on macroscopically flat feldspar surfaces, can accelerate water condensation followed by freezing, thus promoting glaciation and dehydration of mixed - phase clouds. Motivated by the fact that current climate simulations do not properly account for feldspar's extreme efficiency to nucleate ice, we modified DOE's climate model, the Energy Exascale Earth System Model (E3SM), to increase the activation of ice nucleation on feldspar dust. This included adding a new aerosol tracer into the model and updating the ice nucleation parameterization, based on Classical Nucleation Theory, for multiple mineral dust tracers. Although t he se modifications have little impact on global averages , predictions of regional averages can be strongly affected .
The earliest stages of soot formation in flames are believed to involve the formation of small, nanoscale clusters of polycyclic aromatic hydrocarbon molecules. The structure of these clusters is still highly uncertain, however, impeding the construction of quantitative models of soot inception and growth. To provide insight into the structure of incipient soot, we produced nanoclusters of hydrocarbon molecules by annealing coronene films deposited on Pt(111), and examined them with scanning tunneling microcopy. We find that clusters containing ∼20–100 molecules, are disordered agglomerations of stacks that are ∼5–6 molecules tall. These structures are quite distinct from crystalline coronene, but bear a striking resemblance to recently proposed models for the equilibrium structure of similarly-sized clusters that are assumed to initiate soot formation. In contrast to mature soot, the surfaces of these clusters contain very few molecules with graphitic planes oriented parallel to the surface.
Next generation metal-ion conducting membranes are key to developing energy storage and utilization technologies like batteries and fuel ce lls. Sodium super-ionic conductors (aka NaSICON) are a class of compounds with AM 1 M 2 (PO 4 ) 3 stoichiometry where the choice of "A" and "M" cation varies widely. This report, which de scribes substitutional derivatives of NZP (NaZr 2 P 3 O 12 ), summarizes the accomplishments of a Laboratory D irected Research and Development (LDRD) project to analyze transport mec hanisms using a combination of in situ studies of structure, composition, and bonding, com bined with first principles theory and modeling. We developed an experimental platform and applied methods, such as synchrotron- based X-ray spectroscopies, to probe the electronic structure of compositionally well-controlled NaSICON films while in operation ( i.e ., conducting Na ions exposed to oxygen or water va por atmospheres). First principles theory and modeling were used to interpret the experimental observations and develop an enhanced understanding of atomistic processes that give rise to, and affect, ion conduction.
Carbon nanostructures, such as nanotubes and graphene, are of considerable interest due to their unique mechanical and electrical properties. The materials exhibit extremely high strength and conductivity when defects created during synthesis are minimized. Atomistic modeling is one technique for high resolution studies of defect formation and mitigation. To enable simulations of the mechanical behavior and growth mechanisms of C nanostructures, a high-fidelity analytical bond-order potential for the C is needed. To generate inputs for developing such a potential, we performed quantum mechanical calculations of various C structures.
The electrochemical reactions of solid oxide fuel cells occur in the region where gas-phase species, electrode, and electrolyte coincide. When the electrode is an ionic insulator and the electrolyte is an electronic insulator, this triple phase boundary is assumed to have atomic dimensions. For this study, we use photoemission electron microscopy to show that the reduced surface of the electrolyte yttria-stabilized zirconia (YSZ) has a sharp electronic metal-insulator boundary near Pt negative electrodes. The electronic conductivity of the reduced YSZ allows for oxygen reduction, allowing the reduced YSZ to behave as an extended triple phase boundary. This extended triple phase boundary can be many microns in size, depending on oxygen pressure, temperature, applied voltage, and time.
In this project we developed t he atomistic models needed to predict how graphene grows when carbon is deposited on metal and semiconductor surfaces. We first calculated energies of many carbon configurations using first principles electronic structure calculations and then used these energies to construct an empirical bond order potentials that enable s comprehensive molecular dynamics simulation of growth. We validated our approach by comparing our predictions to experiments of graphene growth on Ir, Cu and Ge. The robustness of ou r understanding of graphene growth will enable high quality graphene to be grown on novel substrates which will expand the number of potential types of graphene electronic devices.
Graphene films grown by vapour deposition tend to be polycrystalline due to the nucleation and growth of islands with different in-plane orientations. Here, using low-energy electron microscopy, we find that micron-sized graphene islands on Ir(111) rotate to a preferred orientation during thermal annealing. We observe three alignment mechanisms: the simultaneous growth of aligned domains and dissolution of rotated domains, that is, â €- ripeningâ €™; domain boundary motion within islands; and continuous lattice rotation of entire domains. By measuring the relative growth velocity of domains during ripening, we estimate that the driving force for alignment is on the order of 0.1â €‰meV per C atom and increases with rotation angle. A simple model of the orientation-dependent energy associated with the moiré corrugation of the graphene sheet due to local variations in the graphene-substrate interaction reproduces the results. This work suggests new strategies for improving the van der Waals epitaxy of 2D materials.
Increasingly, basic models such as density functional theory and molecular dynamics are being used to simulate different aspects of hydrogen recycling from plasma facing materials. These models provide valuable insight into hydrogen diffusion, trapping, and recombination from surfaces, but their validation relies on knowledge of the detailed behavior of hydrogen at an atomic scale. Despite being the first wall material for ITER, basic single crystal beryllium surfaces have been studied only sparsely from an experimental standpoint. In prior cases researchers used electron spectroscopy to examine surface reconstruction or adsorption kinetics during exposure to a hydrogen atmosphere. While valuable, these approaches lack the ability to directly detect the positioning of hydrogen on the surface. Ion beam techniques, such as low energy ion scattering (LEIS) and direct recoil spectroscopy (DRS), are two of the only experimental approaches capable of providing this information. In this study, we applied both LEIS and DRS to examine how hydrogen binds to the Be(0001) surface. Our measurements were performed using an angle-resolved ion energy spectrometer (ARIES) to probe the surface with low energy ions (500 eV - 3 keV He+ and Ne+). We were able to obtain a 'scattering maps' of the crystal surface, providing insight on how low energy ions are focused along open surface channels. Once we completed a characterization of the clean surface, we dosed the sample with atomic hydrogen using a heated tungsten capillary. A distinct signal associated with adsorbed hydrogen emerged that was consistent with hydrogen residing between atom rows. To aid in the interpretation of the experimental results, we developed a computational model to simulate ion scattering at grazing incidence. For this purpose, we incorporated a simplified surface model into the Kalypso molecular dynamics code. This approach allowed us to understand how the incident ions interacted with the surface hydrogen, providing confirmation of the preferred binding site.
Growth of high quality graphene films on SiC is regarded as one of the more viable pathways toward graphene-based electronics. Graphitic films form on SiC at elevated temperature because of preferential sublimation of Si. Little is known, however, about the atomistic processes of interrelated SiC decomposition and graphene growth. We have observed the formation of graphene on SiC by Si sublimation in an Ar atmosphere using low energy electron microscopy, scanning tunneling microcopy and atomic force microscopy. This work reveals that the growth mechanism depends strongly on the initial surface morphology, and that carbon diffusion governs the spatial relationship between SiC decomposition and graphene growth. Isolated bilayer SiC steps generate narrow ribbons of graphene, whereas triple bilayer steps allow large graphene sheets to grow by step flow. We demonstrate how graphene quality can be improved by controlling the initial surface morphology specifically by avoiding the instabilities inherent in diffusion-limited growth.
We have studied the dislocation structures that occur in films of Ag, Au, and Ag{sub 0.5}Au{sub 0.5} alloy on a Ru(0001) substrate. Monolayer (ML) films form herringbone phases while films two or more layers thick contain triangular patterns of dislocations. We use scanning tunneling microscopy (STM) and low-energy electron diffraction (LEED) to determine how the film composition affects the structure and periodicity of these ordered structures. One layer of Ag forms two different herringbone phases depending on the exact Ag coverage and temperature. Low-energy electron microscopy (LEEM) establishes that a reversible, first-order phase transition occurs between these two phases at a certain temperature. We critically compare our 1 ML Ag structures to conflicting results from an X-ray scattering study [H. Zajonz et al., Phys. Rev. B 67 (2003) 155417]. Unlike Ag, the herringbone phases of Au and AgAu alloy are independent of the exact film coverage. For two layer films in all three systems, none of the dislocations in the triangular networks thread into the second film layer. In all three systems, the in-plane atomic spacing of the second film layer is nearly the same as in the bulk. Film composition does, however, affect the details of the two layer structures. Ag and Au films form interconnected networks of dislocations, which we refer to as 'trigons.' In 2 ML AgAu alloy, the dislocations form a different triangular network that shares features of both trigon and moire structures. Yet another well-ordered structure, with square symmetry, forms at the boundaries of translational trigon domains in 2 ML Ag films but not in Au films.