Senanayake, Hasini S.; Wimalasiri, Pubudu N.; Godahewa, Sahan M.; Thompson, Ward H.; Greathouse, Jeffery A.
Here, we present a classical interatomic force field, silica-DDEC, to describe the interactions of amorphous and crystalline silica surfaces, parametrized using density functional theory-based charges. Charge schemes for silica surfaces were developed using the density-derived electrostatic and chemical (DDEC) method, which reproduces atomic charges of the periodic models as well as the electrostatic potential away from the atom sites. Lennard–Jones parameters were determined by requiring the correct description of (i) the amorphous silica density, coordination defects, and local coordination geometry, relative to experimental measurements, and (ii) water-silica interatomic distances compared with ab initio results. Deprotonated surface silanol sites are also described within the model based on DDEC charges. The result is a general electronic structure-derived model for describing fully flexible amorphous and crystalline silica surfaces and interactions of liquids with silica surfaces of varying structure and protonation state.
This report summarizes the activities performed by Sandia National Laboratories in FY23 to identify and test coating materials for the prevention, mitigation, and/or repair of potential chloride-induced stress corrosion cracking in spent nuclear fuel dry storage canisters. This work continues efforts by Sandia National Laboratories that are summarized in previous reports from FY20 through FY22 on the same topic. In FY23, Sandia National Laboratories, in collaboration with five industry partners through a memorandum of understanding, evaluated the physical, mechanical, and corrosion-resistance properties of eight different coating systems. The evaluation included thermal and radiation environments relevant to various time periods of storage for spent nuclear fuel canisters. The coating systems include polymeric (polyetherketoneketone, modified polyimide/polyurea, modified phenolic resin, epoxy), organic/inorganic ceramic hybrids (silane-based polyurethane hybrid and a quasi-ceramic sol-gel polyurethane hybrid), and coatings utilizing a Zn-rich primer applied to stainless steel coupons. The results and implications of these tests are summarized in this report. These analyses will be used to identify the most effective coatings for potential use on spent nuclear fuel dry storage canisters and to identify specific needs for further optimization of coating technologies for application on spent nuclear fuel canisters.
White, Zachary K.; Gott, Ryan P.; Bentz, Brian Z.; Xu, Kunning G.
Here we have observed the behavior of striations caused by ionization waves propagating in low-pressure helium DC discharges using the non-invasive laser-collision induced fluorescence (LCIF) diagnostic. To achieve this, we developed an analytic fit of collisional radiative model (CRM) predictions to interpret the LCIF data and recover quantitative two-dimensional spatial maps of the electron density, ne, and the ratios of LCIF emission states that can be correlated with Te with the use of accurate distribution functions at localized positions within striated helium discharges at 500 mTorr, 750 mTorr, and 1 Torr. To our knowledge, these are the first spatiotemporal, laser-based, experimental measurements of ne in DC striations. The ne and 447:588 ratio distributions align closely with striation theory. Constriction of the positive column appears to occur with decreased gas pressure, as shown by the radial ne distribution. We identify a transition from a slow ionization wave to a fast ionization wave between 750 mTorr and 1 Torr. These experiments validate our analytic fit of ne, allowing the implementation of an LCIF diagnostic in helium without the need to develop a CRM.
Chemically robust, low-power sensors are needed for the direct electrical detection of toxic gases. Metal-organic frameworks (MOFs) offer exceptional chemical and structural tunability to meet this challenge, though further understanding is needed regarding how coadsorbed gases influence or interfere with the electrical response. To probe the influence of competitive gases on trace NO2 detection in a simulated flue gas stream, a combined structure-property study integrating synchrotron powder diffraction and pair distribution function analyses was undertaken, to elucidate how structural changes associated with gas binding inside Ni-MOF-74 pores correlate with the electrical response from Ni-MOF-74-based sensors. Data were evaluated for 16 gas combinations of N2, NO2, SO2, CO2, and H2O at 50 °C. Fourier difference maps from a rigid-body Rietveld analysis showed that additional electron density localized around the Ni-MOF-74 lattice correlated with large decreases in Ni-MOF-74 film resistance of up to a factor of 6 × 103, observed only when NO2 was present. These changes in resistance were significantly amplified by the presence of competing gases, except for CO2. Without NO2, H2O rapidly (<120 s) produced small (1-3×) decreases in resistance, though this effect could be differentiated from the slower adsorption of NO2 by the evaluation of the MOF’s capacitance. Furthermore, samples exposed to H2O displayed a significant shift in lattice parameters toward a larger lattice and more diffuse charge density in the MOF pore. Evaluating the Ni-MOF-74 impedance in real time, NO2 adsorption was associated with two electrically distinct processes, the faster of which was inhibited by competitive adsorption of CO2. Together, this work points to the unique interaction of NO2 and other specific gases (e.g., H2O, SO2) with the MOF’s surface, leading to orders of magnitude decrease in MOF resistance and enhanced NO2 detection. Understanding and leveraging these coadsorbed gases will further improve the gas detection properties of MOF materials.
Over the past few decades, inorganic nitride materials have grown in importance in part due to their potential as catalysts for the synthesis of NH3, a key ingredient in fertilizer and precursor to industrial chemicals. Of particular interest are the ternary (ABN) or higher-order nitrides with high metal-to-nitrogen ratios that show promise in enhancing NH3 synthesis reaction rates and yields via heterogeneous catalysis or chemical looping. Although metal nitrides are predicted to be numerous, the stability of nitrogen triple bonds found in N2, especially in comparison to the metal-nitrogen bonds, has considerably hindered synthetic efforts to produce complex nitride compounds. In this study, we present an exhaustive down-selection process to identify ternary nitrides for a promising chemical looping NH3 production mechanism. We also report on a facile and efficient two-step synthesis method that can produce well-characterized η-carbide Co3Mo3N/Fe3Mo3N or filled β-manganese Ni2Mo3N ternaries, as well as their associated quaternary, (Co,Fe)3Mo3N, (Fe,Ni)2Mo3N, and (Co,Ni)2Mo3N, solid solutions. To further explore the quaternary space, syntheses of (Co,Ni)3Mo3N (Ni ≤ 10 mol %) and Co3(Mo,W)3N (W ≤ 10 mol %) were also investigated. The structures of the nitrides were characterized via X-ray powder diffraction. The morphology and compositions were characterized with scanning electron microscopy. The multitude of chemically unique, but structurally related, nitrides suggests that properties such as nitrogen activity may be tunable, making the materials of great interest for NH3 synthesis schemes.
We present an experimental study on a terahertz quantum cascade laser (THz QCL) design that combines both two-well injector and direct-phonon scattering schemes, i.e., a so-called two-well injector direct-phonon design. As a result of the two-well injector direct-phonon scheme presented here, the lasers benefit from both a direct phonon scattering scheme for the lower laser level depopulation and a setback for the doping profile that reduces the overlap of the doped region with active laser states. Additionally, our design also has efficient isolation of the active laser levels from excited and continuum states as indicated by negative differential resistance behavior all the way up to room temperature. This scheme serves as a good platform for improving the temperature performance of THz QCLs as indicated by the encouraging temperature performance results of the device with a relatively high doping level of 7.56 × 1010 cm−2 and Tmax ∼ 167 K. With the right optimization of the molecular beam epitaxy growth and interface quality, the injection coupling strength, and the doping density and its profile, the device could potentially reach higher temperatures than the latest records reached for the maximum operating temperature (Tmax) of THz QCLs.
Short-ranged and line-gapped non-Hermitian Hamiltonians have strong topological invariants given by an index of an associated Fredholm operator. It is shown how these invariants can be accessed via the signature of a suitable spectral localizer. Here, this numerical technique is implemented in an example with relevance to the design of topological photonic systems, such as topological lasers.
Traditional point-to-point communication sends data only after the entirety of the data is available. This includes situations where multiple actors (e.g., threads) contribute to the send buffer. As a result, cases where the completion times of these actors are widely distributed may be lost opportunities for optimization because data ready to be sent is waiting to be transmitted. Fine-grained communication exposes these opportunities by allowing buffers to be divided into elements that can then be sent independently (see e.g., Partitioned Communication in Message Passing Interface v4.0). While some research has been directed at exploring the utility of such 'early-bird' transmission, the overall search space for finding the best performing actor completion timings and element counts is large. In this work, we present an abstract model of fine-grained communication based on the LogGP model and a complementary benchmark. We use the model to explore actor completion timing scenarios and identify trends in communication behavior based on factors such as overall message size and delay between actor completions. We evaluate the benchmarks on three systems utilizing distinct network technologies and show that: (i) smaller numbers of elements are able to exploit most of the benefit of early-bird communication, (ii) performance benefit will depend non-trivially on application behavior, and (iii) benefits are highly network-dependent.
Early-bird communication is a communication/computation overlap technique that combines fine-grained communication with partitioned communication to improve application run-time. Communication is divided among the compute threads such that each individual thread can initiate transmission of its portion of the data as soon as it is complete rather than waiting for all of the threads. However, the benefit of early-bird communication depends on the completion timing of the individual threads. In this paper, we measure and evaluate the potential overlap, the idle time each thread experiences between finishing their computation and the final thread finishing. These measurements help us understand whether a given application could benefit from early-bird communication. We present our technique for gathering this data and evaluate data collected from three proxy applications: MiniFE, MiniMD, and MiniQMC. To characterize the behavior of these workloads, we study the thread timings at both a macro level, i.e., across all threads across all runs of an application, and a micro level, i.e., within a single process of a single run. We observe that these applications exhibit significantly different behavior. While MiniFE and MiniQMC appear to be well-suited for early-bird communication because of their wider thread distribution and more frequent laggard threads, the behavior of MiniMD may limit its ability to leverage early-bird communication.
UV photofragment spectroscopy and IR-UV double resonance methods are used to determine the structure and spectroscopic responses of a three-dimensional [2.2.2]-benzocryptand cage to the incorporation of a single K+ or Ba2+ imbedded inside it (labeled as K+-BzCrypt, Ba2+-BzCrypt). We studied the isolated ion-cryptand complex under cryo-cooled conditions, brought into the gas phase by nano-electrospray ionization. Incorporation of a phenyl ring in place of the central ethyl group in one of the three N-CH2-CH2-O-CH2-CH2-O-CH2-CH2-N chains provides a UV chromophore whose S0-S1 transition we probe. K+-BzCrypt and Ba2+-BzCrypt have their S0-S1 origin transitions at 35,925 and 36,446 cm-1, respectively, blue-shifted by 174 and 695 cm-1 from that of 1,2-dimethoxybenzene. These origins are used to excite a single conformation of each complex selectively and record their IR spectra using IR-UV dip spectroscopy. The alkyl CH stretch region (2800-3000 cm-1) is surprisingly sensitive to the presence and nature of the encapsulated ion. We carried out an exhaustive conformational search of cage conformations for K+-BzCrypt and Ba2+-BzCrypt, identifying two conformations (A and B) that lie below all others in energy. We extend our local mode anharmonic model of the CH stretch region to these strongly bound ion-cage complexes to predict conformation-specific alkyl CH stretch spectra, obtaining quantitative agreement with experiment for conformer A, the gas-phase global minimum. The large electrostatic effect of the charge on the O- and N-lone pairs affects the local mode frequencies of the CH2 groups adjacent to these atoms. The localized CH2 scissors modes are pushed up in frequency by the adjacent O/N-atoms so that their overtones have little effect on the alkyl CH stretch region. However, the localized CH2 wags are nearly degenerate and strongly coupled to one another, producing an array of delocalized wag normal modes, whose highest frequency members reach up above 1400 cm-1. As such, their overtones mix significantly with the CH stretch modes, most notably involving the CH2 symmetric stretch fundamentals of the central ethyl groups in the all-alkyl chains and the CH stretches adjacent to the N-atoms and antiperiplanar to the nitrogen lone pair.
Metamaterial resonators have become an efficient and versatile platform in the terahertz frequency range, finding applications in integrated optical devices, such as active modulators and detectors, and in fundamental research, e.g., ultrastrong light–matter investigations. Despite their growing use, characterization of modes supported by these subwavelength elements has proven to be challenging and it still relies on indirect observation of the collective far-field transmission/reflection properties of resonator arrays. Here, we present a broadband time-domain spectroscopic investigation of individual metamaterial resonators via a THz aperture scanning near-field microscope (a-SNOM). The time-domain a-SNOM allows the mapping and quantitative analysis of strongly confined modes supported by the resonators. In particular, a cross-polarized configuration presented here allows an investigation of weakly radiative modes. These results hold great potential to advance future metamaterial-based optoelectronic platforms for fundamental research in THz photonics.
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.
Evolutionary algorithms have been shown to be an effective method for training (or configuring) spiking neural networks. There are, however, challenges to developing accessible, scalable, and portable solutions. We present an extension to the Fugu framework that wraps the NEAT framework, bringing evolutionary algorithms to Fugu. This approach provides a flexible and customizable platform for optimizing network architectures, independent of fitness functions and input data structures. We leverage Fugu's computational graph approach to evaluate all members of a population in parallel. Additionally, as Fugu is platform-agnostic, this population can be evaluated in simulation or on neuromorphic hardware. We demonstrate our extension using several classification and agent-based tasks. One task illustrates how Fugu integration allows for spiking pre-processing to lower the search space dimensionality. We also provide some benchmark results using the Intel Loihi platform.
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.
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.
Wang, Qian; Guillaume, Joseph H.A.; Jakeman, John D.; Bennett, Frederick R.; Croke, Barry F.W.; Fu, Baihua; Yang, Tao; Jakeman, Anthony J.
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.
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.
In this paper, we highlight how computational properties of biological dendrites can be leveraged for neuromorphic applications. Specifically, we demonstrate analog silicon dendrites that support multiplication mediated by conductance-based input in an interception model inspired by the biological dragonfly. We also demonstrate spatiotemporal pattern recognition and direction selectivity using dendrites on the Loihi neuromorphic platform. These dendritic circuits can be assembled hierarchically as building blocks for classifying complex spatiotemporal patterns.
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.
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.
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.
Dannemann Dugick, Fransiska K.; Bishop, Jordan W.; Martire, Leo; Iezzi, Alexandra M.; Assink, Jelle D.; Brissaud, Quentin; Arrowsmith, Stephen
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.
IT Service Management (ITSM) is an inimitable, ever-changing practice that is critical to all far-reaching organizations, including Sandia. With recent developments in technology and society – spurred by major events like the advent of ChatGPT and the pandemic – it is important 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.
Recent work in neuromorphic computing has proposed a range of new architectures for Spiking Neural Network (SNN)-based systems. However, neuromorphic design lacks a framework to facilitate exploration of different SNN-based architectures and aid with early design decisions. While there are various SNN simulators, none can be used to rapidly estimate latency and energy of different spiking architectures. We show that while current spiking designs differ in implementation, they have common features which can be represented as a generic architecture template. We describe an initial version of a framework that simulates a range of neuromorphic architectures at an abstract time-step granularity. We demonstrate our simulator by modeling Intel's Loihi platform, estimating time-varying energy and latency with less than 10% mean error for various sizes of a two-layer SNN.
Deceglie, Michael G.; Anderson, Kevin; Fregosi, Daniel; Hobbs, William B.; Mikofski, Mark A.; Theristis, Marios; Meyers, Bennet E.
Photovoltaic systems may underperform expectations for several reasons, including inaccurate initial estimates, suboptimal operations and maintenance, or component degradation. Accurate assessment of these loss factors aids in addressing root causes of underperformance and in realizing accurate expectations and models. The performance loss rate (PLR) is a commonly cited high-level metric for the change in system output over time, but there is no precise, standard definition. Herein, an annualized definition of PLR that is inclusive of all loss factors and that can capture nonlinear changes to performance over time is proposed. The importance of distinguishing between recoverable and nonrecoverable losses which underly PLR is highlighted.
Self-healing or self-assembling power systems that rely on local measurements for decision making can provide significant resilience benefits, but they also must include safeguards that prevent the system from self-assembling into an undesirable configuration. One potential undesirable configuration would be the formation of closed loops for which the system was not designed, a situation that can arise any time that two intentional-island systems can be connected in more than one place, e.g., if tie-line breakers are included in the self-assembling system. This paper discusses the unintentional loop formation problem in self-assembling systems and presents a method for mitigating it. This method involves calculating the correlation or the mean absolute error (MAE) between the two local frequency measurements made on either side of a line relay. The correlation and MAE between these frequencies changes significantly between the loop and non-loop cases, and this difference can be used for loop detection. This article presents and explains the method in detail, presents evidence that the method's underlying assumptions are valid, and demonstrates in PSCAD two implementations of the method. The paper concludes with a discussion of the strengths and weaknesses of the proposed method.
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.
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.
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.
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.
Perhaps no single radar component has a more profound effect on Synthetic Aperture Radar (SAR) performance than the antenna. Especially for spaceborne SAR, one particular common design constraint for the antenna is the minimum antenna area constraint. While useful, it relies on a number of assumptions and approximations that may not always be valid or applicable. Indeed, useful operational systems have been built and flown that do not strictly adhere to this constraint. A closer examination
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.
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.
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.
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
Machine learning methodologies can provide insight into Brønsted-Guggenheim-Scatchard specific ion interaction theory (SIT) parameter values where experimental data availability may be limited. This study develops and executes machine learning frameworks to model the SIT interaction coefficient, ε. Key findings include successful estimations of ε via artificial neural networks using clustering and value prediction approaches. Applicability to other chemical parameters is also assessed briefly. Models developed here provide support for a use-case of machine learning in geologic nuclear waste disposal research applications, namely in predictions of chemical behaviors of high ionic strength solutions (i.e., subsurface brines).
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.
Obtaining lightweight and accurate approximations of discretized objective functional Hessians in inverse problems governed by partial differential equations (PDEs) is essential to make both deterministic and Bayesian statistical large-scale inverse problems computationally tractable. The cubic computational complexity of dense linear algebraic tasks, such as Cholesky factorization, that provide a means to sample Gaussian distributions and determine solutions of Newton linear systems is a computational bottleneck at large-scale. These tasks can be reduced to log-linear complexity by utilizing hierarchical off-diagonal low-rank (HODLR) matrix approximations. In this work, we show that a class of Hessians that arise from inverse problems governed by PDEs are well approximated by the HODLR matrix format. In particular, we study inverse problems governed by PDEs that model the instantaneous viscous flow of ice sheets. In these problems, we seek a spatially distributed basal sliding parameter field such that the flow predicted by the ice sheet model is consistent with ice sheet surface velocity observations. We demonstrate the use of HODLR Hessian approximation to efficiently sample the Laplace approximation of the posterior distribution with covariance further approximated by HODLR matrix compression. Computational studies are performed which illustrate ice sheet problem regimes for which the Gauss-Newton data-misfit Hessian is more efficiently approximated by the HODLR matrix format than the low-rank (LR) format. We then demonstrate that HODLR approximations can be favorable, when compared to global LR approximations, for large-scale problems by studying the data-misfit Hessian associated with inverse problems governed by the first-order Stokes flow model on the Humboldt glacier and Greenland ice sheet.
Oldfield, Ron A.; Allan, Benjamin A.; Doutriaux, Charles; Lewis, Katherine; Ahrens, James; Sims, Benjamin; Sweeney, Christine; Banesh, Divya; Wofford, Quincy
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.”
We report flow statistics and visualizations from molecular-gas-dynamics simulations using the direct simulation Monte Carlo (DSMC) method for turbulent Couette flow in a minimal domain where the lower wall is replaced by an idealized permeable fibrous substrate representative of thermal-protection-system materials for which the Knudsen number is O(10-1). Comparisons are made with smooth-wall DSMC simulations and smooth-wall direct numerical simulations (DNS) of the Navier-Stokes equations for the same conditions. Roughness, permeability, and noncontinuum effects are assessed. In the range of Reynolds numbers considered herein, the scalings of the skin friction on the permeable substrate and of the mean flow within the substrate suggest that they are dominated by viscous effects. While the regenerative cycle characteristic of smooth-wall turbulence remains intact for all cases considered, we observe that the near-wall velocity fluctuations are modulated by the permeable substrate with a wavelength equal to the pore spacing. Additionally, the flow within the substrate shows significant rarefaction effects, resulting in an apparent permeability that is 13% larger than the intrinsic permeability. In contrast, the smooth-wall DSMC and DNS simulations exhibit remarkably good agreement for the statistics examined, despite the Knudsen number based on the viscous length scale being as large as O(10-1). This latter result is at variance with classical estimates for the breakdown of the continuum assumption and calls for further investigations into the interaction of noncontinuum effects and turbulence.
Motion primitives allow for application of discrete search algorithms to rapidly produce trajectories in complex continuous space. The maneuver automaton (MA) provides an elegant formulation for creating a primitive library based on trims and maneuvers. However, performance is fundamentally limited by the contents of the primitive library. If the library is too sparse, performance can be poor in terms of path cost, whereas a library that is too large can increase run time. This work outlines new methods for using genetic algorithms to prune a primitive library. The proposed methods balance the path cost and planning time while maintaining the reachability of the MA. The genetic algorithm in this paper evaluates and mutates populations of motion primitive libraries to optimize both objectives. Here, we illustrate the performance of these methods with a simulated study using a nonlinear medium-fidelity F-16 model. We optimize a library with the presented algorithm for obstacle-free navigation and a nap-of-the-Earth navigation task. In the obstacle-free navigation task, we show a tradeoff of a 10.16% higher planning cost for a 96.63% improvement in run time. In the nap-of-the-Earth task, we show a tradeoff of a 9.712% higher planning cost for a 92.06% improvement in run time.
Understanding and controlling nanoscale interface phenomena, such as band bending and secondary phase formation, is crucial for electronic device optimization. In granular metal (GM) studies, where metal nanoparticles are embedded in an insulating matrix, the importance of interface phenomena is frequently neglected. Here, we demonstrate that GMs can serve as an exemplar system for evaluating the role of secondary phases at interfaces through a combination of x-ray photoemission spectroscopy (XPS) and electrical transport studies. We investigated SiNx as an alternative to more commonly used oxide-insulators, as SiNx-based GMs may enable high temperature applications when paired with refractory metals. Comparing Co-SiNx and Mo-SiNx GMs, we found that, in the tunneling-dominated insulating regime, Mo-SiNx had reduced metal-silicide formation and orders-of-magnitude lower conductivity. XPS measurements indicate that metal-silicide and metal-nitride formation are mitigatable concerns in Mo-SiNx. Given the metal-oxide formation seen in other GMs, SiNx is an appealing alternative for metals that readily oxidize. Furthermore, SiNx provides a path to metal-nitride nanostructures, potentially useful for various applications in plasmonics, optics, and sensing.
Recent work on atomic-precision dopant incorporation technologies has led to the creation of both boron and aluminum δ -doped layers in silicon with densities above the solid solubility limit. We use density functional theory to predict the band structure and effective mass values of such δ layers, first modeling them as ordered supercells. Structural relaxation is found to have a significant impact on the impurity band energies and effective masses of the boron layers, but not the aluminum layers. However, disorder in the δ layers is found to lead to a significant flattening of the bands in both cases. We calculate the local density of states and doping potential for these δ -doped layers, demonstrating that their influence is highly localized with spatial extents at most 4 nm. We conclude that acceptor δ -doped layers exhibit different electronic structure features dependent on both the dopant atom and spatial ordering. This suggests prospects for controlling the electronic properties of these layers if the local details of the incorporation chemistry can be fine-tuned.
Abu Rasel, Mdjafar; Schoell, Ryan; Al-Mamun, Nahid S.; Hattar, Khalid; Harris, Charles T.; Haque, Aman; Wolfe, Douglas E.; Ren, Fan; Pearton, Stephen J.
While radiation is known to degrade AlGaN/GaN high-electron-mobility transistors (HEMTs), the question remains on the extent of damage governed by the presence of an electrical field in the device. In this study, we induced displacement damage in HEMTs in both ON and OFF states by irradiating with 2.8 MeV Au4+ ion to fluence levels ranging from 1.72 × 10 10 to 3.745 × 10 13 ions cm−2, or 0.001-2 displacement per atom (dpa). Electrical measurement is done in situ, and high-resolution transmission electron microscopy (HRTEM), energy dispersive x-ray (EDX), geometrical phase analysis (GPA), and micro-Raman are performed on the highest fluence of Au4+ irradiated devices. The selected heavy ion irradiation causes cascade damage in the passivation, AlGaN, and GaN layers and at all associated interfaces. After just 0.1 dpa, the current density in the ON-mode device deteriorates by two orders of magnitude, whereas the OFF-mode device totally ceases to operate. Moreover, six orders of magnitude increase in leakage current and loss of gate control over the 2-dimensional electron gas channel are observed. GPA and Raman analysis reveal strain relaxation after a 2 dpa damage level in devices. Significant defects and intermixing of atoms near AlGaN/GaN interfaces and GaN layer are found from HRTEM and EDX analyses, which can substantially alter device characteristics and result in complete failure.
Experiments with trapped ions and neutral atoms typically employ optical modulators in order to control the phase, frequency, and amplitude of light directed to individual atoms. These elements are expensive, bulky, consume substantial power, and often rely on free-space I/O channels, all of which pose scaling challenges. To support many-ion systems like trapped-ion quantum computers or miniaturized deployable devices like clocks and sensors, these elements must ultimately be microfabricated, ideally monolithically with the trap to avoid losses associated with optical coupling between physically separate components. In this work we design, fabricate, and test an optical modulator capable of monolithic integration with a surface-electrode ion trap. These devices consist of piezo-optomechanical photonic integrated circuits configured as multi-stage Mach-Zehnder modulators that are used to control the intensity of light delivered to a single trapped ion on a separate chip. We use quantum tomography employing hundreds of multi-gate sequences to enhance the sensitivity of the fidelity to the types and magnitudes of gate errors relevant to quantum computing and better characterize the performance of the modulators, ultimately measuring single qubit gate fidelities that exceed 99.7%.
Material extrusion additive manufacturing (AM) has enabled an elegant fabrication pathway for a vast material library. Nonetheless, each material requires optimization of printing parameters generally determined through significant trial-and-error testing. To eliminate arduous, iteration-based optimization approaches, many researchers have used machine learning (ML) algorithms which provide opportunities for automated process optimization. In this work, we demonstrate the use of an ML-driven approach for real-time material extrusion print-parameter optimization through in-situ monitoring of printed line geometry. To do this, we use deep invertible neural networks (INNs) which can solve both forward and inverse, or optimization, problems using a single network. By combining in-situ computer vision and deep INNs, the printing parameters can be autonomously optimized to print a target line width in 1.2 s. Furthermore, defects that occur during printing can be rapidly identified and corrected autonomously. The methods developed and presented in this work eliminate user-intensive, time-consuming, and iterative parameter discovery approaches that currently limit accelerated implementation of extrusion-based AM processes. Furthermore, the presented approach can be generalized to provide real-time monitoring and optimization pathways for increasingly complex AM environments.