Publications

Abstract not provided.

The Disposal Research & Development (Disposal R&D) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Spent Fuel & High-Level Waste Disposition is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). A high priority for Disposal R&D is disposal system modeling (Sassani et al. 2023). The Geologic Disposal Safety Assessment (GDSA) work package is charged with developing a disposal system modeling and analysis capability for evaluating generic disposal system performance for nuclear waste in geologic media.

Reducing lifecycle carbon-dioxide (CO₂) and toxic emissions via electrification or switching to carbon-free fuels is not currently feasible for many off-road, rail, and marine applications. This is due to factors including excessive cost, weight, or size of a battery of sufficient capacity to meet the application requirements, lack of infrastructure, insufficient time for recharging, demanding duty cycles, and severe ambient conditions. The guiding vision for the activities described herein is to enable rapid, cost-effective reductions of the environmental impacts of such applications by using improved, high efficiency engine combustion strategies with currently available and emerging low lifecycle-CO₂ fuels (LLCFs). This report summarizes progress toward achieving this vision in two project areas. The first is a Technology Commercialization Fund (TCF) project focused on facilitating the commercialization of ducted fuel injection (DFI) with LLCFs. The second is a more fundamental, Advanced Combustion Engines (ACE) research project focused on elucidating a new strategy called lean mixing-controlled combustion (LMCC) for use with emerging LLCFs.

Sandia National Laboratories has tested and evaluated three Colt broadband seismometers designed and manufactured by Reftek. The purpose of this seismometer evaluation is to measure performance characteristics in areas such as power consumption, sensitivity, frequency response, full scale, self-noise, dynamic range, calibration system response, and passband. The Colt model of sensors are being evaluated to explore the potential for a future seismometer Type Approval process in the International Monitoring System (IMS) of the Comprehensive Nuclear-Test-Ban Treaty (CTBT).

The Burrowing Owl Survey Report for Sandia National Laboratories: 2024 provides data and analysis on occupancy surveys that were conducted during the 2024 survey season. The surveys were conducted in partnership with Kirtland Air Force Base’s Natural Resource Program to maintain and expand on the long-term dataset on burrowing owls on base.

The Z line VISAR system (ZLV) is a spatially-resolved velocimeter that measures surface velocities in high-energy density experiments on the Z Pulsed Power Facility to facilitate the investigation of fusion, power flow, and dynamic material physics. The data measured in these experiments are analyzed with the LineVISAR SMASH class, a MATLAB software suite that provides tools for data importation, streak image correction, spatiotemporal registration, wrapped phase computation, phase unwrapping, shock handling, and velocity calculation. This report overviews the LineVISAR class, discusses its use, provides example implementations, and supplies analysis specifics not typically recorded in journal publications.

The Strategic Petroleum Reserve (SPR) is the world’s largest supply of emergency crude oil. The reserve consists of four sites in Louisiana and Texas. Each site stores crude in deep, underground salt caverns. It is the mission of the SPR’s Enhanced Monitoring Program to examine available sensing data to inform our understanding of each site. This report discusses the monitoring data, processes, and results for each of the four sites for fiscal year 2024.

Herein, we report on the ultrafast photodissociation of nickel tetracarbonyl─a prototypical metal-ligand model system─at 197 nm. Using mid-infrared transient absorption spectroscopy to probe the bound C≡O stretching modes, we find evidence for the picosecond time scale production of highly vibronically excited nickel dicarbonyl and nickel monocarbonyl, in marked contrast with a prior investigation at 193 nm. Further spectral evolution with a 50 ps time constant suggests an additional dissociation step; the absence of any corresponding growth in signal strongly indicates the production of bare Ni, a heretofore unreported product from single-photon excitation of nickel tetracarbonyl. Thus, by probing the deep UV-induced photodynamics of a prototypical metal carbonyl, this Letter adds time-resolved spectroscopic signatures of these dynamics to the sparse literature at high excitation energies.

