Publications

Abstract not provided.

Sandia National Laboratories (SNL) is advancing technical capabilities used in passive loop seals. The “Puck” seal used a set of International Atomic Energy Agency (IAEA) requirements for new passive loop seals published in 2020 as a design guide. The seal is based on an oxygen-sensitive inner mixture encased in an oxygen-impermeable shell, is monolithic rather than two-part, incorporates self-capturing wire features, contains colored water beads and bubbles formed during processing as unique identifiers (UIDs), and visually indicates tamper (whether from seal body penetration or from seal wire removal) by irreversibly changing the seal body from multi-colored to black. This paper will provide details on the design, development, and testing of Puck seal prototypes.

Downhole logging tools are commonly used to characterize multi-thousand-foot geothermal wells. The elevated temperatures, pressures, and harsh chemical environments present significant challenges for the long-term operation of these tools, especially when real-time data transmission to the surface is required via data cable lines. Teflon-based single or multi-conductor cables with grease-filled cable heads are typically used for downhole tools. However, over extended periods of operation, the grease used to seal the conductors can slowly dissolve into the well fluid, creating electrical shorts and disabling data transmission. Additionally, when temperatures exceed 260 °C, Teflon can soften, potentially allowing parallel conductors to make contact and cause shorts. Between 2009 and 2015, Draka Cableteq USA, now part of the Prysmian Group, developed a multi-conductor/fiber cable and a four-conductor cable capable of operating above 300 °C. While a full study was conducted on the conductor/fiber cable, the evaluation of the four-conductor cable remained incomplete. With the increasing need for long-term high-temperature (HT) operation of logging tools, Sandia National Laboratories is now completing the evaluation of the four-conductor cable. The four-conductor cable has two major novel aspects. Firstly, its glass braid insulation can operate above 300 °C, eliminating the potential for shorts. Secondly, the insulated conductors are encased in metal tubing along the full length of the cable, creating a high-pressure seal between the cable and the tool. This metal tubing eliminates the need for a grease seal, a major limiting factor in the operation time of common cable lines. Sandia National Laboratories will conduct multiple tests to characterize the cable at temperatures above 300 °C and pressures up to 5,000 psi. This cable would enable tools to operate continuously at elevated temperatures, pressures, and in harsh fluids for extended periods, potentially lasting months.

Abstract not provided.

Abstract not provided.

This report presents development of a new meso-scale geomechanics capability using an 8.84-m radius centrifuge for physically simulating underground chemical explosions and stress cage processes in scaled models.

Abstract not provided.

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.

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.

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.

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.

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.

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.

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 in the University of Texas at Austin Nuclear Engineering Teaching Laboratory (NETL) TRIGA reactor for the beam port (BP) 1/5 free-field environment at the 128-inch location adjacent to the core centerline has been accomplished. NETL is being explored as an auxiliary neutron test facility for the Sandia National Laboratories radiation effects sciences research and development campaigns. The NETL reactor is a TRIGA Mark-II pulse and steady-state, above-ground pool-type reactor. NETL is intended as a university research reactor typically used to perform irradiation experiments for students and customers, radioisotope production, as well as a training reactor. Initial criticality of the NETL TRIGA reactor was achieved on March 12, 1992, making it one of the newest test reactor facilities in the US. The neutron energy spectra, uncertainties, and covariance matrices are presented as well as a neutron fluence map of the experiment area of the cavity. For an unmoderated condition, the neutron fluence at the center of BP 1/5, at the adjacent core axial centerline, is about 8.2×1012 n/cm2 per MJ of reactor energy. About 67% of the neutron fluence is below 1 keV and 22% above 100 keV. The 1-MeV Damage-Equivalent Silicon (DES) fluence is roughly 1.6×1012 n/cm2 per MJ of reactor energy.

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.

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.

A neutron fluence map and a total ionizing dose map of the Los Alamos National Laboratory Godiva IV fast burst critical assembly was generated using passive reactor dosimetry, comprised of sulfur pellets and thermoluminescent dosimeters. Godiva IV is an unmoderated, fast burst, critical assembly constructed of approximately 65 kg of highly enriched uranium fuel alloyed with 1.5 % molybdenum for strength. [1] The mapping was performed during a single 75.6 ºC temperature rise burst operation, with the top and sides of the cylindrical Godiva-IV Top Hat covered in passive dosimetry. Dosimetry was placed in a symmetric pattern around the Top Hat, with higher concentrations near the control rods and burst rod. A specific portion of the lower quadrant of the burst rod was mapped to confirm a testing region where the neutron fluence varied by no more than ± 5%. The results will be used to assess the neutron, gamma, and total ionizing dose environment in three-dimensional space around the assembly for higher fidelity experiment placement, active dosimetry positioning, and radiation field characterization.

