Systems consisting of a massive black hole (MBH) and a stellar-origin compact object (CO), known as extreme mass-ratio inspirals (EMRIs), are of great significance for space-based gravitational-wave detectors, as they will allow for testing gravitational theories in the strong field regime, and for checking the validity of the black hole no-hair theorem. In this work, we present a calculation of the EMRI rate and parameter estimation capabilities of the TianQin observatory, for various astrophysical models for these sources. We find that TianQin can observe EMRIs involving COs with mass of 10$M_\odot$ up to redshift $\sim2$. We also find that detections could reach tens or hundreds per year in the most optimistic astrophysical scenarios. Intrinsic parameters are expected to be recovered to within fractional errors of $\sim 10^{-6}$, while typical errors on the luminosity distance and sky localization are 10% and 10 deg$^2$, respectively. TianQin observation of EMRIs can also constrain possible deviations from the Kerr quadrupole moment to within fractional errors $\lesssim10^{-4}$. We also find that a network of multiple detectors would allow for improvements in both detection rates (by a factor $\sim 1.5$--$3$) and in parameter estimation precision (20-fold improvement for the sky localization and 5-fold improvement for the other parameters.)

We study the impact of gas accretion on the orbital evolution of black-hole binaries initially at large separation in the band of the planned Laser Interferometer Space Antenna (LISA). We focus on two sources: (i)~stellar-origin black-hole binaries~(SOBHBs) that can migrate from the LISA band to the band of ground-based gravitational-wave observatories within weeks/months; and (ii) intermediate-mass black-hole binaries~(IMBHBs) in the LISA band only. Because of the large number of observable gravitational-wave cycles, the phase evolution of these systems needs to be modeled to great accuracy to avoid biasing the estimation of the source parameters. Accretion affects the gravitational-wave phase at negative ($-4$) post-Newtonian order, and is therefore dominant for binaries at large separations. If accretion takes place at the Eddington or at super-Eddington rate, it will leave a detectable imprint on the dynamics of SOBHBs. In optimistic astrophysical scenarios, a multiwavelength strategy with LISA and a ground-based interferometer can detect about $10$ (a few) SOBHB events for which the accretion rate can be measured at $50\%$ ($10\%$) level. In all cases the sky position can be identified within much less than $0.4\,{\rm deg}^2$ uncertainty. Likewise, accretion at $\gtrsim 10\%$ ($\gtrsim 100\%$) of the Eddington rate can be measured in IMBHBs up to redshift $z\approx 0.1$ ($z\approx 0.5$), and the position of these sources can be identified within less than $0.01\,{\rm deg}^2$ uncertainty. Altogether, a detection of SOBHBs or IMBHBs would allow for targeted searches of electromagnetic counterparts to black-hole mergers in gas-rich environments with future X-ray detectors (such as Athena) and radio observatories (such as SKA).

We study scalar fields in a black hole background and show that, when the scalar is suitably coupled to curvature, rapid rotation can induce a tachyonic instability. This instability, which is the hallmark of spontaneous scalarization in the linearized regime, is expected to be quenched by nonlinearities and endow the black hole with scalar hair. Hence, our results demonstrate the existence of a broad class of theories that share the same stationary black hole solutions with general relativity at low spins, but which exhibit black hole hair at sufficiently high spins ($a/M\gtrsim 0.64$). This result has clear implications for tests of general relativity and the nature of black holes with gravitational and electromagnetic observations.

We revisit the predictions for the merger rate of massive black hole binaries detectable by the Laser Interferometer Space Antenna (LISA) and their background signal for pulsar-timing arrays. We focus on the effect of the delays between the merger of galaxies and the final coalescence of black hole binaries, and on the effect of supernova feedback on the growth of black holes. By utilizing a semi-analytic galaxy formation model, not only do we account for the processes that drive the evolution of binaries at separations $\lesssim 1$ pc (gas-driven migration, stellar hardening and triple/quadruple massive black hole systems), but we also improve on previous studies by accounting for the time spent by massive black hole pairs from kpc down to a few pc separation. We also include the effect of supernova feedback, which may eject a substantial amount of gas from the nuclear region of low-mass galaxies, thus hampering the growth of black holes via accretion and suppressing their orbital migration in circumbinary disks. In spite of the inclusion of these novel physical effects, we predict that the LISA detection rate should still be in excess of $\sim 2 \mbox{yr}^{-1}$, irrespective of the model for the seeds of the black hole population at high redshifts. However, scenarios in which black holes form from $\sim100 M_\odot$ seeds are more significantly impacted by the feedback from supernovae. We also present predictions for the mass ratio distribution of the merger population, and find that binaries typically have mass ratios between $\sim 0.1$ and $1$. Predictions for the stochastic background in the band of pulsar-timing array experiments are instead rather robust, and show only a mild dependence on the model.

We consider the observation of stellar-mass black holes binaries with the Laser Interferometer Space Antenna (LISA). Preliminary results based on Fisher information matrix analyses have suggested that gravitational waves from those sources could be very sensitive to possible deviations from the theory of general relativity and from the strong equivalence principle during the low-frequency binary inspiral. We perform a full Markov Chain Monte Carlo Bayesian analysis to quantify the sensitivity of these signals to two phenomenological modifications of general relativity, namely a putative gravitational dipole emission and a non-zero mass for the graviton, properly accounting for the detector's response. Moreover, we consider a scenario where those sources could be observed also with Earth-based detectors, which should measure the coalescence time with precision better than $1 \ {\rm ms}$. This constraint on the coalescence time further improves the bounds that we can set on those phenomenological deviations from general relativity. We show that tests of dipole radiation and the graviton's mass should improve respectively by seven and half an order(s) of magnitude over current bounds. Finally, we discuss under which conditions one may claim the detection of a modification to general relativity.

Based on the rate of resolved stellar origin black hole and neutron star mergers measured by LIGO and Virgo, it is expected that these detectors will also observe an unresolved Stochastic Gravitational Wave Background (SGWB) by the time they reach design sensitivity. Since the same binaries observed by LIGO and Virgo pass through the LISA mHz frequency band at an earlier stage of their orbital evolution,it is foreseen that their SGWB will also be observable by LISA with Signal to Noise Ratio (SNR) $\sim 350$. Unlike the stochastic signal from Galactic white dwarf binaries, for which a subtraction is expected to be possible by exploiting its yearly modulation (induced by the motion of the LISA constellation), the background from unresolved stellar origin black hole and neutron star binaries acts as a foreground for other stochastic signals of cosmological or astrophysical origin, which may also be present in the LISA band. Here, we employ a principal component analysis to model and extract an additional hypothetical SGWB in the LISA band, without making any a priori assumptions on its spectral shape. At the same time, we account for the presence of the foreground from stellar origin black holes and neutron stars, as well as for possible uncertainties in the LISA noise calibration. We find that our technique leads to a linear problem and is therefore suitable for fast and reliable extraction of SGWBs with SNR up to ten times weaker than the LIGO/Virgo foreground, quite independently of their spectral shape.

Photons propagating in a plasma acquire an effective mass $\mu$, which is given by the plasma frequency and which scales with the square root of the plasma density. As noted previously by Conlon and Herdeiro, for electron number densities $n_e\sim 10^{-3}$ cm$^{-3}$ (such as those measured in the interstellar medium) the effective mass induced by the plasma is $\mu \sim 10^{-12}$ eV. This would cause superradiant instabilities for spinning black holes of a few tens of solar masses. An obvious problem with this picture is that densities in the vicinity of black holes are much higher than in the interstellar medium because of accretion. We have conducted numerical simulations of the superradiant instability in spinning black holes surrounded by a plasma with density increasing closer to the black hole, in order to mimic the effect of accretion. While we confirm that superradiant instabilities appear for plasma densities that are sufficiently low near the black hole, we find that astrophysically realistic accretion disks are unlikely to trigger such instabilities.

We investigate the ability of the Laser Interferometer Space Antenna (LISA) to measure the center of mass acceleration of stellar-origin black hole binaries emitting gravitational waves. Our analysis is based on the idea that the acceleration of the center of mass induces a time variation in the redshift of the gravitational wave, which in turn modifies its waveform. We confirm that while the cosmological acceleration is too small to leave a detectable imprint on the gravitational waveforms observable by LISA, larger peculiar accelerations may be measurable for sufficiently long lived sources. We focus on stellar mass black hole binaries, which will be detectable at low frequencies by LISA and near coalescence by ground based detectors. These sources may have large peculiar accelerations, for instance, if they form in nuclear star clusters or in AGN accretion disks. If that is the case, we find that in an astrophysical population calibrated to the LIGO-Virgo observed merger rate, LISA will be able to measure the peculiar acceleration of a small but significant fraction of the events if the mission lifetime is extended beyond the nominal duration of 4 years. In this scenario LISA will be able to assess whether black hole binaries form close to galactic centers, particularly in AGN disks, and will thus help discriminate between different formation mechanisms. Although for a nominal 4 years LISA mission the peculiar acceleration effect cannot be measured, a consistent fraction of events may be biased by strong peculiar accelerations which, if present, may imprint large systematic errors on some waveform parameters. In particular, estimates of the luminosity distance could be strongly biased and consequently induce large systematic errors on LISA measurements of the Hubble constant with stellar mass black hole binaries.

