Source author record

Jonathan Gair

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

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

22 published item(s)

preprint2022arXiv

Constraining Unmodeled Physics with Compact Binary Mergers from GWTC-1

We present a flexible model to describe the effects of generic deviations of observed gravitational wave signals from modeled waveforms in the LIGO and Virgo gravitational wave detectors. With the detection of 11 gravitational wave events from the GWTC-1 catalog, we are able to constrain possible deviations from our modeled waveforms. In this paper we present our coherent spline model that describes the deviations, then choose to validate our model on an example phenomenological and astrophysically motivated departure in waveforms based on extreme spontaneous scalarization. We find that the model is capable of recovering the simulated deviations. By performing model comparisons we observe that the spline model effectively describes the simulated departures better than a normal compact binary coalescence (CBC) model. We analyze the entire GWTC-1 catalog of events with our model and compare it to a normal CBC model, finding that there are no significant departures from the modeled template gravitational waveforms used.

preprint2022arXiv

Known unknowns: assessing the impact of instrumental calibration uncertainty on LISA science

The primary scientific results of the future space-based gravitational wave interferometer LISA will come from the parameter inference of a large variety of gravitational wave sources. However, the presence of calibration errors could potentially degrade the measurement precision of the system parameters. Here, we assess the impact of calibration uncertainties on parameter estimation for individual sources, focusing on massive black holes, extreme-mass-ratio inspirals (EMRIs), galactic binaries, and stellar origin black hole binaries. Using a Fisher matrix formalism, we investigate how the measurement precision of source parameters degrades as a function of the size of the assumed calibration uncertainties. If we require that parameter measurements are degraded by no more than a factor of two relative to their value in the absence of calibration error, we find that calibration errors should be smaller than a few tenths of a percent in amplitude and $10^{-3}$ in phase. We also investigate the possibility of using verification binaries and EMRIs to constrain calibration uncertainties. Verification binaries can constrain amplitude calibration uncertainties at the level of a few percent, while both source types can provide constrain phase calibration at the level of a few$\times10^{-2}$.

preprint2022arXiv

New Horizons for Fundamental Physics with LISA

The Laser Interferometer Space Antenna (LISA) has the potential to reveal wonders about the fundamental theory of nature at play in the extreme gravity regime, where the gravitational interaction is both strong and dynamical. In this white paper, the Fundamental Physics Working Group of the LISA Consortium summarizes the current topics in fundamental physics where LISA observations of GWs can be expected to provide key input. We provide the briefest of reviews to then delineate avenues for future research directions and to discuss connections between this working group, other working groups and the consortium work package teams. These connections must be developed for LISA to live up to its science potential in these areas.

preprint2021arXiv

Mapping the inhomogeneous Universe with Standard Sirens: Degeneracy between inhomogeneity and modified gravity theories

The detection of gravitational waves (GWs) and an accompanying electromagnetic (E/M) counterpart have been suggested as a future probe for cosmology and theories of gravity. In this paper, we present calculations of the luminosity distance of sources taking into account inhomogeneities in the matter distribution that are predicted in numerical simulations of structure formation. In addition, we show that inhomogeneities resulting from clustering of matter can mimic certain classes of modified gravity theories, or other effects that dampen GW amplitudes, and deviations larger than $δν\sim \mathcal{O}(0.1)\ (99\%\ \rm{C.L.})$ to the extra friction term $ν$, from zero, would be necessary to distinguish them. For these, we assume mock GWs sources, with known redshift, based on binary population synthesis models, between redshifts $z=0$ and $z=5$. We show that future GW detectors, like Einstein Telescope or Cosmic Explorer, will be needed for strong constraints on the inhomogeneity parameters and breaking the degeneracy between modified gravity effects and matter anisotropies by measuring $ν$ at $5 \%$ and $1 \%$ level with $100$ and $350$ events respectively.

preprint2021arXiv

Observing intermediate-mass black holes and the upper--stellar-mass gap with LIGO and Virgo

