Search NASA⌕ Search

SEARCH · Search NASA

Results for “Statistical accuracy”

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 19 records

Combining machine-learned and empirical force fields with the parareal algorithm: application to the diffusion of atomistic defects

We numerically investigate an adaptive version of the parareal algorithm in the context of molecular dynamics. This adaptive variant has been originally introduced in [1]. We focus here on test cases of physical interest where the dynamics of the system is modelled by the Langevin equation and is simulated using the molecular dynamics software LAMMPS. In this work, the parareal algorithm uses a family of machine-learning spectral neighbor analysis potentials (SNAP) as fine, reference, potentials and embedded-atom method potentials (EAM) as coarse potentials. We consider a self-interstitial atom in a tungsten lattice and compute the average residence time of the system in metastable states. Our numerical results demonstrate significant computational gains using the adaptive parareal algorithm in comparison to a sequential integration of the Langevin dynamics. We also identify a large regime of numerical parameters for which statistical accuracy is reached without being a consequence of trajectorial accuracy.

36 MATERIALS SCIENCE↗

Robust measurement of microbial reduction of graphene oxide nanoparticles using image analysis

ABSTRACT Shewanella oneidensis ( S. oneidensis ) has the capacity to reduce electron acceptors within a medium and is thus used frequently in microbial fuel generation, pollutant breakdown, and nanoparticle fabrication. Microbial fuel setups, however, often require costly or labor-intensive components, thus making optimization of their performance onerous. For rapid optimization of setup conditions, a model reduction assay can be employed to allow simultaneous, large-scale experiments at lower cost and effort. Since S. oneidensis uses different extracellular electron transfer pathways depending on the electron acceptor, it is essential to use a reduction assay that mirrors the pathways employed in the microbial fuel system. For microbial fuel setups that use nanoparticles to stimulate electron transfer, reduction of graphene oxide provides a more accurate model than other commonly used assays as it is a bulk material that forms flocculates in solutions with a large ionic component. However, graphene oxide flocculates can interfere with traditional absorbance-based measurement techniques. This study introduces a novel image analysis method for quantifying graphene oxide reduction, showing improved performance and statistical accuracy over traditional methods. A comparative analysis shows that the image analysis method produces smaller errors between replicates and reveals more statistically significant differences between samples than traditional plate reader measurements under conditions causing graphene oxide flocculation. Image analysis can also detect reduction activity at earlier time points due to its use of larger solution volumes, enhancing color detection. These improvements in accuracy make image analysis a promising method for optimizing microbial fuel cells that use nanoparticles or bulk substrates. IMPORTANCE Shewanella oneidensis ( S. oneidensis ) is widely used in reduction processes such as microbial fuel generation due to its capacity to reduce electron acceptors. Often, these setups are labor-intensive to operate and require days to produce results, so use of a model assay would reduce the time and expenses needed for optimization. Our research developed a novel digital analysis method for analysis of graphene oxide flocculates that may be utilized as a model assay for reduction platforms featuring nanoparticles. Use of this model reduction assay will enable rapid optimization and drive improvements in the microbial fuel generation sector.

Bennett, Danielle T. (ORCID:0009000188748827)↗

Observation of Cosmic-Ray Anisotropy in the Southern Hemisphere with 12 yr of Data Collected by the IceCube Neutrino Observatory

We analyzed the 7.92 × 10 11 cosmic-ray-induced muon events collected by the IceCube Neutrino Observatory from 2011 May 13, when the fully constructed experiment started to take data, to 2023 May 12. This data set provides an up-to-date cosmic-ray arrival direction distribution in the Southern Hemisphere with unprecedented statistical accuracy covering more than a full period length of a solar cycle. Improvements in Monte Carlo event simulation and better handling of year-to-year differences in data processing significantly reduce systematic uncertainties below the level of statistical fluctuations compared to the previously published results. We confirm the observation of a change in the angular structure of the cosmic-ray anisotropy between 10 TeV and 1 PeV, more specifically in the 100–300 TeV energy range. For the first time, we analyzed the angular power spectrum at different energies. The observed variations of the power spectra with energy suggest relatively reduced large-scale features at high energy compared to those of medium and small scales. The large volume of data enhances the statistical significance at higher energies, up to the PeV scale, and smaller angular scales, down to approximately 6° compared to previous findings.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of the Challenges in Developing Sample-Based Multi-fidelity Estimators for Non-deterministic Models