Tidal disruptions are extremely powerful phenomena that have been designated as candidate sources of ultra-high-energy cosmic rays. The disruption of a star by a black hole can naturally provide protons and heavier nuclei, which can be injected and accelerated to ultra-high energies within a jet. Inside the jet, accelerated nuclei are likely to interact with a dense photon field, leading to a significant production of neutrinos and secondary particles. We model numerically the propagation and interactions of high-energy nuclei in jetted tidal disruption events in order to evaluate consistently their signatures in cosmic rays and neutrinos. We propose a simple model of the light curve of tidal disruption events, consisting of two stages: a high state with bright luminosity and short duration and a medium state, less bright and longer lasting. These two states have different impacts on the production of cosmic rays and neutrinos. In order to calculate the diffuse fluxes of cosmic rays and neutrinos, we model the luminosity function and redshift evolution of jetted tidal disruption events. We find that we can fit the latest ultra-high-energy cosmic-ray spectrum and composition results of the Auger experiment for a range of reasonable parameters. The diffuse neutrino flux associated with this scenario is found to be subdominant, but nearby events can be detected by IceCube or next-generation detectors such as IceCube-Gen2.

Ho\v{r}ava gravity breaks boost invariance in the gravitational sector by introducing a preferred time foliation. The dynamics of this preferred slicing is governed, in the low-energy limit suitable for most astrophysical applications, by three dimensionless parameters $\alpha$, $\beta$ and $\lambda$. The first two of these parameters are tightly bound by solar system and gravitational wave propagation experiments, but $\lambda$ remains relatively unconstrained ($0\leq\lambda\lesssim 0.01-0.1$). We restrict here to the parameter space region defined by $\alpha=\beta=0$ (with $\lambda$ kept generic), which in a previous paper we showed to be the only one where black hole solutions are non-pathological at the universal horizon, and we focus on possible violations of the strong equivalence principle in systems involving neutron stars. We compute neutron star 'sensitivities', which parametrize violations of the strong equivalence principle at the leading post-Newtonian order, and find that they vanish identically, like in the black hole case, for $\alpha=\beta=0$ and generic $\lambda\neq0$. This implies that no violations of the strong equivalence principle (neither in the conservative sector nor in gravitational wave fluxes) can occur at the leading post-Newtonian order in binaries of compact objects, and that data from binary pulsars and gravitational interferometers are unlikely to further constrain $\lambda$.

Future gravitational-wave observations will enable unprecedented and unique science in extreme gravity and fundamental physics answering questions about the nature of dynamical spacetimes, the nature of dark matter and the nature of compact objects.

We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $\mu$-Hz frequency range. By the 2040s', the $\mu$-Hz frequency band, bracketed in between the Laser Interferometer Space Antenna (LISA) and pulsar timing arrays, will constitute the largest gap in the coverage of the astrophysically relevant GW spectrum. Yet many outstanding questions related to astrophysics and cosmology are best answered by GW observations in this band. We show that a $\mu$-Hz GW detector will be a truly overarching observatory for the scientific community at large, greatly extending the potential of LISA. Conceived to detect massive black hole binaries from their early inspiral with high signal-to-noise ratio, and low-frequency stellar binaries in the Galaxy, this instrument will be a cornerstone for multimessenger astronomy from the solar neighbourhood to the high-redshift Universe.

We investigate the prospects of detecting gravitational waves from coalescing massive black hole binaries in the Universe with the TianQin observatory, a space-based gravitational wave interferometer proposed to be launched in the 2030s. To frame the scientific scope of the mission, in this paper, we carry out a preliminary estimation of the signal-to-noise ratio, detection rate and parameter estimation precision of massive black hole binaries detectable by TianQin. In order to make our results as robust as possible, we consider several models of the growth history of massive black holes, exploring the effect of some key astrophysical prescriptions as well the impact of the employed computational methods. In the most optimistic model, TianQin can detect as many as approximately 60 mergers per year. If TianQin detects a merger at redshift of 15, it will be capable of estimating its luminosity distance to within an accuracy of 10%; for a nearby event at redshift approximately 2, TianQin can issue early warnings 24 hours before coalescence, with a timing accuracy of around three hours and a sky localization ability of approximately 80 deg$^2$, thus enabling multi-messenger observations.

Modifications of General Relativity leave their imprint both on the cosmic expansion history through a non-trivial dark energy equation of state, and on the evolution of cosmological perturbations in the scalar and in the tensor sectors. In particular, the modification in the tensor sector gives rise to a notion of gravitational-wave (GW) luminosity distance, different from the standard electromagnetic luminosity distance, that can be studied with standard sirens at GW detectors such as LISA or third-generation ground based experiments. We discuss the predictions for modified GW propagation from some of the best studied theories of modified gravity, such as Horndeski or the more general degenerate higher order scalar-tensor (DHOST) theories, non-local infrared modifications of gravity, bigravity theories and the corresponding phenomenon of GW oscillation, as well as theories with extra or varying dimensions. We show that modified GW propagation is a completely generic phenomenon in modified gravity. We then use a simple parametrization of the effect in terms of two parameters $(\Xi_0,n)$, that is shown to fit well the results from a large class of models, to study the prospects of observing modified GW propagation using supermassive black hole binaries as standard sirens with LISA. We construct mock source catalogs and perform detailed Markov Chain Monte Carlo studies of the likelihood obtained from LISA standard sirens alone, as well as by combining them with CMB, BAO and SNe data to reduce the degeneracies between cosmological parameters. We find that the combination of LISA with the other cosmological datasets allows one to measure the parameter $\Xi_0$ that characterizes modified GW propagation to the percent level accuracy, sufficient to test several modified gravity theories. [Abridged]

The recent detections of gravitational waves from binary systems of black holes are in remarkable agreement with the predictions of General Relativity. In this pedagogical mini-review, I will go through the physics of the different phases of the evolution of black hole binary systems, providing a qualitative physical interpretation of each one of them. I will also briefly describe how these phases would be modified if gravitation were described by a theory extending or deforming General Relativity, or if the binary components turned out to be more exotic compact objects than black holes.

Coalescing massive black hole binaries (MBHBs) of $10^{4-7} \rm M_{\odot}$, forming in the aftermath of galaxy mergers, are primary targets of the space mission LISA, the {\it Laser Interferometer Space Antenna}. An assessment of LISA detection prospects requires an estimate of the abundance and properties of MBHBs that form and evolve during the assembly of cosmic structures. To this aim, we employ a semi-analytic model to follow the co-evolution of MBHBs within their host galaxies. We identify three major evolutionary channels driving the binaries to coalescence: two standard paths along which the binary evolution is driven by interactions with the stellar and/or gaseous environment, and a novel channel where MBHB coalescence occurs during the interaction with a third black hole. For each channel, we follow the orbital evolution of MBHBs with physically motivated models that include a self-consistent treatment of the orbital eccentricity. We find that LISA will detect between $\approx 25$ and $\approx 75$ events per year depending on the seed model. We show that triple-induced coalescences can range from a few detected events up to $\sim 30\%$ of the total detected mergers. Moreover, even if the standard gas/stars-driven evolutionary channels should fail and MBHBs were to stall, triple interactions would still occur as a result of the hierarchical nature of galaxy formation, resulting in about $\approx 10$ to $\approx 20$ LISA detections per year. Remarkably, triple interactions among the black holes can produce coalescing binaries with large eccentricities ($\gtrsim 0.9$) upon entrance into the LISA band. This eccentricity will remain significant ($\sim 0.1$) also at merger, requiring suitable templates for parameter estimation.

Coalescing, massive black-hole (MBH) binaries are the most powerful sources of gravitational waves (GWs) in the Universe, which makes MBH science a prime focus for ongoing and upcoming GW observatories. The Laser Interferometer Space Antenna (LISA) -- a gigameter scale space-based GW observatory -- will grant us access to an immense cosmological volume, revealing MBHs merging when the first cosmic structures assembled in the Dark Ages. LISA will unveil the yet unknown origin of the first quasars, and detect the teeming population of MBHs of $10^4 - 10^7$ solar masses. forming within protogalactic halos. The Pulsar Timing Array, a galactic-scale GW survey, can access the largest MBHs the Universe, detecting the cosmic GW foreground from inspiraling MBH binaries of about 10^9 solar masses. LISA can measure MBH spins and masses with precision far exceeding that from electromagnetic (EM) probes, and together, both GW observatories will provide the first full census of binary MBHs, and their orbital dynamics, across cosmic time. Detecting the loud gravitational signal of these MBH binaries will also trigger alerts for EM counterpart searches, from decades (PTAs) to hours (LISA) prior to the final merger. By witnessing both the GW and EM signals of MBH mergers, precious information will be gathered about the rich and complex environment in the aftermath of a galaxy collision. The unique GW characterization of MBHs will shed light on the deep link between MBHs of $10^4-10^{10}$ solar masses and the grand design of galaxy assembly, as well as on the complex dynamics that drive MBHs to coalescence.

