Publications

Various laboratory-focused tools and methodologies for completing a safety risk assessment have been published, yet few similar resources to address chemical security exist.Herein, we describe a chemical security risk assessment case study at a university in a developing country.In this case study, we demonstrate a chemical security risk assessment for a university chemistry department, using an original inventory of 645 entries which was condensed to 295 chemicals after removing duplicates and erroneous entries.We then prioritized to highlight 83 chemicals of interest based on hazardous or dual-use properties that could lead to unacceptable consequences.We further refined to a list of 34 high-risk chemicals that required action, 48 chemicals that may need further justification and consideration for additional protection, and 1 chemical that did not need further consideration for additional protection.

Time-dependent visible and near-ultraviolet emission spectroscopy was used to track the presence of excited and ionized states induced as an electron beam transited a cavity filled with air at low pressures. An electron beam was produced in vacuum using a Febetron pulsed-power generator modified to produce a peak voltage of 100 kV, a peak current of 4.5 kA, and a pulse width of 100 ns. The beam electrons then passed through thin anode and pressure foils and transited a cavity filled with dry air at either 0.1 Torr or 1 Torr, exciting and ionizing the air along the way. The spectral measurements were combined with information from electrical, interferometric, and imaging diagnostics to reveal different ionization dynamics for the two cases under study.

A well-posed physics-based compact model for a three-terminal silicon–oxide–nitride–oxide–silicon (SONOS) synaptic circuit element is presented for use by neuromorphic circuit/system engineers. Based on technology computer aided design (TCAD) simulations of a SONOS device, the model contains a nonvolatile memristor with the state variable QM representing the memristor charge under the gate of the three-terminal element. By incorporating the exponential dependence of the memristance on QM and the applied bias V for the gate, the compact model agrees quantitatively with the results from TCAD simulations as well as experimental measurements for the drain current. The compact model is implemented through VerilogA in the circuit simulation package Cadence Spectre and reproduces the experimental training behavior for the source–drain conductance of a SONOS device after applying writing pulses ranging from –12 V to +11 V, with an accuracy higher than 90%.

We report pure-rotational N2-N2, N2-air, and O2-air S-branch linewidths for temperatures of 80-200 K by measuring the time-dependent decay of rotational Raman coherences in an isentropic free-jet expansion from a sonic nozzle. We recorded pure-rotational hybrid femtosecond/picosecond coherent anti-Stokes Raman scattering (fs/ps CARS) spectra along the axial centerline of the underexpanded jet, within the barrel shock region upstream of the Mach disk. The dephasing of the pure-rotational Raman coherence was monitored using probe-time-delay scans at different axial positions in the jet, corresponding to varying local temperatures and pressures. The local temperature was obtained by fitting CARS spectra acquired at zero probe time delay, where the impact of collisions was minimal. The measured decay of each available Raman transition was fit to a dephasing constant and corrected for the local pressure, which was obtained from the CARS-measured static temperature and thermodynamic relationships for isentropic expansion from the known stagnation state. Nitrogen self-broadened transitions decayed more rapidly than those broadened in air for all temperatures, corresponding to higher Raman linewidths. In general, the measured S-branch linewidths deviated significantly in absolute and relative magnitudes from those predicted by extrapolating the modified exponential gap model to low temperatures. The temperature dependence of the Raman linewidth for each measured rotational state of nitrogen (J ≤ 10) and oxygen (N ≤ 11) was fit to a temperature-dependent power law over the measurable temperature domain (80-200 K) and extrapolated to both higher rotational states and room temperature. The measured and modeled low-temperature linewidth data provided here will aid low temperature gas-phase pressure measurements with fs/ps CARS.

Inverse prediction models have commonly been developed to handle scalar data from physical experiments. However, it is not uncommon for data to be collected in functional form. When data are collected in functional form, it must be aggregated to fit the form of traditional methods, which often results in a loss of information. For expensive experiments, this loss of information can be costly. In this study, we introduce the functional inverse prediction (FIP) framework, a general approach which uses the full information in functional response data to provide inverse predictions with probabilistic prediction uncertainties obtained with the bootstrap. The FIP framework is a general methodology that can be modified by practitioners to accommodate many different applications and types of data. We demonstrate the framework, highlighting points of flexibility, with a simulation example and applications to weather data and to nuclear forensics. Results show how functional models can improve the accuracy and precision of predictions.

