ترغب بنشر مسار تعليمي؟ اضغط هنا

${tt bajes}$: Bayesian inference of multimessenger astrophysical data, methods and application to gravitational-waves

133   0   0.0 ( 0 )
 نشر من قبل Matteo Breschi
 تاريخ النشر 2021
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

We present ${tt bajes}$, a parallel and lightweight framework for Bayesian inference of multimessenger transients. ${tt bajes}$ is a Python modular package with minimal dependencies on external libraries adaptable to the majority of the Bayesian models and to various sampling methods. We describe the general workflow and the parameter estimation pipeline for compact-binary-coalescence gravitational-wave transients. The latter is validated against injections of binary black hole and binary neutron star waveforms, including confidence interval tests that demonstrates the inference is well-calibrated. Binary neutron star postmerger injections are also studied using a network of five detectors made of LIGO, Virgo, KAGRA and Einstein Telescope. Postmerger signals will be detectable for sources at ${lesssim}80,$Mpc, with Einstein Telescope contributing over 90% of the total signal-to-noise ratio. As a full scale application, we re-analyze the GWTC-1 black hole transients using the effective-one-body ${tt TEOBResumS}$ approximant, and reproduce selected results with other approximants. ${tt bajes}$ inferences are consistent with previous results; the direct comparison of ${tt bajes}$ and ${tt bilby}$ analyses of GW150914 shows a maximum Jensen-Shannon divergence of $5.2{times}10^{-4}$. GW170817 is re-analyzed using ${tt TaylorF2}$ with 5.5PN point-mass and 7.5PN tides, ${tt TEOBResumSPA}$, and ${tt IMRPhenomPv2_NRTidal}$ with different cutoff-frequencies of $1024,$Hz and $2048,$Hz. We find that the former choice minimizes systematics on the reduced tidal parameter, while a larger amount of tidal information is gained with the latter choice. ${tt bajes}$ can perform these analyses in about 1~day using 128 CPUs.

قيم البحث

اقرأ أيضاً

Third-generation (3G) gravitational-wave detectors will observe thousands of coalescing neutron star binaries with unprecedented fidelity. Extracting the highest precision science from these signals is expected to be challenging owing to both high si gnal-to-noise ratios and long-duration signals. We demonstrate that current Bayesian inference paradigms can be extended to the analysis of binary neutron star signals without breaking the computational bank. We construct reduced order models for $sim 90,mathrm{minute}$ long gravitational-wave signals, covering the observing band ($5-2048,mathrm{Hz}$), speeding up inference by a factor of $sim 1.3times 10^4$ compared to the calculation times without reduced order models. The reduced order models incorporate key physics including the effects of tidal deformability, amplitude modulation due to the Earths rotation, and spin-induced orbital precession. We show how reduced order modeling can accelerate inference on data containing multiple, overlapping gravitational-wave signals, and determine the speedup as a function of the number of overlapping signals. Thus, we conclude that Bayesian inference is computationally tractable for the long-lived, overlapping, high signal-to-noise-ratio events present in 3G observatories.
A central challenge in Gravitational Wave Astronomy is identifying weak signals in the presence of non-stationary and non-Gaussian noise. The separation of gravitational wave signals from noise requires good models for both. When accurate signal mode ls are available, such as for binary Neutron star systems, it is possible to make robust detection statements even when the noise is poorly understood. In contrast, searches for un-modeled transient signals are strongly impacted by the methods used to characterize the noise. Here we take a Bayesian approach and introduce a multi-component, variable dimension, parameterized noise model that explicitly accounts for non-stationarity and non-Gaussianity in data from interferometric gravitational wave detectors. Instrumental transients (glitches) and burst sources of gravitational waves are modeled using a Morlet-Gabor continuous wavelet frame. The number and placement of the wavelets is determined by a trans-dimensional Reversible Jump Markov Chain Monte Carlo algorithm. The Gaussian component of the noise and sharp line features in the noise spectrum are modeled using the BayesLine algorithm, which operates in concert with the wavelet model.
The Neutron Star Interior Composition Explorer (NICER) recently measured the mass and equatorial radius of the isolated neutron star PSR J0030+0451. We use these measurements to infer the moment of inertia, the quadrupole moment, and the surface ecce ntricity of an isolated neutron star for the first time, using relations between these quantities that are insensitive to the unknown equation of state of supranuclear matter. We also use these results to forecast the moment of inertia of neutron star A in the double pulsar binary J0737-3039, a quantity anticipated to be directly measured in the coming decade with radio observations. Combining this information with the measurement of the tidal Love number with LIGO/Virgo observations, we propose and implement the first theory-agnostic and equation-of-state-insensitive test of general relativity. Specializing these constraints to a particular modified theory, we find that consistency with general relativity places the most stringent constraint on gravitational parity violation to date, surpassing all other previously reported bounds by seven orders of magnitude and opens the path for future a test of general relativity with multimessenger neutron star observations.
Gravitational wave data from ground-based detectors is dominated by instrument noise. Signals will be comparatively weak, and our understanding of the noise will influence detection confidence and signal characterization. Mis-modeled noise can produc e large systematic biases in both model selection and parameter estimation. Here we introduce a multi-component, variable dimension, parameterized model to describe the Gaussian-noise power spectrum for data from ground-based gravitational wave interferometers. Called BayesLine, the algorithm models the noise power spectral density using cubic splines for smoothly varying broad-band noise and Lorentzians for narrow-band line features in the spectrum. We describe the algorithm and demonstrate its performance on data from the fifth and sixth LIGO science runs. Once fully integrated into LIGO/Virgo data analysis software, BayesLine will produce accurate spectral estimation and provide a means for marginalizing inferences drawn from the data over all plausible noise spectra.
The RIT numerical relativity group is releasing the second public catalog of black-hole-binary waveforms url{http://ccrg.rit.edu/~RITCatalog}. This release consists of 320 accurate simulations that include 46 precessing and 274 nonprecessing binary s ystems with mass ratios $q=m_1/m_2$ in the range $1/6leq qleq1$ and individual spins up to $s/m^2=0.95$. The new catalog contains search and ordering tools for the waveforms based on initial parameters of the binary, trajectory information, peak radiation, and final remnant black hole properties. The final black hole remnant properties provided here can be used to model the merger of black-hole binaries from its initial configurations. The waveforms are extrapolated to infinite observer location and can be used to independently interpret gravitational wave signals from laser interferometric detectors. As an application of this waveform catalog we reanalyze the signal of GW150914 implementing parameter estimation techniques that make use of only numerical waveforms without any reference to information from phenomenological models.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا