Publications

Abstract not provided.

A dry etching process to transfer the pattern of a photonic integrated circuit design for high-speed laser communications is described. The laser stack under consideration is a 3.2-µm-thick InGaAs/InAlAs/InAlGaAs epitaxial structure grown by molecular beam epitaxy. The etching was performed using Cl2-based inductively-coupled-plasma and reactive-ion-etching (ICP-RIE) reactors. Four different recipes are presented in two similar ICP-RIE reactors, with special attention paid to the etched features formed with various hard mask compositions, in-situ passivations, and process temperatures. The results indicate that it is possible to produce high-aspect-ratio features with sub-micron separation on this multilayer structure. Additionally, the results of the etching highlight the tradeoffs involved with the corresponding recipes.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Quantifying the radioactive sources present in gamma spectra is an ever-present and growing national security mission and a time-consuming process for human analysts. While machine learning models exist that are trained to estimate radioisotope proportions in gamma spectra, few address the eventual need to provide explanatory outputs beyond the estimation task. In this work, we develop two machine learning models for a NaI detector measurements: one to perform the estimation task, and the other to characterize the first model’s ability to provide reasonable estimates. To ensure the first model exhibits a behavior that can be characterized by the second model, the first model is trained using a custom, semi-supervised loss function which constrains proportion estimates to be explainable in terms of a spectral reconstruction. The second auxiliary model is an out-of-distribution detection function (a type of meta-model) leveraging the proportion estimates of the first model to identify when a spectrum is sufficiently unique from the training domain and thus is out-of-scope for the model. In demonstrating the efficacy of this approach, we encourage the use of meta-models to better explain ML outputs used in radiation detection and increase trust.

In the framework of SFERA-III WP10 Task3, ENEA has organized the 3D-shape round-robin (RR); the purpose is to compare the main geometrical parameters of 3D shape measurement of parabolic-trough (PT) reflective panels evaluated with the instruments adopted by each participant among: ENEA, DLR, F-ISE, NREL, and SANDIA. The last two institutions are outside of the EU, but benefited from the Transnational Access institute to visit several European laboratories, including the ENEA Casaccia research center where they accomplished some measurements with a portable experimental set-up. RR is based on the inter-laboratory circulation of 3 inner plus 3 outer PT panels. The start of the RR was delayed by the covid pandemic, then the circulation of the specimen-set and their measurement took more than one year. At the time of drafting this deliverable at the end of SFERA-III project, NREL has not yet completed the analysis of the measurements, making available only the deviations of the slopes. Therefore here will be reported only the preliminary results. The full comparison will be published as soon as possible, maybe in the open access venue Open Research Europe.

Developing atomically synergistic bifunctional catalysts relies on the creation of colocalized active atoms to facilitate distinct elementary steps in catalytic cycles. Herein, we show that the atomically-synergistic binuclear-site catalyst (ABC) consisting of Zn^δ+ -O-Cr⁶⁺ on zeolite SSZ-13 displays unique catalytic properties for iso-stoichiometric co-conversion of ethane and CO₂. Ethylene selectivity and utilization of converted CO₂ can reach 100 % and 99.0% under 500 °C at ethane conversion of 9.6%, respectively. In-situ/ex-situ spectroscopic studies and DFT calculations reveal atomic synergies between acidic Zn and redox Cr sites. Zn^δ+ (0 < δ < 2) sites facilitate β-C-H bond cleavage in ethane and the formation of Zn-H^δ- hydride, thereby the enhanced basicity promotes CO₂ adsorption/activation and prevents ethane C-C bond scission. The redox Cr site accelerates CO₂ dissociation by replenishing lattice oxygen and facilitates H₂O formation/desorption. This study presents the advantages of the ABC concept, paving the way for the rational design of novel advanced catalysts.

The high-pressure compaction of three dimensional granular packings is simulated using a bonded particle model (BPM) to capture linear elastic deformation. In the model, grains are represented by a collection of point particles connected by bonds. A simple multibody interaction is introduced to control Poisson's ratio and the arrangement of particles on the surface of a grain is varied to model both high- and low-frictional grains. At low pressures, the growth in packing fraction and coordination number follow the expected behavior near jamming and exhibit friction dependence. As the pressure increases, deviations from the low-pressure power-law scaling emerge after the packing fraction grows by approximately 0.1 and results from simulations with different friction coefficients converge. These results are compared to predictions from traditional discrete element method simulations which, depending on the definition of packing fraction and coordination number, may only differ by a factor of two. As grains deform under compaction, the average volumetric strain and asphericity, a measure of the change in the shape of grains, are found to grow as power laws and depend heavily on the Poisson's ratio of the constituent solid. Larger Poisson's ratios are associated with less volumetric strain and more asphericity and the apparent power-law exponent of the asphericity may vary. The elastic properties of the packed grains are also calculated as a function of packing fraction. In particular, we find the Poisson's ratio near jamming is 1/2 but decreases to around 1/4 before rising again as systems densify.

Granular matter takes many paths to pack in natural and industrial processes. The path influences the packing microstructure, particularly for frictional grains. We perform discrete element modeling simulations of different paths to construct packings of frictional spheres. Specifically, we explore four stress-controlled protocols implementing packing expansions and compressions in various combinations thereof. We characterize the eventual packed states through their dependence of the packing fraction and coordination number on packing pressure, identifying non-monotonicities with pressure that correlate with the fraction of frictional contacts. These stress-controlled, bulk-like particle simulations access very low-pressure packings, namely, the marginally stable limit, and demonstrate the strong protocol dependence of frictional granular matter.

