Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The classical Drude model provides an accurate description of the plasma resonance of three-dimensional materials, but only partially explains two-dimensional systems where quantum mechanical effects dominate such as P:δ layers - atomically thin sheets of phosphorus dopants in silicon that induce electronic properties beyond traditional doping. Previously it was shown that P:δ layers produce a distinct Drude tail feature in ellipsometry measurements. However, the ellipsometric spectra could not be properly fit by modeling the δ layer as a discrete layer of classical Drude metal. In particular, even for large broadening corresponding to extremely short relaxation times, a plasma resonance feature was anticipated but not evident in the experimental data. In this work, we develop a physically accurate description of this system, which reveals a general approach to designing thin films with intentionally suppressed plasma resonances. Our model takes into account the strong charge-density confinement and resulting quantum mechanical description of a P:δ layer. We show that the absence of a plasma resonance feature results from a combination of two factors: (i) the sharply varying charge-density profile due to strong confinement in the direction of growth; and (ii) the effective mass and relaxation time anisotropy due to valley degeneracy. The plasma resonance reappears when the atoms composing the δ layer are allowed to diffuse out from the plane of the layer, destroying its well-confined two-dimensional character that is critical to its distinctive electronic properties.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Reversible logic schemes using flux solitons (fluxons) on long Josephson junctions (LJJs) have recently been proposed. The attraction of the fluxon is that it propagates ballistically along an LJJ until it encounters a change in the character of the LJJ, often a designed circuit element. Logic gates involve fluxons interacting with circuit elements and with other fluxons. However, testing of ballistic fluxon circuits requires other circuits outside the logic family to direct and control fluxon motion. We discuss two such non-reversible fluxon control circuits. First, the polarity filter gate is a simple non-reversible gate that allows one polarity of fluxon to pass, while reflecting the other polarity. In the off state both polarities reflect. Second, the polarity separator generalizes on the polarity filter concept and allows separation of the two fluxon polarities into different LJJs. We discuss simulations of these structures and possible applications.

A robust data-management infrastructure is a key enabler for National Security Enterprise (NSE) capabilities in artificial intelligence and machine learning. This document describes efforts from a team of researchers at Sandia National Laboratories, Los Alamos National Laboratory, and Livermore National Laboratory to complete ASC Level II milestone #8854 “Assessment of Data-Management Infrastructure Needs for Production use of Advanced Machine learning and Artificial Intelligence.”

Abstract not provided.

Abstract not provided.

Abstract not provided.

The Xyce XDM Netlist Translator is a software solution developed by Sandia National Laboratories for translating SPICE netlists into a format readable by Xyce, Sandia’s internally developed SPICE-compatible simulation tool.

Abstract not provided.

Abstract not provided.

Abstract not provided.

An x-ray imaging scheme using spherically bent crystals was implemented on the Z-machine to image x rays emitted by the hot, dense plasma generated by a Magnetized Liner Inertial Fusion (MagLIF) target. This diagnostic relies on a spherically bent crystal to capture x-ray emission over a narrow spectral range (<15 eV), which is established by a limiting aperture placed on the Rowland circle. The spherical crystal optic provides the necessary high-throughput and large field-of-view required to produce a bright image over the entire, one-cm length of the emitting column of a plasma. The average spatial resolution was measured and determined to be 18 µm for the highest resolution configuration. With this resolution, the radial size of the stagnation column can be accurately determined and radial structures, such as bifurcations in the column, are clearly resolved. The success of the spherical-crystal imager has motivated the implementation of a new, two-crystal configuration for identifying sources of spectral line emission using a differential imaging technique.

Abstract not provided.

Abstract not provided.

Thermal-Hydrologic (TH) modeling of DECOVALEX 2023, Task C has continued in FY23. This report summarizes progress in TH modeling of Step 1c, with calibration modeling and the addition of shotcrete. The work involves 3-D modeling of the full-scale emplacement experiment at the Mont Terri Underground Rock Laboratory (Nagra, 2019). While Step 1 is focused on modeling the heating phase of the FE experiment with changes in pore pressure in the Opalinus clay resulting from heating, Step 1c is focused on calibration of models using available data.

The table presented below suggests the basic information that should be covered in a facility NMAC Plan for an NMAC program that is designed for nuclear security. The topics are appropriate for and should be addressed by all facilities in their NMAC Plans. They are appropriate for NMAC Plans for nuclear power plants, research reactors, fuel manufacturing facilities, facilities that produce medical isotopes, and other facilities. The difference is in the intensity with which the various measures are applied and the thoroughness of the description of the application (i.e., the program requirements). The robustness of a facility NMAC program and the content of its NMAC Plan should be graded in accordance with the type of facility and the category of its nuclear material.