High-entropy ceramics have garnered interest due to their remarkable hardness, compressive strength, thermal stability, and fracture toughness; yet the discovery of new high-entropy ceramics (out of a tremendous number of possible elemental permutations) still largely requires costly, inefficient, trial-and-error experimental and computational approaches. The entropy forming ability (EFA) factor was recently proposed as a computational descriptor that positively correlates with the likelihood that a 5-metal high-entropy carbide (HECs) will form the desired single phase, homogeneous solid solution; however, discovery of new compositions is computationally expensive. If you consider 8 candidate metals, the HEC EFA approach uses 49 optimizations for each of the 56 unique 5-metal carbides, requiring a total of 2744 costly density functional theory calculations. Here, we describe an orders-of-magnitude more efficient active learning (AL) approach for identifying novel HECs. To begin, we compared numerous methods for generating composition-based feature vectors (e.g., magpie and mat2vec), deployed an ensemble of machine learning (ML) models to generate an average and distribution of predictions, and then utilized the distribution as an uncertainty. We then deployed an AL approach to extract new training data points where the ensemble of ML models predicted a high EFA value or was uncertain of the prediction. Our approach has the combined benefit of decreasing the amount of training data required to reach acceptable prediction qualities and biases the predictions toward identifying HECs with the desired high EFA values, which are tentatively correlated with the formation of single phase HECs. Using this approach, we increased the number of 5-metal carbides screened from 56 to 15,504, revealing 4 compositions with record-high EFA values that were previously unreported in the literature. Our AL framework is also generalizable and could be modified to rationally predict optimized candidate materials/combinations with a wide range of desired properties (e.g., mechanical stability, thermal conductivity).

Photoinitiated polymerization enables spatiotemporal control of reaction conditions and can thereby generate materials with high complexity while consuming minimal energy. Where ring opening metathesis polymerization (ROMP) is concerned, photo-activated processes are typically enabled by chemical inhibition of ruthenium carbenes via the careful design of complexed ligands such that photoactivation can proceed through an isomerization or ligand dissociation event. In this contribution, we have explored a new approach to photoinitiation of ROMP based on physical inhibition through microencapsulation and controlled release of metathesis catalysts. Micron-sized particles of poly(phthalaldehyde) (PPA), catalyst, and photoacid generator were fabricated by spray drying. The particles were dispersed in dicyclopentadiene monomer, after which polymerization was initiated through temperature or UV exposure, both inducing depolymerization of the PPA particles and in situ catalyst release. The monomer/particle dispersions were found to be stable and reproducibly polymerizable with 3 weeks of storage at room temperature. Furthermore, the dispersions can be used for both photo- and thermal-initiated frontal ROMP, yielding a polymerized thermoset of equivalent properties to conventional bulk- and frontally-polymerized analogues. In conclusion, this work will ultimately enable new manufacturing techniques for ROMP-based materials, due to the modular, easily tunable nature of the underlying initiating system and its unparalleled stability.

Recently, the spectral localizer framework has emerged as an efficient approach for classifying topology in photonic systems featuring local nonlinearities and radiative environments. In nonlinear systems, this framework provides rigorous definitions for concepts such as topological solitons and topological dynamics, where a system’s occupation induces a local change in its topology due to nonlinearity. For systems embedded in radiative environments that do not possess a shared bulk spectral gap, this framework enables the identification of local topology and shows that local topological protection is preserved despite the lack of a common gap. However, as the spectral localizer framework is rooted in the mathematics of C*-algebras, and not vector bundles, understanding and using this framework requires developing intuition for a somewhat different set of underlying concepts than those that appear in traditional approaches for classifying material topology. In this tutorial, we introduce the spectral localizer framework from a ground-up perspective and provide physically motivated arguments for understanding its local topological markers and associated local measure of topological protection. In doing so, we provide numerous examples of the framework’s application to a variety of topological classes, including crystalline and higher-order topology. We then show how Maxwell’s equations can be reformulated to be compatible with the spectral localizer framework, including the possibility of radiative boundary conditions. To aid in this introduction, we also provide a physics-oriented introduction to multi-operator pseudospectral methods and numerical K-theory, two mathematical concepts that form the foundation for the spectral localizer framework. Finally, we provide some mathematically oriented comments on the C*-algebraic origins of this framework, including a discussion of real C*-algebras and graded C*-algebras that are necessary for incorporating physical symmetries. Looking forward, we hope that this tutorial will serve as an approachable starting point for learning the foundations of the spectral localizer framework.

