Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Uncertainty”

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 883 records · Page 49

Measurements of the thermal Sunyaev-Zel’dovich effect with ACT and DESI luminous red galaxies

Cosmic Microwave Background (CMB) photons scatter off the free-electron gas in galaxies and clusters, allowing us to use the CMB as a backlight to probe the gas in and around low-redshift galaxies. The thermal Sunyaev-Zel’dovich effect, sourced by hot electrons in high-density environments, measures the thermal pressure of the target objects, shedding light on halo thermodynamics and galaxy formation, and providing a path toward understanding the baryon distribution around cosmic structures. We use a combination of high-resolution CMB maps from the Atacama Cosmology Telescope and photometric luminous red galaxy catalogs from the Dark Energy Spectroscopic Instrument to measure the thermal Sunyaev-Zel’dovich signal in four redshift bins from 𝑧 = 0.4 to 𝑧 = 1.2, with a combined detection significance of 19⁢𝜎 when stacking on the fiducial CMB Compton-𝑦 map. We discuss possible sources of contamination, finding that residual dust emission associated with the target galaxies is important and limits current analyses. We discuss several mitigation strategies and quantify the residual modeling uncertainty. Furthermore, this work complements closely related measurements of the kinematic Sunyaev-Zel’dovich and weak lensing of the same galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Enhanced antineutrino emission from 𝛽 decay in core-collapse supernovae with self-consistent weak decay rates

Nuclear weak-interaction rates are known to exert a prominent effect in the late stages of stellar collapse. Despite their importance, most studies to date on core-collapse supernovae (CCSNe) have focused primarily on the effects of electron captures, neglecting 𝛽 decay contributions. In this work, we present the first CCSNe simulation incorporating global 𝛽 decay rates from a microscopic theory. These are enabled by a large-scale evaluation of both electron capture and 𝛽 decay rates, obtained self-consistently utilizing the relativistic energy density functional theory and finite-temperature quasiparticle random-phase approximation. Including 𝛽 decay leads to a dramatic enhancement of the pre-bounce antineutrino signal as the antineutrino emissivity increases by more than two orders of magnitude and the luminosity by a factor of 50 relative to thermal emission alone, while the average antineutrino energy increases by over 1 MeV. Furthermore, it is expected that these new rates could help us constrain the model uncertainties related to weak-interaction processes, improving the prediction of antineutrino signal during the final stages of stellar death.

Beta decay↗

Global Bayesian analysis of J/𝜓 photoproduction on proton and lead targets

We perform a global Bayesian analysis of diffractive J/𝜓 production in 𝛾 + 𝑝 and 𝛾 + Pb collisions using a color glass condensate (CGC) based calculation framework. As past calculations have shown that CGC-based models typically overpredict the J/𝜓 production in 𝛾 + Pb collisions at high center of mass energy, we address the question of whether it is possible to describe coherent and incoherent diffractive J/𝜓 data from 𝛾 + 𝑝 collisions at Hadron–Electron Ring Accelerator (HERA) and the LHC, and from 𝛾 + Pb collisions at the LHC simultaneously. Our results indicate that a simultaneous description of 𝛾 + 𝑝 and 𝛾 + Pb data is challenging, with results improving when an overall 𝐾-factor—scaling 𝛾 + 𝑝 and 𝛾 + Pb cross sections to absorb model uncertainties—is introduced.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Proton pairing in neutron stars from chiral effective field theory

