Source author record

Jonathan R. Gair

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

ResearcherUnclaimed source record

gr-qc astro-ph.CO astro-ph.HE astro-ph.IM astro-ph.GA hep-th astro-ph Computation Applications Machine Learning Methodology

Catalog footprint

What is connected

56works

11topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2026arXiv

Trust the process: mapping data-driven reconstructions to informed models using stochastic processes

Gravitational-wave astronomy has entered a regime where it can extract information about the population properties of the observed binary black holes. The steep increase in the number of detections will offer deeper insights, but it will also significantly raise the computational cost of testing multiple models. To address this challenge, we propose a procedure that first performs a non-parametric (data-driven) reconstruction of the underlying distribution, and then remaps these results onto a posterior for the parameters of a parametric (informed) model. The computational cost is primarily absorbed by the initial non-parametric step, while the remapping procedure is both significantly easier to perform and computationally cheaper. In addition to yielding the posterior distribution of the model parameters, this method also provides a measure of the model's goodness-of-fit, opening for a new quantitative comparison across models.

preprint2020arXiv

2 Fast 2 Fiducial: Gaussian processes for the interpolation and marginalization of waveform error in extreme-mass-ratio-inspiral parameter estimation

A number of open problems hinder our present ability to extract scientific information from data that will be gathered by the near-future gravitational-wave mission LISA. Many of these relate to the modeling, detection and characterization of signals from binary inspirals with an extreme component-mass ratio of $\lesssim10^{-4}$. In this paper, we draw attention to the issue of systematic error in parameter estimation due to the use of fast but approximate waveform models; this is found to be relevant for extreme-mass-ratio inspirals even in the case of waveforms with $\gtrsim90\%$ overlap accuracy and moderate ($\gtrsim30$) signal-to-noise ratios. A scheme that uses Gaussian processes to interpolate and marginalize over waveform error is adapted and investigated as a possible precursor solution to this problem. Several new methodological results are obtained, and the viability of the technique is successfully demonstrated on a three-parameter example in the setting of the LISA Data Challenge.

preprint2020arXiv

Constraining the spin parameter of near-extremal black holes using LISA

We describe a model that generates first order adiabatic EMRI waveforms for quasi-circular equatorial inspirals of compact objects into rapidly rotating (near-extremal) black holes. Using our model, we show that LISA could measure the spin parameter of near-extremal black holes (for $a \gtrsim 0.9999$) with extraordinary precision, $\sim$ 3-4 orders of magnitude better than for moderate spins, $a \sim 0.9$. Such spin measurements would be one of the tightest measurements of an astrophysical parameter within a gravitational wave context. Our results are primarily based off a Fisher matrix analysis, but are verified using both frequentest and Bayesian techniques. We present analytical arguments that explain these high spin precision measurements. The high precision arises from the spin dependence of the radial inspiral evolution, which is dominated by geodesic properties of the secondary orbit, rather than radiation reaction. High precision measurements are only possible if we observe the exponential damping of the signal that is characteristic of the near-horizon regime of near-extremal inspirals. Our results demonstrate that, if such black holes exist, LISA would be able to successfully identify rapidly rotating black holes up to $a = 1-10^{-9}$ , far past the Thorne limit of $a = 0.998$.

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.

preprint2020arXiv

Transition from Inspiral to Plunge: A Complete Near-Extremal Trajectory and Associated Waveform

We extend the Ori and Thorne (OT) procedure to compute the transition from an adiabatic inspiral into a geodesic plunge for any spin, with emphasis on near-extremal ones. Our analysis revisits the validity of the approximations made in OT. In particular, we discuss possible effects coming from eccentricity and non-geodesic past-history of the orbital evolution. We find three different scaling regimes according to whether the mass ratio is much smaller, of the same order or much larger than the near extremal parameter describing how fast the primary black hole rotates. Eccentricity and non-geodesic past-history corrections are always sub-leading, indicating that the quasi-circular approximation applies throughout the transition regime. However, we show that the OT assumption that the energy and angular momentum evolve linearly with proper time must be modified in the near-extremal regime. Using our transition equations, we describe an algorithm to compute the full worldline in proper time for an extreme mass ratio inspiral (EMRI) and the resultant gravitational waveform in the high spin limit.

preprint2016arXiv

Fast methods for training Gaussian processes on large data sets

Gaussian process regression (GPR) is a non-parametric Bayesian technique for interpolating or fitting data. The main barrier to further uptake of this powerful tool rests in the computational costs associated with the matrices which arise when dealing with large data sets. Here, we derive some simple results which we have found useful for speeding up the learning stage in the GPR algorithm, and especially for performing Bayesian model comparison between different covariance functions. We apply our techniques to both synthetic and real data and quantify the speed-up relative to using nested sampling to numerically evaluate model evidences.

preprint2016arXiv

Improving gravitational-wave parameter estimation using Gaussian process regression

Folding uncertainty in theoretical models into Bayesian parameter estimation is necessary in order to make reliable inferences. A general means of achieving this is by marginalizing over model uncertainty using a prior distribution constructed using Gaussian process regression (GPR). As an example, we apply this technique to the measurement of chirp mass using (simulated) gravitational-wave signals from binary black holes that could be observed using advanced-era gravitational-wave detectors. Unless properly accounted for, uncertainty in the gravitational-wave templates could be the dominant source of error in studies of these systems. We explain our approach in detail and provide proofs of various features of the method, including the limiting behavior for high signal-to-noise, where systematic model uncertainties dominate over noise errors. We find that the marginalized likelihood constructed via GPR offers a significant improvement in parameter estimation over the standard, uncorrected likelihood both in our simple one-dimensional study, and theoretically in general. We also examine the dependence of the method on the size of training set used in the GPR; on the form of covariance function adopted for the GPR, and on changes to the detector noise power spectral density.

preprint2016arXiv

Tunable compression of template banks for fast gravitational-wave detection and localisation

One strategy for reducing the online computational cost of matched-filter searches for gravitational waves is to introduce a compressed basis for the waveform template bank in a grid-based search. In this paper, we propose and investigate several tunable compression schemes for a general template bank. Through offline compression, such schemes are shown to yield faster detection and localisation of signals, along with moderately improved sensitivity and accuracy over coarsened banks at the same level of computational cost. This is potentially useful for any search involving template banks, and especially in the analysis of data from future space-based detectors such as eLISA, for which online grid searches are difficult due to the long-duration waveforms and large parameter spaces.

preprint2015arXiv

Detection of eccentric supermassive black hole binaries with pulsar timing arrays: Signal-to-noise ratio calculations

We present a detailed analysis of the expected signal-to-noise ratios of supermassive black hole binaries on eccentric orbits observed by pulsar timing arrays. We derive several analytical relations that extend the results of Peters and Mathews [Phys. Rev. D 131, 435 (1963)] to quantify the impact of eccentricity in the detection of single resolvable binaries in the pulsar timing array band. We present ready-to-use expressions to compute the increase/loss in signal-to-noise ratio of eccentric single resolvable sources whose dominant harmonic is located in the low/high frequency sensitivity regime of pulsar timing arrays. Building upon the work of Phinney (arXiv:astro-ph/0108028) and Enoki and Nagashima [Prog. Theor. Phys. 117, 241 (2007)], we present an analytical framework that enables the construction of rapid spectra for a stochastic gravitational-wave background generated by a cosmological population of eccentric sources. We confirm previous findings which indicate that, relative to a population of quasicircular binaries, the strain of a stochastic, isotropic gravitational-wave background generated by a cosmological population of eccentric binaries will be suppressed in the frequency band of pulsar timing arrays. We quantify this effect in terms of signal-to-noise ratios in a pulsar timing array.

preprint2015arXiv

Estimating the sensitivity of pulsar timing arrays

The sensitivity curve of a canonical pulsar timing array is calculated for two types of source: a monochromatic wave and a stochastic background. These calculations are performed in both a Bayesian and frequentist framework, using both analytical and numerical methods. These calculations are used to clarify the interpretation of the sensitivity curves and to illustrate the sometimes overlooked fact that the sensitivity curve depends not only on the properties of the pulse time-of-arrival data set but also on the properties of the source being observed. The Bayesian and frequentist frameworks were found to give consistent results and the analytic and numerical calculations were also found to be in good agreement.

preprint2015arXiv

