Search NASA⌕ Search

SEARCH · Search NASA

Results for “error bars”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Putting error bars on density functional theory dataset

This dataset contains submission files and raw output files from high-throughput DFT simulations to analyze the systemic errors in lattice constant, bulk moduli and formation energy predictions for a range of binary and ternary oxides using four exchange correlation functionals (LDA, PBE, PBEsol and vdW-DF-C09). This data was then used as the basis for employing materials informatics methods to predict the expected errors in the lattice constants of the studied compounds. Predicted errors were also used to better the DFT-predicted lattice parameters. Our results emphasize the link between the computed errors and the electron density and hybridization errors of a functional. In essence, these results provide “error bars” for choosing a functional for the creation of high-accuracy, high-throughput datasets as well as avenues for the development of XC functionals with enhanced performance, thereby enabling the accelerated discovery and design of new materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Measuring and unbiasing the BAO shift in the Ly α forest with AbacusSummit

ABSTRACT The Dark Energy Spectroscopic Instrument (DESI) places sub- per cent constraints on measurements of the Baryon Acoustic Oscillation (BAO) scaling parameters from the Ly $\alpha$ forest. However, no systematic error budget stemming from non-linearities in the three-dimensional clustering of the Ly $\alpha$ forest is included in the DESI-Ly $\alpha$ analysis. In this work, we measure the size of the shift of the BAO peak using large Ly $\alpha$ forest mocks produced on the N-body simulation suite AbacusSummit, which adopt the Fluctuating–Gunn–Peterson Approximation (FGPA). Specifically, we measure the Ly $\alpha$ autocorrelation and the Ly $\alpha$-quasar cross-correlation functions. To mitigate the noise, we adopt a linear control variates technique, reducing the error bars by a factor of up to $\sim \sqrt{50}$ on large scales. From the autocorrelation, we detect a small positive shift in radial direction of $\Delta \alpha _{\parallel }= 0.35~{{\ \rm per\ cent}}$ at the 3$\sigma$ level and virtually no shift in the transverse direction, $\alpha _\perp$. From the cross-correlation, we see a similar shift to $\Delta \alpha _\parallel$, albeit with larger error bars, and a small negative shift, $\Delta \alpha _{\perp }=\sim$0.25 per cent, at the 2$\sigma$ level. We also make a connection with the Ly $\alpha$ forest effective field theory (EFT) framework and find that the one-loop EFT power spectrum yields unbiased measurements of the BAO shift parameters in radial and transverse direction for Ly $\alpha$ auto- and the Ly $\alpha$-quasar cross-correlation measurements. When using the one-loop EFT framework, we find that we can recover the BAO parameters without a shift, which has important implications for future Ly $\alpha$ forest analyses based on EFT. This work paves the way for novel full-shape analyses of the currently observing DESI and future surveys such as the PFS, WEAVE-QSO, and 4MOST.

Hadzhiyska, Boryana↗
DOE OSTI2026

Predicting Li-Ion Battery Capacity Fade Using Early-Life Data and a Hybrid Data-Driven Gaussian Process-Bayesian Regression Approach