The finite element method can be used to compute accurate stress intensity factors (SIFs) for cracks with complex geometries and boundary conditions. In contrast, handbook solutions act as surrogate SIF models that provide significantly faster evaluation times. However, the development of conventional surrogate SIF models relies on manual development based on low-order parameterizations. This limits surrogate model accuracy and generalizability. In this paper, we develop a framework for the automated development of mechanics-guided handbook SIF solutions by using interpretable machine learning via genetic programming for symbolic regression (GPSR). Formalizing the mechanics-based approach of Raju and Newman, SIF training data is decomposed into multiple subsets. This decomposition enables parallel GPSR model development of subfunctions, each of which accounts for specific geometrical corrections with respect to a known analytical model. Using this mechanics-based approach with GPSR allows for equations to be learned with improved accuracy and reduced complexity relative to the Raju Newman equations while maintaining the inherent interpretability of mathematical expressions. In this paper, we present equations that match the complexity of the Raju Newman equations while having reduced error, as well as equations with similar errors and reduced complexity.

High-entropy materials (HEMs) emerged as promising candidates for a diverse array of chemical transformations, including CO2 utilization. However, traditional HEMs catalysts are nonporous, limiting their activity to surface sites. Designing HEMs with intrinsic porosity can open the door toward enhanced reactivity while maintaining the many benefits of high configurational entropy. Here, a synergistic experimental, analytical, and theoretical approach to design the first high-entropy metal-organic frameworks (HEMOFs) derived from polynuclear metal clusters is implemented, a novel class of porous HEMs that is highly active for CO2 fixation under mild conditions and short reaction times, outperforming existing heterogeneous catalysts. HEMOFs with up to 15 distinct metals are synthesized (the highest number of metals ever incorporated into a single MOF) and, for the first time, homogenous metal mixing within individual clusters is directly observed via high-resolution scanning transmission electron microscopy. Importantly, density functional theory studies provide unprecedented insight into the electronic structures of HEMOFs, demonstrating that the density of states in heterometallic clusters is highly sensitive to metal composition. This work dramatically advances HEMOF materials design, paving the way for further exploration of HEMs and offers new avenues for the development of multifunctional materials with tailored properties for a wide range of applications.

The cells in battery energy storage systems are monitored, protected, and controlled by battery management systems whose sensors are susceptible to cyberattacks. False data injection attacks (FDIAs) targeting batteries’ voltage sensors affect cell protection functions and the estimation of critical battery states like the state of charge (SoC). Inaccurate SoC estimation could result in battery overcharging and over discharging, which can have disastrous consequences on grid operations. This paper proposes a three-pronged online and offline method to detect, identify, and classify FDIAs corrupting the voltage sensors of a battery stack. To accurately model the dynamics of the series-connected cells a single particle model is used and to estimate the SoC, the unscented Kalman filter is employed. FDIA detection, identification, and classification was accomplished using a tuned cumulative sum (CUSUM) algorithm, which was compared with a baseline method, the chi-squared error detector. Online simulations and offline batch simulations were performed to determine the effectiveness of the proposed approach. Throughout the batch simulations, the CUSUM algorithm detected attacks, with no false positives, in 99.83% of cases, identified the corrupted sensor in 97% of cases, and determined if the attack was positively or negatively biased in 97% of cases.

Vertically aligned nanocomposites (VANs) are unique thin films with vertical nanostructures embedded in a matrix material, allowing for the integration of two distinct materials. These nanocomposites offer novel combined physical properties, such as nanocomposite-based multiferroics and strongly coupled physical properties, such as magneto-optic coupling. Much work has been conducted in exploring different two-phase combinations and various processing conditions to achieve novel tunable properties that cannot be obtained by any single-phase material alone. In this work, the target configuration effects are explored for the growth of LaFeO3-CoNi2O4 VANs. Both mixed and pie-shaped targets are utilized to compare the target configuration effects on the phase separation, morphology tuning, and their resulting physical properties, including optical and magnetic properties. This work suggests that the target configuration is another important parameter for achieving the desired VAN morphology and can be used to design different two-phase VANs with tailorable properties.