Characterization techniques for powder feedstocks used in additive manufacturing (AM) have long been relied upon to describe the inputs to an AM workflow. However, functional gaps remain between tests to measure intrinsic and extrinsic properties with the direct performance within AM equipment. Furthermore, the common practice of reusing powder through multiple build cycles introduces effects and changes to feedstock performance that are otherwise difficult to measure quantitatively. Here, standardization and the development of new test methods have not kept pace with the rapid evolution of the AM industry and its reliance on highly coupled process-structure–property-performance relationships.

Here, this paper presents the conceptual design of a tension leg platform (TLP) for the ARCUS “towerless” vertical-axis wind turbine (VAWT). VAWTs are ideal for floating offshore sites and have several advantages over horizontal-axis wind turbines (HAWT) including reduced top mass, lower center of gravity, increased energy capture, and in turn lower cost. The towerless ARCUS VAWT drives these advantages further through increased structural efficiency and by enabling more optimized TLP designs with simplified installation procedures. For hull sizing, we have studied three turbine sizes with corresponding power ratings of 5.1 MW, 10.4 MW and 22.3 MW. The largest turbine was identified as having the greatest potential to reduce the levelized cost of energy (LCOE) and is the reference size used for the further detailed design process. The conceptual design of the VAWT TLP has been awarded with an ABS Approval in Principle Certificate. This paper contains brief analysis results and design findings for a TLP designed to house a VAWT, including the following topics: • Applicable Design Codes • Metocean Conditions • ARCUS Turbine Loads • Design Load Cases and Requirements - Pre-service TLP Stability - In-place TLP Global Performance • Platform Configurations, Hull Structure Scantling Design, Weight and CG Estimation, and General Arrangement Drawings • Hull Ballast Plan for both Pre-service and In-place Conditions • Pre-service Quayside Integration, Transportation and Wet Tow Stability Analysis • Global Performance Analysis for Motions and Tendon tensions • Summary of cost components and system levelized cost of energy

3D integration of multiple microelectronic devices improves size, weight, and power while increasing the number of interconnections between components. One integration method involves the use of metal bump bonds to connect devices and components on a common interposer platform. Significant variations in the coefficient of thermal expansion in such systems lead to stresses that can cause thermomechanical and electrical failures. More advanced characterization and failure analysis techniques are necessary to assess the bond quality between components. Frequency domain thermoreflectance (FDTR) is a nondestructive, noncontact testing method used to determine thermal properties in a sample by fitting the phase lag between an applied heat flux and the surface temperature response. The typical use of FDTR data involves fitting for thermal properties in geometries with a high degree of symmetry. In this work, finite element method simulations are performed using high performance computing codes to facilitate the modeling of samples with arbitrary geometric complexity. A gradient-based optimization technique is also presented to determine unknown thermal properties in a discretized domain. Using experimental FDTR data from a GaN-diamond sample, thermal conductivity is then determined in an unknown layer to provide a spatial map of bond quality at various points in the sample.

The epitaxial development and characterization of metamorphic “GaSb-on-silicon” buffers as substrates for antimonide devices is presented. The approach involves the growth of a spontaneously and fully relaxed GaSb metamorphic buffer in a primary epitaxial reactor, and use of the resulting “GaSb-on-silicon” wafer to grow subsequent layers in a secondary epitaxial reactor. The buffer growth involves four steps—silicon substrate preparation for oxide removal, nucleation of AlSb on silicon, growth of the GaSb buffer, and finally capping of the buffer to prevent oxidation. This approach on miscut silicon substrates leads to a buffer with negligible antiphase domain density. The growth of this buffer is based on inducing interfacial misfit dislocations between an AlSb nucleation layer and the underlying silicon substrate, which results in a fully relaxed GaSb buffer. A 1 μm thick GaSb layer buffer grown on silicon has ~9.2 × 10⁷ dislocations/cm². The complete lack of strain in the epitaxial structure allows subsequent growths to be accurately lattice matched, thus making the approach ideal for use as a substrate. Here we characterize the GaSb-on-silicon wafer using high-resolution x-ray diffraction and transmission electron microscopy. The concept’s feasibility is demonstrated by growing interband cascade light emitting devices on the GaSb-on-silicon wafer. The performance of the resulting LEDs on silicon approaches that of counterparts grown lattice matched on GaSb.

Alkali metals are among the most desirable negative electrodes for long duration energy storage due to their extremely high capacities. Currently, only high-temperature (>250 °C) batteries have successfully used alkali electrodes in commercial applications, due to limitations imposed by solid electrolytes, such as low conductivity at moderate temperatures and susceptibility to dendrites. Toward enabling the next generation of grid-scale, long duration batteries, we aim to develop molten sodium (Na) systems that operate with commercially attractive performance metrics including high current density (>100 mA cm-2), low temperature (<200 °C), and long discharge times (>12 h). In this work, we focus on the performance of NaSICON solid electrolytes in sodium symmetric cells at 110 °C. Specifically, we use a tin (Sn) coating on NaSICON to reduce interfacial resistance by a factor of 10, enabling molten Na symmetric cell operation with “discharge” durations up to 23 h at 100 mA cm-2 and 110 °C. Unidirectional galvanostatic testing shows a 70% overpotential reduction, and electrochemical impedance spectroscopy (EIS) highlights the reduction in interfacial resistance due to the Sn coating. Detailed scanning electron microscopy (SEM) and energy-dispersive spectroscopy (EDS) show that Sn-coated NaSICON enables current densities of up to 500 mA cm-2 at 110 °C by suppressing dendrite formation at the plating interface (Mode I). This analysis also provides a mechanistic understanding of dendrite formation at current densities up to 1000 mA cm-2, highlighting the importance of effective coatings that will enable advanced battery technologies for long-term energy storage.

