Search NASA⌕ Search

SEARCH · Search NASA

Results for “Bayesian Analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Methods for Probabilistic Uncertainty Analysis and Bayesian Analysis with Examples of Statistically Analyzing Data to Revise MMOD Risk Estimates and Compare Models

Probabilistic methods are presented for characterizing and quantifying uncertainties in models and model predictions. Techniques are given for constructing specific uncertainty distributions based on available information. Alternative techniques are given for propagating uncertainties in model inputs to obtain the uncertainty in the model result or prediction. Bayesian techniques are also described for utilizing data and information to update and revise model results and predictions. The focus is on applications with numerous specific examples given. The use of data to revise Micrometeoroid and Orbital Debris (MMOD) risk prediction models are among the examples given.

Risk↗

Bayesian analysis of inflation: Parameter estimation for single field models

Future astrophysical data sets promise to strengthen constraints on models of inflation, and extracting these constraints requires methods and tools commensurate with the quality of the data. In this paper we describe ModeCode, a new, publicly available code that computes the primordial scalar and tensor power spectra for single-field inflationary models. ModeCode solves the inflationary mode equations numerically, avoiding the slow roll approximation. It is interfaced with CAMB and CosmoMC to compute cosmic microwave background angular power spectra and perform likelihood analysis and parameter estimation. ModeCode is easily extendable to additional models of inflation, and future updates will include Bayesian model comparison. Errors from ModeCode contribute negligibly to the error budget for analyses of data from Planck or other next generation experiments. We constrain representative single-field models ({phi}{sup n} with n=2/3, 1, 2, and 4, natural inflation, and 'hilltop' inflation) using current data, and provide forecasts for Planck. From current data, we obtain weak but nontrivial limits on the post-inflationary physics, which is a significant source of uncertainty in the predictions of inflationary models, while we find that Planck will dramatically improve these constraints. In particular, Planck will link the inflationary dynamics with the post-inflationary growth of the horizon, and thus begin to probe the ''primordial dark ages'' between TeV and grand unified theory scale energies.

79 ASTRONOMY AND ASTROPHYSICS↗

Radiative properties of quantum emitters in boron nitride from excited state calculations and Bayesian analysis

Abstract Point defects in hexagonal boron nitride (hBN) have attracted growing attention as bright single-photon emitters. However, understanding of their atomic structure and radiative properties remains incomplete. Here we study the excited states and radiative lifetimes of over 20 native defects and carbon or oxygen impurities in hBN using ab initio density functional theory and GW plus Bethe-Salpeter equation calculations, generating a large data set of their emission energy, polarization and lifetime. We find a wide variability across quantum emitters, with exciton energies ranging from 0.3 to 4 eV and radiative lifetimes from ns to ms for different defect structures. Through a Bayesian statistical analysis, we identify various high-likelihood charge-neutral defect emitters, among which the native V N N B defect is predicted to possess emission energy and radiative lifetime in agreement with experiments. Our work advances the microscopic understanding of hBN single-photon emitters and introduces a computational framework to characterize and identify quantum emitters in 2D materials.

Chemistry↗

Transverse mode analysis for free-space laser beams using Bayesian analysis

Reliable and rapid assessment of the transverse mode quality of a free-space laser beam has a wide range of applications in laser development, research, and utilization. It has become even more important with recent advances in developing orbital angular momentum photon beams across a broad spectral region. In this work, a general modal analysis method for a free-space multimode laser beam has been developed based on Bayesian analysis. After transforming mode decomposition into a linear system problem, a Gaussian probabilistic model is used to find a closed-form solution. The method is found to be robust with the presence of Gaussian noise. Prior knowledge about the mode content can be incorporated into the method to improve the solution for situations when coherent disturbances or contamination are present in the laser beam. This method can be used to analyze the mode content for laser beams in different bases, such as Hermite–Gaussian (HG) modes and Laguerre–Gaussian (LG) modes. Three applications of this method are presented: a detailed modal analysis of the beam image from the incoherent intensity addition of HG modes and two examples of mode decomposition using the complex wavefront from the coherent superposition of HG and LG modes. The feasibility of this method is demonstrated using various simulation results. Based on digital images of a laser beam recorded without complex wavelength-limiting optics, in principle, this method can be used in a wide wavelength range from infrared to ultraviolet, and possibly x ray.

Liu, Peifan (ORCID:0000000273816153)↗

Bayesian analysis of rare events