Accurately predicting Li-ion battery capacity trajectories using early-life data can dramatically improve battery-life understandings and be used to rapidly evaluate design/cost/performance trade-offs when developing new battery materials. Accurate early-life predictions enable researchers to quickly iterate over cell designs and material precursor properties without consistently cycling cells to failure. To this end, we present a toolbox that uses a combined Gaussian Process and Bayesian regression approach that capitalizes on signals other than just capacity (e.g., dQ/dV, voltage drops) to rapidly predict capacity-fade trajectories. The prediction tool uses Bayesian regression to fit functional forms, e.g., power law, sigmoids, etc., to predict capacity-fade dynamics. By fitting functional forms, the capacity fade can be interrogated at any point in the future, allowing for early cell-failure prediction. Additionally, Bayesian regression allows for accurate uncertainty estimates that account for cell-to-cell variability (aleatoric uncertainty) and the lack of observation data (epistemic uncertainty). By only using early cycle data to predict the capacity fade trajectory, uncertainty bounds at end-of-life can be extremely large. The large uncertainty bounds are further exacerbated because there is no systematic way to define the prior distribution of the functional forms' parameters. We improve our the predicted trajectory confidence interval of our predicted trajectory using two methods. First, we shows that a small amount of held-out cycling data is sufficientuse some train cells, that have been cycled to failure to derive information regarding the appropriate prior distributions for the functional forms' parameters of the functional form, effectively leading to data-driven priors.. We propose constructing the data-driven priors by first running a Bayesian regression starting with uninformed priors to generate intermediate cell-specific posterior parameter distributions. These posterior distributions are combined using a Ggaussian mixture model for each parameter to create the data-driven priors. These mixture models serve as the data-driven prior distributions for the parameters for. Second, we derive multiple features, e.g., C_dchg 0.5 DoD 0.5, log (|mean(dQ/dV_(w_3-w_0 ) (V)|), etc., from the train cellsheld-out cycling data, identify which the features are that best predicting capacity at early/mid-life cycles, and then create Ggaussian process regression models that are used for predicting capacity at early/mid-life cycles for the test cells (see blue dots with error bars in Fig 1b). Finally, these predicted data-points are used in addition to the actual early cycle data capacity fade to construct the Bayesian regression trajectory for the test cell s. Notably. We note that these two methods are complementary and can be combined with each other. We evaluate the performance of our proposed method on an testing open-source dataset from Iowa State University and Iowa Lakes Community College (ISU-ILCC). This dataset comprises of 251 nickel-manganese-cobalt/graphite Lithium-ion cells that are cycled under 63 different conditions. We compute the mean average percentage error (MAPE) and negative log predictive density (NLPD) to quantify the efficacy of our method. Our initial findings suggest that, when only few observations are available, for test cells, when using only Bayesian regression with uninformed priors, a power law functional provides the most accurate predictions. with very few data points. However, asHowever, a the number of data points increases, a twin sigmoidal function becomes more accurate as the number of observations further increases. We also find that using as little as 10% of the data set towards generating data-driven priors can lead to significant improvement in prediction accuracy when using early cycle data. Lastly, we found that augmenting early-cycle data with Gaussian process-predicted capacity data for Bayesian regression greatly improves the prediction accuracy. We will present a comprehensive comparison of our methods to other methods available in the literature and apply this method to additional battery datasets.

42 ENGINEERING↗
DOE OSTI2024

Ensemble Monte Carlo calculations with five novel moves

We introduce five novel types of Monte Carlo (MC) moves that brings the number of moves of ensemble MC calculations from three to eight. So far such calculations have relied on affine invariant stretch moves that were originally introduced by Christen (2007), walk moves by Goodman and Weare (2010) and quadratic moves by Militzer (2023). Ensemble MC methods have been very popular because they harness information about the fitness landscape from a population of walkers rather than relying on expert knowledge. Here we modified the affine method and employed a simplex of points to set the stretch direction. We adopt the simplex concept to quadratic moves. We also generalize quadratic moves to arbitrary order. Finally, we introduce directed moves that employ the values of the probability density while all other types of moves rely solely on the location of the walkers. We apply all algorithms to the Rosenbrock density in 2 and 20 dimensions and to the ring potential in 12 and 24 dimensions. We evaluate their efficiency by comparing error bars, autocorrelation time, travel time, and the level of cohesion that measures whether any walkers were left behind. Our code is open source.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Generalized framework for likelihood-based field-level inference of growth rate from velocity and density fields

Measuring the growth rate of large-scale structures ( f ) as a function of redshift has the potential to break degeneracies between modified gravity and dark energy models, when combined with expansion-rate probes. Direct estimates of peculiar velocities of galaxies have attracted interest as a means of estimating fσ 8 . In particular, field-level methods can be used to fit the field nuisance parameter along with cosmological parameters simultaneously. This article aims to provide the community with a unified framework for the theoretical modeling of the likelihood-based field-level inference by performing fast field covariance calculations for velocity and density fields. Our purpose is to lay the foundations for a nonlinear extension of the likelihood-based method at the field level. We have developed a generalized framework, implemented in the dedicated software flip to perform a likelihood-based inference of fσ 8 . We derived a new field covariance model, which includes wide-angle corrections. We also included the models previously described in the literature inside our framework. We compared their performance against ours, and we validated our model by comparing it with the two-point statistics of a recent N-body simulation. The tests we performed have allowed us to validate our software and determine the appropriate wavenumber range to integrate our covariance model and its validity in terms of separation. Our framework allows for a wider wavenumber coverage to be used in our calculations than in previous works, which is particularly interesting for nonlinear model extensions. Finally, our generalized framework allows us to efficiently perform a survey geometry-dependent Fisher forecast of the fσ 8 parameter. We show that the Fisher forecast method we developed gives an error bar that is 30% closer to a full likelihood-based estimation than a standard volume Fisher forecast.

Ravoux, Corentin↗
DOE OSTI2024

X-ray imaging and electron temperature evolution in laser-driven magnetic reconnection experiments at the national ignition facility

We present results from x-ray imaging of high-aspect-ratio magnetic reconnection experiments driven at the National Ignition Facility. Two parallel, self-magnetized, elongated laser-driven plumes are produced by tiling 40 laser beams. A magnetic reconnection layer is formed by the collision of the plumes. A gated x-ray framing pinhole camera with micro-channel plate detector produces multiple images through various filters of the formation and evolution of both the plumes and current sheet. As the diagnostic integrates plasma self-emission along the line of sight, two-dimensional electron temperature maps ⟨Te⟩Y are constructed by taking the ratio of intensity of these images obtained with different filters. The plumes have a characteristic temperature ⟨Te⟩Y=240 ± 20 eV at 2 ns after the initial laser irradiation and exhibit a slow cooling up to 4 ns. The reconnection layer forms at 3 ns with a temperature ⟨Te⟩Y=280 ± 50 eV as the result of the collision of the plumes. The error bars of the plumes and current sheet temperatures separate at 4 ns, showing the heating of the current sheet from colder inflows. Using a semi-analytical model, we survey various heating mechanisms in the current sheet. We find that reconnection energy conversion would dominate at low density (ne≲7×1018 cm−3) and electron-ion collisional drag at high-density (≳1019 cm−3).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

High-n Rydberg transition spectroscopy for heavy impurity transport studies in W7-X (invited)

Here, we present a novel spectroscopy approach to investigate impurity transport by analyzing line-radiation following high-n Rydberg transitions. While high-n Rydberg states of impurity ions are unlikely to be populated via impact excitation, they can be accessed by charge exchange (CX) reactions along the neutral beams in high-temperature plasmas. Hence, localized radiation of highly ionized impurities, free of passive contributions, can be observed at multiple wavelengths in the visible range. For the analysis and modeling of the observed Rydberg transitions, a technique for calculating effective emission coefficients is presented that can well reproduce the energy dependence seen in datasets available on the OPEN-ADAS database. By using the rate coefficients and comparing modeling results with the new high-n Rydberg CX measurements, impurity transport coefficients are determined with well-documented 2σ confidence intervals for the first time. This demonstrates that high-n Rydberg spectroscopy provides important constraints on the determination of impurity transport coefficients. By additionally considering Bolometer measurements, which provide constraints on the overall impurity emissivity and, therefore, impurity densities, error bars can be reduced even further.

Instruments & Instrumentation↗
DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Accuracy of Kohn–Sham density functional theory for warm- and hot-dense matter equation of state

We study the accuracy of Kohn–Sham density functional theory (DFT) for warm- and hot-dense matter (WDM and HDM). Specifically, considering a wide range of systems, we perform accurate ab initio molecular dynamics simulations with temperature-independent local/semilocal density functionals to determine the equations of state at compression ratios of 3x–7x and temperatures near 1 MK. We find very good agreement with path integral Monte Carlo benchmarks, while having significantly smaller error bars and smoother data, demonstrating the accuracy of DFT for the study of WDM and HDM at such conditions. In addition, using a Δ-machine learned force field scheme, we confirm that the DFT results are insensitive to the choice of exchange-correlation functional, whether local, semilocal, or nonlocal.

Suryanarayana, Phanish (ORCID:0000000151720049)↗
DOE OSTI2025

The Atacama Cosmology Telescope: DR6 constraints on extended cosmological models

We use new cosmic microwave background (CMB) primary temperature and polarization anisotropy measurements from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) to test foundational assumptions of the standard cosmological model, ΛCDM, and set constraints on extensions to it. We derive constraints from the ACT DR6 power spectra alone, as well as in combination with legacy data from the Planck mission. To break geometric degeneracies, we include ACT and Planck CMB lensing data and baryon acoustic oscillation data from DESI Year-1. To test the dependence of our results on non-ACT data, we also explore combinations replacing Planck with WMAP and DESI with BOSS, and further add supernovae measurements from Pantheon+ for models that affect the late-time expansion history. We verify the near-scale-invariance (running of the spectral index dn s /d ln k = 0.0062 ± 0.0052) and adiabaticity of the primordial perturbations. Neutrino properties are consistent with Standard Model predictions: we find no evidence for new light, relativistic species that are free-streaming (N eff = 2.86 ± 0.13, which combined with astrophysical measurements of primordial helium and deuterium abundances becomes N eff = 2.89 ± 0.11), for non-zero neutrino masses (∑m ν < 0.089 eV at 95% CL), or for neutrino self-interactions. We also find no evidence for self-interacting dark radiation (N idr < 0.134), or for early-universe variation of fundamental constants, including the fine-structure constant (α EM /α EM,0 = 1.0043 ± 0.0017) and the electron mass (m e /m e,0 = 1.0063 ± 0.0056). Our data are consistent with standard big bang nucleosynthesis (we find Y p = 0.2312 ± 0.0092), the COBE/FIRAS-inferred CMB temperature (we find T CMB = 2.698 ± 0.016 K), a dark matter component that is collisionless and with only a small fraction allowed as axion-like particles, a cosmological constant (w = -0.986 ± 0.025), and the late-time growth rate predicted by general relativity (γ = 0.663 ± 0.052). We find no statistically significant preference for a departure from the baseline ΛCDM model. In fits to models invoking early dark energy, primordial magnetic fields, or an arbitrary modified recombination history, we find H 0 = 69.9 +0.8 -1.5 , 69.1 ± 0.5, or 69.6 ± 1.0 km/s/Mpc, respectively; using BOSS instead of DESI BAO data reduces the central values of these constraints by 1–1.5 km/s/Mpc while only slightly increasing the error bars. In general, models introduced to increase the Hubble constant or to decrease the amplitude of density fluctuations inferred from the primary CMB are not favored over ΛCDM by our data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Toroidal injection angle dependence of EC assisted plasma initiation at DIII-D

