Publications

Abstract not provided.

This report details the heat source impact testing performed at the Shock Thermodynamics Applied Research Facility (STAR). The purpose of these tests were to measure the impact behavior of a heated Ta-10W tantalum/tungsten alloy heat source exemplar. Each exemplar resembled a cylinder with a rectangular through-hole orthogonal to the cylinder axis-of-symmetry. In each test, the exemplar was impacted by a hardened steel impactor accelerated using the STAR air gun. The exemplars were impacted ...

A comprehensive investigation into the pooling and vaporization of liquid hydrogen spills onto concrete and steel surfaces in a steady cross-wind is presented in this work. This is the first instance of liquid hydrogen pooling and dispersion in a steady environment. A high-capacity fan in a large tunnel was used to generate the steady cross-winds while spilling roughly 10-20, or 40 l/min of liquid hydrogen onto substrates. Temperature and extractive concentration measurements were made, in ad

Abstract not provided.

Abstract not provided.

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.

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.

Poster for Intern Symposium

Abstract not provided.

Abstract not provided.

Abstract not provided.

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

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

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.

Phenomenological CALPHAD (CALculation of PHAse Diagrams) models, widely used for multicomponent materials, often contain a considerable number of parameters and require fitting using data from a relatively small number of experimental measurements or theoretical calculations. Sometimes these parameters are introduced for the purpose of improving model fits but without clear physical justification, which leads to overparametrized models with poor generalization performance. Automated approaches for optimal model selection based on the available data therefore become critical. In this work, a least absolute shrinkage and selection operator (LASSO)-based approach is developed for model selection by leveraging the linearity of the CALPHAD model with respect to its parameters to convert the model selection and fitting to a LASSO minimization problem. We demonstrate its utility for thermodynamic modeling of thermochemical hydrogen (TCH) production materials using lanthanum strontium manganite (LSM) as an example. Various TCH-relevant properties, including oxygen stoichiometry as a function of oxygen partial pressure, enthalpy of reduction, and entropy of reduction, are successfully predicted with reasonable accuracy using a minimal set of model parameters. Importantly, the model selection and fitting involve minimal human decision; it can therefore be applied to high-throughput DFT defect calculations and yield efficient workflows for TCH material modeling and optimization.

We explore a nonvariational quantum state preparation approach combined with the ADAPT operator selection strategy in the application of preparing the ground state of a desired target Hamiltonian. In this algorithm, energy gradient measurements determine both the operators and the gate parameters in the quantum circuit construction. We compare this nonvariational algorithm with ADAPT-VQE and with feedback-based quantum algorithms in terms of the rate of energy reduction, the circuit depth, and the measurement cost in molecular simulation. We find that, despite using deeper circuits, this new algorithm reaches chemical accuracy at a similar measurement cost to ADAPT-VQE. Since it does not rely on a classical optimization subroutine, it may provide robustness against circuit parameter errors due to imperfect control or gate synthesis.

Recent studies have applied variational calculus, conformal mapping, and point transformations to generalize the one-dimensional (1D) space-charge limited current density (SCLCD) and electron emission mechanisms to nonplanar geometries; however, these assessments have focused on extending the Child–Langmuir law (CLL) for SCLCD in vacuum. Since the charge in the diode is independent of the coordinate system (i.e., covariant), we apply bijective point transformations to extend the Mott–Gurney law (MGL) for the SCLCD in a collisional or semiconductor gap to nonplanar 1D geometries. This yields a modified MGL that replaces the Cartesian gap distance with a canonical gap distance that may be written generally in terms of geometric scale factors that are known for multiple geometries. We tabulate results for common geometries. Such an approach may be applied to any current density, including non-space-charge limited gaps and SCLCD that may fall between the CLL and MGL.

One-dimensional Coulomb drag has been an essential tool to probe the physics of interacting Tomonaga-Luttinger liquids. To date, most experimental work has focused on the linear regime while the predictions for Luttinger liquids beyond the linear response theory remain largely untested. In this Letter, we report measurements of reciprocal momentum transfer induced Coulomb drag between vertically coupled quasi-one-dimensional quantum wires in the nonlinear regime. Measurements were performed at ultralow temperatures between wires only 15 nm apart. Our results reveal a nonlinear dependence of the drag voltage as a function of the drive current superimposed with an oscillatory contribution, in agreement with theoretical predictions for Coulomb drag between Tomonaga-Luttinger liquids. Additionally, the observed current-voltage characteristics exhibit a nonmonotonic temperature dependence, further corroborating the presence of non-Fermi-liquid behavior in our system. In conclusion, these findings are observed both in the single and in the multiple subband regimes and in the presence of disorder, extending the onset of this behavior beyond the clean single channel Tomonaga-Luttinger regime where the predictions were originally formulated.

This study investigates the fatigue crack growth rate (FCGR) behavior of pipeline and low-alloy pressure vessel steels in high-pressure gaseous hydrogen. Despite a broad range of yield strengths and microstructures ranging from ferrite/pearlite, acicular ferrite, bainite, and martensite, the FCGR in gaseous hydrogen remained consistent (falling within a factor of 2–3). Steels with higher fractions of pearlite, typical of older vintage pipeline steels, exhibited modestly lower crack growth rates in gaseous hydrogen compared to steels with lower fractions of pearlite. Crack growth rates in these materials exhibit a systematic dependence on stress ratio and partial pressure of hydrogen, as captured in the recently published fatigue design curves in ASME B31 code case 220 for pipeline steels and ASME BPVC code case 2938 for pressure-vessel steels.

Mesoscale modeling of shock waves in Ni+Al multilayers poses significant challenges that are due, in part, to shock-induced chemical reactions. Current modeling approaches utilize reactive molecular dynamics (MD), but they are limited to resolving domains of only a few hundred nanometers. In contrast, actual multilayer superlattices can be tens of micrometers thick, and they exhibit non-ideal (i.e., wavy) interfaces. The second part of our research builds upon previous work developing physically based, thermodynamically complete equations of state for various Ni and Al intermetallic compositions. Here, we introduce a novel workflow for high-fidelity mesoscale simulations of Ni+Al multilayers using a continuum hydrocode. By increasing the simulation domain size beyond MD limitations (e.g., 2 × 6 μm²) and incorporating explicit interfacial roughness, we investigate the shock response of Ni+Al multilayers at previously unexplored scales. Our experimental design encompasses nine multilayer geometries with varying roughness amplitudes and tilt angles (θ = 15°, 30°, and 45°), alongside 19 flyer impact velocities ranging from 0.3 to 3.0 km/s, resulting in a total of 171 high-fidelity simulations. The bulk shock state from inert 2D mesoscale simulations aligns with the law of mixtures, while temperature and pressure fluctuations strongly correlate with multilayer geometry types. A new metric dubbed the “hot spot probability integral” shows a greater dependence on a tilt angle than interfacial roughness.

Geogenic Helium-4 (⁴He) in-situ increases locally in regions of large deformation generated naturally or anthropogenically. This gas release by deformation is a potential geochemical precursor signal for subsurface deformation. To evaluate the applicability of ⁴He degassing for correlating deformation in different lithologies, we conducted high force crush tests, up to 97,800 N axial load, to assess the total ⁴He released during fragmentation of the rocks. We observed that the highest ⁴He released occurred in the sedimentary rocks and that release correlated strongly with lithologic age and U/Th content. Microstructural changes of the pre- and post-test rocks indicate that the degree of grain size reduction relates directly to the total ⁴He released during crushing. The range of in-place ⁴He was calculated based XRF measurements of uranium and thorium in each lithology, with the results indicating that the majority of the trapped ⁴He was not released. However, the ⁴He released by deformation depended upon how the each rock deformed during deformation and the degree of grain size reduction. We postulate that ⁴He precursor signals can be used to understand subsurface deformation only if geomechanical and geochemical conditions for ⁴He enrichment in a lithology are met.

Time Reversal (TR) is a signal processing technique that can be used to focus acoustic waves to a specific location in space, with most applications aiming to create an impulsive focus. This study instead aims to focus long-duration noise signals using TR. This paper seeks to generate higher amplitude noise at a desired location over an existing method of broadcasting equalized noise. Additionally, this paper explores various characteristics associated with focusing long duration noise using TR. The dependence of the focal amplitude on the duration of the focused signal is explored as well as the implications of using multiple sources when focusing noise. The focal amplitude decreases with longer duration and then levels off when the duration exceeds a few seconds. Coherent addition of focused noise is observed if all loudspeakers have coherent noise signals convolved with their reversed impulse responses. Lastly, focusing noise with a desired spectrum is explored.

This paper proposes a new approach for the calibration of material parameters in local elastoplastic constitutive models. The calibration is posed as a constrained optimization problem, where the constitutive model evolution equations for a single material point serve as constraints. The objective function quantifies the mismatch between the stress predicted by the model and corresponding experimental measurements. To improve calibration efficiency, a novel direct-adjoint approach is presented to compute the Hessian of the objective function, which enables the use of second-order optimization algorithms. Automatic differentiation is used for gradient and Hessian computations. Two numerical examples are employed to validate the Hessian matrices and to demonstrate that the Newton–Raphson algorithm consistently outperforms gradient-based algorithms such as L-BFGS-B.

We introduce a grain boundary (GB) translation vector, tWS, to describe and quantify the domain of the microscopic degrees of freedom of GBs. It has long been recognized that for fixed macroscopic degrees of freedom of a GB there exists a large multiplicity of states characterized by different relative grain translations. More recently another degree of freedom, [n], the number of GB atoms, has emerged and is now recognized as an equally important component of GB structural multiplicity. In this work, we show that all GB microstates can be uniquely characterized by their value of tWS, which is located within the Wigner–Seitz (WS) cell of the Displacement-Shift-Complete lattice (DSCL) of the GB. The GB translation vector captures information about both the translation state and the number of GB atoms. We show that the density of GB microstates inside the WS cell of the DSCL is not uniform and can form clusters that correspond to different GB phases. The vectors connecting the centers of the clusters correspond to the Burgers vectors of GB phase junctions, which can be predicted without building the junctions. Using tWS, we quantify GB excess shear and argue that it is defined up to a DSCL vector, which has implications for thermodynamic equilibrium conditions. Additionally, this work generalizes the definition of the number of GB atoms [n] to asymmetric boundaries.

Abstract not provided.

The creep behavior of 316 L stainless steel at room temperature was evaluated as a function of time and applied stress using a new high-Throughput approach. Several common creep models were evaluated against the observations, leading to deeper analysis of a stress-dependent modified logarithmic creep model. Within this model, multiple sources of uncertainty were compared. Aleatoric stochastic variation between samples under nominally identical conditions was identified as the primary contributor to uncertainty in creep response. Under any particular set of conditions, the sample-To-sample variability in creep strain was as high as a factor of two, highlighting the engineering importance of characterizing large statistical datasets. The model's extrapolation capabilities were assessed by comparing predictions derived from calibration on partial, shorter-duration subsets of the data. These findings underscore the importance of accounting for stochastic effects in predictive modeling of aging phenomena.

A variational approach is developed with a meshless discretization to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. We propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework, providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5× improvement in accuracy and allows for taking an 8× larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

Many overarching standards, regulations, or other requirements necessitate that the critical temperature (Tcrit) of an energetic material be known or estimate prior to operations such as heating. However, they rarely, if ever, provide a specific method for doing so. While other methods exist to calculate Tcrit, such as the Frank-Kamenetskii equation, computer simulation, etc., these cannot be readily utilized for energetics that are physical mixtures (i.e. most pyrotechnics) or materials where detailed material properties required for such calculations are lacking. In this study, a COTS (commercial off-the-shelf) SBAT (Simulated Bulk Auto-ignition Test) apparatus is modified to perform Henkin cookoff tests. This creates a simple, efficient, and cost-effective solution to estimate required Tcrit values for energetic materials. After a historical overview of the evolution of the Henkin test, several common energetic materials (PETN, RDX, etc.) for which Tcrit has been readily calculated, and historical data is available, were analyzed for comparison and verification purposes. This was followed by a variety of pyrotechnic mixtures and other materials to where calculation methods cannot be readily used. The modified apparatus, as well as an updated sample shell sealing method, produced results that aligned well with historical data. It also produced reasonable Tcrit estimates for those materials where the Frank-Kamenetskii equation cannot be applied.

Bentonite-based engineered barrier system (EBS) is a key component of many repository designs for the geological disposal of high-level radioactive waste. Given the complexity and interaction of the phenomena affecting the barrier system, coupled thermo-hydro-mechanical (THM) numerical analyses are a potentially useful tool for a better understanding of their behaviour. In this context, a Task (the Horonobe EBS experiment) was undertaken to study, using numerical analyses, the thermo-hydro-mechanical (and thermo-hydro) interactions in bentonite based engineered barriers within the international cooperative project DECOVALEX 2023. One full-scale in-situ experiment and four laboratory experiments, largely complementary, were selected for modelling. The Horonobe EBS experiment is a temperature-controlled non-isothermal experiment combined with artificial groundwater injection. The Horonobe EBS experiment consists of the heating and cooling phases. Six research teams performed the THM or TH (depended on research team approach) numerical analyses using a variety of computer codes, formulations and constitutive laws. For each experiment, the basic features of the analyses are described and the comparison between calculations and laboratory experiments and field observations are presented and discussed.

Defect reduction remains a critical objective in the integrated circuit manufacturing process, particularly within the highly re-entrant lithography modules where minimizing defects is crucial. Defects at the wafer edge can contaminate lithography modules and downstream processing equipment, leading to redistribution onto the wafer surface and adversely affecting overall device yield. A persistent challenge in the resist coating process is the formation of resist edge beads, driven by the strong Van der Waals attraction of excess photoresist (PR) to itself and the underlying substrate. The edge bead removal (EBR) process is a standard cleaning step designed to eliminate these edge beads and prevent potential contamination. In this study, we identify the sources of EBR induced defects and additional EBR process encroachment toward edge patterning during the EBR cleaning process. This study provides a comprehensive study aimed at optimizing the EBR cleaning process to effectively eliminate EBR-induced defects, thereby enhancing overall device yield. Specifically, we identify three primary defects induced by the EBR cleaning process: rainbow-type, finger-shaped, and teardrop-type defects. Our experimental study reveals that in addition to EBR rinse time, PR cast time is crucial parameters contributing to the formation of these defects. By properly optimizing the PR cast time and EBR rinse time, we were able to remove nearly 100 % of dense clusters of defects that were easily visible even at low magnification optical microscopy throughout the wafer edge. We observed that shorter PR casting times shows edge defects caused by inefficient EBR process because of insufficient time for PR to fully settle causing superfluous PR to continue flowing toward wafer edge during EBR clearing step, leading to partial removal of PR at the wafer edge and the formation of rainbow defects. Proper optimization of both PR casting time and EBR chemistries dispense time is essential to resolve these defects, ensuring efficient EBR cleaning process and improved overall device yield.

This study investigates the application of microbeads as an innovative encapsulation technique to protect electronic components from harsh mechanical strain. Traditional encapsulation methods using hard epoxy provide substantial mechanical support but create thermal expansion mismatch issues, potentially leading to electronic component failure. We explore the use of finely powdered microbeads to achieve protective structures combining stiffness and energy absorption. The research focuses on key variables, including microbead size, microbead roughness, compaction of microbeads, and circuit board mounting in the encapsulation, all of which influence the encapsulation’s effectiveness. Experimental setups and testing protocols were developed to assess the performance of various microbead materials under different impact conditions. Results demonstrate that microbead encapsulation significantly reduces strain on circuit boards, minimizing the risk of damage during mechanical shocks. However, challenges remain, such as optimizing microbead characteristics and modeling their behavior within large-scale circuit board assemblies. Despite these challenges, the findings suggest that microbead encapsulation offers a promising alternative to conventional methods, enhancing the durability and reliability of electronic components in high-stress environments.

The Ghareb formation, a shallowly buried porous chalk in Israel, is currently a candidate for nuclear waste disposal. The potential repository is somewhat unique for its host rock and emplacement in shallow (500 m) engineered large diameter boreholes. Herein, the thermal properties of the Ghareb are determined to support design and performance assessment; the relevant properties measured are thermal conductivity, specific heat capacity, thermal diffusivity, and thermal expansion coefficient along with their relationship with varying temperature. For the temperature range of 40 to 275 °C, the thermal conductivity ranges from 0.30 to 1.10 W/m·K, the thermal diffusivity ranges from 0.20 to 0.72 mm2/s, and volumetric heat capacity ranges from 0.86 to 2.00 MJ/m3·K. Thermal strain measurements were used to estimate the linear thermal expansion coefficient to be 6·10–4–9·10–2 °C−1 from 40 to 300 °C. These measured properties were used in a thermomechanical model to estimate near-field stresses an hour and 10 years after waste emplacement; the borehole was found to be stable. Thermal loading after 10 years was predicted to elevate local pore pressures by 1–1.5 MPa. The laboratory measurements coupled with analyses are the first attempts at performance assessment characterization for this first of its kind potential repository setting with this chosen host rock.

The pharmaceutical drug product development process can be greatly accelerated through the use of modeling and simulation techniques to predict the manufacturability and performance of a given formulation. The anticipation and possible mitigation of tablet damage due to manufacturing stresses represents a specific area of interest in the pharmaceutical industry for predicting formulation and tableting performance. While the finite element method (FEM) has been extensively used for predicting the mechanical behavior of powder material in the compaction processes, a shortcoming of the approach is the inherent difficulty to predict discontinuities (e.g., damage or cracking) within a tablet as FEM is a continuum-based approach. In this work, we propose a novel method utilizing peridynamics (PD), a numerical method that can capture discontinuities such as tablet fracture, to predict the evolution of damage and breakage in pharmaceutical tablets. The approach links (1) the finite element method – to elucidate the behavior of powders during die compaction – with (2) the peridynamics modeling technique – to model the discontinuous nature of damage and predict tablet breakage during the critical stages of unloading and ejection from the compression die. This short communication presents a proof of concept including a workflow to calibrate the linked FEM-PD simulation models. It demonstrates promising results from a preliminary experimental validation of the approach. Following further development, this approach could be used to guide the optimization of compression processes through targeted changes to formulation material properties, compression process conditions, and/or tooling geometries to deliver improved process efficiency and tablet robustness.