Multifidelity (MF) uncertainty quantification (UQ) seeks to leverage and fuse information from a collection of models to achieve greater statistical accuracy with respect to a single-fidelity counterpart, while maintaining an efficient use of computational resources. Despite many recent advancements in MF UQ, several challenges remain and these often limit its practical impact in certain application areas. In this manuscript, we focus on the challenges introduced by nondeterministic models to sampling MF UQ estimators. Nondeterministic models produce different responses for the same inputs, which means their outputs are effectively noisy. MF UQ is complicated by this noise since many state-of-the-art approaches rely on statistics, e.g., the correlation among models, to optimally fuse information and allocate computational resources. Here, we demonstrate how the statistics of the quantities of interest, which impact the design, effectiveness, and use of existing MF UQ techniques, change as functions of the noise. With this in hand, we extend the unifying approximate control variate framework to account for nondeterminism, providing for the first time a rigorous means of comparing the effect of nondeterminism on different multifidelity estimators and analyzing their performance with respect to one another. Numerical examples are presented throughout the manuscript to illustrate and discuss the consequences of the presented theoretical results.

97 MATHEMATICS AND COMPUTING↗

PROTOCALC, a W -band Polarized Calibrator for Cosmic Microwave Background Telescopes: Application to Simons Observatory and CLASS

Current- and next-generation cosmic microwave background (CMB) experiments will measure polarization anisotropies with unprecedented sensitivities. The need for high precision in these measurements underscores the importance of gaining a comprehensive understanding of instrument properties, with a particular emphasis on the study of the beam properties, and especially their polarization characteristics and the measurement of the polarization angle. In this context, a major challenge lies in the scarcity of millimeter polarized astrophysical sources with sufficient brightness and calibration knowledge to meet the stringent accuracy requirements of future CMB missions. This led to the development of a drone-borne calibration source designed for the frequency band centered on approximately 90 GHz, matching a commonly used channel in ground-based CMB measurements. The Prototype Calibrator for Cosmology, PROTOCALC, has undergone thorough in-lab testing, and its properties have been subsequently modeled through simulation software integrated into the standard Simons Observatory analysis pipeline. Moreover, the PROTOCALC system has been tested in the field, having been deployed twice on calibration campaigns with CMB telescopes in the Atacama Desert. The data collected constrain the roll angle of the source with a statistical accuracy of 0$^°_•$045.

79 ASTRONOMY AND ASTROPHYSICS↗

Hierarchical Gaussian Random Field Sampling for Multilevel Markov Chain Monte Carlo: Coupling Stochastic Partial Differential Equation and the Karhunen–Loève Decomposition

This work introduces structure preserving hierarchical decompositions for sampling Gaussian random fields (GRFs) within the context of multilevel Bayesian inference in high-dimensional space. Existing scalable hierarchical sampling methods, such as those based on stochastic partial differential equations (SPDEs), often reduce the dimensionality of the sample space at the cost of accuracy of inference. Other approaches, such that those based on Karhunen-Loève (KL) expansions, offer sample space dimensionality reduction but sacrifice GRF representation accuracy and ergodicity of the Markov chain Monte Carlo (MCMC) sampler and are computationally expensive for high-dimensional problems. The proposed method integrates the dimensionality reduction capabilities of KL expansions with the scalability of SPDE-based sampling, thereby providing a robust, unified framework for high-dimensional uncertainty quantification (UQ) that is scalable and accurate, preserves ergodicity, and offers dimensionality reduction of the sample space. The hierarchy in our multilevel algorithm is derived from the geometric multigrid hierarchy. By constructing a hierarchical decomposition that maintains the covariance structure across the levels in the hierarchy, the approach enables efficient coarse-to-fine sampling while ensuring that all samples are drawn from the desired distribution. The effectiveness of the proposed method is demonstrated on a benchmark subsurface flow problem, demonstrating its effectiveness in improving computational efficiency and statistical accuracy. Furthermore, our proposed technique is more efficient and accurate and displays better convergence properties than existing methods for high-dimensional Bayesian inference problems.

