Search NASA⌕ Search

SEARCH · Search NASA

Results for “likelihood”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

The DESI DR1 peculiar velocity survey: growth rate measurements from the maximum likelihood fields method

We present the constraint on the growth rate of structure from the combination of DESI DR1 BGS sample, Fundamental Plane, and Tully-Fisher peculiar velocity catalogues using the maximum likelihood fields method. The combined catalogue contains 415,523 galaxy redshifts and 76,616 peculiar velocity measurements. To handle the large amount of data in the DESI DR1 peculiar velocity catalogue, we significantly improve the computational efficiency by rewriting the algorithm with JAX. After removing outliers and Tully-Fisher galaxies that are affected by systematics, we find fσ 8 = 0.483 -0.043 +0.080 (stat) ± 0.018(sys), consistent within 1σ with the power spectrum and correlation function analysis using the same dataset. Combining all three measurements with appropriate correlations, the consensus measurement is fσ 8 (z eff = 0.07) = 0.450±0.055, consistent with Planck +ΛCDM cosmology (fσ 8 = 0.449±0.008). Combining with the high redshift growth rate of structure measurements from DESI ShapeFit, the constraint on the growth index is γ = 0.58±0.11, consistent with GR.

cosmic flows↗
DOE OSTI2024

Dark energy survey year 3 results: likelihood-free, simulation-based w CDM inference with neural compression of weak-lensing map statistics

We present simulation-based cosmological wcold dark matter (wCDM) inference using dark energy survey year 3 weak-lensing maps, via neural data compression of weak-lensing map summary statistics: power spectra, peak counts, and direct map-level compression/inference with convolutional neural networks (CNN). Using simulation-based inference, also known as likelihood-free or implicit inference, we use forward-modelled mock data to estimate posterior probability distributions of unknown parameters. This approach allows all statistical assumptions and uncertainties to be propagated through the forward-modelled mock data; these include sky masks, non-Gaussian shape noise, shape measurement bias, source galaxy clustering, photometric redshift uncertainty, intrinsic galaxy alignments, non-Gaussian density fields, neutrinos, and non-linear summary statistics. We include a series of tests to validate our inference results. This paper also describes the Gower Street simulation suite: 791 full-sky pkdgrav3 dark matter simulations, with cosmological model parameters sampled with a mixed active-learning strategy, from which we construct over 3000 mock dark energy survey lensing data sets. For wCDM inference, for which we allow –1 < w < –$\frac{1}{3}$⁠, our most constraining result uses power spectra combined with map-level (CNN) inference. Using gravitational lensing data only, this map-level combination gives Ω m = 0.283$^{+0.020}_{–0.027}$⁠, S 8 = 0.804$^{+0.025}_{–0.017⁠}$, and w < –0.80 (with a 68 per cent credible interval); compared to the power spectrum inference, this is more than a factor of two improvement in dark energy parameter (Ω⁠ DE , w⁠) precision.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A Latent-Variable Formulation of the Poisson Canonical Polyadic Tensor Model: Maximum Likelihood Estimation and Fisher Information

We establish parameter inference for the Poisson canonical polyadic (PCP) tensor model through a latent-variable formulation. Our approach exploits the observation that any random PCP tensor can be derived by marginalizing an unobservable random tensor of one dimension larger. The loglikelihood of this larger dimensional tensor, referred to as the “complete” loglikelihood, is comprised of multiple rank one PCP loglikelihoods. Using this methodology, we first derive maximum likelihood estimators for the PCP model and demonstrate that several existing algorithms for fitting non-negative matrix and tensor factorizations are Expectation-Maximization algorithms. Next, we derive the observed and expected Fisher information matrices for the PCP model. The Fisher information provides us crucial insights into the well-posedness of the tensor model, such as the role that tensor rank plays in identifiability and indeterminacy. For the special case of rank one PCP models, we demonstrate that these results are greatly simplified.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Likelihood-Based Particle Identification in SBND