Low-frequency gravitational-wave astronomy can perform precision tests of general relativity and probe fundamental physics in a regime previously inaccessible. A space-based detector will be a formidable tool to explore gravity's role in the cosmos, potentially telling us if and where Einstein's theory fails and providing clues about some of the greatest mysteries in physics and astronomy, such as dark matter and the origin of the Universe.

The upcoming LISA mission offers the unique opportunity to study the Milky Way through gravitational wave radiation from Galactic binaries. Among the variety of Galactic gravitational wave sources, LISA is expected to individually resolve signals from $\sim 10^5$ ultra-compact double white dwarf (DWD) binaries. DWDs detected by LISA will be distributed across the Galaxy, including regions that are hardly accessible to electromagnetic observations such as the inner part of the Galactic disc, the bulge and beyond. We quantitatively show that the large number of DWD detections will allow us to use these systems as tracers of the Milky Way potential. We demonstrate that density profiles of DWDs detected by LISA may provide constraints on the scale length parameters of the baryonic components that are both accurate and precise, with statistical errors of a few percent to $10$ percent level. Furthermore, the LISA sample is found to be sufficient to disentangle between different (commonly used) disc profiles, by well covering the disc out to sufficiently large radii. Finally, up to $\sim 80$ DWDs can be detected through both electromagnetic and gravitational wave radiation. This enables multi-messenger astronomy with DWD binaries and allows one to extract their physical properties using both probes. We show that fitting the Galactic rotation curve constructed using distances inferred from gravitational waves {\it and} proper motions from optical observations yield a unique and competitive estimate of the bulge mass. Instead robust results for the stellar disc mass are contingent upon knowledge of the Dark Matter content.

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 Laser Interferometer Space Antenna (LISA) will observe black hole binaries of stellar origin during their gravitational wave inspiral, months to years before coalescence. Due to the long duration of the signal in the LISA band, a faithful waveform is necessary in order to keep track of the binary phase. This is crucial to extract the signal from the data and to perform an unbiased estimation of the source parameters. We consider Post-Newtonian (PN) waveforms, and analyze the PN order needed to keep the bias caused by the PN approximation negligible relative to the statistical parameter estimation error, as a function of the source parameters. By considering realistic population models, we conclude that for $\sim 90\%$ of the stellar black hole binaries detectable by LISA, waveforms at low Post-Newtonian (PN) order (PN $\le 2$) are sufficiently accurate for an unbiased recovery of the source parameters. Our results provide a first estimate of the trade-off between waveform accuracy and information recovery for this class of LISA sources.

The upcoming LISA mission is the only experiment that will allow us to study the Milky Way's structure using gravitational wave signals from Galactic double white dwarfs (DWDs). The total number of expected detections exceeds $10^5$. Furthermore, up to a hundred DWDs can be simultaneously detected in both gravitational and optical radiation (e.g. with Gaia and LSST as eclipsing), making DWDs ideal sources for performing a multi-messenger tomography of the Galaxy. We show that LISA will detect DWDs everywhere, mapping also the opposite side of the Galaxy. This complete coverage will: (1) provide precise and unbiased constraints on the scale radii of the Milky Way's bulge and disc, and (2) allow us to compute the rotation curve and derive competitive estimates for the bulge and disc masses, when combining gravitational wave and optical observations.

Gravitational-wave astronomy has the potential to explore one of the deepest and most puzzling aspects of Einstein's theory: the existence of black holes. A plethora of ultracompact, horizonless objects have been proposed to arise in models inspired by quantum gravity. These objects may solve Hawking's information-loss paradox and the singularity problem associated with black holes, while mimicking almost all of their classical properties. They are, however, generically unstable on relatively short timescales. Here, we show that this "ergoregion instability" leads to a strong stochastic background of gravitational waves, at a level detectable by current and future gravitational-wave detectors. The absence of such background in the first observation run of Advanced LIGO already imposes the most stringent limits to date on black-hole alternatives, showing that certain models of "quantum-dressed" stellar black holes can be at most a small percentage of the total population. The future LISA mission will allow for similar constraints on supermassive black-hole mimickers.

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.

Massive black hole binaries (MBHBs) are expected to form at the centre of merging galaxies during the hierarchical assembly of the cosmic structure, and are expected to be the loudest sources of gravitational waves (GWs) in the low frequency domain. However, because of the dearth of energy exchanges with background stars and gas, many of these MBHBs may stall at separations too large for GW emission to drive them to coalescence in less than a Hubble time. Triple MBH systems are then bound to form after a further galaxy merger, triggering a complex and rich dynamics that can eventually lead to MBH coalescence. Here we report on the results of a large set of numerical simulations, where MBH triplets are set in spherical stellar potentials and MBH dynamics is followed through 2.5 post-Newtonian order in the equations of motion. From our full suite of simulated systems we find that a fraction $\simeq 20-30$% of the MBH binaries that would otherwise stall are led to coalesce within a Hubble time. The corresponding coalescence timescale peaks around 300 Myr, while the eccentricity close to the plunge, albeit small, is non-negligible ($\lesssim 0.1$). We construct and discuss marginalised probability distributions of the main parameters involved and, in a companion paper of the series, we will use the results presented here to forecast the contribution of MBH triplets to the GW signal in the nHz regime probed by Pulsar Timing Array experiments.

Inspiraling massive black-hole binaries (MBHBs) forming in the aftermath of galaxy mergers are expected to be the loudest gravitational-wave (GW) sources relevant for pulsar-timing arrays (PTAs) at nHz frequencies. The incoherent overlap of signals from a cosmic population of MBHBs gives rise to a stochastic GW background (GWB) with characteristic strain around $h_c\sim10^{-15}$ at a reference frequency of 1 yr$^{-1}$, although uncertainties around this value are large. Current PTAs are piercing into the GW amplitude range predicted by MBHB-population models, but no detection has been reported so far. To assess the future success prospects of PTA experiments, it is therefore important to estimate the minimum GWB level consistent with our current understanding of the formation and evolution of galaxies and massive black holes (MBHs). To this purpose, we couple a semianalytic model of galaxy evolution and an extensive study of the statistical outcome of triple MBH interactions. We show that even in the most pessimistic scenario where all MBHBs stall before entering the GW-dominated regime, triple interactions resulting from subsequent galaxy mergers inevitably drive a considerable fraction of the MBHB population to coalescence. At frequencies relevant for PTA, the resulting GWB is only a factor of 2-to-3 suppressed compared to a fiducial model where binaries are allowed to merge over Gyr timescales. Coupled with current estimates of the expected GWB amplitude range, our findings suggest that the minimum GWB from cosmic MBHBs is unlikely to be lower than $h_c\sim 10^{-16}$ (at $f = 1$ yr$^{-1}$), well within the expected sensitivity of projected PTAs based on future observations with FAST, MeerKAT and SKA.

We highlight some subtleties that affect naive implementations of quadrupolar and octupolar gravitational waveforms from numerically-integrated trajectories of three-body systems. Some of those subtleties arise from the requirement that the source be contained in its "coordinate near zone" when applying the standard PN formulae for gravitational-wave emission, and from the need to use the non-linear Einstein equations to correctly derive the quadrupole emission formula. We show that some of these subtleties were occasionally overlooked in the literature, with consequences for published results. We also provide prescriptions that lead to correct and robust predictions for the waveforms computed from numerically-integrated orbits.

Ultralight bosons can induce superradiant instabilities in spinning black holes, tapping their rotational energy to trigger the growth of a bosonic condensate. Possible observational imprints of these boson clouds include (i) direct detection of the nearly monochromatic (resolvable or stochastic) gravitational waves emitted by the condensate, and (ii) statistically significant evidence for the formation of "holes" at large spins in the spin versus mass plane (sometimes also referred to as "Regge plane") of astrophysical black holes. In this work, we focus on the prospects of LISA and LIGO detecting or constraining scalars with mass in the range $m_s\in [10^{-19},\,10^{-15}]$ eV and $m_s\in [10^{-14},\,10^{-11}]$ eV, respectively. Using astrophysical models of black-hole populations calibrated to observations and black-hole perturbation theory calculations of the gravitational emission, we find that, in optimistic scenarios, LIGO could observe a stochastic background of gravitational radiation in the range $m_s\in [2\times 10^{-13}, 10^{-12}]$ eV, and up to $10^4$ resolvable events in a $4$-year search if $m_s\sim 3\times 10^{-13}\,{\rm eV}$. LISA could observe a stochastic background for boson masses in the range $m_s\in [5\times 10^{-19}, 5\times 10^{-16}]$, and up to $\sim 10^3$ resolvable events in a $4$-year search if $m_s\sim 10^{-17}\,{\rm eV}$. LISA could further measure spins for black-hole binaries with component masses in the range $[10^3, 10^7]~M_\odot$, which is not probed by traditional spin-measurement techniques. A statistical analysis of the spin distribution of these binaries could either rule out scalar fields in the mass range $\sim [4 \times 10^{-18}, 10^{-14}]$ eV, or measure $m_s$ with ten percent accuracy if light scalars in the mass range $\sim [10^{-17}, 10^{-13}]$ eV exist.

