Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parameter estimation using data”

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 1,027 records · Page 57

Data-driven discovery of dynamics from time-resolved coherent scattering

Coherent X-ray scattering (CXS) techniques are capable of interrogating dynamics of nano- to mesoscale materials systems at time scales spanning several orders of magnitude. However, obtaining accurate theoretical descriptions of complex dynamics is often limited by one or more factors—the ability to visualize dynamics in real space, computational cost of high-fidelity simulations, and effectiveness of approximate or phenomenological models. In this work, we develop a data-driven framework to uncover mechanistic models of dynamics directly from time-resolved CXS measurements without solving the phase reconstruction problem for the entire time series of diffraction patterns. Our approach uses neural differential equations to parameterize unknown real-space dynamics and implements a computational scattering forward model to relate real-space predictions to reciprocal-space observations. This method is shown to recover the dynamics of several computational model systems under various simulated conditions of measurement resolution and noise. Moreover, the trained model enables estimation of long-term dynamics well beyond the maximum observation time, which can be used to inform and refine experimental parameters in practice. Finally, we demonstrate an experimental proof-of-concept by applying our framework to recover the probe trajectory from a ptychographic scan. Our proposed framework bridges the wide existing gap between approximate models and complex data.

36 MATERIALS SCIENCE↗

Solving high-dimensional inverse problems using amortized likelihood-free inference with noisy and incomplete data

Here, we present a likelihood-free probabilistic inversion method based on normalizing flows for high-dimensional inverse problems. The proposed method is composed of two complementary networks: a summary network for data compression and an inference network for parameter estimation. The summary network encodes raw observations into a fixed-size vector of summary features, while the inference network generates samples of the approximate posterior distribution of the model parameters based on these summary features. The posterior samples are produced in a deep generative fashion by sampling from a latent Gaussian distribution and passing these samples through an invertible transformation. We construct this invertible transformation by sequentially alternating conditional invertible neural network and conditional neural spline flow layers. The summary and inference networks are trained simultaneously. We apply the proposed method to an inversion problem in groundwater hydrology to estimate the posterior distribution of the log-conductivity field conditioned on spatially sparse time-series observations of the system’s hydraulic head responses. The conductivity field is represented with 706 degrees of freedom in the considered problem. Comparison with the likelihood-based iterative ensemble smoother PEST-IES method demonstrates that the proposed method accurately estimates the parameter posterior distribution and the observations’ predictive posterior distribution at a fraction of the inference time of PEST-IES.

conditional invertible neural network↗

Applying Bayesian parameter estimation to relativistic heavy-ion collisions: Simultaneous characterization of the initial state and quark-gluon plasma medium

We quantitatively estimate properties of the quark-gluon plasma created in ultrarelativistic heavy-ion collisions utilizing Bayesian statistics and a multiparameter model-to-data comparison. The study is performed using a recently developed parametric initial condition model, T R ENTo, which interpolates among a general class of particle production schemes, and a modern hybrid model which couples viscous hydrodynamics to a hadronic cascade.We calibrate the model to multiplicity, transverse momentum, and flow data and report constraints on the parametrized initial conditions and the temperature-dependent transport coefficients of the quark-gluon plasma. We show that initial entropy deposition is consistent with a saturation-based picture, extract a relation between the minimum value and slope of the temperature-dependent specific shear viscosity, and find a clear signal for a nonzero bulk viscosity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: theoretical systematics and Baryon Acoustic Oscillations in the galaxy correlation function

