Search NASA⌕ Search

SEARCH · Search NASA

Results for “parameter inference”

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 235 records · Page 13

Modular Spectral Inference Framework Applied to Young Stars and Brown Dwarfs

In practice, synthetic spectral models are imperfect, causing inaccurate estimates of stellar parameters. Using forward modeling and statistical inference, we derive accurate stellar parameters for a given observed spectrum by emulating a grid of precomputed spectra to track uncertainties. Spectral inference as applied to brown dwarfs re: Synthetic spectral models (Marley et al 1996 and 2014) via the newest grid spans a massive multi-dimensional grid applied to IGRINS spectra, improving atmospheric models for JWST. When applied to young stars(10Myr) with large starpots, they can be measured spectroscopically, especially in the near-IR with IGRINS.

Validation of the DESI DR2 Ly$α$ forest full-shape analysis

We present the validation of the Dark Energy Spectroscopic Instrument (DESI) Data Release 2 (DR2) Lyman-$α$ (Ly$α$) forest full-shape analysis. This analysis combines three-dimensional Ly$α$ forest auto-correlations and cross-correlations with quasars to extract information from both the baryon acoustic oscillation (BAO) feature and the broadband clustering signal, with primary emphasis on the Alcock-Paczynski (AP) measurement. Compared to the DESI DR1 analysis, the DR2 validation uses substantially larger and more realistic mock datasets, including CoLoRe 2LPT and AbacusSummit Ly$α$ forest simulations. The modeling framework is also improved through analytic marginalization over small scales ($<10$$h^{-1}$Mpc) and the impact of ultraviolet background fluctuations. The validation program was completed prior to unblinding and defines quantitative requirements for the cosmological parameters of interest, which are evaluated using hundreds of mock realizations. We further test the analysis through independent fits to the auto- and cross-correlations, multiple catalog splits, and a broad suite of analysis and modeling variations applied to both mocks and blinded observational data. We find that the BAO and AP parameters satisfy all validation requirements and remain stable across all tests. In contrast, mock studies reveal a significant bias in the inferred growth-rate parameter $fσ_8$, leading us to exclude this measurement from the final analysis. The consistency across mocks, data splits, and robustness tests demonstrates that the DR2 Ly$α$ full-shape analysis provides a reliable and substantially improved broadband AP measurement over previous Ly$α$ forest studies.