Here, we study the 1 S 0 proton pairing gap in beta-equilibrated neutron star matter within the framework of chiral effective field theory. We focus on the role of three-body forces, which strongly modify the effective proton-proton spin-singlet interaction in dense matter. We find that three-body forces generically reduce both the size of the pairing gap and the maximum density at which proton pairing may occur. The pairing gap is computed within BCS theory using a single-particle dispersion relation calculated up to second order in perturbation theory. Model uncertainties are estimated by varying the nuclear potential (its order in the chiral expansion and high-momentum cutoff) and the choice of single-particle spectrum in the gap equation. We find that a second-order perturbative treatment of the single-particle spectrum suppresses the proton 1 S 0 pairing gap relative to the use of a free spectrum. We estimate the critical temperature for the onset of proton superconductivity to be T c = (3:2–5:1) 10 9 K, which is consistent with previous theoretical results in the literature and marginally within the range deduced from a recent Bayesian analysis of neutron star cooling observations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Prehydrodynamic evolution and its impact on quark-gluon plasma signatures

State-of-the-art hydrodynamic models of heavy-ion collisions have considerable theoretical model uncertainties in the description of the very early prehydrodynamic stage. Here we add a new computational module, K T Iso, that describes the prehydrodynamic evolution kinetically, based on the relativistic Boltzmann equation with collisions treated in the isotropization time approximation. As a novelty, K T Iso allows for the inclusion and evolution of initial-state momentum anisotropies. To maintain computational efficiency K T Iso assumes strict longitudinal boost invariance and allows collisions to isotropize only the transverse momenta. We use it to explore the sensitivity of hadronic observables measured in relativistic heavy-ion collisions to initial-state momentum anisotropies and microscopic scattering during the prehydrodynamic stage.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gluon distributions and their applications to Ioffe-time distributions

We investigate unpolarized and polarized gluon distributions and their applications to the Ioffe-time distributions, which are related to lattice QCD calculations of parton distribution functions. Guided by the counting rules based on the perturbative QCD at large momentum fraction x and the color coherence of gluon couplings at small x, we parametrize gluon distributions in the helicity basis. By fitting the unpolarized gluon distribution, the inferred polarized gluon distribution from our parametrization agrees with the one from global analysis. A simultaneous fit to both unpolarized and polarized gluon distributions is also performed to explore the model uncertainty. The agreement with the global analysis supports the (1 - x) power suppression of the helicity-antialigned distribution relative to the helicity-aligned distribution. The corresponding Ioffe-time distributions and their asymptotic expansions are calculated from the gluon distributions. Our results of the Ioffe-time distributions can provide guidance to the extrapolation of lattice QCD data to the region lacking precise gluonic matrix elements. Therefore, they can help regulate the ill-posed inverse problem associated with extracting the gluon distributions from discrete data from first-principle calculations, which are available in a limited range of the nucleon momentum and the spatial separation between the gluonic currents. Given various limitations in obtaining lattice QCD data at large Ioffe time, phenomenological approaches can provide important complementary information to extract the gluon distributions in the entire momentum fraction region, especially at small x. The possibility of investigating higher-twist effects and other systematic uncertainties in the contemporary first-principle calculations of parton distributions from phenomenologically well-determined Ioffe-time distributions in the large Ioffe-time region is also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dark Energy Survey Year 3 Results: Three-point shear correlations and mass aperture moments

We present high signal-to-noise measurements of three-point shear correlations and the third moment of the mass aperture statistic using the first 3 years of data from the Dark Energy Survey. We additionally obtain the first measurements of the configuration and scale dependence of the four three-point shear correlations which carry cosmological information. With the third-order mass aperture statistic, we present tomographic measurements over angular scales of 4 to 60 arcminutes with a combined statistical significance of 15.0 σ . Using the tomographic information and measuring also the second-order mass aperture, we additionally obtain a skewness parameter and its redshift evolution. We find that the amplitudes and scale-dependence of these shear 3pt functions are in qualitative agreement with measurements in a mock galaxy catalog based on N-body simulations, indicating promise for including them in future cosmological analyses. We validate our measurements by showing that B-modes, parity-violating contributions and PSF modeling uncertainties are negligible, and determine that the measured signals are likely to be of astrophysical and gravitational origin.

79 ASTRONOMY AND ASTROPHYSICS↗

Muon spin force