In many areas of engineering and science there is an interest in predicting the probability of rare events, in particular in applications related to safety and security. Increasingly, such predictions are made through computer models of physical systems in an uncertainty quantification framework. Additionally, with advances in IT, monitoring and sensor technology, an increasing amount of data on the performance of the systems is collected. This data can be used to reduce uncertainty, improve the probability estimates and consequently enhance the management of rare events and associated risks. Bayesian analysis is the ideal method to include the data into the probabilistic model. It ensures a consistent probabilistic treatment of uncertainty, which is central in the prediction of rare events, where extrapolation from the domain of observation is common. We present a framework for performing Bayesian updating of rare event probabilities, termed BUS. It is based on a reinterpretation of the classical rejection-sampling approach to Bayesian analysis, which enables the use of established methods for estimating probabilities of rare events. By drawing upon these methods, the framework makes use of their computational efficiency. These methods include the First-Order Reliability Method (FORM), tailored importance sampling (IS) methods and Subset Simulation (SuS). In this contribution, we briefly review these methods in the context of the BUS framework and investigate their applicability to Bayesian analysis of rare events in different settings. We find that, for some applications, FORM can be highly efficient and is surprisingly accurate, enabling Bayesian analysis of rare events with just a few model evaluations. In a general setting, BUS implemented through IS and SuS is more robust and flexible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Determination of Quark-Gluon-Plasma Parameters from a Global Bayesian Analysis

The quality of data taken at RHIC and LHC as well as the success and sophistication of computational models for the description of ultra-relativistic heavy-ion collisions have advanced to a level that allows for the quantitative extraction of the transport properties of the Quark-Gluon-Plasma. However, the complexity of this task as well as the computational effort associated with it can only be overcome by developing novel methodologies: in this paper we outline such an analysis based on Bayesian Statistics and systematically compare an event-by-event heavy-ion collision model to data from the Large Hadron Collider. We simultaneously probe multiple model parameters including fundamental quark-gluon plasma properties such as the temperature-dependence of the specific shear viscosity η/s, calibrate the model to optimally reproduce experimental data, and extract quantitative constraints for all parameters simultaneously. As a result, the method is universal and easily extensible to other data and collision models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

KMT-2019-BLG-0371 and the Limits of Bayesian Analysis

We show that the perturbation at the peak of the light curve of microlensing event KMT-2019-BLG-0371 is explained by a model with a mass ratio between the host star and planet of q ∼ 0.08. Due to the short event duration (t(sub E) ∼ 6.5 days), the secondary object in this system could potentially be a massive giant planet. A Bayesian analysis shows that the system most likely consists of a host star with a mass M(sub h) = 0.09(+0.14/-0.05) M⨀ and a massive giant planet with a mass = M(sub h) = 7.70(+11.34/-3.90) M(sub Jup). However, the interpretation of the secondary as a planet (i.e., as having M(sub p) < 13M(sub Jup)) rests entirely on the Bayesian analysis. Motivated by this event, we conduct an investigation to determine which constraints meaningfully affect Bayesian analyses for microlensing events. We find that the masses inferred from such a Bayesian analysis are determined almost entirely by the measured value of θ(sub E) and are relatively insensitive to other factors such as the direction of the event (l,b), the lens–source relative proper motion μ(sub rel), or the specific Galactic model prior.

Binary lens microlensing↗

KMT-2019-BLG-0371 and the Limits of Bayesian Analysis

We show that the perturbation at the peak of the light curve of microlensing event KMT-2019-BLG-0371 is explained by a model with a mass ratio between the host star and planet of q ∼ 0.08. Due to the short event duration (t {sub E} ∼ 6.5 days), the secondary object in this system could potentially be a massive giant planet. A Bayesian analysis shows that the system most likely consists of a host star with a mass M{sub h}=0.09{sub −0.05}{sup +0.14}M{sub ⊙} and a massive giant planet with a mass M{sub p}=7.70{sub −3.90}{sup +11.34}M{sub Jup}. However, the interpretation of the secondary as a planet (i.e., as having M {sub p} < 13M {sub Jup}) rests entirely on the Bayesian analysis. Motivated by this event, we conduct an investigation to determine which constraints meaningfully affect Bayesian analyses for microlensing events. We find that the masses inferred from such a Bayesian analysis are determined almost entirely by the measured value of θ {sub E} and are relatively insensitive to other factors such as the direction of the event (ℓ, b), the lens–source relative proper motion μ {sub rel}, or the specific Galactic model prior.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ockham's razor and Bayesian analysis

