LISA Pathfinder (LPF) has been a space-based mission designed to test new technologies that will be required for a gravitational wave observatory in space. Magnetically driven forces play a key role in the instrument sensitivity in the low-frequency regime (mHz and below), the measurement band of interest for a space-based observatory. The magnetic field can couple to the magnetic susceptibility and remanent magnetic moment from the test masses and disturb them from their geodesic movement. LISA Pathfinder carried on-board a dedicated magnetic measurement subsystem with noise levels of 10 $ \rm nT \ Hz^{-1/2}$ from 1 Hz down to 1 mHz. In this paper we report on the magnetic measurements throughout LISA Pathfinder operations. We characterise the magnetic environment within the spacecraft, study the time evolution of the magnetic field and its stability down to 20 $\mu$Hz, where we measure values around 200 $ \rm nT \ Hz^{-1/2}$ and identify two different frequency regimes, one related to the interplanetary magnetic field and the other to the magnetic field originating inside the spacecraft. Finally, we characterise the non-stationary component of the fluctuations of the magnetic field below the mHz and relate them to the dynamics of the solar wind.

An idealized "test" object in general relativity moves along a geodesic. However, if the object has a finite mass, this will create additional curvature in the spacetime, causing it to deviate from geodesic motion. If the mass is nonetheless sufficiently small, such an effect is usually treated perturbatively and is known as the gravitational self-force due to the object. This issue is still an open problem in gravitational physics today, motivated not only by basic foundational interest, but also by the need for its direct application in gravitational-wave astronomy. In particular, the observation of extreme-mass-ratio inspirals by the future space-based detector LISA will rely crucially on an accurate modeling of the self-force driving the orbital evolution and gravitational wave emission of such systems. In this paper, we present a novel derivation, based on conservation laws, of the basic equations of motion for this problem. They are formulated with the use of a quasilocal (rather than matter) stress-energy-momentum tensor---in particular, the Brown-York tensor---so as to capture gravitational effects in the momentum flux of the object, including the self-force. Our formulation and resulting equations of motion are independent of the choice of the perturbative gauge. We show that, in addition to the usual gravitational self-force term, they also lead to an additional "self-pressure" force not found in previous analyses, and also that our results correctly recover known formulas under appropriate conditions. Our approach thus offers a fresh geometrical picture from which to understand the self-force fundamentally, and potentially useful new avenues for computing it practically.

We present the scientific motivation for future space tests of the equivalence principle, and in particular the universality of free fall, at the $10^{-17}$ level or better. Two possible mission scenarios, one based on quantum technologies, the other on electrostatic accelerometers, that could reach that goal are briefly discussed.

We report on the results of the LISA Pathfinder (LPF) free-fall mode experiment, in which the control force needed to compensate the quasistatic differential force acting on two test masses is applied intermittently as a series of "impulse" forces lasting a few seconds and separated by roughly 350 s periods of true free fall. This represents an alternative to the normal LPF mode of operation in which this balancing force is applied continuously, with the advantage that the acceleration noise during free fall is measured in the absence of the actuation force, thus eliminating associated noise and force calibration errors. The differential acceleration noise measurement presented here with the free-fall mode agrees with noise measured with the continuous actuation scheme, representing an important and independent confirmation of the LPF result. An additional measurement with larger actuation forces also shows that the technique can be used to eliminate actuation noise when this is a dominant factor.

LISA Pathfinder (LPF) was a technology pioneering mission designed to test key technologies required for gravitational wave detection in space. In the low frequency regime (milli-Hertz and below), where space-based gravitational wave observatories will operate, temperature fluctuations play a crucial role since they can couple into the interferometric measurement and the test masses' free-fall accuracy in many ways. A dedicated temperature measurement subsystem, with noise levels in 10$\,\mu$K$\,$Hz$^{-1/2}$ down to $1\,$mHz was part of the diagnostics unit on board LPF. In this paper we report on the temperature measurements throughout mission operations, characterize the thermal environment, estimate transfer functions between different locations and report temperature stability (and its time evolution) at frequencies as low as 10$\,\mu$Hz, where typically values around $1\,$K$\,$Hz$^{-1/2}$ were measured.

Non-recurrent short term variations of the galactic cosmic-ray (GCR) flux above 70 MeV n$^{-1}$ were observed between 2016 February 18 and 2017 July 3 aboard the European Space Agency LISA Pathfinder (LPF) mission orbiting around the Lagrange point L1 at 1.5$\times$10$^6$ km from Earth. The energy dependence of three Forbush decreases (FDs) is studied and reported here. A comparison of these observations with others carried out in space down to the energy of a few tens of MeV n$^{-1}$ shows that the same GCR flux parameterization applies to events of different intensity during the main phase. FD observations in L1 with LPF and geomagnetic storm occurrence is also presented. Finally, the characteristics of GCR flux non-recurrent variations (peaks and depressions) of duration $<$ 2 days and their association with interplanetary structures are investigated. It is found that, most likely, plasma compression regions between subsequent corotating high-speed streams cause peaks, while heliospheric current sheet crossing cause the majority of the depressions.

Since the 2017 Nobel Prize in Physics was awarded for the observation of gravitational waves, it is fair to say that the epoch of gravitational wave astronomy (GWs) has begun. However, a number of interesting sources of GWs can only be observed from space. To demonstrate the feasibility of the Laser Interferometer Space Antenna (LISA), a future gravitational wave observatory in space, the LISA Pathfinder satellite was launched on December, 3rd 2015. Measurements of the spurious forces accelerating an otherwise free-falling test mass, and detailed investigations of the individual subsystems needed to achieve the free-fall, have been conducted throughout the mission. This overview article starts with the purpose and aim of the mission, explains satellite hardware and mission operations and ends with a summary of selected important results and an outlook towards LISA. From the LISA Pathfinder experience, we can conclude that the proposed LISA mission is feasible.

The inspiral of a stellar-mass compact object into a massive ($\sim 10^{4}$-$10^{7} M_{\odot}$) black hole produces an intricate gravitational-wave signal. Due to the extreme-mass ratios involved, these systems complete $\sim 10^{4}$-$10^{5}$ orbits, most of them in the strong-field region of the massive black hole, emitting in the frequency range $\sim10^{-4}-1~$Hz. This makes them prime sources for the space-based observatory LISA (Laser Interferometer Space Antenna). LISA observations will enable high-precision measurements of the physical characteristics of these extreme-mass-ratio inspirals (EMRIs): redshifted masses, massive black hole spin and orbital eccentricity can be determined with fractional errors $\sim 10^{-4}$-$10^{-6}$, the luminosity distance with better than $\sim 10\%$ precision, and the sky localization to within a few square degrees. EMRIs will provide valuable information about stellar dynamics in galactic nuclei, as well as precise data about massive black hole populations, including the distribution of masses and spins. They will enable percent-level measurements of the multipolar structure of massive black holes, precisely testing the strong-gravity properties of their spacetimes. EMRIs may also provide cosmographical data regarding the expansion of the Universe if inferred source locations can be correlated with galaxy catalogs.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

