Publications

Abstract not provided.

We consider estimating the discretization error in a nonlinear functional J(u) in the setting of an abstract variational problem: find u∈V such that B(u,φ)=L(φ)∀φ∈V, as approximated by a Galerkin finite element method. Here, V is a Hilbert space, B(⋅,⋅) is a bilinear form, and L(⋅) is a linear functional. We consider well-known error estimates η of the form J(u)−J(uh)≈η=L(z)−B(uh,z), where uh denotes a finite element approximation to u, and z denotes the solution to an auxiliary adjoint variational problem. We show that there exist nonlinear functionals for which error estimates of this form are not reliable, even in the presence of an exact adjoint solution z. An estimate η is said to be reliable if there exists a constant C∈R>0 independent of uh such that |J(u)−J(uh)|≤C|η|. We present several example pairs of bilinear forms and nonlinear functionals where reliability of η is not achieved.

Binary Code Similarity Analysis (BCSA) has a wide spectrum of applications, including plagiarism detection, vulnerability discovery, and malware analysis, thus drawing significant attention from the security community. However, conventional techniques often face challenges in balancing both accuracy and scalability simultaneously. To overcome these existing problems, a surge of deep learning-based work has been recently proposed. Unfortunately, many researchers still find it extremely difficult to conduct relevant studies or extend existing approaches. First, prior work typically relies on proprietary benchmark without making the entire dataset publicly accessible. Consequently, a large-scale, well-labeled dataset for binary code similarity analysis remains precious and scarce. Moreover, previous work has primarily focused on comparing at the function level, rather than exploring other finer granularities. Therefore, we argue that the lack of a fine-grained dataset for BCSA leaves a critical gap in current research. To address these challenges, we construct a benchmark dataset for fine-grained binary code similarity analysis called BinSimDB, which contains equivalent pairs of smaller binary code snippets, such as basic blocks. Specifically, we propose BMerge and BPair algorithms to bridge the discrepancies between two binary code snippets caused by different optimization levels or platforms. Furthermore, we empirically study the properties of our dataset and evaluate its effectiveness for the BCSA research. The experimental results demonstrate that BinSimDB significantly improves the performance of binary code similarity comparison.

Scripta Materialia

Journal of Computational and Applied Mathematics

Nano-hydrides have been predicted to precipitate at the core of edge dislocations in the Ni-H system, a mechanism that may promote hydrogen embrittlement. However, nano-hydride nucleation, growth, and effects on dislocation behavior have seldom been explored. This work combines molecular dynamics grand canonical Monte Carlo (MD-GCMC) simulations and continuum modeling to uncover a wide range of phenomena linked to dislocation nano-hydrides. Simulations reveal that nano-hydrides can be stabilized at dislocation cores with all character angles, including screw segments, due to the hydrostatic stresses around the cores of the Shockley partials. Nano-hydride nucleation takes place in these regions, and growth is dictated by the character angle θ of the perfect dislocation. The equilibrium stacking fault width deq varies dynamically to increase the local hydrostatic stress field and facilitate the formation of the nano-hydride, forming a constriction-like feature and leading to three distinct behaviors: deq decreases for θ>30°, deq remains unchanged for θ=30°, and deq increases for θ<30°. Remote hydrostatic and Escaig stresses are also shown to influence the nucleation stage, implying stress concentrations such as those ahead of crack tips may facilitate nano-hydride precipitation. Moreover, we identify a new hydrogen-induced 60° dislocation reaction that emits a Shockley partial on a conjugate plane, with potential implications for twin nucleation. Testable predictions from this study are then used to reinterpret previous results from the literature. These findings provide a comprehensive framework to assess nano-hydride formation and evolution at dislocations in nickel and other face-centered cubic metals, with important implications to hydrogen embrittlement.

Although synthetic biology can produce valuable chemicals in a renewable manner, its progress is still hindered by a lack of predictive capabilities. Media optimization is a critical, and often overlooked, process which is essential to obtain the titers, rates and yields needed for commercial viability. Here, we present a molecule- and host-agnostic active learning process for media optimization that is enabled by a fast and highly repeatable semi-automated pipeline. Its application yielded 60% and 70% increases in titer, and 350% increase in process yield in three different campaigns for flaviolin production in Pseudomonas putida KT2440. Explainable Artificial Intelligence techniques pinpointed that, surprisingly, common salt (NaCl) is the most important component influencing production. The optimal salt concentration is very high, comparable to seawater and close to the limits that P. putida can tolerate. The availability of fast Design-Build-Test-Learn (DBTL) cycles allowed us to show that performance improvements for active learning are rarely monotonous. This work illustrates how machine learning and automation can change the paradigm of current synthetic biology research to make it more effective and informative, and suggests a cost-effective and underexploited strategy to facilitate the high titers, rates and yields essential for commercial viability.

A simple co-design example in a reduced parameter space is presented for an oscillating flap device. Initially, the WEC geometry and mass properties are considered along with drivetrain gear ratio, inertia, motor constant and stiffness under both PI and optimal control. This parameter space is reduced to those to which performance is most sensitive for a fixed geometry. The gear ratio, drivetrain stiffness, and flap mass are found to be the most impactful design criteria as they can create orders of magnitude variations in power performance. The performance of the optimized system is compared with several sub-optimal variants in terms of electrical and mechanical power capture, transmission coefficients, and transducer power gain. Notably, though substantial power capture improvements are demonstrated when an optimal controller is employed, this power capture remains sensitive to appropriate selections of drivetrain and flap design parameters, implying that control co-design procedures remain necessary for high-performing WECs. A number of practical caveats and extensions to the presented co-design methodology are suggested, including the characterization of system static friction, especially in the presence of high gear ratios.

Magnetic anisotropy is essential for many applications of ferromagnetic/ferrimagnetic materials, including permanent magnets and magnetic recording media. Attempts have been made recently to build up 3-D nanoparticle and quantum dot assemblies, however, it is not understood yet if a nanoparticle assembly can possess high magnetic anisotropy with low anisotropic materials. In this article, we report our discovery of high magnetic anisotropy resulted from Fe3O4 nanoparticle chains. We started with closely-packed nanoparticle assemblies of spherical Fe3O4 nanoparticles that exhibit low magnetocrystalline anisotropy and shape anisotropy, and corresponding negligible coercivity. When the nanoparticle assemblies are compressed under pressure, they form bundles or arrays that consist of Fe3O4 chains with a length scale of several hundred nanometers. Magnetic measurements show that these Fe3O4 chain arrays possess a high uniaxial magnetic anisotropy (Keff ~ 2.9×105J/m³) and significant magnetic coercivity. Our simulations reveal that interparticle magnetic dipolar interactions contribute to this type of superstructure magnetic anisotropy. This study demonstrates the feasibility and approaches to create “patterned” high magnetic anisotropy in nanoparticle superstructures/assemblies.

Viral hemorrhagic Lassa fever (LF), caused by Lassa virus (LASV), is a significant public health concern endemic in West Africa with high morbidity and mortality rates, limited treatment options, and potential for international spread. Despite advances in interrogating its epidemiology and clinical manifestations, the molecular mechanisms driving pathogenesis of LASV and other arenaviruses remain incompletely understood. This review synthesizes current knowledge regarding the role of LASV host-virus interactions in mediating the pathogenesis of LF, with emphasis on interactions between viral and host proteins. Through investigation of these critical protein–protein interactions, we identify potential therapeutic targets and discuss their implications for development of medical countermeasures including antiviral drugs. This review provides an update in recent literature of significant LASV host-virus interactions important in informing the development of targeted therapies and improving clinical outcomes for LF patients. Knowledge gaps are highlighted as opportunities for future research efforts that would advance the field of LASV and arenavirus pathogenesis.

Material design and accessible manufacturing are often at odds with each other, calling for creative solutions to adapt high-performance materials to available processes. This challenge is represented well by in-mold electronics, an innovative approach to the manufacture of 3D circuitry and electronic components that offers game-changing advantages. In-mold electronics relies on vacuum forming processes, which are historically limited to thermoplastics. Extending these methods to include thermosets would enable manufacturing of robust components with desirable properties. Here, we provide a solution to make thermoset materials amenable to vacuum forming. Specifically, an ambient polymerization is used to transition a liquid monomeric solution to an elastomeric gel. These free-standing gels can then be vacuum formed, and the reaction can be completed via frontal polymerization. Thermoset materials produced with this method have properties that provide benefits over traditionally employed thermoplastic substrates and enable 3D device integration into environmentally demanding architectural, automotive, and extraterrestrial structures.

Recent trends within computational and data sciences show an increasing recognition and adoption of computational workflows as tools for productivity and reproducibility that also democratize access to platforms and processing know-how. As digital objects to be shared, discovered, and reused, computational workflows benefit from the FAIR principles, which stand for Findable, Accessible, Interoperable, and Reusable. The Workflows Community Initiative’s FAIR Workflows Working Group (WCI-FW), a global and open community of researchers and developers working with computational workflows across disciplines and domains, has systematically addressed the application of both FAIR data and software principles to computational workflows. We present recommendations with commentary that reflects our discussions and justifies our choices and adaptations. These are offered to workflow users and authors, workflow management system developers, and providers of workflow services as guidelines for adoption and fodder for discussion. The FAIR recommendations for workflows that we propose in this paper will maximize their value as research assets and facilitate their adoption by the wider community.

Ship tracks, long thin artificial cloud features formed from the pollutants in ship exhaust, are satellite-observable examples of aerosol-cloud interactions (ACI) that can lead to increased cloud albedo and thus increased solar reflectivity, phenomena of interest in solar radiation management. In addition to ship tracks being of interest to meteorologists and policy makers, their observed cloud perturbations provide benchmark evidence of ACI that remain poorly captured by climate models. To broadly analyze the effects of ship tracks, high-resolution satellite imagery data highlighting their presence are required. To support this, we provide a hand labelled dataset to serve as a benchmark for a variety of subsequent analyses. Established from a previous dataset that identified ship track presence using NASA’s MODIS Aqua satellite imager, our first-of-its-kind dataset is comprised of image masks: capturing full ship track regions, including their contours, emission points and dispersive patterns. In total, 300 images, or around 2,500 masked ship tracks, observed under varying conditions are provided, and may facilitate training of machine learning algorithms to automate extraction.

International Journal of Engineering Sciences

Abstract not provided.

Nuclear Materials and Energy

Low-velocity impact of hybrid metal-composite structures was investigated experimentally and computationally. Composite laminates consisting of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) were joined with a 6061-T6 aluminum plate using an epoxy adhesive. Two variations of the structure were studied; one consisting of all plies oriented at 0° and one consisting of all plies oriented at 45°. A drop tower was used to impact structures at a range of energies, including energies above and below the threshold at which the aluminum layer was perforated. Numerical simulations were implemented using Sierra/SM, an in-house transient dynamics finite element code developed at Sandia National Laboratories. A Hosford plasticity model was used to describe the response of the aluminum layer. A newly implemented orthotropic continuum damage mechanics (CDM) constitutive model was used to represent the composite laminate. This 3D-CDM model was compared to a cohesive zone model (2D-CDM/CZM) to investigate efficacy of aluminum perforation energy prediction, delamination prediction, and computational cost. Accuracy of each model was evaluated using the experimental results. Each showed good agreement with the tests for both the force and velocity histories, as well as the observed damage mechanisms. The 2D-CDM/CZM model was marginally more accurate in capturing both the composite and aluminum behavior — this model averaged error percentages of −11.2% and 10.8% for residual velocity and peak force, respectively. Meanwhile, the 3D-CDM model predictions yielded average error percentages of −35.5% (velocity) and 22.6% (force). However, the 3D-CDM model generally resulted in a decreased computational cost; the average run time was 14% shorter than the 2D-CDM/CZM model and 3x as many timesteps per hour were computed using the same computational resources. New experimental data on the impact and perforation resistance of metal-composite laminates is presented in addition to numerical predictions of the impact behavior.