Pulsed-power generators can produce well-controlled continuous ramp compression of condensed matter for high-pressure equation-of-state studies using the magnetic loading technique. X-ray diffraction (XRD) data from dynamically compressed samples provide direct measurements of the elastic compression of the crystal lattice, onset of plastic flow, strength–strain rate dependence, structural phase transitions, and density of crystal defects, such as dislocations. Here, we present a cost-effective, compact, pulsed x-ray source for XRD measurements on pulsed-power-driven ramp-loaded samples. This combination of magnetically driven ramp compression of materials with a single, short-pulse XRD diagnostic will be a powerful capability for the dynamic materials’ community to investigate in situ dynamic phase transitions critical to equation of states. Finally, we present results using this new diagnostic to evaluate lattice compression in Zr and Al and to capture signatures of phase transitions in CdS.

We have developed a strategy for distinguishing between small-angle X-ray scattering (SAXS) from gas-phase species and newly formed nanoparticles in mixed gas- and particle-phase reacting flows. This methodology explicitly accounts for temperature-dependent scattering from gases. We measured SAXS in situ in a sooting linear laminar partially premixed co-flow ethylene/air diffusion flame. The scattering signal demonstrates a downward curvature as a function of the momentum transfer (q) at q values of 0.2-0.57 Å-1. The q-dependent curvature is consistent with the Debye equation and the independent-atom model for gas-phase scattering. This behavior can also be modeled using the Guinier approximation and could be characterized as a Guinier knee for gas-phase scattering. The Guinier functional form can be fit to the scattering signal in this q range without a priori knowledge of the gas-phase composition, enabling estimation of the gas-phase contribution to the scattering signal while accounting for changes in the gas-phase composition and temperature. We coupled the SAXS measurements with in situ temperature measurements using coherent anti-Stokes Raman spectroscopy. This approach to characterizing the gas-phase SAXS signal provides a physical basis for distinguishing among the contributions to the scattering signal from the instrument function, flame gases, and nanoparticles. The results are particularly important for the analysis of the SAXS signal in the q range associated with particles in the size range of 1-6 nm.

We present an in-depth study of metal-insulator interfaces within granular metal (GM) films and correlate their interfacial interactions with structural and electrical transport properties. Nominally 100 nm thick GM films of Co and Mo dispersed within yttria-stabilized zirconia (YSZ), with volumetric metal fractions (φ) from 0.2-0.8, were grown by radio frequency co-sputtering from individual metal and YSZ targets. Scanning transmission electron microscopy and DC transport measurements find that the resulting metal islands are well-defined with 1.7-2.6 nm average diameters and percolation thresholds between φ = 0.4-0.5. The room temperature conductivities for the φ = 0.2 samples are several orders of magnitude larger than previously-reported for GMs. X-ray photoemission spectroscopy indicates both oxygen vacancy formation within the YSZ and band-bending at metal-insulator interfaces. The higher-than-predicted conductivity is largely attributed to these interface interactions. In agreement with recent theory, interactions that reduce the change in conductivity across the metal-insulator interface are seen to prevent sharp conductivity drops when the metal concentration decreases below the percolation threshold. These interface interactions help interpret the broad range of conductivities reported throughout the literature and can be used to tune the conductivities of future GMs.

Here we present a new analysis methodology that allows for the self-consistent integration of multiple diagnostics including nuclear measurements, x-ray imaging, and x-ray power detectors to determine the primary stagnation parameters, such as temperature, pressure, stagnation volume, and mix fraction in magnetized liner inertial fusion (MagLIF) experiments. The analysis uses a simplified model of the stagnation plasma in conjunction with a Bayesian inference framework to determine the most probable configuration that describes the experimental observations while simultaneously revealing the principal uncertainties in the analysis. We validate the approach by using a range of tests including analytic and three-dimensional MHD models. An ensemble of MagLIF experiments is analyzed, and the generalized Lawson criterion χ is estimated for all experiments.

Aqueous zinc-metal batteries are plagued by poor Zn reversibility owing to zinc dendrite and layered double hydroxide (LDH) formation. Here, we introduce a novel additive—N,N-dimethylformamidium trifluoromethanesulfonate (DOTf)—in a low-cost aqueous electrolyte that can very effectively address these issues. The initial water-assisted dissociation of DOTf into triflic superacid creates a robust nanostructured solid-electrolyte interface (SEI)—revealed by operando spectroscopy and cryomicroscopy—which excludes water and enables dense Zn deposition. We demonstrate excellent Zn plating/stripping in a Zn||Cu asymmetric cell for more than 3,500 cycles. Furthermore, near 100% CE is realized at a combined high current density of 4 mA cm−2 and an areal capacity of 4 mAh cm−2 over long-term cycling. Zn||Zn0.25V2O5·nH2O full cells retain ∼83% of their capacity after 1,000 cycles with mass-limited Zn anodes. By restricting the depth of discharge, the cathodes exhibit less proton intercalation and LDH formation with an extended lifetime of 2,000 cycles.

