Extreme-mass-ratio-inspiral observations from future space-based gravitational-wave detectors such as LISA will enable strong-field tests of general relativity with unprecedented precision, but at prohibitive computational cost if existing statistical techniques are used. In one such test that is currently employed for LIGO black-hole binary mergers, generic deviations from relativity are represented by $N$ deformation parameters in a generalised waveform model; the Bayesian evidence for each of its $2^N$ combinatorial submodels is then combined into a posterior odds ratio for modified gravity over relativity in a null-hypothesis test. We adapt and apply this test to a generalised model for extreme-mass-ratio inspirals constructed on deformed black-hole spacetimes, and focus our investigation on how computational efficiency can be increased through an evidence-free method of model selection. This method is akin to the algorithm known as product-space Markov chain Monte Carlo, but uses nested sampling and improved error estimates from a rethreading technique. We perform benchmarking and robustness checks for the method, and find order-of-magnitude computational gains over regular nested sampling in the case of synthetic data generated from the null model.

The Gaia mission offers a new opportunity to search for the low frequency gravitational wave background using astrometric measurements. In this paper, the astrometric effect of gravitational waves is reviewed, with a particular focus on the effect of non-Einsteinian gravitational wave polarizations. A stochastic gravitational wave background generates a correlated vector field of astrometric deflections on the sky. A convenient decomposition for the correlation matrix is introduced, enabling it to be calculated for all possible gravitational wave polarizations and compared to the redshift correlations from the pulsar-timing literature; in the case of a GR background of transverse traceless gravitational waves, this also allows us to identify an astrometric analog of the famous Hellings-Downs curve. Finally, the cross-correlation between the redshift and astrometric signal is also calculated; this may form the basis for future joint pulsar-timing and astrometry searches for arbitrarily polarized gravitational wave backgrounds.

The observation of the electromagnetic counterpart of gravitational-wave (GW) transient GW170817 demonstrated the potential in extracting astrophysical information from multimessenger discoveries. The forthcoming deployment of the first telescopes of the Cherenkov Telescope Array (CTA) observatory will coincide with Advanced LIGO/Virgo's next observing run, O3, enabling the monitoring of gamma-ray emission at E > 20 GeV, and thus particle acceleration, from GW sources. CTA will not be greatly limited by the precision of GW localization as it will be be capable of rapidly covering the GW error region with sufficient sensitivity. We examine the current status of GW searches and their follow-up effort, as well as the status of CTA, in order to identify some of the general strategies that will enhance CTA's contribution to multimessenger discoveries.

We present $\texttt{ENIGMA}$, a time domain, inspiral-merger-ringdown waveform model that describes non-spinning binary black holes systems that evolve on moderately eccentric orbits. The inspiral evolution is described using a consistent combination of post-Newtonian theory, self-force and black hole perturbation theory. Assuming eccentric binaries that circularize prior to coalescence, we smoothly match the eccentric inspiral with a stand-alone, quasi-circular merger, which is constructed using machine learning algorithms that are trained with quasi-circular numerical relativity waveforms. We show that $\texttt{ENIGMA}$ reproduces with excellent accuracy the dynamics of quasi-circular compact binaries. We validate $\texttt{ENIGMA}$ using a set of $\texttt{Einstein Toolkit}$ eccentric numerical relativity waveforms, which describe eccentric binary black hole mergers with mass-ratios between $1 \leq q \leq 5.5$, and eccentricities $e_0 \lesssim 0.2$ ten orbits before merger. We use this model to explore in detail the physics that can be extracted with moderately eccentric, non-spinning binary black hole mergers. We use $\texttt{ENIGMA}$ to show that GW150914, GW151226, GW170104, GW170814 and GW170608 can be effectively recovered with spinning, quasi-circular templates if the eccentricity of these events at a gravitational wave frequency of 10Hz satisfies $e_0\leq \{0.175,\, 0.125,\,0.175,\,0.175,\, 0.125\}$, respectively. We show that if these systems have eccentricities $e_0\sim 0.1$ at a gravitational wave frequency of 10Hz, they can be misclassified as quasi-circular binaries due to parameter space degeneracies between eccentricity and spin corrections. Using our catalog of eccentric numerical relativity simulations, we discuss the importance of including higher-order waveform multipoles in gravitational wave searches of eccentric binary black hole mergers.

The space based interferometer LISA will be capable of detecting the gravitational waves emitted by stellar mass black holes or neutron stars slowly inspiralling into the supermassive black holes found in the centre of most galaxies. The gravitational wave signal from such an extreme mass ratio inspiral (EMRI) event will provide a unique opportunity to test whether the spacetime metric around the central black hole is well described by the Kerr solution. In this paper a well studied model for EMRIs around Kerr black holes is extended to a family of parametrically deformed bumpy black holes which preserve the basic symmetries of the Kerr metric. The new EMRI model is then used to quantify the constraints that LISA observations of EMRIs may be able to place on the deviations, or bumps, on the Kerr metric.

The extreme-mass-ratio inspirals (EMRIs) of stellar-mass compact objects into massive black holes are an important class of source for the future space-based gravitational-wave detector LISA. Detecting signals from EMRIs will require waveform models that are both accurate and computationally efficient. In this paper, we present the latest implementation of an augmented analytic kludge (AAK) model, publicly available at github.com/alvincjk/EMRI_Kludge_Suite as part of an EMRI waveform software suite. This version of the AAK model has improved accuracy compared to its predecessors, with two-month waveform overlaps against a more accurate fiducial model exceeding 0.97 for a generic range of sources; it also generates waveforms 5-15 times faster than the fiducial model. The AAK model is well suited for scoping out data analysis issues in the upcoming round of mock LISA data challenges. A simple analytic argument shows that it might even be viable for detecting EMRIs with LISA through a semi-coherent template bank method, while the use of the original analytic kludge in the same approach will result in around 90% fewer detections.

The inspiral of stellar-mass compact objects, like neutron stars or stellar-mass black holes, into supermassive black holes provides a wealth of information about the strong gravitational-field regime via the emission of gravitational waves. In order to detect and analyse these signals, accurate waveform templates which include the effects of the compact object's gravitational self-force are required. For computational efficiency, adiabatic templates are often used. These accurately reproduce orbit-averaged trajectories arising from the first-order self-force, but neglect other effects, such as transient resonances, where the radial and poloidal fundamental frequencies become commensurate. During such resonances the flux of gravitational waves can be diminished or enhanced, leading to a shift in the compact object's trajectory and the phase of the waveform. We present an evolution scheme for studying the effects of transient resonances and apply this to an astrophysically motivated population. We find that a large proportion of systems encounter a low-order resonance in the later stages of inspiral; however, the resulting effect on signal-to-noise recovery is small as a consequence of the low eccentricity of the inspirals. Neglecting the effects of transient resonances leads to a loss of 4% of detectable signals.

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.

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.

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.

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.

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.

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.

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 $\Omega_\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\mu/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\mu/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 $\Omega^\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.

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.

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.

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.

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.

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.

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.

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.

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.

In this article we consider prospects for detecting extreme mass ratio inspirals (EMRIs) using gravitational wave (GW) observations by a future space borne interferometric observatory eLISA. We start with a description of EMRI formation channels. Different formation scenarios lead to variations in the expected event rate and predict different distributions of the orbital parameters when the GW signal enters the eLISA sensitivity band. Then we will briefly overview the available theoretical models describing the GW signal from EMRIs and describe proposed methods for their detection.

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.

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.