Background: Non-invasive simulations of coronary hemodynamics have improved clinical risk stratification and treatment outcomes for coronary artery disease, compared to relying on anatomical imaging alone. However, simulations typically use empirical approaches to distribute total coronary flow amongst the arteries in the coronary tree, which ignores patient variability, the presence of disease, and other clinical factors. Further, uncertainty in the clinical data often remains unaccounted for in the modeling pipeline. Objective: We present an end-to-end uncertainty-aware pipeline to (1) personalize coronary flow simulations by incorporating vessel-specific coronary flows as well as cardiac function; and (2) predict clinical and biomechanical quantities of interest with improved precision, while accounting for uncertainty in the clinical data. Methods: We assimilate patient-specific measurements of myocardial blood flow from clinical CT myocardial perfusion imaging to estimate branch-specific coronary artery flows. Simulated noise in the clinical data is used to estimate the joint posterior distributions of the model parameters using adaptive Markov Chain Monte Carlo sampling. Additionally, the posterior predictive distribution for the relevant quantities of interest is determined using a new approach combining multi-fidelity Monte Carlo estimation with non-linear, data-driven dimensionality reduction. This leads to improved correlations between high- and low-fidelity model outputs. Results: Our framework accurately recapitulates clinically measured cardiac function as well as branch-specific coronary flows under measurement noise uncertainty. We observe substantial reductions in confidence intervals for estimated quantities of interest compared to single-fidelity Monte Carlo estimation and state-of-the-art multi-fidelity Monte Carlo methods. This holds especially true for quantities of interest that showed limited correlation between the low- and high-fidelity model predictions. In addition, the proposed multi-fidelity Monte Carlo estimators are significantly cheaper to compute than traditional estimators, under a specified confidence level or variance. Conclusions: The proposed pipeline for personalized and uncertainty-aware predictions of coronary hemodynamics is based on routine clinical measurements and recently developed techniques for CT myocardial perfusion imaging. The proposed pipeline offers significant improvements in precision and reduction in computational cost.

This paper introduces the topics that need consideration to allow an effective feasibility analysis of electron beams and/or X-rays technologies prior to potential investment and implementation of these technologies for industrial applications. It also highlights considerations for planning, developing and construction of an electron beam or X-ray facility The focus of this paper is on business topics and it provides details of where additional technical information can be found. The paper also describes current applications for electron beam and X-ray technologies to demonstrate their potential.

The demand for low-cost, low-energy, and highly selective gas capture and separations is an ongoing driver of porous material development. Porous liquids have been identified as a promising gas separation material by creating permanent porosity in inorganic solvents through inclusion of nanoporous materials that sterically exclude solvent from their internal porosity. Among the nanoporous materials that can be used to form porous liquids, porous-organic cages (POCs) have been one of the most popular due to the inherent tunability of POCs. “Scrambled” POCs with varying functionalities on the POC vertices have been developed and incorporated into porous liquid compositions, increasing their gas adsorption capacity. An unexplored avenue to tailor the properties of porous liquids is through scrambling the functionality of the core of the POC. Therefore, we have synthesized a new POC, a CC3-OH derivative with scrambled hydroxides on the core and evaluated the impact on the CO2 uptake capacity in silicon oil-based porous liquids. Core scrambling of the POC resulted in a twofold increase CO2 adsorption capacity in the porous liquid, an emergent property that is a dramatic increase beyond a linear combination of the gas adsorption capacity of the neat solvent and the POC. Density functional theory modeling of the CC3 POC and its hydroxide-based derivatives identified that free rotation of the linker hydroxide allowed for forced interaction between the CO2 molecule and the hydroxide in the pore window. Solvation of the POC may release scrambled core hydroxides from intramolecular bonding with a neighboring imine, allowing for increased gas uptake in the porous liquid over the neat POC. These results identify a key structural relationship of POCs that enables emergent properties in porous liquids and can guide future development of liquid phase gas capture and separation materials for environmental and industrial applications.

Distributed acoustic sensing (DAS) systems, which use existing telecommunication fibers, offer high-resolution capabilities ideal for recording anthropogenic sources. However, the complexity of urban environments and the large amount of data recorded by DAS require automated methods to efficiently detect and categorize anthropogenic sources. We evaluate how well three machine learning models (k-nearest neighbor [k-NN], convolutional neural networks, and recurrent-convolutional neural networks) can identify various anthropogenic sources recorded by DAS. Our findings reveal that both k-NN and neural network methods perform well in high signal-to-noise ratio (SNR) settings. However, their accuracy decreases at SNRs < 4. We also use Kalman filtering, a form of Bayesian inference, on backprojected locations of these sources to recover locations that generally fall within standard smartphone Global Positioning System errors. By combining machine learning and Kalman filter results, we calculate a multidimensional model of moving anthropogenic sources. These results demonstrate the potential of DAS data in urban seismology for accurately identifying and locating such sources. Depending on the research objectives, these sources can be further studied or filtered out to improve the quality of seismic data for earthquake studies. Such methods provide a valuable tool for urban seismology and seismic hazard analysis.

A predictive analytical model is presented for the electrical conductivity of multi-principal-element alloys (MPEAs), including those containing aluminum, transition metals, and refractory metals. Given that the lattice parameter of the Wigner-Seitz cell of an MPEA is similarly variable to a bulk metallic glass, it is postulated that electron scattering can be approximated by a series of two-level systems. The resulting reduced-order model enabled an accurate determination of electrical resistivity and electron thermal conductivity based on the scattering of electrons in a two-level system across a Bloch-potential-based virtual crystal approximation. Model results are compared to experimental four-point probe electrical resistivity measurements between 300 K and 700 K for Al0.3CoCrCuFeNi, CoCrFeMnNi, (CoCrFeMnNi)0.98W0.02, (CoCrFeMnNi)0.95W0.05, and Nb4Ta4V3Ti, for model validation.

Remote coastal communities, which could be early adopters of wave energy projects, have concerns over costs, conflicts, and potential risks of development. Designers and developers are challenged to address these community concerns as they continue to develop wave energy technologies. One potential means of reducing costs, conflicts, and risks, especially for demonstration and pilot-scale projects, could be planning a deployment that operates for only a portion of the year—a seasonal deployment. In this paper we examine the impacts of a seasonal deployment in terms of cost, electricity production, operations and maintenance, environmental impacts, and community benefits. We take a holistic, comparative approach to feasibility that can be replicated for other comparative studies. We estimate electricity production using a point absorber WEC modeled near Sitka, AK, USA and optimized for the given sea conditions. We determine that, for remote community sized projects, seasonal deployments could result in small cost savings (less than 10 %), but larger decreases in annual energy production (around 30 % for our case study area). Seasonal deployments could be preferable in places with seasonal energy needs, if failures and device access become a major hindrance to wave energy technology development, or as a cautionary approach to introducing new technology to the oceans. We also determine that a highly seasonal wave resource is not necessarily a requirement for seasonal deployments to be considered. Seasonal deployments are an alternative to year-round deployments that can be considered in places where marine spatial conflict is a seasonal concern.

Low-temperature soldered wire interconnection (LTSWI) is a technology utilizing many interconnect wires carried on a polymer foil to form electrical connections against cell gridlines without a separate soldering process. In this work, LTSWI module samples were characterized for material properties and assembly dimensions and subjected to accelerated aging experiments to induce degradation. A finite element analysis model was developed based on characterization results, to analyze internal stressors during environmental exposures. The polymer foil contains polyethylene terephthalate and low-density polyethylene layers, and solder composition was tin bismuth, which notably was not metallurgically bonded to cell gridlines. High temperature accelerated exposures created power loss up to 9% in minimodule samples, with fill factor losses implicating contact degradation. Posttest characterization identified solder-gridline cracking and wire-cell separation as contributing mechanisms. Finite element modeling demonstrated that wire-to-cell contact is maintained by polymer contraction post lamination but is reversible, resulting in contact loss and wire separation during high temperature exposure. Simulations also detected in-plane wire-to-cell displacements, driven by surrounding polymer motion in response to high temperatures and mechanical load. We hypothesize that the propensity for wire movement during environmental exposure damages the not-metallurgically bonded wire-gridline interface and contributes to LTSWI contact degradation. Because distinct from thermal expansion mismatches which damage traditionally soldered modules, current test protocols are likely not applying the intended acceleration factors to LTSWI modules. This work highlights how construction-specific accelerated testing may be needed for nontraditional module designs and provides a starting point for accurate LTSWI life assessment.

Lost circulation material (LCM) selection is critical to effectively and efficiently treating wellbore fluid losses in geothermal drilling where costs of treatment can be as much as 30% of the total drilling cost. We conducted several uniaxial compaction experiments on 10 different materials and several material mixtures to identify critical mechanical parameters of each. Materials degraded at 200°C were also investigated to understand how elevated temperatures in geothermal wells would degrade their compaction behavior. Granular materials tended to have lower compressibility and higher compression resistance, while more elongated and softer materials had less mechanical stiffness. Mixing materials tended to moderate the mechanical behaviors while heating universally increased the compaction of materials. Microscopy showed that particle strength tended to correlate positively with roundness and circularity and negatively with elongation of a material. Convexity of the degraded and undegraded materials showed heating may have increased the convexity or roughness of the individual particles. We concluded that granular materials are likely to provide the best seals in wells but that a mixture of size distribution, mechanical rigidity, and elongation is more likely to form a better seal for geothermal wells.

International Journal of Hydrogen Energy

Additive manufacturing (AM) processes, like 3D printing, help to facilitate complex and customizable battery geometries which can provide design freedom and enhance volumetric energy density within electronic devices. AM materials must have the thermal and mechanical properties that enable printability, and when used in batteries, AM materials must also be chemically and electrochemically compatible with the battery chemistry. The compatibility between AM materials and the battery is of particular importance for the cell packaging materials which must be inert and are often overlooked. This study systematically studies AM-compatible polymeric materials for use as gaskets in lithium-ion cells. The materials investigated include three thermoplastics suitable for material extrusion printing: polylactic acid (PLA), polycarbonate, and polypropylene/polyethylene copolymer (PPPEC); and two photoresins suitable for vat photopolymerization (VPP) printing: an acrylate-based photoresin and a polyethylene glycol diacrylate photoresin. The AM gasket materials were tested in comparison to a conventional commercial polypropylene gasket. Mechanical testing (swell measurements and material stiffness) and electrochemical testing (linear sweep voltammetry and galvanostatic cycling of full cells) demonstrated that PLA and the VPP polymers were the least compatible with the lithium-ion battery chemistry, despite their prevalent use in studies of AM batteries, and that PPPEC was the most compatible.

Abstract not provided.

Current methods for modeling non-adiabatic molecular dynamics face fundamental limitations when treating geometric phase effects, quantum mechanical phenomena where nuclear wavepackets acquire phase shifts when encircling conical intersections. Existing approaches either neglect these effects entirely or rely on potential energy surfaces arising from the Born-Oppenheimer approximation, which introduce artificial singularities and can overestimate geometric phase contributions. This project deve

Fluid-fluid interfacial instability and subsequent fluid mixing are ubiquitous in nature and engineering. The hydrodynamic instability of fluid interfaces has long centered on the pressure gradient-driven long-wavelength Rayleigh-Taylor instability and the resonance-induced short-wavelength Faraday instability. However, neither instability alone can explain the dynamics when both mechanisms are present. We identify a previously unseen multi-modal instability emerging from their coexistence. When the denser fluid is polydimethylsiloxane, the mixed region at a high density contrast (Atwood number = 0.9) spans a vibration amplitude range approximately twice the gravitational acceleration. Using Floquet stability analysis, we show how vibrations govern transitions between the RT and Faraday instabilities, leading to contention between these instabilities rather than resonant enhancement. The initial transient growth is represented by the exponential modal growth of the most unstable Floquet exponent, along with its accompanying periodic behavior. Direct numerical simulations validate these findings and track interface breakup into the multiscale and nonlinear regimes. Specifically, we show that growing RT modes nonlinearly suppresses Faraday responses even when the initial growth rate of the Faraday instability is 3.63 times that of RT, so a bidirectional competition hinders their sustained coexistence.