The science operations of the LISA Pathfinder mission has demonstrated the feasibility of sub-femto-g free-fall of macroscopic test masses necessary to build a LISA-like gravitational wave observatory in space. While the main focus of interest, i.e. the optical axis or the $x$-axis, has been extensively studied, it is also of interest to evaluate the stability of the spacecraft with respect to all the other degrees of freedom. The current paper is dedicated to such a study, with a focus set on an exhaustive and quantitative evaluation of the imperfections and dynamical effects that impact the stability with respect to its local geodesic. A model of the complete closed-loop system provides a comprehensive understanding of each part of the in-loop coordinates spectra. As will be presented, this model gives very good agreements with LISA Pathfinder flight data. It allows one to identify the physical noise source at the origin and the physical phenomena underlying the couplings. From this, the performances of the stability of the spacecraft, with respect to its geodesic, are extracted as a function of frequency. Close to $1 mHz$, the stability of the spacecraft on the $X_{SC}$, $Y_{SC}$ and $Z_{SC}$ degrees of freedom is shown to be of the order of $5.0\ 10^{-15} m\ s^{-2}/\sqrt{Hz}$ for X and $4.0 \ 10^{-14} m\ s^{-2}/\sqrt{Hz}$ for Y and Z. For the angular degrees of freedom, the values are of the order $3\ 10^{-12} rad\ s^{-2}/\sqrt{Hz}$ for $\Theta_{SC}$ and $3\ 10^{-13} rad\ s^{-2}/\sqrt{Hz}$ for $H_{SC}$ and $\Phi_{SC}$.

The Space Technology 7 Disturbance Reduction System (ST7-DRS) is a NASA technology demonstration payload that operated from January 2016 through July of 2017 on the European Space Agency's LISA Pathfinder spacecraft. The joint goal of the NASA and ESA missions was to validate key technologies for a future space-based gravitational wave observatory targeting the source-rich milliHertz band. The two primary components of ST7-DRS are a micropropulsion system based on colloidal micro-Newton thrusters (CMNTs) and a control system that simultaneously controls the attitude and position of the spacecraft and the two free-flying test masses (TMs). This paper presents our main experimental results and summarizes the overall the performance of the CMNTs and control laws. We find that the CMNT performance to be consistent with pre-flight predictions, with a measured system thrust noise on the order of $100\,\textrm{nN}/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band. The control system maintained the TM-spacecraft separation with an RMS error of less than 2$\,$nm and a noise spectral density of less than $3\,\textrm{nm}/\sqrt{\textrm{Hz}}$ in the same band. Thruster calibration measurements yield thrust values consistent with the performance model and ground-based thrust-stand measurements, to within a few percent. We also report a differential acceleration noise between the two test masses with a spectral density of roughly $3\,\textrm{fm}/\textrm{s}^2/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band, slightly less than twice as large as the best performance reported with the baseline LISA Pathfinder configuration and below the current requirements for the Laser Interferometer Space Antenna (LISA) mission.

LISA Pathfinder (LPF) was a European Space Agency mission with the aim to test key technologies for future space-borne gravitational-wave observatories like LISA. The main scientific goal of LPF was to demonstrate measurements of differential acceleration between free-falling test masses at the sub-femto-g level, and to understand the residual acceleration in terms of a physical model of stray forces, and displacement readout noise. A key step toward reaching the LPF goals was the correct calibration of the dynamics of LPF, which was a three-body system composed by two test-masses enclosed in a single spacecraft, and subject to control laws for system stability. In this work, we report on the calibration procedures adopted to calculate the residual differential stray force per unit mass acting on the two test-masses in their nominal positions. The physical parameters of the adopted dynamical model are presented, together with their role on LPF performance. The analysis and results of these experiments show that the dynamics of the system was accurately modeled and the dynamical parameters were stationary throughout the mission. Finally, the impact and importance of calibrating system dynamics for future space-based gravitational wave observatories is discussed.

Galactic cosmic-ray (GCR) energy spectra observed in the inner heliosphere are modulated by the solar activity, the solar polarity and structures of solar and interplanetary origin. A high counting rate particle detector (PD) aboard LISA Pathfinder (LPF), meant for subsystems diagnostics, was devoted to the measurement of galactic cosmic-ray and solar energetic particle integral fluxes above 70 MeV n$^{-1}$ up to 6500 counts s$^{-1}$. PD data were gathered with a sampling time of 15 s. Characteristics and energy-dependence of GCR flux recurrent depressions and of a Forbush decrease dated August 2, 2016 are reported here. The capability of interplanetary missions, carrying PDs for instrument performance purposes, in monitoring the passage of interplanetary coronal mass ejections (ICMEs) is also discussed.

The space-based Laser Interferometer Space Antenna (LISA) will be able to observe the gravitational-wave signals from systems comprised of a massive black hole and a stellar-mass compact object. These systems are known as extreme-mass-ratio inspirals (EMRIs) and are expected to complete $\sim 10^4$-$10^5$ cycles in band, thus allowing exquisite measurements of their parameters. In this work, we attempt to quantify the astrophysical uncertainties affecting the predictions for the number of EMRIs detectable by LISA, and find that competing astrophysical assumptions produce a variance of about three orders of magnitude in the expected intrinsic EMRI rate. However, we find that irrespective of the astrophysical model, at least a few EMRIs per year should be detectable by the LISA mission, with up to a few thousands per year under the most optimistic astrophysical assumptions. We also investigate the precision with which LISA will be able to extract the parameters of these sources. We find that typical fractional statistical errors with which the intrinsic parameters (redshifted masses, massive black hole spin and orbital eccentricity) can be recovered are $\sim 10^{-6}$-$10^{-4}$. Luminosity distance (which is required to infer true masses) is inferred to about $10\%$ precision and sky position is localized to a few square degrees, while tests of the multipolar structure of the Kerr metric can be performed to percent-level precision or better.

Test mass charging caused by cosmic rays will be a significant source of acceleration noise for space-based gravitational wave detectors like LISA. Operating between December 2015 and July 2017, the technology demonstration mission LISA Pathfinder included a bespoke monitor to help characterise the relationship between test mass charging and the local radiation environment. The radiation monitor made in situ measurements of the cosmic ray flux while also providing information about its energy spectrum. We describe the monitor and present measurements which show a gradual 40% increase in count rate coinciding with the declining phase of the solar cycle. Modulations of up to 10% were also observed with periods of 13 and 26 days that are associated with co-rotating interaction regions and heliospheric current sheet crossings. These variations in the flux above the monitor detection threshold (approximately 70 MeV) are shown to be coherent with measurements made by the IREM monitor on-board the Earth orbiting INTEGRAL spacecraft. Finally we use the measured deposited energy spectra, in combination with a GEANT4 model, to estimate the galactic cosmic ray differential energy spectrum over the course of the mission.

