BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

John Veitch

John Veitch appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-ph.IMastro-ph.COphysics.data-anastro-ph.SRComputationcond-mat.mtrl-scihep-ph

Catalog footprint

What is connected

31works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2023arXiv

Rapid pre-merger localization of binary neutron stars in third generation gravitational wave detectors

Pre-merger localization of binary neutron stars (BNSs) is one of the most important scientific goals for the third generation (3G) gravitational wave (GW) detectors. It will enable the electromagnetic observation of the whole process of BNS coalescence, especially for the pre-merger and merger phases which have not been observed yet, opening a window for deeper understandings of compact objects. To reach this goal, we describe a novel combination of multi-band matched filtering and semi-analytical localization algorithms to achieve early-warning localization of long BNS signals in 3G detectors. Using our method we are able to efficiently simulate one month of observations with a three-detector 3G network, and show that it is possible to provide accurate sky localizations more than 30 minutes before the merger. Our simulation shows that there could be ~ 10 (~ 100) BNS events localized within 100 deg2, 20 (6) minutes before merger, per month of observation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A Pixelated Approach to Galaxy Catalogue Incompleteness: Improving the Dark Siren Measurement of the Hubble Constant

The use of gravitational wave standard sirens for cosmological analyses is becoming well known, with particular interest in measuring the Hubble constant, $H_0$, and in shedding light on the current tension between early- and late-time measurements. The current tension is over $4σ$ and standard sirens will be able to provide a completely independent measurement. Dark sirens (binary black hole or neutron star mergers with no electromagnetic counterparts) can be informative if the missing redshift information is provided through the use of galaxy catalogues to identify potential host galaxies of the merger. However, galaxy catalogue incompleteness affects this analysis, and accurate modelling of it is essential for obtaining an unbiased measurement of $H_0$. Previously this has been done by assuming uniform completeness within the sky area of a gravitational wave event. This paper presents an updated methodology in which the completeness of the galaxy catalogue is estimated in a directionally-dependent matter, by pixelating the sky and computing the completeness of the galaxy catalogue along each line of sight. The $H_0$ inference for a single event is carried out on a pixel-by-pixel basis, and the pixels are combined for the final result. A reanalysis of the events in the first gravitational wave transient catalogue (GWTC-1) leads to an improvement on the measured value of $H_0$ of approximately 5% compared to the 68.3% highest-density interval of the equivalent LIGO and Virgo result, with $H_0= 68.8^{+15.9}_{-7.8}$ km s$^{-1}$ Mpc$^{-1}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Asimov: A framework for coordinating parameter estimation workflows

Since the first detection in 2015 of gravitational waves from compact binary coalescence, improvements to the Advanced LIGO and Advanced Virgo detectors have expanded our view into the universe for these signals. Searches of the of the latest observing run (O3) have increased the number of detected signals to 90, at a rate of approximately 1 per week. Future observing runs are expected to increase this even further. Bayesian analysis of the signals can reveal the properties of the coalescing black holes and neutron stars by comparing predicted waveforms to the observed data. The proliferating number of detected signals, the increasing number of methods that have been deployed, and the variety of waveform models create an ever-expanding number of analyses that can be considered. Asimov is a python package which is designed to simplify and standardise the process of configuring these analyses for a large number of events. It has already been used in developing analyses in three major gravitational wave catalog publications.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Assessing the model waveform accuracy of gravitational waves

With the improvement in sensitivity of gravitational wave (GW) detectors and the increasing diversity of GW sources, there is a strong need for accurate GW waveform models for data analysis. While the current model accuracy assessments require waveforms generated by numerical relativity (NR) simulations as the "true waveforms", in this paper we propose an assessment approach that does not require NR simulations, which enables us to assess model accuracy everywhere in the parameter space. By measuring the difference between two waveform models, we derive a necessary condition for a pair of waveform models to both be accurate, for a particular set of parameters. We then apply this method to the parameter estimation samples of the Gravitational-Wave Transient Catalogs GWTC-3 and GWTC-2.1, and find that the waveform accuracy for high signal-to-noise ratio events in some cases fails our assessment criterion. Based on analysis of real events' posterior samples, we discuss the correlation between our quantified accuracy assessments and systematic errors in parameter estimation. We find waveform models that perform worse in our assessment are more likely to give inconsistent estimations. We also investigate waveform accuracy in different parameter regions, and find the accuracy degrades as the spin effects go up, the mass ratio deviates from one, or the orbital plane is near-aligned to the line of sight. Furthermore, we make predictions of waveform accuracy requirements for future detectors and find the accuracy of current waveform models should be improved by at least 3 orders of magnitude, which is consistent with previous works.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nested sampling for physical scientists

