Search NASA⌕ Search

SEARCH · Search NASA

Results for “Credible interval”

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.

17 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Relativistic Shapiro Delay Measurements of an Extremely Massive Millisecond Pulsar

Despite its importance to our understanding of physics at supranuclear densities, the equation of state (EoS) of matter deep within neutron stars remains poorly understood. Millisecond pulsars (MSPs) are among the most useful astrophysical objects in the Universe for testing fundamental physics, and place some of the most stringent constraints on this high-density EoS. Pulsar timing—the process of accounting for every rotation of a pulsar over long time periods—can precisely measure a wide variety of physical phenomena, including those that allow the measurement of the masses of the components of a pulsar binary system1. One of these, called relativistic Shapiro delay2, can yield precise masses for both an MSP and its companion; however, it is only easily observed in a small subset of high-precision, highly inclined (nearly edge-on) binary pulsar systems. By combining data from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 12.5-yr data set with recent orbitalphase- specific observations using the Green Bank Telescope, we have measured the mass of the MSP J0740+6620 to be 2.14+0:10-0:09 M⊙ (68.3% credibility interval; the 95.4% credibility interval is 2.14+0:20-0:18 M⊙). It is highly likely to be the most massive neutron star yet observed, and serves as a strong constraint on the neutron star interior EoS.

H T Cromartie↗
NASA NTRS2016

A Comparison of Metamodeling Techniques via Numerical Experiments

This paper presents a comparative analysis of a few metamodeling techniques using numerical experiments for the single input-single output case. These experiments enable comparing the models' predictions with the phenomenon they are aiming to describe as more data is made available. These techniques include (i) prediction intervals associated with a least squares parameter estimate, (ii) Bayesian credible intervals, (iii) Gaussian process models, and (iv) interval predictor models. Aspects being compared are computational complexity, accuracy (i.e., the degree to which the resulting prediction conforms to the actual Data Generating Mechanism), reliability (i.e., the probability that new observations will fall inside the predicted interval), sensitivity to outliers, extrapolation properties, ease of use, and asymptotic behavior. The numerical experiments describe typical application scenarios that challenge the underlying assumptions supporting most metamodeling techniques.

Crespo, Luis G.↗
NASA NTRS2018

Variability Timescale and Spectral Index of Sgr A* in the Near Infrared: Approximate Bayesian Computation Analysis of the Variability of the Closest Supermassive Black Hole

Sagittarius A* (Sgr A* ) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability data set of Sgr A* to date: eight 24 hr epochs of continuous monitoring of Sgr A* at 4.5 μm with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 μm data from Naos Conica at the Very Large Telescope, and 30 epochs of 2.12 μm data from the NIRC2 camera at the Keck Observatory, in total 94,929 measurements. A new approximate Bayesian computation method for fitting the first-order structure function extracts information beyond current fast Fourier transformation (FFT) methods of power spectral density (PSD) estimation. With a combined fit of the data of all three observatories, the characteristic coherence timescale of Sgr A* is b 243 57 82 t = - + minutes (90% credible interval). The PSD has no detectable features on timescales down to 8.5 minutes (95% credible level), which is the ISCO orbital frequency for a dimensionless spin parameter a = 0.92. One light curve measured simultaneously at 2.12 and 4.5 μm during a low flux-density phase gave a spectral index αs = 1.6 ± 0.1 (F µ n s n -a ). This value implies that the Sgr A* NIR color becomes bluer during higher flux-density phases. The probability densities of flux densities of the combined data sets are best fit by log-normal distributions. Based on these distributions, the Sgr A* spectral energy distribution is consistent with synchrotron radiation from a nonthermal electron population from below 20 GHz through the NIR.

accretion↗
NASA NTRS2021

Population Properties of Compact Objects from the Second LIGO-Virgo Gravitational-Wave Transient Catalog

