Parametrized artificial neural networks (ANNs) can be very expressive ansatzes for variational algorithms, reaching state-of-the-art energies on many quantum many-body Hamiltonians. Nevertheless, the training of the ANN can be slow and stymied by the presence of local minima in the parameter landscape. One approach to mitigate this issue is to use parallel tempering methods, and in this work, we focus on the role played by the temperature distribution of the parallel tempering replicas. Using an adaptive method that adjusts the temperatures in order to equate the exchange probability between neighboring replicas, we show that this temperature optimization can significantly increase the success rate of the variational algorithm with negligible computational cost by eliminating bottlenecks in the replicas' random walk. We demonstrate this using two different neural networks, a restricted Boltzmann machine and a feedforward network, which we use to study a toy problem based on a permutation invariant Hamiltonian with a pernicious local minimum and the J1-J2 model on a rectangular lattice.
This work presents a data-driven method for learning low-dimensional time-dependent physics-based surrogate models whose predictions are endowed with uncertainty estimates. We use the operator inference approach to model reduction that poses the problem of learning low-dimensional model terms as a regression of state space data and corresponding time derivatives by minimizing the residual of reduced system equations. Standard operator inference models perform well with accurate training data that are dense in time, but producing stable and accurate models when the state data are noisy and/or sparse in time remains a challenge. Another challenge is the lack of uncertainty estimation for the predictions from the operator inference models. Our approach addresses these challenges by incorporating Gaussian process surrogates into the operator inference framework to (1) probabilistically describe uncertainties in the state predictions and (2) procure analytical time derivative estimates with quantified uncertainties. The formulation leads to a generalized least-squares regression and, ultimately, reduced-order models that are described probabilistically with a closed-form expression for the posterior distribution of the operators. The resulting probabilistic surrogate model propagates uncertainties from the observed state data to reduced-order predictions. We demonstrate the method is effective for constructing low-dimensional models of two nonlinear partial differential equations representing a compressible flow and a nonlinear diffusion–reaction process, as well as for estimating the parameters of a low-dimensional system of nonlinear ordinary differential equations representing compartmental models in epidemiology.
Background/Objectives: Children’s biological age does not always correspond to their chronological age. In the case of BMI trajectories, this can appear as phase variation, which can be seen as shift, stretch, or shrinking between trajectories. With maturation thought of as a process moving towards the final state - adult BMI, we assessed whether children can be divided into latent groups reflecting similar maturational age of BMI. The groups were characterised by early factors and time-related features of the trajectories. Subjects/Methods: We used data from two general population birth cohort studies, Northern Finland Birth Cohorts 1966 and 1986 (NFBC1966 and NFBC1986). Height (n = 6329) and weight (n = 6568) measurements were interpolated in 34 shared time points using B-splines, and BMI values were calculated between 3 months to 16 years. Pairwise phase distances of 2999 females and 3163 males were used as a similarity measure in k-medoids clustering. Results: We identified three clusters of trajectories in females and males (Type 1: females, n = 1566, males, n = 1669; Type 2: females, n = 1028, males, n = 973; Type 3: females, n = 405, males, n = 521). Similar distinct timing patterns were identified in males and females. The clusters did not differ by sex, or early growth determinants studied. Conclusions: Trajectory cluster Type 1 reflected to the shape of what is typically illustrated as the childhood BMI trajectory in literature. However, the other two have not been identified previously. Type 2 pattern was more common in the NFBC1966 suggesting a generational shift in BMI maturational patterns.
Sandia National Laboratories has tested and evaluated a suite of four T120 broadband seismometers designed and manufactured by Nanometrics. Specifically, two T120 Horizon V2 sensors, one T120 Horizon V1 sensor and one T120 Slim Posthole (PH) sensor were evaluated. The purpose of this seismometer evaluation is to measure performance characteristics in areas such as power consumption, sensitivity, frequency response, full scale, self-noise, dynamic range, calibration system response, and passband. The T120 model of sensors are being evaluated to explore the potential for a future seismometer Type Approval process in the International Monitoring System (IMS) of the Comprehensive Nuclear-Test-Ban Treaty (CTBT).
Enblom, Veronica; Clulow, Rebecca; Ha, Tae J.; Witman, Matthew D.; Way, Lauren E.; Han, Sung J.; Brant Carvalho, Paulo H.B.; Stavila, Vitalie; Suh, Jin Y.; Sahlberg, Martin; Fadonougbo, Julien O.
A series of high entropy AB2-type Ti2-xZrxMnCrFeNi alloys (x = 0.6, 0.7, 0.8, 0.9, 1.0, 1.1 and 1.2) were synthesized to investigate their potential for hydrogen storage and chemical compression. The influence of the Ti/Zr ratio was explored in terms of structural, microstructural and thermodynamic properties. The storage capacity together with the reaction enthalpy and entropy changes of the synthesized high entropy alloys were compared to predictions from Machine Learning (ML) to investigate changes in these properties across the explored composition space. The results revealed that a decreasing Zr content consistently lowered the hydride formation enthalpy and increased the plateau pressure from 8 to >90 bar H2 at 25 °C, in good agreement with ML predictions. Selected compositions (x = 1.0 and 1.2) demonstrated reversible hydrogen storage capability over 150 cycles, with capacities of 1.34–1.40 wt % H2 and remarkable reaction kinetics (<4 min) at ambient temperature. These experimental and computational findings highlight the potential of this Laves-HEA system as tuneable, stable, and cost-effective materials suitable for long-term operations in stationary hydrogen storage and compression applications.
Carbon fiber provides opportunity to reduce weight in structural composites, including wind turbine blades, due to the material's superior specific stiffness and specific strength compared to alternatives. Despite these advantages, cost and compressive performance are considered weaknesses for carbon fiber products available today. Studies to produce low-cost carbon fiber alternatives, including the use of textile-derived precursor systems, have shown progress and merit through the DOE/ORNL low-cost carbon fiber initiatives. This work focuses on enabling increases in compressive strength through design of the carbon fiber geometry, applicable to both textile and conventional precursor systems, while also providing opportunities to reduce carbon fiber processing costs. Fiber-resin interface and fiber alignment are among the most frequently cited factors controlling composite compressive performance. However, it is believed that there is opportunity in traditionally unexplored routes to increasing compressive strength through alteration of the carbon fiber geometry by increasing the fiber area moment of inertia and/or the fiber perimeter and interfacial area. This paper presents initial results from manufacturing carbon fiber materials to assess the impacts of carbon fiber size on tested composite compressive performance with projected neutral or even beneficial impact on fiber and composite manufacturing economics. Carbon fiber systems with increasing size illustrate a favorable correlation for compressive performance greater than predicted from a micromechanical failure model. The manufacturing and mechanical test results support the hypothesis of this work that alterations to fiber geometry can be used to produce improvements of the compressive strength of carbon fiber reinforced polymers and provide incentive for related work in designing alternative shapes to further enhance compressive performance.
This work presents a data-driven method for learning low-dimensional time-dependent physics-based surrogate models whose predictions are endowed with uncertainty estimates. We use the operator inference approach to model reduction that poses the problem of learning low-dimensional model terms as a regression of state space data and corresponding time derivatives by minimizing the residual of reduced system equations. Standard operator inference models perform well with accurate training data that are dense in time, but producing stable and accurate models when the state data are noisy and/or sparse in time remains a challenge. Another challenge is the lack of uncertainty estimation for the predictions from the operator inference models. Our approach addresses these challenges by incorporating Gaussian process surrogates into the operator inference framework to (1) probabilistically describe uncertainties in the state predictions and (2) procure analytical time derivative estimates with quantified uncertainties. The formulation leads to a generalized least-squares regression and, ultimately, reduced-order models that are described probabilistically with a closed-form expression for the posterior distribution of the operators. The resulting probabilistic surrogate model propagates uncertainties from the observed state data to reduced-order predictions. We demonstrate the method is effective for constructing low-dimensional models of two nonlinear partial differential equations representing a compressible flow and a nonlinear diffusion–reaction process, as well as for estimating the parameters of a low-dimensional system of nonlinear ordinary differential equations representing compartmental models in epidemiology.
Reconstructing 3D granular microstructures within volumes of arbitrary geometries from limited 2D image data is crucial for predicting the material properties, as well as performances of structural components accounting for material microstructural effects. We present a novel generative learning framework that enables exascale reconstruction of granular microstructures within complex 3D geometric volumes. Building upon existing transfer learning techniques using pre-trained convolutional neural networks (CNN), we introduce several key innovations to overcome the difficulties inherent in arbitrary geometries. Our framework incorporates periodic boundary conditions using circular padding techniques, ensuring continuity and representativeness of the reconstructed microstructures. We also introduce a novel seamless transition reconstruction (STR) method that creates statistically equivalent transition zones to integrate multiple pre-existing 3D microstructure volumes. Based on STR, we propose a cost-effective strategy for reconstructing microstructures within complex geometric volumes, minimizing computational waste. Validation through numerical experiments using kinetic Monte Carlo simulations demonstrates accurate reproduction of grain statistics, including grain size distributions and morphology. A case study involving the reconstruction of a 4-blade propeller microstructure illustrates the method's capability to efficiently handle complex geometries. The proposed framework significantly reduces computational demands while maintaining high reconstruction quality, paving the way for scalable microstructure reconstruction in materials design and analysis.
We present the first fully curated, publicly accessible archive of infrasonic records from ten large bolide events documented by the U.S. Air Force Technical Applications Center’s global microbarometer network between 1960 and 1972. Captured on analog strip-chart paper, these waveforms predate modern digital arrays and space-based sensors, making them a unique window on meteoroid activity in the mid-twentieth century. Prior studies drew important scientific conclusions from the records but released only limited artifacts, chiefly period–amplitude tables and unprocessed scans, leaving the underlying data inaccessible for independent study. The present release transforms those limited excerpts into a research-ready resource. By capturing ten large events in the mid-20th century, the dataset constitutes a critical reference point for assessing bolide activity before the advent of modern space-based and digital ground-based monitoring. The multi-year coverage and worldwide distribution of events provide a valuable reference for comparing past and more recent detections, facilitating assessments of long-term flux and the dynamics of acoustic wave propagation in Earth’s atmosphere. The dataset’s availability in a consolidated format ensures straightforward access to waveforms and derived measurements, supporting a wide range of scientific inquiries into bolide physics and infrasound monitoring. By preserving these historical acoustic observations, the collection maintains a significant record of mid-20th-century meteoroid entries. It thereby establishes a basis for further refinement of impact hazard evaluations, contributes to historical continuity in atmospheric observation, and enriches the study of meteoroid-generated infrasound signals on a global scale. Dataset: DOI data identification number: 10.7910/DVN/VGSN7Q (direct URL to data: https://doi.org/10.7910/DVN/VGSN7Q). Dataset License: CC-BY-NC
We present theories and techniques to support safe, predictable, and efficient impact testing at the Sandia National Laboratories TA-III Mechanical Shock Complex (MSC), located in Building 6570. The two testing instruments used at the MSC are the 20-inch actuator, which is a launcher that propels sleds down an indoor track, and an indoor 6-in bore compressed-gas powered cannon. The theories and techniques presented derive from basic mechanical engineering principles but are significantly informe
We introduce the Poisson tensor completion (PTC) estimator, a non-parametric differential entropy estimator. The PTC estimator leverages inter-sample relationships to compute a low-rank Poisson tensor decomposition of the frequency histogram. Our crucial observation is that the histogram bins are an instance of a space partitioning of counts and thus can be identified with a spatial Poisson process. The Poisson tensor decomposition leads to a completion of the intensity measure over all bins—including those containing few to no samples—and leads to our proposed PTC differential entropy estimator. A Poisson tensor decomposition models the underlying distribution of the count data and guarantees non-negative estimated values and so can be safely used directly in entropy estimation. Our estimator is the first tensor-based estimator that exploits the underlying spatial Poisson process related to the histogram explicitly when estimating the probability density with low-rank tensor decompositions for the purpose of tensor completion. Furthermore, we demonstrate that our PTC estimator is a substantial improvement over standard histogram-based estimators for sub-Gaussian probability distributions because of the concentration of norm phenomenon.
DOE maintains an up-to-date documentation of the number of available full drawdowns of each of the caverns at the Strategic Petroleum Reserve (SPR), where an available full drawdown is defined as the ability to move all the oil out from a cavern at full drawdown rates. This report covers impacts on drawdown availability due to SPR operations during Calendar Year 2024.
On October 16-17, 2024, in Long Beach, California, the U.S. DOE Bioenergy Technologies Office (BETO), in collaboration with Sandia National Laboratories, NREL, PNNL, and LANL, convened experts from industry and academia to explore the potential for utilizing existing petroleum refinery infrastructure to expand the bioeconomy. Although refinery integration is not a new concept, recent trends have revitalized interest in this topic.
Cheung, Lawrence; Kuhn, Michael B.; Henry De Frahan, Marc T.; Mohan, Prakash; Deskos, Georgios; Churchfield, Matt; Sharma, Ashesh; Almgren, Ann; Ananthan, Shreyas; Brazell, Michael J.; Martinez-Tossas, Luis A.; Thedin, Regis; Rood, Jon; Sakievich, Philip; Vijayakumar, Ganesh; Zhang, Weiqun; Sprague, Michael
We present AMR-Wind, a verified and validated high-fidelity computational-fluid-dynamics code for wind farm flows. AMR-Wind is a block-structured, adaptive-mesh, incompressible-flow solver that enables predictive simulations of the atmospheric boundary layer and wind plants. It is a highly scalable code designed for parallel high-performance computing with a specific focus on performance portability for current and future computing architectures, including graphical processing units (GPUs). In this paper, we detail the governing equations, the numerical methods, and the turbine models. Establishing a foundation for the correctness of the code, we present the results of formal verification and validation. The verification studies, which include a novel actuator line test case, indicate that AMR-Wind is spatially and temporally second-order accurate. The validation studies demonstrate that the key physics capabilities implemented in the code, including actuator disk models, actuator line models, turbulence models, and large eddy simulation (LES) models for atmospheric boundary layers, perform well in comparison to reference data from established computational tools and theory. We conclude with a demonstration simulation of a 12-turbine wind farm operating in a turbulent atmospheric boundary layer, detailing computational performance and realistic wake interactions.
Sputter-deposited, nanocrystalline Cu-Ag thin films produced across a broad compositional and deposition-parameter space were evaluated to unravel the process-structure-property relationships important for creating hard, conductive electrical contacts and coatings. Combinatorial deposition involving pulsed direct current magnetron sputtering of elemental targets enabled swift examination of nearly the full range of alloy compositions and a relevant portion of deposition atomistics. Several high-throughput characterization modalities were employed to evaluate the chemistry, structure, and properties of the films. The resultant hardness, modulus, film density, crystal texture, and resistivity were analyzed in terms of key deposition characteristics (incident atom kinetic energy and incidence angle) predicted by binary-collision, kinematic Monte Carlo simulations. The study revealed improved hardness, parabolic resistivity dependence on composition, and compositional and process dependencies of film tarnishing. The results are discussed in the context of variations in microstructure and film density. Transmission electron microscopy and X-ray diffraction demonstrate several forms of compositional variation including solute segregation to grain boundaries as well as periodic, intragranular compositional modulations. Annealing of a Cu-rich alloy film exhibiting grain boundary segregation showed that this as-deposited, compositional variation is not stable above 100 °C. Finally, the Cu-Ag system is shown to have potential for hard, conductive, tarnish-resistant and room temperature-stable nanocrystalline thin films across the composition space.
Morag, Eilam; Li, Andrey; Kim, Anthony D.; Addamane, Sadhvikas J.; Williams, Benjamin S.
Quantum-cascade vertical-external-cavity surface-emitting-lasers (VECSELs) based on disordered amplifying metasurfaces are demonstrated and explored as potential broadband, multi-mode THz sources. The disorder is introduced along one spatial axis of the metasurface by pseudo-randomly varying the width of its resonant ridge antennas. Compared to a quantum-cascade (QC) VECSEL based on a uniform metasurface, the disordered structure supports much more localized transverse modes with reduced spatial overlap within the QC gain material. This localization is hypothesized to facilitate the spatial hole burning of the gain material and, therefore, enable multi-mode lasing, particularly for short cavities on the order of a few wavelengths. Several devices have been fabricated and shown to differ from uniform QC-VECSELs in a few key ways, possessing highly nonlinear light-current characteristics, angle-dependent emission spectra and broadband multi-mode lasing. At most, 17 modes are simultaneously observed, spanning 680 GHz. The number of VECSEL modes is shown to have an inverse relationship with the cavity length, which is attributed to increased diffractive losses in the open plano-plano cavity. As the cavity length is tuned, the device emits over a quasi-continuous band from 3.15 to 3.97 THz or up to 4.24 THz with a longitudinal mode hop.
Asymmetric energy exchange interactions, also known as Hatano-Nelson type couplings, enable the study of non-Hermitian physics and associated phenomena like the non-Hermitian skin effect and exceptional points (EPs). Since these interactions are by definition nonreciprocal, there have been very few options for implementations in integrated photonics. In this work, we show that asymmetric couplings are readily achievable in integrated photonic systems through time-domain dynamic modulation. We experimentally study this concept using a two-resonator photonic molecule produced in a lithium niobate on insulator platform that is electro-optically modulated by rf stimuli. We demonstrate the dynamic tuning of the Hatano-Nelson coupling between the resonators, surpassing the asymmetry that has been achieved in previous work, to reach an EP for the first time. We are additionally able to flip the relative sign of the couplings for opposite directions by going past the EP. Using this capability, we show that the through-chain transport can be configured to exhibit both giant (∼60 dB) optical contrast as well as photonic gyration or nonreciprocal π phase contrast.
The influence of the external cathodic area, stress intensity, and electrochemical polarization level on the resultant crack growth rate and crack tip pH and potential was studied for stainless steel 304L in a sodium chloride environment. Before testing, the validity of in situ crack tip property measurements using microprobes and drilled holes was supported through modeling efforts and confirmed with experimental observations. Here, it was found that increasing the R-ratio during corrosion fatigue increased the crack growth rate and decreased the measured crack tip pH of the 304L sample. Decreasing the electrochemical potential (i.e., cathodic polarization) of the bulk sample increased the crack tip pH and stopped measured crack extension under constant K conditions. Conversely, increasing the potential of the bulk sample decreased the crack tip pH and induced cracking under constant K conditions. Decreasing the cathodic area external to the crack increased the measured pH at the crack tip and decreased the crack growth rate under fatigue conditions. For the given configuration, environment, and alloy studied, fatigue crack growth is not sustained by cathodic reactions solely occurring in the crack wake and notch. These results have implications for the ability to predict atmospheric crack growth rates from laboratory, full immersion experiments.
Simulation-based approaches to microstructure generation can suffer from a variety of limitations, such as high memory usage, long computational times, and difficulties in generating complex geometries. Generative machine learning models present a way around these issues, but they have previously been limited by the fixed size of their generation area. We present a new microstructure generation methodology leveraging advances in inpainting using denoising diffusion models to overcome this generation area limitation. We show that microstructures generated with the presented methodology are statistically similar to grain structures generated with a kinetic Monte Carlo simulator, SPPARKS.
Baker, Andrew J.; Vishnubhotla, Sai B.; Karpe, Sanjana; Yang, Yahui; Veser, Gotz; Jacobs, Tevis B.
The adhesion of nanoparticles to their supports is key to their performance and stability. However, scientific advances in this area have been hampered by the difficulty of experimentally probing adhesion. To date, only a single technique has been developed that can directly measure nanoparticle adhesion, and this technique is inherently limited to monometallic systems. We present a versatile technique for the direct measurement of adhesion for bimetallic nanoparticle systems. This technique combines the spatial resolution of transmission electron microscopy with the force resolution of an atomic force microscope to probe individual, well-characterized nanoparticles. A first study of supported bimetallic nanoparticles provides new insights into the complex impact of alloying on nanoparticle adhesion, explained by charge transfer between constituent metals. The new experimental technique is readily extensible to study other multimetallic nanoparticle systems, including the effects of particle size, shape, and orientation, thus enabling advances in our understanding of nanoparticle physics.
Optimization-based coupling (OBC) is an attractive alternative to traditional Lagrange multiplier approaches in multiple modeling and simulation contexts. However, application of OBC to time-dependent problems has been hindered by the computational cost of finding the stationary points of the associated Lagrangian, which requires primal and adjoint solves. This issue can be mitigated by using OBC in conjunction with computationally efficient reduced order models (ROMs). To demonstrate the potential of this combination, in this paper, we develop an optimization-based ROM-ROM coupling for a transient advection-diffusion transmission problem. We pursue the “optimize-then-reduce” path toward solving the minimization problem at each time step and solve reduced space adjoint system of equations, where the main challenge in this formulation is the generation of adjoint snapshots and reduced bases for the adjoint systems required by the optimizer. One of the main contributions of the paper is a new technique for an efficient adjoint snapshot collection for gradient-based optimizers in the context of optimization-based ROM-ROM couplings. In conclusion, we present numerical studies demonstrating the accuracy of the approach along with comparison between various approaches for selecting a reduced order basis for the adjoint systems, including decay of snapshot energy, average iteration counts, and timings.
Oxygen vacancies in H f x Z r ( 1 − x ) O 2 (HZO) both contribute to stabilization of the ferroelectric orthorhombic phase and promote leakage pathways that limit the endurance of devices based on the material. For this reason, the defect states of oxygen vacancies were investigated using photoemission electron microscopy (PEEM) and photoluminescence spectroscopy (PL), as their concentration was varied via ex situ laser exposure. Following a controlled oxygen vacancy reduction via visible (2.54 eV) laser dosing of HZO, deep-ultraviolet (DUV, 5.82 eV) PEEM was used to spatially probe the resulting mid-gap defect states and work function. Work function was found to increase monotonically with the laser-induced reduction in oxygen vacancy concentration culminating in a total increase near 70 meV. The change implies a Fermi level shift toward the valence band as the total available electron-filled charge states are reduced with the removal of oxygen vacancies. A reduction in charge states is corroborated by the observed lessening of both photoemission and photoluminescence intensities after laser dosing. The deduced position of the Fermi level is within a band of near-conduction band defect states produced by oxygen vacancies that are linked to endurance limiting leakage currents. Together, these results directly identify the primary role of oxygen vacancies on the defect states in HZO while demonstrating that laser exposure can be used for their modification.
Five alternative metals are investigated as ohmic contacts to n-GaN including Cr/Au, Mo/Au, Pt/Au, Pd/Au, and Ge/Au. Ti-based contacts are traditionally used for ohmic contacts on n-GaN. However, conventional Ti/Al/Ni/Au metallization is found to be incompatible with a self-aligned process for GaN trench MOSFETs due to wet etch restrictions. Therefore, an alternative metallization is needed that is unreactive to the etch chemistry used in the self-aligned process. Additionally, the contact should remain ohmic following anneal at 900 °C so that contact formation can precede the anneal required for p-dopant activation. Here, in the present work, an n-GaN bilayer, consisting of a thin heavily doped contact layer (n0 = 1 × 1020 cm−3) atop a thick lesser doped layer, is used to demonstrate ohmic contacts of alternative metals with low specific contact resistance and extended thermal budget. Cr/Au ohmic contacts are demonstrated up to anneal temperatures of 800 °C, an increase of 200 °C compared to the highest known reports for Cr/Au contacts on n-GaN. Pt/Au metallization is demonstrated as an ohmic contact to n-GaN for the first time and exhibits true temperature-agnostic behavior up to anneal temperatures of 900 °C with specific contact resistance that is near parity with Ti/Al/Ni/Au. The temperature-agnostic behavior of Pt/Au ohmic contacts on the n-GaN bilayer, in addition to chemical compatibility with the self-aligned process, positions Pt/Au contacts as a key enabling element for self-aligned trench MOSFETs on GaN.
A finite element model of a 60-cell monocrystalline silicon glass-polymer photovoltaic module was simulated with ±1.0 kPa and ±2.4 kPa loads applied to the glass to calculate the deformation under load. Cell-to-cell displacements were used to approximate interconnect strain and stress. A mathematical fatigue cycle life relation was fitted to data for the interconnect material (copper), to generate a life prediction at each interconnect location based on the local stress means, reversal extents, and amplitudes. Interconnect stress was found to be significantly asymmetric about zero despite symmetric positive and negative module loads due to laminate thickness offsets about the neutral plane and the effects of module framing. Cycle life results indicated that interconnect fatigue failure was unlikely to occur over a 30-year lifetime of conservative wind and snow load cycles since the typical cell design feature of leaving some unconstrained length between the cell edge and first solder pad increases the effective gauge length and decreases the stress levels below the material endurance limit. Follow-up analyses found that 3.6 mm and 6.4 mm were the minimum unconstrained lengths required to survive the assumed lifetime of wind and snow cycles, respectively, confirming that typical industrial module constructions with 8–15 mm unconstrained lengths should survive conservatively. Notably, large magnitude, low-cycle snow loading was consistently the limiting factor requiring a longer unconstrained interconnect length. Finally, insights and workflows from this study inform module interconnection design limits for survival against mechanical fatigue in deployment environments.
Miniature atomic clocks based on the interrogation of the ground state hyperfine splitting of buffer gas cooled ions confined in radio frequency Paul traps have shown great promise as high precision prototype clocks. We report on the performance of two miniature ion trap vacuum packages after being sealed for as much as 10 years. We find the lifetime of the ions within the trap has increased over time for both traps and can be as long as 50 days. We form two clocks using the two traps and compare their relative frequency instability one with another to demonstrate a short-term instability of 5×10-13$τ$-1/2 integrating down to 1×10-14 after 2 ks of integration. The trapped ion lifetime and clock instability demonstrated by these miniature devices despite only being passively pumped for many years represents a critical advance toward their proliferation in the clock community.
NasGen provides a path for migration of structural models from Nastran bulk data format (BDF) into both an Exodus mesh file and an ASCII input file for Sierra Structural Dynamics (Salinas) and Solid Mechanics (Adagio). Many tools at Sandia National Labs (SNL) use the Exodus format. This document describes capabilities and limitations of the NasGen translation software.
Chapare virus (CHAPV) is an emerging New World arenavirus that is the causative agent of Chapare hemorrhagic fever (CHHF) responsible for recent outbreaks with alarmingly high case fatality rates in Bolivia near the Brazilian border. Here, we describe a nonhuman primate (NHP) model of CHHF infection which represents an essential tool to understand this emerging biological threat agent. Cynomolgus macaques challenged intravenously with CHAPV develop clinical disease, which recapitulates several key features of human CHHF. All subjects lost weight and had clinical scores following CHAPV challenge. Notably, one of four NHPs developed lethal disease with viral hepatitis and hemorrhagic features. Clinical chemistry and hematology revealed leukopenia, anemia, thrombocytopenia, and increased transaminase levels. In all four subjects, viremia was detectable for the first week following challenge and viral RNA was detectable in serum and many tissues persisting 35 days-post challenge. Several medical countermeasures (MCM) have efficacy against CHAPV infection in vitro, but the current model for MCM testing and approval of new drugs is reliant on the availability of animal models. This work lays the foundation for future CHHF MCM development.
Brittle behavior of metal alloys is often critical to modeling ballistic impact and penetration. The ALEGRA multiphysics finite element software incorporates calibrated models for the equation of state, elasticity, yield stress, plasticity and fracture, but simulations do not always capture expected metal fracture. Here we report concerted efforts to do so for one important case where experiments clearly show shear fractures: a tungsten sphere impacting a steel plate at various angles. Our best simulations show fractures that are qualitatively similar to experiments, but there are significant differences in quantitative metrics. Specifically, velocities of tracers used to quantify simulated plug parameters consistently fall short of measured plug velocities. Also, simulated plugs break apart more than expected from experimental evidence. We attribute these shortfalls to the lack of an explicit shear fracture mechanism in the material models, leading to over-estimated resistance to plug formation and movement.
Mageeney, Catherine M.; Touceda-Suarez, Maria; Perry, Matthew A.; Frizzo, Riccardo; Lotz-Mcmillen, John H.; Gilmore, Ruby A.; Bennett, Shauna M.; Basso, Jonelle T.R.; Fudyma, Jane D.; Geonczy, Sara E.; Gittrich, Marissa; Gogul, Grant; Hazard, Christina; Jameson, Ellie; Jiraska, Lucie; Johnson, Sarah S.; Kosmopoulos, James C.; Leleiwi, Ikaia; Bin MaBin; Millard, Andrew; Neri, Uri; Rodriguez-Ramos, Josue; Roux, Simon; Di TongDi; Wang, Yiling; Williamson, Kurt; Wu, Ruonan; Martins, Paula D.; Sapkota, Rumakanta; Emerson, Joanne B.; Trubl, Gareth
The research field of soil viral ecology continues to advance rapidly as the roles of viruses in the functioning of soil ecosystems are increasingly recognized. To address recent developments in the field, the second International Soil Virus Conference was held in Livermore, California, USA, from June 25 to 27th, 2024, providing soil viral ecologists the opportunity to share new findings and suggest guidelines for future research, while encouraging international scientific discussion and collaboration. The meeting was held in person with sessions simultaneously streamed online. Fifty researchers attended from ten different countries and spanned a wide range of subfields and career stages. A total of 21 oral presentations were presented, followed by discussions covering key themes in soil viral research. This report summarizes the main takeaways and recommendations from the talks and discussions.
Photovoltaic (PV) systems are essential for the transition to sustainable energy, reducing fossil fuel dependence and mitigating climate change. Although PV requires minimal land area — PV can meet the European Union's energy needs using only 0.26% of its land — space for deployment is often scarce in densely populated regions. Floating photovoltaics (FPV) offer an effective solution to land-use challenges by installing PV systems on floating structures in water bodies. FPV is a growing niche within PV with a cumulative installed capacity reaching 7.7 GW globally by 2023. Almost 90% of the installed FPV capacity is in Asia, with close to 50% of in China alone, while the Netherlands and France are the largest markets outside Asia. FPV shows strong potential to support climate targets, but still faces challenges like regulatory barriers, cost competitiveness compared to ground-based PV (GPV), and uncertainties about environmental impacts and system reliability. FPV systems are currently installed mainly on sheltered inland waters, such as quarry lakes, irrigation ponds and reservoirs. FPV technical standards are still being developed. Guidelines have been published by the World Bank, DNV, and Solar Power Europe, and emerging national standards from South Korea, China, and Singapore address design, components, and safety. The International Electrotechnical Commission (IEC) is working on formal standards for floats, mooring systems, and electrical connectors. However, the published best practices lack quantitative guidance for yield modelling and reliability, which this report aims to address. It provides data-driven insights, models, and parameters essential for accurate energy yield, reliability, and maintenance predictions over FPV systems' lifetimes.
Public-facing solar hosting capacity (HC) maps, which show the maximum amount of solar energy that can be installed at a location without adverse effects, have proven to be a key driver of solar soft cost reductions through a variety of pathways (e.g., streamlining interconnection, siting, and customer acquisition processes). However, current methods for generating HC maps require detailed grid models and time-consuming simulations that limit both their accuracy and scalability—today, only a handful out of almost 2,000 utilities provide these maps. This project developed and validated data-driven algorithms for calculating solar HC using data from AMI without the need of detailed grid models or simulations. The algorithms were validated on utility datasets and incorporated as an application into NRECA’s Open Modeling Framework (OMF.coop) for the over 260 coops and vendors throughout the US to use. The OMF is free and open-source for everyone.
Development of a defensible source-term model (STM), usual ly a thermodynamical model for radionuclide solubility calculations, is critical to a performance assessment (PA) of a geologic repository for nuclear waste disposal. Such a model is generally subjected to rigorous regulatory scrutiny. In this article, we highlight key guiding principles for STM model development and validation in nuclear waste management. We illustrate these principles by closely examining three recently developed thermodynamic models with the Pitzer formulism for aqueous H+—Nd3+—NO3−(—oxalate) systems in a reverse alphabetical order of the authors: the XW model developed by Xiong and Wang, the OWC model developed by Oakes et al., and the GLC model developed by Guignot et al., among which the XW model deals with trace activity coefficients for Nd(III), while the OWC and GLC models are for concentrated Nd(NO3)3 electrolyte solutions. The principles highlighted include the following: (1) Principle 1. Validation against independent experimental data: A model should be validated against experimental data or field observations that have not been used in the original model parameterization. We tested the XW model against multiple independent experimental data sets including electromotive force (EMF), solubility, water vapor, and water activity measurements. The results show that the XW model is accurate and valid for its intended use for predicting trace activity coefficients and therefore Nd solubility in repository environments. (2) Principle 2. Testing for relevant and sensitive variables: Solution pH is such a variable for an STM and easily acquirable. All three models are checked for their ability to predict pH conditions in Nd(NO3)3 electrolyte solutions. The OWC model fails to provide a reasonable estimate for solution pH conditions, thus casting serious doubt on its validity for a source-term calculation. In contrast, both the XW and GLC models predict close-to-neutral pH values, in agreement with experimental measurements. (3) Principle 3. Honoring physical constraints: Upon close examination, it is found that the Nd(III)-NO3 association schema in the OWC model suffers from two shortcomings. Firstly, its second stepwise stability constant for Nd(NO3)2+ (log K2) is much higher than the first stepwise stability constant for NdNO32+ (log K1), thus violating the general rule of (log K2–log K1) < 0, or (Formula presented.). Secondly, the OWC model predicts abnormally high activity coefficients for Nd(NO3)2+ (up to ~900) as the concentration increases. (4) Principle 4. Minimizing degrees of freedom for model fitting: The OWC model with nine fitted parameters is compared with the GLC model with five fitted parameters, as both models apply to the concentrated region for Nd(NO3)3 electrolyte solutions. The latter appears superior to the former because the latter can fit osmotic coefficient data equally well with fewer model parameters. The work presented here thus illustrates the salient points of geochemical model development, selection, and validation in nuclear waste management.
Tamper-indicating devices (TIDs), also known as seals, play a crucial role in various sectors including international nuclear safeguards, arms control, domestic security, and commercial products, by ensuring that monitored or high-value items are not accessed undetected. These devices do not block access but alert to unauthorized tampering. With adversaries' capabilities evolving, there's a pressing need for seals to advance in terms of effectiveness (e.g., better tamper indication and unique identification), and new technology can improve the efficiency of installation and verification. Passive loop seals, widely used in international nuclear safeguards to ensure that continuity of knowledge is maintained on declared items, face stringent International Atomic Energy Agency (IAEA) requirements that surpass those met by commercial products. The metal cup seal (Figure 1, left), a staple IAEA seal, is robust but requires significant resources for post-use verification – specifically, the seal’s unique identity can only be verified at IAEA headquarters after removal from facilities. Further, the seal has been in use for decades and seal types should periodically be replaced to counter adversarial efforts for defeating seals. In 2020, the IAEA outlined about 40 requirements for a new passive loop seal, aiming for in-situ verification, minimal external tool use, unique identification (UID), and clear tamper indication. In response, research and development efforts focused on creating a new passive loop seal that meets these criteria and in 2022 the IAEA announced the completion of the Field Verifiable Passive Loop Seal (FVPS) (Figure 1, right). Concurrently to the IAEA’s efforts, Sandia National Laboratories (SNL) and Oak Ridge National Laboratory (ORNL) designed, developed, and tested two seal versions – Puck and Puck/SAW, with Puck based on the IAEA’s requirements and including a novel visually-obvious tamper response, and Puck/SAW adding additional beneficial capabilities like the ability to receive a unique identifier from a standoff distance and monitoring the wire integrity. Puck/SAW was specifically designed and developed to address sealing applications in dry spent fuel storage facilities, where the number of sealed spent fuel containers results in heavy verification burden and inspector safety issues related to radiation exposure. These efforts are described in this Executive Summary.
Recent successes in the exploration, drilling, and discovery of geologic hydrogen have generated notable excitement. This new energy resource has the potential to make an important contribution to our nation’s energy supply, resiliency, and security. Contemporary studies of geologic hydrogen have a common theme of suggesting places where it might be found or even more specifically, what rocks in what geologic formations may contribute to its formation — either naturally or via artificially induced means. This vital ongoing body of work sets the stage for imagining what may be possible with vast available quantities of naturally occurring hydrogen in the subsurface. While acknowledging current approaches to characterizing geologic hydrogen, this report advances the discussion by suggesting next steps, including the critical science and engineering necessary to make geologic hydrogen an affordable and reliable part of the U.S. energy portfolio.
Autonomous manipulation is a challenging problem in field robotics due to uncertainty in object properties, constraints, and coupling phenomenon with robot control systems. Humans learn motion primitives over time to effectively interact with the environment. We postulate that autonomous manipulation can be enabled by basic sets of motion primitives as well, but do not necessitate mimicking human motion primitives. This work presents an approach to generalized optimal motion primitives using physics-informed neural networks. Our simulated and experimental results demonstrate that optimality is notionally maintained where the mean maximum observed final position percent error was 0.564% and the average mean error for all the trajectories was 1.53%. These results indicate that notional generalization is attained using a physics-informed neural network approach that enables near optimal real-time adaptation of primitive motion profiles.
This document details a data mining exercise that resulted in an exploratory dataset of publicly reported foreign (non-US) hypersonic vehicle test events. Using a combination of targeted English language searches and country-specific queries, the study aggregates information from digital news media, official press releases, and social media posts. The resulting list of events captures the publicly available accounts of foreign hypersonic tests, although it does not represent an exhaustive record. Limitations such as inconsistent reporting, translation challenges, and the inherently provisional nature of open-source data are acknowledged. This dataset serves as an initial reference point for further inquiries into high-speed atmospheric phenomena and may facilitate future efforts to correlate these events with geophysical measurements.
The Fusion of Simulation, Experiment, and Data (FuSED) team provides a set of tools for solving inverse problems in structural dynamics and thermal physics, and also sensor placement optimization via Optimal Experimental Design (OED). These methods are used for designing experiments, model calibration, and verification/validation analysis of systems. This document provides a user’s guide to the input for the three apps that are supported for these methods. Details of input specifications, output options, and optimization parameters are included.
A mesoscale model to predict helium bubble evolution is needed for tritium applications. Such a model requires that the conventional kinetic Monte Carlo (kMC) simulations be significantly accelerated. The objective of this report is to (a) highlight the concepts and mathematical expressions of the accelerated method for defect implementation that have not been published, (b) show an example input file to run the kMC code, and (c) provide suggestions on future improvement following my retirement.