We demonstrate an order of magnitude reduction in the sensitivity to optical crosstalk for neighboring trapped-ion qubits during simultaneous single-qubit gates driven with individual addressing beams. Gates are implemented via two-photon Raman transitions, where crosstalk is mitigated by offsetting the drive frequencies for each qubit to avoid first-order crosstalk effects from inter-beam two-photon resonance. The technique is simple to implement, and we find that phase-dependent crosstalk due to optical interference is reduced on the most impacted neighbor from a maximal fractional rotation error of 0.185 ( 4 ) without crosstalk mitigation to ≤ 0.006 with the mitigation strategy. Furthermore, we characterize first-order crosstalk in the two-qubit gate and avoid the resulting rotation errors for the arbitrary-axis Mølmer-Sørensen gate via a phase-agnostic composite gate. Finally, we demonstrate holistic system performance by constructing a composite CNOT gate using the improved single-qubit gates and phase-agnostic two-qubit gate. This work is done on the Quantum Scientific Computing Open User Testbed; however, our methods are widely applicable for individual addressing Raman gates and impose no significant overhead, enabling immediate improvement for quantum processors that incorporate this technique.

Because of the high-risk nature of emergencies and illegal activities at sea, it is critical that algorithms designed to detect anomalies from maritime traffic data be robust. However, there exist no publicly available maritime traffic data sets with real-world expert-labeled anomalies. As a result, most anomaly detection algorithms for maritime traffic are validated without ground truth. We introduce the HawaiiCoast_GT data set, the first ever publicly available automatic identification system (AIS) data set with a large corresponding set of true anomalous incidents. This data set—cleaned and curated from raw Bureau of Ocean Energy Management (BOEM) and National Oceanic and Atmospheric Administration (NOAA) automatic identification system (AIS) data—covers Hawaii’s coastal waters for four years (2017–2020) and contains 88,749,176 AIS points for a total of 2622 unique vessels. This includes 208 labeled tracks corresponding to 154 rigorously documented real-world incidents.

Here, we show that a laser at threshold can be utilized to generate the class of coherent and transform-limited waveforms (vt — z)^me^i(kz—ωt) at optical frequencies. We derive these properties analytically and demonstrate them in semiclassical time-domain laser simulations. We then utilize these waveforms to expand other waveforms with high modulation frequencies and demonstrate theoretically the feasibility of complex-frequency coherent absorption at optical frequencies, with efficient energy transduction and cavity loading. This approach has potential applications in quantum computing, photonic circuits, and biomedicine.

As global temperatures continue to rise, climate mitigation strategies such as stratospheric aerosol injections (SAI) are increasingly discussed, but the downstream effects of these strategies are not well understood. As such, there is interest in developing statistical methods to quantify the evolution of climate variable relationships during the time period surrounding an SAI. Feature importance applied to echo state network (ESN) models has been proposed as a way to understand the effects of SAI using a data-driven model. This approach depends on the ESN fitting the data well. If not, the feature importance may place importance on features that are not representative of the underlying relationships. Typically, time series prediction models such as ESNs are assessed using out-of-sample performance metrics that divide the times series into separate training and testing sets. However, this model assessment approach is geared towards forecasting applications and not scenarios such as the motivating SAI example where the objective is using a data driven model to capture variable relationships. Here, in this paper, we demonstrate a novel use of climate model replicates to investigate the applicability of the commonly used repeated hold-out model assessment approach for the SAI application. Simulations of an SAI are generated using a simplified climate model, and different initialization conditions are used to provide independent training and testing sets containing the same SAI event. The climate model replicates enable out-of-sample measures of model performance, which are compared to the single time series hold-out validation approach. For our case study, it is found that the repeated hold-out sample performance is comparable, but conservative, to the replicate out-of-sample performance when the training set contains enough time after the aerosol injection.

We investigate the interplay between the quantum Hall (QH) effect and superconductivity in InAs surface quantum well (SQW)/NbTiN heterostructures using a quantum point contact (QPC). We use QPC to control the proximity of the edge states to the superconductor. By measuring the upstream and downstream resistances of the device, we investigate the efficiency of Andreev conversion at the InAs/NbTiN interface. Our experimental data is analyzed using the Landauer-Büttiker formalism, generalized to allow for Andreev reflection processes. We show that by varying the voltage of the QPC, VQPC, the average Andreev reflection, A, at the QH-SC interface can be tuned from 50% to ∼10%. The evolution of A with VQPC extracted from the measurements exhibits plateaus separated by regions for which A varies continuously with VQPC. The presence of plateaus suggests that for some ranges of VQPC the QPC might be pinching off almost completely from the QH-SC interface some of the edge modes. Our work shows an experimental setup to control and advance the understanding of the complex interplay between superconductivity and QH effect in two-dimensional gas systems.