A Machine and Deep Learning (MLDL) methodology is developed and applied to give a high fidelity, fast surrogate for 2D resistive MagnetoHydroDynamic (MHD) simulations of Magnetic Liner Inertial Fusion (MagLIF) implosions. The resistive MHD code GORGON is used to generate an ensemble of implosions with different liner aspect ratios, initial gas preheat temperatures (that is, different adiabats), and different liner perturbations. The liner density and magnetic field as functions of x, y, and z were generated. The Mallat Scattering Transformation (MST) is taken of the logarithm of both fields and a Principal Components Analysis (PCA) is done on the logarithm of the MST of both fields. The fields are projected onto the PCA vectors and a small number of these PCA vector components are kept. Singular Value Decompositions of the cross correlation of the input parameters to the output logarithm of the MST of the fields, and of the cross correlation of the SVD vector components to the PCA vector components are done. This allows the identification of the PCA vectors vis-a-vis the input parameters. Finally, a Multi Layer Perceptron (MLP) neural network with ReLU activation and a simple three layer encoder/decoder architecture is trained on this dataset to predict the PCA vector components of the fields as a function of time. Details of the implosion, stagnation, and the disassembly are well captured. Examination of the PCA vectors and a permutation importance analysis of the MLP show definitive evidence of an inverse turbulent cascade into a dipole emergent behavior. The orientation of the dipole is set by the initial liner perturbation. The analysis is repeated with a version of the MST which includes phase, called Wavelet Phase Harmonics (WPH). While WPH do not give the physical insight of the MST, they can and are inverted to give field configurations as a function of time, including field-to-field correlations.

Abstract not provided.

This paper describes a summer internship project undertaken at Sandia National Labs (SNL), both current status and future work. The project was to explore various machine learning approaches for use on turbulent flow data. Specifically, unsupervised classification of turbulent flow data was explored. First, the usage of models in this field is discussed, and several issues in the common usage of the models are identified. Solutions to these issues are then proposed, in the form of a Bayesian filtering approach which probabilistically incorporates multiple sources of data to improve confidence in a result. Several types of sensors are suggested for this method, the incorporation of which range from semi-supervised learning approaches to fully unsupervised. These approaches are then tested on several turbulent flow cases.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Several sources of interest often generate both low-frequency acoustic and seismic signals due to energy propagation through the atmosphere and the solid Earth. Seismic and acoustic observations are associated with a wide range of sources, including earthquakes, volcanoes, bolides, chemical and nuclear explosions, ocean noise, and others. The fusion of seismic and acoustic observations contributes to a better understanding of the source, both in terms of constraining source location and physics, as well as the seismic to acoustic coupling of energy. In this review, we summarize progress in seismoacoustic data processing, including recent developments in open-source data availability, low-cost seismic and acoustic sensors, and large-scale deployments of collocated sensors from 2010 to 2022. Similarly, we outline the recent advancements in modeling efforts for both source characteristics and propagation dynamics. Finally, we highlight the advantages of fusing multiphenomenological signals, focusing on current and future techniques to improve source detection, localization, and characterization efforts. This review aims to serve as a reference for seismologists, acousticians, and others within the growing field of seismoacoustics and multiphenomenology research.

An aircraft commander needs to be aware of weather phenomena that might be hazardous to his aircraft and mission. An important tool for this is airborne weather (WX) detection radar. The airborne WX radar needs to map weather for the aircraft commander that might be relevant to the safety of the aircraft, which involves both detecting a weather phenomenon, and to some extent seeing through it to detect weather phenomena behind it. Many factors influence the performance of an airborne WX radar

Abstract not provided.

Dominic D’Onofrio is currently a Junior studying Information and Technology at New Mexico Institute of Mining and Technology. He recently secured an Internship with NMCCoE, where he is involved with the TracerFIRE 12 project. Additionally, he is contributing to the load and security testing team by researching ways to implement pipelining and DevSecOps; this is his main project while he is at part time capacity for TracerFIRE 12. He is doing these projects to enhance his knowledge as a system administrator and gain a deeper understating of cybersecurity practices within national labs.

Abstract not provided.

Liquefied petroleum gas (LPG) is used in heating, cooking, and as a vehicle fuel (called autogas). A safety risk assessment may be needed to assess potential hazard scenarios and inform the regulations, codes, and standards that apply to LPG facilities, such as autogas refueling facilities. The frequency of unintended releases in an LPG system is an important aspect of a system quantitative risk assessment. This report documents estimation of leakage frequencies for individual components of LPG systems. These frequencies are described using uncertainty distributions obtained with Bayesian statistical methods, generic data, and LPG data which were publicly available. There was a lack of LPG data in the literature, so frequencies for most components were developed with generic data and should be used cautiously; without additional information about component leak frequencies in LPG systems, it is not known whether these generic frequencies may be conservative or non-conservative.

