Search NASA⌕ Search

SEARCH · Search NASA

Results for “Statistical Methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Antenna arrays for neutrino mass measurements with cyclotron radiation emission spectroscopy

Cyclotron Radiation Emission Spectroscopy (CRES) is a technique for precision measurements of kinetic energies of charged particles, pioneered by the Project 8 experiment to measure the neutrino mass using the tritium end-point method. It was recently employed for the first time to measure the molecular tritium spectrum and place a limit on the neutrino mass using a cubic-centimeter-scale detector. Future direct neutrino mass experiments are developing the technique to overcome the systematic and statistical limitations of current detectors. Here, this paper describes one such approach, namely the use of antenna arrays for CRES in free space. Phenomenology, detector design, simulation, and performance estimates are discussed, culminating with an example design with a projected sensitivity of 𝑚 𝛽 < ⁢0.04 eV/𝑐 2 . Prototype antenna array measurements are also shown for a demonstrator-scale setup as a benchmark for the simulation. By consolidating these results, this paper serves as a comprehensive reference for the development and performance of antenna arrays for CRES.

Physics - Nuclear physics and radiation physics↗
DOE OSTI2024

Testing hadronic-model predictions of depth of maximum of air-shower profiles and ground-particle signals using hybrid data of the Pierre Auger Observatory

We test the predictions of hadronic interaction models regarding the depth of maximum of air-shower profiles, X max , and ground-particle signals in water-Cherenkov detectors at 1000 m from the shower core, S ( 1000 ) , using the data from the fluorescence and surface detectors of the Pierre Auger Observatory. The test consists of fitting the measured two-dimensional ( S ( 1000 ) , X max ) distributions using templates for simulated air showers produced with hadronic interaction models pos-, et--04, 2.3d and leaving the scales of predicted X max and the signals from hadronic component at ground as free-fit parameters. The method relies on the assumption that the mass composition remains the same at all zenith angles, while the longitudinal shower development and attenuation of ground signal depend on the mass composition in a correlated way. The analysis was applied to 2239 events detected by both the fluorescence and surface detectors of the Pierre Auger Observatory with energies between 10 18.5 eV to 10 19.0 eV and zenith angles below 60°. We found, that within the assumptions of the method, the best description of the data is achieved if the predictions of the hadronic interaction models are shifted to deeper X max values and larger hadronic signals at all zenith angles. Given the magnitude of the shifts and the data sample size, the statistical significance of the improvement of data description using the modifications considered in the paper is larger than 5 σ even for any linear combination of experimental systematic uncertainties. Published by the American Physical Society 2024

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Machine-learning-enabled on-the-fly analysis of RHEED patterns during thin film deposition by molecular beam epitaxy

Thin film deposition is a fundamental technology for the discovery, optimization, and manufacturing of functional materials. Deposition by molecular beam epitaxy (MBE) typically employs reflection high-energy electron diffraction (RHEED) as a real-time in situ probe of the growing film. However, the state-of-the-art for RHEED analysis during deposition requires human observation. Here, we present an approach using machine learning (ML) methods to monitor, analyze, and interpret RHEED images on-the-fly during thin film deposition. In the analysis workflow, RHEED pattern images are collected at one frame per second and featurized using a pretrained deep convolutional neural network. The feature vectors are then statistically analyzed to identify changepoints; these changepoints can be related to changes in the deposition mode from initial film nucleation to a transition regime, smooth film deposition, and in some cases, an additional transition to a rough, islanded deposition regime. The feature vectors are additionally analyzed via graph analysis and community classification. The graph is quantified as a stabilization plot, and we show that inflection points in the stabilization plot correspond to changes in the growth regime. The full RHEED analysis workflow is termed RHAAPsody and includes data transfer and output to a visual dashboard. We demonstrate the functionality of RHAAPsody by analyzing the precaptured RHEED images from epitaxial depositions of anatase TiO2 on SrTiO3(001) and show that the analysis workflow can be executed in less than 1 s. Our approach shows promise as one component of ML-enabled real-time feedback control of the MBE deposition process.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Beyond Contrast Transfer: Spectral SNR as a Finite-Dose Metric for STEM Phase Retrieval