One of the key astrophysical sources for the Laser Interferometer Space Antenna (LISA) are the inspirals of stellar-origin compact objects into massive black holes in the centres of galaxies. These extreme-mass-ratio inspirals (EMRIs) have great potential for astrophysics, cosmology and fundamental physics. In this paper we describe the likely numbers and properties of EMRI events that LISA will observe. We present the first results computed for the 2.5 Gm interferometer that was the new baseline mission submitted in January 2017 in response to the ESA L3 mission call. In addition, we attempt to quantify the astrophysical uncertainties in EMRI event rate estimates by considering a range of different models for the astrophysical population. We present both likely event rates and estimates for the precision with which the parameters of the observed sources could be measured. We finish by discussing the implications of these results for science using EMRIs.

The gravitational waves emitted by binary systems with extreme mass ratios carry unique astrophysical information expected to be probed by the next generation of gravitational wave detectors such as LISA. The detection of these binaries rely on an accurate modeling of the gravitational self-force that drives their orbital evolution. Although the theoretical formalism to compute the self-force has been largely established, the mathematical tools needed to implement it are still under development, and the self-force computation remains an open problem. We present here a frequency-domain implementation of the particle-without-particle (PwP) technique previously developed for the computation of the scalar self-force -- a helpful testbed for the gravitational self-force.

Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA 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 infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz.

We report on electrostatic measurements made on board the European Space Agency mission LISA Pathfinder. Detailed measurements of the charge-induced electrostatic forces exerted on free-falling test masses (TMs) inside the capacitive gravitational reference sensor are the first made in a relevant environment for a space-based gravitational wave detector. Employing a combination of charge control and electric-field compensation, we show that the level of charge-induced acceleration noise on a single TM can be maintained at a level close to 1.0 fm/s^2/sqrt(Hz) across the 0.1-100 mHz frequency band that is crucial to an observatory such as LISA. Using dedicated measurements that detect these effects in the differential acceleration between the two test masses, we resolve the stochastic nature of the TM charge build up due to interplanetary cosmic rays and the TM charge-to-force coupling through stray electric fields in the sensor. All our measurements are in good agreement with predictions based on a relatively simple electrostatic model of the LISA Pathfinder instrument.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

The LISA Pathfinder mission will demonstrate the technology of drag-free test masses for use as inertial references in future space-based gravitational wave detectors. To accomplish this, the Pathfinder spacecraft will perform drag-free flight about a test mass while measuring the acceleration of this primary test mass relative to a second reference test mass. Because the reference test mass is contained within the same spacecraft, it is necessary to apply forces on it to maintain its position and attitude relative to the spacecraft. These forces are a potential source of acceleration noise in the LISA Pathfinder system that are not present in the full LISA configuration. While LISA Pathfinder has been designed to meet it's primary mission requirements in the presence of this noise, recent estimates suggest that the on-orbit performance may be limited by this `suspension noise'. The drift-mode or free-flight experiments provide an opportunity to mitigate this noise source and further characterize the underlying disturbances that are of interest to the designers of LISA-like instruments. This article provides a high-level overview of these experiments and the methods under development to analyze the resulting data.

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.

The gravitational waves emitted by binary systems with extreme-mass ratios carry unique astrophysical information that can only be detected by space-based detectors like eLISA. To that end, a very accurate modelling of the system is required. The gravitational self-force program, which has been fully developed in the Lorenz gauge, is the best approach we have so far. However, the computations required would be done more efficiently if we could work in other gauges, like the Regge-Wheeler (RW) one in the case of Schwarzschild black holes. In this letter we present a new scheme, based on the Particle-without-Particle formulation of the field equations, where the gravitational self-force can be obtained from just solving individual wave-type equations like the master equations of the RW gauge. This approach can help to tackle the yet unsolved Kerr case.

The main goal of the LISA Pathfinder (LPF) mission is to fully characterize the acceleration noise models and to test key technologies for future space-based gravitational-wave observatories similar to the eLISA concept. The data analysis team has developed complex three-dimensional models of the LISA Technology Package (LTP) experiment on-board LPF. These models are used for simulations, but more importantly, they will be used for parameter estimation purposes during flight operations. One of the tasks of the data analysis team is to identify the physical effects that contribute significantly to the properties of the instrument noise. A way of approaching this problem is to recover the essential parameters of a LTP model fitting the data. Thus, we want to define the simplest model that efficiently explains the observations. To do so, adopting a Bayesian framework, one has to estimate the so-called Bayes Factor between two competing models. In our analysis, we use three main different methods to estimate it: The Reversible Jump Markov Chain Monte Carlo method, the Schwarz criterion, and the Laplace approximation. They are applied to simulated LPF experiments where the most probable LTP model that explains the observations is recovered. The same type of analysis presented in this paper is expected to be followed during flight operations. Moreover, the correlation of the output of the aforementioned methods with the design of the experiment is explored.

Neutron stars and black holes are the most compact astrophysical objects we can think of and as a consequence they are the main sources of gravitational waves. There are many astrophysically relevant scenarios in which these objects are immersed in or endowed with strong magnetic fields, in such a way that gravitational perturbations can couple to electromagnetic ones and can potentially trigger synergistic electromagnetic signatures. In a recent paper we derived the main equations for gravito-electromagnetic perturbations and studied in detail the case of polar electromagnetic perturbations driven by axial gravitational perturbations. In this paper we deal with the case of axial electromagnetic perturbations driven by polar black-hole or neutron stars oscillations, in which the energy emitted in case is considerably larger than in the previous case. In the case of neutron stars the phenomenon lasts considerably longer since the fluid acts as an energy reservoir that shakes the magnetic field for a timescale of the order of secs.

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete binary black hole hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a "blind injection challenge" similar to that conducted in recent LIGO and Virgo science runs, we added 7 hybrid waveforms to two months of data recolored to predictions of Advanced LIGO and Advanced Virgo sensitivity curves during their first observing runs. The resulting data was analyzed by gravitational-wave detection algorithms and 6 of the waveforms were recovered with false alarm rates smaller than 1 in a thousand years. Parameter estimation algorithms were run on each of these waveforms to explore the ability to constrain the masses, component angular momenta and sky position of these waveforms. We also perform a large-scale monte-carlo study to assess the ability to recover each of the 60 hybrid waveforms with early Advanced LIGO and Advanced Virgo sensitivity curves. Our results predict that early Advanced LIGO and Advanced Virgo will have a volume-weighted average sensitive distance of 300Mpc (1Gpc) for $10M_{\odot}+10M_{\odot}$ ($50M_{\odot}+50M_{\odot}$) binary black hole coalescences. We demonstrate that neglecting the component angular momenta in the waveform models used in matched-filtering will result in a reduction in sensitivity for systems with large component angular momenta. [Abstract abridged for ArXiv, full version in PDF]