An experimental scan of the electron cyclotron waves (EC) toroidal injection angle in plasma breakdown is performed at the DIII-D tokamak. The second harmonic, extraordinary mode EC is used for the study. The dependence of n e and T e on the EC injection angle cannot be conclusively drawn from this study due to the large error bars in the n e and T e measurements. On the other hand, high T e data points are observed in some discharges which can be explained by nonlinear heating. The D α emission measurement shows a clear relation between the breakdown time and the injection angle. An experimental investigation of the cause of the dependence of breakdown delay on the EC injection angle suggests that when the injection angle is oblique, the EC heating after the reflection at the inboard wall may become ineffective and cause the breakdown delay even when the EC heating before and directly upon injection remains effective. A preliminary run of the heat and transport balance code DYON indicates that the obtained dataset is suitable for a quantitative validation of EC absorption models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Machine learning materials properties with accurate predictions, uncertainty estimates, domain guidance, and persistent online accessibility

One compelling vision of the future of materials discovery and design involves the use of machine learning (ML) models to predict materials properties and then rapidly find materials tailored for specific applications. However, realizing this vision requires both providing detailed uncertainty quantification (model prediction errors and domain of applicability) and making models readily usable. At present, it is common practice in the community to assess ML model performance only in terms of prediction accuracy (e.g. mean absolute error), while neglecting detailed uncertainty quantification and robust model accessibility and usability. Here, we demonstrate a practical method for realizing both uncertainty and accessibility features with a large set of models. We develop random forest ML models for 33 materials properties spanning an array of data sources (computational and experimental) and property types (electrical, mechanical, thermodynamic, etc). All models have calibrated ensemble error bars to quantify prediction uncertainty and domain of applicability guidance enabled by kernel-density-estimate-based feature distance measures. All data and models are publicly hosted on the Garden-AI infrastructure, which provides an easy-to-use, persistent interface for model dissemination that permits models to be invoked with only a few lines of Python code. We demonstrate the power of this approach by using our models to conduct a fully ML-based materials discovery exercise to search for new stable, highly active perovskite oxide catalyst materials.