Efficient solution of the Vlasov equation, which can be up to six-dimensional, is key to the simulation of many difficult problems in plasma physics. The discontinuous Petrov-Galerkin (DPG) finite element methodology provides a framework for the development of stable (in the sense of Ladyzhenskaya–Babuška–Brezzi conditions) finite element formulations, with built-in mechanisms for adaptivity. While DPG has been studied extensively in the context of steady-state problems and to a lesser extent with space-time discretizations of transient problems, relatively little attention has been paid to time-marching approaches. In the present work, we study a first application of time-marching DPG to the Vlasov equation, using backward Euler for a Vlasov-Poisson discretization. We demonstrate adaptive mesh refinement for two problems: the two-stream instability problem, and a cold diode problem. We believe the present work is novel both in its application of unstructured adaptive mesh refinement (as opposed to block-structured adaptivity, which has been studied previously) in the context of Vlasov-Poisson, as well as in its application of DPG to the Vlasov-Poisson system. We also discuss extensive additions to the Camellia library in support of both the present formulation as well as extensions to higher dimensions, Maxwell equations, and space-time formulations.

Thermochemical air separation to produce high-purity N2 was demonstrated in a vertical tube reactor via a two-step reduction–oxidation cycle with an A-site substituted perovskite Ba0.15Sr0.85FeO3–δ (BSF1585). BSF1585 particles were synthesized and characterized in terms of their chemical, morphological, and thermophysical properties. A thermodynamic cycle model and sensitivity analysis using computational heat and mass transfer models of the reactor were used to select the system operating parameters for a concentrating solar thermal-driven process. Thermal reduction up to 800 °C in air and temperature-swing air separation from 800 °C to minimum temperatures between 400 and 600 °C were performed in the reactor containing a 35 g packed bed of BSF1585. The reactor was characterized for dispersion, and air separation was characterized via mass spectrometry. Gas measurements indicated that the reactor produced N2 with O2 impurity concentrations as low as 0.02 % for > 30 min of operation. A parametric study of air flow rates suggested that differences in observed and thermodynamically predicted O2 impurities were due to imperfect gas transport in the bed. Temperature swing reduction/oxidation cycling experiments between 800 and 400 °C in air were conducted with no statistically significant degradation in N2 purity over 50 cycles.

Optimization of the radiation pattern from a Bremsstrahlung target for a given application is possible by controlling the electron beam that impacts the high-atomic-number target. In this work, the electron beam is generated by a 13MV vacuum diode that terminates a coaxial magnetically insulted transmission line (MITL) on the HERMES-III machine at Sandia National Labs. Work by Sanford introduced a geometry for vacuum diodes that can control the flow within bounds. The "indented anode", as coined by Sanford, can straighten out the electron beam in a high-current diode that would otherwise be prone to beam pinching. A straighter beam will produce a more forwardly directed radiation pattern while a pinching electron beam will yield a focal point or hot spot on axis and a more diffuse radiation pattern. Either one of these may be desirable depending on the application. This work serves as a first attempt to optimize the radiation pattern in the former sense of collimating the radiation pattern given a limited parameter space. The optimization is attempted first using electromagnetic particle-in-cell simulations in the EMPIRE code suite. The setup of the models used in EMPIRE is discussed along with some basic theory behind some of the models used in the simulations such as anode heating and secondary ions. Theoretical work performed by Allen Garner and his students at Purdue is included here, which concerns the impact of collisions in these vacuum diodes. The EMPIRE simulations consider both an aggressive and a conservative design. The aggressive design is inherently riskier while the conservative design is chosen as something that, while still a risk, is more likely to perform as expected. The ultimate goal of this work was to validate the EMPIRE code results with experimental data. While the experiment that tested the diode designs proposed by the simulation results fell outside of the fiscal boundaries of this project (and for that reason the results of which are not included in this report), the hardware for the experiment was designed and drafted within those same fiscal boundaries, and is thus included in this report. However, there was yet another experiment performed in this project that tested a key feature of the diode: the hemispherical cathode. Those results are documented here as well, which show that the cathode tip is an important aspect to controlling the diode flow. A short series of simulations on this diode were also performed after the experiment in order to gain a better understanding of the effect of ions. on the flow pattern and faceplate dose profile.

Calcite (CaCO3) is one of the most common minerals in geologic and engineered systems. It is often in contact with aqueous solutions, causing chemically assisted fracture that is critical to understanding the stability of subsurface systems and manmade structures. Calcite fracture was evaluated with reactive molecular dynamics simulations, including the impacts of crack tip geometry (notch), the presence of water, and surface hydroxyl groups. Chemo-mechanical weakening was assessed by comparing the loads where fracture began to propagate. Our analyses show that in the presence of a notch, the load at which crack growth begins is lower, compared to the effect of water or surface hydroxyls. Additionally, the breaking of two adjacent Ca-O bonds is the kinetic limitation for crack initiation, since transiently broken bonds can reform, not resulting in crack growth. In aqueous environments, fresh (not hydroxylated) calcite surfaces exhibited water strengthening. Manual addition of H+ and/or OH- species on the (104) calcite surface resulted in chemo-mechanical weakening of calcite by 9%. Achieving full hydroxylation of the calcite surface was thermodynamically and kinetically limited, with only 0.17-0.01 OH/nm2 surface hydroxylation observed on the (104) surface at the end of the simulations. The limited reactivity of pure water with the calcite surface restricts the chemo-mechanical effects and suggests that reactions between physiosorbed water and localized structural defects may be dominating the chemo-mechanical process in the studies where water weakening has been reported.