Using ground-based gravitational-wave detectors, we probe the mass function of intermediate-mass black holes (IMBHs) wherein we also include BHs in the upper mass gap $\sim 60-130~M_\odot$. Employing the projected sensitivity of the upcoming LIGO and Virgo fourth observing (O4) run, we perform Bayesian analysis on quasi-circular non-precessing, spinning IMBH binaries (IMBHBs) with total masses $50\mbox{--} 500\, M_\odot$, mass ratios 1.25, 4, and 10, and dimensionless spins up to 0.95, and estimate the precision with which the source-frame parameters can be measured. We find that, at $2σ$, the mass of the heavier component of IMBHBs can be constrained with an uncertainty of $\sim 10-40\%$ at a signal-to-noise ratio of $20$. Focusing on the stellar-mass gap with new tabulations of the $^{12}\text{C}(α, γ)^{16} \text{O}$ reaction rate and its uncertanties, we evolve massive helium core stars using \MESA\, to establish the lower and upper edge of the mass gap as $\simeq$\,59$^{+34}_{-13}$\,$M_{\odot}$ and $\simeq$\,139$^{+30}_{-14}$\,$M_{\odot}$ respectively, where the error bars give the mass range that follows from the $\pm 3σ$ uncertainty in the $^{12}\text{C}(α, γ) ^{16} \text{O}$ nuclear reaction rate. We find that high resolution of the tabulated reaction rate and fine temporal resolution are necessary to resolve the peak of the BH mass spectrum. We then study IMBHBs with components lying in the mass gap and show that the O4 run will be able to robustly identify most such systems. Finally, we re-analyse GW190521 with a state-of-the-art aligned-spin waveform model, finding that the primary mass lies in the mass gap with 90\% credibility.

preprint2020arXiv

A Precessing Numerical Relativity Waveform Surrogate Model for Binary Black Holes: A Gaussian Process Regression Approach

Gravitational wave astrophysics relies heavily on the use of matched filtering both to detect signals in noisy data from detectors, and to perform parameter estimation on those signals. Matched filtering relies upon prior knowledge of the signals expected to be produced by a range of astrophysical systems, such as binary black holes. These waveform signals can be computed using numerical relativity techniques, where the Einstein field equations are solved numerically, and the signal is extracted from the simulation. Numerical relativity simulations are, however, computationally expensive, leading to the need for a surrogate model which can predict waveform signals in regions of the physical parameter space which have not been probed directly by simulation. We present a method for producing such a surrogate using Gaussian process regression which is trained directly on waveforms generated by numerical relativity. This model returns not just a single interpolated value for the waveform at a new point, but a full posterior probability distribution on the predicted value. This model is therefore an ideal component in a Bayesian analysis framework, through which the uncertainty in the interpolation can be taken into account when performing parameter estimation of signals.

preprint2020arXiv

Complete parameter inference for GW150914 using deep learning

The LIGO and Virgo gravitational-wave observatories have detected many exciting events over the past five years. As the rate of detections grows with detector sensitivity, this poses a growing computational challenge for data analysis. With this in mind, in this work we apply deep learning techniques to perform fast likelihood-free Bayesian inference for gravitational waves. We train a neural-network conditional density estimator to model posterior probability distributions over the full 15-dimensional space of binary black hole system parameters, given detector strain data from multiple detectors. We use the method of normalizing flows---specifically, a neural spline normalizing flow---which allows for rapid sampling and density estimation. Training the network is likelihood-free, requiring samples from the data generative process, but no likelihood evaluations. Through training, the network learns a global set of posteriors: it can generate thousands of independent posterior samples per second for any strain data consistent with the prior and detector noise characteristics used for training. By training with the detector noise power spectral density estimated at the time of GW150914, and conditioning on the event strain data, we use the neural network to generate accurate posterior samples consistent with analyses using conventional sampling techniques.

preprint2019arXiv

Astrometric Effects of Gravitational Wave Backgrounds with non-Luminal Propagation Speeds

A passing gravitational wave causes a deflection in the apparent astrometric positions of distant stars. The effect of the speed of the gravitational wave on this astrometric shift is discussed. A stochastic background of gravitational waves would result in a pattern of astrometric deflections which are correlated on large angular scales. These correlations are quantified and investigated for backgrounds of gravitational waves with sub- and super-luminal group velocities. The statistical properties of the correlations are depicted in two equivalent and related ways: as correlation curves and as angular power spectra. Sub-(super-)luminal gravitational wave backgrounds have the effect of enhancing (suppressing) the power in low-order angular modes. Analytical representations of the redshift-redshift and redshift-astrometry correlations are also derived. The potential for using this effect for constraining the speed of gravity is discussed.

preprint2016arXiv

Science with the space-based interferometer eLISA. I: Supermassive black hole binaries