Carbon dots have attracted widespread interest for sensing applications based on their low cost, ease of synthesis, and robust optical properties. We investigate structure-function evolution on multiemitter fluorescence patterns for model carbon-nitride dots (CNDs) and their implications on trace-level sensing. Hydrothermally synthesized CNDs with different reaction times were used to determine how specific functionalities and their corresponding fluorescence signatures respond upon the addition of trace-level analytes. Archetype explosives molecules were chosen as a testbed due to similarities in substituent groups or inductive properties (i.e., electron withdrawing), and solution-based assays were performed using ratiometric fluorescence excitation-emission mapping (EEM). Analyte-specific quenching and enhancement responses were observed in EEM landscapes that varied with the CND reaction time. We then used self-organizing map models to examine EEM feature clustering with specific analytes. The results reveal that interactions between carbon-nitride frameworks and molecular-like species dictate response characteristics that may be harnessed to tailor sensor development for specific applications.

There is growing interest in material candidates with properties that can be engineered beyond traditional design limits. Compositionally complex oxides (CCO), often called high entropy oxides, are excellent candidates, wherein a lattice site shares more than four cations, forming single-phase solid solutions with unique properties. However, the nature of compositional complexity in dictating properties remains unclear, with characteristics that are difficult to calculate from first principles. Here, compositional complexity is demonstrated as a tunable parameter in a spin-transition oxide semiconductor La1− x(Nd, Sm, Gd, Y)x/4CoO3, by varying the population x of rare earth cations over 0.00≤ x≤ 0.80. Across the series, increasing complexity is revealed to systematically improve crystallinity, increase the amount of electron versus hole carriers, and tune the spin transition temperature and on-off ratio. At high a population (x = 0.8), Seebeck measurements indicate a crossover from hole-majority to electron-majority conduction without the introduction of conventional electron donors, and tunable complexity is proposed as new method to dope semiconductors. First principles calculations combined with angle resolved photoemission reveal an unconventional doping mechanism of lattice distortions leading to asymmetric hole localization over electrons. Thus, tunable complexity is demonstrated as a facile knob to improve crystallinity, tune electronic transitions, and to dope semiconductors beyond traditional means.

The tunability of metal-organic frameworks (MOFs) makes them exceptional materials for the development of highly selective, low-power sensors for toxic gas detection. Herein, we demonstrate enhanced detection of NO2 gas by a MOF-based electrical impedance sensor made using a unique mixed metal MOF-on-MOF synthesis. A combined experimental and computational study was performed using the exemplar NixMg1-x-MOF-74 to understand the fundamental structure-property relationships behind metal mixing and MOF film synthesis methods on sensor performance. Density functional theory results indicated that the presence of Ni in Mg-MOF-74 increased framework stability and increased the electron density of states at lower energies near the HOMO, as well as enhanced the NO2-Mg adsorption interaction. Impedance data of the NixMg1-x-MOF-74 films with larger Ni contents showed greater impedance change after exposure to 1 ppm of NO2 gas. Furthermore, when synthesized through either a drop-cast or direct solvothermal film growth approach, the monometallic Ni-based sensors had the best performance. However, the mixed metal NixMg1-x-MOF-74 sensors synthesized through a MOF-on-MOF approach resulted in the highest impedance change, outperforming all monometallic Ni-based sensors. In particular, the mixed metal Ni-on-Mg-MOF-74 film was the best-performing sensor with an impedance change of 309 upon trace NO2 exposure. Change in impedance response after NO2 exposure was improved by 52% compared to the best monometallic Ni-on-Ni-MOF-74 sensor. Structural analysis of the Ni-on-Mg film showed that the first Mg-MOF-74 layer acts as a structural template controlling the structural features of the final film after metal exchange with Ni. This led to improved film quality, evidenced by the greater crystallinity and larger MOF grain sizes, and resulted in enhanced sensor performance which was not achievable through other metal mixing methods. Altogether, this study identifies structure-property relationships and synthetic templating methods that inform MOF-based sensor design, allowing for improved detection of toxic compounds.

Harmonic and subharmonic RF injection locking is demonstrated in a terahertz (THz) quantum-cascade vertical-external-cavity surface-emitting laser (QC-VECSEL). By tuning the RF injection frequency around integer multiples and submultiples of the cavity round-trip frequency, different harmonic and subharmonic orders can be excited in the same device. Modulation-dependent behavior of the device has been studied with recorded lasing spectral broadening and locking bandwidths in each case. In particular, harmonic injection locking results in the observation of harmonic spectra with bandwidths over 200 GHz. A semiclassical Maxwell-density matrix formalism has been applied to interpret QC-VECSEL dynamics, which aligns well with experimental observations.