LISA Pathfinder is a mission planned by the European Space Agency to test the key technologies that will allow the detection of gravitational waves in space. The instrument on-board, the LISA Technology package, will undergo an exhaustive campaign of calibrations and noise characterisation campaigns in order to fully describe the noise model. Data analysis plays an important role in the mission and for that reason the data analysis team has been developing a toolbox which contains all the functionalities required during operations. In this contribution we give an overview of recent activities, focusing on the improvements in the modelling of the instrument and in the data analysis campaigns performed both with real and simulated data.

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.

Gravitational perturbations of neutron stars and black holes are well known sources of gravitational radiation. If the compact object is immersed in or endowed with a magnetic field, the gravitational perturbations would couple to electromagnetic perturbations and potentially trigger synergistic electromagnetic signatures. We present a detailed analytic calculation of the dynamics of coupled gravitational and electromagnetic perturbations for both neutron stars and black holes. We discuss the prospects for detecting the electromagnetic waves in these scenarios and the potential that these waves have for providing information about their source.

A common result to all EMRI investigations on rates is that the possibility that a compact object merges with the MBH after only one intense burst of GWs is much more likely than a slow adiabatic inspiral, an EMRI. The later is referred to as a "plunge" because the compact object dives into the MBH. The event rates for plunges are orders of magnitude larger than slow inspirals. On the other hand, nature MBH's are most likely Kerr and the magnitude of the spin has been sized up to be high. We calculate the number of periapsis passages that a compact object set on to an extremely radial orbit goes through before being actually swallowed by the Kerr MBH and we then translate it into an event rate for a LISA-like observatory, such as the proposed ESA mission eLISA/NGO. We prove that a "plunging" compact object is conceptually indistinguishable from an adiabatic, slow inspiral; plunges spend on average up to hundred of thousands of cycles in the bandwidth of the detector for a two years mission. This has an important impact on the event rate, enhancing in some cases significantly, depending on the spin of the MBH and the inclination. Moreover, it has been recently proved that the production of low-eccentricity EMRIs is severely blocked by the presence of a blockade in the rate at which orbital angular momenta change takes place. This is the result of relativistic precession on to the stellar potential torques and hence affects EMRIs originating via resonant relaxation at distances of about $\sim 10^{-2}$ pc from the MBH. Since high-eccentricity EMRIs are a result of two-body relaxation, they are not affected by this phenomenon. Therefore we predict that eLISA EMRI event rates will be dominated by high-eccentricity binaries, as we present here.

The capture of a compact object in a galactic nucleus by a massive black hole (MBH), an extreme-mass ratio inspiral (EMRI), is the best way to map space and time around it. Recent work on stellar dynamics has demonstrated that there seems to be a complot in phase space acting on low-eccentricity captures, since their rates decrease significantly by the presence of a blockade in the rate at which orbital angular momenta change takes place. This so-called "Schwarzschild barrier" is a result of the impact of relativistic precession on to the stellar potential torques, and thus it affects the enhancement on lower-eccentricity EMRIs that one would expect from resonant relaxation. We confirm and quantify the existence of this barrier using a statistical sample of 2,500 direct-summation N-body simulations using both a post-Newtonian and also for the first time in a direct-summation code a geodesic approximation for the relativistic orbits. The existence of the barrier prevents low-eccentricity EMRIs from approaching the central MBH, but high-eccentricity EMRIs, which have been wrongly classified as "direct plunges" until recently, ignore the presence of the barrier, because they are driven by two-body relaxation. Hence, since the rates are significantly affected in the case of low-eccentricity EMRIs, we predict that a LISA-like observatory such as eLISA will predominantly detect high-eccentricity EMRIs.

The capture of a compact object in a galactic nucleus by a massive black hole (MBH) is the best way to map space and time around it. Compact objects such as stellar black holes on a capture orbit with a very high eccentricity have been wrongly assumed to be lost for the system after an intense burst of radiation, which has been described as a "direct plunge". We prove that these very eccentric capture orbits spend actually a similar number of cycles in a LISA-like detector as those with lower eccentricities if the central MBH is spinning. Although the rates are higher for high-eccentricity EMRIs, the spin also enhances the rates of lower-eccentricity EMRIs. This last kind have received more attention because of the fact that high-eccentricity EMRIs were thought to be direct plunges and thus negligible. On the other hand, recent work on stellar dynamics has demonstrated that there seems to be a complot in phase space acting on these lower-eccentricity captures, since their rates decrease significantly by the presence of a blockade in the rate at which orbital angular momenta change takes place. This so-called "Schwarzschild barrier" is a result of the impact of relativistic precession on to the stellar potential torques, and thus it affects the enhancement on lower-eccentricity EMRIs that one would expect from resonant relaxation. We confirm and quantify the existence of this barrier using a statitical sample of 2,500 direct-summation N-body simulations using both a post-Newtonian but also, and for the first time, a geodesic approximation for the relativistic orbits. The existence of the barrier prevents "traditional EMRIs" from approaching the central MBH, but if the central MBH is spinning the rate will be anyway dominated by highly-eccentric extreme-mass ratio inspirals, which insolently ignore the presence of the barrier, because they are driven by two-body relaxation.

The capture of a stellar-mass compact object by a supermassive black hole and the subsequent inspiral (driven by gravitational radiation emission) constitute one of the most important sources of gravitational waves for space-based observatories like eLISA/NGO. In this article we describe their potential as high-precision tools that can be used to perform tests of the geometry of black holes and also of the strong field regime of gravity.

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

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

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.

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

We describe a new kludge scheme to model the dynamics of generic extreme-mass-ratio inspirals (EMRIs; stellar compact objects spiraling into a spinning supermassive black hole) and their gravitational-wave emission. The Chimera scheme is a hybrid method that combines tools from different approximation techniques in General Relativity: (i) A multipolar, post-Minkowskian expansion for the far-zone metric perturbation (the gravitational waveforms) and for the local prescription of the self-force; (ii) a post-Newtonian expansion for the computation of the multipole moments in terms of the trajectories; and (iii) a BH perturbation theory expansion when treating the trajectories as a sequence of self-adjusting Kerr geodesics. The EMRI trajectory is made out of Kerr geodesic fragments joined via the method of osculating elements as dictated by the multipolar post-Minkowskian radiation-reaction prescription. We implemented the proper coordinate mapping between Boyer-Lindquist coordinates, associated with the Kerr geodesics, and harmonic coordinates, associated with the multipolar post-Minkowskian decomposition. The Chimera scheme is thus a combination of approximations that can be used to model generic inspirals of systems with extreme to intermediate mass ratios, and hence, it can provide valuable information for future space-based gravitational-wave observatories, like LISA, and even for advanced ground detectors. The local character in time of our multipolar post-Minkowskian self-force makes this scheme amenable to study the possible appearance of transient resonances in generic inspirals.

Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011.