Electrical resistivity tests can be used to evaluate the transport properties of concrete and provide a durability assessment. However, the electrical resistivity is largely dependent on the pore solution composition and recent work suggests that some aggregates have the capacity for cation uptake. This study first aims to provide further evidence for adsorption of cations on aggregate surfaces, without formation of reaction products (e.g. alkali-silica reaction). Secondly, hardened mortar samples were prepared using a fine aggregate with a high alkali affinity and a non-reactive fine aggregate as a control. The electrical resistivity of mortars was measured, and the pore solution of these mortars was obtained through high-pressure extraction. The effect of aggregate moisture dilution on the pore solution was decoupled by using a pore partitioning model. The results indicate that aggregate minerology can influence the pore solution composition through cation uptake. Specific minerals of minor quantity, like biotite, may be responsible for cation exchange. While adsorbed cations strongly affected pore solution and formation factor measurement, the bulk resistivity measurements on hardened mortar were only marginally influenced. Research on other implications of similar aggregate interactions with pore solutions are an intriguing area for future research.

This research proposes a sensitivity-based framework for selecting the optimal prescribed loading path for a biaxial cruciform specimen. Optimality here is determined by the direction and magnitude of the prescribed displacement that minimizes the influence of random noise on the material model parameter identification. Using simulated experimental data based on finite element simulation, in this work, we identify the material model parameters of a Ludwik hardening model and plane stress implementation of the Hill-48 yield criterion using finite element model updating (FEMU). Our analysis reveals that the identification (or estimator) uncertainty of model parameters depends on the displacement boundary conditions (i.e., loading sequence) and the ground-truth value of the individual parameters. Optimal experimental design (OED) criteria based on the Fisher information matrix were investigated to mitigate indecision in the choice of optimal load path when the identification uncertainty of different material model parameters optimized at different load paths. The determinant of the Fisher information matrix was chosen here as the more useful metric due to its ability to capture uncertainty of the most influential material model parameters. The proposed framework demonstrates potential for real-time automated load step selection using scalar criteria derived prior to mechanical loading. The framework can be generalized to other geometries, boundary conditions and material models, allowing this procedure to be utilized for different experimental configurations and materials.

The explosive growth in data collection and the need to process it efficiently, as well as the desire to automate increasingly complex tasks in transportation, medical care, manufacturing, security and many other fields have motivated a growing interest in neuromorphic computing. Unlike the binary, transistorbased ON/OFF logic gates and separate logic and memory functionalities employed in digital computing, neuromorphic computing is inspired by animal brains that use interconnected synapses and neurons to perform processing, storage and transmission of information at the same location, while only consuming ~20 W or less of power. Motivated by the brain’s efficiency, adaptability, self-learning and resiliency qualities, neuromorphic computing can be broadly defined as an approach to processing and storing information using hardware and algorithms inspired by models of biological neural systems. Present research in neuromorphic computing encompasses approaches that vary significantly in their degree of neuro-inspiration, from systems that only incorporate features such as asynchronous, event-driven operation or use crossbar arrays of non-volatile memory (NVM) elements to accelerate deep neural networks (DNNs), to designs that embrace the extreme parallelism, sparsity, reconfigurability, adaptability, complexity and stochasticity observed in nervous systems. The term ‘neuromorphic’ computing is often credited to Carver Mead, who in the 1980s investigated Si-based analog electronics to replicate functions of the animal retina. Earlier important advances in this field include the work of Frank Rosenblatt, who proposed the concept of the perceptron, Bernard Widrow, who used this concept to build one of the first analog neural networks, the Adaline and many other researchers (see ref. 6 for an historical perspective on neuromorphic computing). With the recent increase in the use of artificial intelligence and large language models, and rising concerns over the associated energy costs, interest in neuromorphic hardware has expanded rapidly. According to some estimates, driven largely by the drastic growth in the training use of artificial intelligence (AI) models using the current computing architectures, the energy cost of computing is projected to reach the energy supply worldwide by 2045. Furthermore, while this is not a realistic outcome, it means that, if more efficient computing technologies are not developed -- soon -- the world will soon become one where demand for energy and market constraints limit the continued increase of societal access to AI and cloud services from data centers. Data centers used for training and use of these models consume hundreds of terawatt hours of electricity, already past 4% of the US electricity demand.

As the power grid is undergoing rapid transformations, numerous questions are emerging about its vulnerability to wide-area extreme events (WAEE), which could influence its operations. Relatively few analyses have been conducted to date regarding the impact of correlated outages during WAEEs on the grid’s ability to balance resources with demand. This study addresses this gap by conducting a resource adequacy analysis for a hurricane-inspired WAEE on a 2035 synthetic power grid system for the United States. A sensitivity analysis was also conducted to characterize the relative impact of weather on unserved energy. Our results indicate that although the magnitude and duration of the shortfalls vary depending on weather conditions, persistent shortfalls are observed in some regions. Initial explorations indicate a strong correlation between transmission-constrained regions and regions with persistent shortfalls. Future work could generate empirically-grounded representations for generator outages as well as conduct causal analyses of these shortfalls to improve understanding of drivers as well as possible mitigation strategies. Continued exploration of extreme weather impacts on the grid is important to develop more robust understanding of the reliability and resilience of our power systems, especially as they undergo rapid transformations.

The daily conduct of high-risk, high-consequence work, by its very nature, can be mentally demanding. Research demonstrates that failure to manage or mitigate these demands can degrade psychological well-being and mental health. Degradations in employee well-being impact individual performance, jeopardizing safety and mission success in the workplace. This Commentary describes efforts taken at Sandia National Laboratories over the past eight years to evaluate, understand, learn from, and mitigate factors such as stress, burnout, work-life imbalances, and disengagement in the workplace that can degrade employee well-being. After evaluating these four facets of employee well-being, Sandia National Laboratories initiated several programs to address observations and outcomes, including the Thrive program, the Take 10 Initiative, work-life balance resources, employee resource groups such as the Sandia Parents Group, and Workplace Improvement Networks. Collectively, the intent of these programs is to ensure employee mental readiness to conduct hazardous high-risk work effectively and safely. Preliminary data suggest that these programs are succeeding. Other national laboratories and organizations, regardless of size, may wish to apply similar approaches to improve employee well-being and thereby increase the likelihood of mission success.

Photonic integrated circuits incorporating intersubband transitions are ideal for mid-infrared and terahertz nanophotonics. However, the design of epitaxies has long been inhibited by two factors: the modest predictivity of ab initio theory and the absence of absolute intersubband gain and loss measurements under operating conditions. Existing measurements either yield inaccurate gain profiles or do not accurately assess dependence on frequency, bias, and temperature. Here, we present a delay-resolved absolute-referencing method for accurate gain evaluation without these limitations, addressing a long-standing challenge. By creating a photonic circuit that allows broadband pulses to traverse different lengths of a gain medium, we measure the absolute transmission of intersubband structures. Gain profiles match theoretical predictions at lower temperatures, with gain and dispersion clamping after lasing, and faster-than-expected degradation occurs at higher temperatures. Our approach provides a precise experimental evaluation of temperature-dependent gain performance and gives insight into optimizing temperature performance and frequency comb designs.

This study challenges the assumption of the non-flammability of lithium metal all-solid-state batteries (LiSSBs) and other lithium metal batteries without flammable electrolytes. Through thermodynamic calculations and ex situ experiments, we reveal for the first time the risk of thermite reactions between lithium metal and LiFePO4 in both charged and discharged states. Reactivity is worsened by excess lithium metal in the cell, reaching final maximum adiabatic temperatures of 2,500°C in the charged state, which is hot enough to boil lithium. The thermite reaction triggers spontaneously at 500°C, with poor surface contact, while increasing surface contact through mixing initiates the reaction at room temperature in an inert environment. Despite its fast kinetics, this reaction is transport limited due to lithium passivation, leading to long burn times and reignition risks. Given the risk of lithium metal contacting the cathode during failure, understanding these reactions is crucial for ensuring the safe deployment of LiSSBs.

Silicon processing techniques such as atomic precision advanced manufacturing (APAM) and epitaxial growth require surface preparations that activate oxide desorption (typically >1000 °C) and promote surface reconstruction toward atomically clean, flat, and ordered Si(100)-2 × 1. We compare the aqueous and vapor phase cleaning of Si and Si/SiGe surfaces to prepare APAM-ready and epitaxy-ready surfaces at lower temperatures. Angle resolved X-ray photoelectron spectroscopy (ARXPS) and Fourier transform infrared (FTIR) spectroscopy indicate that vapor hydrogen fluoride (VHF) cleans dramatically reduce carbon surface contamination and allow the chemically prepared surface to reconstruct at lower temperatures, 600 °C for Si and 580 °C for a Si/Si0.7Ge0.3 heterostructure, into an ordered atomic terrace structure indicated by scanning tunneling microscopy (STM). After thermal treatment and vacuum hydrogen termination, we demonstrate STM hydrogen desorption lithography (HDL) on VHF-treated Si samples, creating reactive zones that enable area-selective chemistry by using a thermal budget similar to CMOS process flows. We anticipate that these results will establish new pathways to integrate APAM with Si foundry processing.

We quantify the effect of a nuclear-reactor environment on the hydrogen isotope equilibrium vapor pressure over pure zirconium and zirconium hydride. A vacuum-sealed capsule containing a zirconium foil with 6 atom% deuterium was irradiated at a neutron flux of ~10¹⁴ cm^-2 s^-1 at the University of Missouri Research Reactor (MURR). The internal stainless-steel (SS) sample holder acted as the heat source via gamma absorption. To measure low desorption pressures in a high-flux environment, we developed a method to transduce pressure from the measured sample temperature during irradiation, calibrating with known deuterium pressures in unirradiated capsules at various heating powers using an internal filament-heated system designed to mimic irradiation-induced heating. Our temperature-pressure transduction method operates similarly to a Pirani or thermocouple pressure gauge. The in-reactor measurements revealed a roughly 4-fold enhancement in desorption pressure after only 6 h of irradiation (~2 × 10¹⁸ cm^-2 neutron fluence) compared to thermal desorption in control experiments, indicating a nonthermal contribution from neutron irradiation. The slower temperature/pressure stabilization rate in the reactor suggests that desorption pressure enhancement increases with neutron fluence. Further, this enhancement signifies increased solubility of hydrogen isotopes in zirconium during irradiation. We propose that high-energy neutron collisions with hydrogen isotopes in hydrides lead to their decomposition at lower temperatures, supersaturating the surrounding αZr lattice and resulting in higher desorption pressure, which continues to rise as more hydrides dissolve with increasing neutron fluence.

Ongoing research in isolating high-level nuclear waste and spent fuel has highlighted compacted bentonite as a suitable material for engineered barrier systems in deep geological repositories due to its extraordinary swelling and retention properties. This research focuses on the chemo-mechanical behavior of compacted bentonite exposed to different pore fluids with different concentrations and loading conditions. The study involves swelling pressure and compressibility experiments along with mineralogy analysis employing X-ray diffraction (XRD) and Cation exchange. The tests were conducted on BCV (a Mg/Ca- bentonite) compacted at a dry density of 1.48 ± .02 Mg/m³. An advanced chemical-mechanical constitutive model for unsaturated highly expansive clays was adopted to simulate the material response and better understand its behavior. The model is able to account for the main phenomena at both macro and microstructural levels and the interactions between them. The model successfully replicated experimental observations. The XRD analyses support the macroscopic observation, indicating that salinity impacts crystalline swelling as demonstrated by the reduction of basal spacing from 19.27 Å to 15.68 Å when the osmotic suction increases from 0 MPa to 33 MPa. The results suggested that the osmotic pressure generated by the concentration in the pore fluids promotes a reduction in swelling pressures, swelling strains, and crystalline swelling of clay minerals. Also, it affects the pre-consolidation stress and the compressibility of the compacted samples. In conclusion, it was also observed that both solution type and solution concentration impact the clay swelling pressure.

Here, a novel approach is presented for computing general rigid body motion based on a few known linear accelerations. This method utilizes linear acceleration data obtained from three distinct points on the body, all within a body-fixed reference frame. The only requirement is that the three chosen points must not be collinear. A system of differential-algebraic equations is derived, combining principles of rigid body kinematics with theory of the rotation group SO(3). These equations provide a framework for numerically computing various motion parameters, including angular velocity, angular acceleration, body orientation, velocity field, acceleration field, and displacement field. By numerically solving this system of equations, we can fully characterize rigid body motion in three-dimensional space. A numerical example is provided to demonstrate the practical implementation and efficacy of the proposed technique, illustrating its potential for accurate motion computation in various applications.

Sandia National Laboratories is currently investigating scalable architectural simulation capabilities, with a focus on simulating and evaluating highly scalable supercomputers for high-performance computing applications. This exploration is driven by the shift toward more specialized forms of compute and the need for a more diverse set of accurate models. This project will explore the use of General-Purpose Graphical Processing Units (GPGPUs) in high-performance computing using both physical systems and new simulator models – traditional GPUs as well as tightly-coupled SIMT accelerators.

Advancement of Earth System Models (ESMs) is becoming increasingly challenging due to a confluence of factors including increasing model complexity – to more fully represent the earth system, increasing spatial resolution - to achieve higher accuracy by resolving fine-scale dynamical to physical, biological, and chemical processes and their interaction, increasing ensemble size - to more accurately represent predictive uncertainty, and increased computing requirements – to enable more accurate and timely weather predictions and climate projections for societal benefit. The belief by many that computing will take care of itself is no longer valid given the disruptive changes in HPC that are driving up the cost of computing, increasing the difficulty of using emerging HPC effectively, and exposing limits in parallelism, portability and scalability of the ESM applications themselves.