The contrast transfer function (CTF) is widely used to evaluate phase retrieval methods in scanning transmission electron microscopy (STEM), including center-of-mass imaging, parallax imaging, direct ptychography, and iterative ptychography. However, the CTF reflects only the maximum usable signal, neglecting the effects of finite electron fluence and the Poisson-limited nature of detection. As a result, it can significantly overestimate practical performance, especially in low-dose regimes. Here, we employ the spectral signal-to-noise ratio (SSNR), as a finite-dose statistical framework to evaluate the recoverable signal as a function of spatial frequency. Using numerical reconstructions of white-noise objects, we show that center-of-mass, parallax, and direct ptychography exhibit dose-independent SSNRs, with close-form analytic expressions. In contrast, iterative ptychography exhibits a surprising dose dependence: at low fluence, its SSNR converges to that of direct ptychography; at high fluence, it saturates at a value consistent with the maximum detective quantum efficiency predicted by recent quantum Fisher information bounds. The results highlight the limitations of CTF-based evaluation and motivate SSNR as a more accurate, finite-dose metric for assessing STEM phase retrieval methods.

STEM phase retrieval↗
DOE OSTI2024

Quantum Simulations of Radiation Damage in a Molecular Polyethylene Analog

Abstract An atomic‐level understanding of radiation‐induced damage in simple polymers like polyethylene is essential for determining how these chemical changes can alter the physical and mechanical properties of important technological materials such as plastics. Ensembles of quantum simulations of radiation damage in a polyethylene analog are performed using the Density Functional Tight Binding method to help bind its radiolysis and subsequent degradation as a function of radiation dose. Chemical degradation products are categorized with a graph theory approach, and occurrence rates of unsaturated carbon bond formation, crosslinking, cycle formation, chain scission reactions, and out‐gassing products are computed. Statistical correlations between product pairs show significant correlations between chain scission reactions, unsaturated carbon bond formation, and out‐gassing products, though these correlations decrease with increasing atom recoil energy. The results present relatively simple chemical descriptors as possible indications of network rearrangements in the middle range of excitation energies. Ultimately, the work provides a computational framework for determining the coupling between nonequilibrium chemistry in polymers and potential changes to macro‐scale properties that can aid in the interpretation of future radiation damage experiments on plastic materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Constitutive model development of aluminum alloy 1100 for elevated temperature forming process

Commercially pure aluminum alloy, AA1100, presents good electrical and thermal conductivity, high formability, and low cost. Those favorable characteristics have the potential to enable bipolar plates with improved economics and enhanced performance compared to current stainless steel bipolar plates for proton exchange membrane fuel cells. An accurate constitutive model is essential to develop and optimize processing parameters and effectively control the forming process. Here, the objective of this work is to develop a constitutive model of AA1100 that is able to simulate stress-strain relation, formed geometry, and predict the onset of fracture strain to avoid forming failure. Initially, a set of tensile tests at temperature between 300 and 500°C and strain rate between 0.005 and 1.0/s were conducted to examine the deformation behavior. Then, a set of damage-based unified visco-plastic constitutive equations is proposed and calibrated based on the results of stress-strain data. A genetic algorithm optimization method is applied to search for best fitting material constants in constitutive equations. The proposed model shows good predictability of both the stress-strain relation and fracture strain at low strain rate and high temperature conditions. The accuracy of proposed model is also evaluated statistically. A comparison of the proposed model with three popular models (Arrhenius-type mode, Johnson-Cook model and Zerilli-Armstrong model) was made. The proposed model shows the best experimental agreement with correlation coefficient of 0.96 in contrast to 0.25, 0.38 and 0.75 for the popular models, respectively. The proposed model can help to optimize the elevated temperature forming process and guide die design to enable optimal geometric features in the formed components.

08 HYDROGEN↗
DOE OSTI2025

