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.

The Regge-Wheeler-Zerilli (RWZ) wave-equation describes Schwarzschild-Droste black hole perturbations. The source term contains a Dirac distribution and its derivative. We have previously designed a method of integration in time domain. It consists of a finite difference scheme where analytic expressions, dealing with the wave-function discontinuity through the jump conditions, replace the direct integration of the source and the potential. Herein, we successfully apply the same method to the geodesic generic orbits of EMRI (Extreme Mass Ratio Inspiral) sources, at second order. An EMRI is a Compact Star (CS) captured by a Super Massive Black Hole (SMBH). These are considered the best probes for testing gravitation in strong regime. The gravitational wave-forms, the radiated energy and angular momentum at infinity are computed and extensively compared with other methods, for different orbits (circular, elliptic, parabolic, including zoom-whirl).

We apply our method of indirect integration, described in Part I, at fourth order, to the radial fall affected by the self-force. The Mode-Sum regularisation is performed in the Regge-Wheeler gauge using the equivalence with the harmonic gauge for this orbit. We consider also the motion subjected to a self-consistent and iterative correction determined by the self-force through osculating stretches of geodesics. The convergence of the results confirms the validity of the integration method. This work complements and justifies the analysis and the results appeared in Int. J. Geom. Meth. Mod. Phys., 11, 1450090 (2014).

We adopt the Dirac-Detweiler-Whiting radiative and regular effective field in curved spacetime. Thereby, we derive straightforwardly the first order perturbative correction to the geodesic of the background in a covariant form, for the extreme mass ratio two-body problem. The correction contains the self-force contribution and a background metric dependent term.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

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.

We intend to develop part of the theoretical tools needed for the detection of gravitational waves coming from the capture of a compact object, 1-100 solar masses, by a Supermassive Black Hole, up to a 10 billion solar masses, located at the centre of most galaxies. The analysis of the accretion activity unveils the star population around the galactic nuclei, and tests the physics of black holes and general relativity. The captured small mass is considered a probe of the gravitational field of the massive body, allowing a precise measurement of the particle motion up to the final absorption. The knowledge of the gravitational signal, strongly affected by the self-force - the orbital displacement due to the captured mass and the emitted radiation - is imperative for a successful detection. The results include a strategy for wave equations with a singular source term for all type of orbits. We are now tackling the evolution problem, first for radial fall in Regge- Wheeler gauge, and later for generic orbits in the harmonic or de Donder gauge for Schwarzschild-Droste black holes. In the Extreme Mass Ratio Inspiral, the determination of the orbital evolution demands that the motion of the small mass be continuously corrected by the self-force, i.e. the self-consistent evolution. At each of the integration steps, the self-force must be computed over an adequate number of modes; further, a differential-integral system of general relativistic equations is to be solved and the outputs regularised for suppressing divergences. Finally, for the provision of the computational power, parallelisation is under examination.

The advanced interferometer network will herald a new era in observational astronomy. There is a very strong science case to go beyond the advanced detector network and build detectors that operate in a frequency range from 1 Hz-10 kHz, with sensitivity a factor ten better in amplitude. Such detectors will be able to probe a range of topics in nuclear physics, astronomy, cosmology and fundamental physics, providing insights into many unsolved problems in these areas.

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.

On the basis of a recently proposed strategy of finite element integration in time domain for partial differential equations with a singular source term, we present a fourth order algorithm for non-rotating black hole perturbations in the Regge-Wheeler gauge. Herein, we address even perturbations induced by a particle plunging in. The forward time value at the upper node of the $(r^*,t)$ grid cell is obtained by an algebraic sum of i) the preceding node values of the same cell, ii) analytic expressions, related to the jump conditions on the wave function and its derivatives, iii) the values of the wave function at adjacent cells. In this approach, the numerical integration does not deal with the source and potential terms directly, for cells crossed by the particle world line. This scheme has also been applied to circular and eccentric orbits and it will be object of a forthcoming publication.

Perturbations of Schwarzschild-Droste black holes in the Regge-Wheeler gauge benefit from the availability of a wave equation and from the gauge invariance of the wave function, but lack smoothness. Nevertheless, the even perturbations belong to the C\textsuperscript{0} continuity class, if the wave function and its derivatives satisfy specific conditions on the discontinuities, known as jump conditions, at the particle position. These conditions suggest a new way for dealing with finite element integration in time domain. The forward time value in the upper node of the $(t, r^*$) grid cell is obtained by the linear combination of the three preceding node values and of analytic expressions based on the jump conditions. The numerical integration does not deal directly with the source term, the associated singularities and the potential. This amounts to an indirect integration of the wave equation. The known wave forms at infinity are recovered and the wave function at the particle position is shown. In this series of papers, the radial trajectory is dealt with first, being this method of integration applicable to generic orbits of EMRI (Extreme Mass Ratio Inspiral).

We establish the jump conditions for the wavefunction and its derivatives through the formal solutions of the wave equation. These conditions respond to the requirement of continuity of the perturbations at the position of the particle and they are given for any mode at first order. Using these jump conditions, we then propose a new method for computing the radiated waveform without direct integration of the source term. We consider this approach potentially applicable to generic orbits.