While first-principles density functional theory modeling has become a vital tool to investigate defect properties in semiconductors, the lack of crystalline periodicity in pseudobinary random composition alloys, such as In1−xGaxAs, complicates such analyses. We present a simulation strategy to systematically take into account the variability in the local defect environment in order to predict statistical properties of neutral intrinsic defects in In1−xGaxAs. We use a comprehensive sampling from a modest-sized 64-atom special quasirandom structure (SQS) to define a statistically representative set of defects, and use a 512-atom hypercell, a 2×2×2 supercell of SQS supercells, to achieve cell-size convergence. We articulate an equivalent site principle and describe how it constrains atomic chemical reference energies in computation of defect formation energies in pseudobinary alloys. A simple protocol for estimating reference energies for the Ga and In atoms sharing the III site succeeds in obtaining the equivalence of defects at Ga-sites and In sites in the SQS supercell, (<30 meV differences in average formation energies). For III-site defects, such as the As antisite AsIII, the statistical variability in formation energies is modest, ≈0.1–0.2 eV. The variability in formation energy at As-site defects, such as the As vacancy vAs, can be much larger, >1 eV. The As antisite is shown to be a low-energy defect and the most likely to be present in as-grown materials, just as in GaAs. All other defects are higher-energy defects unlikely to be important in native material, but potentially important in radiation-damaged material. With a strong variability in defect energies, especially on the As-site, explicit consideration of statistical variability due to compositional randomness will be imperative for meaningful and quantitative comparisons to experiment.

The Wiedemann-Franz (WF) law correlates heat and charge transport in metals. However, the validity of this correlation remains an open-ended question, especially in the context of inelastic scattering at room temperature. To address this gap in knowledge, we perform independent measurements of the in-plane thermal and electrical conductivities across four AlCu (0.5% Cu) films [thickness ( h ) ≈ 174, 98, 53, and 24 nm] using optical pump-probe metrologies and four-point probe techniques, respectively. For in-plane thermal conductivity measurements, we utilize time-domain thermoreflectance, in both concentric and beam-offset configurations, and the time-resolved magneto-optic Kerr effect. Our results show that the WF law overpredicts the thermal conductivity by at least ∼10% in all films, thus demonstrating modest deviations in predicted thermal conductivity when applying the WF law to dilute AlCu films. Using infrared variable angle spectroscopic ellipsometry, we demonstrate increased electron scattering rates in the thinnest film ( h ≈ 24 nm), indicating electron-boundary scattering drives the reduction in in-plane thermal conductivity. This is generally an elastic scattering process, which is supported by our thermal conductivity measurements and analysis.

This paper studies Tikhonov regularization (ridge regression) parameter selection for problems in vibrations and acoustics. The selection method is based on a popular Bayesian method, but it incorporates measurements of sensor noise. The regularization parameter is closely related to the ratio of system input energy to noise energy, so noise measurements inform the inference procedure and improve parameter identification. In cases where standard Bayesian regularization identifies zero as the optimal regularization parameter, noise measurements guarantee a unique nonzero optimum. Sufficient theoretical criteria are developed for this guarantee. The method is verified in even-determined and under-determined configurations in an acoustic source localization simulation and a vibration load identification experiment. It is shown to yield significant improvements over existing empirical Bayesian regularization. Improvements are larger in the even-determined case and smaller in the under-determined case, wherein the inverse solution is less sensitive to the regularization parameter.

Abstract not provided.

The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysics scenarios. The code suite is composed of several specialized applications which can operate either in standalone mode or coupled with each other. Arpeggio is a supported utility that enables loose coupling of the various Sierra Mechanics applications by providing access to Framework services that facilitate the coupling.

memo on outgassing of via fill material

Nanoguide optical waveguides are thermally processed materials based on polymers that have different refractive indices. Scintillating versions of these materials have been under development since 2019 and feature PMMA as the low refractive index constituent and Organic Glass Scintillator-polymer blends as the high index component. Characterization efforts to date have focused on the functional characteristics of nanoguide pertaining to its envisioned use in high-resolution transmission radiogra

Abstract not provided.

The National Nuclear Security Administration (NNSA) seeks to assist and support all partners in the fields of radiobiology, health physics, radiation physics, and related areas to transition from cesium-137 chloride-based technologies that can potentially be used in an act of terrorism to X-ray technologies. The absence of information on experimental procedures, equipment, and irradiation parameters directly and negatively impacts the reproducibility and translatability of a sizable portion of radiation biology studies. It also discourages researchers from adopting new tools that eliminate the risks of a radiological dispersal device.

Abstract not provided.

Journal of Marine Science and Engineering

Abstract not provided.

Abstract not provided.

GPS technology is used in a multitude of applications around the world for navigation. However, GPS can be contested or denied and thus alternatives to GPS are needed for positioning in high-consequence systems. This project investigates a non-RF based method to produce an estimated geolocation using stellar/celestial measurements in combination with measurements from inclinometers and a real time clock.

Final report for the swept lookback Exploratory Express project.

Abstract not provided.

We present the Materials Learning Algorithms (MALA) package, a scalable machine learning framework designed to accelerate density functional theory (DFT) calculations suitable for large-scale atomistic simulations. Using local descriptors of the atomic environment, MALA models efficiently predict key electronic observables, including local density of states, electronic density, density of states, and total energy. The package integrates data sampling, model training and scalable inference into a unified library, while ensuring compatibility with standard DFT and molecular dynamics codes. We demonstrate MALA's capabilities with examples including boron clusters, aluminum across its solid-liquid phase boundary, and predicting the electronic structure of a stacking fault in a large beryllium slab. Scaling analyses reveal MALA's computational efficiency and identify bottlenecks for future optimization. With its ability to model electronic structures at scales far beyond standard DFT, MALA is well suited for modeling complex material systems, making it a versatile tool for advanced materials research.

Abstract not provided.

Abstract not provided.

Presented in this document is a small portion of the tests that exist in the Sierra/SolidMechanics (Sierra/SM) verification test suite. Most of these tests are run nightly with the Sierra/SM code suite, and the results of the test are checked versus the correct analytical result. For each of the tests presented in this document, the test setup, a description of the analytic solution, and comparison of the Sierra/SM code results to the analytic solution is provided. Mesh convergence is also ...

This is the level 3 milestone report (M3) for the workpackage SF-25SN01030407 in the Geologic Disposal Safety Assessment (GDSA) program.

High-fidelity large-eddy simulations (LES) and wall-modeled large-eddy simulations (WMLES) are performed on the interaction of a Mach 9 turbulent boundary layer with a 34-degree compression ramp. Our primary interest is the sensitivity of mean and unsteady wall quantities to grid density. Additionally within the context of WMLES, we seek to explore solution sensitivity to the wall-model exchange location.

The report covers the RS105 Wire Array Pulser testbed. The setup of the wire array at Reapplication, building 996, is laid out. Setup instructions are provided. Details of the commissioning test are given. Measurement data from the commissioning test shows the field output of the array exceeds 50 kV/m and that tuning the pulser output with a peaking switch can provide a rise time of ~2.3 ns and a pulse width of ~23 ns. The field levels and timeframes shown are within the tolerance conditions of

Abstract not provided.

LDRD Small project report

The influence of the internal structure at micrometer length scales on the deformation of polycrystalline materials can be effectively captured using crystal plasticity finite element methods (CPFEM). However, the complexity and nonlinearity of the deformation equations CPFEM solves demand significant computational power and resources to achieve accurate predictions, limiting its broader application. To address this challenge, we have identified a reduced-order representation of the complex data in order to establish a computationally efficient reduced-order models (ROM) and drastically reduce the computational expense of CPFEM. Specifically, in this work, we developed a parametric, data-driven, and non-intrusive ROM framework for CPFEM using proper orthogonal decomposition (POD) and sparse variational Gaussian process (SVGP) regression for single-crystal microstructures under tensile loading conditions. The developed protocol enables one to compress field into a latent/low-dimensional space described by principal component analysis (PCA) via the singular value decomposition (SVD) algorithm. As a result, the high-dimensional data are reduced to a significantly smaller amount of dimensions with POD bases and POD coefficients. Furthermore, we deployed an ensemble of SVGPs—extended from the classical Gaussian process (GP) regression for scalability and handling big data—in a massively parallel manner to train and predict latent POD coefficients using known POD bases from a set of previously obtained simulations results. Lastly, using the predicted POD coefficients, we reconstructed the full-field results and showed reasonable agreement compared with the true values obtained from running CPFEM. The developed framework is validated with a set of CPFEM simulations of a single embedded void in single-crystal aluminum alloy. While the framework is broadly applicable, this work specifically focuses on single-crystal microstructures, a single load case (e.g., tensile), and a specific void geometry (spherical).

Presented in this document is a portion of the tests that exist in the Sierra Thermal/Fluids verification test suite. Each of these tests is run nightly with the Sierra/TF code suite and the results of the test checked under mesh refinement against the correct analytic result. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the code results to the analytic solution is provided.

AGU Advances

Abstract not provided.

This report reviews regulations, codes, and standards to be considered for FCEB operations on airports governing vehicle safety, hydrogen infrastructure, and airport operations. Overall, the review found the existing regulatory framework largely suitable for FCEB operations on airports, with only minor gaps identified pertaining to heavy-duty vehicle fueling, operation and fire safety on airport aprons, and authorization for emergency power use.

Abstract not provided.

Accurate identification of material parameters is crucial for predictive modeling in computational mechanics. The two primary approaches in the experimental mechanics community for calibration from full-field digital image correlation data are known as finite element model updating (FEMU) and the virtual fields method (VFM). In VFM, the objective function is a squared mismatch between internal and external virtual work or power. In FEMU, the objective function quantifies the weighted mismatch between model predictions and corresponding experimentally measured quantities of interest. It is minimized by iteratively updating the parameters of an FE model. While FEMU is seen as more flexible, VFM is commonly used instead of FEMU due to its considerably greater computational expense. However, comparisons between the two methods usually involve approximations of gradients or sensitivities with finite difference schemes, thereby making direct assessments difficult. Hence, in this study, we compare VFM and FEMU in the context of numerically-exact sensitivities obtained through local sensitivity analyses and the application of automatic differentiation software. To this end, we conduct a series of test cases to assess both methods under practical challenges using a finite strain elastoplasticity model.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A series of experiments where deuterium was released in trace amounts into a room with a fire were performed at Sandia. This report describes the corresponding effort to model the test series using SIERRA/FUEGO. The objective of this modeling effort was to produce a simulation test matrix that can be utilized to help interpret the corresponding experiments and be used to assess the credibility of using the SIERRA/FUEGO simulations as a surrogate for real tritium reaction data.

This report describes challenges associated with the hierarchical Bayesian approach to inform model-form uncertainty (MFU) representations, which are parameterized modifications to a mathematical models' governing equations to express uncertainty in form of the equations.

IEEE Transactions on Quantum Engineering

Abstract not provided.

This report evaluates leakage behavior from tritium containment vessels under thermal abuse and combined thermal-mechanical abuse conditions to better understand safety implications for releases occurring in a fire scenario. Surrogate gases were used for all tests in this report. Leakage through the valves from thermal pressurization was observed when heating rates >0.8°C/s were sustained to >300°C. Gas plumes were visualized from vessels that were heated above 260°C and then dropped.

We present a curated, openly accessible dataset of 71 regional meteor events simultaneously recorded by optical and infrasound instrumentation between 2006 and 2011. These events were captured during an observational campaign using the all-sky cameras of the Southern Ontario Meteor Network and the co-located Elginfield Infrasound Array. Each entry provides optical trajectory measurements, infrasound waveforms, and atmospheric specification profiles. The integration of optical and acoustic data enables robust linkage between observed acoustic signals and specific points along meteor trajectories, offering new opportunities to examine shock wave generation, propagation, and energy deposition processes. This release fills a critical observational gap by providing the first validated, openly accessible archive of simultaneous optical–infrasound meteor observations that supports trajectory reconstruction, acoustic propagation modeling, and energy deposition analyses. By making these data openly available in a structured format, this work establishes a durable reference resource that advances reproducibility, fosters cross-disciplinary research, and underpins future developments in meteor physics, atmospheric acoustics, and planetary defense. Dataset: https://doi.org/10.5281/zenodo.15868512. Dataset License: CC-BY-NC

