Publications

Li metal anodes are enticing for batteries due to high theoretical charge storage capacity, but commercialization is plagued by dendritic Li growth and short circuits when cycled at high currents. Applied pressure has been suggested to improve morphology, and therefore performance. We hypothesized that increasing pressure would suppress dendritic growth at high currents. To test this hypothesis, here, we extensively use cryogenic scanning electron microscopy to show that varying the applied pressure from 0.01 to 1 MPa has little impact on Li morphology after one deposition. We show that pressure improves Li density and preserves Li inventory after 50 cycles. However, contrary to our hypothesis, pressure exacerbates dendritic growth through the separator, promoting short circuits. Therefore, we suspect Li inventory is better preserved in cells cycled at high pressure only because the shorts carry a larger portion of the current, with less being carried by electrochemical reactions that slowly consume Li inventory.

The canonical beam splitter - a fundamental building block of quantum optical systems - is a reciprocal element. It operates on forward- and backward-propagating modes in the same way, regardless of direction. The concept of nonreciprocal quantum photonic operations, by contrast, could be used to transform quantum states in a momentum- and direction-selective fashion. Here we demonstrate the basis for such a nonreciprocal transformation in the frequency domain through intermodal Bragg scattering four-wave mixing (BSFWM). Since the total number of idler and signal photons is conserved, the process can preserve coherence of quantum optical states, functioning as a nonreciprocal frequency beam splitter. We explore the origin of this nonreciprocity and find that the phase-matching requirements of intermodal BSFWM produce an enormous asymmetry (76×) in the conversion bandwidths for forward and backward configurations, yielding ∼25 dB of nonreciprocal contrast over several hundred GHz. We also outline how the demonstrated efficiencies (∼10-4) may be scaled to near-unity values with readily accessible powers and pumping configurations for applications in integrated quantum photonics.

Natural and anthropogenic events may create low frequency sound waves, or infrasound, that can travel for vast distances in planetary atmospheres. They permit the remote monitoring of geophysical activity over local to global scales. Most studies have utilized ground-based recorders, but it is possible to deploy acoustic sensors to altitudes of over 50 km. Such elevated platforms can capture sounds that their surface analogs cannot access. High altitude balloons and low altitude aerostats are filling this observation gap, but key environments remain out of reach of both of these. Recent work by den Ouden, Smets et al. (2021) addressed this with a new instrumentation platform—a large seabird flying just above the ocean's surface. Their work demonstrates that, infrasound sensing using heavier-than-air platforms in windy environments is possible, which has implications both terrestrially (e.g., extending sensor networks over the oceans) and extraterrestrially (proposed or planned missions to Venus and Titan).

We present a minimally invasive method for forward propagation of material property uncertainty to full-field quantities of interest in solid dynamics. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known a priori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the evolution equations based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the driving term of the evolution equations. The implementation is model independent, storage and communication efficient, and straightforward. We demonstrate its efficiency, accuracy, scaling with dimension of the parameter space, and convergence in distribution with two problems: a quasi-one-dimensional bar impact, and a two material notched plate impact. For the bar impact problem, we provide an analytical solution to PDF of the solution fields for method validation. With the notched plate problem, we also demonstrate good parallel efficiency and scaling of the method.

Despite widespread use of factor fixing in environmental modeling, its effect on model predictions has received little attention and is instead commonly presumed to be negligible. We propose a proof-of-concept adaptive method for systematically investigating the impact of factor fixing. The method uses Global Sensitivity Analysis methods to identify groups of sensitive parameters, then quantifies which groups can be safely fixed at nominal values without exceeding a maximum acceptable error, demonstrated using the 21-dimensional Sobol’ G-function. Furthermore, three error measures are considered for quantities of interest, namely Relative Mean Absolute Error, Pearson Product-Moment Correlation and Relative Variance. Results demonstrate that factor fixing may cause large errors in the model results unexpectedly, when preliminary analysis suggests otherwise, and that the default value selected affects the number of factors to fix. To improve the applicability and methodological development of factor fixing, a new research agenda encompassing five opportunities is discussed for further attention.

Motivation: Crucial aspect of mechanical design is joining methodology of parts. Ability to analyze joint and fasteners in system for structural integrity is fundamental. Different modeling representations of fasteners include spring, beam, and solid elements. Various methods compared for linear system to decide method appropriate for design study. New method for modeling fastener joint is explored from full system perspective. Analysis results match well with published experimental data for new method.