We propose a method to extract the upper laser level’s (ULL’s) excess electronic temperature from the analysis of the maximum light output power (P_max) and current dynamic range ΔJ_d = (J_max – J_th) of terahertz quantum cascade lasers (THz QCLs). We validated this method, both through simulation and experiment, by applying it on THz QCLs supporting a clean three-level system. Detailed knowledge of electronic excess temperatures is of utmost importance in order to achieve high temperature performance of THz QCLs. Our method is simple and can be easily implemented, meaning an extraction of the excess electron temperature can be achieved without intensive experimental effort. This knowledge should pave the way toward improvement of the temperature performance of THz QCLs beyond the state-of-the-art.

Color centers in diamond are one of the most promising tools for quantum information science. Of particular interest is the use of single-crystal diamond membranes with nanoscale-thickness as hosts for color centers. Indeed, such structures guarantee a better integration with a variety of other quantum materials or devices, which can aid the development of diamond-based quantum technologies, from nanophotonics to quantum sensing. A common approach for membrane production is what is known as “smart-cut”, a process where membranes are exfoliated from a diamond substrate after the creation of a thin sub-surface amorphous carbon layer by He⁺ implantation. Due to the high ion fluence required, this process can be time-consuming. In this work, we demonstrated the production of thin diamond membranes by neon implantation of diamond substrates. With the target of obtaining membranes of ~200 nm thickness and finding the critical damage threshold, we implanted different diamonds with 300 keV Ne⁺ ions at different fluences. We characterized the structural properties of the implanted diamonds and the resulting membranes through SEM, Raman spectroscopy, and photoluminescence spectroscopy. We also found that a SRIM model based on a two-layer diamond/sp² -carbon target better describes ion implantation, allowing us to estimate the diamond critical damage threshold for Ne⁺ implantation. Compared to He⁺ smart-cut, the use of a heavier ion like Ne⁺ results in a ten-fold decrease in the ion fluence required to obtain diamond membranes and allows to obtain shallower smart-cuts, i.e. thinner membranes, at the same ion energy.

We hereby offer a comprehensive analysis of various factors that could potentially enable terahertz quantum cascade lasers (THz QCLs) to achieve room temperature performance. We thoroughly examine and integrate the latest findings from recent studies in the field. Our work goes beyond a mere analysis; it represents a nuanced and comprehensive exploration of the intricate factors influencing the performance of THz QCLs. Through a comprehensive and holistic approach, we propose novel insights that significantly contribute to advancing strategies for improving the temperature performance of THz QCLs. This all-encompassing perspective allows us not only to present a synthesis of existing knowledge but also to offer a fresh and nuanced strategy to improve the temperature performance of THz QCLs. We draw new conclusions from prior works, demonstrating that the key to enhancing THz QCL temperature performance involves not only optimizing interface quality but also strategically managing doping density, its spatial distribution, and profile. This is based on our results from different structures, such as two experimentally demonstrated devices: the spit-well resonant-phonon and the two-well injector direct-phonon schemes for THz QCLs, which allow efficient isolation of the laser levels from excited and continuum states. In these schemes, the doping profile has a setback that lessens the overlap of the doped region with the active laser states. Our work stands as a valuable resource for researchers seeking to gain a deeper understanding of the evolving landscape of THz technology. Furthermore, we present a novel strategy for future endeavors, providing an enhanced framework for continued exploration in this dynamic field. This strategy should pave the way to potentially reach higher temperatures than the latest records reached for T_max of THz QCLs.

Intermediate verification languages like Why3 and Boogie have made it much easier to build program verifiers, transforming the process into a logic compilation problem rather than a proof automation one. Why3 in particular implements a rich logic for program specification with polymorphism, algebraic data types, recursive functions and predicates, and inductive predicates; it translates this logic to over a dozen solvers and proof assistants. Accordingly, it serves as a backend for many tools, including Frama-C, EasyCrypt, and GNATProve for Ada SPARK. But how can we be sure that these tools are correct? The alternate foundational approach, taken by tools like VST and CakeML, provides strong guarantees by implementing the entire toolchain in a proof assistant, but these tools are harder to build and cannot directly take advantage of SMT solver automation. As a first step toward enabling automated tools with similar foundational guarantees, we give a formal semantics in Coq for the logic fragment of Why3. We show that our semantics are useful by giving a correct-by-construction natural deduction proof system for this logic, using this proof system to verify parts of Why3's standard library, and proving sound two of Why3's transformations used to convert terms and formulas into the simpler logics supported by the backend solvers.