Ultralight scalar fields around spinning black holes can trigger superradiant instabilities, forming a long-lived bosonic condensate outside the horizon. We use numerical solutions of the perturbed field equations and astrophysical models of massive and stellar-mass black hole populations to compute, for the first time, the stochastic gravitational-wave background from these sources. In optimistic scenarios the background is observable by Advanced LIGO and LISA for field masses $m_s$ in the range $\sim [2\times 10^{-13}, 10^{-12}]\,{\rm eV}$ and $\sim 5\times[ 10^{-19}, 10^{-16}]\,{\rm eV}$, respectively, and it can affect the detectability of resolvable sources. Our estimates suggest that an analysis of the stochastic background limits from LIGO O1 might already be used to marginally exclude axions with mass $\sim 10^{-12.5}{\rm eV}$. Semicoherent searches with Advanced LIGO (LISA) should detect $\sim 15~(5)$ to $200~(40)$ resolvable sources for scalar field masses $3\times 10^{-13}$ ($10^{-17}$) eV. LISA measurements of massive BH spins could either rule out bosons in the range $\sim [10^{-18}, 2\times 10^{-13}]$ eV, or measure $m_s$ with ten percent accuracy in the range $\sim[10^{-17}, 10^{-13}]$ eV.

Binary pulsars are ideal to test the foundations of General Relativity, such as Lorentz symmetry, which requires that experiments produce the same results in all free-falling (i.e.inertial) frames. We here break this symmetry in the gravitational sector by specifying a preferred time direction, and thus a preferred frame, at each spacetime point. We then examine the consequences of this gravitational Lorentz symmetry breaking in the orbital evolution of binary pulsars, focusing on the dissipative effects. We find that Lorentz symmetry breaking modifies these effects, and thus the orbital dynamics, in two different ways. First, it generically causes the emission of dipolar radiation, which makes the orbital separation decrease faster than in General Relativity. Second, the quadrupole component of the emission is also modified. The orbital evolution depends critically on the sensitivities of the stars, which measure how their binding energies depend on the motion relative to the preferred frame. We calculate the sensitivities numerically and compute the predicted orbital decay rate of binary pulsars in Lorentz-violating gravity. By testing these predictions against observations, we place very stringent constraints on gravitational Lorentz violation.

Massive black-hole binaries, formed when galaxies merge, are among the primary sources of gravitational waves targeted by ongoing Pulsar Timing Array (PTA) experiments and the upcoming space-based LISA interferometer. However, their formation and merger rates are still highly uncertain. Recent upper limits on the stochastic gravitational-wave background obtained by PTAs are starting being in marginal tension with theoretical models for the pairing and orbital evolution of these systems. This tension can be resolved by assuming that these binaries are more eccentric or interact more strongly with the environment (gas and stars) than expected, or by accounting for possible selection biases in the construction of the theoretical models. However, another (pessimistic) possibility is that these binaries do not merge at all, but stall at large ($\sim$ pc) separations. We explore this extreme scenario by using a galaxy-formation semi-analytic model including massive black holes (isolated and in binaries), and show that future generations of PTAs will detect the stochastic gravitational-wave background from the massive black-hole binary population within $10-15$ years of observations, even in the "nightmare scenario" in which all binaries stall at the hardening radius. Moreover, we argue that this scenario is too pessimistic, because our model predicts the existence of a sub-population of binaries with small mass ratios ($q \lesssim 10^{-3}$) that should merge within a Hubble time simply as a result of gravitational-wave emission. This sub-population will be observable with large signal-to-noise ratios by future PTAs thanks to next-generation radio telescopes such as SKA or FAST, and possibly by LISA.

Recent work has confirmed that the masses of supermassive black holes, estimated from scaling relations with global properties such as the stellar masses of their host galaxies, may be biased high. Much of this may be caused by the requirement that the gravitational sphere of influence of the black hole must be resolved for the black-hole mass to be reliably estimated. We revisit this issue by using a comprehensive galaxy evolution semi-analytic model, which self-consistently evolves supermassive black holes from high-redshift seeds via gas accretion and mergers, and also includes AGN feedback. Once tuned to reproduce the (mean) correlation of black-hole mass with velocity dispersion, the model is unable to also account for the correlation with stellar mass. This behaviour is independent of the model's parameters, thus suggesting an internal inconsistency in the data. The predicted distributions, especially at the low-mass end, are also much broader than observed. However, if selection effects are included, the model's predictions tend to align with the observations. We also demonstrate that the correlations between the residuals of the local scaling relations are more effective than the scaling relations themselves at constraining AGN feedback models. In fact, we find that our semi-analytic model, while in apparent broad agreement with the scaling relations when accounting for selection biases, yields very weak correlations between their residuals at fixed stellar mass, in stark contrast with observations. This problem persists when changing the AGN feedback strength, and is also present in the $z\sim 0$ outputs of the hydrodynamic cosmological simulation Horizon-AGN, which includes state-of-the-art treatments of AGN feedback. This suggests that current AGN feedback models may be too weak or are simply not capturing the effect of the black hole on the stellar velocity dispersion.

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.

We investigate the capability of various configurations of the space interferometer eLISA to probe the late-time background expansion of the universe using gravitational wave standard sirens. We simulate catalogues of standard sirens composed by massive black hole binaries whose gravitational radiation is detectable by eLISA, and which are likely to produce an electromagnetic counterpart observable by future surveys. The main issue for the identification of a counterpart resides in the capability of obtaining an accurate enough sky localisation with eLISA. This seriously challenges the capability of four-link (2 arm) configurations to successfully constrain the cosmological parameters. Conversely, six-link (3 arm) configurations have the potential to provide a test of the expansion of the universe up to $z\sim 8$ which is complementary to other cosmological probes based on electromagnetic observations only. In particular, in the most favourable scenarios, they can provide a significant constraint on $H_0$ at the level of 0.5%. Furthermore, $(\Omega_M, \Omega_\Lambda)$ can be constrained to a level competitive with present SNIa results. On the other hand, the lack of massive black hole binary standard sirens at low redshift allows to constrain dark energy only at the level of few percent.

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 estimate the potential of present and future interferometric gravitational-wave detectors to test the Kerr nature of black holes through "gravitational spectroscopy," i.e. the measurement of multiple quasinormal mode frequencies from the remnant of a black hole merger. Using population synthesis models of the formation and evolution of stellar-mass black hole binaries, we find that Voyager-class interferometers will be necessary to perform these tests. Gravitational spectroscopy in the local Universe may become routine with the Einstein Telescope, but a 40-km facility like Cosmic Explorer is necessary to go beyond $z\sim 3$. In contrast, eLISA-like detectors should carry out a few - or even hundreds - of these tests every year, depending on uncertainties in massive black hole formation models. Many space-based spectroscopical measurements will occur at high redshift, testing the strong gravity dynamics of Kerr black holes in domains where cosmological corrections to general relativity (if they occur in nature) must be significant.

Massive black-hole binaries (MBHBs) are thought to be the main source of gravitational waves (GWs) in the low-frequency domain surveyed by ongoing and forthcoming Pulsar Timing Array campaigns and future space-borne missions, such as {\it eLISA}. However, many low-redshift MBHBs in realistic astrophysical environments may not reach separations small enough to allow significant GW emission, but rather stall on (sub)pc-scale orbits. This "last-parsec problem" can be eased by the appearance of a third massive black hole (MBH) -- the "intruder" -- whose action can force, under certain conditions, the inner MBHB on a very eccentric orbit, hence allowing intense GW emission eventually leading to coalescence. A detailed assessment of the process, ultimately driven by the induced Kozai-Lidov oscillations of the MBHB orbit, requires a general relativistic treatment and the inclusion of external factors, such as the Newtonian precession of the intruder orbit in the galactic potential and its hardening by scattering off background stars. In order to tackle this problem, we developed a three-body Post-Newtonian (PN) code framed in a realistic galactic potential, including both non-dissipative 1PN and 2PN terms, and dissipative terms such as 2.5PN effects, orbital hardening of the outer binary, and the effect of the dynamical friction on the early stages of the intruder dynamics. In this first paper of a series devoted at studing the dynamics of MBH triplets from a cosmological perspective, we describe, test and validate our code.