Current discrepancy between the measurement and the prediction of the muon anomalous magnetic moment can be resolved in the presence of a long-range force created by ordinary atoms acting on the muon spin via axial-vector and/or pseudoscalar coupling, and requiring a tiny O⁡(10 –13 eV) spin energy splitting between muon state polarized in the vertical direction. We suggest that an extension of the muon spin resonance (μ⁢SR) experiments can provide a definitive test of this class of models. We also derive indirect constraints on the strength of the muon spin force by considering the muon-loop-induced interactions between nuclear spin and external directions. The limits on the muon spin force extracted from the comparison of 199 Hg/ 201 Hg and 129 Xe/ 131 Xe spin precession are strong for the pseudoscalar coupling but are significantly relaxed for the axial-vector one. These limits suffer from significant model uncertainties, poorly known proton/neutron spin content of these nuclei, and therefore do not exclude the possibility of a muon spin force relevant for the muon g – 2.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Machine learning assisted unfolding for neutrino cross-section measurements with the OmniFold technique

The choice of unfolding method for a cross-section measurement is tightly coupled to the model dependence of the efficiency correction and the overall impact of cross-section modeling uncertainties in the analysis. A key issue is the dimensionality used in unfolding, as the kinematics of all outgoing particles in an event typically affect the reconstruction performance in a neutrino detector. OmniFold is an unfolding method that iteratively reweights a simulated dataset, using machine learning to utilize arbitrarily high-dimensional information, that has previously been applied to proton-proton and proton-electron datasets. This paper demonstrates OmniFold’s application to a neutrino cross-section measurement for the first time using a public T2K near detector simulated dataset, comparing its performance with traditional approaches using a mock data study.

Machine learning↗

Observer-Based Nonlinear Control Scheme to Reduce Oscillations and Zero Crossing in Skid-Steer Vehicles

Motion sickness is a common condition experienced by drivers of skid-steer vehicles, primarily caused by zero crossing and oscillations in undamped systems. This study proposes an observer-based nonlinear control scheme to reduce transient oscillations and zero-crossing phenomena in skid-steer vehicles, thereby potentially alleviating motion sickness. Reducing transient oscillations and zero crossing in the transient response may alleviate motion sickness. A nonlinear damping controller is designed to improve transient response by reducing oscillations and zero-crossing. To design the controller, a reduced-order kinematic model based on coordinate transformation is developed. This transformation not only converts the system modeling into a controllable form but also enhances control performance. Modeling error is addressed by considering the distance between the center of the vehicle and the sensor location. Despite these improvements, model uncertainties and external disturbances remain, which may degrade control performance. To ensure robustness and estimate such disturbances, a high-order sliding mode observer (HOSMO) is incorporated. The effectiveness of the proposed method is validated through MATLAB/Simulink and TruckMaker simulations. From the simulation results, it was shown that the proposed method reduced the mean squared error of the tracking error to within 10 % compared to the state feedback controller with the HOSMO.

Seo, Jiwon [Chung-Ang University, Seoul (Korea, Re↗

Change Detection in Constellations of Buried Objects Extracted From Ground-Penetrating Radar Data

Detection of deliberately buried objects in ground-penetrating radar (GPR) data acquired along a path is a clutter-limited problem. Detection-false alarm rate performance can be improved by replacing the detection statistic with a change statistic that incorporates information from previous path traversals. A constellation matching approach is developed for buried-object change detection in GPR data. Network topologies of buried objects detected in GPR data from previous path traversals are maintained in a constellation database. Localized groups of buried objects newly detected on the latest path traversal are matched to the constellation. Buried objects from the latest path traversal whose locations or strengths cannot be reconciled with the constellation are identified as changes. The system has one component that generates constellation databases offline and another component suitable for change detection in real time. It can tolerate paths with significant translational misalignments. The system uses the following: 1) a customized translational relaxation algorithm for point pattern matching that incorporates detection strength and a probabilistic uncertainty model for buried-object location into the objective function and 2) a change statistic that accounts for the magnitude of change relative to predicted detection strength. A constellation database can typically be generated offline from a single path traversal roughly two orders of magnitude faster than the time typically required for a vehicle to travel the extent of the path. Database sizes are typically four to five orders of magnitude smaller than the data sets of GPR signal scans or focused 3-D GPR images that they were generated from. On bumpy dirt roads buried exclusively with nonmetallic objects at various depths, detection-false alarm rate performance is shown to be significantly better for our change statistics than for our detection statistics.