Abstract not provided.

Coherent anti-Stokes Raman scattering (CARS) and nitric oxide molecular tagging velocimetry (NO-MTV) were used to characterize the freestream in Sandia’s Hypersonic Shock Tunnel (HST) using a burst-mode laser operated at 100 kHz. Experiments were performed at nominal freestream velocities of 3 and 4 km/s using both air and N2 test gas. The CARS diagnostic provides nonequilibrium characterization of the flow by measuring vibrational and rotational temperatures of N2 and O2, which are compared to NO temperatures from separate laser absorption experiments. Simultaneous, collinear freestream velocities were measured using NO MTV along with pitot pressure measurements. This extensive freestream dataset is compared to nonequilibrium CFD that is capable of modeling species-specific, vibrational temperatures throughout the nozzle expansion. Significant nonequilibrium between vibrational and rotational temperatures was measured at each flow condition. N2 exhibits the most nonequilibrium, followed by O2 and NO. The CFD model captures this trend, although it consistently overpredicts N2 and O2 vibrational temperatures. At 3 km/s, the modeled NO nonequilibrium is underpredicted, whereas it is overpredicted at 4 km/s. Good agreement is seen between CFD and the velocity and rotational temperature measurements. Experiments with water added to the test gas yielded no discernable difference in vibrational relaxation.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The sputter deposition of alternating layers of Ni(V) and Al forms a reactive multilayer known to undergo self-propagating formation reactions when ignited. The sequential deposition process leads to nanometer-scale premixing of reactants at each included interface, which ultimately affects multilayer exothermicity. This work performs the direct measurement of a disordered face-centered cubic (FCC) solid solution premixed phase at the interfaces of Ni(V)/Al multilayers via scanning transmission electron microscopy. The crystallinity of the observed phase differs from previously reported a priori predictions of an amorphous interlayer. The disordered FCC phase retains its symmetry after annealing for 16 h at 135 ± 5 °C, but the lattice parameter shifts consistently with an Al-rich composition. The existence of a crystalline premix in Ni(V)/Al is attributed to the electronic contribution to the entropy of crystallization. The importance of electronic entropy to the phase formation of energetic materials motivates its inclusion when constructing digital twins for atomistic kinetics and ignition sensitivity.

Different numerical modeling methods have been developed and applied to evaluate a variety of performance indicators of wave energy converters (WECs), including the power performance, structural loads, levelized cost of energy, etc. Based on the modeling fidelity, the commonly used numerical modeling approaches can be classified as linear modeling, weakly nonlinear modeling and fully nonlinear modeling approaches. Each method differs in accuracy and computational efficiency, making them suitable for different stages of WEC design. However, the selection of modeling approach could significantly impact evaluation outcomes. For instance, simplified linear models may underestimate structural loads or overestimate energy production in some operational conditions, potentially leading to less cost-effective designs. Given the widespread utilization of these models, it is essential to understand the uncertainties brought by them in performance evaluations. This work is dedicated to benchmarking different linear-potential-flow-based numerical models for evaluating the systematic performance of WECs. Three representative numerical modeling approaches are considered in this work, including linear frequency-domain modeling, statistically linearized spectral-domain modeling and Cummins equation-based nonlinear time-domain modeling. A generic point absorber WEC is considered as the research reference in this work, and different sea sites are taken into account. The numerical models are utilized to predict critical performance indicators, including power performance, the annual energy production, the capacity factor, the levelized cost of energy and the PTO fatigue loads. By comparing the results, this work identifies the uncertainties associated with different modeling approaches in evaluating WEC performance.

The abstract is too long to fit inside this text box.

This report summarizes accomplishments in the development and maintenance of modeling and simulation tools to support material accountancy of bulk nuclear facilities. In FY25, we added new capabilities to MAPIT (new statistical test, improved error handling), launched the open source F3M modeling library, and added three new facility models to the SSPM-L model library.

Abstract not provided.

In previous years, SGTech funded enhancements to the isotopics routine in the software called Gamma Detector Response and Analysis Software-Detector Response Function (GADRAS-DRF), including the addition of peak fit customization capabilities [1][4][3]. A project was also funded that focused on implementing a peak-based model fitting routine, allowing model fitting to be performed without dependence on export-controlled cross-sections. In FY25 significant improvements were made to the custom pea

Abstract not provided.

The re-entry of the Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer (OSIRIS-REx) sample return capsule (SRC) on 24 September 2023 presented a rare opportunity to study atmospheric entry dynamics through a dense network of ground-based infrasound sensors. As the first interplanetary capsule to re-enter over the United States since Stardust in 2006, this event allowed for unprecedented observations of infrasound signals generated during hypersonic descent. We deployed 39 single-sensor stations across Nevada and Utah, strategically distributed to capture signals from distinct trajectory points. Infrasound data were analyzed to examine how signal amplitude and period vary with altitude and propagation path for a nonablating hypersonic object with welldefined physical and aerodynamic properties. Raytracing simulations incorporated atmospheric specifications from the ground-2-space model to estimate source altitudes for observed signals. Results confirmed ballistic arrivals at all stations, with source altitudes ranging from 44 to 62 km along the trajectory. Signal period and amplitude exhibited strong dependence on source altitude, with higher altitudes corresponding to lower amplitudes, longer periods, and reduced high-frequency content. Regression analysis demonstrated strong correlations between signal characteristics and both altitude and propagation geometry. Our results suggest, when attenuation is considered, the amplitude is primarily determined by the source, with the propagation path playing a secondary role over the distances examined. These findings emphasize the utility of controlled SRC re-entries for advancing our understanding of natural meteoroid dynamics, refining atmospheric entry models, and improving methodologies for planetary defense. The OSIRIS-REx SRC campaign represents the most comprehensive infrasound study of a hypersonic re-entry to date, showcasing the potential of coordinated geophysical observational networks for high-energy atmospheric phenomena, including space debris re-entries.

The trusted inertial terrain-aided navigation (TITAN) algorithm leverages an airborne vertical synthetic aperture radar to measure the range to the closest ground points along several prescribed iso-Doppler contours. These TITAN minimum-range, prescribed-Doppler measurements are the result of a constrained nonlinear optimization problem whose optimization function and constraints both depend on the radar position and velocity. Owing to the complexity of this measurement definition, analysis of the TITAN algorithm is lacking in prior work. This publication offers such an analysis, making the following three contributions: (1) an analytical solution to the TITAN constrained optimization measurement problem, (2) a derivation of the TITAN measurement function Jacobian, and (3) a derivation of the Cramér–Rao lower bound on the estimated position and velocity error covariance. These three contributions are verified via Monte Carlo simulations over synthetic terrain, which further reveal two remarkable properties of the TITAN algorithm: (1) the along-track positioning errors tend to be smaller than the cross-track positioning errors, and (2) the cross-track positioning errors are independent of the terrain roughness.

The SNL-developed F3M and MAPIT tools have the capability to analyze MC&A approaches for nuclear facilities via facility process flow simulation and statistical tests. Improvements on the application of F3M and MAPIT in simulating a generic TRISO fuel fabrication facility were successfully completed. This modeling framework can support the U.S. DOE and industry stakeholders in developing MC&A approaches for fuel fabrication facilities via demonstration of regulatory compliance.

Abstract not provided.

This project was motivated primarily by the opportunity of turning photocatalysis into a cost efficient hydrogen generation method for domestic energy and chemical production, and by the need to provide detailed mechanisms of photocatalytic reactions of broad scientific and technological interest. Specifically, we focused on the development of a new technique for analysis of (photo)catalytic reaction mechanisms to guide photocatalyst design. Our approach was based on applying velocity map ima

Abstract not provided.

Abstract not provided.

​​This white paper describes ongoing work and portfolios at Sandia National Laboratories that could be leveraged in AI for electric grid applications. This document highlights several areas where Sandia has developed capabilities that can be used in future work. These areas are human factors, uncertainty quantification, explainability, and trust maturity frameworks.

Abstract not provided.

This project extends the state of the art in formal verification modeling with modules and automatically checkable data-sharing patterns such that component modules can retain their assurance case when composed within a larger system. For users, smaller models make reasoning easier and help to ensure they accurately reflect text specifications. For automated methods, smaller models give exponential benefits for verification algorithm execution time.

In recent years, experiments involving micofabricated surface ion traps have grown in complexity, and as this complexity grows, a common design has emerged in the form of quantum charge-coupled device architecture. This architecture, already utilized on multiple systems, supports multiple zones on a device for areas such as memory or computation. The shutting process between these zones is a process often seen to be minimized. An additional component to scalable surface trapped-ion experiments

Abstract not provided.

This research project investigates the fundamental mechanisms of silicon nitride (SiN) crystallization, aiming to enhance the understanding of this critical material in microelectronics manufacturing. Through a collaborative effort between Sandia National Laboratories, the University of Tennessee, and the University of Florida, we developed a comprehensive framework that integrates experimental techniques, atomistic modeling, meso-scale simulations, and an integrated multi-scale model to capture

The SIERRA Low Mach Module: Fuego, henceforth referred to as Fuego, is the key element of the ASC fire environment simulation project. The fire environment simulation project is directed at characterizing both open large-scale pool fires and building enclosure fires. Fuego represents the turbulent, buoyantly-driven incompressible flow, heat transfer, mass transfer, combustion, soot, and absorption coefficient model portion of the simulation software.

Electrical characterization of Al doped SrTiO3 nanoparticles by undergraduate intern Aaron Coll

Fireballs (bolides) are high-energy luminous phenomena produced when meteoroids and small asteroids enter Earth’s atmosphere at hypersonic speeds, often resulting in fragmentation or complete disintegration accompanied by significant energy release. The resulting bolide light curves capture temporal brightness variations as these objects traverse increasingly dense atmospheric layers, providing essential information on meteoroid entry dynamics, fragmentation behavior, and atmospheric energy deposition processes. The Center for Near-Earth Object Studies’ (CNEOS) continuously expanding fireball database offers a globally comprehensive archive of bolide events, including light curves and associated metadata. Events associated with infrasound detections allow direct correlations between acoustic signatures and light curve features, therefore enabling detailed analyses of fragmentation dynamics and energy deposition. Here, we introduce Bolide Light-curve Analysis and Discrimination Explorer (BLADE), a robust and high-fidelity framework specifically designed to analyze bolide light curves for objects detected from space. BLADE incorporates a processing pipeline integrating Savitzky-Golay filtering, prominence-based peak detection, and gradient analysis, enabling systematic identification and classification of fragmentation events and their associated energy release characteristics. Preliminary results demonstrate that BLADE reliably distinguishes distinct bolide behaviors, providing an objective, scalable methodology for characterization and analysis of large bolide light curve data sets. This foundational work establishes a novel pathway for advanced bolide research, with promising applications in planetary defense and global atmospheric monitoring. Future research should adopt an integrative approach combining CNEOS optical data with complementary infrasound measurements, further clarifying relationships between bolide energy deposition and acoustic signatures, thus refining our understanding of meteoroid and asteroid atmospheric entry processes.

The motivation for the experiments reported here pertains to the siting of Liquefied Natural Gas (LNG) facilities which requires assessing the potential adverse radiant thermal impacts of accidental fires on the public. The objective is to obtain data on jet fires, pool fires, fireballs, and concrete walls that could serve as thermal barriers for model validation. The fuels tested include ethane, ethylene, propane, and isopentane.

This document details the findings from the FY25 RAD-Tech LDRD titled “Radiation Effects on NoC (Network on Chips).” We utilized the Versal FPGA from AMD as an exemplar platform for NoC. We conducted two radiation tests, one at Texas A&M University (TAMU) in June 2025 and another at Lawrence Berkeley National Laboratory (LBNL) in August 2025. These experiments showed that radiation could upset the NoC and that it experiences a variety of failures, that we could detect those upsets,