We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the application of NS in three scientific fields in which the algorithm has proved to be useful: cosmology, gravitational-wave astronomy, and materials science. We close by making recommendations for best practice when using NS and by summarizing potential limitations and optimizations of NS.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Comparing Short Gamma-Ray Burst Jet Structure Models

A structured gamma-ray burst jet could explain the dimness of the prompt emission observed from GRB$\,170817$A but the exact form of this structure is still ambiguous. However, with the promise of future joint gravitational wave and gamma-ray burst observations, we shall be able to examine populations of binary neutron star mergers rather than a case-by-case basis. We present an analysis that considers gravitational wave triggered binary neutron star events both with and without short gamma-ray burst counterparts assuming that events without a counterpart were observed off-axis. This allows for Bayes factors to be calculated to compare different jet structure models. We perform model comparison between a Gaussian and power-law apparent jet structure on simulated data to demonstrate that the correct model can be distinguished with a log Bayes factor of $>5$ after less than 100 events. Constraints on the apparent structure jet model parameters are also made. After 25(100) events the angular width of the core of a power-law jet structure can be constrained within a $90\%$ credible interval of width $ \sim9.1(4.4)^{\circ} $, and the outer beaming angle to be within $\sim19.9(8.5)^{\circ}$. Similarly we show the width of a Gaussian jet structure to be constrained to $\sim2.8(1.6)^{\circ}$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cosmological Inference using Gravitational Wave Standard Sirens: A Mock Data Challenge

The observation of binary neutron star merger GW170817, along with its optical counterpart, provided the first constraint on the Hubble constant $H_0$ using gravitational wave standard sirens. When no counterpart is identified, a galaxy catalog can be used to provide the necessary redshift information. However, the true host might not be contained in a catalog which is not complete out to the limit of gravitational-wave detectability. These electromagnetic and gravitational-wave selection effects must be accounted for. We describe and implement a method to estimate $H_0$ using both the counterpart and the galaxy catalog standard siren methods. We perform a series of mock data analyses using binary neutron star mergers to confirm our ability to recover an unbiased estimate of $H_0$. Our simulations used a simplified universe with no redshift uncertainties or galaxy clustering, but with different magnitude-limited catalogs and assumed host galaxy properties, to test our treatment of both selection effects. We explore how the incompleteness of catalogs affects the final measurement of $H_0$, as well as the effect of weighting each galaxy's likelihood of being a host by its luminosity. In our most realistic simulation, where the simulated catalog is about three times denser than the density of galaxies in the local universe, we find that a 4.4\% measurement precision can be reached using galaxy catalogs with 50\% completeness and $\sim 250$ binary neutron star detections with sensitivity similar to that of Advanced LIGO's second observing run.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Parameter estimation on gravitational waves from neutron-star binaries with spinning components

Inspiraling binary neutron stars are expected to be one of the most significant sources of gravitational-wave signals for the new generation of advanced ground-based detectors. We investigate how well we could hope to measure properties of these binaries using the Advanced LIGO detectors, which began operation in September 2015. We study an astrophysically motivated population of sources (binary components with masses $1.2~\mathrm{M}_\odot$--$1.6~\mathrm{M}_\odot$ and spins of less than $0.05$) using the full LIGO analysis pipeline. While this simulated population covers the observed range of potential binary neutron-star sources, we do not exclude the possibility of sources with parameters outside these ranges; given the existing uncertainty in distributions of mass and spin, it is critical that analyses account for the full range of possible mass and spin configurations. We find that conservative prior assumptions on neutron-star mass and spin lead to average fractional uncertainties in component masses of $\sim 16\%$, with little constraint on spins (the median $90\%$ upper limit on the spin of the more massive component is $\sim 0.7$). Stronger prior constraints on neutron-star spins can further constrain mass estimates, but only marginally. However, we find that the sky position and luminosity distance for these sources are not influenced by the inclusion of spin; therefore, if LIGO detects a low-spin population of BNS sources, less computationally expensive results calculated neglecting spin will be sufficient for guiding electromagnetic follow-up.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constraining the neutron star equation of state with gravitational wave signals from coalescing binary neutron stars