47 OTHER INSTRUMENTATION↗

Global distribution, formation and fate of mineral‐associated soil organic matter under a changing climate: A trait‐based perspective

Abstract Soil organic matter (SOM) is the largest actively cycling reservoir of terrestrial carbon (C), and the majority of SOM in Earth's mineral soils (~65%) is mineral‐associated organic matter (MAOM). Thus, the formation and fate of MAOM can exert substantial influence on the global C cycle. To predict future changes to Earth's climate, it is critical to mechanistically understand the processes by which MAOM is formed and decomposed, and to accurately represent this process‐based understanding in biogeochemical and Earth system models. In this review, we use a trait‐based framework to synthesize the interacting roles of plants, soil micro‐organisms, and the mineral matrix in regulating MAOM formation and decomposition. Our proposed framework differentiates between plant and microbial traits that influence total OM inputs to the soil (‘feedstock traits’) versus traits that influence the proportion of OM inputs that are ultimately incorporated into MAOM (‘MAOM formation traits’). We discuss how these feedstock and MAOM formation traits may be altered by warming, altered precipitation and elevated carbon dioxide. At a planetary scale, these feedstock and MAOM formation traits help shape the distribution of MAOM across Earth's biomes, and modulate biome‐specific responses of MAOM to climate change. We leverage a global synthesis of MAOM measurements to provide estimates of the total amount of MAOM‐C globally (~840–1540 Pg C; 34%–51% of total terrestrial organic C), and its distribution across Earth's biomes. We show that MAOM‐C concentration is highest in temperate forests and grasslands, and lowest in shrublands and savannas. Grasslands and croplands have the highest proportion of soil organic carbon (SOC) in the MAOM fraction (i.e. the MAOM‐C:SOC ratio), while boreal forests and tundra have the lowest MAOM‐C:SOC ratio. Drawing on our trait framework, we then review experimental data and posit the effects of climate change on MAOM pools in different biomes. We conclude by discussing how MAOM is integrated into soil C models, and how feedstock and MAOM formation traits may be included in these models. We also summarize the projected fate of MAOM under climate change scenarios (Representative Concentration Pathways 4.5 and 8.5) and discuss key model uncertainties. Read the free Plain Language Summary for this article on the Journal blog.

54 ENVIRONMENTAL SCIENCES↗

Data-Efficient Dimensionality Reduction and Surrogate Modeling of High-Dimensional Stress Fields

Tensor datatypes representing field variables like stress, displacement, velocity, etc., have increasingly become a common occurrence in data-driven modeling and analysis of simulations. Numerous methods [such as convolutional neural networks (CNNs)] exist to address the meta-modeling of field data from simulations. As the complexity of the simulation increases, so does the cost of acquisition, leading to limited data scenarios. Modeling of tensor datatypes under limited data scenarios remains a hindrance for engineering applications. Here, in this article, we introduce a direct image-to-image modeling framework of convolutional autoencoders enhanced by information bottleneck loss function to tackle the tensor data types with limited data. The information bottleneck method penalizes the nuisance information in the latent space while maximizing relevant information making it robust for limited data scenarios. The entire neural network framework is further combined with robust hyperparameter optimization. We perform numerical studies to compare the predictive performance of the proposed method with a dimensionality reduction-based surrogate modeling framework on a representative linear elastic ellipsoidal void problem with uniaxial loading. The data structure focuses on the low-data regime (fewer than 100 data points) and includes the parameterized geometry of the ellipsoidal void as the input and the predicted stress field as the output. The results of the numerical studies show that the information bottleneck approach yields improved overall accuracy and more precise prediction of the extremes of the stress field. Additionally, an in-depth analysis is carried out to elucidate the information compression behavior of the proposed framework.