We investigate the potential sources of theoretical systematics in the anisotropic Baryon Acoustic Oscillation (BAO) distance scale measurements from the clustering of galaxies in configuration space using the final Data Release (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS). We perform a detailed study of the impact on BAO measurements from choices in the methodology such as fiducial cosmology, clustering estimators, random catalogues, fitting templates, and covariance matrices. The theoretical systematic uncertainties in BAO parameters are found to be 0.002 in the isotropic dilation α and 0.003 in the quadrupolar dilation ε. The leading source of systematic uncertainty is related to the reconstruction techniques. Theoretical uncertainties are sub-dominant compared with the statistical uncertainties for BOSS survey, accounting 0.2σ stat for α and 0.25σ stat for ε (σ α , stat ~ 0.010 and σ ε , stat ~ 0.012, respectively). We also present BAO-only distance scale constraints from the anisotropic analysis of the correlation function. Our constraints on the angular diameter distance D A (z) and the Hubble parameter H(z), including both statistical and theoretical systematic uncertainties, are 1.5 percent and 2.8 percent at z eff = 0.38, 1.4 percent and 2.4 percent at z eff = 0.51, and 1.7 percent and 2.6 percent at z eff = 0.61. This paper is part of a set that analyses the final galaxy clustering data set from BOSS. The measurements and likelihoods presented here are cross-checked with other BAO analysis in Alam et al. The systematic error budget concerning the methodology on post-reconstruction BAO analysis presented here is used in Alam et al. to produce the final cosmological constraints from BOSS.

79 ASTRONOMY AND ASTROPHYSICS↗

Non-stationary Geostatistical Modeling: A Case Study Comparing LVA Estimation Frameworks

Incorporating locally varying anisotropy (LVA) in geostatistical modeling improves estimates for structurally complex domains where a single set of anisotropic parameters modeled globally do not account for all geological features. In this work, the properties of two LVA-geostatistical modeling frameworks are explored through application to a complexly folded gold deposit in Ghana. The inference of necessary parameters is a significant requirement of geostatistical modeling with LVA; this work focuses on the case where LVA orientations, derived from expert geological interpretation, are used to improve the grade estimates. The different methodologies for inferring the required parameters in this context are explored. The results of considering different estimation frameworks and alternate methods of parameterization are evaluated with a cross-validation study, as well as visual inspection of grade continuity along select cross sections. Results show that stationary methodologies are outperformed by all LVA techniques, even when the LVA framework has minimal guidance on parameterization. Findings also show that additional improvements are gained by considering parameter inference where the LVA orientations and point data are used to infer the local range of anisotropy. Considering LVA for geostatistical modeling of the deposit considered in this work results in better reproduction of curvilinear geological features.

58 GEOSCIENCES↗

SPT-3G D1: Quadratic-Estimator CMB Lensing Reconstruction and Cosmology

We present a map of the cosmic microwave background (CMB) lensing potential reconstructed from observations taken during the 2019 and 2020 seasons with the third-generation camera on the South Pole Telescope (SPT), covering the $1500\,{\rm deg}^{2}$ SPT-3G Main field, referred to as the SPT-3G D1 dataset. From the multi-frequency temperature and polarization data, we reconstruct the CMB lensing field using a quadratic estimator that jointly accounts for the $T$, $E$, and $B$ fields and their covariance. The resulting lensing map is dominated by polarization information for $L \lesssim 600$ and provides the highest signal-to-noise measurement per mode reported to date. With nuisance parameters fixed to their best-fit values, we measure a lensing amplitude consistent with unity at $2\%$ precision relative to the $Λ$CDM model that best fits the combined Planck, ACT DR6, and SPT-3G D1 $TT/TE/EE$ likelihoods (${\rm CMB}_{\rm SPA}$). We further measure the structure-growth parameter $σ_{8}Ω_{\rm m}^{0.25}$ to be $0.6046\pm0.0096$ from the SPT-3G D1 lensing spectrum alone and $0.6020\pm0.0084$ when combined with ACT DR6 and Planck PR4 CMB lensing. By further combining this with ${\rm CMB}_{\rm SPA}$ and the latest DESI DR2 BAO data, we obtain $\sum m_ν < 0.072\,\mathrm{eV}$ (95% C.L.) when allowing the neutrino mass to vary within $Λ$CDM. Compared with previous work, the better agreement of our measurement with DESI DR2 BAO yields both this relaxed upper bound and reduced ($\mathord{\sim}2σ$) preferences for nonzero spatial curvature and for deviations of $(w_0,w_a)$ from the $Λ$CDM expectation. When we combine CMB lensing with the DES Y3 3$\times$2pt analysis, we obtain $S_{8}=0.811\pm0.011$, corresponding to a $1.4\%$ constraint on the late-time clustering amplitude. This precision is competitive with that obtained from the primary CMB within $Λ$CDM.