Gaussian random fields↗

Reconstructing neutrinoless double beta decay event kinematics in a xenon gas detector with vertex tagging

If neutrinoless double beta decay is discovered, the next natural step would be understanding the lepton number violating physics responsible for it. Several alternatives exist beyond the exchange of light neutrinos. Some of these mechanisms can be distinguished by measuring phase-space observables, namely the opening angle cos θ among the two decay electrons, and the electron energy spectra, T 1 and T 2 . In this work, we study the statistical accuracy and precision in measuring these kinematic observables in a future xenon gas detector with the added capability to precisely locate the decay vertex. For realistic detector conditions (a gas pressure of 10 bar and spatial resolution of 4 mm), we find that the average $\overline{cos θ}$ and $\overline{T_1}$ values can be reconstructed with a precision of 0.19 and 110 keV, respectively, assuming that only 10 neutrinoless double beta decay events are detected.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Determining reference standard strength for neutron-irradiated reduced activation ferritic/martensitic steel F82H by Bayesian method

The deterministic approach widely adopted in the design of structural components relies on systematically defined design limits using empirically determined safety factors. However, this approach is not always appropriate because structures are subjected to a variety of loads in the practical environment, which may result in excessively conservative design limits. In recent years, a more rigorous probabilistic approach that incorporates material strength distributions has become an important solution. In the probabilistic approach, the probability density functions of material strength properties underpin the design criteria. Here, the objective of this study is to identify the density distribution functions that best describe tensile properties of irradiated F82H to define a reference strength for DEMO design. Due to the limited number of existing data, this study specifically employs a Bayesian prediction method based on Monte Carlo simulations to determine a material reference value with statistical reliability and to investigate its effectiveness. For example, the dependence of tensile properties of 300 °C irradiated materials on irradiation damage and the range predicted by 95% Bayesian estimation was evaluated. As a statistical model for the dose dependence of statistical parameters, the normal distribution exhibited a better fit for 0.2% proof strength and tensile strength, whereas the distribution of total elongation data gave comparable reference values for both the normal and Weibull distribution models. Both models gave comparable criteria for the distribution of total elongation data. The Weibull model also gave better results for uniform elongation. The function best describing the model was a logarithmic law for both 0.2% proof strength and tensile strength, while a power law for both total and uniform elongation, which allowed for more comprehensive data prediction of irradiation data with statistical accuracy for DEMO reactor design.

36 MATERIALS SCIENCE↗

Measurements of short-lived fission product yields from photofission of 238 U using 13.0 MeV monoenergetic photons

Photon-induced fission product yield (FPY) measurements were conducted on the isotope 238U. Fission was induced using Eγ = 13.0 MeV monoenergetic photons produced by the Triangle Universities Nuclear Laboratory’s (TUNL’s) High Intensity γ-ray Source (HIγS) facility. Short-lived FPYs were measured by performing cyclic activation of the sample using a rapid target transfer system. Following activation, the 238U target was rapidly (0.4 s) transferred to a counting station consisting of two well-shielded high-purity germanium (HPGe) detectors. The irradiation-counting cycle was repeated until the summed data had sufficient statistical accuracy. Twenty-eight unique fission products with half-lives ranging from 1 s to 450 s were identified, and their cumulative FPYs determined. Furthermore, the results are compared with previous independent FPY measurements using inverse kinematics. Good agreement between the data sets is found despite the different excitation energy distributions of the fissioning nucleus in the experiments.