Understanding the mechanisms of charge transport in batteries is important for the rational design of new electrolyte formulations. Persistent questions about ion transport mechanisms in battery electrolytes are often framed in terms of vehicular diffusion by persistent ion-solvent complexes versus structural diffusion through the breaking and reformation of ion-solvent contacts, i.e., solvent exchange events. Ultrafast two-dimensional (2D) IR spectroscopy can probe exchange processes directly via the evolution of the cross-peaks on picosecond time scales. However, vibrational energy transfer in the absence of solvent exchange gives rise to the same spectral signatures, hiding the desired processes. We employ 2D IR on solvent resonances of a mixture of acetonitrile isotopologues to differentiate chemical exchange and energy-transfer dynamics in a comprehensive series of Li+, Mg2+, Zn2+, Ca2+, and Ba2+ bis(trifluoromethylsulfonyl)imide electrolytes from the dilute to the superconcentrated regime. No exchange phenomena occur within at least 100 ps, regardless of the ion identity, salt concentration, and presence of water. All of the observed spectral dynamics originate from the intermolecular energy transfer. These results place the lower experimental boundary on the ion-solvent residence times to several hundred picoseconds, much slower than previously suggested. With the help of MD simulations and conductivity measurements on the Li+ and Zn2+ systems, we discuss these results as a continuum of vehicular and structural modalities that vary with concentration and emphasize the importance of collective electrolyte motions to ion transport. These results hold broadly applicable to many battery-relevant ions and solvents.

Constitutive model parameterizations for the General Plastics EF4003 low density 3 pound per cubic foot are needed for design and qualification purposes in normal and abnormal mechanical simulations. The material is expected to be deformed in two ways: first during preloading, and second under impact conditions of the system (transient dynamic). All analyses are to be performed at room temperature. The goal is to provide the analysis community a robust constitutive model parameterization to represent the compression behavior of the EF4003 foam from small deformations up to massive compressive deformations when the foam is densifying. It is worth noting the EF4003 exhibits anisotropy in its stress-strain behavior between the rise and transverse directions (See figure 2.8c-d) as well as plateau behavior that is very likely to cause material stability issues, due to the buckling transition, (and has historically done so) when using Sandia’s current workhorse models for flexible foams, Hyperfoam and Flex Foam. A Stability-informed Hyperfoam parameterization procedure is developed and executed to calibrate a hyperfoam model for the EF4003 room temperature, transversely loaded data. A rise orientation parameterization was not attempted due to localization in the experiments.

Over the past 50 years, the Renewable Energy Program at Sandia has advanced research in the field with a focus on three key goals; 1) reduce the cost, 2) improve resilience and reliability and, 3) decrease the regulatory burden of renewable energy. Sandia’s expertise, coupled with the Village of Questa’s expanding renewable energy portfolio, presents the opportunity to deploy the Labs’ deep science and engineering capabilities towards the energy goals of KCEC and the Village of Questa. Preliminary research efforts by Sandia technical staff has broadly identified early opportunities for further research, development, and demonstration in the emerging renewable energy segment of agrivoltaics. Agrivoltaics is an emerging and promising area of photovoltaics which entails land use considerations as well as concerns regarding landscape transformation, biodiversity, and ecosystem well-being. In recent years, agrivoltaics systems have been the subject of numerous studies due to their potential in the food-energy (and water) nexus. This document is a preliminary evaluation of the projects performance opportunities of agrivoltaics as a renewable energy technology strategy in the region of Questa, NM.