We present the results of a search for gravitational-wave bursts associated with 137 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray experiments during the fifth LIGO science run and first Virgo science run. The data used in this analysis were collected from 2005 November 4 to 2007 October 1, and most of the GRB triggers were from the Swift satellite. The search uses a coherent network analysis method that takes into account the different locations and orientations of the interferometers at the three LIGO-Virgo sites. We find no evidence for gravitational-wave burst signals associated with this sample of GRBs. Using simulated short-duration (<1 s) waveforms, we set upper limits on the amplitude of gravitational waves associated with each GRB. We also place lower bounds on the distance to each GRB under the assumption of a fixed energy emission in gravitational waves, with typical limits of D ~ 15 Mpc (E_GW^iso / 0.01 M_o c^2)^1/2 for emission at frequencies around 150 Hz, where the LIGO-Virgo detector network has best sensitivity. We present astrophysical interpretations and implications of these results, and prospects for corresponding searches during future LIGO-Virgo runs.

We present an up-to-date, comprehensive summary of the rates for all types of compact binary coalescence sources detectable by the Initial and Advanced versions of the ground-based gravitational-wave detectors LIGO and Virgo. Astrophysical estimates for compact-binary coalescence rates depend on a number of assumptions and unknown model parameters, and are still uncertain. The most confident among these estimates are the rate predictions for coalescing binary neutron stars which are based on extrapolations from observed binary pulsars in our Galaxy. These yield a likely coalescence rate of 100 per Myr per Milky Way Equivalent Galaxy (MWEG), although the rate could plausibly range from 1 per Myr per MWEG to 1000 per Myr per MWEG. We convert coalescence rates into detection rates based on data from the LIGO S5 and Virgo VSR2 science runs and projected sensitivities for our Advanced detectors. Using the detector sensitivities derived from these data, we find a likely detection rate of 0.02 per year for Initial LIGO-Virgo interferometers, with a plausible range between 0.0002 and 0.2 per year. The likely binary neutron-star detection rate for the Advanced LIGO-Virgo network increases to 40 events per year, with a range between 0.4 and 400 per year.

Progenitor scenarios for short gamma-ray bursts (short GRBs) include coalescenses of two neutron stars or a neutron star and black hole, which would necessarily be accompanied by the emission of strong gravitational waves. We present a search for these known gravitational-wave signatures in temporal and directional coincidence with 22 GRBs that had sufficient gravitational-wave data available in multiple instruments during LIGO's fifth science run, S5, and Virgo's first science run, VSR1. We find no statistically significant gravitational-wave candidates within a [-5, +1) s window around the trigger time of any GRB. Using the Wilcoxon-Mann-Whitney U test, we find no evidence for an excess of weak gravitational-wave signals in our sample of GRBs. We exclude neutron star-black hole progenitors to a median 90% CL exclusion distance of 6.7 Mpc.

We present a search for gravitational waves from 116 known millisecond and young pulsars using data from the fifth science run of the LIGO detectors. For this search ephemerides overlapping the run period were obtained for all pulsars using radio and X-ray observations. We demonstrate an updated search method that allows for small uncertainties in the pulsar phase parameters to be included in the search. We report no signal detection from any of the targets and therefore interpret our results as upper limits on the gravitational wave signal strength. The most interesting limits are those for young pulsars. We present updated limits on gravitational radiation from the Crab pulsar, where the measured limit is now a factor of seven below the spin-down limit. This limits the power radiated via gravitational waves to be less than ~2% of the available spin-down power. For the X-ray pulsar J0537-6910 we reach the spin-down limit under the assumption that any gravitational wave signal from it stays phase locked to the X-ray pulses over timing glitches, and for pulsars J1913+1011 and J1952+3252 we are only a factor of a few above the spin-down limit. Of the recycled millisecond pulsars several of the measured upper limits are only about an order of magnitude above their spin-down limits. For these our best (lowest) upper limit on gravitational wave amplitude is 2.3x10^-26 for J1603-7202 and our best (lowest) limit on the inferred pulsar ellipticity is 7.0x10^-8 for J2124-3358.

Through detection by low gravitational wave space interferometers, the capture of stars by supermassive black holes will constitute a giant step forward in the understanding of gravitation in strong field. The impact of the perturbations on the motion of the star is computed via the tail, the back-scattered part of the perturbations, or via a radiative Green function. In the former approach, the self-force acts upon the background geodesic, while in the latter, the geodesic is conceived in the total (background plus perturbations) field. Regularisations (mode-sum and Riemann-Hurwitz $\zeta$ function) intervene to cancel divergencies coming from the infinitesimal size of the particle. The non-adiabatic trajectories require the most sophisticated techniques for studying the evolution of the motion, like the self-consistent approach.

This work deals with the motion of a radially falling star in Schwarzschild geometry and correctly identifies radiation reaction terms by the perturbative method. The results are: i) identification of all terms up to first order in perturbations, second in trajectory deviation, and mixed terms including lowest order radiation reaction terms; ii) renormalisation of all divergent terms by the $\zeta$ Riemann and Hurwitz functions. The work implements a method previously identified by one of the authors and corrects some current misconceptions and results.