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

Parameter estimation for gravitational-wave bursts with the BayesWave pipeline

122   0   0.0 ( 0 )
 نشر من قبل Bence B\\'ecsy
 تاريخ النشر 2016
  مجال البحث فيزياء
والبحث باللغة English




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

We provide a comprehensive multi-aspect study on the performance of a pipeline used by the LIGO-Virgo Collaboration for estimating parameters of gravitational-wave bursts. We add simulated signals with four different morphologies (sine-Gaussians, Gaussians, white-noise bursts, and binary black hole signals) to simulated noise samples representing noise of the two Advanced LIGO detectors during their first observing run. We recover them with the BayesWave (BW) pipeline to study its accuracy in sky localization, waveform reconstruction, and estimation of model-independent waveform parameters. BW localizes sources with a level of accuracy comparable for all four morphologies, with the median separation of actual and estimated sky locations ranging from 25.1$^{circ}$ to 30.3$^{circ}$. This is a reasonable accuracy in the two-detector case, and is comparable to accuracies of other localization methods studied previously. As BW reconstructs generic transient signals with sine-Gaussian wavelets, it is unsurprising that BW performs the best in reconstructing sine-Gaussian and Gaussian waveforms. BWs accuracy in waveform reconstruction increases steeply with network signal-to-noise ratio (SNR$_{rm net}$), reaching a $85%$ and $95%$ match between the reconstructed and actual waveform below SNR$_{rm net} approx 20$ and SNR$_{rm net} approx 50$, respectively, for all morphologies. BWs accuracy in estimating central moments of waveforms is only limited by statistical errors in the frequency domain, and is affected by systematic errors too in the time domain as BW cannot reconstruct low-amplitude parts of signals overwhelmed by noise. The figures of merit we introduce can be used in future characterizations of parameter estimation pipelines.

قيم البحث

اقرأ أيضاً

We describe updates and improvements to the BayesWave gravitational wave transient analysis pipeline, and provide examples of how the algorithm is used to analyze data from ground-based gravitational wave detectors. BayesWave models gravitational wav e signals in a morphology-independent manner through a sum of frame functions, such as Morlet-Gabor wavelets or chirplets. BayesWave models the instrument noise using a combination of a parametrized Gaussian noise component and non-stationary and non-Gaussian noise transients. Both the signal model and noise model employ trans-dimensional sampling, with the complexity of the model adapting to the requirements of the data. The flexibility of the algorithm makes it suitable for a variety of analyses, including reconstructing generic unmodeled signals; cross checks against modeled analyses for compact binaries; as well as separating coherent signals from incoherent instrumental noise transients (glitches). The BayesWave model has been extended to account for gravitational wave signals with generic polarization content and the simultaneous presence of signals and glitches in the data. We describe updates in the BayesWave prior distributions, sampling proposals, and burn-in stage that provide significantly improved sampling efficiency. We present standard review checks indicating the robustness and convergence of the BayesWave trans-dimensional sampler.
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.
One of the main bottlenecks in gravitational wave (GW) astronomy is the high cost of performing parameter estimation and GW searches on the fly. We propose a novel technique based on Reduced Order Quadratures (ROQs), an application and data-specific quadrature rule, to perform fast and accurate likelihood evaluations. These are the dominant cost in Markov chain Monte Carlo (MCMC) algorithms, which are widely employed in parameter estimation studies, and so ROQs offer a new way to accelerate GW parameter estimation. We illustrate our approach using a four dimensional GW burst model embedded in noise. We build an ROQ for this model, and perform four dimensional MCMC searches with both the standard and ROQs quadrature rules, showing that, for this model, the ROQ approach is around 25 times faster than the standard approach with essentially no loss of accuracy. The speed-up from using ROQs is expected to increase for more complex GW signal models and therefore has significant potential to accelerate parameter estimation of GW sources such as compact binary coalescences.
The detection of electromagnetic counterparts to gravitational waves has great promise for the investigation of many scientific questions. It has long been hoped that in addition to providing extra, non-gravitational information about the sources of these signals, the detection of an electromagnetic signal in conjunction with a gravitational wave could aid in the analysis of the gravitational signal itself. That is, knowledge of the sky location, inclination, and redshift of a binary could break degeneracies between these extrinsic, coordinate-dependent parameters and the physical parameters, such as mass and spin, that are intrinsic to the binary. In this paper, we investigate this issue by assuming a perfect knowledge of extrinsic parameters, and assessing the maximal impact of this knowledge on our ability to extract intrinsic parameters. However, we find only modest improvements in a few parameters --- namely the primary components spin --- and conclude that, even in the best case, the use of additional information from electromagnetic observations does not improve the measurement of the intrinsic parameters significantly.
We performed a search for event bursts in the XMASS-I detector associated with 11 gravitational-wave events detected during LIGO/Virgos O1 and O2 periods. Simple and loose cuts were applied to the data collected in the full 832 kg xenon volume around the detection time of each gravitational-wave event. The data were divided into four energy regions ranging from keV to MeV. Without assuming any particular burst models, we looked for event bursts in sliding windows with various time width from 0.02 to 10 s. The search was conducted in a time window between $-$400 and $+$10,000 s from each gravitational-wave event. For the binary neutron star merger GW170817, no significant event burst was observed in the XMASS-I detector and we set 90% confidence level upper limits on neutrino fluence for the sum of all the neutrino flavors via coherent elastic neutrino-nucleus scattering. The obtained upper limit was (1.3-2.1)$times 10^{11}$ cm$^{-2}$ under the assumption of a Fermi-Dirac spectrum with average neutrino energy of 20 MeV. The neutrino fluence limits for mono-energetic neutrinos in the energy range between 14 and 100 MeV were also calculated. Among the other 10 gravitational wave events detected as the binary black hole mergers, a burst candidate with a 3.0$sigma$ significance was found at 1801.95-1803.95 s in the analysis for GW151012. However, no significant deviation from the background in the reconstructed energy and position distributions was found. Considering the additional look-elsewhere effect of analyzing the 11 GW events, the significance of finding such a burst candidate associated with any of them is 2.1$sigma$.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

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