European Pulsar Timing Array Limits On An Isotropic Stochastic Gravitational-Wave Background

We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar dataset spanning 18 yr of observations from the 2015 European Pulsar Timing Array data release. Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar System ephemeris errors, obtaining a robust 95$\%$ upper limit on the dimensionless strain amplitude $A$ of the background of $A<3.0\times 10^{-15}$ at a reference frequency of $1\mathrm{yr^{-1}}$ and a spectral index of $13/3$, corresponding to a background from inspiralling super-massive black hole binaries, constraining the GW energy density to $Ω_\mathrm{gw}(f)h^2 < 1.1\times10^{-9}$ at 2.8 nHz. We also present limits on the correlated power spectrum at a series of discrete frequencies, and show that our sensitivity to a fiducial isotropic GWB is highest at a frequency of $\sim 5\times10^{-9}$~Hz. Finally we discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95$\%$ upper limits on the string tension, $Gμ/c^2$, characterising a background produced by a cosmic string network for a set of possible scenarios, and for a stochastic relic GWB. For a Nambu-Goto field theory cosmic string network, we set a limit $Gμ/c^2<1.3\times10^{-7}$, identical to that set by the {\it Planck} Collaboration, when combining {\it Planck} and high-$\ell$ Cosmic Microwave Background data from other experiments. For a stochastic relic background we set a limit of $Ω^\mathrm{relic}_\mathrm{gw}(f)h^2<1.2 \times10^{-9}$, a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array.

preprint2015arXiv

European Pulsar Timing Array Limits on Continuous Gravitational Waves from Individual Supermassive Black Hole Binaries

We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest EPTA dataset, which consists of ultra-precise timing data on 41 millisecond pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algorithms show evidence for the presence of such a CGW signal, indicating that the data are best described by pulsar and radiometer noise only. Depending on the adopted detection algorithm, the 95\% upper limit on the sky-averaged strain amplitude lies in the range $6\times 10^{-15}<A<1.5\times10^{-14}$ at $5{\rm nHz}<f<7{\rm nHz}$. This limit varies by a factor of five, depending on the assumed source position, and the most constraining limit is achieved towards the positions of the most sensitive pulsars in the timing array. The most robust upper limit -- obtained via a full Bayesian analysis searching simultaneously over the signal and pulsar noise on the subset of ours six best pulsars -- is $A\approx10^{-14}$. These limits, the most stringent to date at $f<10{\rm nHz}$, exclude the presence of sub-centiparsec binaries with chirp mass $\cal{M}_c>10^9$M$_\odot$ out to a distance of about 25Mpc, and with $\cal{M}_c>10^{10}$M$_\odot$ out to a distance of about 1Gpc ($z\approx0.2$). We show that state-of-the-art SMBHB population models predict $<1\%$ probability of detecting a CGW with the current EPTA dataset, consistent with the reported non-detection. We stress, however, that PTA limits on individual CGW have improved by almost an order of magnitude in the last five years. The continuing advances in pulsar timing data acquisition and analysis techniques will allow for strong astrophysical constraints on the population of nearby SMBHBs in the coming years.

preprint2015arXiv

Improved analytic extreme-mass-ratio inspiral model for scoping out eLISA data analysis

The space-based gravitational-wave detector eLISA has been selected as the ESA L3 mission, and the mission design will be finalised by the end of this decade. To prepare for mission formulation over the next few years, several outstanding and urgent questions in data analysis will be addressed using mock data challenges, informed by instrument measurements from the LISA Pathfinder satellite launching at the end of 2015. These data challenges will require accurate and computationally affordable waveform models for anticipated sources such as the extreme-mass-ratio inspirals (EMRIs) of stellar-mass compact objects into massive black holes. Previous data challenges have made use of the well-known analytic EMRI waveforms of Barack and Cutler, which are extremely quick to generate but dephase relative to more accurate waveforms within hours, due to their mismatched radial, polar and azimuthal frequencies. In this paper, we describe an augmented Barack-Cutler model that uses a frequency map to the correct Kerr frequencies, along with updated evolution equations and a simple fit to a more accurate model. The augmented waveforms stay in phase for months and may be generated with virtually no additional computational cost.

preprint2015arXiv

Likelihood smoothing using gravitational wave surrogate models

Likelihood surfaces in the parameter space of gravitational wave signals can contain many secondary maxima, which can prevent search algorithms from finding the global peak and correctly mapping the distribution. Traditional schemes to mitigate this problem maintain the number of secondary maxima and thus retain the possibility that the global maximum will remain undiscovered. By contrast, the recently proposed technique of likelihood transform can modify the structure of the likelihood surface to reduce its complexity. We present a practical method to carry out a likelihood transform using a Gaussian smoothing kernel, utilising gravitational wave surrogate models to perform the smoothing operation analytically. We demonstrate the approach with Newtonian and post-Newtonian waveform models for an inspiralling circular compact binary.

preprint2015arXiv

Mapping gravitational-wave backgrounds in modified theories of gravity using pulsar timing arrays

We extend our previous work on applying CMB techniques to the mapping of gravitational-wave backgrounds to backgrounds which have non-GR polarisations. Our analysis and results are presented in the context of pulsar-timing array observations, but the overarching methods are general, and can be easily applied to LIGO or eLISA observations using appropriately modified response functions. Analytic expressions for the pulsar-timing response to gravitational waves with non-GR polarisation are given for each mode of a spin-weighted spherical-harmonic decomposition of the background, which permit the signal to be mapped across the sky to any desired resolution. We also derive the pulsar-timing overlap reduction functions for the various non-GR polarisations, finding analytic forms for anisotropic backgrounds with scalar-transverse ("breathing") and vector-longitudinal polarisations, and a semi-analytic form for scalar-longitudinal backgrounds. Our results indicate that pulsar-timing observations will be completely insensitive to scalar-transverse mode anisotropies in the polarisation amplitude beyond dipole, and anisotropies in the power beyond quadrupole. Analogously to our previous findings that pulsar-timing observations lack sensitivity to tensor-curl modes for a transverse-traceless tensor background, we also find insensitivity to vector-curl modes for a vector-longitudinal background.

preprint2015arXiv

Quantifying and mitigating bias in inference on gravitational wave source populations

When using incorrect or inaccurate signal models to perform parameter estimation on a gravitational wave signal, biased parameter estimates will in general be obtained. For a single event this bias may be consistent with the posterior, but when considering a population of events this bias becomes evident as a sag below the expected diagonal line of the P-P plot showing the fraction of signals found within a certain significance level versus that significance level. It would be hoped that recently proposed techniques for accounting for model uncertainties in parameter estimation would, to some extent, alleviate this problem. Here we demonstrate that this is indeed the case. We derive an analytic approximation to the P-P plot obtained when using an incorrect signal model to perform parameter estimation. This approximation is valid in the limit of high signal-to-noise ratio and nearly correct waveform models. We show how the P-P plot changes if a Gaussian process likelihood that allows for model errors is used to analyse the data. We demonstrate analytically and using numerical simulations that the bias is always reduced in this way. These results provide a way to quantify bias in inference on populations and demonstrate the importance of utilising methods to mitigate this bias.

preprint2015arXiv

Testing the "no-hair" property of black holes with X-ray observations of accretion disks

Accretion disks around black holes radiate a significant fraction of the rest mass of the accreting material in the form of thermal radiation from within a few gravitational radii of the black hole ($ r \lesssim 20 G M / c^{2}$). In addition, the accreting matter may also be illuminated by hard X-rays from the surrounding plasma which adds fluorescent transition lines to the emission. This radiation is emitted by matter moving along geodesics in the metric, therefore the strong Doppler and gravitational redshifts observed in the emission encode information about the strong gravitational field around the black hole. In this paper the possibility of using the X-ray emission as a strong field test of General Relativity is explored by calculating the spectra for both the transition line and thermal emission from a thin accretion disk in a series of parametrically deformed Kerr metrics. In addition the possibility of constraining a number of known black hole spacetimes in alternative theories of gravity is considered.

preprint2014arXiv

Counting And Confusion: Bayesian Rate Estimation With Multiple Populations