Ultrawide-bandgap (UWBG) semiconductor technology is presently going through a renaissance exemplified by advances in material-level understanding, extensions of known concepts to new materials, novel device concepts, and new applications. This focus issue presents a timely selection of papers spanning the current state of the art in UWBG materials and applications, including both experimental results and theoretical developments. It covers broad research subtopics on UWBG bulk crystals and substrate technologies, UWBG defect science and doping, UWBG epitaxy, UWBG electronic and optoelectronic properties, and UWBG power devices and emitters. In this overview article, we consolidate the fundamentals and background of key UWBG semiconductors including aluminum gallium nitride alloys (AlxGa1–xN), boron nitride (BN), diamond, β-phase gallium oxide (β-Ga2O3), and a number of other UWBG binary and ternary oxides. Graphical Abstract: [Figure not available: see fulltext.]

This report documents the analysis conducted by Sandia National Laboratories on the effect of various alterations to published methodologies to calculate Derived Response Levels for Environmental Protection Agency Drinking Water Protective Action Guides. Specifically, this study sought to assess and provide recommendations on calculation of the Derived Response Level accounting for decay during the consumption period, assess the impact of decay on laboratory Minimal Detectable Concentration in water samples as compared to the Derived Response Level, make a recommendation on the calculation of Derived Response Level consistent with existing Public Protection Methods, and make a recommendation on the use of six age groups versus eight age groups based on available dose coefficients for calculation of the Derived Response Level. The authors analyzed these various factors using nominal radionuclide mixes from four scenarios and compared calculation of the Derived Response Level accounting for decay and no decay and then compared those results to the laboratory Minimal Detectable Concentrations. The authors concluded that decay should be included in the calculation of the Derived Response Level, existing Public Protection Methods should be employed to calculate the Derived Response Level, six age groups should be used versus eight, and the use of both decay and Public Protection Methods result in little to no concern for water samples meeting the Minimum Detectable Concentration requirements. The results of this study may be used in further developing and implementing a method for the Environmental Protection Agency Water Derived Response Level calculation in the Federal Radiological Monitoring and Assessment Center Assessment Manual.

The construction of deep geological repositories (DGR) in salt formations requires penetrating through naturally sealing geosphere layers. While the emplaced nuclear waste is primarily protected by the containment-providing rock zone (CRZ), technical barriers are required, for example during handling. For closure geotechnical barriers seal the repository along the accesses against water or solutions from outside and the possible emission paths for radionuclides contained inside. As these barriers must ensure maintenance-free function on a long-term basis, they typically comprise a set of specialized elements with diversified functions that may be used redundantly. The effects of the individual elements are coordinated so that they are collectively referred to as the Engineered Barrier System (EBS).

Targeting host factors for anti-viral development offers several potential advantages over traditional countermeasures that include broad-spectrum activity and prevention of resistance. Characterization of host factors in animal models provides strong evidence of their involvement in disease pathogenesis, but the feasibility of performing high-throughput in vivo analyses on lists of genes is problematic. To begin addressing the challenges of screening candidate host factors in vivo, we combined advances in CRISPR-Cas9 genome editing with an immunocompromised mouse model used to study highly pathogenic viruses. Transgenic mice harboring a constitutively expressed Cas9 allele (Cas9tg/tg) with or without knockout of type I interferon receptors served to optimize in vivo delivery of CRISPR single-guide RNA (sgRNA) using Invivofectamine 3.0, a simple and easy-to-use lipid nanoparticle reagent. Invivofectamine 3.0-mediated liver-specific editing to remove activity of the critical Ebola virus host factor Niemann-Pick disease type C1 in an average of 74% of liver cells protected immunocompromised Cas9tg/tg mice from lethal surrogate Ebola virus infection. We envision that immunocompromised Cas9tg/tg mice combined with straightforward sgRNA in vivo delivery will enable efficient host factor loss-of-function screening in the liver and other organs to rapidly study their effects on viral pathogenesis and help initiate development of broad-spectrum, host-directed therapies against emerging pathogens.