Herbold, M. [Chicago U., KICP; Ohio State U.] (ORC↗

Transition regions in solar system and astrophysical plasmas

A brief review is presented of basic particle and field characteristics of plasmas observed within the solar system, especially near transition regions, and their parameter ranges are compared with those inferred for stellar winds and the interstellar medium. Parameter ranges for solar system and astrophysical plasmas are found to have considerable overlap. In addition, astrophysics provides unique, global perspectives of large-scale systems, whereas solar-system space physics provides for direct quantitative testing of physical processes. Astrophysics and solar-system space physics studies thus have complementary and synergistic roles.

Eastman, Timothy E.↗

Use of aerosol microphysical measurements to model IR backscatter in support of GLOBE

Data on the concentration and composition of free tropospheric aerosol over the Pacific Ocean, collected during the GAMETAG program in 1977-1978 (Davis, 1980 and Patterson et al., 1980) are used to model values of aerosol optical extinction coefficients (sigma) at two wavelengths (0.55 and 1 micron) and values of volume backscatter coefficients (beta) at four wavelengths (1 micron, 9.11 microns, 9.25 microns, and 10.6 microns) and to investigate the relationship between these parameters. The mass concentrations inferred from the GAMETAG measurements with optical particle spectrometers agreed with the results of simultaneous chemical measurements. The study of the relationships among the optical parameters indicates that visible and near-visible values of beta and sigma may be useful in predicting 9.11- and 10.6-micron backscatter.

Patterson, Edward M.↗

Sequential spectral line analysis for accurate density and temperature diagnosis of laboratory opacity measurements

The accuracy of iron opacity calculated in stellar interiors has been questioned since the discovery of the “solar problem” and the discrepancies between the measured and modeled iron opacity reported in 2015. Experimental opacity benchmarks require accurate temperature and density measurements, which were inferred by analyzing tracer magnesium spectra in those experiments. Could the observed discrepancy be explained by insufficient accuracy in the inferred temperature, density, and their uncertainties? Previous analyses may have yielded biased results due to three limitations: (1) simultaneous multi-line fitting, (2) approximations in line-shape models, and (3) exclusion of certain spectral lines due to insufficient background characterization. Notably, the first issue is a common concern for many inversion methods, including Bayesian inferences. We present a refined analysis method that overcomes these limitations, applied to three categories of iron opacity experiments (Anchor 1, 2, and 3). In particular, the sequential fitting method yields unbiased results with more realistic uncertainties by accounting for line inconsistencies in the parameter uncertainties. The average electron temperature and density values are 162 ± 6 eV and (7.0 ± 1.9) × 10 21 cm −3 for six Anchor 1 experiments, 189 ± 7 eV and (3.4 ± 0.3) × 10 22 cm −3 for 21 Anchor 2 experiments, and 201 ± 6 eV and (4.8 ± 1.1) × 10 22 cm −3 for nine Anchor 3 experiments. These results show ∼4% temperature and ∼20% density reproducibility over a decade, which also aligns with the inferred parameter uncertainties. In conclusion, the resulting temperature and density uncertainties lead to a quasi-continuum iron opacity variation of ±4%–7% for wavelengths below 9.5 Å, which is insufficient to explain the significant model-data discrepancies reported in 2015.

Absorption spectroscopy↗

A non-intrusive framework using acoustic signals and deep learning for boiling diagnostics in visual-limited environments

Accurate monitoring of boiling heat transfer is critical for safeguarding high-power systems operating in environments where conventional optical diagnostics are hindered by radiation fields or restricted visual accessibility. This study presents a non-intrusive framework that integrates hydroacoustic sensing with deep learning to infer near-wall boiling characteristics and enable predictive thermal assessment without visual access. In a prototypical subcooled flow-boiling facility representative of the Isotope Production Facility (IPF) at Los Alamos, hydrophones capture boiling-induced acoustic emissions that are transformed into background-removed Short-Time Fourier Transform (STFT) spectrograms. A convolutional neural network (CNN) then regresses heat flux, wall superheat, and key bubble parameters directly from these spectrograms. The CNN achieved predictive accuracy under nominal conditions and demonstrated robustness and generalization under acoustic noise for Signal-to-Noise Ratios (SNRs) down to approximately 0 dB. When integrated into an ANSYS CFX wall-boiling model, the acoustically inferred parameters reproduced boiling curve and critical heat flux (CHF) values consistent with image-based benchmarks. Furthermore, the model retained reliable performance under moderate variations in bulk temperature, flow rate, and hydrophone placement, confirming its generalizability across practical boundary conditions. These results demonstrate the feasibility of hydroacoustic-based deep learning as a viable path toward real-time, radiation-tolerant boiling diagnostics and predictive thermal safety assessment in inaccessible systems such as the IPF.

42 ENGINEERING↗

Neutrinos from SN 1987A - Implications for cooling of the nascent neutron star and the mass of the electron antineutrino

Data on neutrinos from SN 1987A are compared here with parameterized models of the neutrino emission using a consistent and straightforward statistical methodology. The empirically measured detector background spectra are included in the analysis, and the data are compared with a much wider variety of neutrino emission models than was explored previously. It is shown that the inferred neutrino emission model parameters are strongly correlated. The analysis confirms that simple models of the neutrino cooling of the nascent neutron star formed by the SN adequately explain the data. The inferred radius and binding energy of the neutron star are in excellent agreement with model calculations based on a wide range of equations of state. The results also raise the upper limit of the electron antineutrino rest mass to roughly 25 eV at the 95 percent confidence level, roughly 1.5-5 times higher than found previously.

Loredo, Thomas J.↗

TMI Rain Rate Estimation Over Land and Ocean Utilizing Convective and Stratiform Discrimination

Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) radiometer brightness temperature data in the 85 GHz channel (T85) reveal distinct local minima in a regional map containing a Mesoscale Convective System (MCS). This is because of relatively small footprint size (approximately 5.5 km) and strong extinction properties in this channel of the TMI. A map of rain rate for that region, deduced from simultaneous measurements made by the Precipitation Radar (PR) on board the TRMM satellite, reveals that these T85 minima, produced by scattering, correspond to local PR rain maxima. Utilizing the PR rain rate map as a guide, we infer from TMI data the presence of three different kinds of thunderstorms or Cbs. They are young, mature, and decaying Cbs that have a scale of about 20 km on the average. Two parameters enable us to infer these three kinds of Cbs objectively: a) the magnitude of scattering depression deduced from local T85 minima and b) the mean horizontal gradient of T85 around such minima. Knowing the category of a given Cb, we can estimate the rain rate associated with it. Such estimation is done with the help of relationships linking T85 minimum to rain rate in each Cb type. Similarly, a weak background rain rate in all the areas where T85 is less than 260 K is deduced with another relationship linking T85 to rain rate. In our rain retrieval model, this background rain constitutes the stratiform rain where the Cbs are absent. Initially, these relationships are optimized or tuned utilizing the PR and TMI data of a few MCS events. After such tuning, the model is applied to independent MCS cases. The areal distribution of light (1-10 mm/hr), moderate (10-20 mm/hr), and intense (> 20 mm/hr) rain rates are retrieved satisfactorally. Accuracy in the estimates of the light, moderate and intense rain areas and the mean rain rates associated with such areas in these independent MCS cases is on the average about 15%. Taking advantage of this ability of our retrieval method, one could derive the latent heat input into the atmosphere over the 760 km wide swath of the TMI radiometer in the tropics.