We introduce a new kludge scheme to model the dynamics of generic extreme mass-ratio inspirals (stellar compact objects spiraling into a spinning supermassive black hole) and to produce the gravitational waveforms that describe the gravitational-wave emission of these systems. This scheme combines tools from different techniques in General Relativity: It uses a multipolar, post-Minkowskian (MPM) expansion for the far-zone metric perturbation (which provides the gravitational waveforms, here taken up to mass hexadecapole and current octopole order) and for the local prescription of the self-force (since we are lacking a general prescription for it); a post-Newtonian expansion for the computation of the multipole moments in terms of the trajectories; and a BH perturbation theory expansion when treating the trajectories as a sequence of self-adjusting Kerr geodesics. The orbital evolution is thus equivalent to solving the geodesic equations with time-dependent orbital elements, as dictated by the MPM radiation-reaction prescription. To complete the scheme, both the orbital evolution and wave generation require to map the Boyer-Lindquist coordinates of the orbits to the harmonic coordinates in which the different MPM quantities have been derived, a mapping that we provide explicitly in this paper. This new kludge scheme is thus a combination of approximations that can be used to model generic inspirals of systems with extreme mass ratios to systems with more moderate mass ratios, and hence can provide valuable information for future space-based gravitational-wave observatories like LISA and even for advanced ground detectors. Finally, due to the local character in time of our MPM self-force, this scheme can be used to perform studies of the possible appearance of transient resonances in generic inspirals.

The calculation of the self force in the modeling of the gravitational-wave emission from extreme-mass-ratio binaries is a challenging task. Here we address the question of the possible emergence of a persistent spurious solution in time-domain schemes, referred to as a {\em Jost junk solution} in the literature, that may contaminate self force calculations. Previous studies suggested that Jost solutions are due to the use of zero initial data, which is inconsistent with the singular sources associated with the small object, described as a point mass. However, in this work we show that the specific origin is an inconsistency in the translation of the singular sources into jump conditions. More importantly, we identify the correct implementation of the sources at late times as the sufficient condition guaranteeing the absence of Jost junk solutions.

When a stellar-mass compact object is captured by a supermassive black hole located in a galactic centre, the system losses energy and angular momentum by the emission of gravitational waves. Subsequently, the stellar compact object evolves inspiraling until plunging onto the massive black hole. These EMRI systems are expected to be one of the main sources of gravitational waves for the future space-based Laser Interferometer Space Antenna (LISA). However, the detection of EMRI signals will require of very accurate theoretical templates taking into account the gravitational self-force, which is the responsible of the stellar-compact object inspiral. Due to its potential applicability on EMRIs, the obtention of an efficient method to compute the scalar self-force acting on a point-like particle orbiting around a massive black hole is being object of increasing interest. We present here a review of our time-domain numerical technique to compute the self-force acting on a point-like particle and we show its suitability to deal with both circular and eccentric orbits.

The computation of the self-force constitutes one of the main challenges for the construction of precise theoretical waveform templates in order to detect and analyze extreme-mass-ratio inspirals with the future space-based gravitational-wave observatory LISA. Since the number of templates required is quite high, it is important to develop fast algorithms both for the computation of the self-force and the production of waveforms. In this article we show how to tune a recent time-domain technique for the computation of the self-force, what we call the Particle without Particle scheme, in order to make it very precise and at the same time very efficient. We also extend this technique in order to allow for highly eccentric orbits.

The emergent area of gravitational wave astronomy promises to provide revolutionary discoveries in the areas of astrophysics, cosmology, and fundamental physics. One of the most exciting possibilities is to use gravitational-wave observations to test alternative theories of gravity. In this contribution we describe how to use observations of extreme-mass-ratio inspirals by the future Laser Interferometer Space Antenna to test a particular class of theories: Chern-Simons modified gravity.

The gravitational-wave signals emitted by Extreme-Mass-Ratio Inspirals will be hidden in the instrumental LISA noise and the foreground noise produced by galactic binaries in the LISA band. Then, we need accurate gravitational-wave templates to extract these signals from the noise and obtain the relevant physical parameters. This means that in the modeling of these systems we have to take into account how the orbit of the stellar-mass compact object is modified by the action of its own gravitational field. This effect can be described as the action of a local force, the self-force. We present a time-domain technique to compute the self-force for geodesic eccentric orbits around a non-rotating massive black hole. To illustrate the method we have applied it to a testbed model consisting of scalar charged particle orbiting a non-dynamical black hole. A key feature of our method is that it does not introduce a small scale associated with the stellar-mass compact object. This is achieved by using a multidomain framework where the particle is located at the interface between two subdomains. In this way, we just have to evolve homogeneous wave-like equations with smooth solutions that have to be communicated across the subdomain boundaries using appropriate junction conditions. The numerical technique that we use to implement this scheme is the pseudospectral collocation method. We show the suitability of this technique for the modeling of Extreme-Mass-Ratio Inspirals and show that it can provide accurate results for the self-force.

The LISA PathFinder DMU (Data Management Unit) flight model was formally accepted by ESA and ASD on 11 February 2010, after all hardware and software tests had been successfully completed. The diagnostics items are scheduled to be delivered by the end of 2010. In this paper we review the requirements and performance of this instrumentation, specially focusing on the Radiation Monitor and the DMU, as well as the status of their programmed use during mission operations, on which work is ongoing at the time of writing.

GW Notes was born from the need for a journal where the distinct communities involved in gravitation wave research might gather. While these three communities - Astrophysics, General Relativity and Data Analysis - have made significant collaborative progress over recent years, we believe that it is indispensable to future advancement that they draw closer, and that they speak a common idiom. For this GW Notes issue we have approached Carlos F. Sopuerta to write the highlight article of this issue on his thoughts on fundamental physics with LISA.

[abridged] The inspiral of a stellar compact object into a massive black hole is one of the main sources of gravitational waves for the future space-based Laser Interferometer Space Antenna. We expect to be able to detect and analyze many cycles of these slowly inspiraling systems. To that end, the use of very precise theoretical waveform templates in the data analysis is required. To build them we need to have a deep understanding of the gravitational backreaction mechanism responsible for the inspiral. The self-force approach describes the inspiral as the action of a local force that can be obtained from the regularization of the perturbations created by the stellar compact object on the massive black hole geometry. In this paper we extend a new time-domain technique for the computation of the self-force from the circular case to the case of eccentric orbits around a non-rotating black hole. The main idea behind our scheme is to use a multidomain framework in which the small compact object, described as a particle, is located at the interface between two subdomains. Then, the equations at each subdomain are homogeneous wave-type equations, without distributional sources. In this particle-without-particle formulation, the solution of the equations is smooth enough to provide good convergence properties for the numerical computations. This formulation is implemented by using a pseudospectral collocation method for the spatial discretization, combined with a Runge Kutta algorithm for the time evolution. We present results from several simulations of eccentric orbits in the case of a scalar charged particle around a Schwarzschild black hole. In particular, we show the convergence of the method and its ability to resolve the field and its derivatives across the particle location. Finally, we provide numerical values of the self-force for different orbital parameters.