A stable solid electrolyte interphase (SEI) layer is key to high performing lithium ion and lithium metal batteries for metrics such as calendar and cycle life. The SEI must be mechanically robust to withstand large volumetric changes in anode materials such as lithium and silicon, so understanding the mechanical properties and behavior of the SEI is essential for the rational design of artificial SEI and anode form factors. The mechanical properties and mechanical failure of the SEI are challenging to study, because the SEI is thin at only ∼10-200 nm thick and is air sensitive. Furthermore, the SEI changes as a function of electrode material, electrolyte and additives, temperature, potential, and formation protocols. A variety of in situ and ex situ techniques have been used to study the mechanics of the SEI on a variety of lithium ion battery anode candidates; however, there has not been a succinct review of the findings thus far. Because of the difficulty of isolating the true SEI and its mechanical properties, there have been a limited number of studies that can fully de-convolute the SEI from the anode it forms on. A review of past research will be helpful for culminating current knowledge and helping to inspire new innovations to better quantify and understand the mechanical behavior of the SEI. This review will summarize the different experimental and theoretical techniques used to study the mechanics of SEI on common lithium battery anodes and their strengths and weaknesses.

This report documents the feasibility phase of the Critical Experiment Design (CED) conducted as part of integral experiment request (IER) 523. The purpose of IER-523 is to explore the effects of using 35 weight percent enriched uranium dioxide-beryllium oxide (UO₂-BeO) material on critical configurations using the Seven Percent Critical Experiment (7uPCX) at Sandia National Laboratories (Sandia). Preliminary experiment design concepts, neutronic analysis results, and proposed paths for continuing the CED process are presented.

Understanding how atoms interact with hot dense matter is essential for astrophysical and laboratory plasmas. Interactions in high-density plasmas broaden spectral lines, providing a rare window into interactions that govern, for example, radiation transport in stars. However, up to now, spectral line-shape theories employed at least one of three common approximations: second-order Taylor treatment of broadening operator, dipole-only interactions between atom and plasma, and classical treatment of perturbing electrons. In this Letter, we remove all three approximations simultaneously for the first time and test the importance for two applications: neutral hydrogen and highly ionized magnesium and oxygen. We found 15%-50% change in the spectral line widths, which are sufficient to impact applications including white-dwarf mass determination, stellar-opacity research, and laboratory plasma diagnostics.

The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.

Many innate immune receptors function collaboratively to detect and elicit immune responses to pathogens, but the physical mechanisms that govern the interaction and signaling crosstalk between the receptors are unclear. In this study, we report that the signaling crosstalk between Fc gamma receptor (FcγR) and Toll-like receptor (TLR)2/1 can be overall synergistic or inhibitory depending on the spatial proximity between the receptor pair on phagosome membranes. Using a geometric manipulation strategy, we physically altered the spatial distribution of FcγR and TLR2 on single phagosomes. We demonstrate that the signaling synergy between FcγR and TLR2/1 depends on the proximity of the receptors and decreases as spatial separation between them increases. However, the inhibitory effect from FcγRIIb on TLR2-dependent signaling is always present and independent of receptor proximity. The overall cell responses are an integration from these two mechanisms. This study presents quantitative evidence that the nanoscale proximity between FcγR and TLR2 functions as a key regulatory mechanism in their signaling crosstalk.

Big Data in the area of Remote Sensing has been growing rapidly. Remote sensors are used in surveillance, security, traffic, environmental monitoring, and autonomous sensing. Real-time detection of small moving targets using a remote sensor is an ongoing, challenging problem. Since the object is located far away from the sensor, the object often appears too small. The object’s signal-to-noise-ratio (SNR) is often very low. Occurrences such as camera motion, moving backgrounds (e.g., rustling leaves), low contrast and resolution of foreground objects makes it difficult to segment out the targeted moving objects of interest. Due to the limited appearance of the target, it is tough to obtain the target’s characteristics such as its shape and texture. Without these characteristics, filtering out false detections can be a difficult task. Detecting these targets, would often require the detector to operate under a low detection threshold. However, lowering the detection threshold could lead to an increase of false alarms. In this paper, the author will introduce a new method that improves the probability to detect low SNR objects, while decreasing the number of false alarms as compared to using the traditional baseline detection technique.