Recently exploratory studies were performed on the possibility of constraining the neutron star equation of state (EOS) using signals from coalescing binary neutron stars, or neutron star-black hole systems, as they will be seen in upcoming advanced gravitational wave detectors such as Advanced LIGO and Advanced Virgo. In particular, it was estimated to what extent the combined information from multiple detections would enable one to distinguish between different equations of state through hypothesis ranking or parameter estimation. Under the assumption of zero neutron star spins both in signals and in template waveforms and considering tidal effects to 1 post-Newtonian (1PN) order, it was found that O(20) sources would suffice to distinguish between a hard, moderate, and soft equation of state. Here we revisit these results, this time including neutron star tidal effects to the highest order currently known, termination of gravitational waveforms at the contact frequency, neutron star spins, and the resulting quadrupole-monopole interaction. We also take the masses of neutron stars in simulated sources to be distributed according to a relatively strongly peaked Gaussian, as hinted at by observations, but without assuming that the data analyst will necessarily have accurate knowledge of this distribution for use as a mass prior. We find that especially the effect of the latter is dramatic, necessitating many more detections to distinguish between different EOS and causing systematic biases in parameter estimation, on top of biases due to imperfect understanding of the signal model pointed out in earlier work. This would get mitigated if reliable prior information about the mass distribution could be folded into the analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High-energy electromagnetic offline follow-up of LIGO-Virgo gravitational-wave binary coalescence candidate events

We present two different search methods for electromagnetic counterparts to gravitational-wave (GW) events from ground-based detectors using archival NASA high-energy data from the Fermi-GBM and RXTE-ASM instruments. To demonstrate the methods, we use a limited number of representative GW background noise events produced by a search for binary neutron star coalescence over the last two months of the LIGO-Virgo S6/VSR3 joint science run. Time and sky location provided by the GW data trigger a targeted search in the high-energy photon data. We use two custom pipelines: one to search for prompt gamma-ray counterparts in GBM, and the other to search for a variety of X-ray afterglow model signals in ASM. We measure the efficiency of the joint pipelines to weak gamma-ray burst counterparts, and a family of model X-ray afterglows. By requiring a detectable signal in either electromagnetic instrument coincident with a GW event, we are able to reject a large majority of GW candidates. This reduces the signal-to-noise of the loudest surviving GW background event by around 15-20%.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inference on gravitational waves from coalescences of stellar-mass compact objects and intermediate-mass black holes

Gravitational waves from coalescences of neutron stars or stellar-mass black holes into intermediate-mass black holes (IMBHs) of $\gtrsim 100$ solar masses represent one of the exciting possible sources for advanced gravitational-wave detectors. These sources can provide definitive evidence for the existence of IMBHs, probe globular-cluster dynamics, and potentially serve as tests of general relativity. We analyse the accuracy with which we can measure the masses and spins of the IMBH and its companion in intermediate-mass ratio coalescences. We find that we can identify an IMBH with a mass above $100 ~ M_\odot$ with $95\%$ confidence provided the massive body exceeds $130 ~ M_\odot$. For source masses above $\sim200 ~ M_\odot$, the best measured parameter is the frequency of the quasi-normal ringdown. Consequently, the total mass is measured better than the chirp mass for massive binaries, but the total mass is still partly degenerate with spin, which cannot be accurately measured. Low-frequency detector sensitivity is particularly important for massive sources, since sensitivity to the inspiral phase is critical for measuring the mass of the stellar-mass companion. We show that we can accurately infer source parameters for cosmologically redshifted signals by applying appropriate corrections. We investigate the impact of uncertainty in the model gravitational waveforms and conclude that our main results are likely robust to systematics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Measuring intermediate mass black hole binaries with advanced gravitational wave detectors