We report on the population of 47 compact binary mergers detected with a false-alarm rate of <1yr^(−1) in the second LIGO–Virgo Gravitational-Wave Transient Catalog. We observe several characteristics of the merging binary black hole (BBH) population not discernible until now. First, the primary mass spectrum contains structure beyond a power law with a sharp high-mass cutoff; it is more consistent with a broken power law with a break at 39.7^(+20.3)(-9.1)M or a power law with a Gaussian feature peaking at 33.1^(+4.0)(-5.6)M (90% credible interval). While the primary mass distribution must extend to ~65M or beyond, only 2.9^(+3.5)(-1.7)M% of systems have primary masses greater than 45M. Second, we find that a fraction of BBH systems have component spins misaligned with the orbital angular momentum, giving rise to precession of the orbital plane. Moreover, 12% to 44% of BBH systems have spins tilted by more than 90°, giving rise to a negative effective inspiral spin parameter, χeff. Under the assumption that such systems can only be formed by dynamical interactions, we infer that between 25% and 93% of BBHs with nonvanishing |χeff| > 0.01 are dynamically assembled. Third, we estimate merger rates, finding RBBH = 23.9^(+14.3)(-8.6) Gpc^(-3) yr^(-1) for BBHs and RBNS = 320^(+490)(-240) Gpc^(-3) yr^(-1) for binary neutron stars. We find that the BBH rate likely increases with redshift (85% credibility) but not faster than the star formation rate (86% credibility). Additionally, we examine recent exceptional events in the context of our population models, finding that the asymmetric masses of GW190412 and the high component masses of GW190521 are consistent with our models, but the low secondary mass of GW190814 makes it an outlier.

R Abbott↗
NASA NTRS2015

Parameter Estimation for Compact Binaries with Ground-Based Gravitational-Wave Observations Using the LALInference

The Advanced LIGO and Advanced Virgo gravitational wave (GW) detectors will begin operation in the coming years, with compact binary coalescence events a likely source for the first detections. The gravitational waveforms emitted directly encode information about the sources, including the masses and spins of the compact objects. Recovering the physical parameters of the sources from the GW observations is a key analysis task. This work describes the LALInference software library for Bayesian parameter estimation of compact binary signals, which builds on several previous methods to provide a well-tested toolkit which has already been used for several studies. We show that our implementation is able to correctly recover the parameters of compact binary signals from simulated data from the advanced GW detectors. We demonstrate this with a detailed comparison on three compact binary systems: a binary neutron star (BNS), a neutron star - black hole binary (NSBH) and a binary black hole (BBH), where we show a cross-comparison of results obtained using three independent sampling algorithms. These systems were analysed with non-spinning, aligned spin and generic spin configurations respectively, showing that consistent results can be obtained even with the full 15-dimensional parameter space of the generic spin configurations. We also demonstrate statistically that the Bayesian credible intervals we recover correspond to frequentist confidence intervals under correct prior assumptions by analysing a set of 100 signals drawn from the prior. We discuss the computational cost of these algorithms, and describe the general and problem-specific sampling techniques we have used to improve the efficiency of sampling the compact binary coalescence (CBC) parameter space.

gravitational wave↗
NASA NTRS2016

Observation of Gravitational Waves from a Binary Black Hole Merger

On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-Wave Observatory simultaneously observed a transient gravitational-wave signal. The signal sweeps upwards in frequency from 35 to 250 Hz with a peak gravitational-wave strain of 1.0 x 10(exp -21). It matches the waveform predicted by general relativity for the inspiral and merger of a pair of black holes and the ring down of the resulting single black hole. The signal was observed with a matched-filter signal-to-noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203 000 years, equivalent to a significance greater than 5.1 Sigma. The source lies at a luminosity distance of 410(+160/-180) Mpc corresponding to a redshift z = 0.09(+0.03/-0.04). In the source frame, the initial black hole masses are 36(+5/-4) Mass compared to the sun, and 29(+4/-4) Mass compared to the sun, and the final black hole mass is 62(+4/-4) Mass compared to the sun, with 3.0(+0.5/-0.5)sq c radiated in gravitational waves. All uncertainties define 90% credible intervals. These observations demonstrate the existence of binary stellar-mass black hole systems. This is the first direct detection of gravitational waves and the first observation of a binary black hole merger.

gravitational waves↗
NASA NTRS2016

Properties of the Binary Black Hole Merger GW150914

On September 14, 2015, the Laser Interferometer Gravitational-Wave Observatory (LIGO) detected a gravitational-wave transient (GW150914); we characterize the properties of the source and its parameters. The data around the time of the event were analyzed coherently across the LIGO network using a suite of accurate waveform models that describe gravitational waves from a compact binary system in general relativity. GW150914 was produced by a nearly equal mass binary black hole of masses 36(+5/-4) solar mass and 29(+4/-4) solar mass; for each parameter we report the median value and the range of the 90% credible interval. The dimensionless spin magnitude of the more massive black hole is bound to be less than 0.7 (at 90% probability). The luminosity distance to the source is 410(+160/-180) Mpc, corresponding to a redshift 0.09(+0.03/−0.04) assuming standard cosmology. The source location is constrained to an annulus section of 610 sq deg, primarily in the southern hemisphere. The binary merges into a black hole of mass 62(+4/−4) solar mass and spin 0.67(+0.05/−0.07). This black hole is significantly more massive than any other inferred from electromagnetic observations in the stellar-mass regime.