We revisit the problem of predicting the spin magnitude and direction of the black hole resulting from the merger of two black holes with arbitrary masses and spins inspiralling in quasi-circular orbits. We do this by analyzing a catalog of 619 recent numerical-relativity simulations collected from the literature and spanning a large variety of initial conditions. By combining information from the post-Newtonian approximation, the extreme mass-ratio limit and perturbative calculations, we improve our previously proposed phenomenological formulae for the final remnant spin. In contrast with alternative suggestions in the literature, and in analogy with our previous expressions, the new formula is a simple algebraic function of the initial system parameters and is not restricted to binaries with spins aligned/anti-aligned with the orbital angular momentum, but can be employed for fully generic binaries. The accuracy of the new expression is significantly improved, especially for almost extremal progenitor spins and for small mass ratios, yielding a root-mean-square error $\sigma\approx0.002$ for aligned/anti-aligned binaries and $\sigma\approx0.006$ for generic binaries. Our new formula is suitable for cosmological applications and can be employed robustly in the analysis of the gravitational waveforms from advanced interferometric detectors.

Gravitational waves (GWs) provide a revolutionary tool to investigate yet unobserved astrophysical objects. Especially the first stars, which are believed to be more massive than present-day stars, might be indirectly observable via the merger of their compact remnants. We develop a self-consistent, cosmologically representative, semi-analytical model to simulate the formation of the first stars. By extrapolating binary stellar-evolution models at 10% solar metallicity to metal-free stars, we track the individual systems until the coalescence of the compact remnants. We estimate the contribution of primordial stars to the merger rate density and to the detection rate of the Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO). Owing to their higher masses, the remnants of primordial stars produce strong GW signals, even if their contribution in number is relatively small. We find a probability of $\gtrsim1\%$ that the current detection GW150914 is of primordial origin. We estimate that aLIGO will detect roughly 1 primordial BH-BH merger per year for the final design sensitivity, although this rate depends sensitively on the primordial initial mass function (IMF). Turning this around, the detection of black hole mergers with a total binary mass of $\sim 300\,\mathrm{M}_\odot$ would enable us to constrain the primordial IMF.

The aLIGO detection of the black-hole binary GW150914 opened a new era for probing extreme gravity. Many gravity theories predict the emission of dipole gravitational radiation by binaries. This is excluded to high accuracy in binary pulsars, but entire classes of theories predict this effect predominantly (or only) in binaries involving black holes. Joint observations of GW150914-like systems by aLIGO and eLISA will improve bounds on dipole emission from black-hole binaries by six orders of magnitude relative to current constraints, provided that eLISA is not dramatically descoped.

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.

We study the post-Newtonian expansion of a class of Lorentz-violating gravity theories that reduce to khronometric theory (i.e. the infrared limit of Horava gravity) in high-acceleration regimes, and reproduce the phenomenology of the modified Newtonian dynamics (MOND) in the low-acceleration, non-relativistic limit. Like in khronometric theory, Lorentz symmetry is violated in these theories by introducing a dynamical scalar field (the "khronon") whose gradient is enforced to be timelike. As a result, hypersurfaces of constant khronon define a preferred foliation of the spacetime, and the khronon can be thought of as a physical absolute time. The MOND phenomenology arises as a result of the presence, in the action, of terms depending on the acceleration of the congruence orthogonal to the preferred foliation. We find that if the theory is forced to reduce exactly to General Relativity (rather than to khronometric theory) in the high-acceleration regime, the post-Newtonian expansion breaks down at low accelerations, and the theory becomes strongly coupled. Nevertheless, we identify a sizeable region of the parameter space where the post-Newtonian expansion remains perturbative for all accelerations, and the theory passes both solar-system and pulsar gravity tests, besides producing a MOND phenomenology for the rotation curves of galaxies. We illustrate this explicitly with a toy model of a system containing only baryonic matter but no Dark Matter.

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

Gravity theories beyond General Relativity typically predict dipolar gravitational emission by compact-star binaries. This emission is sourced by "sensitivity" parameters depending on the stellar compactness. We introduce a general formalism to calculate these parameters, and show that in shift-symmetric Horndeski theories stellar sensitivities and dipolar radiation vanish, provided that the binary's dynamics is perturbative (i.e. the post-Newtonian formalism is applicable) and cosmological-expansion effects can be neglected. This allows reproducing the binary-pulsar observed orbital decay.

Studying how nuclear star clusters (NSCs) form and how they are related to the growth of the central massive black holes (MBHs) and their host galaxies is fundamental for our understanding of the evolution of galaxies and the processes that have shaped their central structures. We present the results of a semi-analytical galaxy formation model that follows the evolution of dark matter halos along merger trees, as well as that of the baryonic components. This model allows us to study the evolution of NSCs in a cosmological context, by taking into account the growth of NSCs due to both dynamical friction-driven migration of stellar clusters and star formation triggered by infalling gas, while also accounting for dynamical heating from (binary) MBHs. We find that in-situ star formation contributes a significant fraction (up to ~80%) of the total mass of NSCs in our model. Both NSC growth through in-situ star formation and through star cluster migration are found to generate NSC -- host galaxy scaling correlations that are shallower than the same correlations for MBHs. We explore the role of galaxy mergers on the evolution of NSCs, and show that observational data on NSC -- host galaxy scaling relations provide evidence of partial erosion of NSCs by MBH binaries in luminous galaxies. We show that this observational feature is reproduced by our models, and we make predictions about the NSC and MBH occupation fraction in galaxies. We conclude by discussing several implications for theories of NSC formation.

Black holes in General Relativity are known as Kerr black holes and are characterized solely by two parameters, the mass $M$ and the spin $J$. All the higher multipole moments of the gravitational field are functions of these two parameters. For instance, the quadrupole moment is $Q=-J^2/M$, which implies that a measurement of $M$, $J$, and $Q$ for black hole candidates would allow one to test whether these objects are really black holes as described by General Relativity. While future gravitational-wave experiments will be able to test the Kerr nature of these objects with very high accuracy, in this paper we show that it is possible to put constraints on the quadrupole moment of stellar-mass black hole candidates by using presently available X-ray data of the thermal spectrum of their accretion disk.