We perform a systematic study to explore the accuracy with which the parameters of intermediate-mass black-hole binary systems can be measured from their gravitational wave (GW) signatures using second-generation GW detectors. We make use of the most recent reduced-order models containing inspiral, merger and ringdown signals of aligned-spin effective-one-body waveforms (SEOBNR) to significantly speed up the calculations. We explore the phenomenology of the measurement accuracies for binaries with total masses between 50 and 500 $M_\odot$ and mass ratios between 0.1 and 1. We find that (i) at total masses below ~200 $M_\odot$, where the signal-to-noise-ratio is dominated by the inspiral portion of the signal, the chirp mass parameter can be accurately measured; (ii) at higher masses, the information content is dominated by the ringdown, and total mass is measured more accurately; (iii) the mass of the lower-mass companion is poorly estimated, especially at high total mass and more extreme mass ratios; (iv) spin cannot be accurately measured for our injection set with non-spinning components. Most importantly, we find that for binaries with non-spinning components at all values of the mass ratio in the considered range and at network signal-to-noise ratio of 15, analyzed with spin-aligned templates, the presence of an intermediate-mass black hole with mass >100 $M_\odot$ can be confirmed with 95% confidence in any binary that includes a component with a mass of 130 $M_\odot$ or greater.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Parameter estimation for binary neutron-star coalescences with realistic noise during the Advanced LIGO era

Advanced ground-based gravitational-wave (GW) detectors begin operation imminently. Their intended goal is not only to make the first direct detection of GWs, but also to make inferences about the source systems. Binary neutron-star mergers are among the most promising sources. We investigate the performance of the parameter-estimation \edit{(PE)} pipeline that will be used during the first observing run of the Advanced Laser Interferometer Gravitational-wave Observatory (aLIGO) in 2015: we concentrate on the ability to reconstruct the source location on the sky, but also consider the ability to measure masses and the distance. Accurate, rapid sky-localization is necessary to alert electromagnetic (EM) observatories so that they can perform follow-up searches for counterpart transient events. We consider PE accuracy in the presence of \edit{non-stationary}, non-Gaussian noise. We find that the character of the noise makes negligible difference to the PE performance \edit{at a given signal-to-noise ratio}. The source luminosity distance can only be poorly constrained, the median $90\%$ ($50\%$) credible interval scaled with respect to the true distance is $0.85$ ($0.38$). However, the chirp mass is well measured. Our chirp-mass estimates are subject to systematic error because we used gravitational-waveform templates without component spin to carry out inference on signals with moderate spins, but the total error is typically less than $10^{-3} M_\odot$. The median $90\%$ ($50\%$) credible region for sky localization is $\sim600~\mathrm{deg^{2}}$ ($\sim150~\mathrm{deg^{2}}$), with $3\%$ ($30\%$) of detected events localized within $100~\mathrm{deg^{2}}$. Early aLIGO, with only two detectors, will have a sky-localization accuracy for binary neutron stars of hundreds of square degrees; this makes EM follow-up challenging, but not impossible.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Planned search for LIGO-GBM coincidence in the first advanced LIGO data run

In the fall of 2015 the first scientific observing run (O1) of the advanced LIGO detectors will be conducted. Based on the recent commissioning progress at the LIGO Hanford and Livingston sites, the gravitational wave detector range for a neutron star binary inspiral is expected to be of order 60 Mpc. We describe here our planning for an O1 search for coincidence between a LIGO gravitational wave detection and a gamma-ray signal from the Fermi Gamma-ray Burst Monitor. Such a coincidence would constitute measurement of an electromagnetic counterpart to a gravitational wave signal, with significant corresponding scientific benefits, including revealing the central engine powering the gamma-ray burst, enhanced confidence in the event as a genuine astrophysical detection, and a determination of the relative speed of the photon and graviton.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Robust parameter estimation for compact binaries with ground-based gravitational-wave observations using the LALInference software library

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, a neutron star black hole binary and a binary black hole, 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 parameter space.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Measuring the spin of black holes in binary systems using gravitational waves

Compact binary coalescences are the most promising sources of gravitational waves (GWs) for ground based detectors. Binary systems containing one or two spinning black holes are particularly interesting due to spin-orbit (and eventual spin-spin) interactions, and the opportunity of measuring spins directly through GW observations. In this letter we analyze simulated signals emitted by spinning binaries with several values of masses, spins, orientation, and signal-to-noise ratio. We find that spin magnitudes and tilt angles can be estimated to accuracy of a few percent for neutron star--black hole systems and $\sim$ 5-30% for black hole binaries. In contrast, the difference in the azimuth angles of the spins, which may be used to check if spins are locked into resonant configurations, cannot be constrained. We observe that the best performances are obtained when the line of sight is perpendicular to the system's total angular momentum, and that a sudden change of behavior occurs when a system is observed from angles such that the plane of the orbit can be seen both from above and below during the time the signal is in band. This study suggests that the measurement of black hole spin by means of GWs can be as precise as what can be obtained from X-ray binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Optimizing gravitational-wave searches for a population of coalescing binaries: Intrinsic parameters