domain of applicability↗
DOE OSTI2026

Predicting the dark matter-baryon abundance ratio

We discuss relaxation solutions to the dark matter-baryon coincidence problem in the context of QCD axion dark matter. In relaxation solutions, a moduli dynamically adjusts the mass of dark matter and baryons until their energy densities are 𝒪⁡(1) the same. Because the QCD axion is heavily connected to QCD, scanning the QCD axion mass inherently also scans the proton mass. In the context of relaxation solutions, this implies that the ratio of dark matter to baryon abundances (Ω DM /Ω B ) is a ratio of beta functions showing that these models can only accommodate discrete values of ΩDM/ΩB thereby “predicting” the ratio of the dark matter to baryon abundances. The original composite axion model has only a single integer degree of freedom 𝑁, the size of the gauge group, and we show that when 𝑁 =8 the observed value of Ω DM /Ω B =5.36 is reproduced to within its percent level error bars. Novel tests of this model include more precise measurements of Ω DM /Ω B , a better lattice determination of the dependence of the proton mass on the high energy QCD gauge coupling, as well as more traditional tests such as fifth force experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Electromagnetic observables of open-shell nuclei from coupled-cluster theory

Here we develop a new method to describe electromagnetic observables of open-shell nuclei with two nucleons outside a closed shell. This approach combines the equation-of-motion coupled-cluster method for such systems and the Lorentz integral transform technique, expanding the applicability of coupled-cluster theory for electromagnetic observables beyond closed-(sub)shell nuclei. To validate this new approach, we compute the non-energy-weighted dipole sum rule and the dipole polarizability of 16,24 ⁢O in both the closed-(sub)shell and the new equation-of-motion coupled-cluster frameworks, finding agreement within error bars. We then analyze the evolution of the dipole polarizability along the oxygen and calcium isotopic chains. Our predictions agree well with available experimental data and other available theoretical calculations for the closed-(sub)shell 16,22 ⁢O and the open-shell 18 O. In the calcium isotopes, we observe that our dipole polarizability predictions for open-shell nuclei are lower than those of closed-(sub)shell nuclei. We expect that our predictions for 24 O and 54,56 ⁢Ca will motivate future experimental studies at the drip line.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Systematic study of projection biases in the weak lensing analysis of cosmic shear and the combination of galaxy clustering and galaxy-galaxy lensing