Abstract not provided.

Ice-mass change induces regionally varying patterns of sea-level change due to gravitational, rotational, and deformational (GRD) effects, which in turn influence marine-based ice stability in Antarctica. For improved projection of the Antarctic Ice Sheet (AIS), there is a need for including GRD effects in modeling and improving understanding of basin-by-basin sensitivity of ice evolution to GRD effects under a range of climate scenarios. We couple a high-resolution, higher-order ice-sheet model with a 1D global sea-level model that fully captures GRD effects, and simulate ice evolution in Antarctica under the Ice Sheet Model Intercomparison Project for CMIP6 experiments. We perform two sets of coupled simulations incorporating 1D Maxwell solid Earth structure suitable for West and East Antarctica and show that the Amundsen Sea Embayment (ASE) in West Antarctica has the highest sensitivity to GRD effects—in high-emission scenarios, grounding-line retreat accelerates by hundreds of kilometers by 2300 without GRD effects, but GRD effects delay this retreat on a timescale of decades. However, we find that delay times do not show a clear relationship to the strength of climate forcing alone. Furthermore, GRD effects can influence ice-sheet dynamics more than the choice of climate model for a given emissions scenario. In contrast, East Antarctica exhibits minimal sensitivity to GRD effects throughout the study period. These findings underscore the critical role of GRD effects in shaping future West AIS evolution, highlighting the importance of constraining the regional 3D Earth structure and bed topography in West Antarctica, particularly the ASE.

Abstract not provided.

Prussian blue analogues (PBAs) have attracted increasing interest owing to their potential applications in various fields such as energy storage and conversion, neuromorphic computing, and magnetic switching. With a general formula of AxMN[MC(CN)6], they feature an open framework that provides abundant channels for diffusion of alkali metal ions A and allows flexible compositional control of transition metal ions MNand MC. The oxidation states of transition metal ions can be tuned by adjusting the amount (x) of alkali ions A. Here, we carried out density functional theory calculations combined with experimental measurements to investigate the effects of transition metal ions, alkali ions, and oxidation states on the electronic properties of PBAs. Our calculations found that the band gaps of PBAs can be tuned from close to 0 eV to more than 4 eV. Experimentally, we introduced the synthesis/characterization of five previously unreported PBAs (MN= Ru, Os; MC= Fe, Ru, and Os) to complete the nine stable MN:MCtransition metal combinations in group VIII of the periodic table. The optically measured intervalence charge transfer excitation energies of group VIII PBAs are consistent with calculated band gaps. They demonstrate wide band gap tunability by adjusting transition metals and oxidation states, enabling semiconductor-to-metal transitions for memristor applications and enhancing electronic conductivity for battery applications. This work provides a computational/experimental database of electronic properties versus structural compositions for PBAs.

As with other oscillatory power conversion systems, the design of wave energy converters can be understood as an impedance matching problem. By representing the wave energy converter as a multi-port network, two separate but related impedance matching conditions can be established. Satisfying these conditions maximizes power transfer to the load. In practice, these impedance matching conditions may be used to influence the design of the system (including the hull, power take-off, controller, mooring, etc.). To this end, this paper considers some example applications of wave energy converter design with the help of the impedance matching framework.

The June 1991 Mt. Pinatubo eruption resulted in a massive increase of sulfate aerosols in the atmosphere, absorbing radiation and leading to global changes in surface and stratospheric temperatures. A volcanic eruption of this magnitude serves as a natural analog for stratospheric aerosol injection, a proposed solar radiation modification method to combat a warming climate. The impacts of such an event are multifaceted and region-specific. Our goal is to characterize the multivariate and dynamic nature of the atmospheric impacts following the Mt. Pinatubo eruption. We developed a multivariate space-time dynamic linear model to understand the full extent of the spatially- and temporally-varying impacts. Specifically, spatial variation is modeled using a flexible set of basis functions for which the basis coefficients are allowed to vary in time through a vector autoregressive (VAR) structure. This novel model is cast in a Dynamic Linear Model (DLM) framework and estimated via a customized MCMC approach. We demonstrate how the model quantifies the relationships between key atmospheric parameters prior to and following the Mt. Pinatubo eruption with reanalysis data from MERRA-2 and highlight when such a model is advantageous over univariate models.

Topological insulators are characterized by insulating bulk states and robust metallic surface states. Band inversion is a hallmark of topological insulators. At time-reversal invariant points in the Brillouin zone, spin–orbit coupling (SOC) induces a swapping of orbital character at the bulk band edges. Reliably detecting band inversion in solid-state systems with many-body methods would aid in identifying possible candidates for spintronics and quantum computing applications and improve our understanding of the physics behind topologically nontrivial systems. Density functional theory (DFT) methods are a well-established means of investigating these interesting materials due to their favorable balance of computational cost and accuracy but often struggle to accurately model the electron–electron correlations present in the many materials containing heavier elements. In this work, we develop a novel method to detect band inversion within continuum quantum Monte Carlo (QMC) methods that can accurately treat the electron correlation and spin–orbit coupling that are crucial to the physics of topological insulators. Our approach applies a momentum-space-resolved atomic population analysis throughout the first Brillouin zone utilizing the Löwdin method and the one-body reduced density matrix produced with diffusion Monte Carlo (DMC). We integrate this method into QMCPACK, an open source ab initio QMC package, so that these ground-state methods can be used to complement experimental studies and validate prior DFT work on predicting the band structures of correlated topological insulators. We demonstrate this new technique on the topological insulator bismuth telluride, which displays band inversion between its Bi-p and Te-p states at the Γ-point. We show an increase in charge on the bismuth-p orbital and a decrease in charge on the tellurium-p orbital when comparing band structures with and without SOC. Additionally, we use our method to compare the degree of band inversion present in monolayer Bi2Te3, which has no interlayer van der Waals interactions, to that seen in the bilayer and bulk. The method presented here will enable future many-body studies of band inversion that can shed light on the delicate interplay between correlation and topology in correlated topological materials.

The June 1991 Mt. Pinatubo eruption resulted in a massive increase of sulfate aerosols in the atmosphere, absorbing radiation and leading to global changes in surface and stratospheric temperatures. A volcanic eruption of this magnitude serves as a natural analog for stratospheric aerosol injection, a proposed solar radiation modification method to combat a warming climate. The impacts of such an event are multifaceted and region-specific. Our goal is to characterize the multivariate and dynamic nature of the atmospheric impacts following the Mt. Pinatubo eruption. We developed a multivariate space-time dynamic linear model to understand the full extent of the spatially- and temporally-varying impacts. Specifically, spatial variation is modeled using a flexible set of basis functions for which the basis coefficients are allowed to vary in time through a vector autoregressive (VAR) structure. This novel model is cast in a Dynamic Linear Model (DLM) framework and estimated via a customized MCMC approach. We demonstrate how the model quantifies the relationships between key atmospheric parameters prior to and following the Mt. Pinatubo eruption with reanalysis data from MERRA-2 and highlight when such a model is advantageous over univariate models.

State-of-the-art noisy-intermediate-scale quantum processors are currently implemented across a variety of hardware platforms, each with their own distinct gatesets. As such, circuit compilation should not only be aware of but also deeply connect to the native gateset and noise properties of each. Trapped-ion processors are one such platform that provides a gateset that can be continuously parameterized across both one- and two-qubit gates. Here we use the Quantum Scientific Computing Open User Testbed to study noise-aware compilations focused on continuously parameterized two-qubit ZZ gates (based on the Mølmer-Sørensen interaction) using superstaq, a quantum software platform for hardware-aware circuit compiler optimizations. We discuss the realization of ZZ gates with arbitrary angle on the all-to-all connected trapped-ion system. Then we discuss a variety of different compiler optimizations that innately target these ZZ gates and their noise properties. These optimizations include moving from a restricted maximally entangling gateset to a continuously parameterized one, swap mirroring to further reduce the total entangling angle of the operations, focusing the heaviest ZZ angle participation on the best-performing gate pairs, and circuit approximation to remove the least impactful ZZ gates. We demonstrate these compilation approaches on the hardware with randomized quantum volume circuits, observing the potential to realize a larger quantum volume as a result of these optimizations. Using differing yet complementary analysis techniques, we observe the distinct improvements in system performance provided by these noise-aware compilations and study the role of stochastic and coherent error channels for each compilation choice.

This project advanced the development of an active hybrid mooring system integrating experimental testing with numerical simulation to capture complex mooring dynamics not feasible in existing wave basins, including deep-water and shared mooring interactions. Verification testing of the complete hybrid system yielded good agreement in mooring tension and platform translation responses, though discrepancies in platform pitch response indicate that further refinement may be needed in future work.

Abstract not provided.

This report provides quantitative estimates of the potential economic effects on storage, transport, and disposal of high assay low-enriched uranium spent fuel. It builds on a previous report that described these economic effects qualitatively.

This user guide supports a curated dataset of 71 meteor events recorded between 2006 and 2011 in Southwestern Ontario, Canada. Each event was simultaneously observed by ground-based optical cameras and an infrasound array, providing a rare opportunity to examine meteor trajectories and acoustic signals from the same atmospheric entry events.

ThisreportsupplementstheVerificationofEmpire-CableSANDreport[1]byexpandingontheuse ofXycetosimulatethecouplingtoatransmissionlinecablemodel. WhileEmpire-Cablesolvesits governing equations on a high-order, finite-element mesh with an an implicit-in-time formulation, Xyce must use a first order graph for the circuit and explicit-in-time approach to be compatible with non-linear electrical device models[2]. Thus, given the different solution methodologies in Xyce as compared to Empire-Cable, the co

Journal of Advances in Information Technology

The strength of materials is influenced by a range of external conditions, such as temperature and deformation rate. Consequently, materials that demonstrate substantial variations in their mechanical behavior due to fluctuations in temperature and strain rate require complex strength models to accurately predict material performance in real-world applications. To predict such complex behavior, a robust and flexible strength model is necessary. In this work, we utilize genetic programming-based symbolic regression (GPSR) to develop data-driven strength models that accurately represent the measured stress–strain responses of tin across a wide range of strain, strain rate and temperature regimes. The GPSR models are constrained by physically-informed conditions, which leads to significant improvement in extrapolation. The best model is integrated into a multi-physics code to perform Taylor impact simulations, validating the model's accuracy and robustness. The model predictions showed excellent agreement with experimental results, particularly when compared to predictions using traditional strength models.

Abstract not provided.

Abstract not provided.

This is a single slide that highlights MESA packaging and integration capabilities.

This report documents the work done as part of the “Enhancing Developer Productivity” level 2 milestone. The team surveyed developers about impediments and successes; improved our CI pipeline monitoring and reporting; developed tools for line coverage reporting and analysis; improved compiler warning adherence in SIERRA; prototyped static analysis, AI, and mutation testing tooling in SIERRA; and developed a 3-5 year SIERRA plan document to help these initiatives continue past this milestone.

Abstract not provided.

This lessons learned document provides a detailed analysis of the hydroxy termination titration methods employed for determining hydroxyl group content in various chemical samples. The primary focus is on the evaluation of techniques such as ASTM E222 and their practical applications in laboratory settings. Through examination of method selection, sample preparation, and titration execution, key insights have been identified that enhance the accuracy and reliability of results.