Accurate particle identification is crucial in any high-energy physics experiment, allowing scientists to understand the unique interactions and mechanisms at play in a detector. In this project, I develop and study a likelihood-based particle identification (PID) algorithm for the Short-Baseline Near Detector, which offers a more physically motivated strategy for PID.

Vanderwaal, Sophia [U. Alabama, Huntsville]↗
DOE OSTI2026

Red Noise–based False Alarm Thresholds for Astrophysical Periodograms via Whittle’s Approximation to the Likelihood

Astronomers who search for periodic signals using Lomb–Scargle periodograms rely on false alarm level (FAL) estimates to identify statistically significant peaks. Although FALs are often calculated from white noise models, many astronomical time series suffer from red noise. Prewhitening is a statistical technique in which a continuum model is subtracted from the log power spectrum estimate, after which the observer can proceed with a white-noise treatment. Here we present a prewhitening-based method of calculating frequency-dependent FALs. We fit power laws and autoregressive models of order 1 to each Lomb–Scargle periodogram by minimizing the Whittle approximation to the negative log-likelihood (NLL), then calculate FALs based on the best-fit model power spectrum. Our technique is a novel extension of the Whittle NLL to datasets with uneven time sampling. We demonstrate FAL calculations using observations of α Cen B, GJ 581, HD 192310, synthetic data from the radial velocity (RV) fitting challenge, and Kepler observations of a differential rotator. The Kepler data analysis shows that only true rotation signals are detected by red noise FALs, while white noise FALs suggest all spurious peaks in the low-frequency range are significant. A high-frequency sinusoid injected into α Cen B logR$'$ HK observations exceeds the 1% red noise FAL despite having only 8.9% of the power of the dominant rotation signal. In a periodogram of HD 192310 RVs, peaks associated with differential rotation and planets are detected against the 5% red noise FAL without iterative model fitting or subtraction. The software for calculating red noise–based FALs is available on GitHub.

Astrostatistics (1882)↗
DOE OSTI2025

Forecasting Dark Matter Subhalo Constraints from Stellar Streams using Implicit Likelihood Inference

The evidence for dark matter (DM) remains compelling, although attempts to understand its particle nature remain inconclusive. One promising method to study DM is detecting DM subhalos through their gravitational interactions with stellar streams. In this study, we apply Neural Posterior Estimation (NPE) to constrain subhalo interaction parameters, including mass, scale radius, velocity, and encounter geometry, from stellar stream kinematics. We generate particle spray simulations based on the Lagrange Cloud stripping technique, focusing on the ATLAS-Aliqa Uma stream as a test case. We train multiple NPE models across multiple observational scenarios, quantifying how kinematic completeness affects inference and forecasting constraints from upcoming surveys including LSST, 4MOST, and 10-year Gaia data. Our results demonstrate that NPE can produce accurate and well-calibrated posteriors. In the idealized case with full 6D coordinates, we achieve subhalo mass uncertainties of 15-20% for a $10^7 \, \mathrm{M_\odot}$ subhalo, with 5D coordinates (excluding radial velocities) achieving similar performance. Under realistic observational conditions, mass uncertainties range from 50% (present-day) to 20-40% (future scenarios), with comparable performance between the photometric-only LSST sample and a smaller sample that includes Gaia proper motions and 4MOST radial velocities. Most notably, we find that velocity bimodality emerges when phase space is poorly sampled, whether due to missing kinematic information or limited stellar tracers. Combining large photometric samples with targeted spectroscopic follow-up can effectively resolves this degeneracy. These results demonstrate the power of implicit likelihood inference for optimizing stellar stream observational strategies and forecasting DM subhalo constraints from upcoming surveys.