We introduce a new time-domain method for computing the self-force acting on a scalar particle in a Schwarzschild geometry. The principal feature of our method consists in the division of the spatial domain into several subdomains and locating the particle at the interface betweem two them. In this way, we avoid the need of resolving a small length scale associated with the presence of a particle in the computational domain and, at the same time, we avoid numerical problems due to the low differentiability of solutions of equations with point-like singular behaviour.

Testing deviation of GR is one of the main goals of the proposed {\emph{Laser Interferometer Space Antenna}}, a space-based gravitational-wave observatory. For the first time, we consistently compute the generation of gravitational waves from extreme-mass ratio inspirals (stellar compact objects into supermassive black holes) in a well-motivated alternative theory of gravity, that to date remains weakly constrained by double binary pulsar observations. The theory we concentrate on is Chern-Simons (CS) modified gravity, a 4-D, effective theory that is motivated both from string theory and loop-quantum gravity, and which enhances the Einstein-Hilbert action through the addition of a dynamical scalar field and the parity-violating Pontryagin density. We show that although point particles continue to follow geodesics in the modified theory, the background about which they inspiral is a modification to the Kerr metric, which imprints a CS correction on the gravitational waves emitted. CS modified gravitational waves are sufficiently different from the General Relativistic expectation that they lead to significant dephasing after 3 weeks of evolution, but such dephasing will probably not prevent detection of these signals, but instead lead to a systematic error in the determination of parameters. We end with a study of radiation-reaction in the modified theory and show that, to leading-order, energy-momentum emission is not CS modified, except possibly for the subdominant effect of scalar-field emission. The inclusion of radiation-reaction will allow for tests of CS modified gravity with space-borne detectors that might be two orders of magnitude larger than current binary pulsar bounds.

[abridged] Chern-Simons (CS) modified gravity is a 4D effective theory that descends both from string theory and loop quantum gravity, and that corrects the Einstein-Hilbert action by adding the product of a scalar field and the parity-violating, Pontryagin density. In this theory, the gravitational field of spinning black holes is described by a modified Kerr geometry whose multipole moments deviate from the Kerr ones only at the fourth multipole, l = 4. We investigate possible signatures of this theory in the gravitational wave emission produced in the inspiral of stellar compact objects into massive black holes, both for intermediate- and extreme-mass ratios. We use the semi-relativistic approximation, where the trajectories are geodesics of the massive black hole geometry and the gravitational waveforms are obtained from a multipolar decomposition of the radiative field. The main CS corrections to the waveforms arise from modifications to the geodesic trajectories, due to changes to the massive black hole geometry, and manifest themselves as an accumulating dephasing relative to the general relativistic case. We also explore the propagation and the stress-energy tensor of gravitational waves in this theory. We find that, although this tensor has the same form as in General Relativity, the energy and angular momentum balance laws are indeed modified through the stress-energy tensor of the CS scalar field. These balance laws could be used to describe the inspiral through adiabatic changes in the orbital parameters, which in turn would enhance the dephasing effect. Gravitational-wave observations of intermediate- or extreme-mass ratio inspirals with advanced ground detectors or with LISA could use such dephasing to test the dynamical theory to unprecedented levels.

The description of the inspiral of a stellar-mass compact object into a massive black hole sitting at a galactic centre is a problem of major relevance for the future space-based gravitational-wave observatory LISA (Laser Interferometer Space Antenna), as the signals from these systems will be buried in the data stream and accurate gravitational-wave templates will be needed to extract them. The main difficulty in describing these systems lies in the estimation of the gravitational effects of the stellar-mass compact object on his own trajectory around the massive black hole, which can be modeled as the action of a local force, the self-force. In this paper, we present a new time-domain numerical method for the computation of the self-force in a simplified model consisting of a charged scalar particle orbiting a nonrotating black hole. We use a multi-domain framework in such a way that the particle is located at the interface between two domains so that the presence of the particle and its physical effects appear only through appropriate boundary conditions. In this way we eliminate completely the presence of a small length scale associated with the need of resolving the particle. This technique also avoids the problems associated with the impact of a low differentiability of the solution in the accuracy of the numerical computations. The spatial discretization of the field equations is done by using the pseudospectral collocation method and the time evolution, based on the method of lines, uses a Runge-Kutta solver. We show how this special framework can provide very efficient and accurate computations in the time domain, which makes the technique amenable for the intensive computations required in the astrophysically-relevant scenarios for LISA.

Gravitational Wave Astronomy is becoming a reality as Earth-based interferometric gravitational-wave detectors reach the design sensitivities and move towards advanced configurations that may lead to gravitational-wave detections in the immediate future. In this contribution, I briefly summarize the basic characteristics of this new area, the discovery prospects and the potential for fundamental physics. Then, I present results of some investigations of two different sources of gravitational waves that are potential targets for present and future planned observatories. First, I will discuss the generation of gravitational radiation by non-linear effects arising from the coupling between radial and non-radial oscillations of neutron stars, which may produce distinctive gravitational-wave signatures. The gravitational radiation emitted by these sources is in the frequency band of Earth-based detectors. And second, I will discuss the gravitational-wave emission during the inspiral of extreme-mass-ratio compact binaries. In this case, the gravitational waves have low frequencies, inside the frequency band of space observatories like LISA.

We study perturbations of a Schwarzschild black hole in Chern-Simons modified gravity. We begin by showing that Birkhoff's theorem holds for a wide family of Chern-Simons coupling functions, a scalar field present in the theory that controls the strength of the Chern-Simons correction to the Einstein-Hilbert action. After decomposing the perturbations in spherical harmonics, we study the linearized modified field equations and find that axial and polar modes are coupled, in contrast to general relativity. The divergence of the modified equations leads to the Pontryagin constraint, which forces the vanishing of the Cunningham-Price-Moncrief master function associated with axial modes. We analyze the structure of these equations and find that the appearance of the Pontryagin constraint yields an overconstrained system that does not allow for generic black hole oscillations. We illustrate this situation by studying the case characterized by a canonical choice of the coupling function and pure-parity perturbative modes. We end with a discussion of how to extend Chern-Simons modified gravity to bypass the Pontryagin constraint and the suppression of perturbations.

