Publications

Abstract not provided.

Characterizing the shallow structure of the Rock Valley region of the Nevada National Security Site is a critical component of the Rock Valley Direct Comparison project. Geophysical data of the region is needed for operational decisions, to constrain geologic models used for simulation, and to facilitate the analysis of future explosive source data. Local measurements of gravity are a key piece of geophysical information that helps to resolve the underlying geologic composition, fault structure, and density characteristics, yet, in the Rock Valley region these measurements are sparse on the scale of the testbed. In this report, we present the details of a recent gravity data acquisition survey designed to collect a dense dataset in the region of interest that complements the existing gravity work but greatly enhances our resolution. This dataset will be integrated with a complementary Los Alamos National Laboratory gravity collection and combined with the existing seismic data in a joint inversion. These measurements were conducted over two weeks with a portable gravimeter and high-resolution GPS and include repeat measurements at a USGS base station as well as reoccupation of gravity sites in the regional dataset. This collection of over 100 new dense gravity measurements will facilitate refinement of the existing Geologic Framework Model and directly complement newly acquired dense seismic data, ultimately improving the project’s ability to investigate the direct comparison of shallow earthquake and explosive sources.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The Energy Resilience for Mission Assurance (ERMA) project—a Department of Energy Grid Modernization Lab Consortium effort carried out via a partnership among five national laboratories— seeks to develop metrics to quantify how improvements to energy system resilience translate to improved Department of Defense (DoD) mission assurance (MA) during wide-scale, long-duration outages of the bulk power system. DoD missions are integral to national security and highly dependent on electric power. However, energy system planners—both civilian and military—lack a clear and quantifiable mapping between electric power system resilience and MA, leaving a gap in their ability to understand and consider national security outcomes within their planning efforts. The ERMA project seeks to fill this gap, providing stakeholders with new capabilities to understand the impact of electric power system resilience on MA during hazard scenarios.

Abstract not provided.

We propose primal–dual mesh optimization algorithms that overcome shortcomings of the standard algorithm while retaining some of its desirable features. “Hodge-Optimized Triangulations” defines the “HOT energy” as a bound on the discretization error of the diagonalized Delaunay Hodge star operator. HOT energy is a natural choice for an objective function, but unstable for both mathematical and algorithmic reasons: it has minima for collapsed edges, and its extrapolation to non-regular triangulations is inaccurate and has unbounded minima. We propose a different extrapolation with a stronger theoretical foundation, and avoid extrapolation by recalculating the objective just beyond the flip threshold. We propose new objectives, based on normalizations of the HOT energy, with barriers to edge collapses and other undesirable configurations. We propose mesh improvement algorithms coupling these. When HOT optimization nearly collapses an edge, we actually collapse the edge. Otherwise, we use the barrier objective to update positions and weights and remove vertices. By combining discrete connectivity changes with continuous optimization, we more fully explore the space of possible meshes and obtain higher quality solutions.

Abstract not provided.

Laboratory shear tests were conducted on pulverized Sierra White granite (SWG) to investigate slip mechanisms in naturally occurring faults. Synthetic fault geometries were constructed by sandwiching fine grained SWG powder in between steel forcing blocks. For dry experiments, ~3.5 g of SWG powder was poured onto the face of the lower steel forcing block and leveled. For saturated experiments, enough fluid was added to the ~3.5 g of Sierra White granite powder to form a slurry. This slurry was applied to the lower forcing block and leveled. Inclined forcing blocks with 25.4 mm diameter and 35° faces, which were machined from ground steel rods with fine teeth on the faces, help to hold the gouge in place and prevent delamination at the interface. The top forcing block had a 2.03 mm centered hole to allow pore fluid access to the gouge. A fine steel mesh prevented back flow of the gouge into pore fluid lines. Samples were isolated from the confining medium using three layers of heat shrink polyolefin, as shown in Figure 1. The outer layer was shrunk over the o-rings on the end caps to form an impermeable seal, which was reinforced with steel tie wires on both sides of the o-rings. Hardened steel spacers and copper shim stock was placed between the steel forcing blocks and the end caps to preserve the parallelism of the Hastelloy wetted parts. For dry samples, the end caps were plugged, while the end caps for the saturated samples were connected to pore lines.