Mapping strain and structural heterogeneities around bubbles in amorphous ionically conductive Bi 2 O 3

While amorphous materials are often approximated to have a statistically homogeneous atomic structure, they frequently exhibit localized structural heterogeneity that challenges simplified models. This study uses 4D scanning transmission electron microscopy to investigate the strain and structural modifications around gas bubbles in amorphous Bi 2 O 3 induced by argon irradiation. We present a method for determining strain fields surrounding bubbles that can be used to measure the internal pressure of the gas. Compressive strain is observed around the cavities, with higher-order crystalline symmetries emerging near the cavity interfaces, suggesting paracrystalline ordering as a result of bubble coarsening. This ordering, along with a compressive strain gradient, indicates that gas bubbles induce significant localized changes in atomic packing. By analyzing strain fields with maximum compressive strains of 3%, we estimate a lower bound on the internal pressure of the bubbles at 2.5 GPa. These findings provide insight into the complex structural behavior of amorphous materials under stress, particularly in systems with gas inclusions, and offer new methods for probing the local atomic structure in disordered materials. Although considering structural heterogeneity in amorphous systems is non-trivial, these features have crucial impacts on material functionalities, such as mechanical strength, ionic conductivity, and electronic mobility.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Sensitivity analysis of numerical modeling input parameters on floating offshore wind turbine loads in extreme idling conditions

Abstract. Floating offshore wind turbine (FOWT) systems are subject to complex environmental loads, with significant potential for damage in extreme storm conditions. Design simulations in these conditions are required to assess the survivability of the device with some level of confidence. Aero-hydro-servo-elastic engineering tools can be used with a reasonable balance of accuracy and computational efficiency. The models require many input parameters to describe the air and water conditions, the system properties, and the load calculations. Each of these parameters has some possible range, due to either statistical uncertainty or variations with time. Variation in the input parameters can have important effects on the uncertainty in the resulting loads, but it is not practical to perform detailed assessments of the impact of this uncertainty for every input parameter. This work demonstrates a method to identify the input parameters that have the most impact on the loads to focus further inspection. The process is done specifically for extreme storm load cases defined in the International Electrotechnical Commission design requirements for floating offshore wind turbines. The analysis was performed using the International Energy Agency Wind 15 MW offshore reference wind turbine atop the University of Maine VolturnUS-S reference platform in two US offshore wind regions, the Gulf of Maine and Humboldt Bay. It was found that the direction of incident waves and current, yaw misalignment, and the length of mooring line sections were among the primary sensitivities.

17 WIND ENERGY↗
DOE OSTI2024

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗
DOE OSTI2025

Statistical inference of collision frequencies from x-ray Thomson scattering spectra

Thomson scattering spectra measure the response of plasma particles to incident radiation. In warm dense matter, which is opaque to visible light, x-ray Thomson scattering (XRTS) enables a detailed probe of the electron distribution and has been used as a diagnostic for electron temperature, density, and plasma ionization. In this work, we examine the sensitivities of inelastic XRTS signatures to modeling details, including the dynamic collision frequency and the electronic density of states. Applying verified Monte Carlo inversion methods to dynamic structure factors obtained from time-dependent density functional theory, we assess the utility of XRTS signals as a way to inform the dynamic collision frequency, especially its direct-current limit, which is directly related to the electrical conductivity.

Collision frequency↗
DOE OSTI2026

Direct comparison of gamma, electron beam and X-ray radiation effects on the polymers of a pulsed lavage device