Prabhakara, C.↗

From pole parameters to line shapes and branching ratios

Abstract Resonances are uniquely characterized by their complex pole locations and the corresponding residues. In practice, however, resonances are typically identified experimentally as structures in invariant mass distributions, with branching fractions of resonances determined as ratios of count rates. To make contact between these quantities it is necessary to connect line shapes and resonance parameters. In this work we propose such a connection and illustrate the formalism with detailed studies of the $$\rho (770)$$ ρ ( 770 ) and $$f_0(500)$$ f 0 ( 500 ) resonances. Based on the line shapes inferred from the resonance parameters along these lines, expressions for partial widths and branching ratios are derived and compared to other approaches in the literature.

Physics↗

A design methodology for nonlinear systems containing parameter uncertainty

In the present design methodology for nonlinear systems containing parameter uncertainty, a generalized sensitivity analysis is incorporated which employs parameter space sampling and statistical inference. For the case of a system with j adjustable and k nonadjustable parameters, this methodology (which includes an adaptive random search strategy) is used to determine the combination of j adjustable parameter values which maximize the probability of those performance indices which simultaneously satisfy design criteria in spite of the uncertainty due to k nonadjustable parameters.

Young, G. E.↗

Quantiles, parametric-select density estimation, and bi-information parameter estimators

A quantile-based approach to statistical analysis and probability modeling of data is presented which formulates statistical inference problems as functional inference problems in which the parameters to be estimated are density functions. Density estimators can be non-parametric (computed independently of model identified) or parametric-select (approximated by finite parametric models that can provide standard models whose fit can be tested). Exponential models and autoregressive models are approximating densities which can be justified as maximum entropy for respectively the entropy of a probability density and the entropy of a quantile density. Applications of these ideas are outlined to the problems of modeling: (1) univariate data; (2) bivariate data and tests for independence; and (3) two samples and likelihood ratios. It is proposed that bi-information estimation of a density function can be developed by analogy to the problem of identification of regression models.

Parzen, E.↗

Extracting the Breakout Distance from the ECOT Trajectories: Gaussian Process Regression Approach

Enhanced Corner Turning (ECOT) experiments provide an important metric of performance of high explosive (HE) formulations. The breakout distance is a single scalar value that characterizes the corner turning efficiency of an HE. Extracting the breakout distance from the raw ECOT results, whether experimental or simulated, is a conceptually straightforward procedure which, however, is non-unique, especially in the presence of noise. More specifically, this procedure involves numerical smoothing and selecting particular values for parameters of this smoothing introduces human bias. In this work, we propose to use the Gaussian process regression to analyze ECOT results. This analysis involves the effective smoothing of the data, thus allowing for accurate extraction of the breakout distance. Most importantly, the parameters of this smoothing can be inferred from the ECOT data itself, rendering the approach effectively parameter-free and thus diminishing the human bias. An additional benefit of the Gaussian process regression, being a statistical inference method, is that not just the value of the breakout distance, but also its confidence interval can be extracted from the data. This report introduces the Gaussian process regression, as applied to ECOT, and demonstrates its usefulness by extracting the breakout distances for a selection of experimental and simulated data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

New approaches to Bayesian uncertainty quantification for Nuclear Science (Final Technical Report)