We show how to obtain a Bayesian estimate of the rates or numbers of signal and background events from a set of events when the shapes of the signal and background distributions are known, can be estimated, or approximated; our method works well even if the foreground and background event distributions overlap significantly and the nature of any individual event cannot be determined with any certainty. We give examples of determining the rates of gravitational-wave events in the presence of background triggers from a template bank when noise parameters are known and/or can be fit from the trigger data. We also give an example of determining globular-cluster shape, location, and density from an observation of a stellar field that contains a non-uniform background density of stars superimposed on the cluster stars.

preprint2014arXiv

Electromagnetic signatures of far-field gravitational radiation in the 1+3 approach

Gravitational waves from astrophysical sources can interact with background electromagnetic fields, giving rise to distinctive and potentially detectable electromagnetic signatures. In this paper, we study such interactions for far-field gravitational radiation using the 1+3 approach to relativity. Linearised equations for the electromagnetic field on perturbed Minkowski space are derived and solved analytically. The inverse Gertsenshtein conversion of gravitational waves in a static electromagnetic field is rederived, and the resultant electromagnetic radiation is shown to be significant for highly magnetised pulsars in compact binary systems. We also obtain a variety of nonlinear interference effects for interacting gravitational and electromagnetic waves, although wave-wave resonances previously described in the literature are absent when the electric-magnetic self-interaction is taken into account. The fluctuation and amplification of electromagnetic energy flux as the gravitational wave strength increases towards the gravitational-electromagnetic frequency ratio is a possible signature of gravitational radiation from extended astrophysical sources.

preprint2014arXiv

Mapping gravitational-wave backgrounds using methods from CMB analysis: Application to pulsar timing arrays

We describe an alternative approach to the analysis of gravitational-wave backgrounds, based on the formalism used to characterise the polarisation of the cosmic microwave background. In contrast to standard analyses, this approach makes no assumptions about the nature of the background and so has the potential to reveal much more about the physical processes that generated it. An arbitrary background can be decomposed into modes whose angular dependence on the sky is given by gradients and curls of spherical harmonics. We derive the pulsar timing overlap reduction functions for the individual modes, which are given by simple combinations of spherical harmonics evaluated at the pulsar locations. We show how these can be used to recover the components of an arbitrary background, giving explicit results for both isotropic and anisotropic uncorrelated backgrounds. We also find that the response of a pulsar timing array to curl modes is identically zero, so half of the gravitational-wave sky will never be observed using pulsar timing, no matter how many pulsars are included in the array. An isotropic, unpolarised and uncorrelated background can be accurately represented using only three modes, and so a search of this type will be only slightly more complicated than the standard cross-correlation search using the Hellings and Downs overlap reduction function. However, by measuring the components of individual modes of the background and checking for consistency with isotropy, this approach has the potential to reveal much more information. Each individual mode on its own describes a background that is correlated between different points on the sky. A measurement of the components that indicates the presence of correlations in the background on large angular scales would suggest startling new physics.

preprint2014arXiv

Novel Method for Incorporating Model Uncertainties into Gravitational Wave Parameter Estimates

Posterior distributions on parameters computed from experimental data using Bayesian techniques are only as accurate as the models used to construct them. In many applications these models are incomplete, which both reduces the prospects of detection and leads to a systematic error in the parameter estimates. In the analysis of data from gravitational wave detectors, for example, accurate waveform templates can be computed using numerical methods, but the prohibitive cost of these simulations means this can only be done for a small handful of parameters. In this work a novel method to fold model uncertainties into data analysis is proposed; the waveform uncertainty is analytically marginalised over using with a prior distribution constructed by using Gaussian process regression to interpolate the waveform difference from a small training set of accurate templates. The method is well motivated, easy to implement, and no more computationally expensive than standard techniques. The new method is shown to perform extremely well when applied to a toy problem. While we use the application to gravitational wave data analysis to motivate and illustrate the technique, it can be applied in any context where model uncertainties exist.

preprint2014arXiv

Research Update on Extreme-Mass-Ratio Inspirals

The inspirals of stellar-mass mass compact objects into massive black holes in the centres of galaxies are one of the most important sources of gravitational radiation for space-based detectors like LISA or eLISA. These extreme-mass-ratio inspirals (EMRIs) will enable an ambitious research program with implications for astrophysics, cosmology, and fundamental physics. This article is a summary of the talks delivered at the plenary session on EMRIs at the 10th International LISA Symposium. It contains research updates on the following topics: astrophysics of EMRIs; EMRI science potential; and EMRI modeling.

preprint2014arXiv

Self-forced evolutions of an implicit rotating source: A natural framework to model comparable and intermediate mass-ratio systems from inspiral through ringdown

We develop a waveform model to describe the inspiral, merger and ringdown of binary systems with comparable and intermediate mass-ratios. This model incorporates first-order conservative self-force corrections to the energy and angular momentum, which are valid in the strong-field regime [1]. We model the radiative part of the self-force by deriving second-order radiative corrections to the energy flux. These corrections are obtained by minimizing the phase discrepancy between our self-force model and the effective one body model [2, 3] for a variety of mass-ratios. We show that our model performs substantially better than post-Newtonian approximants currently used to model neutron star-black hole mergers from early inspiral to the innermost stable circular orbit. In order to match the late inspiral evolution onto the plunge regime, we extend the 'transition phase' developed by Ori and Thorne [4] by including finite mass-ratio corrections and modelling the orbital phase evolution using an implicit rotating source [5]. We explicitly show that the implicit rotating source approach provides a natural transition from late-time radiation to ringdown that is equivalent to ringdown waveform modelling based on a sum of quasinormal modes.

preprint2014arXiv

The scientific potential of space-based gravitational wave detectors

The millihertz gravitational wave band can only be accessed with a space-based interferometer, but it is one of the richest in potential sources. Observations in this band have amazing scientific potential. The mergers between massive black holes with mass in the range 10 thousand to 10 million solar masses, which are expected to occur following the mergers of their host galaxies, produce strong millihertz gravitational radiation. Observations of these systems will trace the hierarchical assembly of structure in the Universe in a mass range that is very difficult to probe electromagnetically. Stellar mass compact objects falling into such black holes in the centres of galaxies generate detectable gravitational radiation for several years prior to the final plunge and merger with the central black hole. Measurements of these systems offer an unprecedented opportunity to probe the predictions of general relativity in the strong-field and dynamical regime. Millihertz gravitational waves are also generated by millions of ultra-compact binaries in the Milky Way, providing a new way to probe galactic stellar populations. ESA has recognised this great scientific potential by selecting The Gravitational Universe as its theme for the L3 large satellite mission, scheduled for launch in ~2034. In this article we will review the likely sources for millihertz gravitational wave detectors and describe the wide applications that observations of these sources could have for astrophysics, cosmology and fundamental physics.

preprint2013arXiv

Gravitational wave parameter estimation with compressed likelihood evaluations

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.

preprint2013arXiv

Searching For Anisotropic Gravitational-wave Backgrounds Using Pulsar Timing Arrays

We present the results of simulated injections testing the first Bayesian search-pipeline capable of investigating the angular-structure of a gravitational-wave (GW) background influencing pulsar signals. A stochastic background of GWs from the incoherent superposition of many inspiraling supermassive black hole binaries at nHz frequencies is likely to be the dominant GW signal detectable by pulsar timing arrays (PTAs). Even though one might expect a background composed of a high-redshift cosmological population of sources to be fairly isotropic, deviations from isotropy may be indicative of local GW hotspots or some form of continuous anisotropy in the angular-distribution of GW-power. A GWB induces time-of-arrival deviations in pulsar signals which are correlated between separated pulsars. In an isotropic background this cross-correlation follows a distinctive relationship, known as the Hellings and Downs curve, that depends only on the angular separation of the pulsars. If the background is anisotropic, the cross-correlation is different, but predictable, and also depends on the absolute position of the pulsars. By simulating datasets containing GWBs with various anisotropic configurations, we have explored the prospects for constraining anisotropy using near future data. We find that at moderate to high signal to noise ratio the assumption of isotropy is no longer an appropriate description of the simulated background. Furthermore, we can recover the nature of the injected anisotropy in a Bayesian parameter-estimation search, and propose a prior on the anisotropy search-space motivated by the physicality of the implied distribution of sources.

preprint2013arXiv

Testing General Relativity with Low-Frequency, Space-Based Gravitational-Wave Detectors

