
Extrememassratioinspiral observations from future spacebased
gravitationalwave detectors such as LISA will enable strongfield 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 blackhole 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 nullhypothesis test. We adapt and apply
this test to a generalised model for extrememassratio inspirals constructed
on deformed blackhole spacetimes, and focus our investigation on how
computational efficiency can be increased through an evidencefree method of
model selection. This method is akin to the algorithm known as productspace
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 orderofmagnitude 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 nonEinsteinian 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 pulsartiming 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 HellingsDowns curve. Finally, the crosscorrelation
between the redshift and astrometric signal is also calculated; this may form
the basis for future joint pulsartiming and astrometry searches for
arbitrarily polarized gravitational wave backgrounds.

The observation of the electromagnetic counterpart of gravitationalwave (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 gammaray 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 followup 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, inspiralmergerringdown
waveform model that describes nonspinning binary black holes systems that
evolve on moderately eccentric orbits. The inspiral evolution is described
using a consistent combination of postNewtonian theory, selfforce and black
hole perturbation theory. Assuming eccentric binaries that circularize prior to
coalescence, we smoothly match the eccentric inspiral with a standalone,
quasicircular merger, which is constructed using machine learning algorithms
that are trained with quasicircular numerical relativity waveforms. We show
that $\texttt{ENIGMA}$ reproduces with excellent accuracy the dynamics of
quasicircular 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 massratios 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, nonspinning binary black hole mergers. We use
$\texttt{ENIGMA}$ to show that GW150914, GW151226, GW170104, GW170814 and
GW170608 can be effectively recovered with spinning, quasicircular 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
quasicircular 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 higherorder
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 extrememassratio inspirals (EMRIs) of stellarmass compact objects into
massive black holes are an important class of source for the future spacebased
gravitationalwave 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 twomonth waveform
overlaps against a more accurate fiducial model exceeding 0.97 for a generic
range of sources; it also generates waveforms 515 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 semicoherent 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 stellarmass compact objects, like neutron stars or
stellarmass black holes, into supermassive black holes provides a wealth of
information about the strong gravitationalfield 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 selfforce are required. For computational efficiency, adiabatic
templates are often used. These accurately reproduce orbitaveraged
trajectories arising from the firstorder selfforce, 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 loworder resonance in the later stages of
inspiral; however, the resulting effect on signaltonoise 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 matchedfilter
searches for gravitational waves is to introduce a compressed basis for the
waveform template bank in a gridbased 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 spacebased detectors such as eLISA, for which
online grid searches are difficult due to the longduration waveforms and large
parameter spaces.

Gaussian process regression (GPR) is a nonparametric 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 speedup 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)
gravitationalwave signals from binary black holes that could be observed using
advancedera gravitationalwave detectors. Unless properly accounted for,
uncertainty in the gravitationalwave 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 signaltonoise, 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 onedimensional
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 spacebased gravitationalwave 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 extrememassratio inspirals (EMRIs) of stellarmass compact
objects into massive black holes. Previous data challenges have made use of the
wellknown 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 BarackCutler 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 ultraprecise timing data on 41
millisecond pulsars. We develop frequentist and Bayesian detection algorithms
to search both for monochromatic and frequencyevolving 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
skyaveraged 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
subcentiparsec 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 stateoftheart SMBHB
population models predict $<1\%$ probability of detecting a CGW with the
current EPTA dataset, consistent with the reported nondetection. 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 signaltonoise 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 readytouse expressions to compute the
increase/loss in signaltonoise 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:astroph/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 gravitationalwave 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 gravitationalwave 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
signaltonoise ratios in a pulsar timing array.

We present new limits on an isotropic stochastic gravitationalwave
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 supermassive 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
NambuGoto 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 Xrays 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 Xray 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
gravitationalwave backgrounds to backgrounds which have nonGR polarisations.
Our analysis and results are presented in the context of pulsartiming 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 pulsartiming response to gravitational
waves with nonGR polarisation are given for each mode of a spinweighted
sphericalharmonic decomposition of the background, which permit the signal to
be mapped across the sky to any desired resolution. We also derive the
pulsartiming overlap reduction functions for the various nonGR polarisations,
finding analytic forms for anisotropic backgrounds with scalartransverse
("breathing") and vectorlongitudinal polarisations, and a semianalytic form
for scalarlongitudinal backgrounds. Our results indicate that pulsartiming
observations will be completely insensitive to scalartransverse mode
anisotropies in the polarisation amplitude beyond dipole, and anisotropies in
the power beyond quadrupole. Analogously to our previous findings that
pulsartiming observations lack sensitivity to tensorcurl modes for a
transversetraceless tensor background, we also find insensitivity to
vectorcurl modes for a vectorlongitudinal 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 PP 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 PP plot obtained when using
an incorrect signal model to perform parameter estimation. This approximation
is valid in the limit of high signaltonoise ratio and nearly correct waveform
models. We show how the PP 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 timeofarrival 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
postNewtonian 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
gravitationalwave 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 globularcluster shape,
location, and density from an observation of a stellar field that contains a
nonuniform 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
farfield 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 wavewave resonances previously described in
the literature are absent when the electricmagnetic selfinteraction is taken
into account. The fluctuation and amplification of electromagnetic energy flux
as the gravitational wave strength increases towards the
gravitationalelectromagnetic 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 stellarmass mass compact objects into massive black holes
in the centres of galaxies are one of the most important sources of
gravitational radiation for spacebased detectors like LISA or eLISA. These
extrememassratio 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 gravitationalwave
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
gravitationalwave 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 crosscorrelation 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.