Abbott, B. P.↗
NASA NTRS2016

GW151226: Observation of Gravitational Waves from a 22-Solar-Mass Binary Black Hole Coalescence

We report the observation of a gravitational-wave signal produced by the coalescence of two stellar-mass black holes. The signal, GW151226, was observed by the twin detectors of the Laser Interferometer Gravitational-Wave Observatory (LIGO) on December 26, 2015 at 03:38:53 UTC. The signal was initially identified within 70 s by an online matched-filter search targeting binary coalescences. Subsequent off-line analyses recovered GW151226 with a network signal-to-noise ratio of 13 and a significance greater than 5(sigma). The signal persisted in the LIGO frequency band for approximately 1 s, increasing in frequency and amplitude over about 55 cycles from 35 to 450 Hz, and reached a peak gravitational strain of 3.4(+0.7/-0.9) x 10(exp -22). The inferred source-frame initial black hole masses are 14.2(+8.3/-3.7 Stellar Mass and 7.5(+2.3/-2.3) Stellar Mass, and the final black hole mass is 20.8(+6.1/-1.7) Stellar Mass. We find that at least one of the component black holes has spin greater than 0.2. This source is located at a luminosity distance of 440(+180/-190) Mpc corresponding to a redshift of 0.090(+.030/-0.04). All uncertainties define a 90% credible interval. This second gravitational-wave observation provides improved constraints on stellar populations and on deviations from general relativity.

Abbott, B. P.↗
NASA NTRS2016

Improved Analysis of GW150914 Using a Fully Spin-Precessing Waveform Model

This paper presents updated estimates of source parameters for GW150914, a binary black-hole coalescence event detected by the Laser Interferometer Gravitational-wave Observatory (LIGO) in 2015 [Abbott et al. Phys. Rev. Lett. 116, 061102 (2016).]. Abbott et al. [Phys. Rev. Lett. 116, 241102 (2016).] presented parameter estimation of the source using a 13-dimensional, phenomenological precessing-spin model (precessing IMRPhenom) and an 11-dimensional nonprecessing effective-one-body (EOB) model calibrated to numerical-relativity simulations, which forces spin alignment (nonprecessing EOBNR). Here, we present new results that include a 15-dimensional precessing-spin waveform model (precessing EOBNR) developed within the EOB formalism. We find good agreement with the parameters estimated previously [Abbott et al. Phys. Rev. Lett. 116, 241102 (2016).], and we quote updated component masses of 35(+5)(-3) solar M; and 30(+3)(-4) solar M; (where errors correspond to 90 symmetric credible intervals). We also present slightly tighter constraints on the dimensionless spin magnitudes of the two black holes, with a primary spin estimate is less than 0.65 and a secondary spin estimate is less than 0.75 at 90% probability. Abbott et al. [Phys. Rev. Lett. 116, 241102 (2016).] estimated the systematic parameter-extraction errors due to waveform-model uncertainty by combining the posterior probability densities of precessing IMRPhenom and nonprecessing EOBNR. Here, we find that the two precessing-spin models are in closer agreement, suggesting that these systematic errors are smaller than previously quoted.

Abbott, B. P.↗
NASA NTRS2019

A NICER view of PSR J0030+0541: Millisecond Pulsar Parameter Estimation

We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSRJ0030+0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray spectral timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar’s surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass M and equatorial radius Req are, respectively, -1.34+ M 0.16 0.15 and-12.71+ km 1.19 1.14 , while the compactness = -GM R c 0.156+ eq 2 0.010 0.008 is more tightly constrained; the credible interval bounds reported here are approximately the 16% and 84% quantiles in marginal posterior mass.

T E Riley↗
NASA NTRS2020

GW190521: A Binary Black Hole Merger with a Total Mass of 150 Mꙩ

On May 21, 2019 at 03:02:29 UTC Advanced LIGO and Advanced Virgo observed a short duration gravitational-wave signal, GW190521, with a three-detector network signal-to-noise ratio of 14.7, and an estimated false-alarm rate of 1 in 4900 yr using a search sensitive to generic transients. If GW190521 is from a quasicircular binary inspiral, then the detected signal is consistent with the merger of two black holes with masses of 85 (+21,−14) Mꙩ and 66 (+17,−18) Mꙩ (90% credible intervals). We infer that the primary black hole mass lies within the gap produced by (pulsational) pair-instability supernova processes, with only a 0.32% probability of being below 65 M⊙. We calculate the mass of the remnant to be 142 (+28,−16) Mꙩ, which can be considered an intermediate mass black hole (IMBH). The luminosity distance of the source is 5.3 (+2.4,−2.6) Gpc, corresponding to a redshift of 0.82 (+0.28,−0.34) . The inferred rate of mergers similar to GW190521 is 0.13 (+0.30,−0.11) Gpc^(−3) per yr.

R. Abbott↗
NASA NTRS2021

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be-+12.390.981.30km and-+2.0720.0660.067Me respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is=-+TlogK5.99100.060.05([])for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters(within the prior support)and can thus be computed from the posterior samples.

Millisecond pulsars↗
NASA NTRS2011

The Extravehicular Suit Impact Load Attenuation Study for Use in Astronaut Bone Fracture Prediction

The NASA Integrated Medical Model (IMM) assesses the risk, including likelihood and impact of occurrence, of all credible in-flight medical conditions. Fracture of the proximal femur is a traumatic injury that would likely result in loss of mission if it were to happen during spaceflight. The low gravity exposure causes decreases in bone mineral density which heightens the concern. Researchers at the NASA Glenn Research Center have quantified bone fracture probability during spaceflight with a probabilistic model. It was assumed that a pressurized extravehicular activity (EVA) suit would attenuate load during a fall, but no supporting data was available. The suit impact load attenuation study was performed to collect analogous data. METHODS: A pressurized EVA suit analog test bed was used to study how the offset, defined as the gap between the suit and the astronaut s body, impact load magnitude and suit operating pressure affects the attenuation of impact load. The attenuation data was incorporated into the probabilistic model of bone fracture as a function of these factors, replacing a load attenuation value based on commercial hip protectors. RESULTS: Load attenuation was more dependent on offset than on pressurization or load magnitude, especially at small offsets. Load attenuation factors for offsets between 0.1 - 1.5 cm were 0.69 +/- 0.15, 0.49 +/- 0.22 and 0.35 +/- 0.18 for mean impact forces of 4827, 6400 and 8467 N, respectively. Load attenuation factors for offsets of 2.8 - 5.3 cm were 0.93 +/- 0.2, 0.94 +/- 0.1 and 0.84 +/- 0.5, for the same mean impact forces. Reductions were observed in the 95th percentile confidence interval of the bone fracture probability predictions. CONCLUSIONS: The reduction in uncertainty and improved confidence in bone fracture predictions increased the fidelity and credibility of the fracture risk model and its benefit to mission design and operational decisions.

Lewandowski, Beth E.↗
NASA NTRS2021

The Radius of PSR J0740+6620 from NICER and XMM-Newton Data

PSR J0740+6620 has a gravitational mass of 2.08±0.07Me, which is the highest reliably determined mass of any neutron star. As a result, a measurement of its radius will provide unique insight into the properties of neutron star core matter at high densities. Here we report a radius measurement based on fits of rotating hot spot patterns to Neutron Star Interior Composition Explorer(NICER)and X-ray Multi-Mirror(XMM-Newton)X-ray observations. We find that the equatorial circumferential radius of PSR J0740+6620 is-+13.71.52.6km(68%).We apply our measurement, combined with the previous NICER mass and radius measurement of PSR J0030+0451,the masses of two other∼2Mepulsars, and the tidal deformability constraints from two gravitational wave events, to three different frameworks for equation-of-state modeling, and find consistent results at∼1.5–5 times nuclear saturation density. For a given framework, when all measurements are included, the radius of a 1.4Meneutron star is known to±4%(68% credibility)and the radius of a 2.08Meneutron star is known to±5%. The full radius range that spans the±1σcredible intervals of all the radius estimates in the three frameworks is 12.45±0.65 km for a 1.4Meneutron star and 12.35±0.75 km for a 2.08Me neutron star.

M C Miller↗
NASA NTRS1972

Supplementary notes on sea surface temperature anomalies and model-generated meteorological histories

In seasonal computations, the Mintz-Arakawa two-level model is found to be sensitive to a minor alteration in the computational program. Effects of the program change on monthly mean sea level pressure fields are small in the first month, but large in the second and third months, although the meteorological histories generated by both the original and modified programs are equally credible. The inherited effects of a transient (one month) sea surface temperature (SST) anomaly on the computed monthly mean sea level pressure fields over a period of a season are about as large in absolute magnitude as those generated in the model by a persistent (seasonal) SST anomaly. The effects of a transient SST anomaly in the North Pacific Ocean on monthly and seasonal temperature and precipitation in the eastern United States may be large enough to produce a change of one or two class intervals in these predicted weather elements. The model-generated precipitation in the equatorial region is also found to be sensitive to the sea surface temperature field in the North Pacific.

Spar, J.↗
NASA NTRS2022

Unique Observational Constraints on the Seasonal and Longitudinal Variability of the Earth’s Planetary Albedo and Cloud Distribution Inferred From EPIC Measurements

Thorough comparison to observations is key to developing a credible climate model forecasting capability. Deep Space Climate Observatory (DSCOVR) measurements of Earth’s reflected solar and emitted thermal radiation provide a unique observational perspective that permits a more reliable model/data comparison than is possible with the otherwise available satellite data. The uniqueness is in the DSCOVR satellite’s viewing geometry, which enables continuous viewing of the Earth’s sunlit hemisphere from its Lissajous orbit around the Lagrangian L1 point. The key instrument is the Earth Polychromatic Imaging Camera (EPIC), which views the Earth’s sunlit hemisphere with 1024-by-1024-pixel imagery in 10 narrow spectral bands from 317 to 780 nm, acquiring up to 22 high spatial resolution images per day. The additional feature is that the frequency of EPIC image acquisition is nearly identical to that of the climate GCM data generation scheme where climate data for the entire globe are ‘instantaneously’ calculated at 1-h radiation time-step intervals. Implementation of the SHS (Sunlit Hemisphere Sampling) EPIC-view geometry for the in-line GCM output data sampling establishes a precise self-consistency in the space-time data sampling between EPIC observational and GCM output data generation and sampling. The remaining problem is that the GCM generated data are radiative fluxes, while the EPIC measurements are backscatter-dependent radiances. Radiance to flux conversion is a complex problem with no simple way to convert GCM radiative fluxes into spectral radiances. The more expedient approach is to convert the EPIC spectral radiances into broadband radiances by MODIS/CERES-based regression relationships and then into solar radiative fluxes using the CERES angular distribution models. Averaging over the sunlit hemisphere suppresses the meteorological weather noise, but preserves the intra-seasonal larger scale variability. Longitudinal slicing by the Earth’s rotation permits a self-consistent model/data comparison of the longitudinal model/data differences in the variability of the reflected solar radiation. Ancillary EPIC Composite data provide additional cloud property information for climate model diagnostics. Comparison of EPIC-derived seasonal and longitudinal variability of the Earth’s planetary albedo with the GISS ModelE2 results shows systematic overestimate of cloud reflectivity over the Pacific Ocean with corresponding underestimates over continental land areas.

DSCOVR↗
NASA NTRS2011

Sleep Disruption Medical Intervention Forecasting (SDMIF) Module for the Integrated Medical Model

The NASA Integrated Medical Model (IMM) assesses the risk, including likelihood and impact of occurrence, of all credible in-flight medical conditions. Fatigue due to sleep disruption is a condition that could lead to operational errors, potentially resulting in loss of mission or crew. Pharmacological consumables are mitigation strategies used to manage the risks associated with sleep deficits. The likelihood of medical intervention due to sleep disruption was estimated with a well validated sleep model and a Monte Carlo computer simulation in an effort to optimize the quantity of consumables. METHODS: The key components of the model are the mission parameter program, the calculation of sleep intensity and the diagnosis and decision module. The mission parameter program was used to create simulated daily sleep/wake schedules for an ISS increment. The hypothetical schedules included critical events such as dockings and extravehicular activities and included actual sleep time and sleep quality. The schedules were used as inputs to the Sleep, Activity, Fatigue and Task Effectiveness (SAFTE) Model (IBR Inc., Baltimore MD), which calculated sleep intensity. Sleep data from an ISS study was used to relate calculated sleep intensity to the probability of sleep medication use, using a generalized linear model for binomial regression. A human yes/no decision process using a binomial random number was also factored into sleep medication use probability. RESULTS: These probability calculations were repeated 5000 times resulting in an estimate of the most likely amount of sleep aids used during an ISS mission and a 95% confidence interval. CONCLUSIONS: These results were transferred to the parent IMM for further weighting and integration with other medical conditions, to help inform operational decisions. This model is a potential planning tool for ensuring adequate sleep during sleep disrupted periods of a mission.

Lewandowski, Beth↗