We review the tests of general relativity that will become possible with space-based gravitational-wave detectors operating in the ~0.01mHz - 1Hz low-frequency band. The fundamental aspects of gravitation that can be tested include the presence of additional gravitational fields other than the metric; the number and tensorial nature of gravitational-wave polarization states; the velocity of propagation of gravitational waves; the binding energy and gravitational-wave radiation of binaries, and therefore the time evolution of binary inspirals; the strength and shape of the waves emitted from binary mergers and ringdowns; the true nature of astrophysical black holes; and much more. The strength of this science alone calls for the swift implementation of a space-based detector; the remarkable richness of astrophysics, astronomy, and cosmology in the low-frequency gravitational-wave band make the case even stronger.

preprint2012arXiv

Constraining Gravity with LISA Detections of Binaries

General Relativity (GR) describes gravitation well at the energy scales which we have so far been able to achieve or detect. However, we do not know whether GR is behind the physics governing stronger gravitational field regimes, such as near neutron stars or massive black-holes (MBHs). Gravitational-wave (GW) astronomy is a promising tool to test and validate GR and/or potential alternative theories of gravity. The information that a GW waveform carries not only will allow us to map the strong gravitational field of its source, but also determine the theory of gravity ruling its dynamics. In this work, we explore the extent to which we could distinguish between GR and other theories of gravity through the detection of low-frequency GWs from extreme-mass-ratio inspirals (EMRIs) and, in particular, we focus on dynamical Chern-Simons modified gravity (DCSMG). To that end, we develop a framework that enables us, for the first time, to perform a parameter estimation analysis for EMRIs in DCSMG. Our model is described by a 15-dimensional parameter space, that includes the Chern-Simons (CS) parameter which characterises the deviation between the two theories, and our analysis is based on Fisher information matrix techniques together with a (maximum-mismatch) criterion to assess the validity of our results. In our analysis, we study a 5-dimensional parameter space, finding that a GW detector like the Laser Interferometer Space Antenna (LISA) or eLISA (evolved LISA) should be able to discriminate between GR and DCSMG with fractional errors below 5%, and hence place bounds four orders of magnitude better than current Solar System bounds.

preprint2012arXiv

Cosmology with the lights off: Standard sirens in the Einstein Telescope era

We explore the prospects for constraining cosmology using gravitational-wave (GW) observations of neutron-star binaries by the proposed Einstein Telescope (ET), exploiting the narrowness of the neutron-star mass function. Double neutron-star (DNS) binaries are expected to be one of the first sources detected after "first-light" of Advanced LIGO and are expected to be detected at a rate of a few tens per year in the advanced era. However the proposed ET could catalog tens of thousands per year. Combining the measured source redshift distributions with GW-network distance determinations will permit not only the precision measurement of background cosmological parameters, but will provide an insight into the astrophysical properties of these DNS systems. Of particular interest will be to probe the distribution of delay times between DNS-binary creation and subsequent merger, as well as the evolution of the star-formation rate density within ET's detection horizon. Keeping H_0, Ω_{m,0} and Ω_{Λ,0} fixed and investigating the precision with which the dark-energy equation-of-state parameters could be recovered, we found that with 10^5 detected DNS binaries we could constrain these parameters to an accuracy similar to forecasted constraints from future CMB+BAO+SNIa measurements. Furthermore, modeling the merger delay-time distribution as a power-law, and the star-formation rate (SFR) density as a parametrized version of the Porciani and Madau SF2 model, we find that the associated astrophysical parameters are constrained to within ~ 10%. All parameter precisions scaled as 1/sqrt(N), where N is the number of cataloged detections. We also investigated how precisions varied with the intrinsic underlying properties of the Universe and with the distance reach of the network (which may be affected by the low-frequency cutoff of the detector).

preprint2012arXiv

Extreme-mass-ratio bursts from the Galactic Centre

An extreme-mass-ratio burst (EMRB) is a gravitational wave signal emitted when a compact object passes through periapsis on a highly eccentric orbit about a much more massive body, in our case a stellar mass object about the 4.31 \times 10^6 M_sol massive black hole (MBH) in the Galactic Centre. We investigate how EMRBs could constrain the parameters of the Galaxy's MBH. EMRBs should be detectable if the periapsis is r_p < 65 r_g for a μ= 10 M_sol orbiting object, where r_g = GM/c^2 is the gravitational radius. The signal-to-noise ratio ρscales like log(ρ) = -2.7 log(r_p/r_g) + log(μ/M_sol) + 4.9. For periapses smaller than ~ 10 r_g, EMRBs can be informative, providing good constraints on both the MBH's mass and spin.

preprint2012arXiv

Hubble without the Hubble: cosmology using advanced gravitational-wave detectors alone

We investigate a novel approach to measuring the Hubble constant using gravitational-wave (GW) signals from compact binaries by exploiting the narrowness of the distribution of masses of the underlying neutron-star population. Gravitational-wave observations with a network of detectors will permit a direct, independent measurement of the distance to the source systems. If the redshift of the source is known, these inspiraling double-neutron-star binary systems can be used as standard sirens to extract cosmological information. Unfortunately, the redshift and the system chirp mass are degenerate in GW observations. Thus, most previous work has assumed that the source redshift is obtained from electromagnetic counterparts. In this paper, we explore what we can learn about the background cosmology and the mass distribution of neutron stars from the set of neutron-star (NS) mergers detected by such a network. We use a Bayesian formalism to analyze catalogs of NS-NS inspiral detections. We find that it is possible to constrain the Hubble constant, H_0, and the parameters of the NS mass function using gravitational-wave data alone, without relying on electromagnetic counterparts. Under reasonable assumptions, we will be able to determine H_0 to +/- 10% using ~100 observations, provided the Gaussian half-width of the underlying double NS mass distribution is less than 0.04 solar masses. The expected precision depends linearly on the intrinsic width of the NS mass function, but has only a weak dependence on H_0 near the default parameter values. Finally, we consider what happens if, for some fraction of our data catalog, we have an electromagnetically measured redshift. The detection, and cataloging, of these compact-object mergers will allow precision astronomy, and provide a determination of H_0 which is independent of the local distance scale.

preprint2012arXiv

Importance of including small body spin effects in the modelling of intermediate mass-ratio inspirals. II Accurate parameter extraction of strong sources using higher-order spin effects

We improve the numerical kludge waveform model introduced in [1] in two ways. We extend the equations of motion for spinning black hole binaries derived by Saijo et al. [2] using spin-orbit and spin-spin couplings taken from perturbative and post-Newtonian (PN) calculations at the highest order available. We also include first-order conservative self-force corrections for spin-orbit and spin-spin couplings, which are derived by comparison to PN results. We generate the inspiral evolution using fluxes that include the most recent calculations of small body spin corrections, spin-spin and spin-orbit couplings and higher-order fits to solutions of the Teukolsky equation. Using a simplified version of this model in [1], we found that small body spin effects could be measured through gravitational wave observations from intermediate-mass ratio inspirals (IMRIs) with mass ratio eta ~ 0.001, when both binary components are rapidly rotating. In this paper we study in detail how the spin of the small/big body affects parameter measurement using a variety of mass and spin combinations for typical IMRIs sources. We find that for IMRI events of a moderately rotating intermediate mass black hole (IMBH) of ten thousand solar masses, and a rapidly rotating central supermassive black hole (SMBH) of one million solar masses, gravitational wave observations made with LISA at a fixed signal-to-noise ratio (SNR) of 1000 will be able to determine the inspiralling IMBH mass, the central SMBH mass, the SMBH spin magnitude, and the IMBH spin magnitude to within fractional errors of ~0.001, 0.001, 0.0001, and 9%, respectively. LISA can also determine the location of the source in the sky and the SMBH spin orientation to within ~0.0001 steradians. We show that by including conservative corrections up to 2.5PN order, systematic errors no longer dominate over statistical errors for IMRIs with typical SNR ~1000.

preprint2012arXiv

Linearized f(R) Gravity: Gravitational Radiation & Solar System Tests