Inverse problems play a central role in experimentation and theory/data comparisons for many areas of modern Nuclear Physics (NP) and High-Energy Physics (HEP). Bayes’s Theorem is a powerful tool for solving Inverse Problems, providing conceptually transparent and unbiased constraints on theoretical parameters and their uncertainties (“Bayesian Inference”) and enabling the quantification of agreement or tension between models and data. However, analyses based on Bayesian Inference are often challenging for NP and HEP applications, either because of the large number of parameters in the problem, the high computational cost, or both. We propose a multi-institutional collaboration to develop and deploy novel Bayesian analysis tools that advance the scientific scope of a broad range of current and future NP experiments. This project brings together NP domain scientists working on several high-profile NP projects for which new, high-performance Bayesian Uncertainty Quantification (“Bayesian UQ”) methods are essential to carry out the science, and data scientists who are developing state-of-the-art methods applicable to these problems. The NP projects in this proposal comprise measurements of the mass and fundamental nature of the neutrino; study of the Quark-Gluon Plasma that filled the early universe; and mapping of natural and anthropogenic radiation environments. While these NP projects have very different scientific goals, with datasets and analysis approaches that differ significantly, they share common requirements for improving computationally intensive Bayesian analyses using advanced Machine Learning algorithms and will benefit strongly from a coherent effort to develop general solutions. This proposal brings together these projects and forefront ML-based data science algorithms to develop such general solutions. The methods developed in this project will also be more widely applicable, thereby advancing science in the larger Nuclear Physics portfolio.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The SRG/eROSITA All-Sky Survey: Constraints on the structure growth from cluster number counts

Recent advancements in methods used in wide-area surveys have demonstrated the reliability of the number density of galaxy clusters as a viable tool for precision cosmology. Beyond testing the current cosmological paradigm, cluster number counts can also be used to investigate the discrepancies currently affecting cosmological measurements. In particular, cosmological studies based on cosmic shear and other large-scale structure probes routinely find a value for the amplitude of the fluctuations in the universe S​ 8 = σ 8 (Ω m /0.3) 0.5 smaller than the one inferred from the primary cosmic microwave background. In this work, we investigate this tension by measuring structure evolution across cosmic time as probed by the number counts of massive halos with the first SRG/eROSITA All-Sky Survey cluster catalog in the western Galactic hemisphere, complemented with the overlapping Dark Energy Survey Year-3, Kilo-Degree Survey, and Hyper Suprime-Cam data for weak lensing mass calibration, by implementing two different parameterizations and a model-agnostic method. In the first model, we measured the cosmic linear growth index as γ = 1.19 ± 0.21, which is in tension with the standard value of γ = 0.55 but in good statistical agreement with other large-scale structure probes. The second model is a phenomenological scenario in which we rescale the linear matter power spectrum at low redshift to investigate a potential reduction of structure formation, and it provided similar results. Finally, in a third strategy, we considered a standard ΛCDM cosmology, but we separated the cluster catalog into five redshift bins, measuring the cosmological parameters in each and inferring the evolution of the structure formation, finding hints of a reduction. Interestingly, the S ​8 value inferred from the number counts of the cluster eRASS1 when we add a degree of freedom to the matter power spectrum recovers the value inferred by cosmic shear studies. The observed reduction in the growth rate or systematic uncertainties associated with various measurements may account for the discrepancy in the S 8 ​ values suggested between cosmic shear probes and eROSITA cluster number counts and Planck CMB measurements.

cosmological parameters↗

Informed total-error-minimizing priors: Interpretable cosmological parameter constraints despite complex nuisance effects

While Bayesian inference techniques are standard in cosmological analyses, it is common to interpret resulting parameter constraints with a frequentist intuition. This intuition can fail, for example, when marginalizing high-dimensional parameter spaces onto subsets of parameters, because of what has come to be known as projection effects or prior volume effects. We present the method of informed total-error-minimizing (ITEM) priors to address this problem. An ITEM prior is a prior distribution on a set of nuisance parameters, such as those describing astrophysical or calibration systematics, intended to enforce the validity of a frequentist interpretation of the posterior constraints derived for a set of target parameters (e.g., cosmological parameters). Our method works as follows. For a set of plausible nuisance realizations, we generate target parameter posteriors using several different candidate priors for the nuisance parameters. We reject candidate priors that do not accomplish the minimum requirements of bias (of point estimates) and coverage (of confidence regions among a set of noisy realizations of the data) for the target parameters on one or more of the plausible nuisance realizations. Of the priors that survive this cut, we select the ITEM prior as the one that minimizes the total error of the marginalized posteriors of the target parameters. As a proof of concept, we applied our method to the density split statistics measured in Dark Energy Survey Year 1 data. We demonstrate that the ITEM priors substantially reduce prior volume effects that otherwise arise and that they allow for sharpened yet robust constraints on the parameters of interest.