Layered boron compounds have attracted significant interest in applications from energy storage to electronic materials to device applications, owing in part to a diversity of surface properties tied to specific arrangements of boron atoms. Here we report the energy landscape for surface atomic configurations of MgB2 by combining first-principles calculations, global optimization, material synthesis and characterization. We demonstrate that contrary to previous assumptions, multiple disordered reconstructions are thermodynamically preferred and kinetically accessible within exposed B surfaces in MgB2 and other layered metal diborides at low boron chemical potentials. Such a dynamic environment and intrinsic disordering of the B surface atoms present new opportunities to realize a diverse set of 2D boron structures. We validated the predicted surface disorder by characterizing exfoliated boron-terminated MgB2 nanosheets. We further discuss application-relevant implications, with a particular view towards understanding the impact of boron surface heterogeneity on hydrogen storage performance.

The addition of active, nonlinear, and nonreciprocal functionalities to passive piezoelectric acoustic wave technologies could enable all-acoustic and therefore ultra-compact radiofrequency signal processors. Toward this goal, we present a heterogeneously integrated acoustoelectric material platform consisting of a 50 nm indium gallium arsenide epitaxial semiconductor film in direct contact with a 41° YX lithium niobate piezoelectric substrate. We then demonstrate three of the main components of an all-acoustic radiofrequency signal processor: passive delay line filters, amplifiers, and circulators. Heterogeneous integration allows for simultaneous, independent optimization of the piezoelectric-acoustic and electronic properties, leading to the highest performing surface acoustic wave amplifiers ever developed in terms of gain per unit length and DC power dissipation, as well as the first-ever demonstrated acoustoelectric circulator with an isolation of 46 dB with a pulsed DC bias. Finally, we describe how the remaining components of an all-acoustic radiofrequency signal processor are an extension of this work.

Organophosphorus hydrolase (OPH) is a metalloenzyme that can hydrolyze organophosphorus agents resulting in products that are generally of reduced toxicity. The best OPH substrate found to date is diethyl p-nitrophenyl phosphate (paraoxon). Most structural and kinetic studies assume that the binding orientation of paraoxon is identical to that of diethyl 4-methylbenzylphosphonate, which is the only substrate analog co-crystallized with OPH. In the current work, we used a combined docking and molecular dynamics (MD) approach to predict the likely binding mode of paraoxon. Then, we used the predicted binding mode to run MD simulations on the wild type (WT) OPH complexed with paraoxon, and OPH mutants complexed with paraoxon. Additionally, we identified three hot-spot residues (D253, H254, and I255) involved in the stability of the OPH active site. We then experimentally assayed single and double mutants involving these residues for paraoxon binding affinity. The binding free energy calculations and the experimental kinetics of the reactions between each OPH mutant and paraoxon show that mutated forms D253E, D253E-H254R, and D253E-I255G exhibit enhanced substrate binding affinity over WT OPH. Interestingly, our experimental results show that the substrate binding affinity of the double mutant D253E-H254R increased by 19-fold compared to WT OPH.

The Message Passing Interface (MPI) remains the dominant programming model for scientific applications running on today's high-performance computing (HPC) systems. This dominance stems from MPI's powerful semantics for inter-process communication that has enabled scientists to write applications for simulating important physical phenomena. MPI does not, however, specify how messages and synchronization should be carried out. Those details are typically dependent on low-level architecture details and the message characteristics of the application. Therefore, analyzing an application's MPI resource usage is critical to tuning MPI's performance on a particular platform. The result of this analysis is typically a discussion of the mean message sizes, queue search lengths and message arrival times for a workload or set of workloads. While a discussion of the arithmetic mean in MPI resource usage might be the most intuitive summary statistic, it is not always the most accurate in terms of representing the underlying data. In this paper, we analyze MPI resource usage for a number of key MPI workloads using an existing MPI trace collector and discrete-event simulator. Our analysis demonstrates that the average, while easy and efficient to calculate, is a useful metric for characterizing latency and bandwidth measurements, but may not be a good representation of application message sizes, match list search depths, or MPI inter-operation times. Additionally, we show that the median and mode are superior choices in many cases. We also observe that the arithmetic mean is not the best representation of central tendency for data that are drawn from distributions that are multi-modal or have heavy tails. The results and analysis of our work provide valuable guidance on how we, as a community, should discuss and analyze MPI resource usage data for scientific applications.