We investigate the linearized form of metric f(R)-gravity, assuming that f(R) is analytic about R = 0 so it may be expanded as f(R) = R + a_2 R^2/2 + ... . Gravitational radiation is modified, admitting an extra mode of oscillation, that of the Ricci scalar. We derive an effective energy-momentum tensor for the radiation. We also present weak-field metrics for simple sources. These are distinct from the equivalent Kerr (or Schwarzschild) forms. We apply the metrics to tests that could constrain f(R). We show that light deflection experiments cannot distinguish f(R)-gravity from general relativity as both have an effective post-Newtonian parameter γ= 1. We find that planetary precession rates are enhanced relative to general relativity; from the orbit of Mercury we derive the bound |a_2| < 1.2 \times 10^18 m^2. Gravitational wave astronomy may be more useful: considering the phase of a gravitational waveform we estimate deviations from general relativity could be measurable for an extreme-mass-ratio inspiral about a 10^6 M_sol black hole if |a_2| > 10^17 m^2, assuming that the weak-field metric of the black hole coincides with that of a point mass. However Eot-Wash experiments provide the strictest bound |a_2| < 2 \times 10^-9 m^2. Although the astronomical bounds are weaker, they are still of interest in the case that the effective form of f(R) is modified in different regions, perhaps through the chameleon mechanism. Assuming the laboratory bound is universal, we conclude that the propagating Ricci scalar mode cannot be excited by astrophysical sources.

preprint2012arXiv

Testing Chern-Simons Modified Gravity with Gravitational-Wave Detections of Extreme-Mass-Ratio Binaries

[abridged] The detection of gravitational waves from extreme-mass-ratio (EMRI) binaries, comprising a stellar-mass compact object orbiting around a massive black hole, is one of the main targets for low-frequency gravitational-wave detectors in space, like the Laser Interferometer Space Antenna (LISA or eLISA/NGO). The long-duration gravitational-waveforms emitted by such systems encode the structure of the strong field region of the massive black hole, in which the inspiral occurs. The detection and analysis of EMRIs will therefore allow us to study the geometry of massive black holes and determine whether their nature is as predicted by General Relativity and even to test whether General Relativity is the correct theory to describe the dynamics of these systems. To achieve this, EMRI modeling in alternative theories of gravity is required to describe the generation of gravitational waves. In this paper, we explore to what extent EMRI observations with LISA or eLISA/NGO might be able to distinguish between General Relativity and a particular modification of it, known as Dynamical Chern-Simons Modified Gravity. Our analysis is based on a parameter estimation study that uses approximate gravitational waveforms obtained via a radiative-adiabatic method and is restricted to a five-dimensional subspace of the EMRI configuration space. This includes a Chern-Simons parameter that controls the strength of gravitational deviations from General Relativity. We find that, if Dynamical Chern-Simons Modified Gravity is the correct theory, an observatory like LISA or even eLISA/NGO should be able to measure the Chern-Simons parameter with fractional errors below 5%. If General Relativity is the true theory, these observatories should put bounds on this parameter at the level xi^(1/4) < 10^4 km, which is four orders of magnitude better than current Solar System bounds.

preprint2012arXiv

Testing Chern-Simons modified gravity with observations of extreme-mass-ratio binaries

Extreme-Mass-Ratio Inspirals (EMRIs) are one of the most promising sources of gravitational waves (GWs) for space-based detectors like the Laser Interferometer Space Antenna (LISA). EMRIs consist of a compact stellar object orbiting around a massive black hole (MBH). Since EMRI signals are expected to be long lasting (containing of the order of hundred thousand cycles), they will encode the structure of the MBH gravitational potential in a precise way such that features depending on the theory of gravity governing the system may be distinguished. That is, EMRI signals may be used to test gravity and the geometry of black holes. However, the development of a practical methodology for computing the generation and propagation of GWs from EMRIs in theories of gravity different than General Relativity (GR) has only recently begun. In this paper, we present a parameter estimation study of EMRIs in a particular modification of GR, which is described by a four-dimensional Chern-Simons (CS) gravitational term. We focus on determining to what extent a space-based GW observatory like LISA could distinguish between GR and CS gravity through the detection of GWs from EMRIs.

preprint2012arXiv

Using Swarm Intelligence To Accelerate Pulsar Timing Analysis

We provide brief notes on a particle swarm-optimisation approach to constraining the properties of a stochastic gravitational-wave background in the first International Pulsar Timing Array data-challenge. The technique employs many computational-agents which explore parameter space, remembering their most optimal positions and also sharing this information with all other agents. It is this sharing of information which accelerates the convergence of all agents to the global best-fit location in a very short number of iterations. Error estimates can also be provided by fitting a multivariate Gaussian to the recorded fitness of all visited points.

preprint2012arXiv

Verifying the no-hair property of massive compact objects with intermediate-mass-ratio inspirals in advanced gravitational-wave detectors

The detection of gravitational waves from the inspiral of a neutron star or stellar-mass black hole into an intermediate-mass black hole (IMBH) promises an entirely new look at strong-field gravitational physics. Gravitational waves from these intermediate-mass-ratio inspirals (IMRIs), systems with mass ratios from ~10:1 to ~100:1, may be detectable at rates of up to a few tens per year by Advanced LIGO/Virgo and will encode a signature of the central body's spacetime. Direct observation of the spacetime will allow us to use the "no-hair" theorem of general relativity to determine if the IMBH is a Kerr black hole (or some more exotic object, e.g. a boson star). Using modified post-Newtonian (pN) waveforms, we explore the prospects for constraining the central body's mass-quadrupole moment in the advanced-detector era. We use the Fisher information matrix to estimate the accuracy with which the parameters of the central body can be measured. We find that for favorable mass and spin combinations, the quadrupole moment of a non-Kerr central body can be measured to within a ~15% fractional error or better using 3.5 pN order waveforms; on the other hand, we find the accuracy decreases to ~100% fractional error using 2 pN waveforms, except for a narrow band of values of the best-fit non-Kerr quadrupole moment.

preprint2012arXiv

Weighing The Evidence For A Gravitational-Wave Background In The First International Pulsar Timing Array Data Challenge

We describe an analysis of the First International Pulsar Timing Array Data Challenge. We employ a robust, unbiased Bayesian framework developed by van Haasteren to study the three Open and Closed datasets, testing various models for each dataset and using MultiNest to recover the evidence for the purposes of Bayesian model-selection. The parameter constraints of the favoured model are confirmed using an adaptive MCMC technique. Our results for Closed1 favoured a gravitational-wave background with strain amplitude at f=1 yr-1, A, of (1.1 +/- 0.1) x 10^{-14}, power spectral-index gamma=4.30 +/- 0.15 and no evidence for red-timing noise or single-sources. The evidence for Closed2 favours a gravitational-wave background with A=(6.1 +/- 0.3) x 10^{-14}, gamma=4.34 +/- 0.09 with no red-timing noise or single-sources. Finally, the evidence for Closed3 favours the presence of red-timing noise and a gravitational-wave background, with no single-sources. The properties of the background were A=(5 +/- 1) x 10^{-15} and gamma=4.23 +/- 0.35, while the properties of the red-noise were N_{red}=(12 +/- 4) ns and gamma_{red}=1.5 +/- 0.3. In all cases the redness of the recovered background is consistent with a source-population of inspiraling supermassive black-hole binaries. We also investigate the effect that down-sampling of the datasets has on parameter constraints and run-time. Finally we provide a proof-of-principle study of the ability of the Bayesian framework used in this paper to reconstruct the angular correlation of gravitational-wave background induced timing-residuals, comparing this to the Hellings and Downs curve.

preprint2011arXiv

Approximate Waveforms for Extreme-Mass-Ratio Inspirals in Modified Gravity Spacetimes

Extreme-mass-ratio inspirals, in which a stellar-mass compact object spirals into a supermassive black hole, are prime candidates for detection with space-borne milliHertz gravitational wave detectors, similar to the Laser Interferometer Space Antenna. The gravitational waves generated during such inspirals encode information about the background in which the small object is moving, providing a tracer of the spacetime geometry and a probe of strong-field physics. In this paper, we construct approximate, "analytic-kludge" waveforms for such inspirals with parameterized post-Einsteinian corrections that allow for generic, model-independent deformations of the supermassive black hole background away from the Kerr metric. These approximate waveforms include all of the qualitative features of true waveforms for generic inspirals, including orbital eccentricity and relativistic precession. The deformations of the Kerr metric are modeled using a recently proposed, modified gravity bumpy metric, which parametrically deforms the Kerr spacetime while ensuring that three approximate constants of the motion remain for geodesic orbits: a conserved energy, azimuthal angular momentum and Carter constant. The deformations represent modified gravity effects and have been analytically mapped to several modified gravity black hole solutions in four dimensions. In the analytic kludge waveforms, the conservative motion is modeled by a post-Newtonian expansion of the geodesic equations in the deformed spacetimes, which in turn induce modifications to the radiation-reaction force. These analytic-kludge waveforms serve as a first step toward complete and model-independent tests of General Relativity with extreme mass-ratio inspirals.