Understanding pure H2 and H2/CH4 adsorption and diffusion in earth materials is one vital step toward a successful and safe H2 storage in depleted gas reservoirs. Despite recent research efforts such understanding is far from complete. In this work we first use Nuclear Magnetic Resonance (NMR) experiments to study the NMR response of injected H2 into Duvernay shale and Berea sandstone samples, representing materials in confining and storage zones. Then we use molecular simulations to investigate H2/CH4 competitive adsorption and diffusion in kerogen, a common component of shale. Our results indicate that in shale there are two H2 populations, i.e., free H2 and adsorbed H2, that yield very distinct NMR responses. However, only free gas presents in sandstone that yields a H2 NMR response similar to that of bulk H2. About 10 % of injected H2 can be lost due to adsorption/desorption hysteresis in shale, and no H2 loss (no hysteresis) is observed in sandstone. Our molecular simulation results support our NMR results that there are two H2 populations in nanoporous materials (kerogen). The simulation results also indicate that CH4 outcompetes H2 in adsorption onto kerogen, due to stronger CH4-kerogen interactions than H2-kerogen interactions. Nevertheless, in a depleted gas reservoir with low CH4 gas pressure, about ∼30 % of residual CH4 can be desorbed upon H2 injection. The simulation results also predict that H2 diffusion in porous kerogen is about one order of magnitude higher than that of CH4 and CO2. This work provides an understanding of H2/CH4 behaviors in deleted gas reservoirs upon H2 injection and predictions of H2 loss and CH4 desorption in H2 storage.

I experimentally investigated and modeled a proposed frequency-domain method for detecting and tracking cislunar spacecraft and near-earth asteroids using heliostat fields at night. Unlike imaging, which detects spacecraft and asteroids by their streak in star-fixed long-exposure photographs, the proposed detection method oscillates the orientation of heliostats concentrating light from the stellar field and measures the light's photocurrent power spectrum at sub-milliHertz resolution. If heliostat oscillation repetitively traces out a closed loop fixed to the stars, light from spacecraft or asteroids moving along that loop produce photocurrent at a frequency slightly shifted from starlight. The frequency shift is proportional to the spacecraft or asteroid's apparent angular rate relative to sidereal. Relative phase corresponds to relative angular position, enabling tracking. Since heliostats are inexpensive compared to an astronomical observatory and otherwise unused at night, the proposed method may cost-effectively augment observatory systems such as NASA's Asteroid Terrestrial-impact Last Alert System (ATLAS).

The primary goal of any laboratory test is to expose the unit-under-test to conservative realistic representations of a field environment. Satisfying this objective is not always straightforward due to laboratory equipment constraints. For vibration and shock tests performed on shakers over-testing and unrealistic failures can result because the control is a base acceleration and mechanical shakers have nearly infinite impedance. Force limiting and response limiting are relatively standard practices to reduce over-test risks in random-vibration testing. Shaker controller software generally has response limiting as a built-in capability and it is done without much user intervention since vibration control is a closed loop process. Limiting in shaker shocks is done for the same reasons, but because the duration of a shock is only a few milliseconds, limiting is a pre-planned user in the loop process. Shaker shock response limiting has been used for at least 30 years at Sandia National Laboratories, but it seems to be little known or used in industry. This objective of this paper is to re-introduce response limiting for shaker shocks to the aerospace community. The process is demonstrated on the BARBECUE testbed.

The siting of nuclear waste is a process that requires consideration of concerns of the public. This report demonstrates the significant potential for natural language processing techniques to gain insights into public narratives around “nuclear waste.” Specifically, the report highlights that the general discourse regarding “nuclear waste” within the news media has fluctuated in prevalence compared to “nuclear” topics broadly over recent years, with commonly mentioned entities reflecting a limited variety of geographies and stakeholders. General sentiments within the “nuclear waste” articles appear to use neutral language, suggesting that a scientific or “facts-only” framing of “waste”-related issues dominates coverage; however, the exact nuances should be further evaluated. The implications of a number of these insights about how nuclear waste is framed in traditional media (e.g., regarding emerging technologies, historical events, and specific organizations) are discussed. This report lays the groundwork for larger, more systematic research using, for example, transformer-based techniques and covariance analysis to better understand relationships among “nuclear waste” and other nuclear topics, sentiments of specific entities, and patterns across space and time (including in a particular region). By identifying priorities and knowledge needs, these data-driven methods can complement and inform engagement strategies that promote dialogue and mutual learning regarding nuclear waste.

Hail poses a significant threat to photovoltaic (PV) systems due to the potential for both cell and glass cracking. This work experimentally investigates hail-related failures in Glass/Backsheet and Glass/Glass PV modules with varying ice ball diameters and velocities. Post-impact Electroluminescence (EL) imaging revealed the damage extent and location, while high-speed Digital Image Correlation (DIC) measured the out-of-plane module displacements. The findings indicate that impacts of 20 J or less result in negligible damage to the modules tested. The thinner glass in Glass/Glass modules cracked at lower impact energies (-25 J) than Glass/Backsheet modules (-40 J). Furthermore, both module types showed cell and glass cracking at lower energies when impacted at the module's edges compared to central impacts. At the time of presentation, we will use DIC to determine if out-of-plane displacements are responsible for the impact location discrepancy and provide more insights into the mechanical response of hail impacted modules. This study provides essential insights into the correlation between impact energy, impact location, displacements, and resulting damage. The findings may inform critical decisions regarding module type, site selection, and module design to contribute to more reliable PV systems.

More than 90% of utility-scale photovoltaic (PV) power plants in the US use single-axis trackers (SATs) due to their potential for substantially higher power production over fixed-array systems. However, they are subject to software misconfigurations and mechanical failures, leading to suboptimal tracking accuracy. If failures are left undetected, the overall power yield of the PV power plant is reduced significantly. Robust detection and diagnosis of SAT faults is needed to minimize downtime and ensure continuous and efficient operation. This work presents analytic tools based on machine learning to detect deviations in SAT tracking performance and classify SAT faults.