79 ASTRONOMY AND ASTROPHYSICS↗

A simulation survey of galaxy interactions

Many carefully selected samples of interacting galaxies have been observed extensively in attempts to clarify whether interaction produces activity in galaxies. Because the sample members represent a wide range of encounter parameters and times, one can then study whether there are correlations between observable encounter features and, for example, Seyfert activity. On the other hand, in theoretical studies, simulations typically deal with either time-consuming detailed modelling of single galaxy pairs or tracing a few model encounters over time. The authors extend the observational survey approach by combining it with a simulation survey. The authors are conducting a survey of model encounters, covering the most important encounter parameters over a wide range. Some parameters, such as companion structure and initial velocity, are demonstratably less important and can be ignored in a first pass. The parameter range must be richly enough sampled so that the authors can evaluate the uniqueness of the observable morphology and velocity structure of the resulting simulated pairs to diagnose unobservable companion orbit parameters. They are using a self-gravitating polar n-body code run on the Cray X-MP at the Alabama Supercomputer Network. For each simulation, the authors have stellar and gas distributions predicted over, typically, a billion years, along with information on gas motions within the disk and any material captured by the companion or lost to the system. Features of disturbed spiral galaxies are sensitive enough to time and encounter parameters so that a match of the simulation survey results to observations can be applied as starting points to infer unobservable orbital or system parameters in actual sample members. This should enable them to examine whether interesting observed properties (Seyfert activity, nuclear star-formation rate) are functions of unobservable dynamical properties which characterize each encounter. Any correlations (or lack of some expected ones) will provide strong clues as to how or whether these phenomena are related to interactions. Aside from its use with such observed samples, this survey should greatly speed determination of initial orbital parameters for more detailed subsequent simulations of individual systems.

Byrd, Gene G.↗

Multi-Modal Bayesian Neural Network Surrogates with Conjugate Last-Layer Estimation

As data collection and simulation capabilities advance, multi-modal learning, the task of learning from multiple modalities and sources of data, is becoming an increasingly important area of research. Surrogate models that learn from data of multiple auxiliary modalities to support the modeling of a highly expensive quantity of interest have the potential to aid outer loop applications such as optimization, inverse problems, or sensitivity analyses when multi-modal data are available. We develop two multi-modal Bayesian neural network surrogate models and leverage conditionally conjugate distributions in the last layer to estimate model parameters using stochastic variational inference (SVI). We provide a method to perform this conjugate SVI estimation in the presence of partially missing observations. Here, we demonstrate improved prediction accuracy and uncertainty quantification compared to unimodal surrogate models for both scalar and time series data.

97 MATHEMATICS AND COMPUTING↗

Synthetic modeling of soft x-ray emissivity for magnetic island analysis in LTX-β

We present a synthetic soft x-ray (SXR) forward-modeling framework that characterizes the emissivity structure of rotating magnetic islands in the Lithium Tokamak eXperiment-β across space and time. Magnetic islands associated with tearing modes produce modulations in line-integrated SXR brightness. Traditional tomographic methods struggle to resolve these structures in devices with limited sightlines, resulting in an under-determined inversion problem, or otherwise necessitate equilibrium reconstruction data, preventing their use in active control. Our approach removes this challenge by combining a fully three-dimensional ray-tracing model with a time-dependent emissivity prescription of a $m/n$ = 2/1 magnetic island. The model is validated through helical island geometry, rotation measured by magnetic diagnostics, and equilibrium constraints from PSI-Tri reconstructions. Synthetic brightness signals are generated for each photodiode sightline from a modeled emissivity profile and directly compared with experimental data from a tangential midplane SXR array. By fitting the synthetic diagnostic output to the observed brightness across time, we simultaneously infer key geometric parameters-such as island radial location, island width, and rotation frequency-without requiring full tomographic reconstruction, all with an average absolute deviation of less than 3%. This work demonstrates that forward modeling with a single tangential array can extract key island parameters in a spherical tokamak, providing a computationally efficient alternative to conventional SXR tomography and providing a pathway toward real-time magnetic-island characterization in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