It is essential to Sandia National Laboratory’s continued success in scientific and technological advances and mission delivery to embrace a hybrid workforce culture under which current and future employees can thrive. This report focuses on the findings of the Hybrid Work Team for the Center for Computing Research, which met weekly from March to June 2023 and conducted a survey across the Center at Sandia. Conclusions in this report are drawn from the 9 authors of this report, which comprises the Hybrid Work Team, and 15 responses to a center-wide survey, as well as numerous conversations with colleagues. A major finding was widespread dissatisfaction with the quantity, execution, and tooling surrounding formal meetings with remote participants. While there was consensus that remote work enables people to produce high quality individual and technical work, there was also consensus that there was widespread social disconnect, with particular concern about hires that were made after the onset of the Covid-19 pandemic. There were many concerns about tooling and policy to facilitate remote collaboration both within Sandia and with its external collaborators. This report includes recommendations for mitigating these problems. For problems for which obvious recommendations cannot be made, ideas of what a successful solution might look like are presented.

Abstract not provided.

IT Service Management (ITSM) is an inimitable, ever-changing practice that is critical to all far-reaching organizations, including Sandia National Laboratories. With recent developments in technology and society – spurred by major events like the advent of ChatGPT and the pandemic – it is as important as ever to consider the implications and opportunities for ITSM. This report strategically synthesizes existing information about ITSM structures and examines the current state of Sandia’s service management entity: CCHD. Then, it looks at the emerging state of the world, culminating in CCHD-tailored recommendations for continual service improvement (CSI). Ultimately, the biggest matters to address are staff tenure, ticket documentation, and self-service facilitation. The next step might be to introduce AI to service channels in a non-system-invasive manner, namely a chatbot on the main CCHD page. All would serve to enhance end-user experiences, and by proxy Sandia’s output.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract: Advantages of the 2.5D HI (Heterogeneous Integration) electronics packaging of the power electronics compared to PCB packaging will be presented. Current 2.5D packaging effort using TSV (Through Silicon Via) will be presented in terms of fabrication, microstructural analysis, reliability, and thermal simulation.

Abstract not provided.

This work shows that the static random access memory (SRAM) error rate for a commercial 65-nm device in a dose rate environment can be highly dependent upon the integrated dose (dose rate × pulse duration). While the typical metric for such testing is dose rate upset (DRU) level in rad(Si)/s, a series of dose rate experiments at Little Mountain Test Facility (LMTF) shows dependence on the integrated dose. The error rate is also found to be dependent on the core voltage, and the preradiation value of the bits. We believe that these effects are explained by a well charge depletion caused by gamma ray photocurrent.

Abstract not provided.

Abstract not provided.

Abstract not provided.

In many areas of constrained optimization, representing all possible constraints that give rise to an accurate feasible region can be difficult and computationally prohibitive for online use. Satisfying feasibility constraints becomes more challenging in high-dimensional, non-convex regimes which are common in engineering applications. A prominent example that is explored in the manuscript is the security-constrained optimal power flow (SCOPF) problem, which minimizes power generation costs, while enforcing system feasibility under contingency failures in the transmission network. In its full form, this problem has been modeled as a nonlinear two-stage stochastic programming problem. In this work, we propose a hybrid structure that incorporates and takes advantage of both a high-fidelity physical model and fast machine learning surrogates. Neural network (NN) models have been shown to classify highly non-linear functions and can be trained offline but require large training sets. In this work, we present how model-guided sampling can efficiently create datasets that are highly informative to a NN classifier for non-convex functions. We show how the resultant NN surrogates can be integrated into a non-linear program as smooth, continuous functions to simultaneously optimize the objective function and enforce feasibility using existing non-linear solvers. Overall, this allows us to optimize instances of the SCOPF problem with an order of magnitude CPU improvement over existing methods.