The flow and particle deposition patterns on surfaces in an idealized spacer grid for a 17x17 pressurized water reactor (PWR) assembly in a spent fuel canister are modeled using computational fluid dynamics (CFD) with laminar flow. The effects of gravitational settling, non-Stokesian flow, and particle slip are first rigorously analyzed. From the analysis, non-Stokesian effects and slip may be neglected for the particle sizes and conditions expected in a canister. For particles that do not settle out, a swirling flow pattern at the corners of a spacer grid channel directs particles to the leeward side of the flow vanes where much of the deposition occurs. Particle deposition increases with increasing particle diameter. Deposition also increases with decreasing flow velocity as this provides more time for particles to settle and deposit on the leeward side of the flow vanes. The fraction of particles that are transmitted through a spacer grid is determined as a function of inlet gas velocity and particle diameter by running the CFD calculation for each set of conditions and for each particle diameter. Curve fits of the transmission curve as a function of particle diameter for a specified spacer grid and flow velocity are applied to a lognormal particle mass density function for the inlet particles. The resulting mass density function and aerosol mass fraction that passes through the spacer grid can be determined analytically without resorting to numerical iteration. A sample calculation of the analytical solution is demonstrated for a lognormal particle mass density function.

The goal of this work is to develop a wafer level copper electrochemical deposition process for high aspect ratio through silicon vias and test the reliability and functionality of TSVs for integration into MEMS applications.

In typical carbonyl-containing molecules, bond dissociation events follow initial excitation to $nπ_{C=O}$$^*$ states. However, in acetyl iodide, the iodine atom gives rise to electronic states with mixed $nπ_{C=O}$$^*$ and $nπ_{C–I}$$^*$ character, leading to complex excited-state dynamics, ultimately resulting in dissociation. Using ultrafast extreme ultraviolet (XUV) transient absorption spectroscopy and quantum chemical calculations, we present an investigation of the primary photodissociation dynamics of acetyl iodide via time-resolved spectroscopy of core-to-valence transitions of the I atom after 266 nm excitation. The probed I 4d-to-valence transitions show features that evolve on sub-100-fs time scales, reporting on excited-state wavepacket evolution during dissociation. These features subsequently evolve to yield spectral signatures corresponding to free iodine atoms in their spin–orbit ground and excited states with a branching ratio of 1.1:1 following dissociation of the C–I bond. Calculations of the valence excitation spectrum via equation-of-motion coupled cluster with single and double substitutions (EOM-CCSD) show that initial excited states are of spin-mixed character. From the initially pumped spin-mixed state, we use a combination of time-dependent density functional theory (TDDFT)-driven nonadiabatic ab initio molecular dynamics and EOM-CCSD calculations of the N$_{4,5}$ edge to reveal a sharp inflection point in the transient XUV signal that corresponds to rapid C–I homolysis. Here, by examining the molecular orbitals involved in the core-level excitations at and around this inflection point, we are able to piece together a detailed picture of C–I bond photolysis in which d → σ* transitions give way to d → p excitations as the bond dissociates. We also report theoretical predictions of short-lived, weak 4d → 5d transitions in acetyl iodide, validated by weak bleaching in the experimental transient XUV spectra. This joint experimental–theoretical effort has thus unraveled the detailed electronic structure and dynamics of a strongly spin–orbit coupled system.

The mechanical behavior of partial-penetration laser welds exhibits significant variability in engineering quantities such as strength and apparent ductility. Understanding the root cause of this variability is important when using such welds in engineering designs. In Part II of this work, we develop finite element simulations with geometry derived from micro-computed tomography (μCT) scans of partial-penetration 304L stainless steel laser welds that were analyzed in Part I. We use these models to study the effects of the welds’ small-scale geometry, including porosity and weld depth variability, on the structural performance metrics of weld ductility and strength under quasi-static tensile loading. We show that this small-scale geometry is the primary cause of the observed variability for these mechanical response quantities. Additionally, we explore the sensitivity of model results to the conversion of the μCT data to discretized model geometry using different segmentation algorithms, and to the effect of small-scale geometry simplifications for pore shape and weld root texture. The modeling approach outlined and results of this work may be applicable to other material systems with small-scale geometric features and defects, such as additively manufactured materials.

Recent studies of power flow and particle transport in multi-MA pulsed-power accelerators demonstrate that electrode plasmas may reduce accelerator efficiency by shunting current upstream from the load. The detailed generation and evolution of these electrode plasmas are examined here using fully relativistic, Monte Carlo particle-in-cell (PIC) and magnetohydrodynamic (MHD) simulations over a range of peak currents (8–48 MA). The PIC calculations, informed by vacuum science, describe the electrode surface breakdown and particle transport prior to electrode melt. The MHD calculations show the bulk electrode evolution during melt. The physical description provided by this combined study begins with the rising local magnetic field that increases the local electrode surface temperature. This initiates the thermal desorption of contaminants from the electrode surface, with contributions from atoms outgassing from the bulk metal. The contaminants rapidly ionize forming a 10¹⁵-10¹⁸ cm^-3 plasma that is effectively resistive while weakly collisional because it is created within, and rapidly penetrated by, a strong magnetic field (> 30 T). Prior to melting, the density of this surface plasma is limited by the concentration of absorbed contaminants in the bulk (~10¹⁹ cm^-3 for hydrogen), its diffusion, and ionization. Eventually, the melting electrodes form a conducting plasma (10²¹-10²³ cm^-3) that experiences j × B compression and a typical decaying magnetic diffusion profile. This physical sequence ignores the transport of collisional plasmas of 10¹⁹ cm^-3 which may arise from electrode defects and associated instabilities. Nonetheless, this picture of plasma formation and melt may be extrapolated to higher-energy pulsed-power systems.

