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). We find that H-H and H-Pd interactions promote the formation of the well-known 3×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.
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.
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
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 t o 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 i ncrease the activation of ice nucleation on feldspar dust. This included add ing a new aerosol tracer into the model and updat ing the ice nucleation parameterization, based on Classical Nucleation Theory, for multiple mineral dust tracers. Although t he se m odifications 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. Here 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.
Leite, Marina S.; Ruzmetov, Dmitry; Li, Zhipeng; Bendersky, Leonid A.; Bartelt, Norman C.; Kolmakov, Andrei; Talin, A.A.
The atomistic mechanism for lithiation/delithiation in all-solid-state batteries is still an open question, and the 'holy grail' to engineer devices with extended lifetime. Here, by combining real-time scanning electron microscopy in ultra-high vacuum with electrochemical cycling, we quantify the dynamic degradation of Al anodes in Li-ion all-solid-state batteries, a promising alternative for ultra lightweight devices. We find that AlLi alloy mounds are formed on the top surface of the Al anode and that degradation of battery capacity occurs because of Li trapped in them. Our approach establishes a new platform for probing the real-time degradation of electrodes, and can be expanded to other complex systems, allowing for high throughput characterization of batteries with nanoscale resolution.
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{sup +} and Ne{sup +}). 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.
We have used scanning tunneling microscopy and low-energy electron microscopy to measure the thermal decay of two-dimensional Cu, Pb-overlayer, and Pb-Cu alloy islands on Pb-Cu(1 1 1) surface alloys. Decay rates covering 6-7 orders of magnitude are accessible by applying the two techniques to the same system. We find that Cu adatom diffusion across the surface alloy is rate-limiting for the decay of both Pb and Pb-Cu islands on the surface alloy and that this rate decreases monotonically with increasing Pb concentration in the alloy. The decrease is attributed to repulsive interactions between Cu adatoms and embedded Pb atoms in the surface alloy. The measured temperature dependences of island decay rates are consistent with first-principles calculations of the Cu binding and diffusion energies related to this 'site-blocking' effect.
We find that small temperature changes cause steps on the NiAl(110) surface to move. We show that this step motion occurs because mass is transferred between the bulk and the surface as the concentration of bulk thermal defects (i.e., vacancies) changes with temperature. Since the change in an island's area with a temperature change is found to scale strictly with the island's step length, the thermally generated defects are created (annihilated) very near the surface steps. To quantify the bulk/surface exchange, we oscillate the sample temperature and measure the amplitude and phase lag of the system response, i.e., the change in an island's area normalized to its perimeter. Using a one-dimensional model of defect diffusion through the bulk in a direction perpendicular to the surface, we determine the migration and formation energies of the bulk thermal defects. During surface smoothing, we show that there is no flow of material between islands on the same terrace and that all islands in a stack shrink at the same rate. We conclude that smoothing occurs by mass transport through the bulk of the crystal rather than via surface diffusion. Based on the measured relative sizes of the activation energies for island decay, defect migration, and defect formation, we show that attachment/detachment at the steps is the rate-limiting step in smoothing.
The length scale of stress domain patterns formed at solid surfaces is usually calculated using isotropic elasticity theory. Because this length depends exponentially on elastic constants; deviations between isotropic and anisotropic elasticity can lead to large errors. Another inaccuracy of isotropic elasticity theory is that it neglects the dependence of elastic relaxations on stripe orientation. To remove these inaccuracies; we calculate the energy of striped domain patterns using anisotropic elasticity theory for an extensive set of surfaces encountered in experimental studies of self-assembly. We present experimental and theoretical evidence that elastic anisotropy is large enough to determine the stripe orientation when Pb is deposited on Cu(111). Our analytical and numerical results should be useful for analysis of a broad range of experimental systems.
Pb deposition on Cu(111) causes the surface to self-assemble into periodically arranged domains of a Pb-rich phase and a Pb-poor phase. Using low-energy electron microscopy (LEEM) we provide evidence that the observed temperature-dependent periodicity of these self-assembled domain patterns is the result of changing domain-boundary free energy. We determine the free energy of boundaries at different temperatures from a capillary wave analysis of the thermal fluctuations of the boundaries and find that it varies from 22 meV/nm at 600 K to 8 meV/nm at 650 K. Combining this result with previous measurements of the surface stress difference between the two phases we find that the theory of surface-stress-induced domain formation can quantitatively account for the observed periodicities.
We use low-energy electron microscopy to study the mechanisms of thermal smoothing on Rh(001) surfaces at high temperature. By examining the change of areas of two-dimensional islands as a function of time and temperature, we find that smoothing from 1210 K to 1450 K is limited by the rate of surface diffusion on terraces and not by bulk vacancy diffusion as observed in other systems in the same temperature range. However, the activation energy we measure for island decay is inconsistent with previous measurements and calculations of the activation energy of surface diffusion and the adatom formation energy. This inconsistency combined with an unexpectedly large activation entropy suggests a surface transport mechanism other than simple hopping of adatoms across the surface.
The energetics and thermal motion of the self-assembled domain structures of lead on copper were discussed. It was found that the self-assembled patterns arose from a temperature-independent surface stress difference of approximately 1.2 N/m. The domain patterns evolved in a manner consistent with models, when the lead coverage was increased.
Diffusion of Pd in the Pd/Cu(100) surface alloy was studied as a function of temperature using atom tracking scanning tunneling microscopy. It was shown that Pd diffused by a vacany-exchange mechanism, by following the motion of individual Pd atom incorporated in the surface. Activation energy was measured for diffusion of incorporated Pd atom and was found to be in agreement with ab intio calculated energy.