preprint2011arXiv

Forced motion near black holes

We present two methods for integrating forced geodesic equations in the Kerr spacetime, which can accommodate arbitrary forces. As a test case, we compute inspirals under a simple drag force, mimicking the presence of gas. We verify that both methods give the same results for this simple force. We find that drag generally causes eccentricity to increase throughout the inspiral. This is a relativistic effect qualitatively opposite to what is seen in gravitational-radiation-driven inspirals, and similar to what is observed in hydrodynamic simulations of gaseous binaries. We provide an analytic explanation by deriving the leading order relativistic correction to the Newtonian dynamics. If observed, an increasing eccentricity would provide clear evidence that the inspiral was occurring in a non-vacuum environment. Our two methods are especially useful for evolving orbits in the adiabatic regime. Both use the method of osculating orbits, in which each point on the orbit is characterized by the parameters of the geodesic with the same instantaneous position and velocity. Both methods describe the orbit in terms of the geodesic energy, axial angular momentum, Carter constant, azimuthal phase, and two angular variables that increase monotonically and are relativistic generalizations of the eccentric anomaly. The two methods differ in their treatment of the orbital phases and the representation of the force. In one method the geodesic phase and phase constant are evolved together as a single orbital phase parameter, and the force is expressed in terms of its components on the Kinnersley orthonormal tetrad. In the second method, the phase constants of the geodesic motion are evolved separately and the force is expressed in terms of its Boyer-Lindquist components. This second approach is a generalization of earlier work by Pound and Poisson for planar forces in a Schwarzschild background.

preprint2011arXiv

Fundamental physics and cosmology with LISA

In this article we give a brief review of the fundamental physics that can be done with the future space-based gravitational wave detector LISA. This includes detection of gravitational wave bursts coming from cosmic strings, measuring a stochastic gravitational wave background, mapping spacetime around massive compact objects in galactic nuclei with extreme-mass-ratio inspirals and testing the predictions of General Relativity for the strong dynamical fields of inspiralling binaries. We give particular attention to new results which show the capability of LISA to constrain cosmological parameters using observations of coalescing massive Black Hole binaries.

preprint2011arXiv

Importance of including small body spin effects in the modelling of extreme and intermediate mass-ratio inspirals

We explore the ability of future low-frequency gravitational wave detectors to measure the spin of stellar mass and intermediate mass black holes that inspiral onto super-massive Kerr black holes (SMBHs). We develop a kludge waveform model based on the equations of motion derived by Saijo et al. [Phys Rev D 58, 064005, 1998] for spinning BH binaries, augmented with spin-orbit and spin-spin couplings taken from perturbative and post-Newtonian (PN) calculations, and the associated conservative self-force corrections, derived by comparison to PN results. We model the inspiral phase using accurate fluxes which include perturbative corrections for the spin of the inspiralling body, spin-spin couplings and higher-order fits to solutions of the Teukolsky equation. We present results of Monte Carlo simulations of parameter estimation errors and of the model errors that arise when we omit conservative corrections from the waveform template. For a source 5000+10^6 solar mass observed with an SNR of 1000, LISA will be able to determine the two masses to within a fractional error of ~0.001, measure the SMBH spin magnitude, q, and the spin magnitude of the inspiralling BH to 0.0001, 10%, respectively, and determine the location of the source in the sky and the SMBH spin orientation to within 0.0001 steradians. For a 10+10^6 solar mass system observed with SNR of 30, LISA will not be able to determine the spin magnitude of the inspiralling BH, although the measurement of the other waveform parameters is not significantly degraded by the presence of spin. The model errors which arise from ignoring conservative corrections become significant for mass-ratios above 0.0001, but including these corrections up to 2PN order may be sufficient to reduce these systematic errors to an acceptable level.

preprint2011arXiv

Influence of conservative corrections on parameter estimation for extreme-mass-ratio inspirals

We present an improved numerical kludge waveform model for circular, equatorial EMRIs. The model is based on true Kerr geodesics, augmented by radiative self-force corrections derived from perturbative calculations, and in this paper for the first time we include conservative self-force corrections that we derive by comparison to post-Newtonian results. We present results of a Monte Carlo simulation of parameter estimation errors computed using the Fisher Matrix and also assess the theoretical errors that would arise form omitting the conservative correction terms we include here. We present results for three different types of system, namely the inspirals of black holes, neutron stars or white dwarfs into a supermassive black hole (SMBH). The analysis shows that for a typical source (a 10 solar mass compact object captured by a one million solar mass SMBH at signal to noise ratio of 30) we expect to determine the two masses to within a fractional error of ~0.0001, measure the spin parameter q to ~0.00003 and determine the location of the source on the sky and the spin orientation to within 0.001 steradians. We show that, for this kludge model, omitting the conservative corrections leads to a small error over much of the parameter space, i.e., the ratio R of the theoretical model error to the Fisher Matrix error is R<1 for all ten parameters in the model. For the few systems with larger errors typically R<3 and hence the conservative corrections can be marginally ignored. In addition, we use our model and first order self-force results for Schwarzschild black holes to estimate the error that arises from omitting the second-order radiative piece of the self-force. This indicates that it may not be necessary to go beyond first order to recover accurate parameter estimates.

preprint2011arXiv

Intermediate-mass-ratio-inspirals in the Einstein Telescope: I. Signal-to-noise ratio calculations

The Einstein Telescope (ET) is a proposed third generation ground-based interferometer, for which the target is a sensitivity that is a factor of ten better than Advanced LIGO and a frequency range that extends down to about 1Hz. ET will provide opportunities to test Einstein's theory of relativity in the strong field and will realize precision gravitational wave astronomy with a thousandfold increase in the expected number of events over the advanced ground-based detectors. A design study for ET is currently underway, so it is timely to assess the science that could be done with such an instrument. This paper is the first in a series that will carry out a detailed study of intermediate-mass-ratio inspirals (IMRIs) for ET. In the context of ET, an IMRI is the inspiral of a neutron star or stellar-mass black hole into an intermediate mass black hole (IMBH). In this paper we focus on the development of IMRI waveform models for circular and equatorial inspirals. We consider two approximations for the waveforms, which both incorporate the inspiral, merger and ringdown phases in a consistent way. One approximation, valid for IMBHs of arbitrary spin, uses the transition model of Ori and Thorne [1] to describe the merger, and this is then matched smoothly onto a ringdown waveform. The second approximation uses the Effective One Body (EOB) approach to model the merger phase of the waveform and is valid for non-spinning IMBHs. In this paper, we use both waveform models to compute signal-to-noise ratios (SNRs) for IMRI sources detectable by ET. At a redshift of z=1, we find typical SNRs for IMRI systems with masses 1.4+100 solar masses, 10+100 solar masses, 1.4+500 solar masses and 10+500 solar masses of about 10-25, 40-80, 3-15 and 10-60, respectively. We also find that the two models make predictions for non-spinning inspirals that are consistent to about ten percent.

preprint2011arXiv

Intermediate-mass-ratio-inspirals in the Einstein Telescope. II. Parameter estimation errors