We compare the science capabilities of different eLISA mission designs, including four-link (two-arm) and six-link (three-arm) configurations with different arm lengths, low-frequency noise sensitivities and mission durations. For each of these configurations we consider a few representative massive black hole formation scenarios. These scenarios are chosen to explore two physical mechanisms that greatly affect eLISA rates, namely (i) black hole seeding, and (ii) the delays between the merger of two galaxies and the merger of the black holes hosted by those galaxies. We assess the eLISA parameter estimation accuracy using a Fisher matrix analysis with spin-precessing, inspiral-only waveforms. We quantify the information present in the merger and ringdown by rescaling the inspiral-only Fisher matrix estimates using the signal-to-noise ratio from non-precessing inspiral-merger-ringdown phenomenological waveforms, and from a reduced set of precessing numerical relativity/post-Newtonian hybrid waveforms. We find that all of the eLISA configurations considered in our study should detect some massive black hole binaries. However, configurations with six links and better low-frequency noise will provide much more information on the origin of black holes at high redshifts and on their accretion history, and they may allow the identification of electromagnetic counterparts to massive black hole mergers.

preprint2015arXiv

Accelerated gravitational-wave parameter estimation with reduced order modeling

Inferring the astrophysical parameters of coalescing compact binaries is a key science goal of the upcoming advanced LIGO-Virgo gravitational-wave detector network and, more generally, gravitational-wave astronomy. However, current parameter estimation approaches for such scenarios can lead to computationally intractable problems in practice. Therefore there is a pressing need for new, fast and accurate Bayesian inference techniques. In this letter we demonstrate that a reduced order modeling approach enables rapid parameter estimation studies. By implementing a reduced order quadrature scheme within the LIGO Algorithm Library, we show that Bayesian inference on the 9-dimensional parameter space of non-spinning binary neutron star inspirals can be sped up by a factor of 30 for the early advanced detectors' configurations. This speed-up will increase to about $150$ as the detectors improve their low-frequency limit to 10Hz, reducing to hours analyses which would otherwise take months to complete. Although these results focus on gravitational detectors, the techniques are broadly applicable to any experiment where fast Bayesian analysis is desirable.

preprint2015arXiv

Expected properties of the first gravitational wave signal detected with pulsar timing arrays

In this paper we attempt to investigate the nature of the first gravitational wave (GW) signal to be detected by pulsar timing arrays (PTAs): will it be an individual, resolved supermassive black hole binary (SBHB), or a stochastic background made by the superposition of GWs produced by an ensemble of SBHBs? To address this issue, we analyse a broad set of simulations of the cosmological population of SBHBs, that cover the entire parameter space allowed by current electromagnetic observations in an unbiased way. For each simulation, we construct the expected GW signal and identify the loudest individual sources. We then employ appropriate detection statistics to evaluate the relative probability of detecting each type of source as a function of time for a variety of PTAs; we consider the current International PTA, and speculate into the era of the Square Kilometre Array. The main properties of the first detectable individual SBHBs are also investigated. Contrary to previous work, we cast our results in terms of the detection probability (DP), since the commonly adopted criterion based on a signal-to-noise ratio threshold is statistic-dependent and may result in misleading conclusions for the statistics adopted here. Our results confirm quantitatively that a stochastic signal is more likely to be detected first (with between 75 to 93 per cent probability, depending on the array), but the DP of single-sources is not negligible. Our framework is very flexible and can be easily extended to more realistic arrays and to signal models including environmental coupling and SBHB eccentricity.

preprint2015arXiv

Phase-coherent mapping of gravitational-wave backgrounds using ground-based laser interferometers

We extend the formalisms developed in Gair et al. and Cornish and van Haasteren to create maps of gravitational-wave backgrounds using a network of ground-based laser interferometers. We show that in contrast to pulsar timing arrays, which are insensitive to half of the gravitational-wave sky (the curl modes), a network of ground-based interferometers is sensitive to both the gradient and curl components of the background. The spatial separation of a network of interferometers, or of a single interferometer at different times during its rotational and orbital motion around the Sun, allows for recovery of both components. We derive expressions for the response functions of a laser interferometer in the small-antenna limit, and use these expressions to calculate the overlap reduction function for a pair of interferometers. We also construct maximum-likelihood estimates of the + and x-polarization modes of the gravitational-wave sky in terms of the response matrix for a network of ground-based interferometers, evaluated at discrete times during Earth's rotational and orbital motion around the Sun. We demonstrate the feasibility of this approach for some simple simulated backgrounds (a single point source and spatially-extended distributions having only grad or curl components), calculating maximum-likelihood sky maps and uncertainty maps based on the (pseudo)inverse of the response matrix. The distinction between this approach and standard methods for mapping gravitational-wave power is also discussed.

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.