A literature compilation of nuclear star cluster (NSC) masses is used to study the correlation between global and nuclear properties. A comparison of observational data to the predictions of semi-analytical galaxy formation models places constraints on the co-evolution of NSCs, massive black holes (MBHs) and host galaxies. Both data and theoretical predictions show an increased scatter in the NSC scaling correlations at high galaxy masses, and we show that this is due to the progressively more efficient ejection of stars from NSCs caused by MBH binaries in more massive stellar spheroids. Our results provide a natural explanation of why in nucleated galaxies hosting a MBH, the ratio (M_NSC+M_MBH)/M_bulge (with M_bulge the host spheroid's mass) shows significantly less scatter than M_NSC/M_bulge, and suggest that the formation of MBHs and NSCs are not mutually exclusive, as also supported by observations of co-existing systems. Both MBHs and NSCs represent generic products of galaxy formation, with NSCs being destroyed or modified by the merger evolution of their companion MBHs.

As we showed in previous work, the dynamics and gravitational emission of binary neutron star systems in certain scalar-tensor theories can differ significantly from that expected from General Relativity in the coalescing stage. In this work we examine whether the characteristics of the electromagnetic counterparts to these binaries -- driven by magnetosphere interactions prior to the merger event -- can provide an independent way to test gravity in the most strongly dynamical stages of binary mergers. We find that the electromagnetic flux emitted by binaries in these scalar-tensor theories can show deviations from the GR prediction in particular cases. These differences are quite subtle, thus requiring delicate measurements to differentiate between GR and the type of scalar-tensor theories considered in this work using electromagnetic observations alone. However, if coupled with a gravitational-wave detection, electromagnetic measurements might provide a way to increase the confidence with which GR will be confirmed (or ruled out) by gravitational observations.

The evolving Laser Interferometer Space Antenna (eLISA) will revolutionize our understanding of the formation and evolution of massive black holes along cosmic history by probing massive black hole binaries in the $10^3-10^7$ solar mass range out to redshift $z\gtrsim 10$. High signal-to-noise ratio detections of $\sim 10-100$ binary coalescences per year will allow accurate measurements of the parameters of individual binaries (such as their masses, spins and luminosity distance), and a deep understanding of the underlying cosmic massive black hole parent population. This wealth of unprecedented information can lead to breakthroughs in many areas of physics, including astrophysics, cosmology and fundamental physics. We review the current status of the field, recent progress and future challenges.

The upcoming detection of gravitational waves by terrestrial interferometers will usher in the era of gravitational-wave astronomy. This will be particularly true when space-based detectors will come of age and measure the mass and spin of massive black holes with exquisite precision and up to very high redshifts, thus allowing for better understanding of the symbiotic evolution of black holes with galaxies, and for high-precision tests of General Relativity in strong-field, highly dynamical regimes. Such ambitious goals require that astrophysical environmental pollution of gravitational-wave signals be constrained to negligible levels, so that neither detection nor estimation of the source parameters are significantly affected. Here, we consider the main sources for space-based detectors -- the inspiral, merger and ringdown of massive black-hole binaries and extreme mass-ratio inspirals -- and account for various effects on their gravitational waveforms, including electromagnetic fields, cosmological evolution, accretion disks, dark matter, "firewalls" and possible deviations from General Relativity. We discover that the black-hole quasinormal modes are sharply different in the presence of matter, but the ringdown signal observed by interferometers is typically unaffected. The effect of accretion disks and dark matter depends critically on their geometry and density profile, but is negligible for most sources, except for few special extreme mass-ratio inspirals. Electromagnetic fields and cosmological effects are always negligible. We finally explore the implications of our findings for proposed tests of General Relativity with gravitational waves, and conclude that environmental effects will not prevent the development of precision gravitational-wave astronomy.

We present the results of a semianalytical model that evolves the masses and spins of massive black holes together with the properties of their host galaxies along the cosmic history. As a consistency check, our model broadly reproduces a number of observations, e.g. the cosmic star formation history, the black hole mass and luminosity function and the galaxy mass function at low redshift, the black hole to bulge mass relation, and the morphological distribution at low redshift. For the first time in a semianalytical investigation, we relax the simplifying assumptions of perfect coherency or perfect isotropy of the gas fueling the black holes. The dynamics of gas is instead linked to the morphological properties of the host galaxies, resulting in different spin distributions for black holes hosted in different galaxy types. We compare our results with the observed sample of spin measurements obtained through broad K$\alpha$ iron line fitting. The observational data disfavor both accretion along a fixed direction and isotropic fueling. Conversely, when the properties of the accretion flow are anchored to the kinematics of the host galaxy, we obtain a good match between theoretical expectations and observations. A mixture of coherent accretion and phases of activity in which the gas dynamics is similar to that of the stars in bulges (i.e., with a significant velocity dispersion superimposed to a net rotation) best describes the data, adding further evidence in support to the coevolution of massive black holes and their hosts.

[abridged abstract] No, within a broad class of scenarios. With the advent of gravitational-wave (GW) astronomy, environmental effects on the GW signal will eventually have to be quantified. Here we present a wide survey of the corrections due to these effects in two situations of great interest for GW astronomy: the black hole (BH) ringdown emission and the inspiral of two compact objects. We take into account various effects such as: electric charges, magnetic fields, cosmological evolution, possible deviations from General Relativity, firewalls, and various forms of matter such as accretion disks and dark matter halos. Our analysis predicts the existence of resonances dictated by the external mass distribution, which dominate the very late-time behavior of merger/ringdown waveforms. The mode structure can drastically differ from the vacuum case, yet the BH response to external perturbations is unchanged at the time scales relevant for detectors. This is because although the vacuum Schwarzschild resonances are no longer quasinormal modes of the system, they still dominate the response at intermediate times. Our results strongly suggest that both parametrized and ringdown searches should use at least two-mode templates. Our analysis of compact binaries shows that environmental effects are typically negligible for most eLISA sources, with the exception of very few special extreme mass ratio inspirals. We show in particular that accretion and hydrodynamic drag generically dominate over self-force effects for geometrically thin disks, whereas they can be safely neglected for geometrically thick disk environments, which are the most relevant for eLISA. Finally, we discuss how our ignorance of the matter surrounding compact objects implies intrinsic limits on the ability to constrain strong-field deviations from General Relativity.

Binary pulsars are excellent laboratories to test the building blocks of Einstein's theory of General Relativity. One of these is Lorentz symmetry which states that physical phenomena appear the same for all inertially moving observers. We study the effect of violations of Lorentz symmetry in the orbital evolution of binary pulsars and find that it induces a much more rapid decay of the binary's orbital period due to the emission of dipolar radiation. The absence of such behavior in recent observations allows us to place the most stringent constraints on Lorentz violation in gravity, thus verifying one of the cornerstones of Einstein's theory much more accurately than any previous gravitational observation.

[Abridged] [...] We study the mass and spin evolution of massive black holes within a semianalytical galaxy-formation model that follows the evolution of dark-matter halos along merger trees, as well as that of the baryonic components (hot gas, stellar and gaseous bulges, and stellar and gaseous galactic disks). This allows us to study the mass and spin evolution of massive black holes in a self-consistent way, by taking into account the effect of the gas present in galactic nuclei both during the accretion phases and during mergers. Also, we present predictions, as a function of redshift, for the fraction of gas-rich black-hole mergers -- in which the spins prior to the merger are aligned due to the gravito-magnetic torques exerted by the circumbinary disk -- as opposed to gas-poor mergers, in which the orientation of the spins before the merger is roughly isotropic. These predictions may be tested by LISA or similar spaced-based gravitational-wave detectors such as eLISA/NGO or SGO.

We derive an analytic phenomenological expression that predicts the final mass of the black-hole remnant resulting from the merger of a generic binary system of black holes on quasi-circular orbits. Besides recovering the correct test-particle limit for extreme mass-ratio binaries, our formula reproduces well the results of all the numerical-relativity simulations published so far, both when applied at separations of a few gravitational radii, and when applied at separations of tens of thousands of gravitational radii. These validations make our formula a useful tool in a variety of contexts ranging from gravitational-wave physics to cosmology. As representative examples, we first illustrate how it can be used to decrease the phase error of the effective-one-body waveforms during the ringdown phase. Second, we show that, when combined with the recently computed self-force correction to the binding energy of nonspinning black-hole binaries, it provides an estimate of the energy emitted during the merger and ringdown. Finally, we use it to calculate the energy radiated in gravitational waves by massive black-hole binaries as a function of redshift, using different models for the seeds of the black-hole population.

We present a framework to study generic neutron-star binaries in scalar-tensor theories of gravity. Our formalism achieves this goal by suitably interfacing a post-Newtonian orbital evolution (described by a set of ordinary differential equations) with a set of non-linear algebraic equations, which provide a description of the scalar charge of each binary's component along the evolution in terms of isolated-star data. We validate this semi-analytical procedure by comparing its results to those of fully general-relativistic simulations, and use it to investigate the behavior of binary systems in large portions of the parameter space of scalar-tensor theories. This allows us to shed further light on the phenomenon of "dynamical scalarization", which we uncovered in [Barausse, Palenzuela, Ponce and Lehner, Phys. Rev. D 87, 081506(R) (2013)] and which takes place in tight binaries, even for stars that have exactly zero scalar charge in isolation. We also employ our formalism to study representative binary systems, obtain their gravitational-wave signals and discuss the extent to which deviations from General Relativity can be detected. The insights gained by this framework allow us to additionally show that eccentric binaries can undergo scalarization/de-scalarization phenomena.

Lorentz-symmetry and the notion of light cones play a central role in the definition of horizons and the existence of black holes. Current observations provide strong indications that astrophysical black holes do exist in Nature. Here we explore what happens to the notion of a black hole in gravity theories where local Lorentz symmetry is violated, and discuss the relevant astrophysical implications. Einstein-aether theory and Horava gravity are used as the theoretical background for addressing this question. We review earlier results about static, spherically symmetric black holes, which demonstrate that in Lorentz-violating theories there can be a new type of horizon and, hence, a new notion of black hole. We also present both known and new results on slowly rotating black holes in these theories, which provide insights on how generic these new horizons are. Finally, we discuss the differences between black holes in Lorentz-violating theories and in General Relativity, and assess to what extent they can be probed with present and future observations.

This work devises a formalism to obtain the equations of motion for a black hole-fluid configuration. Our approach is based on a Post-Newtonian expansion and adapted to scenarios where obtaining the relevant dynamics requires long time-scale evolutions. These systems are typically studied with Newtonian approaches, which have the advantage that larger time-steps can be employed than in full general-relativistic simulations, but have the downside that important physical effects are not accounted for. The formalism presented here provides a relatively straightforward way to incorporate those effects in existing implementations, up to 2.5PN order, with lower computational costs than fully relativistic simulations.

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.

Scalar-tensor theories of gravity are natural phenomenological alternatives to General Relativity, where the gravitational interaction is mediated by a scalar degree of freedom, besides the usual tensor gravitons. In regions of the parameter space of these theories where constraints from both solar system experiments and binary-pulsar observations are satisfied, we show that binaries of neutron stars present marked differences from General Relativity in both the late-inspiral and merger phases. In particular, phenomena related to the spontaneous scalarization of isolated neutron stars take place in the late stages of the evolution of binary systems, with important effects in the ensuing dynamics. We comment on the relevance of our results for the upcoming Advanced LIGO/Virgo detectors.

Using the main result of a companion paper, in which the binding energy of a circular-orbit non-spinning compact binary system is computed at leading-order beyond the test-particle approximation, the exact expression of the effective-one-body (EOB) metric component g^eff_tt is obtained through first order in the mass ratio. Combining these results with the recent gravitational self-force calculation of the periastron advance for circular orbits in the Schwarzschild geometry, the EOB metric component g^eff_rr is also determined at linear order in the mass ratio. These results assume that the mapping between the real and effective Hamiltonians at the second and third post-Newtonian (PN) orders holds at all PN orders. Our findings also confirm the advantage of resumming the PN dynamics around the test-particle limit if the goal is to obtain a flexible model that can smoothly connect the test-mass and equal-mass limits.

Using the effective-one-body (EOB) formalism and a time-domain Teukolsky code, we generate inspiral, merger, and ringdown waveforms in the small mass-ratio limit. We use EOB inspiral and plunge trajectories to build the Teukolsky equation source term, and compute full coalescence waveforms for a range of black hole spins. By comparing EOB waveforms that were recently developed for comparable mass binary black holes to these Teukolsky waveforms, we improve the EOB model for the (2,2), (2,1), (3,3), and (4,4) modes. Our results can be used to quickly and accurately extract useful information about merger waves for binaries with spin, and should be useful for improving analytic models of such binaries.

In the effective-one-body (EOB) approach the dynamics of two compact objects of masses m1 and m2 and spins S1 and S2 is mapped into the dynamics of one test particle of mass mu = m1 m2/(m1+m2) and spin S* moving in a deformed Kerr metric with mass M = m1+m2 and spin Skerr. In a previous paper we computed an EOB Hamiltonian for spinning black-hole binaries that (i) when expanded in post-Newtonian orders, reproduces the leading order spin-spin coupling and the leading and next-to-leading order spin-orbit couplings for any mass ratio, and (iii) reproduces all spin-orbit couplings in the test-particle limit. Here we extend this EOB Hamiltonian to include next-to-next-to-leading spin-orbit couplings for any mass ratio. We discuss two classes of EOB Hamiltonians that differ by the way the spin variables are mapped between the effective and real descriptions. We also investigate the main features of the dynamics when the motion is equatorial, such as the existence of the innermost stable circular orbit and of a peak in the orbital frequency during the plunge subsequent to the inspiral.

A classical thought-experiment to destroy black holes was envisaged by Wald in 1974: it consists of throwing particles with large angular momentum into an extremal black hole, checking whether their capture can over-spin the black hole past the extremal limit and create a naked singularity. Wald showed that in the test-particle limit, particles that would be otherwise capable of producing naked singularities are simply scattered. Recently Jacobson and Sotiriou showed that if one considers instead a black hole that is almost, but not exactly extremal, then in the absence of backreaction effects particle capture could indeed over-spin the spacetime above the Kerr limit. Here we analyze back-reaction effects and show that for some of the trajectories giving rise to naked singularities, radiative effects can be neglected. However, for these orbits the conservative self-force is important, and seems to have the right sign to prevent the formation of naked singularities.

The $5 - 20 M_\odot$ dark objects in X-ray binary systems and the $10^5 - 10^9 M_\odot$ dark objects in galactic nuclei are currently thought to be the Kerr black holes predicted by General Relativity. However, direct observational evidence for this identification is still elusive, and the only viable approach to confirm the Kerr black hole hypothesis is to explore and rule out any other possibility. Here we investigate the final stages of the accretion process onto generic compact objects. While for Kerr black holes and for more oblate bodies the accreting gas reaches the innermost stable circular orbit (ISCO) and plunges into the compact object, we find that for more prolate bodies several scenarios are possible, depending on the spacetime geometry. In particular, we find examples in which the gas reaches the ISCO, but then gets trapped between the ISCO and the compact object. In this situation, accretion onto the compact object is possible only if the gas loses additional angular momentum, forming torus-like structures inside the ISCO.

We study spherical black-hole solutions in Einstein-aether theory, a Lorentz-violating gravitational theory consisting of General Relativity with a dynamical unit timelike vector (the "aether") that defines a preferred timelike direction. These are also solutions to the infrared limit of Horava-Lifshitz gravity. We explore parameter values of the two theories where all presently known experimental constraints are satisfied, and find that spherical black-hole solutions of the type expected to form by gravitational collapse exist for all those parameters. Outside the metric horizon, the deviations away from the Schwarzschild metric are typically no more than a few percent for most of the explored parameter regions, which makes them difficult to observe with electromagnetic probes, but in principle within reach of future gravitational-wave detectors. Remarkably, we find that the solutions possess a universal horizon, not far inside the metric horizon, that traps waves of any speed relative to the aether. A notion of black hole thus persists in these theories, even in the presence of arbitrarily high propagation speeds.

Using a Legendre transformation, we compute the unconstrained Hamiltonian of a spinning test-particle in a curved spacetime at linear order in the particle spin. The equations of motion of this unconstrained Hamiltonian coincide with the Mathisson-Papapetrou-Pirani equations. We then use the formalism of Dirac brackets to derive the constrained Hamiltonian and the corresponding phase-space algebra in the Newton-Wigner spin supplementary condition (SSC), suitably generalized to curved spacetime, and find that the phase-space algebra (q,p,S) is canonical at linear order in the particle spin. We provide explicit expressions for this Hamiltonian in a spherically symmetric spacetime, both in isotropic and spherical coordinates, and in the Kerr spacetime in Boyer-Lindquist coordinates. Furthermore, we find that our Hamiltonian, when expanded in Post-Newtonian (PN) orders, agrees with the Arnowitt-Deser-Misner (ADM) canonical Hamiltonian computed in PN theory in the test-particle limit. Notably, we recover the known spin-orbit couplings through 2.5PN order and the spin-spin couplings of type S_Kerr S (and S_Kerr^2) through 3PN order, S_Kerr being the spin of the Kerr spacetime. Our method allows one to compute the PN Hamiltonian at any order, in the test-particle limit and at linear order in the particle spin. As an application we compute it at 3.5PN order.

Jacobson and Sotiriou showed that rotating black holes could be spun-up past the extremal limit by the capture of non-spinning test bodies, if one neglects radiative and self-force effects. This would represent a violation of the Cosmic Censorship Conjecture in four-dimensional, asymptotically flat spacetimes. We show that for some of the trajectories giving rise to naked singularities, radiative effects can be neglected. However, for these orbits the conservative self-force is important, and seems to have the right sign to prevent the formation of naked singularities.

Superspinars are ultracompact objects whose mass M and angular momentum J violate the Kerr bound (cJ/GM^2>1). Recent studies analyzed the observable consequences of gravitational lensing and accretion around superspinars in astrophysical scenarios. In this paper we investigate the dynamical stability of superspinars to gravitational perturbations, considering either purely reflecting or perfectly absorbing boundary conditions at the "surface" of the superspinar. We find that these objects are unstable independently of the boundary conditions, and that the instability is strongest for relatively small values of the spin. Also, we give a physical interpretation of the various instabilities that we find. Our results (together with the well-known fact that accretion tends to spin superspinars down) imply that superspinars are very unlikely astrophysical alternatives to black holes.

Building on a recent paper in which we computed the canonical Hamiltonian of a spinning test particle in curved spacetime, at linear order in the particle's spin, we work out an improved effective-one-body (EOB) Hamiltonian for spinning black-hole binaries. As in previous descriptions, we endow the effective particle not only with a mass m, but also with a spin S*. Thus, the effective particle interacts with the effective Kerr background (having spin S_Kerr) through a geodesic-type interaction and an additional spin-dependent interaction proportional to S*. When expanded in post-Newtonian (PN) orders, the EOB Hamiltonian reproduces the leading order spin-spin coupling and the spin-orbit coupling through 2.5PN order, for any mass-ratio. Also, it reproduces all spin-orbit couplings in the test-particle limit. Similarly to the test-particle limit case, when we restrict the EOB dynamics to spins aligned or antialigned with the orbital angular momentum, for which circular orbits exist, the EOB dynamics has several interesting features, such as the existence of an innermost stable circular orbit, a photon circular orbit, and a maximum in the orbital frequency during the plunge subsequent to the inspiral. These properties are crucial for reproducing the dynamics and gravitational-wave emission of spinning black-hole binaries, as calculated in numerical relativity simulations.

The prediction of the spin of the black hole resulting from the merger of a generic black-hole binary system is of great importance to study the cosmological evolution of supermassive black holes. Several attempts have been recently made to model the spin via simple expressions exploiting the results of numerical-relativity simulations. Here, I first review the derivation of a formula, proposed in Barausse & Rezzolla, Apj 704 L40, which accurately predicts the final spin magnitude and direction when applied to binaries with separations of hundred or thousands of gravitational radii. This makes my formula particularly suitable for cosmological merger-trees and N-body simulations, which provide the spins and angular momentum of the two black holes when their separation is of thousands of gravitational radii. More importantly, I investigate the physical reason behind the good agreement between my formula and numerical relativity simulations, and nail it down to the fact that my formula takes into account the post-Newtonian precession of the spins and angular momentum in a consistent manner.

We investigate the atomic (HI) and molecular (H_2) Hydrogen content of normal galaxies by combining observational studies linking galaxy stellar and gas budgets to their host dark matter (DM) properties, with a physically grounded galaxy formation model. This enables us to analyse empirical relationships between the virial, stellar, and gaseous masses of galaxies and explore their physical origins. Utilising a semi-analytic model (SAM) to study the evolution of baryonic material within evolving DM halos, we study the effects of baryonic infall and various star formation and feedback mechanisms on the properties of formed galaxies using the most up-to-date physical recipes. We find that in order to significantly improve agreement with observations of low-mass galaxies we must suppress the infall of baryonic material and exploit a two-phase interstellar medium (ISM), where the ratio of HI to H_2 is determined by the galactic disk structure. Modifying the standard Schmidt-Kennicutt star formation law, which acts upon the total cold gas in galaxy discs and includes a critical density threshold, and employing a star formation law which correlates with the H_2 gas mass results in a lower overall star formation rate. This in turn, allows us to simultaneously reproduce stellar, HI and H_2 mass functions of normal galaxies.

We use our model for the formation and evolution of galaxies within a two-phase galaxy formation scenario, showing that the high-redshift domain typically supports the growth of spheroidal systems, whereas at low redshifts the predominant baryonic growth mechanism is quiescent and may therefore support the growth of a disc structure. Under this framework we investigate the evolving galaxy population by comparing key observations at both low and high-redshifts, finding generally good agreement. By analysing the evolutionary properties of this model, we are able to recreate several features of the evolving galaxy population with redshift, naturally reproducing number counts of massive star-forming galaxies at high redshifts, along with the galaxy scaling relations, star formation rate density and evolution of the stellar mass function. Building upon these encouraging agreements, we make model predictions that can be tested by future observations. In particular, we present the expected evolution to z=2 of the super-massive black hole mass function, and we show that the gas fraction in galaxies should decrease with increasing redshift in a mass, with more and more evolution going to higher and higher masses. Also, the characteristic transition mass from disc to bulge dominated system should decrease with increasing redshift.

Knowledge of the spin of the black hole resulting from the merger of a generic black-hole binary is of great importance for studying the cosmological evolution of supermassive black holes. Several attempts have been made to model the spin via simple expressions exploiting the results of numerical-relativity simulations. While these expressions are in reasonable agreement with the simulations, they neglect the precession of the binary's orbital plane, and cannot therefore be applied directly -- i.e., without evolving the system to small separations using post-Newtonian theory -- to binaries with separations larger than a few hundred gravitational radii. While not a problem in principle, this may be impractical if the formulas are employed in cosmological merger-trees or N-body simulations, which provide the spins and angular momentum of the two black holes when their separation is of hundreds or thousands of gravitational radii. The formula that we propose is instead built on improved assumptions and gives, for any separation, a very accurate prediction both for the norm of the final spin and for its direction. By comparing with the numerical data, we also show that the final-spin direction is very accurately aligned with the binary's total angular momentum at large separation. Hence, observations of the final-spin direction (e.g. via a jet) can provide information on the binary's orbital plane at large separations and could be relevant, for instance, for studying X-shaped radio sources.

Although the Kerr solution is common to many gravity theories, its perturbations are different in different theories. Hence, perturbed Kerr black holes can probe deviations from General Relativity.

We provide a compact analytic formula to compute the spin of the black hole produced by the coalescence of two black holes. The expression, which uses an analytic fit of numerical-relativity data and relies on four assumptions, aims at modelling generic initial spin configurations and mass ratios. A comparison with numerical-relativity simulations already shows very accurate agreements with all of the numerical data available to date, but we also suggest a number of ways in which our predictions can be further improved.

We have studied extreme mass-ratio inspirals (EMRIs) in spacetimes containing a rotating black hole and a non self-gravitating torus with constant specific angular momentum. We have found that the effect of the hydrodynamic drag exerted by the torus on the satellite is much smaller than the corresponding one due to radiation reaction, for systems such as those generically expected in AGNs and at distances from the SMBH which can be probed with LISA. However, given the uncertainty on the parameters of these systems, there exist configurations in which the effect of the hydrodynamic drag can be comparable to the radiation-reaction one in phases of the inspiral which are detectable by LISA. This is the case, for instance, for a 10^6 M_sun SMBH surrounded by a corotating torus of comparable mass and with radius of 10^3-10^4 gravitational radii, or for a 10^5 M_sun SMBH surrounded by a corotating 10^4 M_sun torus with radius of 10^5 gravitational radii. Should these conditions be met in astrophysical systems, EMRI-gravitational waves could provide a characteristic signature of the presence of the torus. In fact, while radiation reaction always increases the inclination of the orbit with respect to the equatorial plane, the hydrodynamic drag from a torus corotating with the SMBH always decreases it. However, even when initially dominating over radiation reaction, the influence of the hydrodynamic drag decays very rapidly as the satellite moves into the very strong-field region of the SMBH (i.e., p <~ 5M). Although our results have been obtained for a specific class of tori, we argue that they will be qualitatively valid also for more generic distributions of the specific angular momentum.

In a previous paper we showed that static spherically symmetric objects which, in the vicinity of their surface, are well-described by a polytropic equation of state with 3/2<Gamma<2 exhibit a curvature singularity in Palatini f(R) gravity. We argued that this casts serious doubt on the validity of Palatini f(R) gravity as a viable alternative to General Relativity. In the present paper we further investigate this characteristic of Palatini f(R) gravity in order to clarify its physical interpretation and consequences.

Non-vacuum static spherically-symmetric solutions in Palatini f(R) gravity are examined. It is shown that for generic choices of f(R), there are commonly-used equations of state for which no satisfactory physical solution of the field equations can be found within this framework, apart from in the special case of General Relativity, casting doubt on whether Palatini f(R) gravity can be considered as giving viable alternatives to General Relativity.

The dynamical friction force experienced by a body moving at relativistic speed in a gaseous medium is examined. This force, which arises due to the gravitational interaction of the body with its own gravitationally-induced wake, is calculated for straight-line and circular motion, generalizing previous results by several authors. Possible applications to the study of extreme mass-ratio inspirals around strongly-accreting supermassive black holes are suggested.

We have performed a detailed analysis of orbital motion in the vicinity of a nearly extremal Kerr black hole. For very rapidly rotating black holes (spin a=J/M>0.9524M) we have found a class of very strong field eccentric orbits whose angular momentum L_z increases with the orbit's inclination with respect to the equatorial plane, while keeping latus rectum and eccentricity fixed. This behavior is in contrast with Newtonian intuition, and is in fact opposite to the "normal" behavior of black hole orbits. Such behavior was noted previously for circular orbits; since it only applies to orbits very close to the black hole, they were named "nearly horizon-skimming orbits". Our analysis generalizes this result, mapping out the full generic (inclined and eccentric) family of nearly horizon-skimming orbits. The earlier work on circular orbits reported that, under gravitational radiation emission, nearly horizon-skimming orbits tend to evolve to smaller orbit inclination, toward prograde equatorial configuration. Normal orbits, by contrast, always demonstrate slowly growing orbit inclination (orbits evolve toward the retrograde equatorial configuration). Using up-to-date Teukolsky-fluxes, we have concluded that the earlier result was incorrect: all circular orbits, including nearly horizon-skimming ones, exhibit growing orbit inclination. Using kludge fluxes based on a Post-Newtonian expansion corrected with fits to circular and to equatorial Teukolsky-fluxes, we argue that the inclination grows also for eccentric nearly horizon-skimming orbits. We also find that the inclination change is, in any case, very small. As such, we conclude that these orbits are not likely to have a clear and peculiar imprint on the gravitational waveforms expected to be measured by the space-based detector LISA.

The Parametrized Post-Newtonian expansion of gravitational theories with a scalar field coupled to the Gauss-Bonnet invariant is performed and confrontation of such theories with Solar system experiments is discussed.

To investigate the imprint on the gravitational-wave emission from extreme mass-ratio inspirals in non-pure Kerr spacetimes, we have studied the ``kludge'' waveforms generated in highly-accurate, numerically-generated spacetimes containing a black hole and a self-gravitating, homogeneous torus with comparable mass and spin. In order to maximize their impact on the produced waveforms, we have considered tori that are compact, massive and close to the central black hole, investigating under what conditions the LISA experiment could detect their presence. Our results show that for a large portion of the space of parameters the waveforms produced by EMRIs in these black hole-torus systems are indistinguishable from pure-Kerr waveforms. Hence, a ``confusion problem'' will be present for observations carried out over a timescale below or comparable to the dephasing time.

The luminosity distance-redshift relation is one of the fundamental tools of modern cosmology. We compute the luminosity distance-redshift relation in a perturbed flat matter-dominated Universe, taking into account the presence of cosmological inhomogeneities up to second order in perturbation theory. Cosmological observations implementing the luminosity distance-redshift relation tell us that the Universe is presently undergoing a phase of accelerated expansion. This seems to call for a mysterious Dark Energy component with negative pressure. Our findings suggest that the need of a Dark Energy fluid may be challenged once a realistic inhomogeneous Universe is considered and that an accelerated expansion may be consistent with a matter-dominated Universe.