We explore the precision with which the Einstein Telescope (ET) will be able to measure the parameters of intermediate-mass-ratio inspirals (IMRIs). We calculate the parameter estimation errors using the Fisher Matrix formalism and present results of a Monte Carlo simulation of these errors over choices for the extrinsic parameters of the source. These results are obtained using two different models for the gravitational waveform which were introduced in paper I of this series. These two waveform models include the inspiral, merger and ringdown phases in a consistent way. One of the models, based on the transition scheme of Ori & Thorne [1], is valid for IMBHs of arbitrary spin, whereas the second model, based on the Effective One Body (EOB) approach, has been developed to cross-check our results in the non-spinning limit. In paper I of this series, we demonstrated the excellent agreement in both phase and amplitude between these two models for non-spinning black holes, and that their predictions for signal-to-noise ratios (SNRs) are consistent to within ten percent. We now use these models to estimate parameter estimation errors for binary systems with masses 1.4+100, 10+100, 1.4+500 and 10+500 solar masses (SMs), and various choices for the spin of the central intermediate-mass black hole (IMBH). Assuming a detector network of three ETs, the analysis shows that for a 10 SM compact object (CO) inspiralling into a 100 SM IMBH with spin q=0.3, detected with an SNR of 30, we should be able to determine the CO and IMBH masses, and the IMBH spin magnitude to fractional accuracies of 0.001, 0.0003, and 0.001, respectively. We also expect to determine the location of the source in the sky and the luminosity distance to within 0.003 steradians, and 10%, respectively. We also assess how the precision of parameter determination depends on the network configuration.

preprint2010arXiv

Classifying LISA gravitational wave burst signals using Bayesian evidence

We consider the problem of characterisation of burst sources detected with the Laser Interferometer Space Antenna (LISA) using the multi-modal nested sampling algorithm, MultiNest. We use MultiNest as a tool to search for modelled bursts from cosmic string cusps, and compute the Bayesian evidence associated with the cosmic string model. As an alternative burst model, we consider sine-Gaussian burst signals, and show how the evidence ratio can be used to choose between these two alternatives. We present results from an application of MultiNest to the last round of the Mock LISA Data Challenge, in which we were able to successfully detect and characterise all three of the cosmic string burst sources present in the release data set. We also present results of independent trials and show that MultiNest can detect cosmic string signals with signal-to-noise ratio (SNR) as low as ~7 and sine-Gaussian signals with SNR as low as ~8. In both cases, we show that the threshold at which the sources become detectable coincides with the SNR at which the evidence ratio begins to favour the correct model over the alternative.

preprint2010arXiv

Detecting coalescences of intermediate-mass black holes in globular clusters with the Einstein Telescope

We discuss the capability of a third-generation ground-based detector such as the Einstein Telescope to detect mergers of intermediate-mass black holes that may have formed through runaway stellar collisions in globular clusters. We find that detection rates of 500 events per year are plausible.

preprint2010arXiv

Exploring intermediate and massive black-hole binaries with the Einstein Telescope

We discuss the capability of a third-generation ground-based detector such as the Einstein Telescope (ET) to enhance our astrophysical knowledge through detections of gravitational waves emitted by binaries including intermediate-mass and massive black holes. The design target for such instruments calls for improved sensitivity at low frequencies, specifically in the ~ 1-10 Hz range. This will allow the detection of gravitational waves generated in binary systems containing black holes of intermediate mass, ~ 100-1000 solar masses. We primarily discuss two different source types -- mergers between two intermediate mass black holes (IMBHs) of comparable mass, and intermediate-mass-ratio inspirals (IMRIs) of smaller compact objects with mass ~ 1-10 solar masses into IMBHs. IMBHs may form via two channels: (i) in dark matter halos at high redshift through direct collapse or the collapse of very massive metal-poor Population III stars, or (ii) via runaway stellar collisions in globular clusters. In this paper, we will discuss both formation channels, and both classes of merger in each case. We review existing rate estimates where these exist in the literature, and provide some new calculations for the approximate numbers of events that will be seen by a detector like the Einstein Telescope. These results indicate that the ET may see a few to a few thousand comparable-mass IMBH mergers and as many as several hundred IMRI events per year. These observations will significantly enhance our understanding of galactic black-hole growth, of the existence and properties of IMBHs and of the astrophysics of globular clusters. We finish our review with a discussion of some more speculative sources of gravitational waves for the ET, including hypermassive white dwarfs and eccentric stellar-mass compact-object binaries.

preprint2010arXiv

Gravitational wave energy spectrum of a parabolic encounter

We derive an analytic expression for the energy spectrum of gravitational waves from a parabolic Keplerian binary by taking the limit of the Peters and Matthews spectrum for eccentric orbits. This demonstrates that the location of the peak of the energy spectrum depends primarily on the orbital periapse rather than the eccentricity. We compare this weak-field result to strong-field calculations and find it is reasonably accurate (~10%) provided that the azimuthal and radial orbital frequencies do not differ by more than ~10%. For equatorial orbits in the Kerr spacetime, this corresponds to periapse radii of rp > 20M. These results can be used to model radiation bursts from compact objects on highly eccentric orbits about massive black holes in the local Universe, which could be detected by LISA.

preprint2010arXiv

LISA extreme-mass-ratio inspiral events as probes of the black hole mass function

One of the sources of gravitational waves for the proposed space-based gravitational wave detector, the Laser Interferometer Space Antenna (LISA), are the inspirals of compact objects into supermassive black holes in the centres of galaxies - extreme-mass-ratio inspirals (EMRIs). Using LISA observations, we will be able to measure the parameters of each EMRI system detected to very high precision. However, the statistics of the set of EMRI events observed by LISA will be more important in constraining astrophysical models than extremely precise measurements for individual systems. The black holes to which LISA is most sensitive are in a mass range that is difficult to probe using other techniques, so LISA provides an almost unique window onto these objects. In this paper we explore, using Bayesian techniques, the constraints that LISA EMRI observations can place on the mass function of black holes at low redshift. We describe a general framework for approaching inference of this type --- using multiple observations in combination to constrain a parameterised source population. Assuming that the scaling of EMRI rate with black hole mass is known and taking a black hole distribution given by a simple power law, dn/d(ln M) = A (M/M_*)^b, we find that LISA could measure the parameters to a precision of D(ln A) ~ 0.08, and D(b) ~ 0.03 for a reference model that predicts ~1000 events. Even with as few as 10 events, LISA should constrain the slope to a precision ~0.3, which is the current level of observational uncertainty in the low-mass slope of the black hole mass function. We also consider a model in which A and b evolve with redshift, but find that EMRI observations alone do not have much power to probe such an evolution.

preprint2010arXiv

Reconstructing the massive black hole cosmic history through gravitational waves

The massive black holes we observe in galaxies today are the natural end-product of a complex evolutionary path, in which black holes seeded in proto-galaxies at high redshift grow through cosmic history via a sequence of mergers and accretion episodes. Electromagnetic observations probe a small subset of the population of massive black holes (namely, those that are active or those that are very close to us), but planned space-based gravitational-wave observatories such as the Laser Interferometer Space Antenna (LISA) can measure the parameters of ``electromagnetically invisible'' massive black holes out to high redshift. In this paper we introduce a Bayesian framework to analyze the information that can be gathered from a set of such measurements. Our goal is to connect a set of massive black hole binary merger observations to the underlying model of massive black hole formation. In other words, given a set of observed massive black hole coalescences, we assess what information can be extracted about the underlying massive black hole population model. For concreteness we consider ten specific models of massive black hole formation, chosen to probe four important (and largely unconstrained) aspects of the input physics used in structure formation simulations: seed formation, metallicity ``feedback'', accretion efficiency and accretion geometry. For the first time we allow for the possibility of ``model mixing'', by drawing the observed population from some combination of the ``pure'' models that have been simulated. A Bayesian analysis allows us to recover a posterior probability distribution for the ``mixing parameters'' that characterize the fractions of each model represented in the observed distribution. Our work shows that LISA has enormous potential to probe the underlying physics of structure formation.

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.

preprint2009arXiv

Cosmic Swarms: A search for Supermassive Black Holes in the LISA data stream with a Hybrid Evolutionary Algorithm

We describe a hybrid evolutionary algorithm that can simultaneously search for multiple supermassive black hole binary (SMBHB) inspirals in LISA data. The algorithm mixes evolutionary computation, Metropolis-Hastings methods and Nested Sampling. The inspiral of SMBHBs presents an interesting problem for gravitational wave data analysis since, due to the LISA response function, the sources have a bi-modal sky solution. We show here that it is possible not only to detect multiple SMBHBs in the data stream, but also to investigate simultaneously all the various modes of the global solution. In all cases, the algorithm returns parameter determinations within $5σ$ (as estimated from the Fisher Matrix) of the true answer, for both the actual and antipodal sky solutions.

preprint2009arXiv

Use of the MultiNest algorithm for gravitational wave data analysis