This is the manual for the TeMPI Shim library, whose goal is to facilitate scientific apps being loosely coupled through MPI. TeMPI Shim can be used to create and manage Message Passing Interface (MPI, see [MPI]) Communicators for Single Program, Multiple Data (SPMD, see [SPMD]) and Multiple Programs, Multiple Data (MPMD, see [MPMD]) workflows. MPI provides the MPI_APPNUM command number to each of the processes within its context. This number, starting at 0, is the application number being executed. In a case where 3 applications are being run in MPMD mode, there would be MPI_APPNUM values of 0, 1, and 2. TeMPI Shim creates intra- and inter-communicators between each pair of MPI_APPNUM values. In the aforementioned case, application 0 would have intra- and inter-communicators to speak to itself, application 1, and application 2. This is replicated for each of the applications. Additionally, TeMPI Shim creates intra- and inter-communicators for the first MPI rank of each application to directly communicate only with each other. Finally, TeMPI Shim creates its own copy of the default world communicator, i.e., MPI_COMM_WORLD. In the case where there is only a single application, it will have the communicators to only speak with itself. TeMPI Shim is useful in this case since it is considered good practice for MPI applications to copy the default world communicator and reference this copy (see [Duplicate World]_) anyways. Ultimately, it provides value independent of the number of applications present.

Thermal spray processes can benefit from cooling to maintain substrate temper, reduce processing times, and manage thermally induced residual stresses. “Plume quenching” is a plume-targeted cooling technique which has been shown to reduce substrate temperatures by redirection of hot plume gases using a lateral argon curtain injected into the plume, while limiting interaction with the substrate or affecting coating properties. Here, this study explores the use of this technique for residual stress management by reducing the thermally driven component in nickel and tantalum coatings on titanium and aluminum substrates. The in-situ residual stress profiles were measured for all substrate and coating pairings during spraying and cooling, and the deposition and thermal stresses recorded. For substrate and coating pairings where the predominant component of residual stress was thermal (driven by a large difference in coefficient of thermal expansion, Δα, between coating and substrate), plume quenching reduced both the thermal stress and the final stress state of the coating. This was seen primarily in tantalum on aluminum coatings where the Δα was -17 × 10^-6 /°C, and thermal stress was reduced by 7.5% and 22.4% for the plume quenching rates of 50 and 100 slpm, respectively.

Thermal spray processes can benefit from cooling to maintain substrate temper, reduce processing times, and manage thermally induced residual stresses. “Plume quenching” is a plume-targeted cooling technique which has been shown to reduce substrate temperatures by redirection of hot plume gases using a lateral argon curtain injected into the plume, while limiting interaction with the substrate or affecting coating properties. Here, this study explores the use of this technique for residual stress management by reducing the thermally driven component in nickel and tantalum coatings on titanium and aluminum substrates. The in-situ residual stress profiles were measured for all substrate and coating pairings during spraying and cooling, and the deposition and thermal stresses recorded. For substrate and coating pairings where the predominant component of residual stress was thermal (driven by a large difference in coefficient of thermal expansion, Δα, between coating and substrate), plume quenching reduced both the thermal stress and the final stress state of the coating. This was seen primarily in tantalum on aluminum coatings where the Δα was -17 × 10^-6 /°C, and thermal stress was reduced by 7.5% and 22.4% for the plume quenching rates of 50 and 100 slpm, respectively.

Abstract not provided.

This report describes computational modeling in the HyRAM+ software for study of hydrogen behavior in two common scenarios. First, models of unignited plumes exiting a vent stack were considered. It was shown that entrainment and vent backpressure were major factors in plume physics. Second, HyRAM+ was extended to include wind effects on plumes by modifying the plume momentum balance and entrainment modeling. Use of the model showed that plume shape and length changed with wind speed and direction; in all cases, wind causes a shortening of the plume along the streamline. While the no-wind case in HyRAM+ has been validated and the newly developed wind model was fitted to very limited experimental data, more controlled experimental configurations would help validate the models and ensure accurate simulation of hydrogen plume behavior for vent stack releases or in wind.

In this paper, we present a Riemannian geometric derivation of the governing equations of motion of nonholonomic dynamic systems. A geometric form of the work-energy principle is first derived. The geometric form can be realized in appropriate generalized quantities, and the independent equations of motion can be obtained if the subspace of generalized speeds allowable by nonholonomic constraints can be determined. We provide a geometric perspective of the governing equations of motion and demonstrate its effectiveness in studying dynamic systems subjected to nonholonomic constraints.