While amorphous materials are often approximated to have a statistically homogeneous atomic structure, they frequently exhibit localized structural heterogeneity that challenges simplified models. This study uses 4D scanning transmission electron microscopy to investigate the strain and structural modifications around gas bubbles in amorphous Bi2O3 induced by argon irradiation. We present a method for determining strain fields surrounding bubbles that can be used to measure the internal pressure of the gas. Compressive strain is observed around the cavities, with higher-order crystalline symmetries emerging near the cavity interfaces, suggesting paracrystalline ordering as a result of bubble coarsening. This ordering, along with a compressive strain gradient, indicates that gas bubbles induce significant localized changes in atomic packing. By analyzing strain fields with maximum compressive strains of 3%, we estimate a lower bound on the internal pressure of the bubbles at 2.5 GPa. These findings provide insight into the complex structural behavior of amorphous materials under stress, particularly in systems with gas inclusions, and offer new methods for probing the local atomic structure in disordered materials. Although considering structural heterogeneity in amorphous systems is non-trivial, these features have crucial impacts on material functionalities, such as mechanical strength, ionic conductivity, and electronic mobility.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Ocean observation buoys require relatively small amounts of power, yet traditionally necessitate costly resupply trips for battery replacement. With the offshore location of the buoys and small power requirements, wave energy may be an effective solution for providing consistent and reliable power to support the buoy instrumentation. The US National Science Foundation Ocean Observatories Initiative (OOI) includes arrays of point absorber-like buoy systems used for ocean observation that have been deployed at multiple locations including the Southern Mid-Atlantic Bight. A study is currently underway to design a pitch resonator wave energy converter to supplement existing renewable energy generation for powering observation instrumentation. This paper details field measurements from surface moorings of the OOI Coastal Pioneer Array, which informs the subsequent development of a numerical model for the moored observation system. The model is developed in Wave Energy Converter Simulator (WEC-Sim), which leverages the Simscape multibody solver within the MATLAB/Simulink framework and linear potential flow theory to simulate the hydrodynamic interactions and multibody dynamics in 6 degrees of freedom. Multiple tuning variables are considered to produce a model for the system that matches well with empirical data (about 8% error). The WEC-Sim model will serve as a platform for integrating the pitch resonator wave energy converter concept and deployment preparation (detailed design including power take-off and control systems, response evaluation, etc.).

Abstract not provided.

New concepts of symmetry related to topological order emerged from the discovery of the fractional quantum Hall effect and high-temperature superconductivity in strongly correlated electron systems. This led to the study of quantum materials-- materials exhibiting emergent quantum phenomena with no classical analogues. While these materials have engendered exciting basic materials science and physics, realizing novel devices is a key challenge in the field. The goal of this proposal is to harnes

Abstract not provided.

A common method of determining elastic material properties is utilizing an ultrasonic pulse receiver. This is a non-destructive test (NDT) causing no plastic deformation of the material. It requires a pulse generator, oscilloscope, and two transducers to measure the sound velocity of a material. Both longitudinal and shear sound velocities may be obtained by utilizing the appropriate transducer. Acoustic nondestructive testing methods are often used in manufacturing certification processes and f

Abstract not provided.

Many systems derive their electromagnetic radiation (EMR) environments from Mil-Std-464C or Mil-Std-464D, Electromagnetic Environmental Effects Requirements for Systems. This document is intended to provide additional clarity on the Mil-Std-464C and Mil-Std-464D EMR environments. This supplementary information may be used to guide tailoring of relevant environments or application to different systems with the additional consideration of any shielding that may reduce the external environments.

Neutrons are produced at the Sandia Ion Beam Laboratory (IBL) for testing effects of neutron irradiation on electronic devices. This report provides detailed information on the methods used for the production and dosimetry of neutrons by the D(D,n)3He, 7Li(p,n)7Be and 12C(D,n)13N nuclear reactions at the IBL.

In order to perform post mortem analysis on materials launched in gas guns and powder guns, various methods have been explored to soft catch these materials. The majority of the deformation profile must be from the initial impact during the experiment, not from the subsequent impact with the energy absorption material in the catch tank. This process becomes increasingly more difficult with increased velocities. This report explores a method to soft catch materials launched at velocities 3000 ms

Deliverables for Summer Undergraduate Laboratory Internship (SULI) Program

Abstract not provided.

Abstract not provided.

We present an improved TIGER solution technique aimed at accurately determining the equilibrium state for systems that can produce a diverse array of products, including both solid and liquid condensed species.

BLADE_LC_processor.py is an optional script within the BLADE (Bolide Light-curve Analysis and Discrimination Explorer) open-source software package.

Physical Review Materials

Abstract not provided.

Abstract not provided.

Recorded seismic data are generally contaminated by noise from different sources, which masks the signals of interest. In the seismology community, frequency filtering (FF) is the standard method for noise suppression. However, when the signal of interest and noise share the same frequency band, the latter cannot be filtered out without infringing on the former. We implemented a noise suppression approach based on the mathematical morphology theorem. The method involves compound operations of dilation and erosion using structuring elements of varying lengths and decomposes an input noisy waveform into several time functions with differing characteristics. The filtered waveform is constructed from the time functions using a self-adaptive reconstruction technique. Application to a data set of > 4700 local waveforms suggests that the implemented mathematical morphological filtering (MMF) approach is efficient for data with low signal-to-noise ratio (SNR) and significantly outperforms FF in that SNR range. For most of the dataset, FF, machine learning (ML) denoising, and continuous wavelet transform (CWT) thresholding result in higher SNR values compared with the MMF method. However, for ∼42% of the waveforms, MMF outperforms FF, and the SNR gain achieved with MMF is as large as ∼23 dB. Compared to ML denoising and CWT thresholding, this proportion drops to only ∼10%–14%. Our results suggests that in an operational setting, MMF cannot replace the other noise suppression methods; however, signal detection can be improved if MMF is used to supplement them in some scenarios. MMF could help detect signals in problematic low-SNR data, which are currently being missed particularly when using FF alone.

Photonic Doppler velocimetry (PDV) is an established technique for measuring the velocities of fast-moving surfaces in high-energy-density experiments. In the standard approach to PDV analysis, the short-time Fourier transform (STFT) is used to generate a spectrogram from which the velocity history of the target is inferred. The user chooses the form, duration, and separation of the window function. Here, we present a Bayesian approach to infer the velocity directly from the PDV oscilloscope trace, without using the spectrogram for analysis. This is clearly a difficult inference problem due to the highly periodic nature of the data, but we find that with carefully chosen prior distributions for the model parameters, we can accurately recover the injected velocity from synthetic data. We validate this method using PDV data collected at the STAR two-stage light gas gun at Sandia National Laboratories, recovering shock-front velocities in quartz that are consistent with those inferred using the STFT-based approach and are interpolated across regions of low signal-to-noise data. Although this method does not rely on the same user choices as the STFT, we caution that it can be prone to misspecification if the chosen model is not sufficient to capture the velocity behavior. Analysis using posterior predictive checks can be used to establish whether a better model is required, although more complex models come with additional computational cost, often taking more than several hours to converge when sampling the Bayesian posterior. We, therefore, recommend it be viewed as a complementary method to that of the STFT-based approach.

New concepts of symmetry related to topological order emerged from the discovery of the fractional quantum Hall effect and high-temperature superconductivity in strongly correlated electron systems. This led to the study of quantum materials-- materials exhibiting emergent quantum phenomena with no classical analogues. While these materials have engendered exciting basic materials science and physics, realizing novel devices is a key challenge in the field. The goal of this proposal is to harnes

Various applications—including brain-like computing and on-chip artificial vision—increasingly demand a combination of electronic and photonic techniques. However, integrating both approaches on a single chip is challenging, and solutions typically rely on disparate components with power-hungry signal conversions. Here we report electro-optical Mott neurons that combine visible light emission with electrical threshold switching, as well as neuron-like oscillations. The devices are based on thin films of sputtered niobium dioxide (NbO2), a Mott insulator–metal transition material, operating at room temperature and emitting light that peaks around 810 nm. Operando measurements reveal an electronic origin to the light emission: charge carrier relaxation initiated by high-field transport in the NbO2. Our devices combine electrical and optical functions within a single material, thereby expanding the options available for future artificial intelligence hardware.

This document provides a step-by-step tutorial on how to set up OpenCSP for both users and developers. It is meant to support novice users and does not require software engineering experience. It is a detailed extension of the OpenCSP getting started on-line documentation, found here: https://opencsp.readthedocs.io/en/main/contributing.html#getting-started For an overview of OpenCSP overall, see the OpenCSP website: https://opencsp.sandia.gov. For details on OpenCSP components and algorithm

Stilbenes are a class of organic compounds with broad-ranging pharmaceutical and agricultural applications, which are typically isolated and purified through recrystallization. We are motivated by reducing experimental waste and optimizing yield via developing predictive simulations for processing-dependent crystal morphologies. Using resveratrol as a model stilbene system, we have developed an approach for simulating crystallization with molecular resolution using on-lattice kinetic Monte Carlo. In this work, we highlight modifications to the Stochastic Parallel PARticle Kinetic Simulator (SPPARKS) software package, which were essential to this application. Key enhancements include the incorporation of non-orthogonal cell shapes and monomer anisotropy approximations using bound hard spheres. This new SPPARKS application has been applied to resveratrol with attachment energy libraries obtained from density functional theory, resulting in excellent agreement with experimental morphology prediction.

Sapphire (Al₂ ⁢O₃), known for its remarkable incompressibility at ambient conditions, plays a pivotal role in both static and dynamic compression research. Accurately characterizing its equation of state (EoS) is essential for these applications. Here, we present a complete Hugoniot of Al₂⁢ O₃ as locus of experimentally assessed, high-precision, pressure, density and temperature states up to 14 Mbar and 43 kK. The Hugoniot is established with single shock experiments using magnetically launched hyper velocity flyers on the Z Accelerator at Sandia National Laboratories. We explore principal Hugoniot states at very high shock 𝑇 and 𝑝 in the solid phase, tracking the solid-liquid boundary and culminating at 2.4-fold compression, where data provides a direct constraint on the liquid phase. Corresponding shock release data probe thermodynamic states complementary to the Hugoniot and place additional constraints on tabular EoS models. Our findings indicate a significant deviation from existing tabular EoS models for Al₂⁢ O₃ dictating a comprehensive overhaul. We develop two advanced EoSs for Al₂⁢ O₃ the SESAME 97412 model, featuring an extensive phase diagram that includes three solid phases and the liquid phase, and the updated LEOS 2200m2 model. EoS development is assisted with Quantum Molecular Dynamics simulations. Our experimental data allows for stringent testing of our EoSs. Both models accurately capture the Hugoniot of Al₂ ⁢O₃ up to the highest pressures and temperatures. Rigorous experimental determination of extreme pressures and temperatures, paired with sophisticated models, advances the frontier of EoS development beyond 1 terapascal.

Coulomb drag is a powerful tool to study interactions in coupled low-dimensional systems. Historically, Coulomb drag has been attributed to a frictional force arising from momentum transfer whose direction is dictated by the current flow. In the absence of electron-electron correlations, treating the Coulomb drag circuit as a rectifier of noise fluctuations yields similar conclusions about the reciprocal nature of Coulomb drag. In contrast, recent findings in one-dimensional systems have identified a nonreciprocal contribution to Coulomb drag that is independent of the current flow direction. In this work, we present Coulomb drag measurements between vertically coupled GaAs/AlGaAs quantum wires separated vertically by a hard barrier only 15 nm wide, where both reciprocal and nonreciprocal contributions to the drag signal are observed simultaneously, and whose relative magnitudes are temperature and gate tunable. Our study opens up the possibility of studying the physical mechanisms behind the onset of both Coulomb drag contributions simultaneously in a single device, ultimately leading to a better understanding of Luttinger liquids in multi-channel wires and paving the way for the creation of energy harvesting devices.

Controlling the island decoration on metal nanoparticle supports is a major opportunity for improving the catalytic activity and an attractive synthetic challenge. The structure of the decorating metal determines how it interacts with the metal support and how it effectively catalyzes the reactants and the intermediates. In this work, we demonstrate that a slow-growth method maximizes the formation of Ru islands on faceted, branched Ni nanoparticles, thereby controlling the number of Ru–Ni atomic interactions and improving the catalytic activity. The Ru islands on branched Ni nanoparticles with the highest loading of Ru (9%) exhibited the highest activity for the electro-oxidation of biomass-derived 5-hydroxymethylfurfural (HMF). In conclusion, these results demonstrate the ability to synthetically control the second metal decoration to tune metal–support interactions, thereby enhancing the catalytic activity.