artificial intelligence↗

SU-E-T-202: Comparison of 4D-Measurement-Guided Dose Reconstructions (MGDR) with COMPASS and OCTAVIUS 4D System

Purpose: To cross-validate the MGDR of COMPASS (IBA dosimetry, GmbH, Germany) and OCTAVIUS 4D system (PTW, Freiburg, Germany). Methods: Volumetric-modulated arc plans (5 head-and-neck and 3 prostate) collapsed to 40° gantry on the OCTAVIUS 4D phantom in QA mode on Monaco v5.0 (Elekta, CMS, Maryland Heights, MO) were delivered on a Elekta Agility linac. This study was divided into two parts: (1) error-free measurements by gantry-mounted EvolutionXX 2D array were reconstructed in COMPASS (IBA dosimetry, GmbH, Germany), and by OCTAVIUS 1500 array in Versoft v6.1 (PTW, Freiburg, Germany) to obtain the 3D doses (COM4D and OCTA4D). COM4D and OCTA4D were compared to the raw measurement (OCTA3D) at the same detector plane for which OCTAVIUS 1500 was perpendicular to 0° gantry axis while the plans were delivered at gantry 40°; (2) beam steering errors of energy (Hump=-2%) and symmetry (2T=+2%) were introduced during the delivery of 5 plans to compare the MGDR doses COM4D-Hump (COM4D-2T), OCTA4D-Hump (OCTA4D-2T), with raw doses OCTA3D-Hump (OCTA3D-2T) and with OCTA3D to assess the error reconstruction and detection ability of MGDR tools. All comparisons used Υ-criteria of 2%(local dose)/2mm and 3%/3mm. Results: Averaged Υ passing rates were 85% and 96% for COM4D,and 94% and 99% for OCTA4D at 2%/2mm and 3%/3mm criteria respectively. For error reconstruction, COM4D-Hump (COM4D-2T) showed 81% (93%) at 2%/2mm and 94% (98%) at 3%/3mm, while OCTA4D-Hump (OCTA4D-2T) showed 96% (96%) at 2%/2mm and 99% (99%) at 3%/3mm. For error detection, OCTA3D doses were compared to COM4D-Hump (COM4D-2T) showing Υ passing rates of 93% (93%) at 2%/2mm and 98% (98%), and to OCTA4D-Hump (OCTA4D -2T) showing 94% (99%) at 2%/2mm and 81% (96%) at 3%/3mm, respectively. Conclusion: OCTAVIUS MGDR showed better agreement to raw measurements in both error- and error-free comparisons. COMPASS MGDR deviated from the raw measurements possibly owing to beam modeling uncertainty.

60 APPLIED LIFE SCIENCES↗

SU-F-J-54: Towards Real-Time Volumetric Imaging Using the Treatment Beam and KV Beam