Factor Fixing (FF) is a common method for reducing the number of model parameters to lower computational cost. FF typically starts with distinguishing the insensitive parameters from the sensitive and pursues uncertainty quantification (UQ) on the resulting reduced-order model, fixing each insensitive parameter at a fixed value. There is a need, however, to expand such a common approach to consider the effects of decision choices in the FF-UQ procedure on metrics of interest. Therefore, to guide the use of FF and increase confidence in the resulting dimension-reduced model, we propose a new adaptive framework consisting of four principles: (a) re-parameterize the model first to reduce obvious non-identifiable parameter combinations, (b) focus on decision relevance especially with respect to errors in quantities of interest (QoI), (c) conduct adaptive evaluation and robustness assessment of errors in the QoI across FF choices as sample size increases, and (d) reconsider whether fixing is warranted. The framework is demonstrated on a spatially-distributed water quality model. The error in estimates of QoI caused by FF can be estimated using a Polynomial Chaos Expansion (PCE) surrogate model. Built with 70 model runs, the surrogate is computationally inexpensive to evaluate and can provide global sensitivity indices for free. For the selected catchment, just two factors may provide an acceptably accurate estimate of model uncertainty in the average annual load of Total Suspended Solids (TSS), suggesting that reducing the uncertainty in these two parameters is a priority for future work before undertaking further formal uncertainty quantification.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Photocathodes based on GaAs and other III-V semiconductors are capable of producing highly spin-polarized electron beams. GaAs/GaAsP superlattice photocathodes exhibit high spin polarization; however, the quantum efficiency (QE) is limited to 1% or less. To increase the QE, we fabricated a GaAs/GaAsP superlattice photocathode with a Distributed Bragg Reflector (DBR) underneath. This configuration creates a Fabry-Pérot cavity between the DBR and GaAs surface, which enhances the absorption of incident light and, consequently, the QE. These photocathode structures were grown using molecular beam epitaxy and achieved record quantum efficiencies exceeding 15% and electron spin polarization of about 75% when illuminated with near-bandgap photon energies.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This special section of the Bulletin of the Seismological Society of America provides a broad overview on recent advances to the understanding of the seismoacoustic wavefield through 19 articles. Leveraging multiphenomenology datasets is instrumental for the continued success of future planetary missions, nuclear test ban treaty verification, and natural hazard monitoring. Progress in our theoretical understanding of mechanical coupling, advancements in coupled-media wave modeling, and developments of efficient multitechnology inversion procedures are key to fully exploiting geophysical datasets on Earth and beyond. We begin by highlighting papers describing experimental setups and instrumentation, followed by characterization of natural and anthropogenic sources of interest, and ending in new open-access datasets. Finally, we conclude with an overview of challenges that remain as well as some potential directions for future investigation within the growing multidisciplinary field of seismoacoustics.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Threaded fastener behavior can be an important aspect of complex component and system behavior, but there is no one-size-fits-all finite element analysis technique. Proper modeling of threaded fastener joints requires careful consideration of many details, from test setup and data acquisition to constitutive modeling and uncertainty quantification approaches. This report details analysis of a “mini-radax” bolted-joint exemplar where a Discrete-Direct uncertainty quantification approach is employed to evaluate margin of the component. The mini-radax geometry is tested to failure on a drop table, and single-coupon tests of individual fasteners serve as foundational data for the analysis. Analysis predictions complement the test data well and provide additional context for engineering decision-making.

The growth of helium bubbles impacts structural integrity of materials in nuclear applications. Understanding helium bubble nucleation and growth mechanisms is critical for improved material applications and aging predictions. Systematic molecular dynamics simulations have been performed to study helium bubble nucleation and growth mechanisms in Fe70Ni11Cr19 stainless steels. First, helium cluster diffusivities are calculated at a variety of helium cluster sizes and temperatures for systems with and without dislocations. Second, the process of diffusion of helium atoms to join existing helium bubbles is not deterministic and is hence studied using ensemble simulations for systems with and without vacancies, interstitials, and dislocations. We find that bubble nucleation depends on diffusion of not only single helium atoms, but also small helium clusters. Defects such as vacancies and dislocations can significantly impact the diffusion kinetics due to the trapping effects. Vacancies always increase the time for helium atoms to join existing bubbles due to the short-range trapping effect. This promotes bubble nucleation as opposed to bubble growth. Interestingly, dislocations can create a long-range trapping effect that reduces the time for helium atoms to join existing bubbles. This can promote bubble growth within a certain region near dislocations.

This report summarizes the fiscal year 2023 (FY23) status of the second phase of a series of borehole heater tests in salt at the Waste Isolation Pilot Plant (WIPP) funded by the Disposal Research and Development (R&D) program of the Spent Fuel & Waste Science and Technology (SFWST) office at the US Department of Energy’s Office of Nuclear Energy’s (DOE-NE) Office in the Spent Fuel and Waste Disposition (SFWD) program.

Abstract not provided.

Abstract not provided.