We revisit the problem of searching for gravitational waves from inspiralling compact binaries in Gaussian coloured noise. For binaries with quasicircular orbits and non-precessing component spins, considering dominant mode emission only, if the intrinsic parameters of the binary are known then the optimal statistic for a single detector is the well-known two-phase matched filter. However, the matched filter signal-to-noise ratio is /not/ in general an optimal statistic for an astrophysical population of signals, since their distribution over the intrinsic parameters will almost certainly not mirror that of noise events, which is determined by the (Fisher) information metric. Instead, the optimal statistic for a given astrophysical distribution will be the Bayes factor, which we approximate using the output of a standard template matched filter search. We then quantify the possible improvement in number of signals detected for various populations of non-spinning binaries: for a distribution of signals uniformly distributed in volume and with component masses distributed uniformly over the range $1\leq m_{1,2}/M_\odot\leq 24$, $(m_1+m_2) /M_\odot\leq 25$ at fixed expected SNR, we find $\gtrsim 20\%$ more signals at a false alarm threshold of $10^{-6}\,$Hz in a single detector. The method may easily be generalized to binaries with non-precessing spins.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The First Two Years of Electromagnetic Follow-Up with Advanced LIGO and Virgo

We anticipate the first direct detections of gravitational waves (GWs) with Advanced LIGO and Virgo later this decade. Though this groundbreaking technical achievement will be its own reward, a still greater prize could be observations of compact binary mergers in both gravitational and electromagnetic channels simultaneously. During Advanced LIGO and Virgo's first two years of operation, 2015 through 2016, we expect the global GW detector array to improve in sensitivity and livetime and expand from two to three detectors. We model the detection rate and the sky localization accuracy for binary neutron star (BNS) mergers across this transition. We have analyzed a large, astrophysically motivated source population using real-time detection and sky localization codes and higher-latency parameter estimation codes that have been expressly built for operation in the Advanced LIGO/Virgo era. We show that for most BNS events the rapid sky localization, available about a minute after a detection, is as accurate as the full parameter estimation. We demonstrate that Advanced Virgo will play an important role in sky localization, even though it is anticipated to come online with only one-third as much sensitivity as the Advanced LIGO detectors. We find that the median 90% confidence region shrinks from ~500 square degrees in 2015 to ~200 square degrees in 2016. A few distinct scenarios for the first LIGO/Virgo detections emerge from our simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

TIGER: A data analysis pipeline for testing the strong-field dynamics of general relativity with gravitational wave signals from coalescing compact binaries

The direct detection of gravitational waves with upcoming second-generation gravitational wave detectors such as Advanced LIGO and Virgo will allow us to probe the genuinely strong-field dynamics of general relativity (GR) for the first time. We present a data analysis pipeline called TIGER (Test Infrastructure for GEneral Relativity), which is designed to utilize detections of compact binary coalescences to test GR in this regime. TIGER is a model-independent test of GR itself, in that it is not necessary to compare with any specific alternative theory. It performs Bayesian inference on two hypotheses: the GR hypothesis $\mathcal{H}_{\rm GR}$, and $\mathcal{H}_{\rm modGR}$, which states that one or more of the post-Newtonian coefficients in the waveform are not as predicted by GR. By the use of multiple sub-hypotheses of $\mathcal{H}_{\rm modGR}$, in each of which a different number of parameterized deformations of the GR phase are allowed, an arbitrarily large number of 'testing parameters' can be used without having to worry about a model being insufficiently parsimonious if the true number of extra parameters is in fact small. TIGER is well-suited to the regime where most sources have low signal-to-noise ratios, again through the use of these sub-hypotheses. Information from multiple sources can trivially be combined, leading to a stronger test. We focus on binary neutron star coalescences, for which sufficiently accurate waveform models are available that can be generated fast enough on a computer to be fit for use in Bayesian inference. We show that the pipeline is robust against a number of fundamental, astrophysical, and instrumental effects, such as differences between waveform approximants, a limited number of post-Newtonian phase contributions being known, the effects of neutron star spins and tidal deformability on the orbital motion, and instrumental calibration errors.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Reconstructing the sky location of gravitational-wave detected compact binary systems: methodology for testing and comparison