Recent studies on off-stoichiometric thermosets reveal unique viscoelastic behavior derived from increased free volume and physical interactions between chain ends. To understand structural characteristics arising from cure and its effect on properties, we developed a Monte Carlo model based on step-growth polymerization. Our model accurately predicted structure-property trends for a two-component system of EPON 828 (EPON) and ethylenediamine. A second epoxy monomer, D.E.R. 732 (DER), was investigated to modulate Tg. Binary mixtures of EPON and DER in off-stoichiometric, amine-rich formulations resulted in nonlinear evolution of thermomechanical properties with respect to initial formulation stoichiometry. Modifying our model with kinetic parameters allowing for differential epoxide/amine reaction kinetics only partially accounted for trends in Tg, suggesting that spatiotemporal contributions─not captured by our model─were significant determinants of material properties compared to polymer architecture for three-component systems. These findings underpin the importance of spatial awareness in modeling to inform the development of dynamic thermosets.

A cloud of very fast, O(km/s), and very fine, O(µm), particles may be ejected when a strong shock impacts and possibly melts the free surface of a solid metal. To quantify these dynamics, this work develops an ultraviolet, long-working distance, two-pulse Digital Holographic Microscopy (DHM) configuration and is the first to replace film recording with digital sensors for this challenging application. A proposed multi-iteration DHM processing algorithm is demonstrated for automated measures of the sizes, velocities, and three-dimensional positions of non-spherical particles. Ejecta as small as 2 µm diameter are successfully tracked, while uncertainty simulations indicate that particle size distributions are accurately quantified for diameters ≥4 µm. These techniques are demonstrated on three explosively driven experiments. Measured ejecta size and velocity statistics are shown to be consistent with prior film-based recording, while also revealing spatial variations in velocities and 3D positions that have yet to be widely investigated. Having eliminated time-consuming analog film processing, the methodologies proposed here are expected to significantly accelerate future experimental investigation of ejecta physics.

Multimodal datasets of materials are rich sources of information which can be leveraged for expedited discovery of process–structure–property relationships and for designing materials with targeted structures and/or properties. For this data descriptor article, we provide a multimodal dataset of magnetron sputter-deposited molybdenum (Mo) thin films, which are used in a variety of industries including high temperature coatings, photovoltaics, and microelectronics. In this dataset we explored a process space consisting of 27 unique combinations of sputter power and Ar deposition pressure. Here, the phase, structure, surface morphology, and composition of the Mo thin films were characterized by x-ray diffraction, scanning electron microscopy, atomic force microscopy, and Rutherford backscattering spectrometry. Physical properties—namely, thickness, film stress and sheet resistance—were also measured to provide additional film characteristics and behaviors. Additionally, nanoindentation was utilized to obtain mechanical load-displacement data. The entire dataset consists of 2072 measurements including scalar values (e.g., film stress values), 2D linescans (e.g., x-ray diffractograms), and 3D imagery (e.g., atomic force microscopy images). An additional 1889 quantities, including film hardness, modulus, electrical resistivity, density, and surface roughness, were derived from the experimental datasets using traditional methods. Minimal analysis and discussion of the results are provided in this data descriptor article to limit the authors’ preconceived interpretations of the data. Overall, the data modalities are consistent with previous reports of refractory metal thin films, ensuring that a high-quality dataset was generated. The entirety of this data is committed to a public repository in the Materials Data Facility.

This paper focuses on the development of multifidelity modeling approaches using neural network surrogates, where training data arising from multiple model forms and resolutions are integrated to predict high-fidelity response quantities of interest at lower cost. We focus on the context of quantum chemistry and the integration of information from multiple levels of theory. Important foundations include the use of symmetry function-based atomic energy vector constructions as feature vectors for representing structures across families of molecules and single-fidelity neural network training capabilities that learn the relationships needed to map feature vectors to potential energy predictions. These foundations are embedded within several multifidelity topologies that decompose the high-fidelity mapping into model-based components, including sequential formulations that admit a general nonlinear mapping across fidelities and discrepancy-based formulations that presume an additive decomposition. Methodologies are first explored and demonstrated on a pair of simple analytical test problems and then deployed for potential energy prediction for C5H5 using B2PLYP-D3/6-311++G(d,p) for high-fidelity simulation data and Hartree-Fock 6-31G for low-fidelity data. For the common case of limited access to high-fidelity data, our computational results demonstrate that multifidelity neural network potential energy surface constructions achieve roughly an order of magnitude improvement, either in terms of test error reduction for equivalent total simulation cost or reduction in total cost for equivalent error.