Gamma radiation used for sterilization of medical devices is challenged by cobalt-60 supply and commercial capacity. To maintain a robust radiation sterilization marketplace for the rapidly growing single-use medical device industry, investigation of potential alternatives to gamma technology, such as electron beam (e-beam) and X-ray technology, is critical. In this work, we directly compare the effects of radiation source and absorbed dose level on the polymeric materials and function of a commercial pulsed lavage device used for wound care. Product functionality, polymer mechanical, and polymer optical properties were evaluated using standard methods and input from the device manufacturer. Test results show that functionality of the product was not inhibited by radiation although the battery in the device exposed to X-ray exhibited greater voltage loss compared to batteries in products exposed to gamma or e-beam. Statistically significant differences between gamma and e-beam exposure and between gamma and X-ray exposure were also observed for product appearance in terms of yellowness index of several of the polymers considered. Overall, the results of this study support the viability of e-beam and X-ray radiation technologies as alternatives to cobalt-60 gamma technology for sterilization of the single-use pulsed lavage medical device investigated.

Electron beam↗
DOE OSTI2025

Block Lanczos algorithm for lattice QCD spectroscopy and matrix elements

Recent work introduced a new framework for analyzing correlation functions with improved convergence and signal-to-noise properties, as well as rigorous quantification of excited-state effects, based on the Lanczos algorithm and spurious eigenvalue filtering with the Cullum-Willoughby test. Here, we extend this framework to the analysis of correlation-function matrices built from multiple interpolating operators in lattice quantum chromodynamics (QCD) by constructing an oblique generalization of the block Lanczos algorithm, as well as a new physically motivated reformulation of the Cullum-Willoughby test that generalizes to block Lanczos straightforwardly. The resulting block Lanczos method directly extends generalized eigenvalue problem (GEVP) methods, which can be viewed as applying a single iteration of block Lanczos. Block Lanczos provides qualitative and quantitative advantages over GEVP methods analogous to the benefits of Lanczos over the standard effective mass, including faster convergence to ground- and excited-state energies, explicitly computable two-sided error bounds, straightforward extraction of matrix elements of external currents, and asymptotically constant signal-to-noise. No fits or statistical inference are required. Proof-of-principle calculations are performed for noiseless mock-data examples as well as two-by-two proton correlation-function matrices in lattice QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Type Ia Supernova Growth-rate Measurement with LSST Simulations: Intrinsic Scatter Systematics

Measurement of the growth rate of structures (fσ 8 ) with Type Ia supernovae (SNe Ia) will improve our understanding of the nature of dark energy and enable tests of general relativity. In this paper, we generate simulations of the 10 yr SN Ia data set of the Rubin-LSST survey, including a correlated velocity field from an N-body simulation and realistic models of SNe Ia properties and their correlations with host-galaxy properties. We find, similar to SN Ia analyses that constrain the dark energy equation-of-state parameters w 0 w a , that constraints on fσ 8 can be biased depending on the intrinsic scatter of SNe Ia. While for the majority of intrinsic scatter models we recover fσ 8 with a precision of ∼13%–14%, for the most realistic dust-based model, we find that the presence of non-Gaussianities in Hubble diagram residuals leads to a bias on fσ 8 of ∼ −20%. When trying to correct for the dust-based intrinsic scatter, we find that the propagation of the uncertainty on the model parameters does not significantly increase the error on fσ 8 . We also find that while the main component of the error budget of fσ 8 is the statistical uncertainty (>75% of the total error budget), the systematic error budget is dominated by the uncertainty on the damping parameter, σ u , that gives an empirical description of the effect of redshift space distortions on the velocity power spectrum. Our results motivate a search for new methods to correct for the non-Gaussian distribution of the Hubble diagram residuals, as well as an improved modeling of the damping parameter.