An additive manufacturing approach combining aerosol jet printing (AJP) and electrodeposition opens a new pathway to the production of lightweight coreless flyback transformer devices for power electronics. AJP of seed layers with resolution on the order of 30μm is combined with electrodeposition of Cu and Ni for decreased resistance. This combined approach addresses known shortcomings of AJP and electrodeposition. Nanoparticle inks used in AJP of metals have low conductivity versus bulk materials due to their high grain boundary resistance. There is a lack of readily available high-resolution patterning techniques for electrodeposition outside of expensive clean-room-based lithography techniques. Combining these two techniques enables the patterning of high-resolution, high-conductivity components. In this manuscript, we report on the construction of coreless flyback transformers consisting of two-layer primary and two-layer secondary spiral inductors separated by layers of a printed UV-curable dielectric. An input voltage of 17 V at 400 kHz was amplified to an output of 1250 V corresponding to a gain of 73.5. COMSOL modeling at the individual inductor level and at the transformer level was used to compare expected inductance, equivalent series resistance (ESR), and coupling with experimentally measured values.

Earthquakes have repeatedly been shown to produce inaudible acoustic signals (< 20 Hz), otherwise known as infrasound. These signals can propagate hundreds to thousands of kilometers and still be detected by ground-based infrasound arrays depending on the source strength, distance between source and receiver, and atmospheric conditions. Another type of signal arrival at infrasound arrays is the seismic induced motion of the sensor itself, or ground-motion-induced sensor noise. Measured acoustic and seismic waves produced by earthquakes can provide insight into properties of the earthquake such as magnitude, depth, and focal mechanism, as well as information about the local lithology and atmospheric conditions. Large earthquakes that produce strong acoustic signals detected at distances greater than 100 km are the most commonly studied; however, more recent studies have found that smaller magnitude earthquakes (Mw < 2:0) can be detected at short ranges. In that vein, this study will investigate the ability for a long-term deployment of infrasound sensors (deployed as part of the Source Physics Experiments [SPE] from 2014 to 2020) to detect both seismic and infrasonic signals from earthquakes at local ranges (< 50 km). Methods used include a combination of spectral analysis and automated array processing, supported by U.S. Geological Survey earthquake bulletins. This investigation revealed no clear acoustic detections for short range earthquakes. However, secondary infrasound from an Mw 7.1 earthquake over 200 km away was detected. Important insights were also made regarding the performance of the SPE networks including detections of other acoustic sources such as bolides and rocket launches. Finally, evaluation of the infrasound arrays is performed to provide insight into optimal deployments for targeting earthquake infrasound.

When high-energy-density materials are subjected to thermal or mechanical insults at extreme conditions (shock loading), a coupled response between the thermo-mechanical and chemical behaviors is systematically induced. Herein we develop a reaction model for the fast chemistry of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) at the mesoscopic scale, where the chemical behavior is determined by underlying microscopic reactive simulations. The slow carbon cluster formation is not discussed in the present work. All-atom reactive molecular dynamics (MD) simulations are performed with the ReaxFF potential, and a reduced-order chemical kinetics model for TATB is fitted to isothermal and adiabatic simulations of single crystal chemical decomposition. Unsupervised machine learning techniques based on non-negative matrix factorization are applied to MD trajectories to model the decomposition kinetics of TATB in terms of a four-component model. The associated heats of reaction are fit to the temperature evolution from adiabatic decomposition trajectories. Using a chemical species analysis, we show that non-negative matrix factorization captures the main chemical decomposition steps of TATB and provides an accurate estimation of their evolution with temperature. The final analytical formulation, coupled to a diffusion term, is incorporated into a continuum formalism, and simulation results are compared one-to-one against MD simulations of 1D reaction propagation along different crystallographic directions and with different initial temperatures. A good agreement is found for both the temporal and spatial evolution of the temperature field.

Abstract not provided.

Abstract not provided.

We present a graph neural network approach that fully automates the prediction of defect formation enthalpies for any crystallographic site from the ideal crystal structure, without the need to create defected atomic structure models as input. Here we used density functional theory reference data for vacancy defects in oxides, to train a defect graph neural network (dGNN) model that replaces the density functional theory supercell relaxations otherwise required for each symmetrically unique crystal site. Interfaced with thermodynamic calculations of reduction entropies and associated free energies, the dGNN model is applied to the screening of oxides in the Materials Project database, connecting the zero-kelvin defect enthalpies to high-temperature process conditions relevant for solar thermochemical hydrogen production and other energy applications. The dGNN approach is applicable to arbitrary structures with an accuracy limited principally by the amount and diversity of the training data, and it is generalizable to other defect types and advanced graph convolution architectures. It will help to tackle future materials discovery problems in clean energy and beyond.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.