Using the thesis of W.R. Nolan (cite) as a guide, Cobalt Iron (CoFe) powders were reacted with 0.1 wt.% and 0.2 wt.% phosphoric acid in a 20:1 ratio of acetone to phosphoric acid. The powders were then dried at room temperature. The resulting phosphate coated CoFe was then mixed with 0.75 wt.% of the lubricant N,N' ethylene bis-stearamide (trade name: Acrawax C) and hot pressed to form a consolidated soft magnetic material referred to as CoFeP. With an avenue of synthesis for CoFeP determined, a proper amount of stock was synthesized for continuous “brick” production. While under current optimization, these 1x1 mm magnetic bricks will ultimately be placed and secured along the inside wall of each MK Magnetics transformer core by an appropriate CoFeP dispersed epoxy. As of now each brick has been produced though a pressing and annealing process via square 2x2 cm die. Before a brick is made a pressure calculation is run to ensure the dies maximum operating pressure is not exceeded. Figure 1. ensures the user’s safety by showing that the tons-on-ram required for a 2x2 cm square die to reach 760 MPa is below the point of die failure.

The nuclear spins of low-density implanted Ga atoms in Ge are interesting candidates for solid state-based qubits. To date, activation studies of implanted Ga in Ge have focused on high densities. Here, we extend activation studies into the low-density regime. We use spreading resistance profiling and secondary ion mass spectrometry to derive electrical activation of Ga ions implanted into Ge as a function of the rapid thermal anneal temperature and implant density. We show that for our implant conditions, the activation is best for anneal temperatures between 400 and 650 °C with a maximum activation of 69% at the highest fluence. Below 400 °C, remaining implant damage results in defects that act as superfluous carriers, and above 650 °C, surface roughening and loss of Ga ions are observed. The activation increased monotonically from 10% to 69% as the implant fluence increased from 6 × 10 10 to 6 × 10 12 cm-2. The results provide thermal anneal conditions to be used for initial studies of using low-density Ga atoms in Ge as nuclear spin qubits.

As concerns about the illicit movement of radioactive materials across international borders increase, so too has the need for increased protection of those borders both foreign and domestic. The challenge is not only to detect hidden radioactive materials, but also to distinguish them from legitimate radionuclides such as radio-pharmaceuticals that are often transported across borders and shipped throughout a country. With more than 600 U.S. border sites to protect, screening imported radioactive material requires a careful balance of high throughput and high search efficiency. However, these requirements are difficult to meet as rapid screening operations leave less time for radiation detectors to efficiently evaluate materials. In support of border security, Sandia developed organic glass scintillators.

Particle fragmentation influences thermochemical energy conversion processes in different ways and is of significance in energy generation technologies. Different reactive material formulations trigger varied thermal response in extreme environments such as high velocity impact. This study investigated optical thermal response of powder gun launched intermetallic (Al:Zr) and thermite (Al:MoO3) projectiles using pyrometry and thermography. Projectiles were launched at 1250 m/s into an air-filled chamber and impacted a steel witness plate to create a dust explosion. The pyrometer was configured to measure temperatures directly at the point of impact, while the thermographic system measured temperatures throughout the explosion chamber. Results show that impact temperatures ranged between 3500 and 4000 K, but that the dynamics of energy conversion were different for the intermetallic and thermite projectiles. The intermetallic exhibited secondary reactions due to fragmented debris impacting the walls of the chamber. The thermite exhibited greater gas generation, propelling the debris field, and producing a more stochastic response with faster spreading and dissipation of thermal energy. Unique features such as microexplosions within fragmented particles were also analyzed. While both reactive materials produce similar temperatures, their mechanisms of energy conversion and release are different, indicating the potential of these materials for different ballistic applications.

Easily measured metrics that could assign quantifiable values to coating batches for quality control have started to be developed. High-density is an attribute of quality films. Increased density results in harder, more wear resistant coatings in inert and humid environments. Denser films are more resistant to oxidation from aging, limiting the severity and depth of oxide into the coating. Future work includes using metrics for quality. The next step is to develop in house deposition capabilities to develop process-structure relationships.

A fast matrix-free formulation of the control volume finite element method is presented, requiring much less memory and computational work than previous efforts. The method is implemented and evaluated as a solver for low-Mach flow, including the evaluation of a preconditioning strategy for the pressure Poisson equation. The efficiency and scaling with polynomial order is evaluated on simple turbulent flows of interest, with appropriate solution quality metrics, and compared with a reference node-centered finite volume discretization. For a turbulent channel flow test, we show improvement in computational work for a given accuracy with the high-order scheme. The performance on a GPU accelerated platform is also investigated, with benefit shown for the matrix-free discretization.