Nguyen, Tri [Northwestern U. (main); SkAI, Chicago↗
DOE OSTI2025

Compounding effects of climate change and WUI expansion quadruple the likelihood of extreme-impact wildfires in California

Previous research has examined individual factors contributing to wildfire risk, but the compounding effects of these factors remain underexplored. Here, we introduce the “Integrated Human-centric Wildfire Risk Index (IHWRI)” to quantify the compounding effects of fire-weather intensification and anthropogenic factors—including ignitions and human settlement into wildland—on wildfire risk. While climatic trends increased the frequency of high-risk fire-weather by 2.5-fold, the combination of this trend with wildland-urban interface expansion led to a 4.1-fold increase in the frequency of conditions conducive to extreme-impact wildfires from 1990 to 2022 across California. More than three-quarters of extreme-impact wildfires—defined as the top 20 largest, most destructive, or deadliest events on record—originated within 1 km from the wildland-urban interface. The deadliest and most destructive wildfires—90% of which were human-caused—primarily occurred in the fall, while the largest wildfires—56% of which were human-caused—mostly took place in the summer. By integrating human activity and climate change impacts, we provide a holistic understanding of human-centric wildfire risk, crucial for policy development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The Atacama Cosmology Telescope: DR6 power spectra, likelihoods and ΛCDM parameters

We present power spectra of the cosmic microwave background (CMB) anisotropy in temperature and polarization, measured from the Data Release 6 maps made from Atacama Cosmology Telescope (ACT) data. These cover 19,000 deg 2 of sky in bands centered at 98, 150 and 220 GHz, with white noise levels three times lower than Planck in polarization. We find that the ACT angular power spectra estimated over 10,000 deg 2 , and measured to arcminute scales in TT, TE and EE, are well fit by the sum of CMB and foregrounds, where the CMB spectra are described by the ΛCDM model. Combining ACT with larger-scale Planck data, the joint P-ACT dataset provides tight limits on the ingredients, expansion rate, and initial conditions of the universe. We find similar constraining power, and consistent results, from either the Planck power spectra or from ACT combined with WMAP data, as well as from either temperature or polarization in the joint P-ACT dataset. When combined with CMB lensing from ACT and Planck, and baryon acoustic oscillation data from the Dark Energy Spectroscopic Instrument (DESI DR1), we measure a baryon density of Ω b h 2 = 0.0226 ± 0.0001, a cold dark matter density of Ω c h 2 = 0.118 ± 0.001, a Hubble constant of H 0 = 68.22 ± 0.36 km/s/Mpc, a spectral index of n s = 0.974 ± 0.003, and an amplitude of density fluctuations of σ 8 = 0.813 ± 0.005. Including the DESI DR2 data tightens the Hubble constant to H 0 = 68.43 ± 0.27 km/s/Mpc; ΛCDM parameters agree between the P-ACT and DESI DR2 data at the 1.6σ level. We find no evidence for excess lensing in the power spectrum, and no departure from spatial flatness. The contribution from Sunyaev-Zel'dovich (SZ) anisotropy is detected at high significance; we find evidence for a tilt with suppressed small-scale power compared to our baseline SZ template spectrum, consistent with hydrodynamical simulations with feedback.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Development of systematic uncertainty-aware neural network trainings for binned-likelihood analyses at the LHC

We propose a neural network training method capable of accounting for the effects of systematic variations of the data model in the training process and describe its extension towards neural network multiclass classification. The procedure is evaluated on the realistic case of the measurement of Higgs boson production via gluon fusion and vector boson fusion in the τ τ decay channel at the CMS experiment. The neural network output functions are used to infer the signal strengths for inclusive production of Higgs bosons as well as for their production via gluon fusion and vector boson fusion. We observe improvements of 12 and 16% in the uncertainty in the signal strengths for gluon and vector-boson fusion, respectively, compared with a conventional neural network training based on cross-entropy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Discriminative versus generative approaches to simulation-based inference

Most of the fundamental, emergent, and phenomenological parameters of particle and nuclear physics are determined through parametric template fits. Simulations are used to populate histograms which are then matched to data. This approach is inherently lossy, since histograms are binned and low-dimensional. Deep learning has enabled unbinned and high-dimensional parameter estimation through neural likelihood(-ratio) estimation. We compare two approaches for neural simulation-based inference (NSBI): one based on discriminative learning (classification) and one based on generative modeling. These two approaches are directly evaluated on the same datasets, with a similar level of hyperparameter optimization in both cases. In addition to a Gaussian dataset, we study NSBI using a Higgs boson dataset from the FAIR Universe Challenge. We find that both the direct likelihood and likelihood ratio estimation are able to effectively extract parameters with reasonable uncertainties. For the numerical examples and within the set of hyperparameters studied, we found that the likelihood ratio method is more accurate and/or precise. Both methods have a significant spread from the network training and would require ensembling or other mitigation strategies in practice.

high energy physics↗
DOE OSTI2025

Persistent Sampling: Enhancing the Efficiency of Sequential Monte Carlo

Sequential Monte Carlo (SMC) samplers are powerful tools for Bayesian inference but suffer from high computational costs due to their reliance on large particle ensembles for accurate estimates. We introduce persistent sampling (PS), an extension of SMC that systematically retains and reuses particles from all prior iterations to construct a growing, weighted ensemble. By leveraging multiple importance sampling and resampling from a mixture of historical distributions, PS mitigates the need for excessively large particle counts, directly addressing key limitations of SMC such as particle impoverishment and mode collapse. Crucially, PS achieves this without additional likelihood evaluations-weights for persistent particles are computed using cached likelihood values. This framework not only yields more accurate posterior approximations but also produces marginal likelihood estimates with significantly lower variance, enhancing reliability in model comparison. Furthermore, the persistent ensemble enables efficient adaptation of transition kernels by leveraging a larger, decorrelated particle pool. Experiments on high-dimensional Gaussian mixtures, hierarchical models, and non-convex targets demonstrate that PS consistently outperforms standard SMC and related variants, including recycled and waste-free SMC, achieving substantial reductions in mean squared error for posterior expectations and evidence estimates, all at reduced computational cost. PS thus establishes itself as a robust, scalable, and efficient alternative for complex Bayesian inference tasks.

Karamanis, Minas↗
DOE OSTI2026

Distinguishing Single and Linked Ruptures in the Laboratory and Nature

Earthquakes can grow either monotonically from a single, stressed patch or through linking multiple stressed regions. The distinction has implications for magnitude predictability with single ruptures requiring knowledge of the local stress state, while linked ruptures require knowing the global stress and energy distribution. Here, we use a laboratory fault that allows direct observation of slip to determine which fault conditions promote linked ruptures, and how to observationally evaluate their likelihood. Higher concentrations of normal stress due to increased normal force, applied stress asperities, or larger heterogeneity between the samples, all lead to significant increases in linked rupture likelihood. The mean radiated energy enhancement factor (REEF) of large events is an excellent proxy for linked event likelihood, and a combination of REEF and duration can identify the rupture style of some but not all events. The results imply that globally observed variations of REEF can be interpreted as variations in stress concentration.

earthquakes↗
DOE OSTI2024

Feasibility of peak temperature targets in light of institutional constraints

Despite faster-than-expected progress in clean energy technology deployment, global annual CO 2 emissions have increased from 2020 to 2023. The feasibility of limiting warming to 1.5 °C is therefore questioned. Here we present a model intercomparison study that accounts for emissions trends until 2023 and compares cost-effective scenarios to alternative scenarios with institutional, geophysical and technological feasibility constraints and enablers informed by previous literature. Our results show that the most ambitious mitigation trajectories with updated climate information still manage to limit peak warming to below 1.6 °C (‘low overshoot’) with around 50% likelihood. However, feasibility constraints, especially in the institutional dimension, decrease this maximum likelihood considerably to 5–45%. Accelerated energy demand transformation can reduce costs for staying below 2 °C but have only a limited impact on further increasing the likelihood of limiting warming to 1.6 °C. Our study helps to establish a new benchmark of mitigation scenarios that goes beyond the dominant cost-effective scenario design.

54 ENVIRONMENTAL SCIENCES↗