Omori, Y. [Chicago U., Astron. Astrophys. Ctr.; Ch↗

BayesMT: A Probabilistic Bayesian Framework for the Seismic Moment Tensor

Moment tensors (MTs) have long been used in earthquake and explosion source analysis, and there has been a renewed interest in how they can inform us about the seismic source, particularly in the geophysical monitoring community due to its application in event identification and yield analysis. However, parameter uncertainties in seismic MT inversion are rarely available. The inverse procedure often does not quantify MT model errors such as event location, data noise and Earth model that are essential for estimating solution robustness. To address this need, we propose to adopt the Bayesian probabilistic framework to incorporate uncertainties in MT inversions. In this study, we present the theoretical background of a probabilistic Bayesian framework for MT inversion accounting for model and measurements errors and illustrate the implementation of the method using a synthetic example.

58 GEOSCIENCES↗

DeepCMB: Lensing reconstruction of the cosmic microwave background with deep neural networks

Next-generation cosmic microwave background (CMB) experiments will have lower noise and therefore increased sensitivity, enabling improved constraints on fundamental physics parameters such as the sum of neutrino masses and the tensor-to-scalar ratio r. Achieving competitive constraints on these parameters requires high signal-to-noise extraction of the projected gravitational potential from the CMB maps. Standard methods for reconstructing the lensing potential employ the quadratic estimator (QE). However, the QE performs suboptimally at the low noise levels expected in upcoming experiments. Other methods, like maximum likelihood estimators (MLE), are under active development. In this work, we demonstrate reconstruction of the CMB lensing potential with deep convolutional neural networks (CNN) - ie, a ResUNet. The network is trained and tested on simulated data, and otherwise has no physical parametrization related to the physical processes of the CMB and gravitational lensing. We show that, over a wide range of angular scales, ResUNets recover the input gravitational potential with a higher signal-to-noise ratio than the QE method, reaching levels comparable to analytic approximations of MLE methods. We demonstrate that the network outputs quantifiably different lensing maps when given input CMB maps generated with different cosmologies. We also show we can use the reconstructed lensing map for cosmological parameter estimation. This application of CNN provides a few innovations at the intersection of cosmology and machine learning. First, while training and regressing on images, we predict a continuous-variable field rather than discrete classes. Second, we are able to establish uncertainty measures for the network output that are analogous to standard methods. We expect this approach to excel in capturing hard-to-model non-Gaussian astrophysical foreground and noise contributions.

79 ASTRONOMY AND ASTROPHYSICS↗

Heisenberg model analysis on inelastic powder neutron scattering data using parent and K doped BaMn 2 As 2 samples

Low temperature powder inelastic neutron scattering measurements were performed on three different powder samples; parent BaMn 2 As 2 ,12.5% K-doped Ba 0.875 K 0.125 Mn 2 As 2 and 25% K-doped Ba (0.75) K 0.25 Mn 2 As 2 . The Heisenberg Model involving J 1 -J 2 -J z coupling constants were compared to the data by a powder integration routine using Monte Carlo integration methods. The best magnetic parameters were selected using a chi-square test where model intensities were compared to the full (q,E) dependence of magnetic scattering. A key step to this analysis is the characterization of the background which is formed mostly by phonon scattering intensities along with other sources including the magnetic impurity scattering events. The calculated powder magnetic intensities added to the estimated background obtained from the non-magnetic high momentum transfer region. The agreement between the simulated and the raw data enabled us to perform quantitative analysis of the unreacted MnO impurities. Overall, this is another confirmation along with earlier studies using this technique, that magnetic exchange constants can be calculated within an acceptable range with a very quick inelastic neutron powder experiment without need for a single crystal sample.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Exploring HOD-dependent systematics for the DESI 2024 Full-Shape galaxy clustering analysis

We analyze the robustness of the DESI 2024 cosmological inference from the full shape of the galaxy power spectrum to uncertainties in the Halo Occupation Distribution (HOD) model of the galaxy-halo connection and the choice of priors on nuisance parameters. We assess variations in the recovered cosmological parameters across a range of mocks populated with different HOD models and find that shifts are often greater than 20% of the expected statistical uncertainties from the DESI data. We encapsulate the effect of such shifts in terms of a systematic covariance term, C HOD , and an additional diagonal contribution quantifying the impact of our choice of nuisance parameter priors on the ability of the effective field theory (EFT) model to correctly recover the cosmological parameters of the simulations. These two covariance contributions are designed to be added to the usual covariance term, C stat , describing the statistical uncertainty in the power spectrum measurement, in order to fairly represent these sources of systematic uncertainty. This novel approach should be more general and robust to the choice of model or additional external datasets used in cosmological fits than the alternative approach of adding systematic uncertainties to the recovered marginalised parameter posteriors. We compare the approaches within the context of a fixed ΛCDM model and demonstrate that our method gives conservative estimates of the systematic uncertainty that nevertheless have little impact on the final posteriors obtained from DESI data.

79 ASTRONOMY AND ASTROPHYSICS↗

Nucleon momentum fraction, helicity and transversity from 2+1-flavor lattice QCD

A detailed analysis of the systematic uncertainties in the calculation of the isovector momentum fraction, ( x ) u - d , helicity moment, ( x ) Δ u - Δ d , and the transversity moment, ( x ) δu - δd , of the nucleon is presented using high-statistics data on seven ensembles of gauge configurations generated by the JLab/W&M/LANL/MIT collaborations using 2 + 1-flavors of dynamical Wilson-clover quarks. The much higher statistics have facilitated better control over all systematics compared to previous lattice calculations. The least understood systematic — excited-state contamination — is quantified by studying the variation of the results as a function of different estimates of the mass gap of the first excited state, obtained from two- and three-point correlation functions, and as a function of the pion mass M π . The final results are obtained using a simultaneous fit in the lattice spacing a , pion mass M π and the finite volume parameter M π L keeping leading order corrections. The data show no significant dependence on the lattice spacing and some evidence for finite-volume corrections. Our final results, in the $ \overline{\mathrm{MS}} $ scheme at 2 GeV, are ( x ) u - d = 0 . 155(17)(20), ( x ) Δ u - Δ d = 0 . 183(14)(20) and ( x ) δu - δd = 0 . 220(18)(20), where the first error is the overall analysis uncertainty assuming excited-state contributions have been removed, and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with phenomenological global fit values.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DeepCMB: Lensing Reconstruction of the CMB with Deep Neural Networks

Next-generation cosmic microwave background (CMB) experiments will have lower noise and therefore increased sensitivity, enabling improved constraints on fundamental physics parameters such as the sum of neutrino masses and the tensor-to-scalar ratio r. Achieving competitive constraints on these parameters requires high signal-to-noise extraction of the projected gravitational potential from the CMB maps. Standard methods for reconstructing the lensing potential employ the quadratic estimator (QE). However, the QE performs suboptimally at the low noise levels expected in upcoming experiments. Other methods, like maximum likelihood estimators (MLE), are under active development. In this work, we demonstrate reconstruction of the CMB lensing potential with deep convolutional neural networks (CNN) - ie, a ResUNet. The network is trained and tested on simulated data, and otherwise has no physical parametrization related to the physical processes of the CMB and gravitati onal lensing. We show that, over a wide range of angular scales, ResUNets recover the input gravitational potential with a higher signal-to-noise ratio than the QE method, reaching levels comparable to analytic approximations of MLE methods. We demonstrate that the network outputs quantifiably different lensing maps when given input CMB maps generated with different cosmologies. We also show we can use the reconstructed lensing map for cosmological parameter estimation. This application of CNN provides a few innovations at the intersection of cosmology and machine learning. First, while training and regressing on images, we predict a continuous-variable field rather than discrete classes. Second, we are able to establish uncertainty measures for the network output that are analogous to standard methods. We expect this approach to excel in capturing hard-to-model non-Gaussian astrophysical foreground and noise contributions.

79 ASTRONOMY AND ASTROPHYSICS↗

The GMRT archive atomic gas survey – II. Mass modelling and dark matter halo properties across late-type spirals

Studying the kinematics and mass modelling of galaxies from H i 21 cm data provides valuable insights into the properties of both the baryonic components and the dark matter halo in nearby galaxies. Despite many observational studies, mass modelling of galaxies remains challenging due to different limitations. For example, most of the previous studies involving mass modelling are based on rotation curves derived from 2D velocity fields from H i or H α spectroscopic observation which are often affected by beam smearing and projection effect. However, kinematic modelling done by fitting the ‘Tilted ring model’ to 3D data cube is not affected by these issues. In this study, we present and compare 3D kinematic modelling of a pilot sample of 11 galaxies from the GMRT archive atomic gas survey (GARCIA) using two different publicly available pipelines. We model the observed H i rotation curve using 3.6-μm infrared data and SDSS r-band data for stellar contribution, H i surface density profile for gas, and Navarro–Frenk–White profile for dark matter halo; and employ the Markov chain Monte Carlo optimization method for parameter estimation. Further, to validate our analysis, we revisit important scaling relations, e.g. the M gas –M star relation, M star –M halo relation, M gas –M halo relation and Baryonic Tully–Fisher relation. Here, the scaling relations from our analysis are broadly consistent with that reported in the literature. A larger sample of galaxies from GARCIA in the near future will allow studying these scaling relations in greater details.

79 ASTRONOMY AND ASTROPHYSICS↗

Exploration of Domain Aware Machine Learning for Grid Analytics: Transfer-Learnt Energy Models to Assist Buildings Control with Sparse Field Data

Buildings are a primary consumer of energy in the United States and are also increasingly being perceived as providers of grid services such as load shifting, shedding and modulation. High fidelity models of building energy consumption are needed to set appropriate baselines for measurement and verification (M&V) of controllers designed for energy efficient operation of buildings and to enable buildings to provide grid services via. participation in demand response programs. State-of-the-art building energy modeling techniques either rely on Physics based models, or extensive instrumentation of the building envelope to gather “big” data to train machine learning based models such as deep neural networks. While Physics based models are often limited by their accuracy, it is not always feasible to gather a significant amount of field data required to train machine learning based models with sufficient accuracy. In this paper, we explore the use of transfer learning-based strategies to address unsatisfactory accuracy of models for estimating building energy consumption when available field data for training is sparse or of unacceptable quality. In particular, we transfer knowledge in the form of data and parameters, from Physics based simulation frameworks to the field to improve the model accuracy, thus resulting in a Physics-informed Machine Learning framework. We evaluated the efficacy of our approach on field data collected from six commercial buildings and our results indicate that the proposed transfer learning based models provide comparative (and in some cases better) accuracy than state-of-the-art machine learning and deep learning solutions, with just one month of field data.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Characterization and Valuation of the Uncertainty of Calibrated Parameters in Microsimulation Decision Models

We evaluated the implications of different approaches to characterize the uncertainty of calibrated parameters of microsimulation decision models (DMs) and quantified the value of such uncertainty in decision making. We calibrated the natural history model of CRC to simulated epidemiological data with different degrees of uncertainty and obtained the joint posterior distribution of the parameters using a Bayesian approach. We conducted a probabilistic sensitivity analysis (PSA) on all the model parameters with different characterizations of the uncertainty of the calibrated parameters. We estimated the value of uncertainty of the various characterizations with a value of information analysis. We conducted all analyses using high-performance computing resources running the Extreme-scale Model Exploration with Swift (EMEWS) framework. The posterior distribution had a high correlation among some parameters. The parameters of the Weibull hazard function for the age of onset of adenomas had the highest posterior correlation of -0.958. When comparing full posterior distributions and the maximum-a-posteriori estimate of the calibrated parameters, there is little difference in the spread of the distribution of the CEA outcomes with a similar expected value of perfect information (EVPI) of $\$$653 and $\$$685, respectively, at a willingness-to-pay (WTP) threshold of $\$$66,000 per quality-adjusted life year (QALY). Ignoring correlation on the calibrated parameters’ posterior distribution produced the broadest distribution of CEA outcomes and the highest EVPI of $\$$809 at the same WTP threshold. Different characterizations of the uncertainty of calibrated parameters affect the expected value of eliminating parametric uncertainty on the CEA. Ignoring inherent correlation among calibrated parameters on a PSA overestimates the value of uncertainty.

97 MATHEMATICS AND COMPUTING↗

Modeling Within-Host Dynamics of SARS-CoV-2 Infection: A Case Study in Ferrets

The pre-clinical development of antiviral agents involves experimental trials in animals and ferrets as an animal model for the study of SARS-CoV-2. Here, we used mathematical models and experimental data to characterize the within-host infection dynamics of SARS-CoV-2 in ferrets. We also performed a global sensitivity analysis of model parameters impacting the characteristics of the viral infection. We provide estimates of the viral dynamic parameters in ferrets, such as the infection rate, the virus production rate, the infectious virus proportion, the infected cell death rate, the virus clearance rate, as well as other related characteristics, including the basic reproduction number, pre-peak infectious viral growth rate, post-peak infectious viral decay rate, pre-peak infectious viral doubling time, post-peak infectious virus half-life, and the target cell loss in the respiratory tract. These parameters and indices are not significantly different between animals infected with viral strains isolated from the environment and isolated from human hosts, indicating a potential for transmission from fomites. While the infection period in ferrets is relatively short, the similarity observed between our results and previous results in humans supports that ferrets can be an appropriate animal model for SARS-CoV-2 dynamics-related studies, and our estimates provide helpful information for such studies.

59 BASIC BIOLOGICAL SCIENCES↗

Tropospheric and stratospheric wildfire smoke profiling with lidar: mass, surface area, CCN, and INP retrieval

We present retrievals of tropospheric and stratospheric height profiles of particle mass, volume, surface area, and number concentrations in the case of wildfire smoke layers as well as estimates of smoke-related cloud condensation nuclei (CCN) and ice-nucleating particle (INP) concentrations from backscatter lidar measurements on the ground and in space. Conversion factors used to convert the optical measurements into microphysical properties play a central role in the data analysis, in addition to estimates of the smoke extinction-to-backscatter ratios required to obtain smoke extinction coefficients. The set of needed conversion parameters for wildfire smoke is derived from AERONET observations of major smoke events, e.g., in western Canada in August 2017, California in September 2020, and southeastern Australia in January–February 2020 as well as from AERONET long-term observations of smoke in the Amazon region, southern Africa, and Southeast Asia. The new smoke analysis scheme is applied to CALIPSO observations of tropospheric smoke plumes over the United States in September 2020 and to ground-based lidar observation in Punta Arenas, in southern Chile, in aged Australian smoke layers in the stratosphere in January 2020. These case studies show the potential of spaceborne and ground-based lidars to document large-scale and long-lasting wildfire smoke events in detail and thus to provide valuable information for climate, cloud, and air chemistry modeling efforts performed to investigate the role of wildfire smoke in the atmospheric system.

54 ENVIRONMENTAL SCIENCES↗

Feature Selection for Monitoring Erosive Cavitation on a Hydroturbine

This paper presents a method for comparing and evaluating cavitation detection features - the first step towards estimating remaining useful life (RUL) of hydroturbine runners that are impacted by erosive cavitation. The method can be used to quickly compare features created from cavitation survey data collected on any type of hydroturbine, sensor type, sensor location, and cavitation sensitivity parameter (CSP). Although manual evaluation and knowledge of hydroturbine cavitation is still required for our feature selection method, the use of principal component analysis greatly reduces the number of plots that require evaluation. We present a case study based on a cavitation survey data collected on a Francis hydroturbine located at a hydroelectric plant and demonstrate the selection of the most advantageous sensor type, sensor location, and CSP to use on this hydroturbine for long-term monitoring of erosive cavitation. Our method provides hydroturbine operators and researchers with a clear and effective means to determine preferred sensors, sensor placements, and CSPs while also laying the groundwork for determining RUL in the future.

13 HYDRO ENERGY↗