Coulomb drag is a powerful tool to study interactions in coupled low-dimensional systems. Historically, Coulomb drag has been attributed to a frictional force arising from momentum transfer whose direction is dictated by the current flow. In the absence of electron-electron correlations, treating the Coulomb drag circuit as a rectifier of noise fluctuations yields similar conclusions about the reciprocal nature of Coulomb drag. In contrast, recent findings in one-dimensional systems have identified a nonreciprocal contribution to Coulomb drag that is independent of the current flow direction. In this work, we present Coulomb drag measurements between vertically coupled GaAs/AlGaAs quantum wires separated vertically by a hard barrier only 15 nm wide, where both reciprocal and nonreciprocal contributions to the drag signal are observed simultaneously, and whose relative magnitudes are temperature and gate tunable. Our study opens up the possibility of studying the physical mechanisms behind the onset of both Coulomb drag contributions simultaneously in a single device, ultimately leading to a better understanding of Luttinger liquids in multi-channel wires and paving the way for the creation of energy harvesting devices.

The adoption of machine learning techniques in the seismology community has led to great performance improvements in several areas, including signal processing. Specifically, the development of deep learning–based seismic waveform denoising models has the potential to yield improvements in signal detection capabilities for networks operating in particularly noisy environments. Recent advancements in the design of these deep learning denoising models have included the incorporation of continuous and discrete wavelet transform functions into the network architecture to improve the learning capabilities and efficiency of said models. These wavelet transform–based seismic denoising models have shown improved denoising capabilities in regions where there is good agreement between the data features present in the training and evaluation datasets. However, questions remain about the overall transportability of these models to other monitoring regions. Here, in this study, we will determine the baseline transportability of a newly developed multilevel wavelet‐transform convolutional neural network (MWCNN) seismic denoising model. We accomplish this by taking a version of the MWCNN denoising model trained on data collected from the Utah region and evaluating its denoising performance on datasets collected from the neighboring Nevada region, which differ with regard to monitoring sensor types and event histories. We find that there is a notable variability in denoising performance related to the degree of similarity between the initial and new target datasets. The most notable difference in denoising performance is the ability of the denoising model to preserve accurate amplitude information associated with the signal energy present in the waveform data. Finally, we evaluate the ability of transfer learning techniques to improve the transportability of the MWCNN denoising model. We find that although there is still a performance gap present in the denoising results of the MWCNN model, transfer learning did yield improved results.

The Stretched-Wire Alignment Technique (SWAT) is one method of magnet alignment for linear induction accelerators. The applications of SWAT have been implemented for aligning solenoid magnets on the Scorpius linear induction accelerator which will be sited at the Nevada National Security Site and the Flash X-Ray (FXR) linear induction accelerator at Lawrence Livermore National Laboratory’s Contained Firing Facility. This article describes both systematic (repeatable) and random sources of background and noise as well as practical ways to eliminate or reduce them to acceptable levels. Systematic sources include reflections from wire ends, rapid sag due to ohmic heating of the wire, magnetic materials, and shot rate. Random sources include air currents, vibration of nearby equipment, mechanical stability of test equipment, and the instruments used to measure the wire motion. Mitigations include curve fitting and adaptive noise signal cancellation, and mechanical damping. Finite Element Analysis (FEA) was used to identify and resolve a repeatable wire vibration frequency interfering with the signal resolution. Two stretched wire alignment technique set ups from Sandia National Labs and Lawrence Livermore National Lab have shown background noise sources and ways of mitigating them by either analysis methods or change of mechanical configuration. Conclusions that were drawn included the severe sensitivity of the deflection to even small external interferences of the SWAT wire such that it requires attention to detail in mechanical set up and analysis.

Cost-effective and reliable hydrogen compression remains a challenging barrier in the widespread adoption of hydrogen as an energy carrier. The prevailing technology of mechanical compression suffers from several drawbacks, some of which can be addressed by nonmechanical compression strategies (e.g., electrochemical or metal hydride-based thermal compression). Thermally driven metal hydride compression strategies typically rely on multistage metal hydride-based compressors; however, discovering or optimizing low-stability metal hydrides that can pressurize hydrogen upward of 1000 bar is difficult, both with respect to computational predictions and experimental validation. Here, we (1) demonstrate that simple machine learning-derived design rules can inform the rational design of alloying strategies yielding low-stability hydrides, (2) validate their experimental pressure–composition–temperature (PCT) isotherms up to 875 bar, and (3) utilize a dynamic system-level model of a metal hydride compressor design to evaluate their performance under realistic operating conditions. Importantly, this analysis yields predicted operational efficiencies of both 2-stage (90–875 bar) and 3-stage (20–875 bar) metal hydride compressors to enable further evaluation of this technology and its techno-economic outlook.

Diazonium compounds are synthetically useful in the production of dyes and textiles, however they are highly explosive under dry conditions. Explosion prevention becomes more difficult when new diazonium compounds are synthesized, because while some syntheses include a counterion to increase their stability, this is not always a reliable method to prevent an explosive incident. Due to the uncertainty surrounding the explosiveness of different diazonium compounds, it is important to understand how to safely clean up after an incident and how to determine when it is safe to return a laboratory to typical operational use, particularly when the incident involves a novel compound where a standard does not exist for instrument calibration. Here, an explosive event is discussed involving the synthesis of 4-bromo-benzenediazonium-2-carboxylate. Following the explosive incident and 3-step cleanup, which involved a precautionary neutralization step, samples were collected from the fume hood where the incident occurred. Because the incident involved an unstable, novel compound that is not commercially available and was deemed unsafe to resynthesize for instrument calibration, we assessed the risk of further explosion by analyzing for the stable decomposition products. Mass spectrometry analysis confirmed that the residue in the fume hood contained 5-bromosalicylic acid, a decomposition product of 4-bromo-benzenediazonium-2-carboxylate. Samples were taken from multiple points in the fume hood and analyzed to estimate the spatial distribution of the decomposition product. Based on this analysis, we inferred that the primary decomposition product was far more abundant than residual energetic, indicating the energetic had been consumed or neutralized to a trace quantity where the risk of further explosion was low. The steps presented here─specifically, initial neutralization and then analyzing the spatial distribution of expected decomposition products to assess risk when a novel explosive material is detonated in a confined space─were our approach to assess further risk following an explosion due to a novel diazonium compound without the need for any further handling or resynthesis of the energetic. Here, we present our approach and critically analyze these steps by discussing retrospective lessons learned and alternative analytical approaches.

Ion implantation emerges as a promising technique to address the persistent challenge of lithium (Li) filament growth in solid-state electrolytes as it can induce compressive stresses inhibiting crack growth and deflect dendrites, de facto mitigating early electrolyte failure. In this study, we examine the potential paradox of ion implantation: while aiming to enhance electrolyte performance, the radiation damage associated with implantation might inadvertently compromise both the ionic conductivity and the intrinsic fracture toughness of the material, rendering the material unsuitable for battery applications. Specifically, we employed molecular dynamics simulations to examine the scope of the downsides of ion implantation, specifically: (i) reduced ionic conductivity (due to radiation-induced amorphization) and (ii) mechanical stability (due to radiation-induced embrittlement) in ion-implanted Li₇La₃Zr₂O₁₂ (LLZO) solid-state electrolytes. We explore how radiation damage impacts LLZO’s crystalline structure, Li-ion diffusion, and fracture properties at various temperatures and radiation damage levels. The study aims to provide insights into the competing effects of ion implantation and suggest potential engineering strategies for developing more robust solid-state electrolytes with improved conductivity and dendrite resistance.

Accurate calibration of irradiance measurement devices, or radiometers, is essential for ensuring the reliability of measurements in high heat applications such as concentrating solar power (CSP), aerospace, defense, and pulsed power systems. Despite the critical need, existing calibration standards and service providers are limited to irradiance levels below 100 kW/m² and specific radiation sources, which is insufficient for many applications. For instance, CSP technologies, particularly those under the Department of Energy’s Solar Energy Technologies Office (SETO) Gen 3 program, require accurate measurements of broadband irradiance at levels exceeding 2000 kW/m². In even more extreme scenarios, such as re-entry vehicles, heat levels can surpass 10000 kW/m². Current ISO standards, specifically ISO 14934–2 and ISO 14934–3, are constrained to lower irradiance levels and dependent on black body heat sources, limiting their applicability for high-intensity broadband irradiance measurements, particularly in concentrated solar applications. Here, to address this shortfall, the National Solar Thermal Test Facility (NSTTF) at Sandia National Laboratories (SNL) proposes a calibration method and facility capable of characterizing radiometers up to 2750 kW/m² using concentrated solar irradiance. Calibrating with concentrated sunlight is important for solar applications as it aligns the calibration process with the solar spectrum. This alignment is crucial for minimizing systematic errors and avoiding the need for additional corrections that may arise when radiometers designed for solar applications are calibrated using black-body or electrical sources. This paper presents the present day NSTTF characterization facility and procedure, detailing the proposed calibration method and uncertainty quantification. The presented method builds upon 1980′s NSTTF methodology and involves both theoretical and empirical methods to establish a robust relationship between gauge voltage output and irradiance intensity, quantifying both measurement and fitting errors. By addressing the limitations of existing standards and extending the characterization range, this work provides an advancement in the field of high-intensity irradiance measurement and instrumentation characterization.

We develop a Stockmayer fluid model for molecular dynamics simulations of ionic liquids that captures molecular polarization, ionic conductivity, viscosity, and glass transition temperature, using ethylammonium nitrate (EAN) as an example. The ions in EAN are treated as spheres interacting via the Lennard-Jones potential with an embedded point charge and a permanent dipole moment. We show that our simulation results for EAN are consistent with experimental data and then explore the effects of the molecular parameters on the viscosity of ionic liquids. Our results indicate that viscosity monotonically increases with ionic charge and dipole moment but non-monotonically changes with ionic diameter (or molar volume). This non-monotonic trend arises from the competition among the electrostatic interactions, molecular packing, and size asymmetry between the cation and anion. Our model also shows that long-lived ion pairs result in higher viscosities.

A photocontrolled isomerization reaction of oxabenzonorbornadiene (OBNBD) derivatives with a photo acid generator (PAG) has been developed. Under irradiation, the PAG releases an acid that catalyzes the heterolytic cleavage of the benzyl ether, triggering an aromatization process accompanied by substantial heat release that accelerates the transformation. Large exotherms can further push low reactivity substrates to achieve full conversion within minutes. This approach not only illuminates practical and environmentally benign advantages in organic synthesis but also opens new avenues for developing photoresponsive molecular architectures and advanced functional materials.

Data collected via seismic and infrasound array deployments are leveraged in the geosciences to detect and characterize a myriad of natural and anthropogenic sources. These deployments consist of numerous sensors placed in a predetermined configuration to amplify signal strength and improve the efficacy of array processing techniques used to measure signal directionality and waveform coherence. High‐fidelity feature extraction is often predicated on interstation distance as well as the frequency content and wavelength of an incident signal. Numerous array processing softwares analyze data in sequential frequency bands to obtain a more detailed characterization of a signal. However, current algorithms are limited in their ability to determine optimal array configuration for each band. We introduce an open‐source Python code, called Cardinal, to process seismic and infrasound array data in discretized time–frequency space with the option of applying an adaptive array design to determine optimal subarray configuration for each frequency band. To reduce computational time, the array processing step can be run in parallel using multithreading. Furthermore, the software has the capability to aggregate array processing results from different time–frequency pixels to produce separate sets of detections, or families, with added utility via the application of an adaptive semblance threshold, which aids in isolating signals‐of‐interest from coherent background noise. Upon appropriate configuration, Cardinal exhibits the potential to combine distinct seismic and infrasound phases into separate families.