For an Energy System to be truly equitable, it should provide affordable and reliable energy services to disadvantaged and underserved populations. Disadvantaged communities often face a combination of economic, social, health, and environmental burdens and may be geographically isolated (e.g., rural communities), which systematically limits their opportunity to fully participate in aspects of economic, social, and civic life.

Abstract not provided.

Visualization of flow structures within post-detonation fireballs has been performed for benchmark validation of numerical simulations. Custom pressed PETN explosives with a 12-mm diameter hemispherical form factor were used to produce a spherically symmetric post-detonation flow with low soot yield. Hydroxyl-radical planar laser induce fluorescence (OH-PLIF) was employed to visualize the structure ranging from approximately 10μs to 35μs after shock breakout from the explosive pellet. Fireball simulations were performed using the HyBurn Computational Fluid Dynamics (CFD) package. Experimental OH-PLIF results were compared to synthetic OH-PLIF from post-processing of CFD simulations. From the comparison of experimental and synthetic OH-PLIF images, CFD is shown to replicate much of the flow structure observed in the experiments, revealing potential differences in turbulent length scales and OH kinetics. Results provide significant advancement in experimental resolution of these harsh turbulent combustion environments and validate physical models thereof.

Multifidelity emulators have found wide-ranging applications in both forward and inverse problems within the computational sciences. Thanks to recent advancements in neural architectures, they provide significant flexibility for integrating information from multiple models, all while retaining substantial efficiency advantages over single-fidelity methods. In this context, existing neural multifidelity emulators operate by separately resolving the linear and nonlinear correlation between equally parameterized high-and low-fidelity approximants. However, many complex models ensembles in science and engineering applications only exhibit a limited degree of linear correlation between models. In such a case, the effectiveness of these approaches is impeded, i.e., larger datasets are needed to obtain satisfactory predictions. In this work, we present a general strategy that seeks to maximize the linear correlation between two models through input encoding. We showcase the effectiveness of our approach through six numerical test problems, and we show the ability of the proposed multifidelity emulator to accurately recover the high-fidelity model response under an increasing number of quasi-random samples. In our experiments, we show that input encoding produces in many cases emulators with significantly simpler nonlinear correlations. Finally, we demonstrate how the input encoding can be leveraged to facilitate the fusion of information between low-and high-fidelity models with dissimilar parametrization, i.e., situations in which the number of inputs is different between low-and high-fidelity models.

Point-source transient events (PSTEs) - optical events that are both extremely fast and extremely small - pose several challenges to an imaging system. Due to their speed, accurately characterizing such events often requires detectors with very high frame rates. Due to their size, accurately detecting such events requires maintaining coverage over an extended field-of-view, often through the use of imaging focal plane arrays (FPA) with a global shutter readout. Traditional imaging systems that meet these requirements are costly in terms of price, size, weight, power consumption, and data bandwidth, and there is a need for cheaper solutions with adequate temporal and spatial coverage. To address these issues, we develop a novel compressed sensing algorithm adapted to the rolling shutter readout of an imaging system. This approach enables reconstruction of a PSTE signature at the sampling rate of the rolling shutter, offering a 1-2 order of magnitude temporal speedup and a proportional reduction in data bandwidth. We present empirical results demonstrating accurate recovery of PSTEs using measurements that are spatially undersampled by a factor of 25, and our simulations show that, relative to other compressed sensing algorithms, our algorithm is both faster and yields higher quality reconstructions. We also present theoretical results characterizing our algorithm and corroborating simulations. The potential impact of our work includes the development of much faster, cheaper sensor solutions for PSTE detection and characterization.

Gallium nitride (GaN)-based nanoscale vacuum electron devices, which offer advantages of both traditional vacuum tube operation and modern solid-state technology, are attractive for radiation-hard applications due to the inherent radiation hardness of vacuum electron devices and the high radiation tolerance of GaN. Here, we investigate the radiation hardness of top-down fabricated n-GaN nanoscale vacuum electron diodes (NVEDs) irradiated with 2.5-MeV protons (p) at various doses. We observe a slight decrease in forward current and a slight increase in reverse leakage current as a function of cumulative protons fluence due to a dopant compensation effect. The NVEDs overall show excellent radiation hardness with no major change in electrical characteristics up to a cumulative fluence of 5E14 p/cm2, which is significantly higher than the existing state-of-the-art radiation-hardened devices to our knowledge. The results show promise for a new class of GaN-based nanoscale vacuum electron devices for use in harsh radiation environments and space applications.

Dendrites enable neurons to perform nonlinear operations. Existing silicon dendrite circuits sufficiently model passive and active characteristics, but do not exploit shunting inhibition as an active mechanism. We present a dendrite circuit implemented on a reconfigurable analog platform that uses active inhibitory conductance signals to modulate the circuit's membrane potential. We explore the potential use of this circuit for direction selectivity by emulating recent observations demonstrating a role for shunting inhibition in a directionally-selective Drosophila (Fruit Fly) neuron.

Anelastic strain recovery, the process of measuring the time dependent recovered strain after a core is cut at depth was utilized to make a measure of the in-situ properties stresses at depth at the FORGE (Frontier Observatory for Research in Geothermal Energy) site in Milford Utah. Core was collected from a region of well 16B at approximately 4860-4870 ft. Core was instrumented with strain gages within 10 hours of the core being cut. The relaxation of the cores was measured for approximately one month, and the results analyzed, which showed that the principal stresses were slightly off vertical, and magnitudes are close to equal.