The problem of reconstructing the sky position of compact binary coalescences detected via gravitational waves is a central one for future observations with the ground-based network of gravitational-wave laser interferometers, such as Advanced LIGO and Advanced Virgo. Different techniques for sky localisation have been independently developed. They can be divided in two broad categories: fully coherent Bayesian techniques, which are high-latency and aimed at in-depth studies of all the parameters of a source, including sky position, and "triangulation-based" techniques, which exploit the data products from the search stage of the analysis to provide an almost real-time approximation of the posterior probability density function of the sky location of a detection candidate. These techniques have previously been applied to data collected during the last science runs of gravitational-wave detectors operating in the so-called initial configuration. Here, we develop and analyse methods for assessing the self-consistency of parameter estimation methods and carrying out fair comparisons between different algorithms, addressing issues of efficiency and optimality. These methods are general, and can be applied to parameter estimation problems other than sky localisation. We apply these methods to two existing sky localisation techniques representing the two above-mentioned categories, using a set of simulated inspiral-only signals from compact binary systems with total mass $\le 20\,M_\odot$ and non-spinning components. We compare the relative advantages and costs of the two techniques and show that sky location uncertainties are on average a factor $\approx 20$ smaller for fully coherent techniques than for the specific variant of the "triangulation-based" technique used during the last science runs, at the expense of a factor $\approx 1000$ longer processing time.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Techniques for targeted Fermi-GBM follow-up of gravitational-wave events

The Advanced LIGO and Advanced Virgo ground-based gravitational-wave detectors are projected to come online 2015-2016, reaching a final sensitivity sufficient to observe dozens of binary neutron star mergers per year by 2018. We present a fully-automated, targeted search strategy for prompt gamma-ray counterparts in offline Fermi-GBM data. The multi-detector method makes use of a detailed model response of the instrument, and benefits from time and sky location information derived from the gravitational-wave signal.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effect of calibration errors on Bayesian parameter estimation for gravitational wave signals from inspiral binary systems in the advanced detectors era

By 2015 the advanced versions of the gravitational-wave detectors Virgo and LIGO will be online. They will collect data in coincidence with enough sensitivity to potentially deliver multiple detections of gravitation waves from inspirals of compact-object binaries. This work is focused on understanding the effects introduced by uncertainties in the calibration of the interferometers. We consider plausible calibration errors based on estimates obtained during LIGO's fifth and Virgo's third science runs, which include frequency-dependent amplitude errors of $\sim 10%$ and frequency-dependent phase errors of $\sim 3$ degrees in each instrument. We quantify the consequences of such errors estimating the parameters of inspiraling binaries. We find that the systematics introduced by calibration errors on the inferred values of the chirp mass and mass ratio are smaller than 20% of the statistical measurement uncertainties in parameter estimation for 90% of signals in our mock catalog. Meanwhile, the calibration-induced systematics in the inferred sky location of the signal are smaller than $\sim 50%$ of the statistical uncertainty. We thus conclude that calibration-induced errors at this level are not a significant detriment to accurate parameter estimation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Estimating parameters of coalescing compact binaries with a detector network including LIGO Australia

One of the goals of gravitational-wave astronomy is simultaneous detection of gravitational-wave signals from merging compact-object binaries and the electromagnetic transients from these mergers. With the next generation of advanced ground-based gravitational wave detectors under construction, we examine the benefits of the proposed extension of the detector network to include a fourth site in Australia in addition to the network of Hanford, Livingston and Cascina sites. Using Bayesian parameter-estimation analyses of simulated gravitational-wave signals from a range of coalescing-binary locations and orientations, we study the improvement in parameter estimation. We find that an Australian detector can break degeneracies in several parameters; in particular, the localization of the source on the sky is improved by a factor of ~4, with more modest improvements in distance and binary inclination estimates. This enhanced ability to localize sources on the sky will be crucial in any search for electromagnetic counterparts to detected gravitational-wave signals.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Testing General Relativity using Bayesian model selection: Applications to observations of gravitational waves from compact binary systems