Physics - Nuclear physics and radiation physics↗

Final SeaQuest Results on the Flavor Asymmetry of the Proton Light-Quark Sea with Proton-Induced Drell-Yan Process

The Fermilab E906/SeaQuest collaboration performed measurements of the Drell-Yan process using 120 GeV proton beams bombarding liquid hydrogen and liquid deuterium targets. A combined analysis of all collected data was performed to obtain the final results for the 𝜎 𝑝⁢𝑑 /2⁢𝜎 𝑝⁢𝑝 Drell-Yan cross section ratio covering the kinematic region of 0.13 < 𝑥 <0.45. The 𝑥 dependencies of $\bar{d}$(𝑥)/$\bar{𝑢}$⁡(𝑥) and $\bar{d}$(𝑥) $-$ $\bar{𝑢}$⁡(𝑥) are extracted from these cross section ratios. It is found that $\bar{d}$(𝑥) is greater than $\bar{𝑢}$⁡(𝑥) over the entire measured 𝑥 range, with improved statistical accuracy compared to previous measurements. The new results on $\bar{d}$(𝑥)/$\bar{𝑢}$⁡(𝑥) and $\bar{d}$(𝑥) $-$ $\bar{𝑢}$⁡(𝑥) are compared to various parton distribution functions and theoretical calculations.