This paper presents the results of a systematic study of projection biases in the weak lensing analysis of cosmic shear and the combination of galaxy clustering and galaxy-galaxy lensing using data collected during the first year of running the Dark Energy Survey experiment. The study uses Lambda cold dark matter ( Λ CDM ) as the cosmological model and two-point correlation functions for the weak lensing (WL) analysis. The results in this paper show that, independent of the WL analysis, projection biases of more than 1 σ exist and are a function of the position of the true values of the parameters h , n s , Ω b , and Ω ν h 2 with respect to their prior probabilities. For cosmic shear, and the combination of galaxy clustering and galaxy-galaxy lensing, this study shows that the coverage probability of the 68.27% credible intervals ranges from as high as 93% to as low as 16% and that these credible intervals are inflated, on average, by 29% for cosmic shear and 20% for the combination of galaxy clustering and galaxy-galaxy lensing. The results of the study also show that, in six out of nine tested cases, the reduction in error bars obtained by transforming credible intervals into confidence intervals is equivalent to an increase in the amount of data by a factor of 3.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Multiscale Physics of Atomic Nuclei from First Principles

Atomic nuclei exhibit multiple energy scales ranging from hundreds of MeV in binding energies to fractions of an MeV for low-lying collective excitations. As the limits of nuclear binding are approached near the neutron and proton drip lines, traditional shell structure starts to melt with an onset of deformation and an emergence of coexisting shapes. It is a long-standing challenge to describe this multiscale physics starting from nuclear forces with roots in quantum chromodynamics. Here, we achieve this within a unified and nonperturbative quantum many-body framework that captures both short- and long-range correlations starting from modern nucleon-nucleon and three-nucleon forces from chiral effective field theory. The short-range (dynamic) correlations which account for the bulk of the binding energy are included within a symmetry-breaking framework, while long-range (static) correlations (and fine details about the collective structure) are included by employing symmetry projection techniques. Our calculations accurately reproduce—within theoretical error bars—available experimental data for low-lying collective states and the electromagnetic quadrupole transitions in 20−30 Ne. In addition, we reveal coexisting spherical and deformed shapes in 30 Ne, which indicates the breakdown of the magic neutron number 𝑁 = 20 as the key nucleus 28 O is approached, and we predict that the drip line nuclei 32,34 Ne are strongly deformed and collective. By developing reduced-order models for symmetry-projected states, we perform a global sensitivity analysis and find that the subleading singlet 𝑆-wave contact and a pion-nucleon coupling strongly impact nuclear deformation in chiral effective field theory. The techniques developed in this work clarify how microscopic nuclear forces generate the multiscale physics of nuclei spanning collective phenomena as well as short-range correlations and allow one to capture emergent and dynamical phenomena in finite fermion systems such as atom clusters, molecules, and atomic nuclei.

74 ATOMIC AND MOLECULAR PHYSICS↗