preprint2014arXiv

Accelerated Bayesian model-selection and parameter-estimation in continuous gravitational-wave searches with pulsar-timing arrays

We describe several new techniques which accelerate Bayesian searches for continuous gravitational-wave emission from supermassive black-hole binaries using pulsar timing arrays. These techniques mitigate the problematic increase of search-dimensionality with the size of the pulsar array which arises from having to include an extra parameter per pulsar as the array is expanded. This extra parameter corresponds to searching over the phase of the gravitational-wave as it propagates past each pulsar so that we can coherently include the pulsar-term in our search strategies. Our techniques make the analysis tractable with powerful evidence-evaluation packages like MultiNest. We find good agreement of our techniques with the parameter-estimation and Bayes factor evaluation performed with full signal templates, and conclude that these techniques make excellent first-cut tools for detection and characterisation of continuous gravitational-wave signals with pulsar timing arrays. Crucially, at low to moderate signal-to-noise ratios the factor by which the analysis is sped up can be > 100, permitting rigorous programs of systematic injection and recovery of signals to establish robust detection criteria within a Bayesian formalism.

preprint2014arXiv

Method for estimation of gravitational-wave transient model parameters in frequency-time maps

A common technique for detection of gravitational-wave signals is searching for excess power in frequency-time maps of gravitational-wave detector data. In the event of a detection, model selection and parameter estimation will be performed in order to explore the properties of the source. In this paper, we develop a Bayesian statistical method for extracting model-dependent parameters from observed gravitational-wave signals in frequency-time maps. We demonstrate the method by recovering the parameters of model gravitational-wave signals added to simulated advanced LIGO noise. We also characterize the performance of the method and discuss prospects for future work.

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.

preprint2012arXiv

eLISA: Astrophysics and cosmology in the millihertz regime

This document introduces the exciting and fundamentally new science and astronomy that the European New Gravitational Wave Observatory (NGO) mission (derived from the previous LISA proposal) will deliver. The mission (which we will refer to by its informal name "eLISA") will survey for the first time the low-frequency gravitational wave band (about 0.1 mHz to 1 Hz), with sufficient sensitivity to detect interesting individual astrophysical sources out to z = 15. The eLISA mission will discover and study a variety of cosmic events and systems with high sensitivity: coalescences of massive black holes binaries, brought together by galaxy mergers; mergers of earlier, less-massive black holes during the epoch of hierarchical galaxy and black-hole growth; stellar-mass black holes and compact stars in orbits just skimming the horizons of massive black holes in galactic nuclei of the present era; extremely compact white dwarf binaries in our Galaxy, a rich source of information about binary evolution and about future Type Ia supernovae; and possibly most interesting of all, the uncertain and unpredicted sources, for example relics of inflation and of the symmetry-breaking epoch directly after the Big Bang. eLISA's measurements will allow detailed studies of these signals with high signal-to-noise ratio, addressing most of the key scientific questions raised by ESA's Cosmic Vision programme in the areas of astrophysics and cosmology. They will also provide stringent tests of general relativity in the strong-field dynamical regime, which cannot be probed in any other way. This document not only describes the science but also gives an overview on the mission design and orbits.

preprint2012arXiv

Low-frequency gravitational-wave science with eLISA/NGO

We review the expected science performance of the New Gravitational-Wave Observatory (NGO, a.k.a. eLISA), a mission under study by the European Space Agency for launch in the early 2020s. eLISA will survey the low-frequency gravitational-wave sky (from 0.1 mHz to 1 Hz), detecting and characterizing a broad variety of systems and events throughout the Universe, including the coalescences of massive black holes brought together by galaxy mergers; the inspirals of stellar-mass black holes and compact stars into central galactic black holes; several millions of ultracompact binaries, both detached and mass transferring, in the Galaxy; and possibly unforeseen sources such as the relic gravitational-wave radiation from the early Universe. eLISA's high signal-to-noise measurements will provide new insight into the structure and history of the Universe, and they will test general relativity in its strong-field dynamical regime.

preprint2011arXiv

Graviton mass bounds from space-based gravitational-wave observations of massive black hole populations