This report investigates the integration of infrasound and seismic data to improve event localization accuracy, specifically focusing on a surface explosion at the Utah Training and Testing Range (UTTR). Utilizing the Seismoacoustic Bayesian Event Locator (SABEL) framework, we incorporated atmospheric specifications derived from Ground to Space (G2S) profiles to enhance celerity-range priors. Our analysis revealed that while the combination of infrasound and seismic observations significantly reduced localization uncertainty, challenges remained, particularly with returns at distances less than 200 km from the source and the influence of specific observations on location estimates. The results indicate that broader celerity distributions, such as those from Blom et al. (2020), facilitate better alignment with ground truth locations compared to narrower models. Overall, this work demonstrates the promise of seismoacoustic approaches in refining event localization and highlights the need for further exploration of celerity-range models to ensure reliable outcomes.

The valorization and dissolution of lignin using ionic liquids (ILs) is critical for developing sustainable biorefineries and a circular bioeconomy. This review aims to critically assess the current state of computational and machine learning methods for understanding and optimizing IL-based lignin dissolution and valorization processes reported since 2022. The paper examines various computational approaches, from quantum chemistry to machine learning, highlighting their strengths, limitations, and recent advances in predicting and optimizing lignin-IL interactions. Key themes include the challenges in accurately modeling lignin’s complex structure, the development of efficient screening methodologies for ionic liquids to enhance lignin dissolution and valorization processes, and the integration of machine learning with quantum calculations. These computational advances will drive progress in IL-based lignin valorization by providing deeper molecular-level insights and facilitating the rapid screening of novel IL-lignin systems.

In recent years, Native tribes across the country have embraced solar energy to power their communities. However, in approximately 25-30 years after installation, the solar panels will enter end-of-life (EOL) status, meaning that they will begin to operate less efficiently than upon initial implementation. At this point, tribes can determine the outcome: they can continue to operate the array, or they can dispose of the panels. Recycling is an option that promotes the use of high-value recyclable materials in solar panels in the circular economy. There are several companies across the US that offer recycling of solar panels, but at various costs. It is important that a hazardous waste expert makes an initial assessment of the solar array to identify any solar panels that cannot be recycled. Currently, recycling solar panels is not cost-effective or environmentally sustainable, but significant research into improving the cost and carbon footprint is ongoing.

Blending hydrogen into existing natural gas pipelines is a current initiative for large-scale hydrogen deployment in the United States. Understanding changes in hazards and risk for fuel blends can help inform regulations. A differential risk assessment was conducted to understand how the safety profile of fires and unignited flammable plumes of leaking fuel would change based on blend composition. The area affected by harmful heat fluxes from jet fires was found to be larger for blends with more natural gas in them, while the distance affected by unignited flammable fuel concentrations was found to be larger for blends with more hydrogen in them. Component leak frequencies and ignition probabilities greatly impact both heat flux risk and the flammable plume risk. The collection of additional blend-specific industry data could help understand which assumptions are most valid to make when calculating risk.

This document serves as a comprehensive final report for the Single-Volume Scatter Camera (SVSC) project. Most of the work over the course of this project has been documented in journal articles, conference papers, and other reports. We therefore reference those materials for detailed presentation of our technical results. The most recent efforts on the project have not been published and are presented in detail here. We present characterizations of two neutron scatter camera prototypes; one using a monolithic geometry, and one using an optically segmented geometry. Both detectors employ plastic scintillator with SiPM-based readout. For the monolithic prototype, we present calibrations of several detector parameters made using a set of dark counts. In particular, we employ a coincidence analysis to characterize the distribution external optical crosstalk, along with the relative timing of the different SiPMs. For the optically segmented prototype we present the results of calibrations made using a tagged ²²Na source, as well as the results of an imaging measurement made using a tagged AmBe source.

Sandia National Labs will provide Technical Assistance and consulting to Requester on a series of Requester’s test structures and devices to characterize and aid in the determination of optimal fabrication parameters via fabrication of functional test structures and cells. Working with previously established processes and test die modules, Sandia will construct a series of test chips with varying materials properties to allow the electrical characterization and component testing of Requester’s test structures and devices. Material variations will be selected for optimization of yield and desired target operation parameters. A deliverable in the form of a report will be written. Sandia intends to deliver testable die to the Requester.

The H3D was characterized using GADRAS 19.3.5 using data sent by the IAEA. Everything related to the measurement setup was provided by the IAEA. The sources measured were detailed and the distance from the detector was also provided. The data sheets of the sources were used to determine the source activity at the time of calibration, as needed for characterization.

Abstract not provided.