This report is the revised (Revision 8) Task F specification for DECOVALEX-2023. Task F is a comparison of the models and methods used in deep geologic repository performance assessment. The task proposes to develop a reference case for a mined repository in a fractured crystalline host rock and a reference case for a mined repository in a salt formation. Teams may choose to participate in the comparison for either or both of the reference cases. For each reference case, a common set of conceptual models and parameters describing features, events, and processes that impact performance will be given, and teams will be responsible for determining how best to implement and couple the models. The comparison will be conducted in stages, beginning with a comparison of key outputs of individual process models, followed by a comparison of a single deterministic simulation of the full reference case, and moving on to uncertainty propagation and uncertainty and sensitivity analysis. This report provides background information, a summary of the proposed reference cases, and a staged plan for the analysis.

Precise temperature determination is a significant challenge in extreme environments of dynamic compression studies. How can radiance measurements taken in high-pressure shock experiments constrain temperature in a meaningful and physically consistent way? Experiments maintaining sample compression against a transparent window can be tailored to present a uniform measurement area with uncertain spectral emissivity. We compare several methods to analyze radiance collected at multiple wavelengths, applying statistical methods and physical principles to improve temperature inference. With proper radiance collection and analysis, dynamic temperature uncertainties become comparable to thermomechanical ambiguities of the emitting surface.

This paper summarizes the evolution of the Gen 3 Particle Pilot Plant (G3P3) receiver design with the goal of reducing heat losses and increasing thermal efficiencies. New features that were investigated included aperture covers and shrouds, active airflow, multistage catch-and-release devices (stairs), and optimization of receiver cavity geometry. Simulations and ground-based testing showed that a reduced receiver volume and aperture shroud could reduce advective heat losses by ∼40 - 50%, and stairs could increase opacity and reduce backwall temperatures. The reduced volume receiver and stairs were selected for on-sun testing, and receiver efficiencies up to 80 - 90% were achieved in the current test campaign. The receiver thermal efficiency generally increased as a function of incident power and particle mass flow rates. In addition, particle outlet temperatures were maintained to within ±10 °C of a prescribed setpoint temperature up to ∼700 °C using a PID controller that adjusted the particle mass flow rate into the receiver in response to the measured particle outlet temperatures.

This paper summarizes the evolution of the Gen 3 Particle Pilot Plant (G3P3) receiver design with the goal of reducing heat losses and increasing thermal efficiencies. New features that were investigated included aperture covers and shrouds, active airflow, multistage catch-and-release devices (stairs), and optimization of receiver cavity geometry. Simulations and ground-based testing showed that a reduced receiver volume and aperture shroud could reduce advective heat losses by ∼40 - 50%, and stairs could increase opacity and reduce backwall temperatures. The reduced volume receiver and stairs were selected for on-sun testing, and receiver efficiencies up to 80 - 90% were achieved in the current test campaign. The receiver thermal efficiency generally increased as a function of incident power and particle mass flow rates. In addition, particle outlet temperatures were maintained to within ±10 °C of a prescribed setpoint temperature up to ∼700 °C using a PID controller that adjusted the particle mass flow rate into the receiver in response to the measured particle outlet temperatures.

Design and construction of a particle-to-sCO2 heat exchanger test station is described in this paper. The purpose of this test station is to make steady-state measurements of thermal performance with sCO2 as the working fluid. While the test station was initially constructed to test a 20 kW heat exchanger developed under the Gen3 Particle pilot plant (G3P3) project, it also was designed to accommodate testing of other heat exchanger configurations. Improvements for this test station design is based on lessons learned from prior heat exchanger testing. Maximum pressure and temperature ratings are based on the desire to use primarily stainless steel in the construction to reduce cost and lead time of components. Construction of the test station was completed and commissioning and initial testing took place during the October to November 2020 timeframe.

Particle-based heat transfer materials used in concentrating solar power systems benefit from gravity-fed arrangements such as vertically integrated components inside the receiver tower which can eliminate the need for conveyance machinery. However, the amount of particles required for commercial scale systems near 100 MWe can require towers with very thick walls that must be built with high-strength concrete. Cost models for particle-based receiver towers with internal particle storage are being developed in this work and compared to well-established cost models that have been used to estimate tower costs for molten salt systems with external storage tanks. New cost models were developed to accommodate the high-temperature applications required for CSP. Further research is needed to directly compare costs between tower-integrated and external storage. For now, a method is proposed to superimpose increased storage costs with existing molten salt CSP towers. For instances where suitable materials are unavailable or do not meet the structural requirements, ground based storage bins must be used in concert with mechanical conveyance systems. Ground based storage vessels have been shown to be consistent with low thermal energy storage cost and heat loss goals. Ground based storage vessels are well-established in industry.