Second generation interferometric gravitational wave detectors, such as Advanced LIGO and Advanced Virgo, are expected to begin operation by 2015. Such instruments plan to reach sensitivities that will offer the unique possibility to test General Relativity in the dynamical, strong field regime and investigate departures from its predictions, in particular using the signal from coalescing binary systems. We introduce a statistical framework based on Bayesian model selection in which the Bayes factor between two competing hypotheses measures which theory is favored by the data. Probability density functions of the model parameters are then used to quantify the inference on individual parameters. We also develop a method to combine the information coming from multiple independent observations of gravitational waves, and show how much stronger inference could be. As an introduction and illustration of this framework - and a practical numerical implementation through the Monte Carlo integration technique of nested sampling - we apply it to gravitational waves from the inspiral phase of coalescing binary systems as predicted by General Relativity and a very simple alternative theory in which the graviton has a non-zero mass. This method can trivially (and should) be extended to more realistic and physically motivated theories.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Bayesian coherent analysis of in-spiral gravitational wave signals with a detector network

The present operation of the ground-based network of gravitational-wave laser interferometers in "enhanced" configuration brings the search for gravitational waves into a regime where detection is highly plausible. The development of techniques that allow us to discriminate a signal of astrophysical origin from instrumental artefacts in the interferometer data and to extract the full range of information are some of the primary goals of the current work. Here we report the details of a Bayesian approach to the problem of inference for gravitational wave observations using a network of instruments, for the computation of the Bayes factor between two hypotheses and the evaluation of the marginalised posterior density functions of the unknown model parameters. The numerical algorithm to tackle the notoriously difficult problem of the evaluation of large multi-dimensional integrals is based on a technique known as Nested Sampling, which provides an attractive alternative to more traditional Markov-chain Monte Carlo (MCMC) methods. We discuss the details of the implementation of this algorithm and its performance against a Gaussian model of the background noise, considering the specific case of the signal produced by the in-spiral of binary systems of black holes and/or neutron stars, although the method is completely general and can be applied to other classes of sources. We also demonstrate the utility of this approach by introducing a new coherence test to distinguish between the presence of a coherent signal of astrophysical origin in the data of multiple instruments and the presence of incoherent accidental artefacts, and the effects on the estimation of the source parameters as a function of the number of instruments in the network.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Mock LISA Data Challenges: from Challenge 3 to Challenge 4

The Mock LISA Data Challenges are a program to demonstrate LISA data-analysis capabilities and to encourage their development. Each round of challenges consists of one or more datasets containing simulated instrument noise and gravitational waves from sources of undisclosed parameters. Participants analyze the datasets and report best-fit solutions for the source parameters. Here we present the results of the third challenge, issued in Apr 2008, which demonstrated the positive recovery of signals from chirping Galactic binaries, from spinning supermassive--black-hole binaries (with optimal SNRs between ~ 10 and 2000), from simultaneous extreme-mass-ratio inspirals (SNRs of 10-50), from cosmic-string-cusp bursts (SNRs of 10-100), and from a relatively loud isotropic background with Omega_gw(f) ~ 10^-11, slightly below the LISA instrument noise.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Status of NINJA: the Numerical INJection Analysis project

The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Studying stellar binary systems with the Laser Interferometer Space Antenna using Delayed Rejection Markov chain Monte Carlo methods

Bayesian analysis of LISA data sets based on Markov chain Monte Carlo methods has been shown to be a challenging problem, in part due to the complicated structure of the likelihood function consisting of several isolated local maxima that dramatically reduces the efficiency of the sampling techniques. Here we introduce a new fully Markovian algorithm, a Delayed Rejection Metropolis-Hastings Markov chain Monte Carlo method, to efficiently explore these kind of structures and we demonstrate its performance on selected LISA data sets containing a known number of stellar-mass binary signals embedded in Gaussian stationary noise.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Testing gravitational-wave searches with numerical relativity waveforms: Results from the first Numerical INJection Analysis (NINJA) project

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Detecting gravitational radiation from neutron stars using a six-parameter adaptive MCMC method

We present a Markov chain Monte Carlo technique for detecting gravitational radiation from a neutron star in laser interferometer data. The algorithm can estimate up to six unknown parameters of the target, including the rotation frequency and frequency derivative, using reparametrization, delayed rejection and simulated annealing. We highlight how a simple extension of the method, distributed over multiple computer processors, will allow for a search over a narrow frequency band. The ultimate goal of this research is to search for sources at a known locations, but uncertain spin parameters, such as may be found in SN1987A.

0 citations0 reviewsNo code linkedOpen access