Focused ion beam implantation is ideally suited for placing defect centers in wide bandgap semiconductors with nanometer spatial resolution. However, the fact that only a few percent of implanted defects can be activated to become efficient single photon emitters prevents this powerful capability to reach its full potential in photonic/electronic integration of quantum defects. Here an industry adaptive scalable technique is demonstrated to deterministically create single defects in commercial grade silicon carbide by performing repeated low ion number implantation and in situ photoluminescence evaluation after each round of implantation. An array of 9 single defects in 13 targeted locations is successfully created—a ≈70% yield which is more than an order of magnitude higher than achieved in a typical single pass ion implantation. The remaining emitters exhibit non-classical photon emission statistics corresponding to the existence of at most two emitters. This approach can be further integrated with other advanced techniques such as in situ annealing and cryogenic operations to extend to other material platforms for various quantum information technologies.

Mineral dust can indirectly impact climate by nucleation of atmospheric solids, for example, by heterogeneously nucleating ice in mixed-phase clouds or by impacting the phase of aerosols and clouds through contact nucleation. The effectiveness toward nucleation of individual components of mineral dust requires further study. Here, the nucleation behavior of metal oxide nanoparticle components of atmospheric mineral dust is investigated. A long-working-distance optical trap is used to study contact and immersion nucleation of ammonium sulfate by transition-metal oxides, and an environmental chamber is used to probe depositional ice nucleation on metal oxide particles. Previous theory dictates that ice nucleation and heterogeneous nucleation of atmospheric salts can be impacted by several factors including morphology, lattice match, and surface area. Here, we observe a correlation between the cationic oxidation states of the metal oxide heterogeneous nuclei and their effectiveness in causing nucleation in both contact efflorescence mode and depositional freezing mode. In contrast to the activity of contact efflorescence, the same metal oxide particles did not cause a significant increase in efflorescence relative humidity when immersed in the droplet. These experiments suggest that metal speciation, possibly as a result of cationic charge sites, may play a role in the effectiveness of nucleation that is initiated at particle surfaces.

Lithium dendrite growth hinders the use of lithium metal anodes in commercial batteries. We present a 3D model to study the mechanical and electrochemical mechanisms that drive microscale plating. With this model, we investigate electrochemical response across a lithium protrusion characteristic of rough anode surfaces, representing the separator as a porous polymer in non-conformal contact with a lithium anode. The impact of pressure on separator morphology and electrochemical response is of particular interest, as external pressure can improve cell performance. We explore the relationships between plating propensity, stack pressure, and material properties. External pressure suppresses lithium plating due to interfacial stress and separator pore closure, leading to inhomogeneous plating rates. For moderate pressures, dendrite growth is completely suppressed, as plating will occur in the electrolyte-filled gaps between anode and separator. In fast-charging conditions and systems with low electrolyte diffusivities, the benefits of pressure are overridden by ion transport limitations.

Lithium metal solid-state batteries (LiSSBs) present new challenges in the measurement of material, component, and cell mechanical behaviors and in the measurement and theory of fundamental mechanical-electrochemical (thermodynamics, transport, and kinetics) couplings. Here, we classify the major mechanical and electrochemical-mechanical (ECM) studies underway and provide an overview of major mechanical testing platforms. We emphasize key distinctions among testing platforms, including tip- vs. platen-based sample compression, surface- vs. volume-based analysis, ease of integration with a vacuum or inert atmosphere environment, the ability to control and measure force/displacement over long periods of time, ranges of force and contact area, and others. Among the techniques we review, nanoindentation platforms offer some unique benefits associated with being able to use both tip-based nanoindentation techniques as well as platen-based compression over areas approaching 1 mm2. Sample design is also important: while most efforts are particle-based (i.e., using particles of solid electrolyte and cathode-active materials and densifying them using sintering or pressure), the resulting electrochemical response is from the overall collection of particles present. In contrast, thin-film (<1 μm) solid-state battery materials (e.g., Li, LiPON, LCO) provide well defined and uniform structures well suited for fundamental electrochemical-mechanical studies and offer an important opportunity to drive underlying scientific advances in LiSSB and other areas. We believe there are exciting opportunities to advance the measurement of both mechanical properties and electrochemical-mechanical couplings through the careful and novel co-design of test structures and experimental approaches for LiSSB materials, components, and cells.