Carreres, Bastien [Duke Univ., Durham, NC (United ↗
DOE OSTI2024

Nuclear level density and γ -decay strength of Sr 93

This work presents the first experimentally determined nuclear level density and γ-ray strength function of the short-lived fission product 93 Sr, accomplished using the β-Oslo method. Direct measurement of the 92 Sr(n, γ) 93 Sr cross section is not currently possible, as the half-life of 2.66 hours is too short; instead, 93 Sr was formed through β decay of 93 Rb to excitation energies around the neutron separation energy. The γ-ray spectra were measured using a total absorption spectrometer at the National Superconducting Cyclotron Laboratory (NSCL) at Michigan State University (MSU). The statistical properties of the 93 Sr nucleus were experimentally determined, including the γ-ray strength function and nuclear level density. At low energies, the γ-ray strength function exhibits a constant γ-decay strength, rather than a slightly increasing strength with decreasing γ-ray energy as had been previously observed for several nuclei in this mid-mass region. Finally, these statistical properties were then implemented in the reaction code TALYS1.95 to calculate the 92 Sr(n, γ) 93 Sr cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Towards revealing intrinsic vortex-core states in Fe-based superconductors through statistical discovery

Abstract In type-II superconductors, electronic states within magnetic vortices hold crucial information about the paring mechanism and can reveal non-trivial topology. While scanning tunneling microscopy/spectroscopy (STM/S) is a powerful tool for imaging superconducting vortices, it is challenging to isolate the intrinsic electronic properties from extrinsic effects like subsurface defects and disorders. Here we combine STM/STS with basic machine learning to develop a method for screening out the vortices pinned by embedded disorder in iron-based superconductors. Through a principal component analysis of large STS data within vortices, we find that the vortex-core states in Ba(Fe 0.96 Ni 0.04 ) 2 As 2 start to split into two categories at certain magnetic field strengths, reflecting vortices with and without pinning by subsurface defects or disorders. Our machine-learning analysis provides an unbiased approach to reveal intrinsic vortex-core states in novel superconductors and shed light on ongoing puzzles in the possible emergence of a Majorana zero mode.

Guo, Yueming↗
DOE OSTI2026

Measurement of low energy characteristic X-rays emitted present in U-10Mo HALEU alloy fuel

Advance high density monolithic alloy HALEU fuels have been scanned with a silicon drift detector scanning system to measure low-energy characteristic X-rays emitted by uranium daughters, as well as zirconium present as a foil diffusion barrier. Four peaks were measured across the surfaces of three plates, and are consistent with other thorium and protactinium peaks present in depleted uranium samples that are induced by the decay of 235 U and 238 U. Of the four peaks detected, three are dominated by the decay of 235 U, and one peak at 15.70 keV covers two characteristic X-rays from the zirconium diffusion barrier (Kα 1 ) and from a 231 Th 15.62 keV X-ray induced by 235 U α-decay. The low energies of these lines suggest that they are mostly emitted from the fuel meat surface, and the poor counting statistics caused by strong self-attenuation make these data a poor candidate for 235 U/ 238 U ratio analysis. We suggest, then, that these X-rays could be used for passive measurement of variations in the cladding thickness of these fuels to complement manual thickness measurement methods.

HALEU↗
DOE OSTI2026

Generation of random geological models using multi-randomization for machine learning

Generating high-fidelity geological models is essential for advancing machine learning (ML) methods in automated seismic interpretation. For instance, seismic images paired with corresponding fault labels are foundational for ML-based fault detection from seismic migration sections. While several open-access datasets of random geological models exist, open-source tools specifically designed to produce large volumes of such models for ML applications remain scarce. To address this gap, we present RGM (Random Geological Model), an open-source software package for efficiently generating 2D and 3D synthetic geological models tailored for ML workflows. RGM supports the creation of diverse model components, including medium property distributions (P-/S-wave velocities and density), seismic reflectivity images (i.e., synthetic migration sections), relative geological time, and discrete fault attributes such as probability, dip, strike, rake, and displacement. It also accommodates the creation of complex geological features such as salt bodies and unconformities. The model generation algorithm employs a multi-randomization strategy, yielding an effectively infinite-dimensional model space that encompasses a wide range of geological scenarios and associated seismic features. Furthermore, RGM incorporates a method to generate synthetic elastic migration images using analytical elastic reflection coefficients combined with frequency-dependent scaling. This functionality enables the creation of training datasets for ML models that leverage elastic seismic images. RGM is implemented in modern object-oriented Fortran, allowing users to flexibly control statistical parameters governing model variability. We demonstrate the capability, performance, and geological realism of the package through comprehensive 2D and 3D examples.

58 GEOSCIENCES↗