Extreme mass ratio bursts (EMRBs) have been proposed as a possible source for future space-borne gravitational wave detectors, such as the Laser Interferometer Space Antenna (LISA). These events are characterized by long-period, nearly-radial orbits of compact objects around a central massive black hole. The gravitational radiation emitted during such events consists of a short burst, corresponding to periapse passage, followed by a longer, silent interval. In this paper we investigate the impact of including relativistic corrections to the description of the compact object's trajectory via a geodesic treatment, as well as including higher-order multipole corrections in the waveform calculation. The degree to which the relativistic corrections are important depends on the EMRB's orbital parameters. We find that relativistic EMRBs (v_{max}}/c > 0.25) are not rare and actually account for approximately half of the events in our astrophysical model. The relativistic corrections tend to significantly change the waveform amplitude and phase relative to a Newtonian description, although some of this dephasing could be mimicked by parameter errors. The dephasing over several bursts could be of particular importance not only to gravitational wave detection, but also to parameter estimation, since it is highly correlated to the spin of the massive black hole. Consequently, we postulate that if a relativistic EMRB is detected, such dephasing might be used to probe the relativistic character of the massive black hole and obtain information about its spin.

[abridged] The coalescence of a binary black hole system is one of the main sources of gravitational waves that present and future detectors will study. Apart from the energy and angular momentum that these waves carry, for unequal-mass binaries there is also a net flux of linear momentum that implies a recoil velocity of the resulting final black hole in the opposite direction. We present a computation of the recoil velocity based on the close-limit approximation scheme, which gives excellent results for head-on and grazing collisions of black holes when compared to full numerical relativistic calculations. We obtain a maximum recoil velocity of ~ 64 km/s for a symmetric mass ratio eta = M_1 M_2/(M_1+M_2)^2 ~ 0.19 and an initial proper separation of 4 M, where M is the total ADM mass of the system. This separation is the maximum at which the close-limit approximation is expected to provide accurate results. If we supplement this estimate with PN calculations up to the innermost stable circular orbit, we obtain a lower bound for the recoil velocity, with a maximum around 84 km/s. This is a lower bound because it neglects the initial merger phase. We can however obtain a rough estimate by using PN methods or the close-limit approximation. Since both methods are known to overestimate the amount of radiation, we obtain in this way an upper bound for the recoil with maxima in the range of 220-265 km/s. We also provide non-linear fits to these estimated upper and lower bounds. These estimates are subject to uncertainties related to issues such as the choice of initial data and higher effects in perturbation theory. Nonetheless, our estimates are consistent with previous results in the literature and suggest a narrower range of possible recoil velocities.

The formation and growth of supermassive black holes is a key issue to unveil the secrets of galaxy formation. In particular, the gravitational recoil produced in the merger of unequal mass black hole binaries could have a number of astrophysical implications, such as the ejection of black holes from the host galaxy or globular cluster. We present estimates of the recoil velocity that include the effect of small eccentricities. The approach is specially suited for the last stage of the merger, where most of the emission of linear momentum in gravitational waves takes place. Supplementing our estimates with post-Newtonian approximations, we obtain lower and upper bounds that constrain previous recoil velocities estimates as well as a best estimate that agrees with numerical simulations in the quasi-circular case. For eccentricities e <= 0.1, the maximum recoil is found for mass ratios of M_1/M_2 ~ 0.38 with a best estimate of ~ 167 (1 + e) km/s and upper and lower bounds of 79 (1 + e) km/s and 216 (1 + e) km/s respectively.

This is a brief report on time-domain numerical simulations of extreme-mass-ratio binaries based on finite element methods. We discuss a new technique for solving the perturbative equations describing a point-like object orbiting a non-rotating massive black hole and the prospects of using it for the evaluation of the gravitational self-force responsible of the inspiral of these binary systems. We also discuss the perspectives of transferring this technology to the more astrophysically relevant case of a central rotating massive black hole.

We introduce a computational framework which avoids solving explicitly hydrodynamic equations and is suitable to study the pre-merger evolution of black hole-neutron star binary systems. The essence of the method consists of constructing a neutron star model with a black hole companion and freezing the internal degrees of freedom of the neutron star during the course of the evolution of the space-time geometry. We present the main ingredients of the framework, from the formulation of the problem to the appropriate computational techniques to study these binary systems. In addition, we present numerical results of the construction of initial data sets and evolutions that demonstrate the feasibility of this approach.

We investigate the non-linear coupling between radial and non-radial oscillations of static spherically symmetric neutron stars as a possible mechanism for the generation of gravitational waves that may lead to observable signatures. In this paper we concentrate on the axial sector of the non-radial perturbations. By using a multi-parameter perturbative framework we introduce a complete description of the non-linear coupling between radial and axial non-radial oscillations; we study the gauge invariant character of the associated perturbative variables and develop a computational scheme to evolve the non-linear coupling perturbations in the time domain. We present results of simulations corresponding to different physical situations and discuss the dynamical behaviour of this non-linear coupling. Of particular interest is the occurrence of signal amplifications in the form of resonance phenomena when a frequency associated with the radial pulsations is close to a frequency associated with one of the axial w-modes of the star. Finally, we mention possible extensions of this work and improvements towards more astrophysically motivated scenarios.

The description of extreme-mass-ratio binary systems in the inspiral phase is a challenging problem in gravitational wave physics with significant relevance for the space interferometer LISA. The main difficulty lies in the evaluation of the effects of the small body's gravitational field on itself. To that end, an accurate computation of the perturbations produced by the small body with respect the background geometry of the large object, a massive black hole, is required. In this paper we present a new computational approach based on Finite Element Methods to solve the master equations describing perturbations of non-rotating black holes due to an orbiting point-like object. The numerical computations are carried out in the time domain by using evolution algorithms for wave-type equations. We show the accuracy of the method by comparing our calculations with previous results in the literature. Finally, we discuss the relevance of this method for achieving accurate descriptions of extreme-mass-ratio binaries.

Linear perturbation theory is appropriate to describe small oscillations of stars, while a mild non-linearity is still tractable perturbatively but requires to consider mode coupling. It is natural to start to look at this problem by considering the coupling between linear radial and non-radial modes. Radial pulsations of a spherical compact objects do not per se emit gravitational waves but, if the coupling is efficient in driving and possibly amplifying the non-radial modes, gravitational radiation could then be produced to a significant level. In this paper we develop the relativistic formalism to study the coupling of radial and non-radial first order perturbations of a compact spherical star. From a mathematical point of view, it is convenient to treat the two sets of perturbations as separately parametrized, using a 2-parameter perturbative expansion of the metric, the energy-momentum tensor and Einstein equations in which $\lambda$ is associated with the radial modes, $\epsilon$ with the non-radial perturbations, and the $\lambda\epsilon$ terms describe the coupling. This approach provides a well-defined framework to consider the gauge dependence of perturbations, allowing us to use $\epsilon$ order gauge-invariant non-radial variables on the static background and to define new second order $\lambda\epsilon$ gauge-invariant variables describing the non-linear coupling. We present the evolution and constraint equations for our variables outlining the setup for numerical computations, and briefly discuss the surface boundary conditions in terms of the second order $\lambda\epsilon$ Lagrangian pressure perturbation.

