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.