Purpose: Existing real-time imaging uses dual (orthogonal) kV beam fluoroscopies and may result in significant amount of extra radiation to patients, especially for prolonged treatment cases. In addition, kV projections only provide 2D information, which is insufficient for in vivo dose reconstruction. We propose real-time volumetric imaging using prior knowledge of pre-treatment 4D images and real-time 2D transit data of treatment beam and kV beam. Methods: The pre-treatment multi-snapshot volumetric images are used to simulate 2D projections of both the treatment beam and kV beam, respectively, for each treatment field defined by the control point. During radiation delivery, the transit signals acquired by the electronic portal image device (EPID) are processed for every projection and compared with pre-calculation by cross-correlation for phase matching and thus 3D snapshot identification or real-time volumetric imaging. The data processing involves taking logarithmic ratios of EPID signals with respect to the air scan to reduce modeling uncertainties in head scatter fluence and EPID response. Simulated 2D projections are also used to pre-calculate confidence levels in phase matching. Treatment beam projections that have a low confidence level either in pre-calculation or real-time acquisition will trigger kV beams so that complementary information can be exploited. In case both the treatment beam and kV beam return low confidence in phase matching, a predicted phase based on linear regression will be generated. Results: Simulation studies indicated treatment beams provide sufficient confidence in phase matching for most cases. At times of low confidence from treatment beams, kV imaging provides sufficient confidence in phase matching due to its complementary configuration. Conclusion: The proposed real-time volumetric imaging utilizes the treatment beam and triggers kV beams for complementary information when the treatment beam along does not provide sufficient confidence for phase matching. This strategy minimizes the use of extra radiation to patients. This project is partially supported by a Varian MRA grant.

60 APPLIED LIFE SCIENCES↗

Enabling Hyper-Differential Sensitivity Analysis for Ill-Posed Inverse Problems

Inverse problems constrained by partial differential equations (PDEs) play a critical role in model development and calibration. In many applications, there are multiple uncertain parameters in a model that must be estimated. However, high dimensionality of the parameters and computational complexity of the PDE solves make such problems challenging. A common approach is to reduce the dimension by fixing some parameters (which we will call auxiliary parameters) to a best estimate and use techniques from PDE-constrained optimization to estimate the other parameters. In this article, hyper-differential sensitivity analysis (HDSA) is used to assess the sensitivity of the solution of the PDE-constrained optimization problem to changes in the auxiliary parameters. Foundational assumptions for HDSA require satisfaction of the optimality conditions which are not always practically feasible as a result of ill-posedness in the inverse problem. Here we introduce novel theoretical and computational approaches to justify and enable HDSA for ill-posed inverse problems by projecting the sensitivities on likelihood informed subspaces and defining a posteriori updates. Our proposed framework is demonstrated on a nonlinear multiphysics inverse problem motivated by estimation of spatially heterogeneous material properties in the presence of spatially distributed parametric modeling uncertainties.

97 MATHEMATICS AND COMPUTING↗

Virtual Compton scattering measurements in the nucleon resonance region

We report on new measurements of the electric Generalized Polarizability (GP) of the proton \( \alpha_{E}\) in a kinematic region where a puzzling dependence on momentum transfer has been observed, and we have found that \( \alpha_E = (5.3 \pm 0.6_{stat} \pm 1.3_{sys}) 10^{-4}\) fm 3 at \( Q^{2}=0.20\) (GeV/ c ) 2 . The new measurements, when considered along with the rest of the world data, suggest that \( \alpha_E\) can be described by either a local plateau or by an enhancement in the region \( Q^2=0.20\) (GeV/ c ) 2 to 0.33 (GeV/ c ) 2 . The experiment also provides the first measurement of the Coulomb quadrupole amplitude in the \( N \rightarrow \Delta\) transition through the exploration of the \( p(e,e'p)\gamma\) reaction. The new measurement gives \( CMR = (-4.4 \pm 0.8_{stat} \pm 0.6_{sys})\%\) at \( Q^2=0.20\) (GeV/ c ) 2 and is consistent with the results from the pion electroproduction world data. It has been obtained using a completely different extraction method, and therefore represents a strong validation test of the world data model uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quasi-Newton Variational Bayes v.3.0

SAND2024-09052O Quasi-Newton Variational Bayes (QNVB) is a training algorithm that performs high-dimensional variational inference on machine learning models. It provides mechanisms to calibrate and control model uncertainty during training. The PyTorch implementation of projective integral updates for Gaussian mean-fields supports the manuscript, "Projective Integral Updates for High-Dimensional Variational Inference." Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Duersch, Jed↗