Operation and control of a galvanically isolated three-phase AC-AC converter for solid state transformer applications is described. The converter regulates bidirectional power transfer by phase shifting voltages applied on either side of a high-frequency transformer. The circuit structure and control system are symmetrical around the transformer. Each side operates independently, enabling conversion between AC systems with differing voltage magnitude, phase angle, and frequency. This is achieved in a single conversion stage with low component count and high efficiency. The modulation strategy is discussed in detail and expressions describing the relationship between phase shift and power transfer are presented. Converter operation is demonstrated in a 3 kW hardware prototype.

The explosive BTF (benzotrifuroxan) is an interesting molecule for sub-millimeter studies of initiation and detonation. It has no hydrogen, thus no water in the detonation products and a subsequently high temperature in the reaction zone. The material has impact sensitivity that is comparable or less than that of PETN (pentaerythritol tetranitrate) and slightly greater than RDX, HMX, and CL-20. Physical vapor deposition (PVD) can be used to grow high-density films of pure explosives with precise control over geometry, and we apply this technique to BTF to study detonation and initiation behavior as a function of sample thickness. The geometrical effects on detonation and corner turning behavior are studied with the critical detonation thickness experiment and the micromushroom test, respectively. Initiation behavior is studied with the high-throughput initiation experiment. Vapor-deposited films of BTF show detonation failure, corner turning, and initiation consistent with a heterogeneous explosive. Scaling of failure thickness to failure diameter shows that BTF has a very small failure diameter.

Underground caverns in salt formations are promising geologic features to store hydrogen (H2) because of salt's extremely low permeability and self-healing behavior.Successful salt-cavern H2 storage schemes must maximize the efficiency of cyclic injection-production while minimizing H2 loss through adjacent damaged salt.The salt cavern storage community, however, has not fully understood the geomechanical behaviors of salt rocks driven by quick operation cycles of H2 injection-production, which may significantly impact the cost-effective storage-recovery performance.Our field-scale generic model captures the impact of combined drag and back stressing on the salt creep behavior corresponding to cycles of compression and extension, which may lead to substantial loss of cavern volumes over time and diminish the cavern performance for H2 storage.Our preliminary findings address that it is essential to develop a new salt constitutive model based on geomechanical tests of site-specific salt rock to probe the cyclic behaviors of salt both beneath and above the dilatancy boundary, including reverse (inverse transient) creep, the Bauschinger effect and fatigue.

Our present electric power grid maximizes spinning inertia of fossil fuel generators (inherent energy storage) to meet stability and performance requirements. Our goal is to begin to investigate the replacement of the large spinning inertia of fossil fuel generators with energy storage systems (ESS) including information flow as a necessary part of the renewable energy sources (RES) and subject to certain criteria. General criteria metrics include: energy storage, information flow, estimation, communication links, central versus decentralized, etc. Our focus is on evaluating the Fisher Information Equivalency (FIE) metric as a multi-criteria trade-off cost function for the minimization of ESS options and information flow. This paper begins with a formal conceptual definition of an infinite bus. Then a simple example of a One Machine Infinite Bus (OMIB) system with a Unified Power Flow Controller (UPFC) to demonstrate the FIE-based approach to minimize the ESS. A second more detailed example of several spinning machines are included with representative power electronic and ESS for RES that are attached to the electric power grid. A simple trade-study begins to highlight requirements to support large penetration of RES. Keep in mind for a large scale high penetration of RES will require large investments in ESS which we want to minimize.

The MACCS code was created by Sandia National Laboratories for the U.S. Nuclear Regulatory Commission and has been used for emergency planning, level 3 probabilistic risk assessments, consequence analyses and other scientific and regulatory research for over half a century. Specializing in modeling the transport of nuclear material into the environment, MACCS accounts for atmospheric transport and dispersion, wet and dry deposition, probabilistic treatment of meteorology, exposure pathways, varying protective actions for the emergency, intermediate and long-term phases, dosimetry, health effects (including but not limited to population dose, acute radiation injury and increased cancer risk), and economic impacts. Routine updates and recent enhancements to the MACCS code, such as the inclusion of a higher fidelity atmospheric transport and dispersion model, the addition of a new economic impact model, and the application of nearfield modeling, have continuously increased the codes capabilities in consequence analysis. Additionally, investigations of MACCS capabilities for advanced reactor applications have shown that MACCS can provide realistic and informative risk assessments for the new generation of reactor designs. Even so, areas of improvement as well as gaps have been identified that if resolved can increase the usefulness of MACCS in any application regarding a release of nuclear material into the environment.

Neuromorphic computing platforms hold the promise to dramatically reduce power requirements for calculations that are computationally intensive. One such application space is scientific machine learning (SciML). Techniques in this space use neural networks to approximate solutions of scientific problems. For instance, the popular physics-informed neural network (PINN) approximates the solution to a partial differential equation by using a trained feed-forward neural network, and injecting the knowledge of the physics through the loss function. Recent efforts have demonstrated how to convert a trained PINN to a spiking network architecture. In this work, we discuss our approach to quantization and implementation required to migrate these spiking PINNs to Intel's Loihi 2 neuromorphic hardware. We explore the effect of quantization on the model accuracy, as well as the energy and throughput characteristics of the implementation. It is our intent that this serve as a starting point for additional SciML implementations on neuromorphic hardware.