Here, this study investigates neutron-induced displacement damage in Bipolar Junction Transistors (BJTs) using TCAD models informed by Deep-Level-Transient-Spectroscopy (DLTS) data. These models are calibrated and validated against experimental measurements performed at various neutron fluences. Both npn and pnp transistor configurations are studied to analyze the effects of individual traps on carrier recombination and base leakage currents. In npn transistors, deep traps (0.42 eV from the conduction band) dominate at low voltages, while shallow traps (0.17 eV from the conduction band) become prominent at higher voltages. Conversely, pnp transistors have base leakage current predominantly due to deep-level traps. The study observes a notable trend in trap density versus fluence, characterized by a linear relationship on a log-log scale. These insights into defect evolution under radiation conditions are crucial for optimizing semiconductor device reliability and performance in radiation-prone environments.

Conceptual models of smectite hydration include planar (flat) clay layers that undergo stepwise expansion as successive monolayers of water molecules fill the interlayer regions. However, X-ray diffraction (XRD) studies indicate the presence of interstratified hydration states, suggesting non-uniform interlayer hydration in smectites. Additionally, recent theoretical studies have shown that clay layers can adopt bent configurations over nanometer-scale lateral dimensions with minimal effect on mechanical properties. Therefore, in this study we used molecular simulations to evaluate structural properties and water adsorption isotherms for montmorillonite models composed of bent clay layers in mixed hydration states. Results are compared with models consisting of planar clay layers with interstratified hydration states (e.g. 1W–2W). The small degree of bending in these models (up to 1.5 Å of vertical displacement over a 1.3 nm lateral dimension) had little or no effect on bond lengths and angle distributions within the clay layers. Except for models that included dry states, porosities and simulated water adsorption isotherms were nearly identical for bent or flat clay layers with the same averaged layer spacing. Similar agreement was seen with Na- and Ca-exchanged clays. In conclusion, while the small bent models did not retain their configurations during unconstrained molecular dynamics simulation with flexible clay layers, we show that bent structures are stable at much larger length scales by simulating a 41.6×7.1 nm² system that included dehydrated and hydrated regions in the same interlayer.

This report details operations and maintenance (O&M) activities performed across Fiscal Years 2022 through 2024 in support of the continued capabilities of the National Solar Thermal Testing Facility (NSTTF) at Sandia National Laboratories. The NSTTF O&M project is funded by the U.S. Department of Energy Solar Energy Technologies Office (SETO) to support research activities and testing on behalf of external customers at the facility under award number CPS 38491. During the project period, the NSTTF made progress in the areas of site metrics, site maintenance and utilization tracking, and customer engagement. The O&M project also supported special initiatives including procurement of a heat exchanger for particle concentrating solar thermal processes and a scoping and cost study for refurbishment and repair of component in the NSTTF heliostat field.

Theoretically and computationally describing the operation of nanodiodes requires characterizing the transitions between multiple electron emission mechanisms for nanodiodes with complicated geometries. This motivates our development of techniques to determine when simplified theories for individual mechanisms suffice compared to more complete, but more computationally expensive, models. Leveraging recent theories that define a canonical gap distance to translate planar theory to nonplanar diodes, we derive the conditions for the transitions among thermal emission, field emission, and space-charge-limited current density (SCLCD) in vacuum and with collisions for non-Cartesian coordinate systems, including spherical, cylindrical, and prolate spheroidal coordinate systems. Particle-in-cell (PIC) simulations of the current density as a function of applied voltage for a tip-to-plate geometry in vacuum agreed qualitatively with the asymptotes for thermal emission at low voltage and SCLCD at higher voltage using the canonical gap distance. As a result, this demonstrates the utility of this approach for guiding system design and suggests future extensions to save simulation time for more realistic geometries that are more computationally expensive.

The stochastic weighted particle method (SWPM) is a generalization of the Direct Simulation Monte Carlo (DSMC) method where particle weights are variable and dynamic. SWPM is backed by a strong theoretical foundation but has not been critically evaluated for problems of practical interest. A thorough assessment of SWPM for boundary-driven flows reveals significant numerical artifacts near the boundary, notably a diverging heat flux. To correct the boundary heat flux, two modifications to SWPM are proposed: separated grouping and a spatially-dependent weight transfer function. To gauge the relative efficiency of SWPM in comparison to DSMC, a high-Mach-number wheel flow which forms a strong density gradient is also simulated.