'Ockham's razor', the ad hoc principle enjoining the greatest possible simplicity in theoretical explanations, is presently shown to be justifiable as a consequence of Bayesian inference; Bayesian analysis can, moreover, clarify the nature of the 'simplest' hypothesis consistent with the given data. By choosing the prior probabilities of hypotheses, it becomes possible to quantify the scientific judgment that simpler hypotheses are more likely to be correct. Bayesian analysis also shows that a hypothesis with fewer adjustable parameters intrinsically possesses an enhanced posterior probability, due to the clarity of its predictions.

Jefferys, William H.↗

A Bayesian analysis of two probability models describing thunderstorm activity at Cape Kennedy, Florida

A Bayesian analysis of the two discrete probability models, the negative binomial and the modified negative binomial distributions, which have been used to describe thunderstorm activity at Cape Kennedy, Florida, is presented. The Bayesian approach with beta prior distributions is compared to the classical approach which uses a moment method of estimation or a maximum-likelihood method. The accuracy and simplicity of the Bayesian method is demonstrated.

Williford, W. O.↗

Characterization of the initial state and QGP medium from a combined Bayesian analysis of LHC data at 2.76 and 5.02 TeV

Here, we perform a global Bayesian analysis of a modern event-by-event heavy-ion collision model and LHC data at √s = 2.76 and 5.02 TeV. After calibration, the model simultaneously describes multiplicity, transverse momentum, and flow data at both beam energies. We report new constraints on the scaling of initial-state entropy deposition and QGP transport coefficients, including a quantitative estimate of the temperature-dependent shear viscosity (η/s)( T ).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

A Bayesian Analysis of SDSS J0914+0853, a Low-mass Dual AGN Candidate

We present the first results from Bayesian AnalYsis of Multiple AGN in X-rays (BAYMAX), a tool that uses a Bayesian framework to quantitatively evaluate whether a given Chandra observation is more likely a single or dual point source. Although the most robust method of determining the presence of dual active galactic nuclei (AGNs) is to use X-ray observations, only sources that are widely separated relative to the instrumentʼs point-spread function are easy to identify. It becomes increasingly difficult to distinguish dual AGNs from single AGNs when the separation is on the order of Chandraʼs angular resolution (<1″). Using likelihood models for single and dual point sources, BAYMAX quantitatively evaluates the likelihood of an AGN for a given source. Specifically, we present results from BAYMAX analyzing the lowest-mass dual AGN candidate to date, SDSS J0914+0853, where archival Chandra data shows a possible secondary AGN ∼ 0"3 from the primary. Analyzing a new 50 ks Chandra observation, results from BAYMAX shows that SDSS J0914+0853 is most likely a single AGN with a Bayes factor of 13.5 in favor of a single point source model. Further, posterior distributions from the dual point source model are consistent with emission from a single AGN. We find a very low probability of SDSS J0914+0853 being a dual AGN system with a flux ratio f>0.3 and separation r>0"3. Overall, BAYMAX will be an important tool for correctly classifying candidate dual AGNs in the literature, as well as studying the dual AGN population where past spatial resolution limits have prevented systematic analyses.

Active galaxies↗

New Bayesian analysis of hybrid EoS constraints with mass-radius data for compact stars

We suggest a new Bayesian analysis using disjunct mass and radius constraints for extracting probability measures for cold, dense nuclear matter equations of state. One of the key issues of such an analysis is the question of a deconfinement transition in compact stars and whether it proceeds as a crossover or rather as a first order transition. The latter question is relevant for the possible existence of a critical endpoint in the QCD phase diagram under scrutiny in present and upcoming heavy-ion collision experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Bayesian analysis of (3 +1)⁢D relativistic nuclear dynamics with the RHIC beam energy scan data

This work presents a Bayesian inference study for relativistic heavy-ion collisions in the beam energy scan program at the BNL Relativistic Heavy-Ion Collider. The theoretical model simulates event-by-event (3+1)-dimensional [(3+1)⁢D] collision dynamics using hydrodynamics and hadronic transport theory. We analyze the model's 20-dimensional posterior distributions obtained using three model emulators with different accuracy and demonstrate the essential role of training an accurate model emulator in the Bayesian analysis. Our analysis provides robust constraints on the quark-gluon plasma's transport properties and various aspects of (3+1)⁢D relativistic nuclear dynamics. By running full model simulations with 100 parameter sets sampled from the posterior distribution, we make predictions for p T -differential observables and estimate their systematic theory uncertainty. Here, a sensitivity analysis is performed to elucidate how individual experimental observables respond to different model parameters, providing useful physics insights into the phenomenological model for heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