We describe an application of the MultiNest algorithm to gravitational wave data analysis. MultiNest is a multimodal nested sampling algorithm designed to efficiently evaluate the Bayesian evidence and return posterior probability densities for likelihood surfaces containing multiple secondary modes. The algorithm employs a set of live points which are updated by partitioning the set into multiple overlapping ellipsoids and sampling uniformly from within them. This set of live points climbs up the likelihood surface through nested iso-likelihood contours and the evidence and posterior distributions can be recovered from the point set evolution. The algorithm is model-independent in the sense that the specific problem being tackled enters only through the likelihood computation, and does not change how the live point set is updated. In this paper, we consider the use of the algorithm for gravitational wave data analysis by searching a simulated LISA data set containing two non-spinning supermassive black hole binary signals. The algorithm is able to rapidly identify all the modes of the solution and recover the true parameters of the sources to high precision.

preprint2008arXiv

"Kludge" gravitational waveforms for a test-body orbiting a Kerr black hole

One of the most exciting potential sources of gravitational waves for low-frequency, space-based gravitational wave (GW) detectors such as the proposed Laser Interferometer Space Antenna (LISA) is the inspiral of compact objects into massive black holes in the centers of galaxies. The detection of waves from such "extreme mass ratio inspiral" systems (EMRIs) and extraction of information from those waves require template waveforms. The systems' extreme mass ratio means that their waveforms can be determined accurately using black hole perturbation theory. Such calculations are computationally very expensive. There is a pressing need for families of approximate waveforms that may be generated cheaply and quickly but which still capture the main features of true waveforms. In this paper, we introduce a family of such "kludge" waveforms and describe ways to generate them. We assess performance of the introduced approximations by comparing "kludge" waveforms to accurate waveforms obtained by solving the Teukolsky equation in the adiabatic limit (neglecting GW backreaction). We find that the kludge waveforms do extremely well at approximating the true gravitational waveform, having overlaps with the Teukolsky waveforms of 95% or higher over most of the parameter space for which comparisons can currently be made. Indeed, we find these kludges to be of such high quality (despite their ease of calculation) that it is possible they may play some role in the final search of LISA data for EMRIs.

preprint2008arXiv

Rates and Characteristics of Intermediate Mass Ratio Inspirals Detectable by Advanced LIGO

Gravitational waves (GWs) from the inspiral of a neutron star (NS) or stellar-mass black hole (BH) into an intermediate-mass black hole (IMBH) with mass between ~50 and ~350 solar masses may be detectable by the planned advanced generation of ground-based GW interferometers. Such intermediate mass ratio inspirals (IMRIs) are most likely to be found in globular clusters. We analyze four possible IMRI formation mechanisms: (1) hardening of an NS-IMBH or BH-IMBH binary via three-body interactions, (2) hardening via Kozai resonance in a hierarchical triple system, (3) direct capture, and (4) inspiral of a compact object from a tidally captured main-sequence star; we also discuss tidal effects when the inspiraling object is an NS. For each mechanism we predict the typical eccentricities of the resulting IMRIs. We find that IMRIs will have largely circularized by the time they enter the sensitivity band of ground-based detectors. Hardening of a binary via three-body interactions, which is likely to be the dominant mechanism for IMRI formation, yields eccentricities under 10^-4 when the GW frequency reaches 10 Hz. Even among IMRIs formed via direct captures, which can have the highest eccentricities, around 90% will circularize to eccentricities under 0.1 before the GW frequency reaches 10 Hz. We estimate the rate of IMRI coalescences in globular clusters and the sensitivity of a network of three Advanced LIGO detectors to the resulting GWs. We show that this detector network may see up to tens of IMRIs per year, although rates of one to a few per year may be more plausible. We also estimate the loss in signal-to-noise ratio that will result from using circular IMRI templates for data analysis and find that, for the eccentricities we expect, this loss is negligible.

Jonathan R. Gair

What is connected

Connect this record

See the researcher in context

Building this map preview

56 published item(s)

Trust the process: mapping data-driven reconstructions to informed models using stochastic processes

2 Fast 2 Fiducial: Gaussian processes for the interpolation and marginalization of waveform error in extreme-mass-ratio-inspiral parameter estimation

Constraining the spin parameter of near-extremal black holes using LISA

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

Transition from Inspiral to Plunge: A Complete Near-Extremal Trajectory and Associated Waveform

Fast methods for training Gaussian processes on large data sets

Improving gravitational-wave parameter estimation using Gaussian process regression

Tunable compression of template banks for fast gravitational-wave detection and localisation

Detection of eccentric supermassive black hole binaries with pulsar timing arrays: Signal-to-noise ratio calculations

Estimating the sensitivity of pulsar timing arrays

European Pulsar Timing Array Limits On An Isotropic Stochastic Gravitational-Wave Background

European Pulsar Timing Array Limits on Continuous Gravitational Waves from Individual Supermassive Black Hole Binaries

Improved analytic extreme-mass-ratio inspiral model for scoping out eLISA data analysis

Likelihood smoothing using gravitational wave surrogate models

Mapping gravitational-wave backgrounds in modified theories of gravity using pulsar timing arrays

Quantifying and mitigating bias in inference on gravitational wave source populations

Testing the "no-hair" property of black holes with X-ray observations of accretion disks

Counting And Confusion: Bayesian Rate Estimation With Multiple Populations

Electromagnetic signatures of far-field gravitational radiation in the 1+3 approach

Mapping gravitational-wave backgrounds using methods from CMB analysis: Application to pulsar timing arrays

Novel Method for Incorporating Model Uncertainties into Gravitational Wave Parameter Estimates

Research Update on Extreme-Mass-Ratio Inspirals

Self-forced evolutions of an implicit rotating source: A natural framework to model comparable and intermediate mass-ratio systems from inspiral through ringdown

The scientific potential of space-based gravitational wave detectors

Gravitational wave parameter estimation with compressed likelihood evaluations

Searching For Anisotropic Gravitational-wave Backgrounds Using Pulsar Timing Arrays

Testing General Relativity with Low-Frequency, Space-Based Gravitational-Wave Detectors

Constraining Gravity with LISA Detections of Binaries

Cosmology with the lights off: Standard sirens in the Einstein Telescope era

Extreme-mass-ratio bursts from the Galactic Centre

Hubble without the Hubble: cosmology using advanced gravitational-wave detectors alone

Importance of including small body spin effects in the modelling of intermediate mass-ratio inspirals. II Accurate parameter extraction of strong sources using higher-order spin effects

Linearized f(R) Gravity: Gravitational Radiation & Solar System Tests

Testing Chern-Simons Modified Gravity with Gravitational-Wave Detections of Extreme-Mass-Ratio Binaries

Testing Chern-Simons modified gravity with observations of extreme-mass-ratio binaries

Using Swarm Intelligence To Accelerate Pulsar Timing Analysis

Verifying the no-hair property of massive compact objects with intermediate-mass-ratio inspirals in advanced gravitational-wave detectors

Weighing The Evidence For A Gravitational-Wave Background In The First International Pulsar Timing Array Data Challenge

Approximate Waveforms for Extreme-Mass-Ratio Inspirals in Modified Gravity Spacetimes

Forced motion near black holes

Fundamental physics and cosmology with LISA

Importance of including small body spin effects in the modelling of extreme and intermediate mass-ratio inspirals

Influence of conservative corrections on parameter estimation for extreme-mass-ratio inspirals

Intermediate-mass-ratio-inspirals in the Einstein Telescope: I. Signal-to-noise ratio calculations

Intermediate-mass-ratio-inspirals in the Einstein Telescope. II. Parameter estimation errors

Classifying LISA gravitational wave burst signals using Bayesian evidence

Detecting coalescences of intermediate-mass black holes in globular clusters with the Einstein Telescope

Exploring intermediate and massive black-hole binaries with the Einstein Telescope

Gravitational wave energy spectrum of a parabolic encounter

LISA extreme-mass-ratio inspiral events as probes of the black hole mass function

Reconstructing the massive black hole cosmic history through gravitational waves

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

Cosmic Swarms: A search for Supermassive Black Holes in the LISA data stream with a Hybrid Evolutionary Algorithm

Use of the MultiNest algorithm for gravitational wave data analysis

"Kludge" gravitational waveforms for a test-body orbiting a Kerr black hole

Rates and Characteristics of Intermediate Mass Ratio Inspirals Detectable by Advanced LIGO