We present the first fully curated, publicly accessible archive of infrasonic records from ten large bolide events documented by the U.S. Air Force Technical Applications Center’s global microbarometer network between 1960 and 1972. Captured on analog strip-chart paper, these waveforms predate modern digital arrays and space-based sensors, making them a unique window on meteoroid activity in the mid-twentieth century. Prior studies drew important scientific conclusions from the records but released only limited artifacts, chiefly period–amplitude tables and unprocessed scans, leaving the underlying data inaccessible for independent study. The present release transforms those limited excerpts into a research-ready resource. By capturing ten large events in the mid-20th century, the dataset constitutes a critical reference point for assessing bolide activity before the advent of modern space-based and digital ground-based monitoring. The multi-year coverage and worldwide distribution of events provide a valuable reference for comparing past and more recent detections, facilitating assessments of long-term flux and the dynamics of acoustic wave propagation in Earth’s atmosphere. The dataset’s availability in a consolidated format ensures straightforward access to waveforms and derived measurements, supporting a wide range of scientific inquiries into bolide physics and infrasound monitoring. By preserving these historical acoustic observations, the collection maintains a significant record of mid-20th-century meteoroid entries. It thereby establishes a basis for further refinement of impact hazard evaluations, contributes to historical continuity in atmospheric observation, and enriches the study of meteoroid-generated infrasound signals on a global scale. Dataset: DOI data identification number: 10.7910/DVN/VGSN7Q (direct URL to data: https://doi.org/10.7910/DVN/VGSN7Q). Dataset License: CC-BY-NC

We introduce the Poisson tensor completion (PTC) estimator, a non-parametric differential entropy estimator. The PTC estimator leverages inter-sample relationships to compute a low-rank Poisson tensor decomposition of the frequency histogram. Our crucial observation is that the histogram bins are an instance of a space partitioning of counts and thus can be identified with a spatial Poisson process. The Poisson tensor decomposition leads to a completion of the intensity measure over all bins—including those containing few to no samples—and leads to our proposed PTC differential entropy estimator. A Poisson tensor decomposition models the underlying distribution of the count data and guarantees non-negative estimated values and so can be safely used directly in entropy estimation. Our estimator is the first tensor-based estimator that exploits the underlying spatial Poisson process related to the histogram explicitly when estimating the probability density with low-rank tensor decompositions for the purpose of tensor completion. Furthermore, we demonstrate that our PTC estimator is a substantial improvement over standard histogram-based estimators for sub-Gaussian probability distributions because of the concentration of norm phenomenon.

A series of high entropy AB2-type Ti2-xZrxMnCrFeNi alloys (x = 0.6, 0.7, 0.8, 0.9, 1.0, 1.1 and 1.2) were synthesized to investigate their potential for hydrogen storage and chemical compression. The influence of the Ti/Zr ratio was explored in terms of structural, microstructural and thermodynamic properties. The storage capacity together with the reaction enthalpy and entropy changes of the synthesized high entropy alloys were compared to predictions from Machine Learning (ML) to investigate changes in these properties across the explored composition space. The results revealed that a decreasing Zr content consistently lowered the hydride formation enthalpy and increased the plateau pressure from 8 to >90 bar H2 at 25 °C, in good agreement with ML predictions. Selected compositions (x = 1.0 and 1.2) demonstrated reversible hydrogen storage capability over 150 cycles, with capacities of 1.34–1.40 wt % H2 and remarkable reaction kinetics (<4 min) at ambient temperature. These experimental and computational findings highlight the potential of this Laves-HEA system as tuneable, stable, and cost-effective materials suitable for long-term operations in stationary hydrogen storage and compression applications.

A structural uncertainty validation study for a large-scale, fire-engulfed, elevated object subjected to crosswind is presented to establish the credibility of a high-fidelity, low-Mach, turbulent reacting flow wall-modeled large-eddy simulation (WMLES) approach that includes multiphysics coupling to participating media radiation and conjugate heat transfer. To establish that WMLES can accurately predict surface quantities including drag and pressure coefficient in the low-Mach crosswind regime, a foundational elevated isothermal cylinder validation case is presented at a similar gap-to-diameter ratio of 0.25, spanning the subcritical to supercritical drag regime (ReD=1.1×105 and 4.3×105, respectively). This study exercised both static and dynamic coefficient LES (Smagorinsky and ksgs) with both local and exchange-based velocity sampling. Results showcase that the drag crisis (or the sudden drop in drag coefficient at increased ReD) is well captured when using an exchange-based dynamic coefficient WMLES methodology, while noting lack of mesh convergence and overall drag and pressure coefficient predictively when using a static coefficient, local velocity sampling WMLES. For the O(10) m JP-8 liquid pool fire crosswind validation study presented, two experimental crosswind configurations (2 m/s and 9.5 m/s) are showcased for a fire-engulfed mock fuselage roughly 4 m in diameter. Using the best model-form practices identified in the isothermal study, dynamic coefficient ksgs exchange-based WMLES fire validation findings demonstrate accurate peak irradiation and skin temperature predictions as a function of crosswind magnitude. Excessive yaw in the low-crosswind fuselage configuration, consistent with experimental findings, captured a significant predicted asymmetry in flame attachment and heat flux toward the downwind cylindrical cap - indicative of axial vortex structures transporting the flame along the upper and lower fuselage leeward surface. All fire mesh resolution simulations captured the experimental finding that as crosswind increased, predicted flame shape and peak irradiation magnitude onto the fuselage transitioned from a windward to a leeward cylinder location due to the migration of the upper- to lower-shear fuel/air mixing layer thereby demonstrating the novelty, significance, and credibility of this high-fidelity WMLES reacting flow framework.

Abstract not provided.

DOE maintains an up-to-date documentation of the number of available full drawdowns of each of the caverns at the Strategic Petroleum Reserve (SPR), where an available full drawdown is defined as the ability to move all the oil out from a cavern at full drawdown rates. This report covers impacts on drawdown availability due to SPR operations during Calendar Year 2024.

We present AMR-Wind, a verified and validated high-fidelity computational-fluid-dynamics code for wind farm flows. AMR-Wind is a block-structured, adaptive-mesh, incompressible-flow solver that enables predictive simulations of the atmospheric boundary layer and wind plants. It is a highly scalable code designed for parallel high-performance computing with a specific focus on performance portability for current and future computing architectures, including graphical processing units (GPUs). In this paper, we detail the governing equations, the numerical methods, and the turbine models. Establishing a foundation for the correctness of the code, we present the results of formal verification and validation. The verification studies, which include a novel actuator line test case, indicate that AMR-Wind is spatially and temporally second-order accurate. The validation studies demonstrate that the key physics capabilities implemented in the code, including actuator disk models, actuator line models, turbulence models, and large eddy simulation (LES) models for atmospheric boundary layers, perform well in comparison to reference data from established computational tools and theory. We conclude with a demonstration simulation of a 12-turbine wind farm operating in a turbulent atmospheric boundary layer, detailing computational performance and realistic wake interactions.

Background/Objectives: Children’s biological age does not always correspond to their chronological age. In the case of BMI trajectories, this can appear as phase variation, which can be seen as shift, stretch, or shrinking between trajectories. With maturation thought of as a process moving towards the final state - adult BMI, we assessed whether children can be divided into latent groups reflecting similar maturational age of BMI. The groups were characterised by early factors and time-related features of the trajectories. Subjects/Methods: We used data from two general population birth cohort studies, Northern Finland Birth Cohorts 1966 and 1986 (NFBC1966 and NFBC1986). Height (n = 6329) and weight (n = 6568) measurements were interpolated in 34 shared time points using B-splines, and BMI values were calculated between 3 months to 16 years. Pairwise phase distances of 2999 females and 3163 males were used as a similarity measure in k-medoids clustering. Results: We identified three clusters of trajectories in females and males (Type 1: females, n = 1566, males, n = 1669; Type 2: females, n = 1028, males, n = 973; Type 3: females, n = 405, males, n = 521). Similar distinct timing patterns were identified in males and females. The clusters did not differ by sex, or early growth determinants studied. Conclusions: Trajectory cluster Type 1 reflected to the shape of what is typically illustrated as the childhood BMI trajectory in literature. However, the other two have not been identified previously. Type 2 pattern was more common in the NFBC1966 suggesting a generational shift in BMI maturational patterns.

Abstract not provided.

This work presents a data-driven method for learning low-dimensional time-dependent physics-based surrogate models whose predictions are endowed with uncertainty estimates. We use the operator inference approach to model reduction that poses the problem of learning low-dimensional model terms as a regression of state space data and corresponding time derivatives by minimizing the residual of reduced system equations. Standard operator inference models perform well with accurate training data that are dense in time, but producing stable and accurate models when the state data are noisy and/or sparse in time remains a challenge. Another challenge is the lack of uncertainty estimation for the predictions from the operator inference models. Our approach addresses these challenges by incorporating Gaussian process surrogates into the operator inference framework to (1) probabilistically describe uncertainties in the state predictions and (2) procure analytical time derivative estimates with quantified uncertainties. The formulation leads to a generalized least-squares regression and, ultimately, reduced-order models that are described probabilistically with a closed-form expression for the posterior distribution of the operators. The resulting probabilistic surrogate model propagates uncertainties from the observed state data to reduced-order predictions. We demonstrate the method is effective for constructing low-dimensional models of two nonlinear partial differential equations representing a compressible flow and a nonlinear diffusion–reaction process, as well as for estimating the parameters of a low-dimensional system of nonlinear ordinary differential equations representing compartmental models in epidemiology.

Update slides for CRADA between SNL and GA

Sandia National Laboratories has tested and evaluated a suite of four T120 broadband seismometers designed and manufactured by Nanometrics. Specifically, two T120 Horizon V2 sensors, one T120 Horizon V1 sensor and one T120 Slim Posthole (PH) sensor were evaluated. The purpose of this seismometer evaluation is to measure performance characteristics in areas such as power consumption, sensitivity, frequency response, full scale, self-noise, dynamic range, calibration system response, and passband. The T120 model of sensors are being evaluated to explore the potential for a future seismometer Type Approval process in the International Monitoring System (IMS) of the Comprehensive Nuclear-Test-Ban Treaty (CTBT).

See document.

Abstract not provided.

NPJ Advanced Manufacturing

Abstract not provided.

On October 16-17, 2024, in Long Beach, California, the U.S. DOE Bioenergy Technologies Office (BETO), in collaboration with Sandia National Laboratories, NREL, PNNL, and LANL, convened experts from industry and academia to explore the potential for utilizing existing petroleum refinery infrastructure to expand the bioeconomy. Although refinery integration is not a new concept, recent trends have revitalized interest in this topic.

Carbon fiber provides opportunity to reduce weight in structural composites, including wind turbine blades, due to the material's superior specific stiffness and specific strength compared to alternatives. Despite these advantages, cost and compressive performance are considered weaknesses for carbon fiber products available today. Studies to produce low-cost carbon fiber alternatives, including the use of textile-derived precursor systems, have shown progress and merit through the DOE/ORNL low-cost carbon fiber initiatives. This work focuses on enabling increases in compressive strength through design of the carbon fiber geometry, applicable to both textile and conventional precursor systems, while also providing opportunities to reduce carbon fiber processing costs. Fiber-resin interface and fiber alignment are among the most frequently cited factors controlling composite compressive performance. However, it is believed that there is opportunity in traditionally unexplored routes to increasing compressive strength through alteration of the carbon fiber geometry by increasing the fiber area moment of inertia and/or the fiber perimeter and interfacial area. This paper presents initial results from manufacturing carbon fiber materials to assess the impacts of carbon fiber size on tested composite compressive performance with projected neutral or even beneficial impact on fiber and composite manufacturing economics. Carbon fiber systems with increasing size illustrate a favorable correlation for compressive performance greater than predicted from a micromechanical failure model. The manufacturing and mechanical test results support the hypothesis of this work that alterations to fiber geometry can be used to produce improvements of the compressive strength of carbon fiber reinforced polymers and provide incentive for related work in designing alternative shapes to further enhance compressive performance.

Abstract not provided.

Reconstructing 3D granular microstructures within volumes of arbitrary geometries from limited 2D image data is crucial for predicting the material properties, as well as performances of structural components accounting for material microstructural effects. We present a novel generative learning framework that enables exascale reconstruction of granular microstructures within complex 3D geometric volumes. Building upon existing transfer learning techniques using pre-trained convolutional neural networks (CNN), we introduce several key innovations to overcome the difficulties inherent in arbitrary geometries. Our framework incorporates periodic boundary conditions using circular padding techniques, ensuring continuity and representativeness of the reconstructed microstructures. We also introduce a novel seamless transition reconstruction (STR) method that creates statistically equivalent transition zones to integrate multiple pre-existing 3D microstructure volumes. Based on STR, we propose a cost-effective strategy for reconstructing microstructures within complex geometric volumes, minimizing computational waste. Validation through numerical experiments using kinetic Monte Carlo simulations demonstrates accurate reproduction of grain statistics, including grain size distributions and morphology. A case study involving the reconstruction of a 4-blade propeller microstructure illustrates the method's capability to efficiently handle complex geometries. The proposed framework significantly reduces computational demands while maintaining high reconstruction quality, paving the way for scalable microstructure reconstruction in materials design and analysis.

Abstract not provided.

Parametrized artificial neural networks (ANNs) can be very expressive ansatzes for variational algorithms, reaching state-of-the-art energies on many quantum many-body Hamiltonians. Nevertheless, the training of the ANN can be slow and stymied by the presence of local minima in the parameter landscape. One approach to mitigate this issue is to use parallel tempering methods, and in this work, we focus on the role played by the temperature distribution of the parallel tempering replicas. Using an adaptive method that adjusts the temperatures in order to equate the exchange probability between neighboring replicas, we show that this temperature optimization can significantly increase the success rate of the variational algorithm with negligible computational cost by eliminating bottlenecks in the replicas' random walk. We demonstrate this using two different neural networks, a restricted Boltzmann machine and a feedforward network, which we use to study a toy problem based on a permutation invariant Hamiltonian with a pernicious local minimum and the J1-J2 model on a rectangular lattice.

We present theories and techniques to support safe, predictable, and efficient impact testing at the Sandia National Laboratories TA-III Mechanical Shock Complex (MSC), located in Building 6570. The two testing instruments used at the MSC are the 20-inch actuator, which is a launcher that propels sleds down an indoor track, and an indoor 6-in bore compressed-gas powered cannon. The theories and techniques presented derive from basic mechanical engineering principles but are significantly informe

This work presents a data-driven method for learning low-dimensional time-dependent physics-based surrogate models whose predictions are endowed with uncertainty estimates. We use the operator inference approach to model reduction that poses the problem of learning low-dimensional model terms as a regression of state space data and corresponding time derivatives by minimizing the residual of reduced system equations. Standard operator inference models perform well with accurate training data that are dense in time, but producing stable and accurate models when the state data are noisy and/or sparse in time remains a challenge. Another challenge is the lack of uncertainty estimation for the predictions from the operator inference models. Our approach addresses these challenges by incorporating Gaussian process surrogates into the operator inference framework to (1) probabilistically describe uncertainties in the state predictions and (2) procure analytical time derivative estimates with quantified uncertainties. The formulation leads to a generalized least-squares regression and, ultimately, reduced-order models that are described probabilistically with a closed-form expression for the posterior distribution of the operators. The resulting probabilistic surrogate model propagates uncertainties from the observed state data to reduced-order predictions. We demonstrate the method is effective for constructing low-dimensional models of two nonlinear partial differential equations representing a compressible flow and a nonlinear diffusion–reaction process, as well as for estimating the parameters of a low-dimensional system of nonlinear ordinary differential equations representing compartmental models in epidemiology.

Sputter-deposited, nanocrystalline Cu-Ag thin films produced across a broad compositional and deposition-parameter space were evaluated to unravel the process-structure-property relationships important for creating hard, conductive electrical contacts and coatings. Combinatorial deposition involving pulsed direct current magnetron sputtering of elemental targets enabled swift examination of nearly the full range of alloy compositions and a relevant portion of deposition atomistics. Several high-throughput characterization modalities were employed to evaluate the chemistry, structure, and properties of the films. The resultant hardness, modulus, film density, crystal texture, and resistivity were analyzed in terms of key deposition characteristics (incident atom kinetic energy and incidence angle) predicted by binary-collision, kinematic Monte Carlo simulations. The study revealed improved hardness, parabolic resistivity dependence on composition, and compositional and process dependencies of film tarnishing. The results are discussed in the context of variations in microstructure and film density. Transmission electron microscopy and X-ray diffraction demonstrate several forms of compositional variation including solute segregation to grain boundaries as well as periodic, intragranular compositional modulations. Annealing of a Cu-rich alloy film exhibiting grain boundary segregation showed that this as-deposited, compositional variation is not stable above 100 °C. Finally, the Cu-Ag system is shown to have potential for hard, conductive, tarnish-resistant and room temperature-stable nanocrystalline thin films across the composition space.

Quantum-cascade vertical-external-cavity surface-emitting-lasers (VECSELs) based on disordered amplifying metasurfaces are demonstrated and explored as potential broadband, multi-mode THz sources. The disorder is introduced along one spatial axis of the metasurface by pseudo-randomly varying the width of its resonant ridge antennas. Compared to a quantum-cascade (QC) VECSEL based on a uniform metasurface, the disordered structure supports much more localized transverse modes with reduced spatial overlap within the QC gain material. This localization is hypothesized to facilitate the spatial hole burning of the gain material and, therefore, enable multi-mode lasing, particularly for short cavities on the order of a few wavelengths. Several devices have been fabricated and shown to differ from uniform QC-VECSELs in a few key ways, possessing highly nonlinear light-current characteristics, angle-dependent emission spectra and broadband multi-mode lasing. At most, 17 modes are simultaneously observed, spanning 680 GHz. The number of VECSEL modes is shown to have an inverse relationship with the cavity length, which is attributed to increased diffractive losses in the open plano-plano cavity. As the cavity length is tuned, the device emits over a quasi-continuous band from 3.15 to 3.97 THz or up to 4.24 THz with a longitudinal mode hop.

Asymmetric energy exchange interactions, also known as Hatano-Nelson type couplings, enable the study of non-Hermitian physics and associated phenomena like the non-Hermitian skin effect and exceptional points (EPs). Since these interactions are by definition nonreciprocal, there have been very few options for implementations in integrated photonics. In this work, we show that asymmetric couplings are readily achievable in integrated photonic systems through time-domain dynamic modulation. We experimentally study this concept using a two-resonator photonic molecule produced in a lithium niobate on insulator platform that is electro-optically modulated by rf stimuli. We demonstrate the dynamic tuning of the Hatano-Nelson coupling between the resonators, surpassing the asymmetry that has been achieved in previous work, to reach an EP for the first time. We are additionally able to flip the relative sign of the couplings for opposite directions by going past the EP. Using this capability, we show that the through-chain transport can be configured to exhibit both giant (∼60 dB) optical contrast as well as photonic gyration or nonreciprocal π phase contrast.

Simulation-based approaches to microstructure generation can suffer from a variety of limitations, such as high memory usage, long computational times, and difficulties in generating complex geometries. Generative machine learning models present a way around these issues, but they have previously been limited by the fixed size of their generation area. We present a new microstructure generation methodology leveraging advances in inpainting using denoising diffusion models to overcome this generation area limitation. We show that microstructures generated with the presented methodology are statistically similar to grain structures generated with a kinetic Monte Carlo simulator, SPPARKS.

The influence of the external cathodic area, stress intensity, and electrochemical polarization level on the resultant crack growth rate and crack tip pH and potential was studied for stainless steel 304L in a sodium chloride environment. Before testing, the validity of in situ crack tip property measurements using microprobes and drilled holes was supported through modeling efforts and confirmed with experimental observations. Here, it was found that increasing the R-ratio during corrosion fatigue increased the crack growth rate and decreased the measured crack tip pH of the 304L sample. Decreasing the electrochemical potential (i.e., cathodic polarization) of the bulk sample increased the crack tip pH and stopped measured crack extension under constant K conditions. Conversely, increasing the potential of the bulk sample decreased the crack tip pH and induced cracking under constant K conditions. Decreasing the cathodic area external to the crack increased the measured pH at the crack tip and decreased the crack growth rate under fatigue conditions. For the given configuration, environment, and alloy studied, fatigue crack growth is not sustained by cathodic reactions solely occurring in the crack wake and notch. These results have implications for the ability to predict atmospheric crack growth rates from laboratory, full immersion experiments.

Optimization-based coupling (OBC) is an attractive alternative to traditional Lagrange multiplier approaches in multiple modeling and simulation contexts. However, application of OBC to time-dependent problems has been hindered by the computational cost of finding the stationary points of the associated Lagrangian, which requires primal and adjoint solves. This issue can be mitigated by using OBC in conjunction with computationally efficient reduced order models (ROMs). To demonstrate the potential of this combination, in this paper, we develop an optimization-based ROM-ROM coupling for a transient advection-diffusion transmission problem. We pursue the “optimize-then-reduce” path toward solving the minimization problem at each time step and solve reduced space adjoint system of equations, where the main challenge in this formulation is the generation of adjoint snapshots and reduced bases for the adjoint systems required by the optimizer. One of the main contributions of the paper is a new technique for an efficient adjoint snapshot collection for gradient-based optimizers in the context of optimization-based ROM-ROM couplings. In conclusion, we present numerical studies demonstrating the accuracy of the approach along with comparison between various approaches for selecting a reduced order basis for the adjoint systems, including decay of snapshot energy, average iteration counts, and timings.

The adhesion of nanoparticles to their supports is key to their performance and stability. However, scientific advances in this area have been hampered by the difficulty of experimentally probing adhesion. To date, only a single technique has been developed that can directly measure nanoparticle adhesion, and this technique is inherently limited to monometallic systems. We present a versatile technique for the direct measurement of adhesion for bimetallic nanoparticle systems. This technique combines the spatial resolution of transmission electron microscopy with the force resolution of an atomic force microscope to probe individual, well-characterized nanoparticles. A first study of supported bimetallic nanoparticles provides new insights into the complex impact of alloying on nanoparticle adhesion, explained by charge transfer between constituent metals. The new experimental technique is readily extensible to study other multimetallic nanoparticle systems, including the effects of particle size, shape, and orientation, thus enabling advances in our understanding of nanoparticle physics.

Oxygen vacancies in H f x Z r ( 1 − x ) O 2 (HZO) both contribute to stabilization of the ferroelectric orthorhombic phase and promote leakage pathways that limit the endurance of devices based on the material. For this reason, the defect states of oxygen vacancies were investigated using photoemission electron microscopy (PEEM) and photoluminescence spectroscopy (PL), as their concentration was varied via ex situ laser exposure. Following a controlled oxygen vacancy reduction via visible (2.54 eV) laser dosing of HZO, deep-ultraviolet (DUV, 5.82 eV) PEEM was used to spatially probe the resulting mid-gap defect states and work function. Work function was found to increase monotonically with the laser-induced reduction in oxygen vacancy concentration culminating in a total increase near 70 meV. The change implies a Fermi level shift toward the valence band as the total available electron-filled charge states are reduced with the removal of oxygen vacancies. A reduction in charge states is corroborated by the observed lessening of both photoemission and photoluminescence intensities after laser dosing. The deduced position of the Fermi level is within a band of near-conduction band defect states produced by oxygen vacancies that are linked to endurance limiting leakage currents. Together, these results directly identify the primary role of oxygen vacancies on the defect states in HZO while demonstrating that laser exposure can be used for their modification.

A finite element model of a 60-cell monocrystalline silicon glass-polymer photovoltaic module was simulated with ±1.0 kPa and ±2.4 kPa loads applied to the glass to calculate the deformation under load. Cell-to-cell displacements were used to approximate interconnect strain and stress. A mathematical fatigue cycle life relation was fitted to data for the interconnect material (copper), to generate a life prediction at each interconnect location based on the local stress means, reversal extents, and amplitudes. Interconnect stress was found to be significantly asymmetric about zero despite symmetric positive and negative module loads due to laminate thickness offsets about the neutral plane and the effects of module framing. Cycle life results indicated that interconnect fatigue failure was unlikely to occur over a 30-year lifetime of conservative wind and snow load cycles since the typical cell design feature of leaving some unconstrained length between the cell edge and first solder pad increases the effective gauge length and decreases the stress levels below the material endurance limit. Follow-up analyses found that 3.6 mm and 6.4 mm were the minimum unconstrained lengths required to survive the assumed lifetime of wind and snow cycles, respectively, confirming that typical industrial module constructions with 8–15 mm unconstrained lengths should survive conservatively. Notably, large magnitude, low-cycle snow loading was consistently the limiting factor requiring a longer unconstrained interconnect length. Finally, insights and workflows from this study inform module interconnection design limits for survival against mechanical fatigue in deployment environments.

Five alternative metals are investigated as ohmic contacts to n-GaN including Cr/Au, Mo/Au, Pt/Au, Pd/Au, and Ge/Au. Ti-based contacts are traditionally used for ohmic contacts on n-GaN. However, conventional Ti/Al/Ni/Au metallization is found to be incompatible with a self-aligned process for GaN trench MOSFETs due to wet etch restrictions. Therefore, an alternative metallization is needed that is unreactive to the etch chemistry used in the self-aligned process. Additionally, the contact should remain ohmic following anneal at 900 °C so that contact formation can precede the anneal required for p-dopant activation. Here, in the present work, an n-GaN bilayer, consisting of a thin heavily doped contact layer (n₀ = 1 × 10²⁰ cm⁻³) atop a thick lesser doped layer, is used to demonstrate ohmic contacts of alternative metals with low specific contact resistance and extended thermal budget. Cr/Au ohmic contacts are demonstrated up to anneal temperatures of 800 °C, an increase of 200 °C compared to the highest known reports for Cr/Au contacts on n-GaN. Pt/Au metallization is demonstrated as an ohmic contact to n-GaN for the first time and exhibits true temperature-agnostic behavior up to anneal temperatures of 900 °C with specific contact resistance that is near parity with Ti/Al/Ni/Au. The temperature-agnostic behavior of Pt/Au ohmic contacts on the n-GaN bilayer, in addition to chemical compatibility with the self-aligned process, positions Pt/Au contacts as a key enabling element for self-aligned trench MOSFETs on GaN.

Miniature atomic clocks based on the interrogation of the ground state hyperfine splitting of buffer gas cooled ions confined in radio frequency Paul traps have shown great promise as high precision prototype clocks. We report on the performance of two miniature ion trap vacuum packages after being sealed for as much as 10 years. We find the lifetime of the ions within the trap has increased over time for both traps and can be as long as 50 days. We form two clocks using the two traps and compare their relative frequency instability one with another to demonstrate a short-term instability of 5×10^-13$τ$^-1/2 integrating down to 1×10^-14 after 2 ks of integration. The trapped ion lifetime and clock instability demonstrated by these miniature devices despite only being passively pumped for many years represents a critical advance toward their proliferation in the clock community.

Chapare virus (CHAPV) is an emerging New World arenavirus that is the causative agent of Chapare hemorrhagic fever (CHHF) responsible for recent outbreaks with alarmingly high case fatality rates in Bolivia near the Brazilian border. Here, we describe a nonhuman primate (NHP) model of CHHF infection which represents an essential tool to understand this emerging biological threat agent. Cynomolgus macaques challenged intravenously with CHAPV develop clinical disease, which recapitulates several key features of human CHHF. All subjects lost weight and had clinical scores following CHAPV challenge. Notably, one of four NHPs developed lethal disease with viral hepatitis and hemorrhagic features. Clinical chemistry and hematology revealed leukopenia, anemia, thrombocytopenia, and increased transaminase levels. In all four subjects, viremia was detectable for the first week following challenge and viral RNA was detectable in serum and many tissues persisting 35 days-post challenge. Several medical countermeasures (MCM) have efficacy against CHAPV infection in vitro, but the current model for MCM testing and approval of new drugs is reliant on the availability of animal models. This work lays the foundation for future CHHF MCM development.

NasGen provides a path for migration of structural models from Nastran bulk data format (BDF) into both an Exodus mesh file and an ASCII input file for Sierra Structural Dynamics (Salinas) and Solid Mechanics (Adagio). Many tools at Sandia National Labs (SNL) use the Exodus format. This document describes capabilities and limitations of the NasGen translation software.

Abstract not provided.

Brittle behavior of metal alloys is often critical to modeling ballistic impact and penetration. The ALEGRA multiphysics finite element software incorporates calibrated models for the equation of state, elasticity, yield stress, plasticity and fracture, but simulations do not always capture expected metal fracture. Here we report concerted efforts to do so for one important case where experiments clearly show shear fractures: a tungsten sphere impacting a steel plate at various angles. Our best simulations show fractures that are qualitatively similar to experiments, but there are significant differences in quantitative metrics. Specifically, velocities of tracers used to quantify simulated plug parameters consistently fall short of measured plug velocities. Also, simulated plugs break apart more than expected from experimental evidence. We attribute these shortfalls to the lack of an explicit shear fracture mechanism in the material models, leading to over-estimated resistance to plug formation and movement.

Task 1 – Extract GADRAS-DRF capabilities for a workflow intuitive to safeguards analysis. This task is almost complete. Significant improvements have been made to the functionality of customizing peak fits for use in the isotopics application, as well as peak-based FSA model analysis. The remaining tasks revolve around bug fixes, testing the application, and adding a density scroll bar to isotopics for real time analysis updates. The density scroll bar values will be reflected in summary tables and the self-shielding form accessed within isotopics.

This document details a data mining exercise that resulted in an exploratory dataset of publicly reported foreign (non-US) hypersonic vehicle test events. Using a combination of targeted English language searches and country-specific queries, the study aggregates information from digital news media, official press releases, and social media posts. The resulting list of events captures the publicly available accounts of foreign hypersonic tests, although it does not represent an exhaustive record. Limitations such as inconsistent reporting, translation challenges, and the inherently provisional nature of open-source data are acknowledged. This dataset serves as an initial reference point for further inquiries into high-speed atmospheric phenomena and may facilitate future efforts to correlate these events with geophysical measurements.

This document summarizes the key processes (thermal, hydrological, mechanical, and chemical; THMC) impacting the features of a deep geological repository for radioactive waste in salt. Some processes are natural and on-going whether the repository is there or not, and other processes are driven by the pertur- bation associated with the repository. The features considered here include both engineered and natural components of the repository system.

Cross-linking a blend of linear and ring polymers creates a new topology-based dual-network elastomer in which the two components differ significantly in their topology. We use molecular simulations and topological analysis to examine key mechanical properties as functions of ring polymer volume fraction Φ_R. For Φ_R < Φ_R*, where the rings begin to overlap, the network shear modulus G and the maximum stretch ratio λ_p are weakly dependent on Φ_R. For Φ_R > Φ_R*, entanglements trapped in the network are diluted as the rings overlap, leading to a significant decrease in G and an increase in λ_p with increasing Φ_R. Here, the peak tensile stress, σ_p, exhibits a maximum around Φ_R*, indicating an enhancement of network strength due to the stronger cohesion from the entanglements between linear and ring polymers.

Autonomous manipulation is a challenging problem in field robotics due to uncertainty in object properties, constraints, and coupling phenomenon with robot control systems. Humans learn motion primitives over time to effectively interact with the environment. We postulate that autonomous manipulation can be enabled by basic sets of motion primitives as well, but do not necessitate mimicking human motion primitives. This work presents an approach to generalized optimal motion primitives using physics-informed neural networks. Our simulated and experimental results demonstrate that optimality is notionally maintained where the mean maximum observed final position percent error was 0.564% and the average mean error for all the trajectories was 1.53%. These results indicate that notional generalization is attained using a physics-informed neural network approach that enables near optimal real-time adaptation of primitive motion profiles.

The research field of soil viral ecology continues to advance rapidly as the roles of viruses in the functioning of soil ecosystems are increasingly recognized. To address recent developments in the field, the second International Soil Virus Conference was held in Livermore, California, USA, from June 25 to 27th, 2024, providing soil viral ecologists the opportunity to share new findings and suggest guidelines for future research, while encouraging international scientific discussion and collaboration. The meeting was held in person with sessions simultaneously streamed online. Fifty researchers attended from ten different countries and spanned a wide range of subfields and career stages. A total of 21 oral presentations were presented, followed by discussions covering key themes in soil viral research. This report summarizes the main takeaways and recommendations from the talks and discussions.

The Fusion of Simulation, Experiment, and Data (FuSED) team provides a set of tools for solving inverse problems in structural dynamics and thermal physics, and also sensor placement optimization via Optimal Experimental Design (OED). These methods are used for designing experiments, model calibration, and verification/validation analysis of systems. This document provides a user’s guide to the input for the three apps that are supported for these methods. Details of input specifications, output options, and optimization parameters are included.

This report delineates an initial systematic data mining investigation aimed at compiling a preliminary database of re-entry events recorded between 1995 and 2024. The dataset encompasses a heterogeneous assortment of occurrences, including controlled space shuttle re-entries, scheduled orbital de-orbits, and the uncontrolled re-entry of space debris, among other unplanned events. By employing targeted English-language queries, country-specific search parameters, and a suite of technical keywords, alongside a methodology analogous to our prior study on foreign hypersonic test events, the investigation extracts pertinent information from digital news outlets, official press releases, and social media platforms. The resulting compilation serves as a robust empirical foundation and establishes critical ground truth for the study of high-speed atmospheric dynamics. Although derived solely from publicly available sources and acknowledged as non-exhaustive, this initial study paves the way for the development of a more comprehensive database. It motivates subsequent efforts in data fusion and multi-modality sensing investigations aimed at correlating geophysical signatures with these events.

Using a belt as a replacement for a rope on a rotary power take-offs (PTOs) system has become more common for wave energy converters, improving cyclic bend over sheave performance with a smaller bending thickness for belts. However, the service life predictions of PTOs are a major concern in design, because belt performance under harsh underwater environments is largely less studied. In this work, the effect of fleet and twist angles on wear life is being investigated both experimentally and numerically. Two three-dimensional equivalent static finite element models are constructed to evaluate the complex stress state of polyurethane-steel belts around steel drums. The first is to capture the response of the experimental investigation performed on the wear life, and the second to predict the wear life of an existing functional PTO. The results show a significant effect for fleet and twist angles on stress concentrations and estimated service life.

Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high fidelity, validated models used in modal, vibration, static and shock analysis of structural systems. This manual describes the theory behind many of the constructs in Sierra/SD. For a more detailed description of how to use Sierra/SD, we refer the reader to User’s Manual. Many of the constructs in Sierra/SD are pulled directly from published material. Where possible, these materials are referenced herein. However, certain functions in Sierra/SD are specific to our implementation. We try to be far more complete in those areas. The theory manual was developed from several sources including general notes, a programmer_notes manual, the user’s notes and of course the material in the open literature.

Development of a defensible source-term model (STM), usual ly a thermodynamical model for radionuclide solubility calculations, is critical to a performance assessment (PA) of a geologic repository for nuclear waste disposal. Such a model is generally subjected to rigorous regulatory scrutiny. In this article, we highlight key guiding principles for STM model development and validation in nuclear waste management. We illustrate these principles by closely examining three recently developed thermodynamic models with the Pitzer formulism for aqueous H+—Nd3+—NO3−(—oxalate) systems in a reverse alphabetical order of the authors: the XW model developed by Xiong and Wang, the OWC model developed by Oakes et al., and the GLC model developed by Guignot et al., among which the XW model deals with trace activity coefficients for Nd(III), while the OWC and GLC models are for concentrated Nd(NO3)3 electrolyte solutions. The principles highlighted include the following: (1) Principle 1. Validation against independent experimental data: A model should be validated against experimental data or field observations that have not been used in the original model parameterization. We tested the XW model against multiple independent experimental data sets including electromotive force (EMF), solubility, water vapor, and water activity measurements. The results show that the XW model is accurate and valid for its intended use for predicting trace activity coefficients and therefore Nd solubility in repository environments. (2) Principle 2. Testing for relevant and sensitive variables: Solution pH is such a variable for an STM and easily acquirable. All three models are checked for their ability to predict pH conditions in Nd(NO3)3 electrolyte solutions. The OWC model fails to provide a reasonable estimate for solution pH conditions, thus casting serious doubt on its validity for a source-term calculation. In contrast, both the XW and GLC models predict close-to-neutral pH values, in agreement with experimental measurements. (3) Principle 3. Honoring physical constraints: Upon close examination, it is found that the Nd(III)-NO3 association schema in the OWC model suffers from two shortcomings. Firstly, its second stepwise stability constant for Nd(NO3)2+ (log K2) is much higher than the first stepwise stability constant for NdNO32+ (log K1), thus violating the general rule of (log K2–log K1) < 0, or (Formula presented.). Secondly, the OWC model predicts abnormally high activity coefficients for Nd(NO3)2+ (up to ~900) as the concentration increases. (4) Principle 4. Minimizing degrees of freedom for model fitting: The OWC model with nine fitted parameters is compared with the GLC model with five fitted parameters, as both models apply to the concentrated region for Nd(NO3)3 electrolyte solutions. The latter appears superior to the former because the latter can fit osmotic coefficient data equally well with fewer model parameters. The work presented here thus illustrates the salient points of geochemical model development, selection, and validation in nuclear waste management.

Lithium primary batteries (LPBs) remain essential in critical applications such as military, aerospace, medical and emergency devices, and portable electronics. Their superior energy density over lithium-ion batteries offers a significant advantage for long-duration use. Therefore, accurate estimation of the state of charge (SoC) is essential for ensuring the reliable and safe operation of these batteries. While extensive research has been conducted on SoC estimation techniques for lithium-ion secondary batteries, LPBs present unique challenges that complicate accurate SoC estimation. Moreover, research on nondestructive testing techniques for SoC estimation in LPBs is significantly lacking. In this review article, it is aimed to provide a comprehensive overview of recent advancements in SoC estimation for LPBs and generates new insights and directions for future research. Herein, existing methods are discussed and their effectiveness and mechanisms are identified, and areas for further optimization are outlined. More theoretical/experimental efforts to advance SoC detection in LPBs is recommended due to challenges identified with existing techniques.

Tamper-indicating devices (TIDs), also known as seals, play a crucial role in various sectors including international nuclear safeguards, arms control, domestic security, and commercial products, by ensuring that monitored or high-value items are not accessed undetected. These devices do not block access but alert to unauthorized tampering. With adversaries' capabilities evolving, there's a pressing need for seals to advance in terms of effectiveness (e.g., better tamper indication and unique identification), and new technology can improve the efficiency of installation and verification. Passive loop seals, widely used in international nuclear safeguards to ensure that continuity of knowledge is maintained on declared items, face stringent International Atomic Energy Agency (IAEA) requirements that surpass those met by commercial products. The metal cup seal (Figure 1, left), a staple IAEA seal, is robust but requires significant resources for post-use verification – specifically, the seal’s unique identity can only be verified at IAEA headquarters after removal from facilities. Further, the seal has been in use for decades and seal types should periodically be replaced to counter adversarial efforts for defeating seals. In 2020, the IAEA outlined about 40 requirements for a new passive loop seal, aiming for in-situ verification, minimal external tool use, unique identification (UID), and clear tamper indication. In response, research and development efforts focused on creating a new passive loop seal that meets these criteria and in 2022 the IAEA announced the completion of the Field Verifiable Passive Loop Seal (FVPS) (Figure 1, right). Concurrently to the IAEA’s efforts, Sandia National Laboratories (SNL) and Oak Ridge National Laboratory (ORNL) designed, developed, and tested two seal versions – Puck and Puck/SAW, with Puck based on the IAEA’s requirements and including a novel visually-obvious tamper response, and Puck/SAW adding additional beneficial capabilities like the ability to receive a unique identifier from a standoff distance and monitoring the wire integrity. Puck/SAW was specifically designed and developed to address sealing applications in dry spent fuel storage facilities, where the number of sealed spent fuel containers results in heavy verification burden and inspector safety issues related to radiation exposure. These efforts are described in this Executive Summary.

Abstract not provided.

Ever-increasing wind turbine size has challenged predictive capabilities on several fronts. To address part of the blade structural modeling uncertainty, a systematic model fidelity comparison study was conducted on commonly used finite elements. pyNuMAD was utilized to create beam, shell, and solid models of a 100 m long blade undergoing large static deflections. The solid model avoided the use of layered-solid elements by resolving core and facesheet layers. An unprecedented model with 73.7 million elements revealed insights that have never been possible from prior experimental and numerical studies. As compared to the solid element model, the tip deflection from the shell and beam model was found to be about 2% and 4.3% too low, respectively. The twist from the beam model was found to be about 5.6% too high, while the twist from shell model was 24% too low, though improvement was demonstrated with mesh refinement. The beam model adhesive stresses were more accurate than the shell model. Out-of-plane stresses were of great significance near geometric and material discontinuities, and neither the shell nor beam model captured these effects well. Failure predictions from beam, shell, or layered-solid models are unlikely to be reliable at trailing edges, adhesives, ply-drops, spar-cap boundaries.

Abstract not provided.

Public-facing solar hosting capacity (HC) maps, which show the maximum amount of solar energy that can be installed at a location without adverse effects, have proven to be a key driver of solar soft cost reductions through a variety of pathways (e.g., streamlining interconnection, siting, and customer acquisition processes). However, current methods for generating HC maps require detailed grid models and time-consuming simulations that limit both their accuracy and scalability—today, only a handful out of almost 2,000 utilities provide these maps. This project developed and validated data-driven algorithms for calculating solar HC using data from AMI without the need of detailed grid models or simulations. The algorithms were validated on utility datasets and incorporated as an application into NRECA’s Open Modeling Framework (OMF.coop) for the over 260 coops and vendors throughout the US to use. The OMF is free and open-source for everyone.

Recent successes in the exploration, drilling, and discovery of geologic hydrogen have generated notable excitement. This new energy resource has the potential to make an important contribution to our nation’s energy supply, resiliency, and security. Contemporary studies of geologic hydrogen have a common theme of suggesting places where it might be found or even more specifically, what rocks in what geologic formations may contribute to its formation — either naturally or via artificially induced means. This vital ongoing body of work sets the stage for imagining what may be possible with vast available quantities of naturally occurring hydrogen in the subsurface. While acknowledging current approaches to characterizing geologic hydrogen, this report advances the discussion by suggesting next steps, including the critical science and engineering necessary to make geologic hydrogen an affordable and reliable part of the U.S. energy portfolio.

Two bolides (2016 June 2 and 2019 April 4) were detected at multiple regional infrasound stations, with many of the locations receiving multiple detections. Analysis of the received signals was used to estimate the yield, location, and trajectory, as well as the type of shock that produced the received signal. The results from the infrasound analysis were compared with ground-truth information that was collected through other sensing modalities. This multimodal framework offers an expanded perspective on the processes governing bolide shock generation and propagation. The majority of signal features showed reasonable agreement between the infrasound-based interpretation and the other observational modalities, though the yield estimate from the 2019 bolide was significantly lower using the infrasound detections. There was also evidence suggesting that one of the detections was from a cylindrical shock that was initially propagating upward, which is unusual though not impossible.

The HyRAM+ software is an open-source toolkit that provides publicly available models and default input values to enable straightforward and consistent safety assessments for hydrogen and other alternative fuel systems, such as natural gas and propane. The HyRAM+ quantitative risk assessment calculation incorporates annual likelihood of leaks or failures for both compressed gaseous and liquefied flammable fuels, as well as probabilistic models for the effects of heat flux and overpressure. HyRAM

This report documents results assessing the presence of Naturally Occurring Radioactive Material (NORM) in illicit drugs using high-resolution gamma spectroscopy. Specifically, high-purity germanium (HPGe) measurements were conducted at the Drug Enforcement Administration (DEA) South Central Laboratory on 11 samples (referred to as exhibits) including cocaine, fentanyl, heroin, ketamine, LSD, marijuana, methamphetamines, psilocybin mushrooms, and synthetic cannabinoids.

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options an

Public-facing solar hosting capacity (HC) maps, which show the maximum amount of solar energy that can be installed at a location without adverse effects, have proven to be a key driver of solar soft cost reductions through a variety of pathways (e.g., streamlining interconnection, siting, and customer acquisition processes). However, current methods for generating HC maps require detailed grid models and time-consuming simulations that limit both their accuracy and scalability—today, only a handful out of almost 2,000 utilities provide these maps. This project developed and validated data-driven algorithms for calculating solar HC using data from AMI without the need of detailed grid models or simulations. The algorithms were validated on utility datasets and incorporated as an application into NRECA’s Open Modeling Framework (OMF.coop) for the over 260 coops and vendors throughout the US to use. The OMF is free and open-source for everyone.

Abstract not provided.

This is the Abstract at UUR Level for the full Airborne Circumvent & Recover Demonstration for Strategic Environments paper which is CUI in a separate release.

A mesoscale model to predict helium bubble evolution is needed for tritium applications. Such a model requires that the conventional kinetic Monte Carlo (kMC) simulations be significantly accelerated. The objective of this report is to (a) highlight the concepts and mathematical expressions of the accelerated method for defect implementation that have not been published, (b) show an example input file to run the kMC code, and (c) provide suggestions on future improvement following my retirement.

The size and attribution of the regional net carbon flux from land-use change (LUC) activities (ELUC) are often highly debated, especially in regions such as China, which has experienced decades-long extensive reforestation activities. Here, using a LUC dataset incorporating remote-sensing and national forest inventory data with two modelling approaches, we show that ELUC in China shifted from a carbon source to a sink in the 1990s, contributing to a net cumulative CO2 removal of 2.0 Pg C during 1981–2020. From 2001 to 2020, the average ELUC was −0.14 Pg C yr−1, accounting for over one-third of the national land carbon sinks. Forest-related LUC activities contributed greatly to national carbon fluxes, while non-forest-related activities played a dominant role in certain areas. Our findings suggest that the carbon sinks from LUC activities in China may be largely underestimated in global assessments, underscoring the need to develop region-specific modelling for evaluation and potential regulation.

Abstract not provided.

This work aimed to apply Sandia’s expertise in metallurgy and modeling to enable the use of hybrid laser arc welding for building nuclear reactor containment structures, via a collaboration with Holtec International. Experimental observations were coupled with finite element analysis to resolve microstructure development, mechanical properties, distortion, and residual stress in welds relevant to the production of the Holtec SMR-160. High residual stresses were observed in welds that were not subjected to preheat. Meanwhile, the microstructure of the welds generally exhibited a narrow heat affected zone relative to conventional arc welds. FEA appeared to be effective in simulating the thermal/mechanical conditions that occur during hybrid laser arc welding of simplified and instrumented test welds. Subsequently, FEA was used to perform sensitivity analyses for various weld geometries that would be prohibitively costly to assess with physical experiments. Insights from the study were used to inform Holtec’s welding process, and successful production welds were performed in 2025.

Tamper-indicating devices (TIDs), also known as seals, play a crucial role in various sectors including international nuclear safeguards, arms control, domestic security, and commercial products, by ensuring that monitored or high-value items are not accessed undetected. These devices do not block access but alert to unauthorized tampering. With adversaries' capabilities evolving, there's a pressing need for seals to advance in terms of effectiveness (e.g., better tamper indication and unique identification), and new technology can improve the efficiency of installation and verification. Passive loop seals, widely used in international nuclear safeguards to ensure that continuity of knowledge is maintained on declared items, face stringent International Atomic Energy Agency (IAEA) requirements that surpass those met by commercial products. The metal cup seal (Figure 1, left), a staple IAEA seal, is robust but requires significant resources for post-use verification – specifically, the seal’s unique identity can only be verified at IAEA headquarters after removal from facilities. Further, the seal has been in use for decades and seal types should periodically be replaced to counter adversarial efforts for defeating seals. In 2020, the IAEA outlined about 40 requirements for a new passive loop seal, aiming for in-situ verification, minimal external tool use, unique identification (UID), and clear tamper indication. In response, research and development efforts focused on creating a new passive loop seal that meets these criteria and in 2022 the IAEA announced the completion of the Field Verifiable Passive Loop Seal (FVPS) (Figure 1, right). Concurrently to the IAEA’s efforts, Sandia National Laboratories (SNL) and Oak Ridge National Laboratory (ORNL) designed, developed, and tested two seal versions – Puck and Puck/SAW, with Puck based on the IAEA’s requirements and including a novel visually-obvious tamper response, and Puck/SAW adding additional beneficial capabilities like the ability to receive a unique identifier from a standoff distance and monitoring the wire integrity. Puck/SAW was specifically designed and developed to address sealing applications in dry spent fuel storage facilities, where the number of sealed spent fuel containers results in heavy verification burden and inspector safety issues related to radiation exposure. These efforts are described in this Executive Summary.

Multiple scattering of light serves as a mechanism for feedback in random lasers. Consequently, internal spatial mode patterns, lasing wavelengths, and output directionality can all be random. Strong mode interaction can occur in such devices due to spatially overlapping modes resulting in nonlinearity with respect to the pump input power. Nevertheless, temporal coherence and lasing mode amplitude can be fixed at a constant pumping rate. This is a property desirable for applications where unique randomness is exploited but expected to be reliable over time, such as physical unclonable functions. Random lasers can also be cheaply and easily fabricated, exhibit relatively low lasing thresholds and high emission intensity. However, the precise scattering properties of such structures and fluctuations in the pump field can make device emission irreproducible, thereby limiting random laser applications. Here we directly compare the random lasing spectra from zinc oxide samples fabricated in four distinct ways: spin-coating, sputtering, solgel deposition, and atomic layer deposition. The particular method of fabrication has a strong impact. Samples made through atomic layer deposition here exhibit both reproducibility and strong nonlinearity desirable for applications. Randomness in emission spectra persists across hundreds of repeated and averaged measurements irrespective of spatial location and is demonstrably nonlinear with respect to input signal intensity.

Photovoltaic (PV) systems are essential for the transition to sustainable energy, reducing fossil fuel dependence and mitigating climate change. Although PV requires minimal land area — PV can meet the European Union's energy needs using only 0.26% of its land — space for deployment is often scarce in densely populated regions. Floating photovoltaics (FPV) offer an effective solution to land-use challenges by installing PV systems on floating structures in water bodies. FPV is a growing niche within PV with a cumulative installed capacity reaching 7.7 GW globally by 2023. Almost 90% of the installed FPV capacity is in Asia, with close to 50% of in China alone, while the Netherlands and France are the largest markets outside Asia. FPV shows strong potential to support climate targets, but still faces challenges like regulatory barriers, cost competitiveness compared to ground-based PV (GPV), and uncertainties about environmental impacts and system reliability. FPV systems are currently installed mainly on sheltered inland waters, such as quarry lakes, irrigation ponds and reservoirs. FPV technical standards are still being developed. Guidelines have been published by the World Bank, DNV, and Solar Power Europe, and emerging national standards from South Korea, China, and Singapore address design, components, and safety. The International Electrotechnical Commission (IEC) is working on formal standards for floats, mooring systems, and electrical connectors. However, the published best practices lack quantitative guidance for yield modelling and reliability, which this report aims to address. It provides data-driven insights, models, and parameters essential for accurate energy yield, reliability, and maintenance predictions over FPV systems' lifetimes.

Tests from the Sierra Structural Dynamics verification test suite are reviewed. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the Sierra code results to the analytic solution is provided. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems.

Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high-fidelity, validated models used in modal, vibration, static and shock analysis of weapons systems. This document provides a user’s guide to the input for Sierra/SD. Details of input specifications for the different solution types, output options, element types and parameters are included. The appendices contain detailed examples, and instructions for running the software on parallel platforms.

Describes the set of example problems, input decks, and meshes that are distributed with Sierra/SD.

One of the most striking measurements taken during DOE’s EGS Collab project at the 4850-foot depth location was the so-called ‘sewer cam’, which enabled direct visualization of the flow of water into the production well through fractures during the stimulation. The ability to see directly which fractures were flowing and (roughly) how much was a breakthrough in understanding the topology of the created fracture network. Achieving this kind of fracture flow imaging at FORGE would be more challenging because of the 225°C temperature, but equally or even more valuable if it could be achieved. In 2017, a joint project between Sandia and Stanford developed a downhole tool concept to measure the enthalpy of multiphase fluid entering a geothermal well from individual fractures (Gao et al., 2017). For the FORGE project, measuring enthalpy is of less interest because the fluid is expected to be single-phase liquid water. However, the foundation of the device was the measurement of chloride ion concentration, which could form the basis for a direct measurement of inflow from fractures. During the 2017 project, this novel chloride sensing system was implemented into a laboratory test instrument, and we confirmed the capability of the system to measure the ion concentration of fluid entering a model wellbore through a small entry port. The wellbore was a 6-inch diameter model well, and the port was approximately 0.08 inch (2mm) in diameter. The device could measure the chloride concentration accurately even when the well was flowing in a bubbly flow. Given its accuracy, the tool should be able to identify locations of water entering the wellbore even if the ion concentration differs only slightly from that of the water in the well. It is likely that different fractures may flow slightly different chloride concentrations, which would make it feasible to detect individual fractures as well as to estimate the volume of their flow. Ultimately, we could also recognize different fractures flowing back significantly different ion concentrations after fracturing in the FORGE wells. This could be realized by adding different ions in the fracturing fluids in different fractures created at different stages of stimulation (and modifying the tool to include different ion specificity). Sandia’s tool was shown during the study to have the capability to withstand the 225°C temperature, and the electrochemical sensing elements were tested in the laboratory to 225°C at 1500 psia for 24 hours. An early implementation of the fully integrated downhole electrochemical tool, including high-temperature electronics, robust housing, and wireline truck interface, had previously been constructed and tested successfully at Sandia; thus, hardware development tasks focused on advancing the technology readiness level (TRL) of this promising technology for FORGE deployment, rather than on developing a new scientific basis for its operation. The data collection electronics in this tool allowed for several other sensors (pressure, temperature, flow spinner) to be implemented in parallel as well. The research was a new collaboration between Stanford and Sandia to modify and refine the tool for FORGE deployment, to make the downhole measurements, and to characterize the evolving fractures.

Brazing and soldering are metallurgical joining techniques that use a wetting molten metal to create a joint between two faying surfaces. The quality of the brazing process depends strongly on the wetting properties of the molten filler metal, namely the surface tension and contact angle, and the resulting joint can be susceptible to various defects, such as run-out and underfill, if the material properties or joining conditions are not suitable. In this work, we implement a finite element simulation to predict the formation of such defects in braze processes. This model incorporates both fluid–structure interaction through an arbitrary Eulerian–Lagrangian technique and free surface wetting through conformal decomposition finite element modeling. Upon validating our numerical simulations against experimental run-out studies on a silver-Kovar system, we then use the model to predict run-out and underfill in systems with variable surface tension, contact angles, and applied pressure. Finally, we consider variable joint/surface geometries and show how different geometrical configurations can help to mitigate run-out. This work aims to understand how brazing defects arise and validate a coupled wetting and fluid–structure interaction simulation that can be used for other industrial problems.

A striking example of frustration in physics is Hofstadter's butterfly, a fractal structure that emerges from the competition between a crystal's lattice periodicity and the magnetic length of an applied field. Current methods for predicting the topological invariants associated with Hofstadter's butterfly are challenging or impossible to apply to a range of materials, including those that are disordered or lack a bulk spectral gap. Here, we demonstrate a framework for predicting a material's local Chern markers using its position-space description and validate it against experimental observations of quantum transport in artificial graphene in a semiconductor heterostructure, inherently accounting for fabrication disorder strong enough to close the bulk spectral gap. By resolving local changes in the system's topology, we reveal the topological origins of antidot-localized states that appear in artificial graphene in the presence of a magnetic field. Moreover, we show the breadth of this framework by simulating how Hofstadter's butterfly emerges from an initially unpatterned 2D electron gas as the system's potential strength is increased and predict that artificial graphene becomes a topological insulator at the critical magnetic field. Overall, we anticipate that a position-space approach to determine a material's Chern invariant without requiring prior knowledge of its occupied states or bulk spectral gaps will enable a broad array of fundamental inquiries and provide a novel route to material discovery, especially in metallic, aperiodic, and disordered systems.

Non-invasive imaging of microwave (MW) magnetic fields with microscale lateral resolution is pivotal for various applications, such as MW technologies and integrated circuit failure analysis. Diamond nitrogen-vacancy (NV) center magnetometry has emerged as an ideal tool, offering micrometer-scale resolution, millimeter-scale field of view, high sensitivity, and non-invasive imaging compatible with diverse samples. However, up until now, it has been predominantly used for imaging of static or low-frequency magnetic fields or, concerning MW field imaging, to directly characterize the same microwave device used to drive the NV spin transitions. In this work, we leverage an NV center ensemble in diamond for wide-field imaging of MW magnetic fields generated by a test device employing a differential measurement protocol. The microscope is equipped with a MW loop to induce Rabi oscillations between NV spin states, and the MW field from the device-under-test is measured through local deviations in the Rabi frequency. This differential protocol yields magnetic field maps of a 2.57 GHz MW field with a sensitivity of ∼ 9 μ T Hz^{− 1/2} for a total measurement duration of T = 357 s, covering a 340 × 340 μ m² field of view with a micrometer-scale spatial resolution and a device-under-test input power dynamic range of 30 dB. This work demonstrates a novel NV magnetometry protocol, based on differential Rabi frequency measurement, that extends NV wide-field imaging capabilities to imaging of weak MW magnetic fields that would be difficult to measure directly through standard NV Rabi magnetometry.

The Integrated Tiger Series (ITS) generates a database containing energy deposition data. This data, when stored on an Exodus file, is not typically suitable for analysis within Sierra Mechanics for finite element analysis. The its2sierra tool maps data from the ITS database to the Sierra database. This document provides information on the usage of its2sierra.

Global wildfire events have had increasingly severe impacts in recent years, particularly in the western USA, driven by extreme fire-weather conditions, fuel accumulation and multiple ignition sources. Wildfires sparked by power lines tend to be larger and more destructive, as they often occur during high winds, which accelerate the spread of fires. Moreover, efforts to contain wildfires frequently result in power outages, causing considerable economic disruption. Here, in this Review, we examine wildfire risks related to power-line-induced ignitions, infrastructure damage, climate-induced environmental impacts, grid operational risks, real-time grid management risks, vegetation management risks, and financial and funding risks in the context of a changing climate and their interdependence with power grid infrastructures. We then explore the resilience of power grids under wildfire threats, looking at risk analysis, prediction and mitigation strategies. The Review also shares practical insights and experiences in the USA to inform researchers, policymakers and industry professionals.

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 codesuite, and the results of the test are checked versus the correct analytical result. For each of thetests presented in this document, the test setup, a description of the analytic solution, andcomparison of the Sierra/SM code results to the analytic solution is provided. Mesh convergenceis also check

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. Sierra/PMR handles the participating-media thermal radiation mechanics. This project is an integral part of the SIERRA multi-mechanics software development project. Fuego depends heavily upon the core architecture developments provided by SIERRA for massively parallel computing, solution adaptivity, and mechanics coupling on unstructured grids.

This is an addendum to the Sierra/SolidMechanics 5.24 User’s Guide that documents additional capabilities available only in alternate versions of the Sierra/SolidMechanics (Sierra/SM) code. These alternate versions are enhanced to provide capabilities that are regulated under the U.S. Department of State’s International Traffic in Arms Regulations (ITAR) export control rules. The ITAR regulated codes are only distributed to entities that comply with the ITAR export control requirements. The ITAR enhancements to Sierra/SM include material models with an energy-dependent pressure response (appropriate for very large deformations and strain rates) and capabilities for blast modeling. This document is an addendum only; the standard Sierra/SolidMechanics 5.24 User’s Guide should be referenced for most general descriptions of code capability and use.

Silicon-based microelectronics are limited to ∼150°C and therefore not suitable for the extremely high temperatures in aerospace, energy, and space applications. While wide-band-gap semiconductors can provide high-temperature logic, nonvolatile memory devices at high temperatures have been challenging. In this work, we develop a nonvolatile electrochemical memory cell that stores and retains analog and digital information at temperatures as high as 600°C. Through correlative scanning transmission electron microscopy, we show that this high-temperature information retention is a result of composition phase separation between the oxidized and reduced forms of amorphous tantalum oxide. This result demonstrates a memory concept that is resilient at extreme temperatures and reveals phase separation as the principal mechanism that enables nonvolatile information storage in these electrochemical memory cells.

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. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems.

The design and realization of highly selective nanoporous materials are necessary to target critical separations across industries. By leveraging pore size, pore shape, and linker functionalization, the design of nanoporous solid adsorbents will enable the rapid production of energy efficient separation materials for high-value gas mixtures. This study uses a combination of modeling, synthesis, and gas adsorption testing to investigate a new class of small-pore isostructural rare-earth (RE) 2,5-dihydroxyterephthalic acid (DOBDC) metal-organic frameworks (MOFs) (RE: Pr-, Gd-, Er-, Yb; DOBDC = 2,5-dihydroxyterephthalic acid) and their adsorption selectivity for acetylene/ethylene mixtures. Density functional theory simulations identified that selective binding of acetylene over ethylene in the Gd-, Er-, and Yb-DOBDC MOFs was due to hydrogen-bonding between acetylene and the linker hydroxyl. Adsorption experiments validated the computational results by identifying mechanisms that control the acetylene/ethylene adsorption selectivity and high acetylene adsorption. Furthermore, dynamic column breakthrough experiments with the Gd-DOBDC MOF validated the simulations and indicated that ethylene can be separated from acetylene in a mixture containing 1 vol % acetylene and 39 vol % ethylene (balance argon). The results highlight the complexity of gas binding in functional porous materials and how combining modeling and experiment enables a fundamental understanding of gas-framework interactions that can be leveraged for the design of future separation materials.

Josephson junctions are the principal circuit element in numerous superconducting quantum information devices and can be readily integrated into large-scale electronics. However, device integration at the wafer scale necessarily depends on having a reliable, high-fidelity, and high-yield fabrication method for creating Josephson junctions. When creating Al/AlO_x based superconducting qubits, the standard Josephson junction fabrication method relies on a sub-micron suspended resist bridge, known as a Dolan bridge, which tends to be particularly fragile and can often times fracture during the resist development process, ultimately resulting in device failure. In this work, we demonstrate a unique Josephson junction lithography mask design that incorporates stress-relief channels. Our simulation results show that the addition of stress-relief channels reduces the lateral stress in the Dolan bridge by more than 70% for all the bridge geometries investigated. In practice, our novel mask design significantly increased the survivability of the bridge during device processing, resulting in 100% yield for over 100 Josephson junctions fabricated.

Under IER-441, critical experiments were done with and without tantalum test rods within a central test region surrounded by 7uPCX fuel rods. The experiments were done in new critical assembly hardware designed to support the 7uPCX fuel in a 1.02 cm triangular-pitched array

Agent-based models (ABM) provide an excellent framework for modeling outbreaks and interventions in epidemiology by explicitly accounting for diverse individual interactions and environments. However, these models are usually stochastic and highly parametrized, requiring precise calibration for predictive performance. When considering realistic numbers of agents and properly accounting for stochasticity, this high-dimensional calibration can be computationally prohibitive. This paper presents a random forest-based surrogate modeling technique to accelerate the evaluation of ABMs and demonstrates its use to calibrate an epidemiological ABM named CityCOVID via Markov chain Monte Carlo (MCMC). The technique is first outlined in the context of CityCOVID's quantities of interest, namely hospitalizations and deaths, by exploring dimensionality reduction via temporal decomposition with principal component analysis (PCA) and via sensitivity analysis. The calibration problem is then presented, and samples are generated to best match COVID-19 hospitalization and death numbers in Chicago from March to June in 2020. These results are compared with previous approximate Bayesian calibration (IMABC) results, and their predictive performance is analyzed, showing improved performance with a reduction in computation.

Qubit performance is often reported in terms of a variety of single-value metrics, each providing a facet of the underlying noise mechanism limiting performance. However, the value of these metrics may drift over long timescales, and reporting a single number for qubit performance fails to account for the low-frequency noise processes that give rise to this drift. In this work, we demonstrate how we can use the distribution of these values to validate or invalidate candidate noise models. We focus on the case of randomized benchmarking (RB), where typically a single error rate is reported but this error rate can drift over time when multiple passes of RB are performed. We show that using a statistical test as simple as the Kolmogorov-Smirnov statistic on the distribution of RB error rates can be used to rule out noise models, assuming the experiment is performed over a long enough time interval to capture relevant low frequency noise. With confidence in a noise model, we show how care must be exercised when performing error attribution using the distribution of drifting RB error rate.

This paper develops a novel method for reconstructing the full-field response of structural dynamic systems using sparse measurements. The singular value decomposition is applied to a frequency response matrix relating the structural response to physical loads, base motion, or modal loads. The left singular vectors form a non-physical reduced basis that can be used for response reconstruction with far fewer sensors than existing methods. The contributions of the singular vectors to measured response are termed singular-vector loads (SVLs) and are used in a regularized Bayesian framework to generate full-field response estimates and confidence intervals. The reconstruction framework is applicable to the estimation of single data records and power spectral densities from multiple records. Reconstruction is successfully performed in configurations where the number of SVLs to identify is less than, equal to, and greater than the number of sensors used for reconstruction. In a simulation featuring a seismically excited shear structure, SVL reconstruction significantly outperforms modal FRF-based reconstruction and successfully estimates full-field responses with as few as two uniaxial accelerometers. SVL reconstruction is further verified in a simulation featuring an acoustically excited cylinder. Finally, response reconstruction and uncertainty quantification are performed on an experimental structure with three shaker inputs and 27 triaxial accelerometer outputs.

Hydrogen geo-storage is attracting substantial interdisciplinary interest as a cost-effective and sustainable option for medium- and long-term storage. Hydrogen can be stored underground in diverse formations, including aquifers, salt caverns, and depleted oil and gas reservoirs. The wetting dynamics of the hydrogen-brine-rock system are critical for assessing both structural and residual storage capacities, and ensuring containment safety. Through molecular dynamics simulations, we explore how varying concentrations of cushion gases (CO2 or CH4) influence the wetting properties of hydrogen-brine-clay systems under geological conditions (15 MPa and 333 K). We employed models of talc and the hydroxylated basal face of kaolinite (kaoOH) as clay substrates. Our findings reveal that the effect of cushion gases on hydrogen-brine-clay wettability is strongly dependent on the clay-brine interactions. Notably, CO2 and CH4 reduce the water wettability of talc in hydrogen-brine-talc systems, while exerting no influence on the wettability of hydrogen-brine-kaoOH systems. Detailed analysis of free energy of cavity formation near clay surfaces, clay-brine interfacial tensions, and the Willard-Chandler surface for gas-brine interfaces elucidate the molecular mechanisms underlying wettability changes. Our simulations identify empirical correlations between wetting properties and the average free energy required to perturb a flat interface when clay-brine interactions are less dominant. Our thorough thermodynamic analysis of rock-fluid and fluid-fluid interactions, aligning with key experimental observations, underscores the utility of simulated interfacial properties in refining contact angle measurements and predicting experimentally relevant properties. These insights significantly enhance the assessment of gas geo-storage potential. Prospectively, the approaches and findings obtained from this study could form a basis for more advanced multiscale simulations that consider a range of geological and operational variables, potentially guiding the development and improvement of geo-storage systems in general, with a particular focus on hydrogen storage.

The presence of defects can significantly improve catalytic activity and stability, as they influence the binding of the reactants, intermediates, and products to the catalyst. Controlling defects in the structures of nanocrystal catalysts is synthetically challenging. In this study, we demonstrate the ability to control the growth of Ir nanocrystals, enabling the tuning of both structural and surface defects. The Ir nanocrystals have unique structures that range from single crystals of a few nanometers to twinned nanoparticles and multiply twinned crystallites with a high density of atomic steps. Further, this approach of defect engineering enables us to understand their roles in enhancing the performance of the OER and producing an Ir catalyst with both high activity and stability. Our results show the importance of the concept of using synthetic control of structural and surface defects in metal nanoparticles as a strategy to improve catalytic performance.

Vat photopolymerization (VP) has emerged as a promising additive manufacturing technique to allow rapid light-based fabrication of 3D objects from a liquid resin. Research in the field of vat photopolymerization spans across multiple disciplines from engineering and materials science to applied chemistry and physics. This perspective brings together early-career researchers from various disciplines in academia and national laboratories around the world to summarize the most recent advancements with special emphasis on the research highlighted as part of the Gordon Research Conference (GRC) 2024 meeting on Additive Manufacturing of Soft Materials. We provide an outlook on next-generation polymer processing methods from synthesis of novel materials to multimodality manufacturing and performance engineering. Further, this article combines the ideas of many of these junior researchers to present a vision for the future of the field by highlighting the challenges and opportunities that lie ahead.

Nitric oxide molecular tagging velocimetry (NO MTV) is used to characterize the wake behind a free-flight spherical model in the test section of a free-piston reflected shock tunnel using a burst-mode laser operated at 100 kHz. A novel (to our knowledge) multi-delay timing scheme was implemented to measure velocity in a varying collisional environment in the wake of the free-flight model. Four simultaneous velocity profiles were measured in the wake of the model from -600 to 3600 m/s for flow enthalpies of 10.3 and 12.0 MJ/kg. Finally, the measured velocity distributions show good agreement when compared to computational fluid dynamics (CFD) modeling.

A unified model for the work hardening Θσ_Θ(σ) and stress-strain σε_σ(ε) behavior is presented that accounts for deformation under tensile loading, from the onset of yielding at the proportional limit up to the ultimate strength as defined at the structural instability. The origin of this approach is based on a negative exponential formulation for an asymptotic-curvilinear work-hardening model that accounts for the rapid strengthening of metals as well as the continuation of steady-state deformation to the instability.

In high-purity polycrystalline metallic materials, voids tend to favor grain boundaries as nucleation sites due to the elevated stress states produced by granular interactions and the weakened grain boundary from the relative atomic disorder. To quantify the key factors of this elevated stress state, simple compression of a small multi-grain cylinder of body-centered cubic tantalum was simulated using a single crystal plasticity model that incorporates non-Schmid effects. Four increasingly complex synthetic microstructures were created to tractably incorporate grain boundary interactions, and a statistically significant number of combinations were performed by varying the initial crystallographic orientations of the microstructure. Most of these simulations produce the maximum von Mises stress on a grain boundary and less frequently at the multi-grain junctions. To build a statistical model for the maximum von Mises stress at the grain boundary, physically based features that could contribute to the elevated stress state were selected. Then, a learning algorithm based on information theory was used to identify which of these features contributed the most information to the data set. The identified features include a grain's propensity to accommodate both elastic and plastic deformations and their directional components. The misalignment of the direction of each grain's mechanical response was found to be strongly correlated to the magnitude of the stress near the grain boundary. For all of the synthetic microstructures, the statistical models produce a residual distribution that is nearly Gaussian with a variance of, at most, 10% of the prior distribution. The successful performance of the statistical model implies the correct identification of the physical features that cause severe stress localization in polycrystalline materials. The statistical models constructed here can be used to formulate a physically motivated void nucleation model which is sensitive to a microstructure's propensity to produce elevated stress states. These statistical models also enable the design of material microstructures, in which the crystallographic orientation is chosen to resist void nucleation.

Significant events, such as volcanic eruptions, can have global and long-lasting impacts on climate. These global impacts, however, are not uniform across space and time. Understanding how the Mt. Pinatubo eruption affects global and regional climate is of great interest for predicting the impact on climate due to similar events as well as understanding the possible effect of the stratospheric aerosol injections proposed to combat climate change. While many studies illustrated the impact of the Pinatubo eruption on a global scale, studies at a fine regional scale are scarce. We propose a Bayesian spatially-varying changepoint detection and estimation method to trace the impact of Mt. Pinatubo eruption on regional climate. Our approach takes into account the diffusing nature and spatial correlation of the climate changes attributed to the volcanic eruption. We illustrate our method and demonstrate its advantages over an existing changepoint detection method through simulations. Finally, we apply our method to monthly stratospheric aerosol optical depth and surface temperature data from 1985 to 1995 to detect and estimate change-points following the 1991 Mt. Pinatubo eruption. Our results quantitatively characterize the spatial pattern of the eruption’s impact on regional climate, complementing the previous studies on the global impact of the Pinatubo eruption.

The Reference Model (RM) project developed six marine energy converter concepts using a sequential design methodology, which, while widely adopted in the industry, often overlooks interactions between system components, resulting in suboptimal designs. One such example is the Reference Model 3 (RM3), a two-body point absorber wave energy converter (WEC). An assessment using the Technology Performance Level (TPL) revealed that RM3’s low power-to-cost ratio, partly due to expensive steel construction, limits its techno-economic performance. This study aims to redesign RM3 by reducing its scale and employing control co-design to integrate WEC and Power Take-Off (PTO) dynamics, constraints, and cost considerations within an optimization framework. We demonstrate the limitations of RM3’s current PTO design and explore the benefits of scaling down to enhance techno-economic viability by lowering material costs. Using WecOptTool, we conduct a parameter sweep over gear ratios and spring stiffnesses for various Commercial Off-The-Shelf generators in irregular wave conditions. Our findings emphasize the importance of aligning PTO components with WEC dynamics, showing that control co-design and strategic scaling can improve RM3’s power-to-cost ratio. This study presents a transferable example of applied control co-design for other WECs, supporting early-stage developers in their design decisions.

An increasing magnetic field perpendicular to an undoped semiconductor surface at low temperature is known to strengthen the binding of localized electrons to stationary ions, as the wavefunction's tails evolve from exponential to Gaussian. It is also known that application of a high bias voltage to a depleted semiconductor can liberate bound charge and induce a large drop in electrical resistance. We connect these established results to experimental electrical transport measurements on off-state germanium Schottky-barrier metal-oxide-semiconductor field-effect transistor (MOSFETs) with an aluminum oxide insulating dielectric and platinum germanide contacts. We make measurements at the three distinct orientations of the magnetic field with respect to the substrate and the current. At 6 K, we observe sharp attenuation of current by more than 2 orders of magnitude, within 60 mT, at a crossover magnetic field perpendicular to the substrate. A 1 T magnetic field attenuates the current by more than 4 orders of magnitude. The strength of the attenuation and the value of the crossover field are controlled by both the gate-source and drain-source voltages. The attenuation is much weaker when the magnetic field is parallel to the current. Finally, we orient the magnetic field parallel to the substrate, but perpendicular to the current, allowing us to distinguish charge hopping at the oxide interface from charge hopping in the bulk. This large off-state magnetoresistance can be exploited for cryogenic magnetic- and photo-detection, and for high-bias, low-leakage MOSFETs.

The computational cost of high-fidelity numerical models makes outer-loop analysis, which requires repeated interrogation of the model such as uncertainty quantification, computationally demanding. Multi-fidelity methods, which construct a surrogate model using data from an ensemble of models of varying cost and accuracy, can substantially reduce the cost of outer-loop analysis. However, these methods can be difficult to apply when the model ensemble does not admit a clear hierarchy a priori and the correlations between models are low. Consequently, in this paper, we present a multi-fidelity method that leverages dimension reduction to enhance the correlation between models, thereby reducing the amount of data needed to train a surrogate from an unordered ensemble of models. Our method utilizes basis adaptation to build low-dimensional polynomial chaos expansions of each model and employs Multi-fidelity Networks to encode the relationships among models. We show that the resulting method exhibit two notable advantages over its counterpart: (1) enhanced accuracy (both reduced bias and variance); and (2) reduced dependency on the graph structure encoding relationships among models. We demonstrate the approach on an analytical test problem and a challenging finite element model for a spent nuclear fuel. Our method produces a surrogate model that is significantly more accurate than either a single-fidelity surrogate or a multi-fidelity surrogate constructed without basis adaptation.

We consider numerical approaches for deterministic, finite-dimensional optimal control problems whose dynamics depend on unknown or uncertain parameters. We seek to amortize the solution over a set of relevant parameters in an offline stage to enable rapid decision-making and be able to react to changes in the parameter in the online stage. To tackle the curse of dimensionality arising when the state and/or parameter are highdimensional, we represent the policy using neural networks. We compare two training paradigms: First, our model-based approach leverages the dynamics and definition of the objective function to learn the value function of the parameterized optimal control problem and obtain the policy using a feedback form. Second, we use actor-critic reinforcement learning to approximate the policy in a data-driven way. Using an example involving a two-dimensional convection-diffusion equation, which features high-dimensional state and parameter spaces, we investigate the accuracy and efficiency of both training paradigms. While both paradigms lead to a reasonable approximation of the policy, the model-based approach is more accurate and considerably reduces the number of PDE solves.

Hydrogen is a key pillar in the global Net Zero strategy. Rapid scaling up of hydrogen production, transport, distribution and utilization is expected. This entails that hydrogen, which is traditionally an industrial gas, will come into proximity of populated urban areas and in some situations handled by the untrained public. To realize all their benefits, hydrogen and its technologies must be safely developed and deployed. The specific properties of hydrogen involving wide flammability range, low ignition energy and fast flame speed implies that any accidental release of hydrogen can be easily ignited. Comparing with conventional fuels, combustion systems fueled by hydrogen are also more prone to flame instability and abnormal combustion. This paper aims to provide a comprehensive review about combustion research related to hydrogen safety. It starts with a brief introduction which includes some overview about risk analysis, codes and standards. The core content covers ignition, fire, explosions and deflagration to detonation transition (DDT). Considering that DDT leads to detonation, and that detonation may also be induced directly under special circumstances, the subject of detonation is also included for completeness. The review covers laboratory, medium and large-scale experiments, as well as theoretical analysis and numerical simulation results. While highlights are provided at the end of each section, the paper closes with some concluding remarks highlighting the achievements and key knowledge gaps.

The magnetostrictive response of a Terfenol-D pellet was measured via a laboratory-based X-ray diffractometer. X-ray diffraction patterns were collected from the pellet sample with and without the presence of an applied magnetic field (~30 mT) generated by placing a large magnet under the pellet. A standard reference material, Silicon 640c, was employed as an internal standard. Magnetostriction values of 323 and 227 ppm Δl/l were determined for the (104) and (110) indexed peaks, respectively, assuming a rhombohedral structure for Terfenol-D. A threshold noise level value of ~20 to 30 ppm Δl/l was suggested based on before/after measurements in the absence of the applied field. No clear evidence of domain wall rotation was detected via changes in relative intensities of diffraction peaks in the presence of the applied magnetic field.

Development of compact photocurrent models is currently dominated by analytical techniques that rely on physical assumptions to render the governing equations solvable in a closed form. Violation of these assumptions can reduce the accuracy of the models and/or limit their scope. In this paper we show that system identification of nonlinear state-space systems can serve as an alternative numerical basis for non-intrusive reduced order modeling of photocurrent effects. To that end we develop a compact gray box photocurrent model (GBPM) by using a state-space representation with a low-dimensional latent state equation that mimics a mathematical model for the response of an idealized class of devices to ionizing radiation. In so doing we obtain a model that learns the dynamics of a quantity of interest directly from its measurements without requiring snapshots of the internal device state or its discretized model, and can be inferred from very small data sets. To demonstrate the approach we train the GBPM using a small experimental data set for a Z5236 Zener diode and a small synthetic data set obtained by simulating a synthetic pn-junction device. We then compare the GBPMs with black box models trained on the same data and show that performance of the latter is limited by the size of the data set, while the former are able to achieve excellent performance in both the reproductive and the predictive regimes.

Chirality plays a crucial role in the helicity mismatch between surface acoustic waves and magnetic spin waves, leading to nonreciprocal transmission of acoustic power for coupled magnetoacoustic modes. Acoustic modes with both longitudinal and shear strain exhibit elliptical particle displacements, making them chiral, and different acoustic modes can exhibit different helicities of this elliptical particle displacement. Here, we study chiral acoustic modes with different helicities supported on the same piezoelectric platform and their interaction with magnetic spin waves. Our study demonstrates that the nonreciprocal transmission of acoustic power is driven by the helicity mismatch effect and, specifically, that the handedness of the nonreciprocity is based on whether the surface acoustic wave has retrograde (Rayleigh mode) or prograde (Sezawa mode) elliptical particle displacement with respect to the propagation direction. We found the transmission nonreciprocity to be significant, with 7.3 dB/mm for the retrograde particle displacement (Rayleigh mode at 2.358 GHz) and 3.3 dB/mm for prograde particle displacement (Sezawa mode at 3.112 GHz). This work highlights that piezoelectric platforms can be engineered to support acoustic modes with opposite helicities to enable frequency-selective nonreciprocal radiofrequency and microwave components, such as isolators and circulators, through coupled acoustic spin wave interactions.

High-dimensional quantum systems are a valuable resource for quantum information processing. They can be used to encode error-correctable logical qubits, which has been demonstrated using continuous-variable states in microwave cavities or the motional modes of trapped ions. For example, high-dimensional systems can be used to realize ‘Schrödinger cat’ states, which are superpositions of widely displaced coherent states that can be used to illustrate quantum effects at large scales. Recent proposals have suggested encoding qubits in high-spin atomic nuclei, which are finite-dimensional systems that can host hardware-efficient versions of continuous-variable codes. Here we demonstrate the creation and manipulation of Schrödinger cat states using the spin-7/2 nucleus of an antimony atom embedded in a silicon nanoelectronic device. We use a multi-frequency control scheme to produce spin rotations that preserve the symmetry of the qudit, and we constitute logical Pauli operations for qubits encoded in the Schrödinger cat states. Our work demonstrates the ability to prepare and control non-classical resource states, which is a prerequisite for applications in quantum information processing and quantum error correction, using our scalable, manufacturable semiconductor platform.

Projection-based reduced order models are effective at approximating parameter-dependent differential equations that are parametrically separable. When parametric separability is not satisfied, which occurs in both linear and nonlinear problems, projection-based methods fail to adequately reduce the computational complexity. Devising alternative reduced order models is crucial for obtaining efficient and accurate approximations to expensive high-fidelity models. In this work, we develop a timestepping procedure for dynamical parameter-dependent problems, in which a neural-network is trained to propagate the coefficients of a reduced basis expansion. This results in an online stage with a computational cost independent of the size of the underlying problem. We demonstrate our method on several parabolic partial differential equations, including a problem that is not parametrically separable.

We introduce physics-informed multimodal autoencoders (PIMA)-a variational inference framework for discovering shared information in multimodal datasets. Individual modalities are embedded into a shared latent space and fused through a product-of-experts formulation, enabling a Gaussian mixture prior to identify shared features. Sampling from clusters allows cross-modal generative modeling, with a mixture-of-experts decoder that imposes inductive biases from prior scientific knowledge and thereby imparts structured disentanglement of the latent space. This approach enables cross-modal inference and the discovery of features in high-dimensional heterogeneous datasets. Consequently, this approach provides a means to discover fingerprints in multimodal scientific datasets and to avoid traditional bottlenecks related to high-fidelity measurement and characterization of scientific datasets.

Abstract not provided.

This document evaluates a hypothetical nuclear powerplant and associated protective force personnel using modern modeling and simulation tools. The facility incorporates security early in the design to consider and integrate methods to resolve security issues and vulnerabilities via the facility’s inherent design characteristics before construction. The evaluation in this document is an example only. It is not intended to recommend or evaluate the effectiveness of existing physical security requirements or identify any method that the U.S. Nuclear Regulatory Commission staff may find acceptable for complying with existing requirements.

Reliably simulating experiments relevant to the National Nuclear Security Administration (NNSA) requires a detailed description of material properties across a wide range of conditions. Such properties include the equations of state, charged-particle transport coefficients, and optical properties like the opacity. Together, these properties make up the material models used in radiation-magnetohydrodynamic simulations of nuclear fusion experiments. Many of these models do not incorporate uncertainties in the data used to produce them. It is unknown whether these uncertainties significantly impact the interpretation of simulation results and diagnostics. The purpose of this work is to quantify how such uncertainties impact simulations of pulsed-power experiments. We accomplished this task by first assessing discrepancies between approaches used to generate the data. This included bringing together members of the high-energy-density community spanning the three NNSA laboratories and multiple universities. Then, using these data, we developed a general framework that systematically incorporates physical uncertainties within the material models suitable for uncertainty quantification analyses. The framework utilizes machine learning, Bayesian inference, and incorporates multi-fidelity datasets. We demonstrated the framework by quantifying the impact that material model uncertainties have on simulations of pulsed-power experiments underway on Z at Sandia National Laboratories. As a result of this work, we discovered that modest uncertainties in material models (roughly 20%) correspond to significant uncertainties in the outputs from simulations. Our framework has enabled rapid construction of material models through an automated procedure and allows for the generation of material models of interest to the NNSA.

Abstract not provided.

The objective of this work was to develop a machine learning ensemble that could assist pebble bed reactor verification by evaluating whether a given pebble circulating through a PBR was normal or anomalous using gamma spectroscopy measurements from a notional PBR burnup measurement system. Using a PBR reference design, data sets of synthetic gamma spectra representative of BUMS measurements of normal and anomalous pebbles that may be used to produce special fissile material were generated to train and test an ML anomaly detection ensemble on two reference scenarios – substitution of normal pebbles with target pebbles for production of Pu or ²³³U. The ML ensemble correctly identified all anomalous pebbles in the testing data set, and while perfect ensemble performance is normally indicative of overfitting, it was concluded that significantly lower photon intensity of target pebbles produced distinctly less intense photon spectra to where perfect ensemble performance was expected.

Abstract not provided.

Dr. Judi E. See, a Systems Analyst and Human Factors Engineer at Sandia National Laboratories, reflects on her experience conducting NDT research in a male-dominated environment. She emphasizes the importance of persistence, flexibility, and persuasive skills in overcoming challenges, ranging from gaining access to test sites and equipment to building trust with inspectors. She shares her personal experience of navigating professional situations where gender disparities were evident, highlighting the need for women to adapt and overcome obstacles in traditionally male-dominated settings. See's journey demonstrates that perseverance and ingenuity can lead to significant contributions, process improvements, and recognition in the NDT field.

Based on planning and reviews conducted in 2020, it is clear today that the Z Facility must be sustained through 2040 or longer. This requires a well-developed and executed maintenance and sustainment effort. Z’s Sustainment Plan outlines the approach that Sandia’s Pulsed Power Sciences Center will pursue to assure that the Z Facility maintains its world class research and development program supporting the NNSA’s Science Based Stockpile Stewardship efforts. In FY24, the Pulsed Power Center at SNL planned and executed many substantial sustainment projects as well as maintenance-driven-sustainment aligned to the Sustainment Plan. These projects are listed in this section with brief descriptions of the work executed.

Abstract not provided.

The objective of this work was to develop a machine learning ensemble that could assist pebble bed reactor verification by evaluating whether a given pebble circulating through a PBR was normal or anomalous using gamma spectroscopy measurements from a notional PBR burnup measurement system. Using a PBR reference design, data sets of synthetic gamma spectra representative of BUMS measurements of normal and anomalous pebbles that may be used to produce special fissile material were generated to train and test an ML anomaly detection ensemble on two reference scenarios – substitution of normal pebbles with target pebbles for production of Pu or ²³³U. The ML ensemble correctly identified all anomalous pebbles in the testing data set, and while perfect ensemble performance is normally indicative of overfitting, it was concluded that significantly lower photon intensity of target pebbles produced distinctly less intense photon spectra to where perfect ensemble performance was expected.

Abstract not provided.

The purpose of this protocol is to define procedures and practices to be used by the PACT center for field testing of metal halide perovskite (MHP) photovoltaic (PV) modules. The protocol defines the physical, electrical, and analytical configuration of the tests and applies equally to mounting systems at a fixed orientation or sun tracking systems. While standards exist for outdoor testing of conventional PV modules, these do not anticipate the unique electrical behavior of perovskite cells. Further, the existing standards are oriented toward mature, relatively stable products with lifetimes that can be measured on the scale of years to decades. The state of the art for MHP modules is still immature with considerable sample to sample variation among nominally identical modules. Version 0.0 of this protocol does not define a minimum test duration, although the intent is for modules to be fielded for periods ranging for weeks to months. This protocol draws from relevant parts of existing standards, and where necessary includes modifications specific to the behavior of perovskites.

Abstract not provided.

Sandia is a federally funded research and development center (FFRDC) focused on developing and applying advanced science and engineering capabilities to mitigate national security threats. This is accomplished through the exceptional staff leading research at the Labs and partnering with universities and companies. Sandia’s LDRD program aims to maintain the scientific and technical vitality of the Labs and to enhance the Labs’ ability to address future national security needs. The program funds foundational, leading-edge discretionary research projects that cultivate and utilize core science, technology, and engineering (ST&E) capabilities. Per Congressional intent (P.L. 101-510) and Department of Energy (DOE) guidance (DOE Order 413.2C, Chg 1), Sandia’s LDRD program is crucial to maintaining the nation’s scientific and technical vitality.

One of the objectives of the United States Department of Energy Office of Nuclear Energy’s Office of Spent Fuel and High-Level Waste Disposition is to better understand the technical bases, risks, and uncertainties associated with the safe and secure disposition of spent nuclear fuel and high-level radioactive waste. Domestic defense and research activities have generated a few thousand metric tons of spent nuclear fuel and high-level radioactive waste, much of which has been or will be processed and vitrified into high-level waste glass. The Nuclear Waste Policy Act 1982 makes the Department of Energy responsible for disposal of these materials.

Disinformation poses a significant and evolving threat to today’s online environment. Individuals encounter challenges in detecting disinformation, subsequently influencing their behavior and decision-making processes. Our research examines the potential use of large language model (LLM) generative agents (LGAs) to replicate human behavior to better understand how disinformation is spread in online environments. Using human subjects research, we first investigate how personality traits, individual differences, and demographic factors relate to decision-making in simulated online disinformation environments. Then, we examine whether LGAs can effectively replicate human responses in the same simulated online environments when assigned personality traits, demographic characteristics and behavioral attributes. Our findings indicate that LGAs can align with human decisions in these scenarios; however, alignment is contingent upon scenario context, persona settings and LLM selection. Results provide valuable insights for methodology refinement in future research and in utilizing LGAs to model complex national security challenges such as disinformation campaigns.

Synthetic aperture radar (SAR) images formed with dechirp-on-receive data collection and rectangular format processing algorithm are the result of a two-dimensional discrete Fourier transform (DFT) applied to sampled data. There are several steps required to compute the range and Doppler values associated with each pixel in a SAR range-Doppler image. This memo walks readers through the process.

Salt formations have long been recognized as a highly favorable host rock for the final disposal of high-level radioactive waste (HLW) in deep geological repositories. Their unique properties, including exceptional impermeability, self-healing capabilities, and thermal conductivity, make them a reliable natural barrier for the deep disposal of radioactive waste. This report focuses on the development and application of a methodology for assessing the integrity and per formance of the Engineered Barrier System (EBS) within salt-based repositories, a critical component of the multi-barrier system ensuring safe radioactive waste disposal.

Abstract not provided.

This report provides a comprehensive assessment of physical security modeling and simulation tools available for use in the vulnerability assessment (VA) process for nuclear facilities. It outlines the historical evolution of VA methodologies, emphasizing the transition from traditional layer-based approaches to a more holistic framework that integrates detection probabilities directly into combat simulations. The document details the critical components of the VA process, including the characterization of targets, threats, and protective measures, as well as the development of adversary scenarios that reflect both insider and outsider threats. It highlights the importance of performance assurance programs, emphasizing the need for continuous evaluation and testing of security systems to ensure their effectiveness against evolving threats. Additionally, the report discusses the significance of utilizing accredited modeling and simulation tools in accredited areas to accurately represent adversary actions and the corresponding responses of protective forces. By establishing a systematic approach to VA, this document aims to enhance the overall security posture of nuclear facilities, ensuring compliance with regulatory standards while effectively mitigating risks associated with potential adversarial actions.

Under direction from the DOE Office of Electricity, Sandia National Laboratories performed testing of the Viziv system to evaluate the quality of the Zenneck surface wave and potential application to long range power transfer. This report documents the test methodology as well as the test results. This includes an analysis of prior test data collected by Viziv.

This report documents a cross-site effort to identify and evaluate materials that are suitable replacements for the polyurethanes known as EN-7 and EN-8. EN-7 and EN-8 are commercial formulations that contain free toluene diisocyanate (TDI) and which are used as adhesives and as encapsulants in many weapon applications. TDI is an OSHA-regulated volatile diisocyanate and has been targeted for elimination from future weapons use, prompting the need for a TDI-free replacement material.

The purpose of this protocol is to define procedures and practices to be used by the PACT center for field testing of metal halide perovskite (MHP) photovoltaic (PV) modules. The protocol defines the physical, electrical, and analytical configuration of the tests and applies equally to mounting systems at a fixed orientation or sun tracking systems. While standards exist for outdoor testing of conventional PV modules, these do not anticipate the unique electrical behavior of perovskite cells. Further, the existing standards are oriented toward mature, relatively stable products with lifetimes that can be measured on the scale of years to decades. The state of the art for MHP modules is still immature with considerable sample to sample variation among nominally identical modules. Version 0.0 of this protocol does not define a minimum test duration, although the intent is for modules to be fielded for periods ranging for weeks to months. This protocol draws from relevant parts of existing standards, and where necessary includes modifications specific to the behavior of perovskites.

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. More importantly Arpeggio orchestrates the execution of applications that participate in the coupling. This document describes the various components of Arpeggio and their operability. The intent of the document is to provide a fast path for analysts interested in coupled applications via simple examples of its usage.

This report examines narrow slot maximum versus average transmission cross sections and the resulting directivity gain. Other aperture shapes at high frequencies are also briefly discussed.

Non-von Neumann computing using neuromorphic systems based on analogue synaptic and neuronal elements has emerged as a potential solution to tackle the growing need for more efficient data processing, but progress toward practical systems has been stymied due to a lack of materials and devices with the appropriate attributes. Recently, solid state electrochemical ion-insertion, also known as electrochemical random access memory (ECRAM) has emerged as a promising approach to realize the needed device characteristics. ECRAM is a three terminal device that operates by tuning electronic conductance in functional materials through solid-state electrochemical redox reactions. This mechanism can be considered as a gate-controlled bulk modulation of dopants and/or phases in the channel. Early work demonstrating that ECRAM can achieve nearly ideal analogue synaptic characteristics has sparked tremendous interest in this approach. More recently, the realization that electrochemical ion insertion can be used to tune the electronic properties of many types of materials including transition metal oxides, layered two-dimensional materials, organic and coordination polymers, and that the changes in conductance can span orders of magnitude has further attracted interest in ECRAM as the basis for analogue synaptic elements for inference accelerators as well as for dynamical devices that can emulate a wide range of neuronal characteristics for implementation in analogue spiking neural networks. At its core, ECRAM shares many fundamental aspects with rechargeable batteries, where ion insertion materials are used extensively for their ability to reversibly store charge and energy. Computing applications, however, present drastically different requirements: systems will require many millions of devices, scaled down to tens of nanometers, all while achieving reliable electronic-state tuning at scaled-up rates and endurances, and with minimal energy dissipation and noise. In this review, we discuss the history, basic concepts, recent progress, as well as the challenges and opportunities for different types of ECRAM, broadly grouped by their primary mobile ionic charge carrier, including Li, protons, and oxygen vacancies.

A 2D phononic pseudocrystal isolator exhibiting cyclic symmetry and radial self-similarity is measured and demonstrated to block a wide range of ultrasonic vibration. Measurements of longitudinal and shear wave blocking effects are made and compared with computational results. The use of the bandgap edge ratio is recommended for quantifying suppression in very-wide-bandgap materials. The upper-to-lower suppression edge frequency ratios of 3-4 are remarkably large for shear waves and even larger for longitudinal waves upper-to-lower suppression ratio (13 at 5 dB), such that 92.5% of frequencies in that range experience ≥ 5 dB of suppression.

Permafrost thaw increases the bioavailability of ancient organic matter, facilitating microbial metabolism of volatile organic compounds (VOCs), carbon dioxide, and methane (CH₄). The formation of thermokarst (thaw) lakes in icy, organic-rich Yedoma permafrost leads to high CH₄ emissions, and subsurface microbes that have the potential to be biogeochemical drivers of organic carbon turnover in these systems. However, to better characterize and quantify rates of permafrost changes, methods that further clarify the relationship between subsurface biogeochemical processes and microbial dynamics are needed. In this study, we investigated four sites (two well-drained thermokarst mounds, a drained thermokarst lake, and the terrestrial margin of a recently formed thermokarst lake) to determine whether biogenic VOCs (1) can be effectively collected during winter, and (2) whether winter sampling provides more biologically significant VOCs correlated with subsurface microbial metabolic potential. During the cold season (March 2023), we drilled boreholes at the four sites and collected cores to simultaneously characterize microbial populations and captured VOCs. VOC analysis of these sites revealed “fingerprints” that were distinct and unique to each site. Total VOCs from the boreholes included > 400 unique VOC features, including > 40 potentially biogenic VOCs related to microbial metabolism. Subsurface microbial community composition was distinct across sites; for example, methanogenic archaea were far more abundant at the thermokarst site characterized by high annual CH₄ emissions. The results obtained from this method strongly suggest that ∼10% of VOCs are potentially biogenic, and that biogenic VOCs can be mapped to subsurface microbial metabolisms. By better revealing the relationship between subsurface biogeochemical processes and microbial dynamics, this work advances our ability to monitor and predict subsurface carbon turnover in Arctic soils.

BiVO4 photoanodes are promising for solar water splitting, with photogenerated electrons and holes preferentially reacting at top {010} and lateral {110} facets, respectively. However, the mechanisms driving this facet-dependent reactivity remain unclear. Here, we investigate facet-dependent photocurrent and material heterogeneity using correlative scanning photoelectrochemical microscopy (SPCM), electron beam induced current (EBIC) mapping, and mid-IR scattering scanning near-field optical microscopy (s-SNOM). SPCM measurements of 62 BiVO4 particles confirmed higher photocurrents at lateral {110} facets compared to top {010} facets, but unexpectedly revealed variations in photocurrent among lateral facets within the same particle. Variations in lateral facet surface termination could explain the intraparticle-level reactivity heterogeneity, consistent with theoretical predictions. Nano-FTIR spectroscopy and Raman microspectroscopy indicated significant materials chemistry heterogeneity within individual particles and facets that could be attributed to variations in lattice vibration distortions that enhance the overlap between Bi 6s and O 2p orbitals. The increased orbital overlap is significant as it potentially increases hole mobility in the valence band and potentially explains the lateral facet-dependent charge separation efficiency observed in photocurrent maps. Facet-dependent electrical and EBIC measurements showed no space charge regions at interfacet junctions or metal-BiVO4 contacts under vacuum, suggesting that photogenerated holes beneath top {010} facets are unlikely to transport to lateral {110} facets to drive water/sulfite oxidation. These findings indicate the potential influence of distinct bulk properties and surface termination chemistries across different particles and facets, highlighting the importance of carefully controlling defects and surface chemistry during sample growth to optimize photocatalytic performance.

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni _xAl _y⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni _xAl _y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

Production of low carbon gasoline-like fuels such as methanol-to-gasoline (MTG) is a promising approach to achieve rapid greenhouse gas emission reduction of the transportation sector. Despite the fact that gasoline that meets the ASTM D4814 standard for automotive spark-ignition engine fuel can be readily produced from these processes, it is unclear how the composition of MTG may affect engine performance and emissions. In this paper, a surrogate for an MTG is used to numerically study the effects of gasoline composition on knock propensity and on the sensitivity of knock to thermal and fuel stratification, to oxygen dilution and to nitric oxide from exhaust gas recirculation of residual gases. Simulations were performed in ANSYS CHEMKIN-PRO using a comprehensive chemical kinetic mechanism for gasoline surrogates, and results of the MTG surrogate were compared against those of a petroleum-based regular E10 gasoline, termed PACE-20. A premium-grade MTG fuel was also formulated by adding ethanol to the MTG surrogate, and results were compared against those of four premium-grade, gasoline-like fuels representative of future alternative gasoline formulations. Surrogates and mechanism were evaluated by comparison against experimental engine data, and the model showed high accuracy at stoichiometric conditions (mean absolute error of ignition timing equal to 1.46 crank angle degrees) but larger deviations at lean conditions (mean absolute error of ignition timing equal to 5.52 crank angle degrees). Despite the fact that the MTG surrogate has a RON 1.1 units higher than that of PACE-20, it may show higher knock propensity at medium temperature conditions due to a less intense NTC behavior. MTG autoignition was more temperature- and equivalence ratio-sensitive than that of PACE20, suggesting that MTG can benefit more from naturally-occurring thermal stratification or from induced fuel stratification of the end gas to mitigate knock intensity. The sensitivity of autoignition reactivity to oxygen dilution and to NO concentration was higher for MTG than for regular gasoline at medium loads, but the opposite trend was observed at high loads due to the effect of pressure on the low-temperature chemistry of regular gasoline. Approximately 14 %vol ethanol content was required to upgrade the octane rating of MTG from regular grade to premium grade. Adding 13.6 %vol ethanol made the fuel autoignition less sensitive to both oxygen dilution and NO content (ignition time varies approx. 17 % and 50 % less with oxygen dilution and NO addition, respectively, when adding ethanol at high engine loads).

An optimization-based design framework is proposed to tune the response of soft metamaterials involving both geometric instabilities and nonlinear viscoelastic material behavior. Designing the response of soft metmaterials to harness instabilities and undergo large, tailored configuration changes will enable advancements in soft robotics, shock and vibration mitigation, and flexible electronics. In line with the metamaterial concept, the response of these materials is governed to a large extent by the geometric and topological makeup of their small-scale features. However, the link between structure and response is less intuitive for soft metamaterials due to their reliance upon highly nonlinear responses triggered by geometric instabilities. This is further complicated by the effects of viscoelastic relaxation, which recent studies have shown to alter the emergence of instabilities in non-intuitive ways. hese effects are accounted for in our framework to achieve various design objectives, including tailored force–displacement response and maximized energy absorption from both geometric and material effects. To fully automate this process, it is essential to have a completely robust equation solver for forward problems involving instabilities and viscoelastic relaxation. We achieve this by casting the search for stable mechanical equilibrium — i.e. the forward problem — as a minimization problem and utilize a trust region algorithm to robustly handle instabilities and follow energetically-favorable equilibrium paths through critical points.

This study investigates the influence of tunneling contact resistances between carbon nanotubes (CNTs) on electron transport and electrical conductivity of macroscopic carbon nanofibers (CNFs), which profoundly impacts the performance of CNT thin film electronics, CNF electron emitters and cathodes, and energy conversion and storage devices. Utilizing a self-consistent electrical contact model coupling a transmission line model with tunneling current, we calculate the contact resistances of a plethora of CNT-CNT contacts within a CNF fiber, which consists of aligned, densely packed CNTs. A statistical analysis is conducted, using Gaussian distributions to account for variations in contact lengths, tunneling gap distances, and single CNT aspect ratios, to calculate the CNT-CNT contact resistance and the overall resistance of CNT fiber. By scaling our model to a macroscopic level, our results are in good agreement with experimental measurements. Our calculation suggests that while increasing the contact overlap length diminishes individual CNT-CNT contact resistance, it could paradoxically increase macroscopic CNT fiber resistance for a given constant CNF mass density, which is due to that fact that a larger overlap length allows more CNTs to pack along an electrical conduction path per unit length, leading to more tunneling contact junctions connected in series and thus less number of parallel conduction paths within the fiber cross section. Increasing tunneling gap distance increases both individual contact and overall fiber resistance. This research provides a simple design tool for tailoring CNT fiber electrical properties to promote real-world applications using CNTs or similar low-dimensional materials.

Underground caverns in a salt dome are promising geologic features to store hydrogen because of salt's extremely low permeability and self-healing behavior. The salt cavern storage community, however, has not fully understood the geomechanical behaviors of salt rock driven by quick operation cycles of injection–production, which may significantly impact the cost-effective storage-recovery performance of multiple caverns. Our field-scale generic model captures the impact of cyclic loading–unloading on the salt creep behavior and deformation under different cycle frequencies, operating pressure, and spatial order of operating cavern(s). This systematic simulation study indicates that the initial operation cycle and arrangement of multiple caverns play a significant role in the creep-driven loss of cavern volumes and cavern deformation. Our future study will develop a new salt constitutive model based on geomechanical tests of site-specific salt rock to probe the cyclic behaviors of salt precisely both beneath and above the dilatancy boundary, including reverse (inverse transient) creep, the Bauschinger effect, and damage-healing mechanism.

At micro- and nanoscales, momentum transfer between surfaces is influenced by various physical mechanisms, including quantum fluctuations, electromagnetic interactions, electric charges, and the dynamics of (rarefied) gases. Under non-isothermal conditions, rarefied gases give rise to thermal Knudsen forces whose magnitudes strongly depend on the gas species and surface characteristics. Knudsen forces are particularly relevant in nanotechnology, optical manipulation, and aerospace systems, where gas rarefaction occurs due to highly confined geometries, sub-micrometer length scales, and reduced particle densities. Despite their significance, predictive modeling of Knudsen forces is limited by a lack of comprehensive experimental data across diverse materials and surface morphologies. In this work, we present a highly sensitive and adaptable measurement platform capable of directly quantifying Knudsen forces using a suspended, interchangeable micro-cantilever within controlled rarefied helium and nitrogen environments. The system integrates optical fiber interferometry to precisely capture out-of-plane displacements at sub-micrometer resolution, driven by Knudsen forces. From the empirical data, we derive a robust correlation linking the magnitudes of Knudsen forces to energy accommodation coefficients, offering deeper insights into the underlying gas-surface interaction mechanisms.

Underground caverns in a salt dome are promising geologic features to store hydrogen because of salt's extremely low permeability and self-healing behavior. The salt cavern storage community, however, has not fully understood the geomechanical behaviors of salt rock driven by quick operation cycles of injection–production, which may significantly impact the cost-effective storage-recovery performance of multiple caverns. Our field-scale generic model captures the impact of cyclic loading–unloading on the salt creep behavior and deformation under different cycle frequencies, operating pressure, and spatial order of operating cavern(s). This systematic simulation study indicates that the initial operation cycle and arrangement of multiple caverns play a significant role in the creep-driven loss of cavern volumes and cavern deformation. Our future study will develop a new salt constitutive model based on geomechanical tests of site-specific salt rock to probe the cyclic behaviors of salt precisely both beneath and above the dilatancy boundary, including reverse (inverse transient) creep, the Bauschinger effect, and damage-healing mechanism.

Conventional methods for extracting rare earth metals (REMs) from mined mineral ores are inefficient, expensive, and environmentally damaging. Recent discovery of lanmodulin (LanM), a protein that coordinates REMs with high-affinity and selectivity over competing ions, provides inspiration for new REM refinement methods. Here, we used quantum mechanical (QM) methods to investigate trivalent lanthanide cation (Ln3+) interactions with coordination systems representing bulk solvent water and protein binding sites. Energy decomposition analysis (EDA) showed differences in the energetic components of Ln3+ interaction with representatives of solvent (water, H2O) and protein binding sites (acetate, CH3COO-), highlighting the importance of accurate description of electrostatics and polarization in computational modeling of REM interactions with biological and bioinspired molecules. Relative binding free energies were obtained for Ln3+ with coordination complexes originating from binding sites in PDB structures of a lanthanum binding peptide (PDB entry 7CCO) and LanM, with explicit consideration of the first hydration shell waters, according to quasi-chemical theory (QCT). Beyond the first shell, the bulk solvent environment was represented with an implicit continuum model. Ln3+ interactions with (H2O)9 and both binding site models became more favorable, moving down the periodic series. This trend was more pronounced with the protein binding site models than with water, resulting in affinity increasing with periodic number, except for the last REM, Lu3+, which bound less favorably than the preceding element, Yb3+. Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln3+ relative binding free energies for the whole 7CCO peptide were reproduced. Conversely, the previously reported experimental data for LanM show a preference for the earlier lanthanides; this is likely due to longer-range interactions and cooperative effects, which are not represented by the reduced models. Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln3+ relative binding free energies for the whole 7CCO peptide were reproduced. In contrast to the previously reported experimental data for LanM, the peptide preferentially binds the earlier lanthanides. This difference likely arises due to longer-range interactions and cooperative effects not represented by the peptide. Further investigation of Ln3+ interactions with whole proteins using polarizable molecular mechanics models with explicit solvent is warranted to understand the influence of longer-ranged interactions, cooperativity, and bulk solvent. Nevertheless, the present work provides new insights into Ln3+ interactions with biomolecules and presents an effective computational platform for designing specific single-site REM binding peptides more efficiently.

This work presents the thermal runaway propagation model LIM1TR (Lithium-ion Modeling with 1-D Thermal Runaway) as an efficient tool to predict different cell-to-cell thermal runaway propagation scenarios. Here, we explored the vent gas volume production and reaction duration highlighting the relationship between these parameters and thermal runaway propagation due to convection by the vented gases. Two metrics based on gas production rate and heating rate are utilized as good indicators of the start and end of thermal runaway. LIM1TR results are compared with and validated by experiments from the literature for single-cell and multicell array experiments of 5 Ah and 10 Ah cells. By accounting for intraparticle diffusion of reacting species in the electrodes, we were able to capture the general dynamics of thermal runaway propagation and estimate acceptable reaction durations compared with the experimental values. Simulation results further demonstrated that varying heating modes lead to distinct reaction durations, consistent with experimental observations. Vent gas volume predictions indicate the need to consider both full and partial oxidation of the electrolyte. The outcomes of this work are building blocks for further investigations of module-to-module propagation by vented gases through convective heat transfer.

The space-charge-limited current (SCLC) in a vacuum diode is given by the Child-Langmuir law (CLL), whose electric potential ϕ(x) ∝ (x/D)^4/3, where x is the spatial coordinate across the gap and D is the gap separation distance. For a collisional diode, SCLC is given by the Mott-Gurney law (MGL) and ϕ(x) ∝ (x/D)^3/2. Here, we apply a capacitance argument for SCLC and use the transit time from a recent exact solution for collisional SCLC to show that ϕ(x) ∝ (x/D)ξ for a general collisional gap, where 4/3 ≤ ξ ≤ 3/2 . Furthermore, ξ is strictly a function of νT, where ν is the collision frequency and T is the electron transit time. Using this definition of ξ, we estimate the spatial dependence of the electron velocity and use the gap capacitance to derive an analytic equation for collisional SCLC that agrees within ~4.5% of the exact solution that requires solving parametrically through T. This analytic equation for general ξ asymptotically recovers the CLL as ν → 0 and the MGL as ν → ∞. As a result, matching these limits shows that ξ ≈ 1.40 and V ∝ D²ν² at the transition from a vacuum to a collisional diode for any device condition.

Low-velocity impact of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) composite laminates was studied experimentally and numerically. Hybrid laminates containing blocked layers of GFRP/CFRP/GFRP with all plies oriented at 0° were investigated. Relatively high impact energies were used to obtain full perforation of the laminate in a low-velocity impact setup. Numerical simulations were carried out using the in-house transient dynamics finite element code, Sierra/SM, developed at Sandia National Laboratories. A three-dimensional continuum damage model was used to describe the response of a woven composite ply. Two methods for handling delamination were considered and compared: (1) cohesive zone modeling and (2) continuum damage mechanics. The reduced model size achieved by omission of the cohesive zone elements produced acceptable results at reduced computational cost. The comparison between different modeling techniques can be used to inform modeling decisions relevant to low velocity impact scenarios. The modeling was validated by comparing with the experimental results and showed good agreement in terms of predicted damage mechanisms and impactor velocity and force histories.