This is a progress report on our study of the coupling of first-order radial and non-radial relativistic perturbations of a static spherical star. Our goal is to investigate the effects of this coupling on the gravitational wave signal of neutron stars. In particular, we are looking for the existence of resonances and parametric amplifications, changes in the damping time of non-radial oscillations, etc. To that end, we have developed a formalism that introduces gauge invariant quantities to describe the coupling. Their equations have the same structure as the equations for first-order non-radial perturbations plus some source terms, which makes them very appealing for time domain studies.

If the observable universe is a braneworld of Randall-Sundrum type, then particle interactions at high energies will produce 5-dimensional gravitons that escape into the bulk. As a result, the Weyl energy density on the brane does not behave like radiation in the early universe, but does so only later, in the low energy regime. Recently a simple model was proposed to describe this modification of the Randall-Sundrum cosmology. We investigate the dynamics of this model, and find the exact solution of the field equations. We use a dynamical systems approach to analyze global features of the phase space of solutions.

We introduce N-parameter perturbation theory as a new tool for the study of non-linear relativistic phenomena. The main ingredient in this formulation is the use of the Baker-Campbell-Hausdorff formula. The associated machinery allows us to prove the main results concerning the consistency of the scheme to any perturbative order. Gauge transformations and conditions for gauge invariance at any required order can then be derived from a generating exponential formula via a simple Taylor expansion. We outline the relation between our novel formulation and previous developments.

In this paper we study thick-shell braneworld models in the presence of a Gauss-Bonnet term. We discuss the peculiarities of the attainment of the thin-shell limit in this case and compare them with the same situation in Einstein gravity. We describe the two simplest families of thick-brane models (parametrized by the shell thickness) one can think of. In the thin-shell limit, one family is characterized by the constancy of its internal density profile (a simple structure for the matter sector) and the other by the constancy of its internal curvature scalar (a simple structure for the geometric sector). We find that these two families are actually equivalent in Einstein gravity and that the presence of the Gauss-Bonnet term breaks this equivalence. In the second case, a shell will always keep some non-trivial internal structure, either on the matter or on the geometric sectors, even in the thin-shell limit.

In this paper we consider the Long-Wavelength Approximation Scheme (LWAS) in the framework of the covariant fluid approach to general relativistic dynamics, specializing to the particular case of irrotational dust matter. We discuss the dynamics of these models during the approach to any spacelike singularity where a BKL-type evolution is expected, studying the validity of this approximation scheme and the role of the magnetic part of the Weyl tensor, $H_{ab}$. Our analytic results confirm a previous numerical analysis: it is $H_{ab}$ that destroys the pure Kasner-like approach to the singularity and eventually produces the bounce to another Kasner phase. Expanding regions evolve as separate universes where inhomogeneities and anisotropies die away.

An implicit fundamental assumption in relativistic perturbation theory is that there exists a parametric family of spacetimes that can be Taylor expanded around a background. The choice of the latter is crucial to obtain a manageable theory, so that it is sometime convenient to construct a perturbative formalism based on two (or more) parameters. The study of perturbations of rotating stars is a good example: in this case one can treat the stationary axisymmetric star using a slow rotation approximation (expansion in the angular velocity Omega), so that the background is spherical. Generic perturbations of the rotating star (say parametrized by lambda) are then built on top of the axisymmetric perturbations in Omega. Clearly, any interesting physics requires non-linear perturbations, as at least terms lambda Omega need to be considered. In this paper we analyse the gauge dependence of non-linear perturbations depending on two parameters, derive explicit higher order gauge transformation rules, and define gauge invariance. The formalism is completely general and can be used in different applications of general relativity or any other spacetime theory.

We report results on the construction of cosmological braneworld models in the context of the Einstein-Gauss-Bonnet gravity, which include the leading correction to the Einstein-Hilbert action suggested by superstring theory. We obtain and study the equations governing the dynamics of the standard cosmological models. We find that they can be written in the same form as in the case of the Randall-Sundrum model but with time-varying four-dimensional gravitational and cosmological constants. Finally, we discuss the cosmological evolution predicted by these models and their compatibility with observational data.

We consider braneworld scenarios including the leading correction to the Einstein-Hilbert action suggested by superstring theory, the Gauss-Bonnet term. We obtain and study the complete set of equations governing the cosmological dynamics. We find they have the same form as those in Randall-Sundrum scenarios but with time-varying four-dimensional gravitational and cosmological constants. Studying the bulk geometry we show that this variation is produced by bulk curvature terms parameterized by the mass of a black hole. Finally, we show there is a coupling between these curvature terms and matter that can be relevant for early universe cosmology.

In this paper we deal with the cosmological dynamics of Randall-Sundrum brane-world type scenarios in which the five-dimensional Weyl tensor has a non-vanishing projection onto the three-brane where matter fields are confined. Using dynamical systems techniques, we study how the state space of Friedmann-Lemaitre-Robertson-Walker (FLRW) and Bianchi type I cosmological models is affected by the bulk Weyl tensor, focusing in the differences that appear with respect to general relativity and also Randall-Sundrum cosmological scenarios without the Weyl tensor contribution.

We present the results of a systematic investigation of the qualitative behaviour of the Friedmann-Lemaitre-Robertson-Walker (FLRW) and Bianchi I and V cosmological models in Randall-Sundrum brane-world type scenarios.

In this work we consider Randall-Sundrum brane-world type scenarios, in which the spacetime is described by a five-dimensional manifold with matter fields confined in a domain wall or three-brane. We present the results of a systematic analysis, using dynamical systems techniques, of the qualitative behaviour of the Friedmann-Lemaitre-Robertson-Walker and the Bianchi I and V cosmological models in these scenarios. We construct the state spaces for these models and discuss how their structure changes with respect to the general-relativistic case, in particular, what new critical points appear and their nature, the occurrence of bifurcations and the dynamics of anisotropy.

The imposition of symmetries or special geometric properties on submanifolds is less restrictive than to impose them in the full space-time. Starting from this idea, in this paper we study irrotational dust cosmological models in which the geometry of the hypersurfaces generated by the fluid velocity is flat, which supposes a relaxation of the restrictions imposed by the Cosmological Principle. The method of study combines covariant and tetrad methods that exploits the geometrical and physical properties of these models. This procedure will allow us to determine all the space-times within this class as well as to study their properties. Some important consequences and applications of this study are also discussed.

In a cold matter universe, the linearized gravito-magnetic tensor field satisfies a transverse condition (vanishing divergence) when it is purely radiative. We show that in the nonlinear theory, it is no longer possible to maintain the transverse condition, since it leads to a non-terminating chain of integrability conditions. These conditions are highly restrictive, and are likely to hold only in models with special symmetries, such as the known Bianchi and $G_2$ examples. In models with realistic inhomogeneity, the gravito-magnetic field is necessarily non-transverse at second and higher order.