Leung, C. H. [University of Illinois at Urbana-Cha↗

Empirical radius formulas for canonical neutron stars from bidirectionally selecting features of equations of state in extended Bayesian analyses of observational data

Significant advancement in Bayesian inference of nuclear equation of state (EOS) from gravitational wave and x-ray observations of neutron stars (NSs) has been made by the nuclear astrophysics community especially since GW170817. By extending the traditional Bayesian analysis which normally ends at presenting the marginalized posterior probability distribution functions (PDFs) of individual EOS parameters and their correlations (or sometimes only the Pearson correlation coefficients which are only reliably useful when the variables are linearly correlated while they are actually often not), we search for a data-driven and robust empirical formula for the radius 𝑅 1.4 of canonical NSs in terms of the characteristic EOS parameters (features). We also identify the single most important but currently poorly known EOS parameter for determining the 𝑅 1.4 . Using three regression-model-building methodologies: bidirectional stepwise feature selection, least absolute shrinkage selection operator (LASSO) regression, and neural network regression on a large set of posterior EOSs and the corresponding 𝑅 1.4 values inferred from earlier comprehensive Bayesian analyses of NS observational data, we systematically and rigorously develop the most probable 𝑅 1.4 formulas with varying statistical accuracy and technical complexity. Here, the most important EOS parameters for determining 𝑅 1.4 are found consistently in each of the feature selection processes to be (in order of decreasing importance): curvature 𝐾 sym , slope 𝐿, skewness 𝐽 sym of nuclear symmetry energy, skewness 𝐽 0 , incompressibility 𝐾 0 of symmetric nuclear matter, and the magnitude 𝐸 sym ⁡(𝜌 0 ) of symmetry energy at the saturation density 𝜌 0 of nuclear matter.

Bayesian methods↗

Beyond PCA: Additional Dimension Reduction Techniques to Consider in the Development of Climate Fingerprints

Abstract Dimension reduction techniques are an essential part of the climate analyst’s toolkit. Due to the enormous scale of climate data, dimension reduction methods are used to identify major patterns of variability within climate dynamics, to create compelling and informative visualizations, and to quantify major named modes such as El Niño–Southern Oscillation. Principal components analysis (PCA), also known as the method of empirical orthogonal functions (EOFs), is the most commonly used form of dimension reduction, characterized by a remarkable confluence of attractive mathematical, statistical, and computational properties. Despite its ubiquity, PCA suffers from several difficulties relevant to climate science: high computational burden with large datasets, decreased statistical accuracy in high dimensions, and difficulties comparing across multiple datasets. In this paper, we introduce several variants of PCA that are likely to be of use in climate sciences and address these problems. Specifically, we introduce non-negative , sparse , and tensor PCA and demonstrate how each approach provides superior pattern recognition in climate data. We also discuss approaches to comparing PCA-family results within and across datasets in a domain-relevant manner. We demonstrate these approaches through an analysis of several runs of the E3SM climate model from 1991 to 1995, focusing on the simulated response to the Mt. Pinatubo eruption; our findings are consistent with a recently identified stratospheric warming fingerprint associated with this type of stratospheric aerosol injection.

Weylandt, Michael↗

Measuring Neutron Polarisation in Deuteron Photo-disintegration with the CLAS Start Counter [Thesis]

Deuteron photo-disintegration (γd → γp) is a reaction that represents the simplest case in which nuclear and hadron physics models can be tested. Despite this, associated polarization analyses are limited in terms of angular coverage and energy ranges, especially in observables related to the recoil neutron. This is largely due to a lack in dedicated polarimetry equipment, and represents a roadblock in global progress to understand high-energy phenomena such as hexaquarks, and quark-gluon degrees of freedom. To address this problem, this PhD thesis pioneers a new methodology for the parasitic measurement of nucleon polarization using kinematic reconstruction of (spin-dependent) nucleon-nucleus scattering of reaction products, prior to their detection in large acceptance particle detector apparatus. Following this novel approach, which requires no dedicated polarimeter, a determination of the double polarization observable, $C^n_{x'}$, from deuteron photo-disintegration is presented, using Jefferson Lab’s CLAS detector. The analysis utilizes the (n,p) charge exchange reaction in CLAS’s "start counter" (plastic scintillator) to determine the final state neutron polarizations. The results present the first ever data for this observable above 0.7 GeV (photon beam energy) and significantly extend the angular range of the world data set. This new data is largely statistically consistent with the previous measurement of $C^n_{x'}$ by Bashkanov et al . in the overlapping energy range of 0.4-0.7 GeV. It is planned for the statistical accuracy of the presented result to be increased by the inclusion of additional data. The analysis herein serves as a key proof of concept for future applications, including a recommended similar analysis to be implemented with data from the more modern CLAS12 detector. This paves the way for a plethora of additional analyses using existing data sets that would provide crucial new constraints for hadron and nuclear physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Transition from Vehicular to Structural Ionic Transport in Electrified Alkali Aqueous Solutions

A molecular understanding of the solvation and dynamics of ions under static electric fields is crucial for modelling a wide range of natural and technological processes. Yet, traditional simulation methods suffer from a trade-off that has to be made between accuracy and statistical convergence. To bridge this gap, herein we extend our recently introduced Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) approach to investigate the solvation structures and ionic transport mechanisms of electrified alkali cationic solutions. We obtain ionic conductivities for Li+, Na+ and Cs+ from the field dependence of the ionic current density in good agreement with experiment. Surprisingly, the migration mechanism is found to be strikingly different for the three ions despite their similar ionic conductivities. While Li+ conducts predominantly through vehicular migration of a stable 4-fold coordinated ion at all field strengths, Cs+ conducts strictly through a structural diffusion mechanism, where 9-12 transient first shell water coordination bonds are broken and renewed. Notably, aqueous Na+ emerges as a “Goldilocks” ion: its ion-water interactions are strong enough to maintain distinct 5-6 fold coordination shells at zero field (unlike Cs+) yet labile enough to be strongly perturbed by electric fields (unlike Li+). As a consequence, we observe an electric field-induced transition from vehicular to structural ionic transport for Na+ that is accompanied by a marked increase in the ionic current density. Our results imply that the conductance mechanism of ions with moderate ion-solvent interactions can be effectively tuned by external electric fields.

Joll, Kit↗