Coastal communities face unique challenges in maintaining continuous service from critical infrastructure. This research advances capabilities for evaluating the impact of using wave energy to desalinate water on the resilience of coastal communities. The study focuses on the feasibility of using wave energy conversion to provide drinking water to communities in need and applying resilience metrics to quantify its impact on the community. To assess the feasibility of wave-powered desalination, this research couples the open-source software Wave Energy Converter SIMulator (WEC-Sim) and Water Network Tool for Resilience (WNTR). This research explores variations in both the wave resource (location, seasonality, and duration) and the ability to maintain drinking water service during a disruption scenario by applying the simulation framework to three case studies, which are based on communities in Puerto Rico. The simulation framework provides a contextualized assessment of the ability of wave-powered desalination to improve the resilience of coastal communities, which can serve as a methodology for future studies seeking the integration of wave-powered desalination with water distribution systems.

The (a)-type screw dislocations are known to be significant mediators of plasticity in hexagonal-close-packed (HCP) metals. These dislocations have polymorphic core structures, and subtle changes in the relative energies of these core structures are known to have a large impact on the dynamics of the dislocations. This work identifies a previously neglected long-range elastic interstitial-solute/dislocation interaction that influences the core structures. Essentially, interstitial solutes induce a change in the dislocation core structure to minimize the energy of interaction between the solutes and the dislocation. Molecular dynamics simulations, continuum linear elasticity, and statistical analysis show that this long-range interaction can locally alter the dislocation cores so that many different polymorphs appear along a single dislocation not only because of direct contact between interstitials and the dislocation core but also because of this long-range elastic interaction.