Space-based gravitational-wave detectors, such as LISA or a similar ESA-led mission, will offer unique opportunities to test general relativity. We study the bounds that space-based detectors could place on the graviton Compton wavelength λ_g=h/(m_g c) by observing multiple inspiralling black hole binaries. We show that while observations of individual inspirals will yield mean bounds λ_g~3x10^15 km, the combined bound from observing ~50 events in a two-year mission is about ten times better: λ_g~3x10^16 km (m_g~4x10^-26 eV). The bound improves faster than the square root of the number of observed events, because typically a few sources provide constraints as much as three times better than the mean. This result is only mildly dependent on details of black hole formation and detector characteristics. The bound achievable in practice should be one order of magnitude better than this figure (and hence almost competitive with the static, model-dependent bounds from gravitational effects on cosmological scales), because our calculations ignore the merger/ringdown portion of the waveform. The observation that an ensemble of events can sensibly improve the bounds that individual binaries set on λ_g applies to any theory whose deviations from general relativity are parametrized by a set of global parameters.

preprint2009arXiv

Massive Black Hole Binary Inspirals: Results from the LISA Parameter Estimation Taskforce

The LISA Parameter Estimation (LISAPE) Taskforce was formed in September 2007 to provide the LISA Project with vetted codes, source distribution models, and results related to parameter estimation. The Taskforce's goal is to be able to quickly calculate the impact of any mission design changes on LISA's science capabilities, based on reasonable estimates of the distribution of astrophysical sources in the universe. This paper describes our Taskforce's work on massive black-hole binaries (MBHBs). Given present uncertainties in the formation history of MBHBs, we adopt four different population models, based on (i) whether the initial black-hole seeds are small or large, and (ii) whether accretion is efficient or inefficient at spinning up the holes. We compare four largely independent codes for calculating LISA's parameter-estimation capabilities. All codes are based on the Fisher-matrix approximation, but in the past they used somewhat different signal models, source parametrizations and noise curves. We show that once these differences are removed, the four codes give results in extremely close agreement with each other. Using a code that includes both spin precession and higher harmonics in the gravitational-wave signal, we carry out Monte Carlo simulations and determine the number of events that can be detected and accurately localized in our four population models.

preprint2009arXiv

Observing gravitational waves from the first generation of black holes

The properties of the first generation of black-hole seeds trace and distinguish different models of formation of cosmic structure in the high-redshift universe. The observational challenge lies in identifying black holes in the mass range ~100-1000 solar masses at redshift z~10. The typical frequencies of gravitational waves produced by the coalescence of the first generation of light seed black-hole binaries fall in the gap between the spectral ranges of low-frequency space-borne detectors (e.g., LISA) and high-frequency ground-based detectors (e.g., LIGO, Virgo and GEO 600). As such, these sources are targets for proposed third-generation ground-based instruments, such as the Einstein Telescope which is currently in design study. Using galaxy merger trees and four different models of black hole accretion - which are meant to illustrate the potential of this new type of source rather than to yield precise event-rate predictions - we find that such detectors could observe a few to a few tens of seed black-hole merger events in three years and provide, possibly unique, information on the evolution of structure in the corresponding era. We show further that a network of detectors may be able to measure the luminosity distance to sources to a precision of ~40%, allowing us to be confident of the high-redshift nature of the sources.

preprint2009arXiv

The massive Dirac field on a rotating black hole spacetime: Angular solutions

The massive Dirac equation on a Kerr-Newman background may be solved by the method of separation of variables. The radial and angular equations are coupled via an angular eigenvalue, which is determined from the Chandrasekhar-Page (CP) equation. Obtaining accurate angular eigenvalues is a key step in studying scattering, absorption and emission of the fermionic field. Here we introduce a new method for finding solutions of the CP equation. First, we introduce a novel representation for the spin-half spherical harmonics. Next, we decompose the angular solutions of the CP equation (the mass-dependent spin-half spheroidal harmonics) in the spherical basis. The method yields a three-term recurrence relation which may be solved numerically via continued-fraction methods, or perturbatively to obtain a series expansion for the eigenvalues. In the case $μ= \pm ω$ (where $ω$ and $μ$ are the frequency and mass of the fermion) we obtain eigenvalues and eigenfunctions in closed form. We study the eigenvalue spectrum, and the zeros of the maximally co-rotating mode. We compare our results with previous studies, and uncover and correct some errors in the literature. We provide series expansions, tables of eigenvalues and numerical fits across a wide parameter range, and present plots of a selection of eigenfunctions. It is hoped this study will be a useful resource for all researchers interested in the Dirac equation on a rotating black hole background.