Pozzolans rich in silica and alumina react with lime to form cementing compounds and are incorporated into portland cement as supplementary cementitious materials (SCMs). However, pozzolanic reactions progress slower than portland cement hydration, limiting their use in modern construction due to insufficient early-age strength. Hence, alternative SCMs that enable faster pozzolanic reactions are necessary including synthetic zeolites, which have high surface areas and compositional purity that indicate the possibility of rapid pozzolanic reactivity. Synthetic zeolites with varying cation composition (Na-zeolite, H-zeolite), SiO2/Al2O3 ratio, and framework type were evaluated for pozzolanic reactivity via Ca(OH)2 consumption using ion exchange and in-situ X-ray diffraction experiments. Na-zeolites exhibited limited exchange reactions with KOH and Ca(OH)2 due to the occupancy of acid sites by Na+ and hydroxyl groups. Meanwhile, H-zeolites readily adsorbed K+ and Ca2+ from a hydroxide solution by exchanging cations with H+ at Brønsted acid sites or cation adsorption at vacant acid sites. By adsorbing cations, the H-zeolite reduced the pH and increased Ca2+ solubility to promote pozzolanic reactions in a system where Ca(OH)2 dissolution/diffusion was a rate limiting factor. High H-zeolite reactivity resulted in 0.8 g of Ca(OH)2 consumed per 1 g of zeolites after 16 h of reaction versus 0.4 g of Ca(OH)2 consumed per 1 g of Na-zeolite. The H-zeolite modulated the pore fluid alkalinity and created a low-density amorphous silicate phase via mechanisms analogous to two-step C-S-H nucleation experiments. Controlling these reaction mechanisms is key to developing next generation pozzolanic cementitious systems with comparable hydration rates to portland cement.

The deliberate introduction of stacking faults has been found to induce properties markedly distinct from those of perfectly stacked particles or sheets, with even minute deviations leading to significant alterations in material characteristics. In this review article, we discuss the effects of stacking faults – both linear and rotational – on surface structures influencing ion adsorption, particle–particle affinity governing crystal growth, friction, and electronic properties. Our comprehensive overview sheds light on the multifaceted impact of stacking faults on these phenomena, offering insights that bridge multiple disciplines, and provides numerous interdisciplinary research questions, paving the way for future innovations.

The rapid expansion of photovoltaic (PV) systems, particularly inverters, has introduced new cybersecurity challenges that threaten both local operations as well as the broader electrical grid’s stability. PV inverters, integrated into critical energy infrastructure are potential targets for cyber attacks due to vulnerabilities in firmware, remote access systems, and communication protocols. The Coordinated Vulnerability Disclosure (CVD) process, as defined by the Cybersecurity and Infrastructure Security Agency (CISA), provides a framework for identifying, reporting, and addressing these vulnerabilities in a transparent and collaborative manner. This report outlines the CVD process as it applies to PV systems, detailing the roles of key stakeholders, such as manufacturers, grid operators, and security researchers. The report also highlights specific challenges in managing vulnerabilities for new and legacy PV systems, which includes those introduced by insecure communications and third-party supply chain components. By adhering to the CVD process, the PV industry can mitigate cybersecurity risks, ensure regulatory compliance, and maintain consumer trust, while safeguarding the operational resilience of the energy grid. Ultimately, the effective coordination of vulnerability management is crucial for securing the future of PV systems within the critical electric grid infrastructure landscape.

Abstract not provided.

Classical short-circuit programs that linearize the power network are no longer applicable for inverter based resources (IBRs), necessitating an iterative approach. Phasor domain programs can model the IBRs using an iterative approach considering nonlinear fault responses. In phasor domain models, the IBR can be modeled as a voltage controlled current source (VCCS) in tabular form with positive, negative, and zero sequence information for balanced and unbalance short-circuit faults. In the VCCS modeling of the IBR, positive and negative incremental reactive current, also known as the k-factor, plays an important role in short circuit program convergence. In this work, a few approaches: conventional VCCS modeling with a k-factor of 2, conventional VCCS modeling with a k-factor of 2 with modified pre-fault voltages, VCCS characteristics based on the power flow solution with a k-factor of 2, and VCCS characteristics based on the power flow solution with a k-factor of 1 are investigated for short circuit program convergence under higher IBR penetration. The IEEE 39 bus New England Test System is taken as the test system, and simulations are carried out in PSS®CAPE 15.0.26 simulation software. Simulation results demonstrate that IBR penetration is higher for the VCCS model, which corresponds to the power flow solution with k-factor 1, compared to other approaches.