A resonant cavity undergoes three distinct behaviors with increasing frequency: 1) fundamental modes, localized in frequency with well defined modal distribution; 2) undermoded region, where modes are still separated, but are sufficiently perturbed by small imperfections that their spectral positions (and distributions) are statistical in nature; and 3) overmoded region, where modes overlap, field distributions follow stochastic distributions, and the slot acts as if in free space. Understanding the penetration through slots in the overmoded region is of great interest, and is the focus of this article. Since full-wave solvers may not be able to provide a timely answer for very high frequencies due to a lack of memory and/or computation resources, we develop bounding methods to estimate worst-case average and maximum fields within the cavity. After discussing the bounding formulation, we compare its results to full-wave simulations at the first, second, and third resonance supported by the slot in the case of a cylindrical cavity. Note that the bounding formulation indicates that results are nearly independent of cavity shape: only the cavity volume, frequency, and cavity quality factor affect the overmoded region, making this formulation a powerful tool to assess electromagnetic interference and electromagnetic compatibility effects within cavities.

Point defects in SiC are an attractive platform for quantum information and sensing applications because they provide relatively long spin coherence times, optical spin initialization, and spin-dependent fluorescence readout in a fabrication-friendly semiconductor. The ability to precisely place these defects at the optimal location in a host material with nano-scale accuracy is desirable for integration of these quantum systems with traditional electronic and photonic structures. Here, we demonstrate the precise spatial patterning of arrays of silicon vacancy (VSi) emitters in an epitaxial 4H-SiC (0001) layer through mask-less focused ion beam implantation of Li+. We characterize these arrays with high-resolution scanning confocal fluorescence microscopy on the Si-face, observing sharp emission lines primarily coming from the V1 ′ zero-phonon line (ZPL). The implantation dose is varied over 3 orders of magnitude, leading to VSi densities from a few per implantation spot to thousands per spot, with a linear dependence between ZPL emission and implantation dose. Optically-detected magnetic resonance (ODMR) is also performed, confirming the presence of V2 VSi. Our investigation reveals scalable and reproducible defect generation.

The isotopic composition of water in Earth’s oceans is challenging to recreate using a plausible mixture of known extraterrestrial sources such as asteroids—an additional isotopically light reservoir is required. The Sun’s solar wind could provide an answer to balance Earth’s water budget. We used atom probe tomography to directly observe an average ~1 mol% enrichment in water and hydroxyls in the solar-wind-irradiated rim of an olivine grain from the S-type asteroid Itokawa. We also experimentally confirm that H+ irradiation of silicate mineral surfaces produces water molecules. These results suggest that the Itokawa regolith could contain ~20 l m−3 of solar-wind-derived water and that such water reservoirs are probably ubiquitous on airless worlds throughout our Galaxy. The production of this isotopically light water reservoir by solar wind implantation into fine-grained silicates may have been a particularly important process in the early Solar System, potentially providing a means to recreate Earth’s current water isotope ratios.

The study of microstructure evolution in additive manufacturing of metals would be aided by knowing the thermal history. Since temperature measurements beneath the surface are difficult, estimates are obtained from computational thermo-mechanical models calibrated against traces left in the sample revealed after etching, such as the trace of the melt pool boundary. Here we examine the question of how reliable thermal histories computed from a model that reproduces the melt pool trace are. To this end, we perform experiments in which one of two different laser beams moves with constant velocity and power over a substrate of 17-4PH SS or Ti-6Al-4V, with low enough power to avoid generating a keyhole. We find that thermal histories appear to be reliably computed provided that (a) the power density distribution of the laser beam over the substrate is well characterized, and (b) convective heat transport effects are accounted for. Poor control of the laser beam leads to potentially multiple three-dimensional melt pool shapes compatible with the melt pool trace, and therefore to multiple potential thermal histories. Ignoring convective effects leads to results that are inconsistent with experiments, even for the mild melt pools here.

The solid-state joining of oxide-dispersion-strengthened (ODS) austenitic steels was achieved using a pulsed electric current joining (PECJ) process. Microstructures of the austenitic grain structures and oxide dispersions in the joint areas were characterized using electron microscopy. Negligible grain growth was observed in austenitic grain structures, while slight coarsening of oxide dispersions occurred at a short holding time. The mechanisms of the PECJ process may involve three steps that occur simultaneously, including the sintering of mechanical alloying powders in the bonding layer, formation of oxide dispersions, and bonding of the mechanical alloying powders with the base alloy. The high hardness and irradiation resistance of ODS alloys were retained in the joint areas. This research revealed the fundamental mechanisms during the PECJ process, which is beneficial for its potential applications during the advanced manufacturing of ODS alloys.