Quantitative Infrared-to-Terahertz Nanospectroscopy of Semiconductors

Semiconductor technology now employs few-nanometer features, necessitating tools probing electronic properties on the same length scale. While the concentration of free charge carriers is routinely measured, the scattering rate remains challenging to access at the nanoscale. Here, we present ultrabroadband (5–50 THz) synchrotron infrared nanospectroscopy as a quantitative metrology tool for semiconductors. This technique can determine both the charge carrier concentration and scattering rate with percent-level accuracy, and it is inherently capable of ∼10 nm spatial resolution. We study silicon with different doping levels and confirm the method’s accuracy by statistical analysis and comparison with established far-field infrared spectroscopy. Near-field measurements systematically reveal charge-carrier concentrations ∼30% lower than far-field values, consistent with increased surface sensitivity and surface depletion. Our work establishes synchrotron infrared nanospectroscopy as a precise tool for quantitative nanoscale semiconductor characterization and paves the way toward all-optical characterization of surface depletion effects.

36 MATERIALS SCIENCE↗

Measurement of neutron production in atmospheric neutrino interactions at Super-Kamiokande

We present measurements of total neutron production from atmospheric neutrino interactions in water, analyzed as a function of electron-equivalent visible energy over a range of 30 MeV to 10 GeV. These results are based on 4,270 days of data collected by Super-Kamiokande, including 564 days with 0.011 wt% gadolinium added to enhance neutron detection. Neutron signal selection is based on a neural network trained on simulation, with its performance validated using an Am/Be neutron point source. The measurements are compared to predictions from neutrino event generators combined with various hadron-nucleus interaction models, which include an intranuclear cascade model and a nuclear deexcitation model. We observe significant variations in the predictions depending on the choice of hadron-nucleus interaction model. We discuss key factors that contribute to describing our data, such as in-medium effects in the intranuclear cascade and the accuracy of statistical evaporation modeling.

Cherenkov detectors↗

Elimination LArTPC Simulation Uncertainty

Liquid Argon Time Projection Chambers (LArTPC) are crucial for measuring muons and neutrinos by capturing the paths of fast-moving particles through argon gas. However, these detectors face challenges such as electron-ion recombination, diffusion, and attenuation, which introduce uncertainties in simulation models. This study, conducted by Ka ren Mkrtchyan at FERMILAB, aims to reduce these uncertainties by adjusting the amplitude and width of signals detected by the TPC wires. Initial findings indicate that the current modification algorithm requires further refinement to better align simulations with observed data. Ongoing work focuses on correcting computational bugs and enhancing the simulation model for improved accuracy and statistical confidence.

Mkrtchyan, Ka'ren↗

Trust Your Gut: Comparing Human and Machine Inference from Noisy Visualizations

People commonly utilize visualizations not only to examine a given dataset, but also to draw generalizable conclusions about the underlying models or phenomena. Prior research has compared human visual inference to that of an optimal Bayesian agent, with deviations from rational analysis viewed as problematic. However, human reliance on non-normative heuristics may prove advantageous in certain circumstances. We investigate scenarios where human intuition might surpass idealized statistical rationality. In two experiments, we examine individuals’ accuracy in characterizing the parameters of known data-generating models from bivariate visualizations. Our findings indicate that, although participants generally exhibited lower accuracy compared to statistical models, they frequently outperformed Bayesian agents, particularly when faced with extreme samples. Participants appeared to rely on their internal models to filter out noisy visualizations, thus improving their resilience against spurious data. However, participants displayed overconfidence and struggled with uncertainty estimation. They also exhibited higher variance than statistical machines. Our findings suggest that analyst gut reactions to visualizations may provide an advantage, even when departing from rationality. These results carry implications for designing visual analytics tools, offering new perspectives on how to integrate statistical models and analyst intuition for improved inference and decision-making. The data and materials for this paper are available at https://osf.io/qmfv6

human-machine collaboration↗