Researchers are exploring adding wave energy converters to existing oceanographic buoys to provide a predictable source of renewable power. A ”pitch resonator” power take-off system has been developed that generates power using a geared flywheel system designed to match resonance with the pitching motion of the buoy. However, the novelty of the concept leaves researchers uncertain about various design aspects of the system. This work presents a novel design study of a pitch resonator to inform design decisions for an upcoming deployment of the system. The assessment uses control co-design via WecOptTool to optimize control trajectories for maximal electrical power production while varying five design parameters of the pitch resonator. Given the large search space of the problem, the control trajectories are optimized within a Monte Carlo analysis to identify optimal designs, followed by parameter sweeps around the optimum to identify trends between the design parameters. The gear ratio between the pitch resonator spring and flywheel are found to be the most sensitive design variables to power performance. The assessment also finds similar power generation for various sizes of resonator components, suggesting that correctly designing for optimal control trajectories at resonance is more critical to the design than component sizing.

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α - Te O 2 paratellurite ( P 4 1 2 1 2 ), Te O 2 cotunnite ( P n m a ), and Te O 2 post-cotunnite ( P 2 1 / m ). Based on the discontinuity found in the calculated U s − u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35 km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of Te O 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α - Te O 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

Accurate material characterization and model calibration are pivotal for simulations used for high-consequence engineering decisions. Current characterization and calibration methods (1) use simplified test specimen geometries and global data, (2) cannot guarantee that sufficient characterization data is collected for a specific model of interest, (3) provide only mean parameter values with no uncertainty quantification, and (4) are sequential, inflexible, and time-consuming. This work developed a new paradigm—coined Interlaced Characterization and Calibration (ICC)—which drives forward the state-of-the-art in model calibration by bringing together recent advancements into one improved workflow. The ICC paradigm (1) employs tools to efficiently use full-field data to calibrate high-fidelity material models, (2) aligns the data needed with the data collected by adopting an optimal experimental design protocol, (3) provides uncertainty metrics on the calibrated model parameters, and (4) incorporates these advances into a quasi real-time feedback loop. The ICC framework was validated synthetically with both low-fidelity and high-fidelity simulations paired with several different elastoplastic material models, and was also demonstrated experimentally with an aluminum 6061 cruciform exemplar specimen. Results showed that the ICC framework—in which Bayesian optimal experimental design actively guided the experiment— resulted in calibrations with similar or better accuracy than predetermined experiments based on subject matter expertise. Moreover, the ICC framework produced a complete model calibration— with quantified uncertainties on model parameters—in 1 week, a 5 - 10× increase in efficiency over traditional approaches. Thus, the ICC paradigm improves both the calibration process and quality, by (1) improving efficiency, which increases agility of solid mechanics modeling and enables utilization of computational simulation (CompSim) at earlier stages of the design cycle and (2) providing quantified, and in some cases reduced, parameter uncertainties, which increases confidence in model predictions and supports credible decision making.

The characterization of the neutron, prompt gamma-ray, and delayed gamma-ray radiation fields for the White Sands Missile Range (WSMR) Fast Burst Reactor, also known as molybdenum-alloy Godiva (Molly-G) has been assessed at the 6-inch irradiation location. The neutron energy spectra, uncertainties, and common radiation metrics are presented. Code-dependent recommended constants are given to facilitate the conversion of various dosimetry readings into radiation metrics desired by experimenters. The Molly-G core was designed and configured similarly to Godiva II, as an unreflected, unmoderated, cylindrical annulus of uranium-molybdenum-alloy fuel with molybdenum loading of 10%. At the 6-inch position, the axial fluence maximum is about 2.4×1013 n/cm2 per MJ of reactor energy; about 0.1% of the neutron fluence is below 1 keV and 96% is above 100 keV. The 1-MeV Damage-Equivalent Silicon (DES) fluence is estimated at 2.2×1013 n/cm2 per MJ of reactor energy. The prompt gamma-ray dose is roughly 2.5E+03 rad(Si) per MJ and the delayed gamma-ray dose is about 1.3E+03 rad(Si) per MJ.

Fully-polarimetric synthetic aperture radar (PolSAR) data contain a rich body of elementary scattering physics information that is critically valuable for a broad range of applications and scientific purposes. However, there is a lack of available high-resolution (< 0.3048-m) data available for PolSAR phenomenology research. This article introduces a high-resolution PolSAR data set collected and provided by Sandia National Laboratories (SNL). The data sets were collected to support studying high-resolution scattering physics from different types of clutter and applications such as polarimetric-based terrain classification.