Causal discovery tools enable scientists to infer meaningful relationships from observational data, spurring advances in fields as diverse as biology, economics, and climate science. Despite these successes, the application of causal discovery to space-time systems remains immensely challenging due to the high-dimensional nature of the data. For example, in climate sciences, modern observational temperature records over the past few decades regularly measure thousands of locations around the globe. To address these challenges, we introduce Causal Space-Time Stencil Learning (CaStLe), a novel meta-algorithm for discovering causal structures in complex space-time systems. CaStLe leverages regularities in local space-time dependencies to learn governing global dynamics. This local perspective eliminates spurious confounding and drastically reduces sample complexity, making space-time causal discovery practical and effective. For causal discovery, CaStLe flexibly accepts any appropriately adapted time series causal discovery algorithm to recover local causal structures. These advances enable causal discovery of geophysical phenomena that were previously unapproachable, including non-periodic, transient phenomena such as volcanic eruption plumes. Regularities in local space-time dependencies are transformed into informative spatial replicates, which actually improve CaStLe's performance when applied to ever-larger spatial grids. We successfully apply CaStLe to discover the atmospheric dynamics governing the climate response to the 1991 Mount Pinatubo volcanic eruption. We provide validation experiments to demonstrate the effectiveness of CaStLe over existing causal-discovery frameworks on a range of geophysics-inspired benchmarks while identifying the method's limitations and domains where its assumptions may not hold.

Pyomo is an open-source optimization modeling software that has undergone significant evolution since its inception in 2008. Pyomo has evolved to enhance flexibility, solver integration, and community engagement. Modern collaborative tools for open-source software have facilitated the development of new Pyomo functionality and improved our development process through automated testing and performance-tracking pipelines. However, Pyomo faces challenges typical of research software, including resource limitations and knowledge retention. The Pyomo team's commitment to better development practices and community engagement reflects a proactive approach to these issues. We describe Pyomo's development journey, highlighting both successes and failures, in the hopes that other open-source research software packages may benefit from our experiences.

Characterization of induced microseismicity at a carbon dioxide (CO₂) storage site is critical for preserving reservoir integrity and mitigating seismic hazards. We apply a multilevel machine learning (ML) approach that combines the nonnegative matrix factorization and hidden Markov model to extract spectral representations of microseismic events and cluster them to identify seismic patterns at the Illinois Basin-Decatur Project. Unlike traditional waveform correlation methods, this approach leverages spectral characteristics of first arrivals to improve event classification and detect previously undetected planes of weakness. By integrating ML-based clustering with focal mechanism analysis, we resolve small-scale fault structures that are below the detection limits of conventional seismic imaging. Our findings reveal temporal bursts of microseismicity associated with brittle failure, providing insights into the spatio-temporal evolution of fault reactivation during CO₂ injection. This approach enhances seismic monitoring capabilities at CO₂ injection sites by improving fault characterization beyond the resolution of standard geophysical surveys.

Frontal ring-opening metathesis polymerization (FROMP) is a promising energy-efficient approach to fabricate polymeric materials. Recent advances have demonstrated FROMP for diverse applications, including additive manufacturing, composites, and foams. However, the characteristic properties of the front are currently controlled primarily by varying the resin composition or the environmental conditions. In this work we present an approach to control FROMP of dicyclopentadiene (DCPD) using photochemical methods. A photobase generator is used to inhibit FROMP of DCPD with UV light while a photosensitizer and co-initiator are used to accelerate FROMP with blue light, enabling orthogonal active photocontrol of front velocity. In addition, photoinhibition-enabled lithographic patterning of frontal polymerizations is demonstrated. Frontal polymerizations are spatially controlled, redirected, and even split into diverging fronts. This work establishes a foundation for advanced control of frontal polymerizations, enabling innovation in traditional and additive manufacturing, as well as emerging processes like morphogenic manufacturing.

Here, a systems governance approach emphasizes a holistic perspective that identifies and navigates the interdependencies and conflicts between security and operational needs. Governance is defined as a collection of metasystems that provide the necessary constraints and processes to support, steer, adapt, transform, and sustain a system (Keating et al. 2022). Utilizing the Cynefin framework, which distinguishes between simple, complicated, complex, and chaotic environments (Snowden and Boone 2007), the article highlights the challenges faced by nuclear power plants in predatory contexts and the importance of integrating security objectives into governance frameworks. By incorporating security as a fundamental aspect of governance, the article underscores its significance for persistence, adaptation, and transformation in the face of uncertainty. Additionally, it introduces key heuristics of systems security, such as the importance of context, knowledge‐based decision‐making, and organization‐specific sociological factors (Williams and Caskey 2024). Ultimately, this work provides valuable insights into enhancing resilient operations in complex environments by reinforcing the connection between effective governance and security in systems engineering.

The widespread implementation of next-generation Li metal anodes is limited, in part, due to the formation of dendritic and/or mossy electrodeposits during cycling. These morphologies can lead to battery failure due to the formation of short circuits and significant volumetric expansion at the anode. One strategy to control the electrodeposition of Li metal is to use lithiophilic materials at the anode. Here, we evaluate the impact of Ag and Au on the early stages of Li metal electrodeposition and cycling. The alloying substrates decrease the voltage for Li reduction and improve Li wetting/adhesion. We probe volumetric expansion directly through dilatometry measurements and find that the degree of volumetric expansion is less when lithium is cycled on an alloying substrate compared to a non-alloying substrate (Cu). Dilatometry experiments reveal that Au has the least amount of volumetric expansion and coin cell cycling experiments indicate that Ag yields more stable cycling compared to Au or Cu. The evaluation of in situ cross-sectional images of cycled coin cells shows that Ag has the lowest volumetric expansion in a coin cell format.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Cements are a critical component in well construction, as they act to prevent well fluid and gas escape, prevent corrosion of the casing, and strengthen the wellbore to prevent deformation. Under the high temperature/pressure conditions common in geothermal systems, the injection of cold water for energy production is expected to induce cyclic damage to the borehole cement through the rapid temperature fluctuations. These “thermal shocks” are expected to cause casing shrinkage, annulus formation, and cement tensile stresses. To understand the effect of cold water injection on the wellbore environment, a set of rock-cement-steel samples were created to simulate the structure of a geothermal well. Lightweight thermally-insulating cement blends were tested under thermal shock conditions in this study. In each test, the samples were pressurized to an effective pressure of ∼3.5 MPa and placed at high temperatures. Thermal shocks were performed by injecting cold water (∼10-15 °C) through the samples at a constant rate while keeping the samples at high temperatures until the sample temperature stopped decreasing and deformation ceased. Eight thermal shock tests were conducted with each sample – two at 100 °C and six at 200 °C. Post-tests analysis was then conducted by cutting open each sample to examine the damage in each component of the simulated wellbore. Experimental results suggest that all samples experienced similar degrees of axial and lateral contraction during cold water injection, but for the most part this contraction is recoverable when injection halts. Post-test analysis revealed that fly ash cenosphere pre-treatment had the best effect on improving thermal shock resistance in the cement blends. Thermomechanical modeling of likely stress paths experienced by the cements during heating/cooling cycles shows that elasto-plastic cement constitutive behavior results in most plastic strain occurring during the initial heating steps, with mostly elastic strain occurring during the thermal shock cycles. This agrees with experimental evidence, suggesting that cement damage from shocking occurs via other mechanisms such as chemical alteration, corrosion, and fatigue.

This manuscript we present a new method for solving PDEs. The method utilizes a mixed formulation of Tensor Train (TT) and Quantized Tensor Train (QTT), designed for the spectral element discretization (Q1-SEM) of the time-dependent convection-diffusion-reaction (CDR) equation. We reformulate the assembly process of the spectral element discretized CDR to enhance its compatibility with tensor operations and introduce a low-rank tensor structure for the spectral element operators.

Abstract not provided.

To date, careful data treatment workflows and statistical detectors are used to perform hyperspectral image (HSI) detection of any gas contained in a spectral library, which is often expanded with physics models to incorporate different spectral characteristics. In general, surrounding evidence or known gas-release parameters are used to provide confidence in or confirm detection capability, respectively. This makes quantifying detection performance difficult as it is nearly impossible to develop an absolute ground truth for gas target pixel presence in collected HSI. Consequently, developing and comparing new detection methods, especially machine learning (ML)-based methods, is susceptible to subjectivity in derived detection map quality. In this work, we demonstrate the first use of transformer-based paired neural networks (PNNs) for one-shot gas target detection for multiple gases while providing quantitative classification and detection metrics for their use on labeled data. Terabytes of training data are generated from a database of long-wave infrared HSI obtained from historical Mako sensor campaigns over Los Angeles. By incorporating labels, singular signature representations, and a model development pipeline, we can tune and select PNNs to detect multiple gas targets that are not seen in training on a quantitative basis. We additionally assess our test set detections using interpretability techniques widely employed with ML-based predictors, but less common with detection methods relying on learned latent spaces.

This manual provides step-by-step instructions for installing, configuring, and using BLADE, an automated framework for analyzing and classifying bolide light curves from NASA CNEOS datasets. BLADE enables efficient, reproducible analysis of atmospheric entry phenomena, supporting planetary defense and atmospheric science research through standardized signal processing and event classification.

This report addresses the need to predict infrasound signal arrival times and back azimuths at monitoring stations, enabling more focused and efficient searches within recorded waveform data.

Abstract not provided.

The authors were tasked with writing a concise memo (now a short report) that adequately addresses the current ``state-of-the-art'' in the field of reactive flow modeling and burn models for mid-year 2025, along with identifying some of the modeling gaps. It is assumed that the reader has experience with burn models and running hydrocode simulations, and thus every effort is made to ensure brevity, clarity, and utility, so that this document will serve as a quick and useful reference. Note that

List of accomplishments from FY2024. All content is UUR.

Collection of news articles on international nuclear developments for the global security community.

The structural frame used for conducting large structural tests was replaced by the ACME (Applied Combined Mechanical Environments) Lab in Building 860 of Sandia National Laboratories in Albuquerque, NM. This new design required simulation and analysis to verify maximum loads and corresponding factors of safety. Two major assemblies which are part of the ACME Lab, the Modular Wall and A-Frame, were designed to operate with a multitude of components in various configurations and with varying lo

Abstract not provided.

Abstract not provided.

ACS Photonics

Sandia National Laboratories (SNL) collaborated with the US Nuclear Regulatory Commission's (NRC) Office of Nuclear Security and Incident Response (NSIR) and Office of Nuclear Regulatory Research (RES) to investigate the reasonable variability in the estimation of dose exceedance distances. SNL performed calculations using the MACCS code to elucidate the sensitivity of dose exceedance distances to variations in user input parameters.

Abstract not provided.

Poster for Intern Symposium

Abstract not provided.

Understanding the safety profile of aged Li-ion batteries is essential for developing effective battery management and hazard mitigation strategies. However, most safety assessments have focused on fresh batteries, with just a few calorimetry studies on aged batteries with metal oxide positive electrodes. This study provides a broad assessment of commercial 18650-type Li-ion batteries with NCA, NMC, and LFP positive electrodes, both uncycled and aged under conditions that promoted solid electrolyte interphase (SEI) growth as the dominant degradation mechanism. The cells underwent mechanical (nail penetration, crush), electrical (overcharge, overdischarge), and thermal (accelerating rate calorimetry) abuse tests. Safety was rated on general characteristics such as mass loss, maximum temperature, and EUCAR (European Council for Automotive R&D) hazard level, as well as characteristics specific to individual abuse tests. Generally, aged cells with SEI growth exhibited similar or improved safety compared to uncycled cells, contrasting with our previous findings on NCA cells with Li plating as the dominant aging mechanism (Part I of this series). Yet, some tests and characteristics indicated reduced aged cell safety, such as earlier triggering of mechanical failure. These results emphasize the need to examine aged battery safety across diverse empirical techniques, degradation modes, and chemistries.

Identification of amorphous phases in nanoparticles by atomic-resolution transmission electron microscopy (TEM) requires analyses such as tilt-angle-dependent TEM imaging and single-nanoparticle electron diffraction, rather than relying on a single TEM image.