This White Paper summarizes the discussions from a five-day workshop, involving 50 people from 18 countries, held in Cape Town, South Africa in February 2020. Convened by the International Astronomical Union's Executive Committee Working Group on Global Coordination of Ground and Space Astrophysics and sponsored by the Kavli Foundation, we discussed existing and potential bottlenecks for transient and multi-messenger astronomy, identifying eight broad areas of concern. Some of these are very similar to the challenges faced by many astronomers engaging in international collaboration, for example, data access policies, funding, theoretical and computational resources and workforce equity. Others, including, alerts, telescope coordination and target-of-opportunity implementation, are strongly linked to the time domain and are particularly challenging as we respond to transients. To address these bottlenecks we offer thirty-five specific recommendations, some of which are simply starting points and require development. These recommendations are not only aimed at collaborative groups and individuals, but also at the various organizations who are essential to making transient collaborations efficient and effective: including the International Astronomical Union, observatories, projects, scientific journals and funding agencies. We hope those involved in transient research will find them constructive and use them to develop collaborations with greater impact and more inclusive teams.

Massive black hole binaries (MBHBs) of $10^5 \, \rm M_\odot - 3 \times 10^7 \, \rm M_\odot $ merging in low redshift galaxies ($z\le4$) are sufficiently loud to be detected weeks before coalescence with the Laser Interferometer Space Antenna (LISA). This allows us to perform the parameter estimation $on$ $the$ $fly$, i.e. as a function of the time to coalescence during the inspiral phase, relevant for early warning of the planned LISA protected periods and for searches of electromagnetic signals. In this work, we study the evolution of the sky position, luminosity distance, chirp mass and mass ratio uncertainties as function of time left before merger. Overall, light systems with total intrinsic mass $\rm M_{\rm tot} = 3 \times 10^5 \, \rm M_\odot$ are characterized by smaller uncertainties than heavy ones ($\rm M_{\rm tot} = 10^7 \, \rm M_\odot$) during the inspiral. Luminosity distance, chirp mass and mass ratio are well constrained at the end of the inspiral. Concerning sky position, at $z=1$, MBHBs with $\rm M_{\rm tot} = 3 \times 10^5 \, \rm M_\odot$ can be localized with a median precision of $\simeq 10^2 \, \rm deg^2 (\simeq 1 \, \rm deg^2)$ at 1 month (1 hour) from merger, while the sky position of heavy MBHBs can be determined to $10 \, \rm deg^2$ only 1 hour before merger. However the uncertainty around the median values broadens with time, ranging in between 0.04 -- 20 $\rm deg^2$ (0.3 -- 3 $\times 10^3 \, \rm deg^2$) for light (heavy) systems at 1 hour before merger. At merger the sky localization improves down to $\simeq 10^{-1} \, \rm deg^2$ for all masses. For the benefit of the observer community, we provide the full set of data from our simulations and simple and ready-to-use analytical fits to describe the time evolution of uncertainties in the aforementioned parameters, valid for systems with total mass between $10^5$--$10^7 \, \rm M_\odot$ and redshift $0.3$--$3$.

We estimate the merger timescale of spectroscopically-selected, subparsec supermassive black hole binary (SMBHB) candidates by comparing their expected contribution to the gravitational wave background (GWB) with the sensitivity of current pulsar timing array (PTAs) experiments and in particular, with the latest upper limit placed by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). We find that the average timescale to coalescence of such SMBHBs is $\langle t_{\rm evol} \rangle > 6\times 10^4\,$yr, assuming that their orbital evolution in the PTA frequency band is driven by emission of gravitational waves. If some fraction of SMBHBs do not reside in spectroscopically detected active galaxies, and their incidence in active and inactive galaxies is similar, then the merger timescale could be $\sim 10$ times longer, $\langle t_{\rm evol} \rangle > 6\times 10^5\,$yr. These limits are consistent with the range of timescales predicted by theoretical models and imply that all the SMBHB candidates in our spectroscopic sample could be binaries without violating the observational constraints on the GWB. This result illustrates the power of the multi-messenger approach, facilitated by the pulsar timing arrays, in providing an independent statistical test of the nature of SMBHB candidates discovered in electromagnetic searches.

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.)

Within the next several years, pulsar-timing array programs will likely usher in the next era of gravitational-wave astronomy through the detection of a stochastic background of nanohertz-frequency gravitational waves, originating from a cosmological population of inspiraling supermassive binary black holes. While the source positions will likely be isotropic to a good approximation, the gravitational-wave angular power distribution will be anisotropic, with the most massive and/or nearby binaries producing signals that may resound above the background. We study such a realistic angular power distribution, developing fast and accurate sky-mapping strategies to localize pixels and extended regions of excess power while simultaneously modeling the background signal from the less massive and more distant ensemble. We find that power anisotropy will be challenging to discriminate from isotropy for realistic gravitational-wave skies, requiring SNR $>10$ in order to favor anisotropy with $10:1$ posterior odds in our case study. Amongst our techniques, modeling the population signal with multiple point sources in addition to an isotropic background provides the most physically-motivated and easily interpreted maps, while spherical-harmonic modeling of the square-root power distribution, $P(\hat\Omega)^{1/2}$, performs best in discriminating from overall isotropy. Our techniques are modular and easily incorporated into existing pulsar-timing array analysis pipelines.

We study the mass assembly and spin evolution of supermassive black holes (BHs) across cosmic time as well as the impact of gravitational recoil on the population of nuclear and wandering black holes (wBHs) by using the semi-analytical model LGalaxies run on top of Millennium merger trees. We track spin changes that BHs experience during both coalescence events and gas accretion phases. For the latter, we assume that spin changes are coupled with the bulge assembly. This assumption leads to predictions for the median spin values of z=0 BHs that depend on whether they are hosted by pseudobulges, classical bulges or ellipticals, being a ~ 0.9, 0.7 and 0.4, respectively. The outcomes of the model display a good consistency with z<4 quasar luminosity functions and the z=0 BH mass function, spin values and black hole-bulge mass correlation. Regarding the wBHs, we assume that they can originate from both the disruption of satellite galaxies (orphan wBH) and ejections due to gravitational recoils (ejected wBH). The model points to a number density of wBHs that increases with decreasing redshift, although this population is always ~2 dex smaller than the one of nuclear black holes. At all redshifts, wBHs are typically hosted in M_{halo}>10^{13} M_{\odot} and M_{stellar}>10^{10} M_{\odot}, being orphan wBHs the dominant type. Besides, independently of redshift and halo mass, \textit{ejected} wBHs inhabit the central regions (<0.3R_{200}) of the host DM halo, while orphan wBHs linger at larger scales (>0.5R_{200}). Finally, we find that gravitational recoils cause a progressive depletion of nuclear BHs with decreasing redshift and stellar mass. Moreover, ejection events lead to changes in the predicted local black hole-bulge relation, in particular for BHs in pseudobulges, for which the relation is flattened at M_{bulge}>10^{10.2} M{\odot} and the scatter increase up to ~3dex.

The gravitational-wave astronomical revolution began in 2015 with LIGO's observation of the coalescence of two stellar-mass black holes. Over the coming decades, ground-based detectors like LIGO will extend their reach, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will enable gravitational-wave observations of the massive black holes in galactic centres. Between LISA and ground-based observatories lies the unexplored decihertz gravitational-wave frequency band. Here, we propose a Decihertz Observatory to cover this band, and complement observations made by other gravitational-wave observatories. The decihertz band is uniquely suited to observation of intermediate-mass ($\sim 10^2$-$10^4 M_\odot$) black holes, which may form the missing link between stellar-mass and massive black holes, offering a unique opportunity to measure their properties. Decihertz observations will be able to detect stellar-mass binaries days to years before they merge and are observed by ground-based detectors, providing early warning of nearby binary neutron star mergers, and enabling measurements of the eccentricity of binary black holes, providing revealing insights into their formation. Observing decihertz gravitational-waves also opens the possibility of testing fundamental physics in a new laboratory, permitting unique tests of general relativity and the Standard Model of particle physics. Overall, a Decihertz Observatory will answer key questions about how black holes form and evolve across cosmic time, open new avenues for multimessenger astronomy, and advance our understanding of gravitation, particle physics and cosmology.

We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrapolation from higher mass ratios, we find that for $q\lesssim 10^{-3}$ the eccentricity growth rate becomes negative, i.e., the binary {\it circularises} as it shrinks. This behaviour is due to the subset of interacting stars captured in metastable counter-rotating orbits; those stars tend to extract angular momentum from the binary, promoting eccentricity growth for $q>10^{-3}$, but tend to inject angular momentum into the binary driving it towards circularisation for $q<10^{-3}$. The physical origin of this behaviour requires a detailed study of the orbits of this subset of stars and is currently under investigation. Our findings might have important consequences for intermediate MBHs (IMBHs) inspiralling onto MBHs (e.g. a putative $10^3\rm M_{\odot}$ black hole inspiralling onto SgrA$^*$).

The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field.

The evolution of cosmic structures, the formation and growth of the first black holes and the connection to their baryonic environment are key unsolved problems in astrophysics. The X-ray Athena mission and the gravitational-wave Laser Interferometer Space Antenna (LISA) offer independent and complementary angles on these problems. We show that up to 10 black hole binaries in the mass range 10^5 - 10^8 Msun discovered by LISA at redshift <~ 3.5 could be detected by Athena in an exposure time up to 100 ks, if prompt X-ray emission of ~ 1% - 10% of the Eddington luminosity is present. Likewise, if any LISA-detected extreme mass ratio inspirals occur in accretion disks, Athena can detect associated electromagnetic emission out to redshift ~ 1. Finally, warned by LISA, Athena can point in advance and stare at stellar-mass binary black hole mergers at redshift <~ 0.1. These science opportunities emphasise the vast discovery space of simultaneous observations from the two observatories, which would be missed if they were operated in different epochs.

We estimate the properties of the double neutron star (DNS) population that will be observable by the planned space-based interferometer LISA. By following the gravitational radiation driven evolution of DNSs generated from rapid population synthesis of massive binary stars, we estimate that around 35 DNSs will accumulate a signal-to-noise ratio above 8 over a four-year LISA mission. The observed population mainly comprises Galactic DNSs (94 per cent), but detections in the LMC (5 per cent) and SMC (1 per cent) may also be expected. The median orbital frequency of detected DNSs is expected to be 0.8 mHz, and many of them will be eccentric (median eccentricity of $0.11$). The orbital properties will provide insights into DNS progenitors and formation channels. LISA is expected to localise these DNSs to a typical angular resolution of $2^\circ$, with best-constrained sources localised to a few arcminutes. These localisations may allow neutron star natal kick magnitudes to be constrained through the Galactic distribution of DNSs, and make it possible to follow up the sources with radio pulsar searches. However, LISA is also expected to resolve $\sim 10^4$ Galactic double white dwarfs, many of which may have binary parameters that resemble DNSs; we discuss how the combined measurement of binary eccentricity, chirp mass, and sky location may aid the identification of a DNS. We expect the best-constrained DNSs to have eccentricities known to a few parts in a thousand, chirp masses measured to better than 1 per cent fractional uncertainty, and sky localisation at the level of a few arcminutes.

We determine the main properties of the Galactic binary black hole (BBH) population detectable by LISA and strategies to distinguish them from the much more numerous white dwarf binaries. We simulate BBH populations based on cosmological simulations of Milky Way-like galaxies and binary evolution models. We then determine their gravitational wave emission as observed by LISA and build mock catalogs. According to our model LISA will detect $\approx4(6)$ binary black holes assuming 4(10) years of operations. Those figures grow to $\approx6(9)$ when models are re-normalized to the inferred LIGO/Virgo merger rates. About 40\%(70\%) of the sources will have a good enough chirp mass measurement to separate them from the much lighter white dwarf and neutron star binaries. Most of the remaining sources should be identifiable by their lack of electromagnetic counterpart within $\approx100$ pc. These results are robust with respect to the current uncertainties of the BBH merger rate as measured by LIGO/Virgo as well as the global mass spectrum of the binaries. We determine there is a 94 per cent chance that LISA finds at least one of these systems, which will allow us to pinpoint the conditions where they were formed and possibly find unique electromagnetic signatures.

One of the central goals of LISA is the detection of gravitational waves from the merger of supermassive black holes. Contrary to stellar-mass black hole mergers, such events are expected to be rich X-ray sources due to the accretion of material from the circumbinary disks onto the black holes. The orbital dynamics before merger is also expected to modulate the resulting X-ray emission via Doppler boosting in a quasi-periodic way, and in a simple phase relation with the gravitational wave from the inspiral of the black holes. Detecting the X-ray source would enable a precise and early localization of the binary, thus allowing many telescopes to observe the very moment of the merger. Although identifying the correct X-ray source in the relatively large LISA sky localization will be challenging due to the presence of many confounding point sources, the quasi-periodic modulation may aid in the identification. We explore the practical feasibility of such idea. We simulate populations of merging supermassive black holes, their detection with LISA and their X-ray lightcurves using a simple model. Taking the parameters of the X-ray Telescope on the proposed NASA Transient Astrophysics Probe, we then design and simulate an observation campaign that searches for the modulated X-ray source while LISA is still observing the inspiral of the black holes. Assuming a fiducial LISA detection rate of $10$ mergers per year at redshift closer than $3.5$, we expect a few detections of modulated X-ray counterparts over the nominal duration of the LISA mission.

In current stellar evolutionary models, the occurrence of pair instability supernovae implies the lack of stellar black holes (BHs) with masses between about $[60, \, 120] \, \rm{M}_\odot$, resulting in the presence of an upper mass gap in the BH mass distribution. In this Letter, we develop a fiducial model describing BHs beyond the pair instability gap, by convolving the initial mass function and star formation rate with the metallicity evolution across cosmic time. Under the ansatz that the underlying physics of binary formation does not change beyond the gap, we then construct the cosmic population of merging BH binaries. The detection rate of BH binaries with both mass components above the gap, is found to range between $\simeq [0.4,\,7] \, \rm{yr}^{-1}$ for LIGO/Virgo at design sensitivity and $[10, \, 460] \, \rm{yr}^{-1}$ for third-generation ground based detectors, considering the most pessimistic and optimistic scenarios. LISA can detect individually these binaries up to thousands of years from coalescence. The number of events merging in less than four years, which enable multi-band observation in sequence, is expected in the range $[1, \, 20]$. While ET will detect all these events, LIGO/Virgo is expected to detect $\lesssim 50\%$ of them. Finally, we estimate that the gravitational-wave background from unresolved sources in the LISA band may be in principle detected with a signal-to-noise ratio between $ \simeq 2.5$ and $\simeq 80$.

We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves.

In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science.

The local black hole mass function (BHMF) is of great interest to a variety of astrophysical problems, ranging from black hole binary merger rates to an indirect census of the dominant seeding mechanism of supermassive black holes. In this Letter, we combine the latest galaxy stellar mass function from the Galaxy And Mass Assembly survey with X-ray-based constraints to the local black hole occupation fraction to probe the BHMF below $10^6$ $M_{\odot}$. Notwithstanding the large uncertainties inherent to the choice of a reliable observational proxy for black hole mass, the resulting range of BHMFs yields a combined normalization uncertainty of $\lesssim$1 dex over the $[10^5-10^6]$ $M_{\odot}$ range, where upcoming, space-based gravitational wave detectors are designed to be most sensitive.

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.

We present an analytic model to describe the supermassive black hole binary (SMBHB) merger rate in the Universe with astrophysical observables: galaxy stellar mass function, pair fraction, merger timescale and black hole - host galaxy relations. We construct observational priors and compute the allowed range of the characteristic spectrum $h_c$ of the gravitational wave background (GWB) to be $10^{-16}<h_c<10^{-15}$ at a frequency of $f=1/{\rm yr}$. We exploit our parametrization to tackle the problem of astrophysical inference from Pulsar Timing Array (PTA) observations. We simulate a series of upper limits and detections and use a nested sampling algorithm to explore the parameter space. Corroborating previous results, we find that the current PTA non-detection does not place significant constraints on any observables; however, either future upper limits or detections will significantly enhance our knowledge of the SMBHB population. If a GWB is not detected at a level of $h_c(f=1/{\rm yr})=10^{-17}$, our current understanding of galaxy and SMBHB mergers is disfavoured at a $5\sigma$ level, indicating a combination of severe binary stalling, over-estimating of the SMBH -- host galaxy relations, and extreme dynamical properties of merging SMBHBs. Conversely, future detections of a Square Kilometre Array (SKA)-type array will allow to constrain the normalization of the SMBHB merger rate in the Universe, the time between galaxy pairing and SMBHB merging, the normalization of the SMBH -- host galaxy relations and the dynamical binary properties, including their eccentricity and density of stellar environment.

Super-massive black hole binaries are expected to produce a GW signal in the nano-Hertz frequency band which may be detected by PTAs in the coming years. The signal is composed of both stochastic and individually resolvable components. Here we develop a generic Bayesian method for the analysis of resolvable sources based on the construction of `null-streams' which cancel the part of the signal held in common for each pulsar (the Earth-term). For an array of $N$ pulsars there are $N-2$ independent null-streams that cancel the GW signal from a particular sky location. This method is applied to the localisation of quasi-circular binaries undergoing adiabatic inspiral. We carry out a systematic investigation of the scaling of the localisation accuracy with signal strength and number of pulsars in the PTA. Additionally, we find that source sky localisation with the International PTA data release one is vastly superior than what is achieved by its constituent regional PTAs.

We compute the expected cosmic rates of tidal disruption events induced by individual massive black holes (MBHs) and by MBH binaries (MBHBs) - with a specific focus on the latter class - to explore the potential of TDEs to probe the cosmic population of sub-pc MBHBs. Rates are computed by combining MBH and MBHB population models derived from large cosmological simulations with estimates of the induced TDE rates for each class of objects. We construct empirical TDE spectra that fit a large number of observations in the optical, UV and X-ray and consider their observability by current and future survey instruments. Consistent with results in the literature, and depending on the detailed assumption of the model, we find that LSST and Gaia in optical and eROSITA in X-ray will observe a total of 3000-6000, 80-180 and 600-900 TDEs per year, respectively. Depending on the survey, one to several percent of these are prompted by MBHBs. In particular both LSST and eROSITA are expected to see 150-450 MBHB induced TDEs in their respective mission lifetimes, including 5-100 repeated flares. The latter provide an observational sample of binary candidates with relatively low contamination and have the potential of unveiling the sub-pc population of MBHBs in the mass range $10^5M_\odot<M<10^7M_\odot$, thus informing future low frequency gravitational wave observatories.

In the local universe, black holes of $10^{5-6}$ M$_{\odot}$ are hosted in galaxies displaying a variety of stellar profiles and morphologies. These black holes are the anticipated targets of LISA, the Laser Interferometer Space Antenna that will detect the low-frequency gravitational-wave signal emitted by binary black holes in this mass interval. In this paper, we infer upper limits on the lifetime of binary black holes of $10^{5-6}$ M$_{\odot}$ and up to $10^8$ M$_{\odot}$, forming in galaxy mergers, exploring two underlying stellar density profiles, by Dehnen and by Prugniel & Simien, and by exploiting local scaling relations between the mass of the black holes and several quantities of their hosts. We focus on the phase of the dynamical evolution when the binary is transitioning from the hardening phase ruled by the interaction with single stars to the phase driven by the emission of gravitational waves. We find that different stellar profiles predict very distinct trends with binary mass, with lifetimes ranging between fractions of a Gyr to more than 10 Gyr, and with a spread of about one order of magnitude, given by the uncertainties in the observed correlations, which are larger in the low-mass tail of the observed black hole population.

We combine sophisticated high precision scattering experiments, together with results from the Millenium-II simulation, to compute the cosmic merger rate of bound compact object (CO) binaries dynamically interacting with supermassive black hole binaries (SMBHBs). We consider binaries composed of white dwarfs (WDs), neutron stars (NSs) and black holes (BHs). The overall merger rates for WD-WD, NS-NS, BH-BH, BH-NS binaries and EBBH (eccentric binaries of black holes) from redshift $\sim 5$ are found to be $4.32\times 10^3\,{\rm yr}^{-1}$($5.93\times10^2\,{\rm yr}^{-1}$ for Type Ia SNe), $82.7\,{\rm yr}^{-1}$, $96.3\,{\rm yr}^{-1}$, $13.1\,{\rm yr}^{-1}$ and $148\,{\rm yr}^{-1}$ , respectively, for a nominal CO binary fraction in the Galactic centre of 0.1. We calculate the distance ($R$) distribution of the merger sites with respect to the host galaxies of the binaries. The distribution shows a wide range of distances up to $\sim$Mpc; this tail is produced by escaped hyper-velocity CO binaries. Due to the differences in the matter density of the surrounding environment, merger events with different $R$ are expected to display significantly different signatures in their EM counterparts. In particular, merger events (and especially NS-NS) producing a relativistic jet but occurring in the intergalactic medium will have very weak afterglow radiation relative to their prompt emission. These events, which we call 'off-center', can only be produced from a close encounter between CO binaries and SMBHBs; hence the detection of such merger events would indicate the existence of nearby SMBHBs, and in particular with high mass ratio, produced in the aftermath of a major galaxy merger.

In this paper we consider a scenario where the currently observed hypervelocity stars in our Galaxy have been ejected from the Galactic center as a result of dynamical interactions with an intermediate-mass black hole (IMBH) orbiting the central supermassive black hole (SMBH). By performing 3-body scattering experiments, we calculate the distribution of the ejected stars' velocities given various parameters of the IMBH-SMBH binary: IMBH mass, semimajor axis and eccentricity. We also calculate the rates of change of the BH binary orbital elements due to those stellar ejections. One of our new findings is that the ejection rate depends (although mildly) on the rotation of the stellar nucleus (its total angular momentum). We also compare the ejection velocity distribution with that produced by the Hills mechanism (stellar binary disruption) and find that the latter produces faster stars on average. Also, the IMBH mechanism produces an ejection velocity distribution which is flattened towards the BH binary plane while the Hills mechanism produces a spherically symmetric one. The results of this paper will allow us in the future to model the ejection of stars by an evolving BH binary and compare both models with \textit{Gaia} observations, for a wide variety of environments (galactic nuclei, globular clusters, the Large Magellanic Clouds, etc.).

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.

Future GW detector networks and EM observatories will provide a unique opportunity to observe the most luminous events in the Universe involving matter in extreme environs. They will address some of the key questions in physics and astronomy: formation and evolution of compact binaries, sites of formation of heavy elements and the physics of jets.

Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanying structure formation at early times, and the runaway collapse of early, ultra-dense stellar clusters. The origin of the early massive BHs remains an intriguing and long-standing unsolved puzzle in astrophysics. Here we discuss strategies for discerning between BH seeding models using electromagnetic observations. We argue that the most direct answers will be obtained through detection of BHs with masses M<10^5 Msun at redshifts z>10, where we expect them to first form. Reaching out to these redshifts and down to these masses is crucial, because BHs are expected to lose the memory of their initial assembly by the time they grow well above 10^5 Msun and are incorporated into higher-mass galaxies. The best way to detect 10^4-10^5 Msun BHs at high redshifts is by a sensitive X-ray survey. Critical constraining power is augmented by establishing the properties and the environments of their host galaxies in deep optical/IR imaging surveys. Required OIR data can be obtained with the JWST and WFIRST missions. The required X-ray flux limits (down to 10^{-19} erg/s/cm^2) are accessible only with a next-generation X-ray observatory which has both high (sub-1") angular resolution and high throughput. A combination of deep X-ray and OIR surveys will be capable of probing several generic markers of the BH seed scenarios, and resolving the long-stanging puzzle of their origin. These electromagnetic observations are also highly synergistic with the information from LISA on high-z BH mergers.

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.

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.

LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center.

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

The tidal disruption of stars by (super-)massive black holes in galactic nuclei has been discussed in theoretical terms for about 30 years but only in the past decade have we been able to detect such events in substantial numbers. Thus, we are now starting to carry out observational tests of models for the disruption. We are also formulating expectations for the inspiral and disruption of white dwarfs by intermediate-mass black holes with masses $< 10^5\;{\rm M}_\odot$. Such events are very rich with information and open a new window to intermediate-mass black holes, thought to live in dwarf galaxies and star clusters. They can inform us of the demographics of intermediate-mass black holes, stellar populations and dynamics in their immediate vicinity, and the physics of accretion of hydrogen-deficient material. The combination of upcoming transient surveys using ground-based, electromagnetic observatories and low-frequency gravitational wave observations is ideal for exploiting tidal disruptions of white dwarfs. The detection rate of gravitational wave signals, optimistically, may reach a few dozen per year in a volume up to $z\approx 0.1$. Gravitational wave observations are particularly useful because they yield the masses of the objects involved and allow determination of the spin of the black hole, affording tests of physical models for black hole formation and growth. They also give us advance warning of the electromagnetic flares by weeks or more. The right computing infrastructure for modern models for the disruption process and event rates will allow us to make the most of the upcoming observing facilities.

We propose a novel methodology to select host galaxy candidates of future pulsar timing array (PTA) detections of resolved gravitational waves (GWs) from massive black hole binaries (MBHBs). The method exploits the physical dependence of the GW amplitude on the MBHB chirp mass and distance to the observer, together with empirical MBH mass-host galaxy correlations, to rank potential host galaxies in the mass-redshift plane. This is coupled to a null-stream based likelihood evaluation of the GW amplitude and sky position in a Bayesian framework that assigns to each galaxy a probability of hosting the MBHB generating the GW signal. We test our algorithm on a set of realistic simulations coupling the likely properties of the first PTA resolved GW signal to synthetic all-sky galaxy maps. For a foreseeable PTA sky-localization precision of 100 squared degrees, we find that the GW source is hosted with 50% (90%) probability within a restricted number of <50 (<500) potential hosts. These figures are orders of magnitude smaller than the total number of galaxies within the PTA sky error-box, enabling extensive electromagnetic follow-up campaigns on a limited number of targets.

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.

Massive black-hole (MBH) binaries, which are expected to form following the merger of their parent galaxies, produce gravitational waves which will be detectable by Pulsar Timing Arrays at nanohertz frequencies (year periods). While no confirmed, compact MBH binary systems have been seen in electromagnetic observations, a large number of candidates have recently been identified in optical surveys of AGN variability. Using a combination of cosmological, hydrodynamic simulations; comprehensive, semi-analytic binary merger models; and analytic AGN spectra and variability prescriptions; we calculate the expected electromagnetic detection rates of MBH binaries as periodically variable AGN. In particular, we consider two independent variability models: (i) Doppler boosting due to large orbital velocities, and (ii) hydrodynamic variability in which the fueling of MBH accretion disks is periodically modulated by the companion. Our models predict that numerous MBH binaries should be present and distinguishable in the existing data. In particular, our fiducial models produce an expectation value of $0.2$ (Doppler) and $5$ (hydrodynamic) binaries to be identifiable in CRTS, while $20$ and $100$ are expected after five years of LSST observations. The brightness variations in most systems are too small to be distinguishable, but almost $1\%$ of AGN at redshifts $z \lesssim 0.6$ could be in massive binaries. We analyze the predicted binary parameters of observable systems and their selection biases, and include an extensive discussion of our model parameters and uncertainties.

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.

Pulsar-timing analyses are sensitive to errors in the solar-system ephemerides (SSEs) that timing models utilise to estimate the location of the solar-system barycentre, the quasi-inertial reference frame to which all recorded pulse times-of-arrival are referred. Any error in the SSE will affect all pulsars, therefore pulsar timing arrays (PTAs) are a suitable tool to search for such errors and impose independent constraints on relevant physical parameters. We employ the first data release of the International Pulsar Timing Array to constrain the masses of the planet-moons systems and to search for possible unmodelled objects (UMOs) in the solar system. We employ ten SSEs from two independent research groups, derive and compare mass constraints of planetary systems, and derive the first PTA mass constraints on asteroid-belt objects. Constraints on planetary-system masses have been improved by factors of up to 20 from the previous relevant study using the same assumptions, with the mass of the Jovian system measured at 9.5479189(3)$\times10^{-4}$ $M_{\odot}$. The mass of the dwarf planet Ceres is measured at 4.7(4)$\times10^{-10}$ $M_{\odot}$. We also present the first sensitivity curves using real data that place generic limits on the masses of UMOs, which can also be used as upper limits on the mass of putative exotic objects. For example, upper limits on dark-matter clumps are comparable to published limits using independent methods. While the constraints on planetary masses derived with all employed SSEs are consistent, we note and discuss differences in the associated timing residuals and UMO sensitivity curves.

Hypervelocity binaries (HVBs) can be produced by a close encounter between stellar binaries and supermassive black hole (SMBH)-intermediate massive black hole (IMBH) binaries. We perform scattering experiments between tight stellar binaries and an SMBH-IMBH binary (within the observationally allowed parameters), using high-precision N-body integrations. Simulation results show that only a few HVBs can be produced in the Galactic center if it harbors a $4\times10^3M_\odot$ IMBH. However, we find that tight HVBs can be efficiently produced in prograde coplanar orbits with respect to the SMBH-IMBH binary. Therefore, the trajectories of the HVBs provide a promising way to constrain the IMBH orbital plane, if one exists. More importantly, we find that the ejected tight binaries are eccentric. For other galaxies hosting equal-mass SMBH binaries for which the HVB ejection rate is much higher relative to SMBH-IMBH binaries, these ejected compact binaries constitute promising extragalactic merger sources that could be detected by gravitational wave detectors.

Blazars are a sub-population of quasars whose jets are nearly aligned with the line-of-sight, which tend to exhibit multi-wavelength variability on a variety of timescales. Quasi-periodic variability on year-like timescales has been detected in a number of bright sources, and has been connected to the orbital motion of a putative massive black hole binary. If this were indeed the case, those blazar binaries would contribute to the nanohertz gravitational-wave stochastic background. We test the binary hypothesis for the blazar population observed by the \textit{Fermi} Gamma-Ray Space Telescope, which consists of BL Lacertae objects and flat-spectrum radio quasars. Using mock populations informed by the luminosity functions for BL Lacertae objects and flat-spectrum radio quasars with redshifts $z \le 2$, we calculate the expected gravitational wave background and compare it to recent pulsar timing array upper limits. The two are consistent only if a fraction $\lesssim 10^{-3}$ of blazars hosts a binary with orbital periods $<5$ years. We therefore conclude that binarity cannot significantly explain year-like quasi-periodicity in blazars.

We describe some key astrophysical processes driving the formation and evolution of black hole binaries of different nature, from stellar-mass to supermassive systems. In the first part, we focus on the mainstream channels proposed for the formation of stellar mass binaries relevant to ground-based gravitational wave detectors, namely the {\it field} and the {\it dynamical} scenarios. For the field scenario, we highlight the relevant steps in the evolution of the binary, including mass transfer, supernovae explosions and kicks, common envelope and gravitational wave emission. For the dynamical scenario, we describe the main physical processes involved in the formation of star clusters and the segregation of black holes in their centres. We then identify the dynamical processes leading to binary formation, including three-body capture, exchanges and hardening. The second part of the notes is devoted to massive black hole formation and evolution, including the physics leading to mass accretion and binary formation. Throughout the notes, we provide several step-by-step pedagogical derivations, that should be particularly suited to undergraduates and PhD students, but also to gravitational wave physicists interested in approaching the subject of gravitational wave sources from an astrophysical perspective.

We review the reports of possible year-long quasi-periodicities of BL Lac objects in the $\gamma$-ray and optical bands, and present a homogeneous time analysis of the light curves of PKS2155$-$304, PG1553+113, and BL Lac. Based on results from a survey covering the entire Fermi $\gamma$-ray sky we have estimated the fraction of possible quasi-periodic BL Lac objects. We compared the cyclical behaviour in BL Lac objects with that derived from the search of possible optical periodicities in quasars, and find that at z$\lesssim$1 the cosmic density of quasi-periodic BL Lac objects is larger than that of quasi-periodic quasars. If the BL Lac quasi-periodicities were due to a supermassive binary black hole (SBBH) scenario, there could be a tension with the upper limits on the gravitational wave background measured by the pulsar timing array. The argument clearly indicates the difficulties of generally associating quasi-periodicities of BL Lac objects with SBBHs.

In currently favoured hierarchical cosmologies, the formation of massive black hole binaries (MBHBs) following galaxy mergers is unavoidable. Still, due the complex physics governing the (hydro)dynamics of the post-merger dense environment of stars and gas in galactic nuclei, the final fate of those MBHBs is still unclear. In gas-rich environments, it is plausible that turbulence and gravitational instabilities feed gas to the nucleus in the form of a series of cold incoherent clumps, thus providing a way to exchange energy and angular momentum between the MBHB and its surroundings. Within this context, we present a suite of smoothed-particle-hydrodynamical models to study the evolution of a sequence of near-radial turbulent gas clouds as they infall towards equal-mass, circular MBHBs. We focus on the dynamical response of the binary orbit to different levels of anisotropy of the incoherent accretion events. Compared to a model extrapolated from a set of individual cloud-MBHB interactions, we find that accretion increases considerably and the binary evolution is faster. This occurs because the continuous infall of clouds drags inwards circumbinary gas left behind by previous accretion events, thus promoting a more effective exchange of angular momentum between the MBHB and the gas. These results suggest that sub-parsec MBHBs efficiently evolve towards coalescence during the interaction with a sequence of individual gas pockets.

We search for continuous gravitational waves (CGWs) produced by individual super-massive black-hole binaries (SMBHBs) in circular orbits using high-cadence timing observations of PSR J1713$+$0747. We observe this millisecond pulsar using the telescopes in the European Pulsar Timing Array (EPTA) with an average cadence of approximately 1.6 days over the period between April 2011 and July 2015, including an approximately daily average between February 2013 and April 2014. The high-cadence observations are used to improve the pulsar timing sensitivity across the GW frequency range of $0.008-5$ $\mu$Hz. We use two algorithms in the analysis, including a spectral fitting method and a Bayesian approach. For an independent comparison, we also use a previously published Bayesian algorithm. We find that the Bayesian approaches provide optimal results and the timing observations of the pulsar place a 95 per cent upper limit on the sky-averaged strain amplitude of CGWs to be $\lesssim3.5\times10^{-13}$ at a reference frequency of 1 $\mu$Hz. We also find a 95 per cent upper limit on the sky-averaged strain amplitude of low-frequency CGWs to be $\lesssim1.4\times10^{-14}$ at a reference frequency of 20~nHz.

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.

Hypervelocity stars (HVSs) are amongst the fastest objects in our Milky Way. These stars are predicted to come from the Galactic center (GC) and travel along unbound orbits across the Galaxy. In the coming years, the ESA satellite Gaia will provide the most complete and accurate catalogue of the Milky Way, with full astrometric parameters for more than $1$ billion stars. In this paper, we present the expected sample size and properties (mass, magnitude, spatial, velocity distributions) of HVSs in the Gaia stellar catalogue. We build three Gaia mock catalogues of HVSs anchored to current observations, exploring different ejection mechanisms and GC stellar population properties. In all cases, we predict hundreds to thousands of HVSs with precise proper motion measurements within a few tens of kpc from us. For stars with a relative error in total proper motion below $10 \%$, the mass range extends to ~$10 M_{\odot}$ but peaks at ~$1$ $M_\odot$. The majority of Gaia HVSs will therefore probe a different mass and distance range compared to the current non-Gaia sample. In addition, a subset of a few hundreds to a few thousands of HVSs with $M$ ~ $3$ $M_\odot$ will be bright enough to have a precise measurement of the three-dimensional velocity from Gaia alone. Finally, we show that Gaia will provide more precise proper motion measurements for the current sample of HVS candidates. This will help identifying their birthplace narrowing down their ejection location, and confirming or rejecting their nature as HVSs. Overall, our forecasts are extremely encouraging in terms of quantity and quality of HVS data that can be exploited to constrain both the Milky Way potential and the GC properties.

Massive black hole binaries (MBHBs) represent an unavoidable outcome of hierarchical galaxy formation, but their dynamical evolution at sub-parsec scales is poorly understood, due to a combination of uncertainties in theoretical models and lack of firm observational evidence. In gas rich environments, it has been shown that a putative extended, steady circumbinary gaseous disc plays an important role in the MBHB evolution, facilitating its coalescence. How gas on galactic scales is transported to the nuclear region to form and maintain such a stable structure is, however, unclear. If, following a galaxy merger, turbulent gas is condenses in cold clumps and filaments that are randomly scattered, gas is naturally transported on parsec scales and interacts with the MBHB in discrete incoherent pockets. The aim of this work is to investigate the gaseous structures arising from this interaction. We employ a suite of smoothed-particle-hydrodynamic simulations to study the formation and evolution of gaseous structures around a MBHB constantly perturbed by the incoherent infall of molecular clouds. We investigate the influence of the infall rate and angular momentum distribution of the clouds on the geometry and stability of the arising structures. We find that the continuous supply of incoherent clouds is a double-edge sword, resulting in the intermittent formation and disruption of circumbinary structures. Anisotropic cloud distributions featuring an excess of co-rotating events generate more prominent co-rotating circumbinary discs. Similar structures are seen when mostly counter-rotating clouds are fed to the binary, even though they are more compact and less stable. In general, our simulations do not show the formation of extended smooth and stable circumbinary discs, typically assumed in analytical and numerical investigations of the the long term evolution of MBHBs. (Abridged)

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.

The advent of time domain astronomy is revolutionizing our understanding of the Universe. Programs such as the Catalina Real-time Transient Survey (CRTS) or the Palomar Transient Factory (PTF) surveyed millions of objects for several years, allowing variability studies on large statistical samples. The inspection of $\approx$250k quasars in CRTS resulted in a catalogue of 111 potentially periodic sources, put forward as supermassive black hole binary (SMBHB) candidates. A similar investigation on PTF data yielded 33 candidates from a sample of $\approx$35k quasars. Working under the SMBHB hypothesis, we compute the implied SMBHB merger rate and we use it to construct the expected gravitational wave background (GWB) at nano-Hz frequencies, probed by pulsar timing arrays (PTAs). After correcting for incompleteness and assuming virial mass estimates, we find that the GWB implied by the CRTS sample exceeds the current most stringent PTA upper limits by almost an order of magnitude. After further correcting for the implicit bias in virial mass measurements, the implied GWB drops significantly but is still in tension with the most stringent PTA upper limits. Similar results hold for the PTF sample. Bayesian model selection shows that the null hypothesis (whereby the candidates are false positives) is preferred over the binary hypothesis at about $2.3\sigma$ and $3.6\sigma$ for the CRTS and PTF samples respectively. Although not decisive, our analysis highlights the potential of PTAs as astrophysical probes of individual SMBHB candidates and indicates that the CRTS and PTF samples are likely contaminated by several false positives.

We calculate the properties, occurrence rates and detection prospects of individually resolvable 'single sources' in the low frequency gravitational wave (GW) spectrum. Our simulations use the population of galaxies and massive black hole binaries from the Illustris cosmological hydrodynamic simulations, coupled to comprehensive semi-analytic models of the binary merger process. Using mock pulsar timing arrays (PTA) with, for the first time, varying red-noise models, we calculate plausible detection prospects for GW single sources and the stochastic GW background (GWB). Contrary to previous results, we find that single sources are at least as detectable as the GW background. Using mock PTA, we find that these 'foreground' sources (also 'deterministic'/'continuous') are likely to be detected with $\sim 20 \, \textrm{yr}$ total observing baselines. Detection prospects, and indeed the overall properties of single sources, are only moderately sensitive to binary evolution parameters---namely eccentricity & environmental coupling, which can lead to differences of $\sim 5 \,\textrm{yr}$ in times to detection. Red noise has a stronger effect, roughly doubling the time to detection of the foreground between a white-noise only model ($\sim 10$ - $15 \, \textrm{yr}$) and severe red noise ($\sim 20$ - $30 \, \textrm{yr}$). The effect of red noise on the GWB is even stronger, suggesting that single source detections may be more robust. We find that typical signal-to-noise ratios for the foreground peak near $f = 0.1 \, \textrm{yr}^{-1}$, and are much less sensitive to the continued addition of new pulsars to PTA.

We discuss the astrophysical science case for a decihertz gravitational-wave mission. We focus on unique opportunities for scientific discovery in this frequency range, including probes of type IA supernova progenitors, mergers in the presence of third bodies, intermediate mass black holes, seeds of massive black holes, improved sky localization, and tracking the population of merging compact binaries.

Pulsars are the most stable macroscopic clocks found in nature. Spinning with periods as short as a few milliseconds, their stability can supersede that of the best atomic clocks on Earth over timescales of a few years. Stable clocks are synonymous with precise measurements, which is why pulsars play a role of paramount importance in testing fundamental physics. As a pulsar rotates, the radio beam emitted along its magnetic axis appears to us as pulses because of the lighthouse effect. Thanks to the extreme regularity of the emitted pulses, minuscule disturbances leave particular fingerprints in the times-of-arrival (TOAs) measured on Earth with the technique of pulsar timing. Tiny deviations from the expected TOAs, predicted according to a theoretical timing model based on known physics, can therefore reveal a plethora of interesting new physical effects. Pulsar timing can be used to measure the dynamics of pulsars in compact binaries, thus probing the post-Newtonian expansion of general relativity beyond the weak field regime, while offering unique possibilities of constraining alternative theories of gravity. Additionally, the correlation of TOAs from an ensemble of millisecond pulsars can be exploited to detect low-frequency gravitational waves of astrophysical and cosmological origins. We present a comprehensive review of the many applications of pulsar timing as a probe of gravity, describing in detail the general principles, current applications and results, as well as future prospects.

Supermassive black hole binaries (SMBHBs) in the 10 million to 10 billion $M_\odot$ range form in galaxy mergers, and live in galactic nuclei with large and poorly constrained concentrations of gas and stars. There are currently no observations of merging SMBHBs--- it is in fact possible that they stall at their final parsec of separation and never merge. While LIGO has detected high frequency GWs, SMBHBs emit GWs in the nanohertz to millihertz band. This is inaccessible to ground-based interferometers, but possible with Pulsar Timing Arrays (PTAs). Using data from local galaxies in the 2 Micron All-Sky Survey, together with galaxy merger rates from Illustris, we find that there are on average $91\pm7$ sources emitting GWs in the PTA band, and $7\pm2$ binaries which will never merge. Local unresolved SMBHBs can contribute to GW background anisotropy at a level of $\sim20\%$, and if the GW background can be successfully isolated, GWs from at least one local SMBHB can be detected in 10 years.

The formation of massive black hole binaries (MBHBs) is an unavoidable outcome of galaxy evolution via successive mergers. However, the mechanism that drives their orbital evolution from parsec separations down to the gravitational wave (GW) dominated regime is poorly understood, and their final fate is still unclear. If such binaries are embedded in gas-rich and turbulent environments, as observed in remnants of galaxy mergers, the interaction with gas clumps (such as molecular clouds) may efficiently drive their orbital evolution. Using numerical simulations, we test this hypothesis by studying the dynamical evolution of an equal-mass, circular MBHB accreting infalling molecular clouds. We investigate different orbital configurations, modelling a total of 13 systems to explore different possible impact parameters and relative inclinations of the cloud-binary encounter. We focus our study on the prompt, transient phase during the first few orbits when the dynamical evolution of the binary is fastest, finding that this evolution is dominated by the exchange of angular momentum through gas capture by the individual black holes and accretion. Building on these results, we construct a simple model for evolving a MBHB interacting with a sequence of clouds, which are randomly drawn from reasonable populations with different levels of anisotropy in their angular momenta distributions. We show that the binary efficiently evolves down to the GW emission regime within a few hundred million years, overcoming the 'final parsec' problem regardless of the stellar distribution.

Pulsar timing arrays (PTAs) are presently the only means to search for the gravitational wave stochastic background from supermassive black hole binary populations, considered to be within the grasp of current or near future observations. However, the stringent upperlimit set by the Parkes PTA (Shannon et al. 2013, 2015) has been interpreted as excluding at $> 90\%$ confidence the current paradigm of binary assembly through galaxy mergers and hardening via stellar interactions, suggesting evolution is accelerated (by stars and/or gas) or stalled. Using Bayesian hierarchical modelling, we consider implications of this upperlimit for a comprehensive range of astrophysical scenarios, without invoking stalling nor more exotic physical processes. We find they are fully consistent with the upperlimit, but (weak) bounds on population parameters can be inferred. Bayes factors between models vary between $\approx 1.03$ -- $5.81$ and Kullback-Leibler divergences between characteristic amplitude prior and posterior lie between $0.37$ -- $0.85$. Considering prior astrophysical information on galaxy merger rates, recent upwards revisions of the black hole-galaxy bulge mass relation (Kormendy & Ho 2013) are disfavoured at $1.6\sigma$ against lighter models (e.g. Shankar et al. 2016). We also show, if no detection is achieved once sensitivity improves by an order of magnitude, the most optimistic scenario is disfavoured at $3.9\sigma$.

We present a fast and versatile method to calculate the characteristic spectrum $h_c$ of the gravitational wave background (GWB) emitted by a population of eccentric massive black hole binaries (MBHBs). We fit the spectrum of a reference MBHB with a simple analytic function and show that the spectrum of any other MBHB can be derived from this reference spectrum via simple scalings of mass, redshift and frequency. We then apply our calculation to a realistic population of MBHBs evolving via 3-body scattering of stars in galactic nuclei. We demonstrate that our analytic prescription satisfactorily describes the signal in the frequency band relevant to pulsar timing array (PTA) observations. Finally we model the high frequency steepening of the GWB to provide a complete description of the features characterizing the spectrum. For typical stellar distributions observed in massive galaxies, our calculation shows that 3-body scattering alone is unlikely to affect the GWB in the PTA band and a low frequency turnover in the spectrum is caused primarily by high eccentricities.

We consider the inverse problem in pulsar timing array (PTA) analysis, investigating what astrophysical information about the underlying massive black hole binary (MBHB) population can be recovered from the detection of a stochastic gravitational wave background (GWB). We employ a physically motivated model that connects the GWB spectrum to a series of parameters describing the underlying redshift evolution of the MBHB mass function and to the typical eccentricity they acquire while interacting with the dense environment of post merger galactic nuclei. This allows the folding in of information about the spectral shape of the GWB into the analysis. The priors on the model parameters are assumed to be uninformative and consistent with the current lack of secure observations of sub-parsec MBHBs. We explore the implications of current upper limits as well as of future detections with a variety of PTA configurations. We confirm our previous finding that current upper limits can only place an upper bound on the overall MBHB merger rate. Depending on the properties of the array, future detections can also constrain several MBHB population models at different degrees of fidelity. In particular, a simultaneous detection of a steepening of the spectrum at high frequency and a bending at low frequency will place strong constraints on both the MBHB mass function and on the typical eccentricity of inspiralling MBHBs, providing insights on MBHB astrophysics unlikely to be achievable by any other means.

Pulsars are natural cosmic clocks. On long timescales they rival the precision of terrestrial atomic clocks. Using a technique called pulsar timing, the exact measurement of pulse arrival times allows a number of applications, ranging from testing theories of gravity to detecting gravitational waves. Also an external reference system suitable for autonomous space navigation can be defined by pulsars, using them as natural navigation beacons, not unlike the use of GPS satellites for navigation on Earth. By comparing pulse arrival times measured on-board a spacecraft with predicted pulse arrivals at a reference location (e.g. the solar system barycenter), the spacecraft position can be determined autonomously and with high accuracy everywhere in the solar system and beyond. We describe the unique properties of pulsars that suggest that such a navigation system will certainly have its application in future astronautics. We also describe the on-going experiments to use the clock-like nature of pulsars to "construct" a galactic-sized gravitational wave detector for low-frequency (f_GW ~1E-9 - 1E-7 Hz) gravitational waves. We present the current status and provide an outlook for the future.

Pulsar Timing Arrays (PTA) around the world are using the incredible consistency of millisecond pulsars to measure low frequency gravitational waves from (super)Massive Black Hole (MBH) binaries. We use comprehensive MBH merger models based on cosmological hydrodynamic simulations to predict the spectrum of the stochastic Gravitational-Wave Background (GWB). We use real Time-of-Arrival (TOA) specifications from the European, NANOGrav, Parkes, and International PTA (IPTA) to calculate realistic times to detection of the GWB across a wide range of model parameters. In addition to exploring the parameter space of environmental hardening processes (in particular: stellar scattering efficiencies), we have expanded our models to include eccentric binary evolution which can have a strong effect on the GWB spectrum. Our models show that strong stellar scattering and high characteristic eccentricities enhance the GWB strain amplitude near the PTA sensitive "sweet-spot" (near the frequency $f = 1 \, \mathrm{yr}^{-1}$), slightly improving detection prospects in these cases. While the GWB $amplitude$ is degenerate between cosmological and environmental parameters, the location of a spectral turnover at low frequencies ($f \lesssim 0.1 \, \mathrm{yr}^{-1}$) is strongly indicative of environmental coupling. At high frequencies ($f\gtrsim 1 \, \mathrm{yr}^{-1}$), the GWB spectral index can be used to infer the number density of sources and possibly their eccentricity distribution. Even with merger models that use pessimistic environmental and eccentricity parameters, if the current rate of PTA expansion continues, we find that the International PTA is highly likely to make a detection within about 10 years.

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.

The recent Advanced LIGO detections of coalescing black hole binaries (BHBs) imply a large population of such systems emitting at milli-Hz frequencies, accessible to the Laser Interferometer Space Antenna (LISA). We show that these systems provide a new class of cosmological standard sirens. Direct LISA luminosity distance -$D_l$- measurements, combined with the inhomogeneous redshift -$z$- distribution of possible host galaxies provide an effective way to populate the $D_l-z$ diagram at $z<0.1$, thus allowing a precise local measurement of the Hubble expansion rate. To be effective, the method requires a sufficiently precise LISA distance determination and sky localization of a sizeable number of BHBs, which is best achieved for a 6-link detector configuration. We find that, for a BHB population consistent with current fiducial LIGO rates, the Hubble constant $H_0$ can be determined at the $\sim$5% and $\sim$2% level (68% confidence) assuming two and five million Km arm-length respectively.

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.

Soon after the observation of the first black hole binary (BHB) by advanced LIGO (aLIGO), GW150914, it was realised that such a massive system would have been observable in the milli-Hz (mHz) band few years prior to coalescence. Operating in the frequency range 0.1-100 mHz, the Laser Interferometer Space Antenna (LISA) can potentially detect up to thousands inspiralling BHBs, based on the coalescence rates inferred from the aLIGO first observing run (O1). The vast majority of them (those emitting at $f<10$ mHz) will experience only a minor frequency drift during LISA lifetime, resulting in signals similar to those emitted by galactic white dwarf binaries. At $f>10$ mHz however, several of them will sweep through the LISA band, eventually producing loud coalescences in the audio-band probed by aLIGO. This contribution reviews the scientific potential of these new class of LISA sources which, in the past few months, has been investigated in several contexts, including multi-messenger and multi-band gravitational wave astronomy, BHB astrophysics, tests of alternative theories of gravity and cosmography.

We present two methods for determining the significance of a stochastic gravitational-wave (GW) background affecting a pulsar-timing array, where detection is based on evidence for quadrupolar spatial correlations between pulsars. Rather than constructing noise simulations, we eliminate the GWB spatial correlations in the true datasets to assess detection significance with all real data features intact. In our first method, we perform random phase shifts in the signal-model basis functions. This phase shifting eliminates signal phase coherence between pulsars, while keeping the statistical properties of the pulsar timing residuals intact. We then explore a method to null correlations between pulsars by using a "scrambled" overlap-reduction function in the signal model for the array. This scrambled function is orthogonal to what we expect of a real GW background signal. We demonstrate the efficacy of these methods using Bayesian model selection on a set of simulated datasets that contain a stochastic GW signal, timing noise, undiagnosed glitches, and uncertainties in the Solar system ephemeris. Finally, we introduce an overarching formalism under which these two techniques are naturally linked. These methods are immediately applicable to all current pulsar-timing array datasets, and should become standard tools for future analyses.

The focus of this Chapter is on describing the prospective sources of the gravitational wave universe accessible to present and future observations, from kHz, to mHz down to nano-Hz frequencies. The multi-frequency gravitational wave universe gives a deep view into the cosmos, inaccessible otherwise. It has as main actors core-collapsing massive stars, neutron stars, coalescing compact object binaries of different flavours and stellar origin, coalescing massive black hole binaries, extreme mass ratio inspirals, and possibly the very early universe itself. Here, we highlight the science aims and describe the gravitational wave signals expected from the sources and the information gathered in it. We show that the observation of gravitational wave sources will play a transformative role in our understanding of the processes ruling the formation and evolution of stars and black holes, galaxy clustering and evolution, the nature of the strong forces in neutron star interiors, and the most mysterious interaction of Nature: gravity. The discovery, by the LIGO Scientific Collaboration and Virgo Collaboration, of the first source of gravitational waves from the cosmos GW150914, and the superb technological achievement of the space mission LISA Pathfinder herald the beginning of the new phase of exploration of the universe.

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.

Up to hundreds of black hole binaries individually resolvable by eLISA will coalesce in the Advanced LIGO/Virgo band within ten years, allowing for multi-band gravitational wave observations. Binaries formed via dynamical interactions in dense star clusters are expected to have eccentricities $e_0\sim 10^{-3}$-$10^{-1}$ at the frequencies $f_0=10^{-2}$ Hz where eLISA is most sensitive, while binaries formed in the field should have negligible eccentricity in both frequency bands. We estimate that eLISA should always be able to detect a nonzero $e_0$ whenever $e_0\gtrsim 10^{-2}$; if $e_0\sim 10^{-3}$, eLISA should detect nonzero eccentricity for a fraction $\sim 90\%$ ($\sim 25\%$) of binaries when the observation time is $T_{\rm obs}=5$ ($2$) years, respectively. Therefore eLISA observations of BH binaries have the potential to distinguish between field and cluster formation scenarios.

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.

A space-based interferometer such as eLISA could observe few to few thousands progenitors of black hole binaries (BHBs) similar to those recently detected by Advanced LIGO. Gravitational radiation circularizes the orbit during inspiral, but some BHBs retain a measurable eccentricity at the low frequencies where eLISA is most sensitive. The eccentricity of a BHB carries precious information about its formation channel: BHBs formed in the field, in globular clusters, or close to a massive black hole (MBH) have distinct eccentricity distributions in the eLISA band. We generate mock eLISA observations, folding in measurement errors, and using Bayesian model selection we study whether eLISA measurements can identify the BHB formation channel. We find that a handful of observations would suffice to tell whether BHBs were formed in the gravitational field of a MBH. Conversely, several tens of observations are needed to tell apart field formation from globular cluster formation. A five-year eLISA mission with the longest possible armlength is desirable to shed light on BHB formation scenarios.

We show that the black hole binary (BHB) coalescence rates inferred from the advanced LIGO (aLIGO) detection of GW150914 imply an unexpectedly loud GW sky at milli-Hz frequencies accessible to the evolving Laser Interferometer Space Antenna (eLISA), with several outstanding consequences. First, up to thousands of BHB will be individually resolvable by eLISA; second, millions of non resolvable BHBs will build a confusion noise detectable with signal-to-noise ratio of few to hundreds; third -- and perhaps most importantly -- up to hundreds of BHBs individually resolvable by eLISA will coalesce in the aLIGO band within ten years. eLISA observations will tell aLIGO and all electromagnetic probes weeks in advance when and where these BHB coalescences are going to occur, with uncertainties of <10s and <1deg^2. This will allow the pre-pointing of telescopes to realize coincident GW and multi-wavelength electromagnetic observations of BHB mergers. Time coincidence is critical because prompt emission associated to a BHB merger will likely have a duration comparable to the dynamical time-scale of the systems, and is only possible with low frequency GW alerts.

Supermassive black hole -- host galaxy relations are key to the computation of the expected gravitational wave background (GWB) in the pulsar timing array (PTA) frequency band. It has been recently pointed out that standard relations adopted in GWB computations are in fact biased-high. We show that when this selection bias is taken into account, the expected GWB in the PTA band is a factor of about three smaller than previously estimated. Compared to other scaling relations recently published in the literature, the median amplitude of the signal at $f=1$yr$^{-1}$ drops from $1.3\times10^{-15}$ to $4\times10^{-16}$. Although this solves any potential tension between theoretical predictions and recent PTA limits without invoking other dynamical effects (such as stalling, eccentricity or strong coupling with the galactic environment), it also makes the GWB detection more challenging.

Recent non-detection of gravitational-wave backgrounds from pulsar timing arrays casts further uncertainty on the evolution of supermassive black hole binaries. We study the capabilities of current gravitational-wave observatories to detect individual binaries and demonstrate that, contrary to conventional wisdom, some are in principle detectable throughout the Universe. In particular, a binary with rest-frame mass $\gtrsim10^{10}\,M_\odot$ can be detected by current timing arrays at arbitrarily high redshifts. The same claim will apply for less massive binaries with more sensitive future arrays. As a consequence, future searches for nanohertz gravitational waves could be expanded to target evolving high-redshift binaries. We calculate the maximum distance at which binaries can be observed with pulsar timing arrays and other detectors, properly accounting for redshift and using realistic binary waveforms.

We report on the high-precision timing of 42 radio millisecond pulsars (MSPs) observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0 extends up to mid-2014 and baselines range from 7-18 years. It forms the basis for the stochastic gravitational-wave background, anisotropic background, and continuous-wave limits recently presented by the EPTA elsewhere. The Bayesian timing analysis performed with TempoNest yields the detection of several new parameters: seven parallaxes, nine proper motions and, in the case of six binary pulsars, an apparent change of the semi-major axis. We find the NE2001 Galactic electron density model to be a better match to our parallax distances (after correction from the Lutz-Kelker bias) than the M2 and M3 models by Schnitzeler (2012). However, we measure an average uncertainty of 80\% (fractional) for NE2001, three times larger than what is typically assumed in the literature. We revisit the transverse velocity distribution for a set of 19 isolated and 57 binary MSPs and find no statistical difference between these two populations. We detect Shapiro delay in the timing residuals of PSRs J1600$-$3053 and J1918$-$0642, implying pulsar and companion masses $m_p=1.22_{-0.35}^{+0.5} \text{M}_{\odot}$, $m_c = 0.21_{-0.04}^{+0.06} \text{M}_{\odot }$ and $m_p=1.25_{-0.4}^{+0.6} \text{M}_{\odot}$, $m_c = 0.23_{-0.05}^{+0.07} \text{M}_{\odot }$, respectively. Finally, we use the measurement of the orbital period derivative to set a stringent constraint on the distance to PSRs J1012$+$5307 and J1909$-$3744, and set limits on the longitude of ascending node through the search of the annual-orbital parallax for PSRs J1600$-$3053 and J1909$-$3744.

We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only.

The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the "International" PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of ~2 and provides a 2-sigma limit on the dimensionless amplitude of a stochastic GWB of 1.7x10^{-15} at a frequency of 1 yr^{-1}. This is 1.7 times less constraining than the limit placed by (Shannon et al. 2015), due mostly to the more recent, high-quality data they used.

The sensitivity of Pulsar Timing Arrays to gravitational waves depends on the noise present in the individual pulsar timing data. Noise may be either intrinsic or extrinsic to the pulsar. Intrinsic sources of noise will include rotational instabilities, for example. Extrinsic sources of noise include contributions from physical processes which are not sufficiently well modelled, for example, dispersion and scattering effects, analysis errors and instrumental instabilities. We present the results from a noise analysis for 42 millisecond pulsars (MSPs) observed with the European Pulsar Timing Array. For characterising the low-frequency, stochastic and achromatic noise component, or "timing noise", we employ two methods, based on Bayesian and frequentist statistics. For 25 MSPs, we achieve statistically significant measurements of their timing noise parameters and find that the two methods give consistent results. For the remaining 17 MSPs, we place upper limits on the timing noise amplitude at the 95% confidence level. We additionally place an upper limit on the contribution to the pulsar noise budget from errors in the reference terrestrial time standards (below 1%), and we find evidence for a noise component which is present only in the data of one of the four used telescopes. Finally, we estimate that the timing noise of individual pulsars reduces the sensitivity of this data set to an isotropic, stochastic GW background by a factor of >9.1 and by a factor of >2.3 for continuous GWs from resolvable, inspiralling supermassive black-hole binaries with circular orbits.

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.

There is compelling evidence that most -if not all- galaxies harbour a super-massive black hole (SMBH) at their nucleus, hence binaries of these massive objects are an inevitable product of the hierarchical evolution of structures in the universe, and represent an important but thus-far elusive phase of galaxy evolution. Gas accretion via a circumbinary disc is thought to be important for the dynamical evolution of SMBH binaries, as well as in producing luminous emission that can be used to infer their properties. One plausible source of the gaseous fuel is clumps of gas formed due to turbulence and gravitational instabilities in the interstellar medium, that later fall toward and interact with the binary. In this context, we model numerically the evolution of turbulent clouds in near-radial infall onto equal-mass SMBH binaries, using a modified version of the SPH code GADGET-3. We present a total of 12 simulations that explore different possible pericentre distances and relative inclinations, and show that the formation of circumbinary discs and discs around each SMBH ('mini-discs') depend on those parameters. We also study the dynamics of the formed discs, and the variability of the feeding rate onto the SMBHs in the different configurations.

We consider the information that can be derived about massive black-hole binary populations and their formation history solely from current and possible future pulsar timing array (PTA) results. We use models of the stochastic gravitational-wave background from circular massive black hole binaries with chirp mass in the range $10^6 - 10^{11} M_\odot$ evolving solely due to radiation reaction. Our parameterised models for the black hole merger history make only weak assumptions about the properties of the black holes merging over cosmic time. We show that current PTA results place an upper limit on the black hole merger density which does not depend on the choice of a particular merger history model, however they provide no information about the redshift or mass distribution. We show that even in the case of a detection resulting from a factor of 10 increase in amplitude sensitivity, PTAs will only put weak constraints on the source merger density as a function of mass, and will not provide any additional information on the redshift distribution. Without additional assumptions or information from other observations, a detection cannot meaningfully bound the massive black hole merger rate above zero or any particular mass.

We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest EPTA dataset, which consists of ultra-precise timing data on 41 millisecond pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algorithms show evidence for the presence of such a CGW signal, indicating that the data are best described by pulsar and radiometer noise only. Depending on the adopted detection algorithm, the 95\% upper limit on the sky-averaged strain amplitude lies in the range $6\times 10^{-15}<A<1.5\times10^{-14}$ at $5{\rm nHz}<f<7{\rm nHz}$. This limit varies by a factor of five, depending on the assumed source position, and the most constraining limit is achieved towards the positions of the most sensitive pulsars in the timing array. The most robust upper limit -- obtained via a full Bayesian analysis searching simultaneously over the signal and pulsar noise on the subset of ours six best pulsars -- is $A\approx10^{-14}$. These limits, the most stringent to date at $f<10{\rm nHz}$, exclude the presence of sub-centiparsec binaries with chirp mass $\cal{M}_c>10^9$M$_\odot$ out to a distance of about 25Mpc, and with $\cal{M}_c>10^{10}$M$_\odot$ out to a distance of about 1Gpc ($z\approx0.2$). We show that state-of-the-art SMBHB population models predict $<1\%$ probability of detecting a CGW with the current EPTA dataset, consistent with the reported non-detection. We stress, however, that PTA limits on individual CGW have improved by almost an order of magnitude in the last five years. The continuing advances in pulsar timing data acquisition and analysis techniques will allow for strong astrophysical constraints on the population of nearby SMBHBs in the coming years.

We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar dataset spanning 18 yr of observations from the 2015 European Pulsar Timing Array data release. Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar System ephemeris errors, obtaining a robust 95$\%$ upper limit on the dimensionless strain amplitude $A$ of the background of $A<3.0\times 10^{-15}$ at a reference frequency of $1\mathrm{yr^{-1}}$ and a spectral index of $13/3$, corresponding to a background from inspiralling super-massive black hole binaries, constraining the GW energy density to $\Omega_\mathrm{gw}(f)h^2 < 1.1\times10^{-9}$ at 2.8 nHz. We also present limits on the correlated power spectrum at a series of discrete frequencies, and show that our sensitivity to a fiducial isotropic GWB is highest at a frequency of $\sim 5\times10^{-9}$~Hz. Finally we discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95$\%$ upper limits on the string tension, $G\mu/c^2$, characterising a background produced by a cosmic string network for a set of possible scenarios, and for a stochastic relic GWB. For a Nambu-Goto field theory cosmic string network, we set a limit $G\mu/c^2<1.3\times10^{-7}$, identical to that set by the {\it Planck} Collaboration, when combining {\it Planck} and high-$\ell$ Cosmic Microwave Background data from other experiments. For a stochastic relic background we set a limit of $\Omega^\mathrm{relic}_\mathrm{gw}(f)h^2<1.2 \times10^{-9}$, a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array.

We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. We set upper limits for a GWB from supermassive black hole binaries under power law, broken power law, and free spectral coefficient GW spectrum models. We place a 95\% upper limit on the strain amplitude (at a frequency of yr$^{-1}$) in the power law model of $A_{\rm gw} < 1.5\times 10^{-15}$. For a broken power law model, we place priors on the strain amplitude derived from simulations of Sesana (2013) and McWilliams et al. (2014). We find that the data favor a broken power law to a pure power law with odds ratios of 22 and 2.2 to one for the McWilliams and Sesana prior models, respectively. The McWilliams model is essentially ruled out by the data, and the Sesana model is in tension with the data under the assumption of a pure power law. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. We then place the most stringent limits so far on the energy density of relic GWs, $\Omega_\mathrm{gw}(f)\,h^2 < 4.2 \times 10^{-10}$, yielding a limit on the Hubble parameter during inflation of $H_*=1.6\times10^{-2}~m_{Pl}$, where $m_{Pl}$ is the Planck mass. Our limit on the cosmic string GWB, $\Omega_\mathrm{gw}(f)\, h^2 < 2.2 \times 10^{-10}$, translates to a conservative limit of $G\mu<3.3\times 10^{-8}$ - a factor of 4 better than the joint Planck and high-$l$ CMB data from other experiments.

The paucity of observed supermassive black hole binaries (SMBHBs) may imply that the gravitational wave background (GWB) from this population is anisotropic, rendering existing analyses sub-optimal. We present the first constraints on the angular distribution of a nanohertz stochastic GWB from circular, inspiral-driven SMBHBs using the $2015$ European Pulsar Timing Array data [Desvignes et al. (in prep.)]. Our analysis of the GWB in the $\sim 2 - 90$ nHz band shows consistency with isotropy, with the strain amplitude in $l>0$ spherical harmonic multipoles $\lesssim 40\%$ of the monopole value. We expect that these more general techniques will become standard tools to probe the angular distribution of source populations.

In this paper we attempt to investigate the nature of the first gravitational wave (GW) signal to be detected by pulsar timing arrays (PTAs): will it be an individual, resolved supermassive black hole binary (SBHB), or a stochastic background made by the superposition of GWs produced by an ensemble of SBHBs? To address this issue, we analyse a broad set of simulations of the cosmological population of SBHBs, that cover the entire parameter space allowed by current electromagnetic observations in an unbiased way. For each simulation, we construct the expected GW signal and identify the loudest individual sources. We then employ appropriate detection statistics to evaluate the relative probability of detecting each type of source as a function of time for a variety of PTAs; we consider the current International PTA, and speculate into the era of the Square Kilometre Array. The main properties of the first detectable individual SBHBs are also investigated. Contrary to previous work, we cast our results in terms of the detection probability (DP), since the commonly adopted criterion based on a signal-to-noise ratio threshold is statistic-dependent and may result in misleading conclusions for the statistics adopted here. Our results confirm quantitatively that a stochastic signal is more likely to be detected first (with between 75 to 93 per cent probability, depending on the array), but the DP of single-sources is not negligible. Our framework is very flexible and can be easily extended to more realistic arrays and to signal models including environmental coupling and SBHB eccentricity.

The N-independence observed in the evolution of massive black hole binaries (MBHBs) in recent simulation of merging stellar bulges suggests a simple interpretation beyond complex time-dependent relaxation processes. We conjecture that the MBHB hardening rate is equivalent to that of a binary immersed in a field of unbound stars with density $\rho$ and typical velocity $\sigma$, provided that $\rho$ and $\sigma$ are the stellar density and the velocity dispersion at the influence radius of the MBHB. By comparing direct N-body simulations to an hybrid model based on 3-body scattering experiments, we verify this hypothesis: when normalized to the stellar density and velocity dispersion at the binary influence radius, the N-body MBHB hardening rate approximately matches that predicted by 3-body scatterings in the investigated cases. The eccentricity evolution obtained with the two techniques is also in reasonable agreement. This result is particularly practical because it allows to estimate the lifetime of MBHBs forming in dry mergers based solely on the stellar density profile of the host galaxy. We briefly discuss some implications of our finding for the gravitational wave signal observable by pulsar timing arrays and for the expected population of MBHBs lurking in massive ellipticals.

We model numerically the evolution of $10^4M_\odot$ turbulent molecular clouds in near-radial infall onto $10^6M_\odot$, equal-mass super-massive black hole binaries, using a modified version of the SPH code GADGET-3. We investigate the different gas structures formed depending on the relative inclination between the binary and the cloud orbits. Our first results indicate that an aligned orbit produces mini-discs around each black hole, almost aligned with the binary; a perpendicular orbit produces misaligned mini-discs; and a counter-aligned orbit produces a circumbinary, counter-rotating ring.

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.

On a time scale of years to decades, gravitational wave (GW) astronomy will become a reality. Low frequency (nanoHz) GWs are detectable through long-term timing observations of the most stable pulsars. Radio observatories worldwide are currently carrying out observing programmes to detect GWs, with data sets being shared through the International Pulsar Timing Array project. One of the most likely sources of low frequency GWs are supermassive black hole binaries (SMBHBs), detectable as a background due to a large number of binaries, or as continuous or burst emission from individual sources. No GW signal has yet been detected, but stringent constraints are already being placed on galaxy evolution models. The SKA will bring this research to fruition. In this chapter, we describe how timing observations using SKA1 will contribute to detecting GWs, or can confirm a detection if a first signal already has been identified when SKA1 commences observations. We describe how SKA observations will identify the source(s) of a GW signal, search for anisotropies in the background, improve models of galaxy evolution, test theories of gravity, and characterise the early inspiral phase of a SMBHB system. We describe the impact of the large number of millisecond pulsars to be discovered by the SKA; and the observing cadence, observation durations, and instrumentation required to reach the necessary sensitivity. We describe the noise processes that will influence the achievable precision with the SKA. We assume a long-term timing programme using the SKA1-MID array and consider the implications of modifications to the current design. We describe the possible benefits from observations using SKA1-LOW. Finally, we describe GW detection prospects with SKA1 and SKA2, and end with a description of the expectations of GW astronomy.

We investigate the consequences of superkicks on the population of supermassive black holes (SMBHs) in the Universe residing in brightest cluster galaxies (BCGs). There is strong observational evidence that BCGs grew prominently at late times (up to a factor 2-4 in mass from z=1), mainly through mergers with satellite galaxies from the cluster, and they are known to host the most massive SMBHs ever observed. Those SMBHs are also expected to grow hierarchically, experiencing a series of mergers with other SMBHs brought in by merging satellites. Because of the net linear momentum taken away from the asymmetric gravitational wave emission, the remnant SMBH experiences a kick in the opposite direction. Kicks may be as large as ~5000 Km/s ("superkicks"), pushing the SMBHs out in the cluster outskirts for a time comparable to galaxy-evolution timescales. We predict, under a number of plausible assumptions, that superkicks can efficiently eject SMBHs from BCGs, bringing their occupation fraction down to a likely range 0.9<f<0.99 in the local Universe. Future thirty-meter-class telescopes like ELT and TMT will be capable of measuring SMBHs in hundreds of BCGs up to z=0.2, testing the occurrence of superkicks in nature and the strong-gravity regime of SMBH mergers.

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.

Pulsar timing arrays (PTAs) are designed to detect gravitational waves (GWs) at nHz frequencies. The expected dominant signal is given by the superposition of all waves emitted by the cosmological population of supermassive black hole (SMBH) binaries. Such superposition creates an incoherent stochastic background, on top of which particularly bright or nearby sources might be individually resolved. In this contribution I describe the properties of the expected GW signal, highlighting its dependence on the overall binary population, the relation between SMBHs and their hosts, and their coupling with the stellar and gaseous environment. I describe the status of current PTA efforts, and prospect of future detection and SMBH binary astrophysics.

In this paper, we address the question: What is the probability of stellar-mass black hole (BH) binaries co-existing in a globular cluster with an intermediate-mass black hole (IMBH)? Our results suggest that the detection of one or more BH binaries can strongly constrain the presence of an IMBH in most Galactic globular clusters. More specifically, the detection of one or more BH binaries could strongly indicate against the presence of an IMBH more massive than $\gtrsim 10^3$ M$_{\rm \odot}$ in roughly 80\% of the clusters in our sample. To illustrate this, we use a combination of N-body simulations and analytic methods to weigh the rate of formation of BH binaries against their ejection and/or disruption rate via strong gravitational interactions with the central (most) massive BH. The eventual fate of a sub-population of stellar-mass BHs (with or without binary companions) is for all BHs to be ejected from the cluster by the central IMBH, leaving only the most massive stellar-mass BH behind to form a close binary with the IMBH. During each phase of evolution, we discuss the rate of inspiral of the central BH-BH pair as a function of both the properties of the binary and its host cluster.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project currently observes 43 pulsars using the Green Bank and Arecibo radio telescopes. In this work we use a subset of 17 pulsars timed for a span of roughly five years (2005--2010). We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Within the timing data, we perform a search for continuous gravitational waves from individual supermassive black hole binaries in circular orbits using robust frequentist and Bayesian techniques. We find that there is no evidence for the presence of a detectable continuous gravitational wave; however, we can use these data to place the most constraining upper limits to date on the strength of such gravitational waves. Using the full 17 pulsar dataset we place a 95% upper limit on the sky-averaged strain amplitude of $h_0\lesssim 3.8\times 10^{-14}$ at a frequency of 10 nHz. Furthermore, we place 95% \emph{all sky} lower limits on the luminosity distance to such gravitational wave sources finding that the $d_L \gtrsim 425$ Mpc for sources at a frequency of 10 nHz and chirp mass $10^{10}{\rm M}_{\odot}$. We find that for gravitational wave sources near our best timed pulsars in the sky, the sensitivity of the pulsar timing array is increased by a factor of $\sim$4 over the sky-averaged sensitivity. Finally we place limits on the coalescence rate of the most massive supermassive black hole binaries.

This paper presents a technique to search for supermassive black hole binaries (MBHBs) in the Sloan Digital Sky Survey (SDSS). The search is based on the peculiar properties of merging galaxies that are found in a mock galaxy catalog from the Millennium Simulation. MBHBs are expected to be the main gravitational wave (GW) sources for pulsar timing arrays (PTAs); however, it is still unclear if the observed GW signal will be produced by a few single MBHBs, or if it will have the properties of a stochastic background. The goal of this work is to produce a map of the sky in which each galaxy is assigned a probability of having suffered a recent merger, and of hosting a MBHB that could be detected by PTAs. This constitutes a step forward in the understanding of the expected PTA signal: the skymap can be used to investigate the clustering properties of PTA sources and the spatial distribution of the observable GW signal power; moreover, galaxies with the highest probabilities could be used as inputs in targeted searches for individual GW sources. We also investigate the distribution of neighboring galaxies around galaxies hosting MBHBs, finding that the most likely detectable PTA sources are located in dense galaxy environments. Different techniques are used in the search, including Bayesian and Machine Learning algorithms, with consistent outputs. Our method generates a list of galaxies classified as MBHB hosts, that can be combined with other searches to effectively reduce the number of misclassifications. The spectral coverage of the SDSS reaches less than a fifth of the sky, and the catalog becomes severely incomplete at large redshifts; however, this technique can be applied in the future to larger catalogs to obtain complete, observationally-based information of the expected GW signal detectable by PTAs.

We study the interplay between mass transfer, accretion and gravitational torques onto a black hole binary migrating in a self-gravitating, retrograde circumbinary disc. A direct comparison with an identical prograde disc shows that: (i) because of the absence of resonances, the cavity size is a factor a(1+e) smaller for retrograde discs; (ii) nonetheless the shrinkage of a circular binary semi--major axis, a, is identical in both cases; (iii) a circular binary in a retrograde disc remains circular while eccentric binaries grow more eccentric. For non-circular binaries, we measure the orbital decay rates and the eccentricity growth rates to be exponential as long as the binary orbits in the plane of its disc. Additionally, for these co-planar systems, we find that interaction (~ non--zero torque) stems only from the cavity edge plus a(1+e) in the disc, i.e. for dynamical purposes, the disc can be treated as a annulus of small radial extent. We find that simple 'dust' models in which the binary- disc interaction is purely gravitational can account for all main numerical results, both for prograde and retrograde discs. Furthermore, we discuss the possibility of an instability occurring for highly eccentric binaries causing it to leave the disc plane, secularly tilt and converge to a prograde system. Our results suggest that there are two stable configurations for binaries in self-gravitating discs: the special circular retrograde case and an eccentric (e~ 0.6) prograde configuration as a stable attractor.

Massive black hole binaries (MBHBs) are unavoidable outcomes of the hierarchical structure formation process, and according to the theory of general relativity are expected to be the loudest gravitational wave (GW) sources in the Universe. In this article I provide a broad overview of MBHBs as GW sources. After reviewing the basics of GW emission from binary systems and of MBHB formation, evolution and dynamics, I describe in some details the connection between binary properties and the emitted gravitational waveform. Direct GW observations will provide an unprecedented wealth of information about the physical nature and the astrophysical properties of these extreme objects, allowing to reconstruct their cosmic history, dynamics and coupling with their dense stellar and gas environment. In this context I describe ongoing and future efforts to make a direct detection with space based interferometry and pulsar timing arrays, highlighting the invaluable scientific payouts of such enterprises.

Pulsar timing arrays (PTAs) are designed to detect the predicted gravitational wave (GW) background produced by a cosmological population of supermassive black hole (SMBH) binaries. In this contribution I review the physics of such GW background, highlighting its dependence on the overall binary population, the relation between SMBHs and their hosts, and their coupling with the stellar and gaseous environment. The latter is particularly relevant when it drives the binaries to extreme eccentricities (e>0.9), which might be the case for stellar-driven systems. This causes a substantial suppression of the low frequency signal, potentially posing a serious threat to the effectiveness of PTA observations. A future PTA detection will allow to directly observe for the first time subparsec SMBH binaries on their way to the GW driven coalescence, providing important answers of the outstanding questions related to the physics underlying the formation and evolution of these spectacular sources.

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.

An astrophysical population of supermassive black hole binaries is thought to be the strongest source of gravitational waves in the frequency range covered by Pulsar Timing Arrays (PTAs). A potential cause for concern is that the standard cross-correlation method used in PTA data analysis assumes that the signals are isotropically distributed and Gaussian random, while the signals from a black hole population are likely to be anisotropic and deterministic. Here we argue that while the conventional analysis is not optimal, it is not hopeless either, as the standard Hellings-Downs correlation curve turns out to hold for point sources, and the small effective number of signal samples blurs the distinction between Gaussian and deterministic signals. Possible improvements to the standard cross-correlation analysis that account for the anisotropy of the signal are discussed.

Within this decade gravitational wave detection will open a new observational window on the Universe. Advanced ground based interferometers covering the kHz frequency range will be online by 2016, and it is foreseeable the announcement of a first detection within five years. At the same time, a worldwide effort of detecting low frequency waves (in the nHz regime) by timing ultra-precise millisecond pulsars is rapidly growing, possibly leading to a positive detection within this decade. The mHz regime, bridging these two windows, is the realm of space based interferometers, which might be launched in the late 20s. I provide here a short overview of the scientific payouts of gravitational wave astronomy, focusing the discussion on the low frequency regime (pulsar timing and space based interferometry). A detailed discussion of advanced ground based interferometer can be found in Patrick Brady's contribution to this proceedings series.

Observations of the Galactic Center (GC) have accumulated a multitude of "forensic" evidence indicating that several million years ago the center of the Milky Way galaxy was teaming with starforming and accretion-powered activity -- this paints a rather different picture from the GC as we understand it today. We examine a possibility that this epoch of activity could have been triggered by the infall of a satellite galaxy into the Milky Way which began at the redshift of 10 and ended few million years ago with a merger of the Galactic supermassive black hole with an intermediate mass black hole brought in by the inspiralling satellite.

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

In this letter we carry out the first systematic investigation of the expected gravitational wave (GW) background generated by supermassive black hole (SMBH) binaries in the nHz frequency band accessible to pulsar timing arrays (PTAs). We take from the literature several estimates of the redshift dependent galaxy mass function and of the fraction of close galaxy pairs to derive a wide range of galaxy merger rates. We then exploit empirical black hole-host relations to populate merging galaxies with SMBHs. The result of our procedure is a collection of a large number of phenomenological SMBH binary merger rates consistent with current observational constraints on the galaxy assembly at z<1.5. For each merger rate we compute the associated GW signal, eventually producing a large set of estimates of the nHz GW background that we use to infer confidence intervals of its expected amplitude. When considering the most recent SMBH-host relations, accounting for ultra-massive black holes in brightest cluster galaxies, we find that the nominal $1\sigma$ interval of the expected GW signal is only a factor of 3-to-10 below current PTA limits, implying a non negligible chance of detection in the next few years.

Pulsar timing arrays (PTAs) might detect gravitational waves (GWs) from massive black hole (MBH) binaries within this decade. The signal is expected to be an incoherent superposition of several nearly-monochromatic waves of different strength. The brightest sources might be individually resolved, and the overall deconvolved, at least partially, in its individual components. In this paper we extend the maximum-likelihood based method developed in Babak & Sesana 2012, to search for individual MBH binaries in PTA data. We model the signal as a collection of circular monochromatic binaries, each characterized by three free parameters: two angles defining the sky location, and the frequency. We marginalize over all other source parameters and we apply an efficient multi-search genetic algorithm to maximize the likelihood function and look for sources in synthetic datasets. On datasets characterized by white Gaussian noise plus few injected sources with signal-to-noise ratio (SNR) in the range 10-60, our search algorithm performs well, recovering all the injections with no false positives. Individual source SNRs are estimated within few % of the injected values, sky locations are recovered within few degrees, and frequencies are determined with sub-Fourier bin precision.

Space based gravitational wave astronomy will open a completely new window on the Universe and massive black holes binaries are expected to be among the primary actors on this upcoming stage. The New Gravitational-wave Observatory (NGO) is a space interferometer proposal derived from the former Laser Interferometer Space Antenna (LISA) concept. We describe here its capabilities of observing massive black hole binaries throughout the Universe, measuring their relevant parameters (masses, spins, distance to the observer) to high precision. The statistical properties of the population of detected systems can be used to constrain the massive black hole cosmic history, providing deep insights into the faint, high redshift Universe.

We review the expected science performance of the New Gravitational-Wave Observatory (NGO, a.k.a. eLISA), a mission under study by the European Space Agency for launch in the early 2020s. eLISA will survey the low-frequency gravitational-wave sky (from 0.1 mHz to 1 Hz), detecting and characterizing a broad variety of systems and events throughout the Universe, including the coalescences of massive black holes brought together by galaxy mergers; the inspirals of stellar-mass black holes and compact stars into central galactic black holes; several millions of ultracompact binaries, both detached and mass transferring, in the Galaxy; and possibly unforeseen sources such as the relic gravitational-wave radiation from the early Universe. eLISA's high signal-to-noise measurements will provide new insight into the structure and history of the Universe, and they will test general relativity in its strong-field dynamical regime.

Intermediate mass black holes (IMBHs) are among the most elusive objects in contemporary astrophysics. Both theoretical and observational evidence of their existence is subject of debate. Conversely, both theory and observations confirm the presence of a large population of millisecond pulsars (MSPs) with low mass companions residing in globular cluster (GC) centers. If IMBHs are common in GC centers as well, then dynamical interactions will inevitably break up many of these binaries, causing the ejection of several fast MSPs in the Galactic halo. Such population of fast halo MSPs, hard to produce with 'standard' MSP generation mechanisms, would provide a strong, albeit indirect, evidence of the presence of a substantial population of IMBHs in GCs. In this paper we study in detail the dynamical formation and evolution of such fast MSPs population, highlighting the relevant observational properties and assessing detection prospects with forthcoming radio surveys.

We study the interplay between gas accretion and gravity torques in changing a binary elements and its total angular momentum (L) budget. Especially, we analyse the physical origin of the gravity torques (T_g) and their location within the disc. We analyse 3D SPH simulations of the evolution of initially quasi-circular massive black hole binaries (BHBs) residing in the central hollow of massive self-gravitating circumbinary discs. We use different thermodynamics within the cavity and for the numerical size of the black holes to show that (i) the BHB eccentricity growth found previously is a general result, independent of the accretion and the adopted thermodynamics; (ii) the semi-major axis decay depends both on the T_g and on the interplay with the disc-binary L-transfer due to accretion; (iii) the spectral structure of the T_g is predominately caused by disc edge overdensities and spiral arms developing in the body of the disc and, in general, does not reflect directly the period of the binary; (iv) the net T_g changes sign across the BHB corotation radius. We quantify the relative importance of the two, which appear to depend on the thermodynamical properties of the instreaming gas, and which is crucial in assessing the disc-binary L-transfer; (v) the net torque manifests as a purely kinematic (non-resonant) effect as it stems from the cavity, where the material flows in and out in highly eccentric orbits. Both accretion onto the black holes and the interaction with gas streams inside the cavity must be taken into account to assess the fate of the BHB. Moreover, the total torque exerted by the disc affects L(BHB) by changing all the elements (mass, mass ratio, eccentricity, semimajor axis) of the BHB. Common prescriptions equating tidal torque to semi-major axis shrinking might therefore be poor approximations for real astrophysical systems.

This document introduces the exciting and fundamentally new science and astronomy that the European New Gravitational Wave Observatory (NGO) mission (derived from the previous LISA proposal) will deliver. The mission (which we will refer to by its informal name "eLISA") will survey for the first time the low-frequency gravitational wave band (about 0.1 mHz to 1 Hz), with sufficient sensitivity to detect interesting individual astrophysical sources out to z = 15. The eLISA mission will discover and study a variety of cosmic events and systems with high sensitivity: coalescences of massive black holes binaries, brought together by galaxy mergers; mergers of earlier, less-massive black holes during the epoch of hierarchical galaxy and black-hole growth; stellar-mass black holes and compact stars in orbits just skimming the horizons of massive black holes in galactic nuclei of the present era; extremely compact white dwarf binaries in our Galaxy, a rich source of information about binary evolution and about future Type Ia supernovae; and possibly most interesting of all, the uncertain and unpredicted sources, for example relics of inflation and of the symmetry-breaking epoch directly after the Big Bang. eLISA's measurements will allow detailed studies of these signals with high signal-to-noise ratio, addressing most of the key scientific questions raised by ESA's Cosmic Vision programme in the areas of astrophysics and cosmology. They will also provide stringent tests of general relativity in the strong-field dynamical regime, which cannot be probed in any other way. This document not only describes the science but also gives an overview on the mission design and orbits.

We study the capability of a pulsar timing array (PTA) to individually resolve and localize in the sky monochromatic gravitational wave (GW) sources. Given a cosmological population of inspiralling massive black hole binaries, their observable signal in the PTA domain is expected to be a superposition of several nearly-monochromatic GWs of different strength. In each frequency bin, the signal is neither a stochastic background nor perfectly resolvable in its individual components. In this context, it is crucial to explore how the information encoded in the spatial distribution of the array of pulsars might help recovering the origin of the GW signal, by resolving individually and locating in the sky the strongest sources. In this paper we develop a maximum-likelihood based method finalized to this purpose. We test the algorithm against noiseless data showing that up to P/3 sources can be resolved and localized in the sky by an array of P pulsars. We validate the code by performing blind searches on both noiseless and noisy data containing an unknown number of signals with different strengths. Even without employing any proper search algorithm, our analysis procedure performs well, recovering and correctly locating in the sky almost all the injected sources.

We outline the eccentricity evolution of sub-parsec massive black hole binaries (MBHBs) forming in galaxy mergers. In both stellar and gaseous environments, MBHBs are expected to grow large orbital eccentricities before they enter the gravitational wave (GW) observational domain. We re--visit the predicted eccentricities detectable by space based laser interferometers (as the proposed ELISA/NGO) for both environments. Close to coalescence, many MBHBs will still maintain detectable eccentricities, spanning a broad range from <10^{-5} up to <~ 0.5. Stellar and gas driven dynamics lead to distinct distributions, with the latter favoring larger eccentricities. At larger binary separations, when emitted GWs will be observed by pulsar timing arrays (PTAs), the expected eccentricities are usually quite large, in the range 0.01-0.7, which poses an important issue for signal modelling and detection algorithms. In this window, large eccentricities also have implications on proposed electromagnetic counterparts to the GW signal, which we briefly review.

We study the evolution of the orientation of the orbital plane of massive black hole binaries (BHBs) in rotating stellar systems in which the total angular momentum of the stellar cusp is misaligned with respect to that of the binary. We compare results from direct summation N-body simulations with predictions from a simple theoretical model. We find that the same encounters between cusp stars and the BHB that are responsible for the hardening and eccentricity evolution of the binary, lead to a reorientation of the binary orbital plane. In particular, binaries whose angular momentum is initially misaligned with respect to that of the stellar cusp tend to realign their orbital planes with the angular momentum of the cusp on a timescale of a few hardening times. This is due to angular momentum exchange between stars and the BHB during close encounters, and may have important implications for the relative orientation of host galaxies and radio jets.

The study of the dynamical evolution of massive black hole pairs in mergers is crucial in the context of a hierarchical galaxy formation scenario. The timescales for the formation and the coalescence of black hole binaries are still poorly constrained, resulting in large uncertainties in the expected rate of massive black hole binaries detectable in the electromagnetic and gravitational wave spectra. Here we review the current theoretical understanding of the black hole pairing in galaxy mergers, with a particular attention to recent developments and open issues. We conclude with a review of the expected observational signatures of massive binaries, and of the candidates discussed in literature to date.

I review our current understanding of massive black hole (MBH) formation and evolution along the cosmic history. After a brief introductory overview of the relevance of MBHs in the hierarchical structure formation paradigm, I discuss the main viable channels for seed BH formation at high redshift and for their subsequent mass growth and spin evolution. The emerging hierarchical picture, where MBHs grow through merger triggered accretion episodes, acquiring their mass while shining as quasars, is overall robust, but too simplistic to explain the diversity observed in MBH phenomenology. I briefly discuss which future observations will help to shed light on the MBH cosmic history in the near future, paying particular attention to the upcoming gravitational wave window.

Space-based gravitational-wave detectors, such as LISA or a similar ESA-led mission, will offer unique opportunities to test general relativity. We study the bounds that space-based detectors could place on the graviton Compton wavelength \lambda_g=h/(m_g c) by observing multiple inspiralling black hole binaries. We show that while observations of individual inspirals will yield mean bounds \lambda_g~3x10^15 km, the combined bound from observing ~50 events in a two-year mission is about ten times better: \lambda_g~3x10^16 km (m_g~4x10^-26 eV). The bound improves faster than the square root of the number of observed events, because typically a few sources provide constraints as much as three times better than the mean. This result is only mildly dependent on details of black hole formation and detector characteristics. The bound achievable in practice should be one order of magnitude better than this figure (and hence almost competitive with the static, model-dependent bounds from gravitational effects on cosmological scales), because our calculations ignore the merger/ringdown portion of the waveform. The observation that an ensemble of events can sensibly improve the bounds that individual binaries set on \lambda_g applies to any theory whose deviations from general relativity are parametrized by a set of global parameters.

In the last few years before merger, supermassive black hole binaries will rapidly inspiral and precess in a magnetic field imposed by a surrounding circumbinary disk. Multiple simulations suggest this relative motion will convert some of the local energy to a Poynting-dominated outflow, with a luminosity 10^{43} erg/s * (B/10^4 G)^2(M/10^8 Msun)^2 (v/0.4 c)^2, some of which may emerge as synchrotron emission at frequencies near 1 GHz where current and planned wide-field radio surveys will operate. On top of a secular increase in power on the gravitational wave inspiral timescale, orbital motion will produce significant, detectable modulations, both on orbital periods and (if black hole spins are not aligned with the binary's total angular momenta) spin-orbit precession timescales. Because the gravitational wave merger time increases rapidly with separation, we find vast numbers of these transients are ubiquitously predicted, unless explicitly ruled out (by low efficiency $\epsilon$) or obscured (by accretion geometry f_{geo}). If the fraction of Poynting flux converted to radio emission times the fraction of lines of sight accessible $f_{geo}$ is sufficiently large (f_{geo} \epsilon > 2\times 10^{-4} for a 1 year orbital period), at least one event is accessible to future blind surveys at a nominal 10^4 {deg}^2 with 0.5 mJy sensitivity. Our procedure generalizes to other flux-limited surveys designed to investigate EM signatures associated with many modulations produced by merging SMBH binaries.

We demonstrate that very massive (>10^8\msun), cosmologically nearby (z<1) black hole binaries (MBHBs), which are primary targets for ongoing and upcoming pulsar timing arrays (PTAs), are particularly appealing multimessenger carriers. According to current models for massive black hole formation and evolution, the planned Square Kilometer Array (SKA) will collect gravitational wave signals from thousands of such massive systems, being able to individually resolve and locate in the sky several of them (maybe up to a hundred). By employing a standard model for the evolution of MBHBs in circumbinary discs, with the aid of dedicated numerical simulations, we characterize the gas-binary interplay, identifying possible electromagnetic signatures of the PTA sources. We concentrate our investigation on two particularly promising scenarios in the high energy domain, namely, the detection of X-ray periodic variability and of double broad K\alpha iron lines. Up to several hundreds of periodic X-ray sources with a flux >10^-13 erg s^-1 cm^-2 will be in the reach of upcoming X-ray observatories. Double relativistic K\alpha lines may be observable in a handful of low redshift (z<0.3) sources by proposed deep X-ray probes, such as Athena. (Abridged)

Supermassive black holes presumably grow through numerous mergers throughout cosmic time. During each merger, supermassive black hole binaries are surrounded by a circumbinary accretion disk that imposes a significant (~1e4 G for a binary of 1e8 Msun) magnetic field. The motion of the binary through that field will convert the field energy to Poynting flux, with a luminosity ~1e43 erg/s (B/1e4 G)^2 (M/1e8 Msun)^2, some of which may emerge as synchrotron emission at frequencies near 1 GHz where current and planned wide-field radio surveys will operate. We find that the short timescales of many mergers will limit their detectability with most planned blind surveys to <1 per year over the whole sky, independent of the details of the emission process and flux distribution. Including an optimistic estimate for the radio flux makes detection even less likely, with <1 mergers per year over the whole sky. However, wide-field radio instruments may be able to localize systems identified in advance of merger by gravitational waves. Further, radio surveys may be able to detect the weaker emission produced by the binary's motion as it is modulated by spin-orbit precession and inspiral well in advance of merger.

In this letter we study the eccentricity evolution of a massive black hole (MBH) binary (MBHB) embedded in a rotating stellar cusp. Following the observation that stars on counter-rotating (with respect to the MBHB) orbits extract angular momentum from the binary more efficiently then their co-rotating counterparts, the eccentricity evolution of the MBHB must depend on the degree of co-rotation (counter-rotation) of the surrounding stellar distribution. Using an hybrid scheme that couples numerical three-body scatterings to an analytical formalism for the cusp-binary interaction, we verify this hypothesis by evolving the MBHB in spherically symmetric cusps with different fractions F of co-rotating stars. Consistently with previous works, binaries in isotropic cusps (F=0.5) tend to increase their eccentricity, and when F approaches zero (counter-rotating cusps) the eccentricity rapidly increases to almost unity. Conversely, binaries in cusps with a significant degree of co-rotation (F>0.7) tend to become less and less eccentric, circularising quite quickly for F approaching unity. Direct N-body integrations performed to test the theory, corroborate the results of the hybrid scheme, at least at a qualitative level. We discuss quantitative differences, ascribing their origin to the oversimplified nature of the hybrid approach.

We study the dynamics of supermassive black hole binaries embedded in circumbinary gaseous discs, with the SPH code Gadget-2. The sub-parsec binary (of total mass M and mass ratio q=1/3) has excavated a gap and transfers its angular momentum to the self--gravitating disc (M_disc=0.2 M). We explore the changes of the binary eccentricity e, by simulating a sequence of binary models that differ in the initial eccentricity e_0, only. In initially low-eccentric binaries, the eccentricity increases with time, while in high-eccentric binaries e declines, indicating the existence of a limiting eccentricity e_crit that is found to fall in the interval [0.6,0.8]. We also present an analytical interpretation for this saturation limit. An important consequence of the existence of e_crit is the detectability of a significant residual eccentricity e_LISA} by the proposed gravitational wave detector LISA. It is found that at the moment of entering the LISA frequency domain e_LISA ~ 10^{-3}-10^{-2}; a signature of its earlier coupling with the massive circumbinary disc. We also observe large periodic inflows across the gap, occurring on the binary and disc dynamical time scales rather than on the viscous time. These periodic changes in the accretion rate (with amplitudes up to ~100%, depending on the binary eccentricity) can be considered a fingerprint of eccentric sub-parsec binaries migrating inside a circumbinary disc.

LISA might detect gravitational waves from mergers of massive black hole binaries strongly lensed by intervening galaxies (Sereno et al. 2010). The detection of multiple gravitational lensing events would provide a new tool for cosmography. Constraints on cosmological parameters could be placed by exploiting either lensing statistics of strongly lensed sources or time delay measurements of lensed gravitational wave signals. These lensing methods do not need the measurement of the redshifts of the sources and the identification of their electromagnetic counterparts. They would extend cosmological probes to redshift z <= 10 and are then complementary to other lower or higher redshift tests, such as type Ia supernovae or cosmic microwave background. The accuracy of lensing tests strongly depends on the formation history of the merging binaries, and the related number of total detectable multiple images. Lensing amplification might also help to find the host galaxies. Any measurement of the source redshifts would allow to exploit the distance-redshift test in combination with lensing methods. Time-delay analyses might measure the Hubble parameter H_0 with accuracy of >= 10 km s^{-1}Mpc^{-1}. With prior knowledge of H_0, lensing statistics and time delays might constrain the dark matter density (delta Omega_M >= 0.08, due to parameter degeneracy). Inclusion of our methods with other available orthogonal techniques might significantly reduce the uncertainty contours for Omega_M and the dark energy equation of state.

In this article we give a brief review of the fundamental physics that can be done with the future space-based gravitational wave detector LISA. This includes detection of gravitational wave bursts coming from cosmic strings, measuring a stochastic gravitational wave background, mapping spacetime around massive compact objects in galactic nuclei with extreme-mass-ratio inspirals and testing the predictions of General Relativity for the strong dynamical fields of inspiralling binaries. We give particular attention to new results which show the capability of LISA to constrain cosmological parameters using observations of coalescing massive Black Hole binaries.

Gravitational wave signals from coalescing Massive Black Hole (MBH) binaries could be used as standard sirens to measure cosmological parameters. The future space based gravitational wave observatory Laser Interferometer Space Antenna (LISA) will detect up to a hundred of those events, providing very accurate measurements of their luminosity distances. To constrain the cosmological parameters we also need to measure the redshift of the galaxy (or cluster of galaxies) hosting the merger. This requires the identification of a distinctive electromagnetic event associated to the binary coalescence. However, putative electromagnetic signatures may be too weak to be observed. Instead, we study here the possibility of constraining the cosmological parameters by enforcing statistical consistency between all the possible hosts detected within the measurement error box of a few dozen of low redshift (z<3) events. We construct MBH populations using merger tree realizations of the dark matter hierarchy in a LambdaCDM Universe, and we use data from the Millennium simulation to model the galaxy distribution in the LISA error box. We show that, assuming that all the other cosmological parameters are known, the parameter w describing the dark energy equation of state can be constrained to a 4-8% level (2sigma error), competitive with current uncertainties obtained by type Ia supernovae measurements, providing an independent test of our cosmological model.

Tidal stellar disruptions have traditionally been discussed as a probe of the single, massive black holes (MBHs) that are dormant in the nuclei of galaxies. In Chen et al. (2009), we used numerical scattering experiments to show that three-body interactions between bound stars in a stellar cusp and a non-evolving "hard" MBH binary will also produce a burst of tidal disruptions, caused by a combination of the secular "Kozai effect" and by close resonant encounters with the secondary hole. Here we derive basic analytical scalings of the stellar disruption rates with the system parameters, assess the relative importance of the Kozai and resonant encounter mechanisms as a function of time, discuss the impact of general relativistic (GR) and extended stellar cusp effects, and develop a hybrid model to self-consistently follow the shrinking of an MBH binary in a stellar background, including slingshot ejections and tidal disruptions. In the case of a fiducial binary with primary hole mass M_1=10^7\msun and mass ratio q=M_2/M_1=1/81, embedded in an isothermal cusp, we derive a stellar disruption rate \dot{N_*} ~ 0.2/yr lasting ~ 3X10^5 yr. This rate is 3 orders of magnitude larger than the corresponding value for a single MBH fed by two-body relaxation, confirming our previous findings. For q<<0.01, the Kozai/chaotic effect could be quenched due to GR/cusp effects by an order of magnitude, but even in this case the stellar-disruption rate is still two orders of magnitude larger than that given by standard relaxation processes around a single MBH. Our results suggest that >~10% of the tidal-disruption events may originate in MBH binaries.

LISA should detect gravitational waves from tens to hundreds of systems containing black holes with mass in the range from 10 thousand to 10 million solar masses. Black holes in this mass range are not well constrained by current electromagnetic observations, so LISA could significantly enhance our understanding of the astrophysics of such systems. In this paper, we describe a framework for combining LISA observations to make statements about massive black hole populations. We summarise the constraints that LISA observations of extreme-mass-ratio inspirals might be able to place on the mass function of black holes in the LISA range. We also describe how LISA observations can be used to choose between different models for the hierarchical growth of structure in the early Universe. We consider four models that differ in their prescription for the initial mass distribution of black hole seeds, and in the efficiency of accretion onto the black holes. We show that with as little as 3 months of LISA data we can clearly distinguish between these models, even under relatively pessimistic assumptions about the performance of the detector and our knowledge of the gravitational waveforms.

The massive black holes we observe in galaxies today are the natural end-product of a complex evolutionary path, in which black holes seeded in proto-galaxies at high redshift grow through cosmic history via a sequence of mergers and accretion episodes. Electromagnetic observations probe a small subset of the population of massive black holes (namely, those that are active or those that are very close to us), but planned space-based gravitational-wave observatories such as the Laser Interferometer Space Antenna (LISA) can measure the parameters of ``electromagnetically invisible'' massive black holes out to high redshift. In this paper we introduce a Bayesian framework to analyze the information that can be gathered from a set of such measurements. Our goal is to connect a set of massive black hole binary merger observations to the underlying model of massive black hole formation. In other words, given a set of observed massive black hole coalescences, we assess what information can be extracted about the underlying massive black hole population model. For concreteness we consider ten specific models of massive black hole formation, chosen to probe four important (and largely unconstrained) aspects of the input physics used in structure formation simulations: seed formation, metallicity ``feedback'', accretion efficiency and accretion geometry. For the first time we allow for the possibility of ``model mixing'', by drawing the observed population from some combination of the ``pure'' models that have been simulated. A Bayesian analysis allows us to recover a posterior probability distribution for the ``mixing parameters'' that characterize the fractions of each model represented in the observed distribution. Our work shows that LISA has enormous potential to probe the underlying physics of structure formation.

We discuss strong gravitational lensing of gravitational waves from merging of massive black hole binaries in the context of the LISA mission. Detection of multiple events would provide invaluable information on competing theories of gravity, evolution and formation of structures and, with complementary observations, constraints on H_0 and other cosmological parameters. Most of the optical depth for lensing is provided by intervening massive galactic halos, for which wave optics effects are negligible. Probabilities to observe multiple events are sizable for a broad range of formation histories. For the most optimistic models, up to 4 multiple events with a signal to noise ratio >= 8 are expected in a 5-year mission. Chances are significant even for conservative models with either light (<= 60%) or heavy (<= 40%) seeds. Due to lensing amplification, some intrinsically too faint signals are brought over threshold (<= 2 per year).

Pulsar timing arrays (PTAs) measure nHz frequency gravitational waves (GWs) generated by orbiting massive black hole binaries (MBHBs) with periods between 0.1-10 yr. Previous studies on the nHz GW background assumed that the inspiral is purely driven by GWs. However, torques generated by a gaseous disk can shrink the binary much more efficiently than GW emission, reducing the number of binaries at these separations. We use simple disk models for the circumbinary gas and for the binary-disk interaction to follow the orbital decay of MBHBs through physically distinct regions of the disk, until GWs take over their evolution. We extract MBHB cosmological merger rates from the Millennium simulation, generate Monte Carlo realizations of a population of gas driven binaries, and calculate the corresponding GW amplitudes of the most luminous individual binaries and the stochastic GW background. For steady state alpha-disks with alpha>0.1 we find that the nHz GW background can be significantly modified. The number of resolvable binaries is however not changed by the presence of gas; we predict 1-10 individually resolvable sources to stand above the noise for a 1-50 ns timing precision. Gas driven migration reduces predominantly the number of small total mass or unequal mass ratio binaries, which leads to the attenuation of the mean stochastic GW--background, but increases the detection significance of individually resolvable binaries. The results are sensitive to the model of binary--disk interaction. The GW background is not attenuated significantly for time-dependent models of Ivanov, Papaloizou, & Polnarev (1999).

We construct evolutionary tracks for massive black hole binaries (MBHBs) embedded in a surrounding distribution of stars. The dynamics of the binary is evolved by taking into account the erosion of the central stellar cusp bound to the massive black holes, the scattering of unbound stars feeding the binary loss cone, and the emission of gravitational waves (GWs). Stellar dynamics is treated in a hybrid fashion by coupling the results of numerical 3-body scattering experiments of bound and unbound stars to an analytical framework for the evolution of the stellar density distribution and for the efficiency of the binary loss cone refilling. Our main focus is on the behaviour of the binary eccentricity, in the attempt of addressing its importance in the merger process and its possible impact for GW detection with the planned Laser Interferometer Space Antenna ({\it LISA}), and ongoing and forthcoming pulsar timing array (PTA) campaigns. We produce a family of evolutionary tracks extensively sampling the relevant parameters of the system which are the binary mass, mass ratio and initial eccentricity, the slope of the stellar density distribution, its normalization and the efficiency of loss cone refilling. We find that, in general, stellar dynamics causes a dramatic increase of the MBHB eccentricity, especially for initially already mildly eccentric and/or unequal mass binaries. When applied to standard MBHB population models, our results predict eccentricities in the ranges $10^{-3}-0.2$ and $0.03-0.3$ for sources detectable by {\it LISA} and PTA respectively. Such figures may have a significant impact on the signal modelling, on source detection, and on the development of parameter estimation algorithms.

The observation of massive black hole binaries (MBHBs) with Pulsar Timing Arrays (PTAs) is one of the goals of gravitational wave astronomy in the coming years. Massive (>10^8 solar masses) and low-redshift (< 1.5) sources are expected to be individually resolved by up-coming PTAs, and our ability to use them as astrophysical probes will depend on the accuracy with which their parameters can be measured. In this paper we estimate the precision of such measurements using the Fisher-information-matrix formalism. We restrict to "monochromatic" sources. In this approximation, the system is described by seven parameters and we determine their expected statistical errors as a function of the number of pulsars in the array, the array sky coverage, and the signal-to-noise ratio (SNR) of the signal. At fixed SNR, the gravitational wave astronomy capability of a PTA is achieved with ~20 pulsars; adding more pulsars (up to 1000) to the array reduces the source error-box in the sky \Delta\Omega by a factor ~5 and has negligible consequences on the statistical errors on the other parameters. \Delta\Omega improves as 1/SNR^2 and the other parameters as 1/SNR. For a fiducial PTA of 100 pulsars uniformly distributed in the sky and a coherent SNR = 10, we find \Delta\Omega~40 deg^2, a fractional error on the signal amplitude of ~30% (which constraints only very poorly the chirp mass - luminosity distance combination M_c^{5/3}/D_L), and the source inclination and polarization angles are recovered at the ~0.3 rad level. The ongoing Parkes PTA is particularly sensitive to systems located in the southern hemisphere, where at SNR = 10 the source position can be determined with \Delta\Omega ~10 deg^2, but has poorer performance for sources in the northern hemisphere. (Abridged)

Massive black holes are key ingredients of the assembly and evolution of cosmic structures. Pulsar Timing Arrays (PTAs) currently provide the only means to observe gravitational radiation from massive black hole binary systems with masses >10^7 solar masses. The whole cosmic population produces a signal consisting of two components: (i) a stochastic background resulting from the incoherent superposition of radiation from the all the sources, and (ii) a handful of individually resolvable signals that raise above the background level and are produced by sources sufficiently close and/or massive. Considering a wide range of massive black hole binary assembly scenarios, we investigate both the level and shape of the background and the statistics of resolvable sources. We predict a characteristic background amplitude in the interval h_c(f = 10^-8 Hz)~5*10^-16 - 5*10^-15, within the detection range of the complete Parkes PTA. We also quantify the capability of PTAs of measuring the parameters of individual sources, focusing on monochromatic signals produced by binaries in circular orbit. We investigate how the results depend on the number and distribution of pulsars in the array, by computing the variance-covariance matrix of the parameter measurements. For plausible Square Kilometre Array (SKA) observations (100 pulsars uniformly distributed in the sky), and assuming a coherent signal-to-noise ratio of 10, the sky position of massive black hole binaries can be located within a ~40deg^2 error box, opening promising prospects for detecting a putative electromagnetic counterpart to the gravitational wave emission. The planned SKA, can plausibly observe these unique systems, although the number of detections is likely to be small. (Abridged)

The International Pulsar Timing Array project combines observations of pulsars from both Northern and Southern hemisphere observatories with the main aim of detecting ultra-low frequency (~10^-9 to 10^-8 Hz) gravitational waves. Here we introduce the project, review the methods used to search for gravitational waves emitted from coalescing supermassive binary black-hole systems in the centres of merging galaxies and discuss the status of the project.

Supermassive black holes (SMBHs) found in the centers of many galaxies have been recognized to play a fundamental active role in the cosmological structure formation process. In hierarchical formation scenarios, SMBHs are expected to form binaries following the merger of their host galaxies. If these binaries do not coalesce before the merger with a third galaxy, the formation of a black hole triple system is possible. Numerical simulations of the dynamics of triples within galaxy cores exhibit phases of very high eccentricity (as high as $e \sim 0.99$). During these phases, intense bursts of gravitational radiation can be emitted at orbital periapsis. This produces a gravitational wave signal at frequencies substantially higher than the orbital frequency. The likelihood of detection of these bursts with pulsar timing and the Laser Interferometer Space Antenna ({\it LISA}) is estimated using several population models of SMBHs with masses $\gtrsim 10^7 {\rm M_\odot}$. Assuming a fraction of binaries $\ge 0.1$ in triple system, we find that few to few dozens of these bursts will produce residuals $>1$ ns, within the sensitivity range of forthcoming pulsar timing arrays (PTAs). However, most of such bursts will be washed out in the underlying confusion noise produced by all the other 'standard' SMBH binaries emitting in the same frequency window. A detailed data analysis study would be required to assess resolvability of such sources. Implementing a basic resolvability criterion, we find that the chance of catching a resolvable burst at a one nanosecond precision level is 2-50%, depending on the adopted SMBH evolution model. On the other hand, the probability of detecting bursts produced by massive binaries (masses $\gtrsim 10^7\msun$) with {\it LISA} is negligible.

Identifying the properties of the first generation of seeds of massive black holes is key to understanding the merger history and growth of galaxies. Mergers between ~100 solar mass seed black holes generate gravitational waves in the 0.1-10Hz band that lies between the sensitivity bands of existing ground-based detectors and the planned space-based gravitational wave detector, the Laser Interferometer Space Antenna (LISA). However, there are proposals for more advanced detectors that will bridge this gap, including the third generation ground-based Einstein Telescope and the space-based detector DECIGO. In this paper we demonstrate that such future detectors should be able to detect gravitational waves produced by the coalescence of the first generation of light seed black-hole binaries and provide information on the evolution of structure in that era. These observations will be complementary to those that LISA will make of subsequent mergers between more massive black holes. We compute the sensitivity of various future detectors to seed black-hole mergers, and use this to explore the number and properties of the events that each detector might see in three years of observation. For this calculation, we make use of galaxy merger trees and two different seed black hole mass distributions in order to construct the astrophysical population of events. We also consider the accuracy with which networks of future ground-based detectors will be able to measure the parameters of seed black hole mergers, in particular the luminosity distance to the source. We show that distance precisions of ~30% are achievable, which should be sufficient for us to say with confidence that the sources are at high redshift.

The properties of the first generation of black-hole seeds trace and distinguish different models of formation of cosmic structure in the high-redshift universe. The observational challenge lies in identifying black holes in the mass range ~100-1000 solar masses at redshift z~10. The typical frequencies of gravitational waves produced by the coalescence of the first generation of light seed black-hole binaries fall in the gap between the spectral ranges of low-frequency space-borne detectors (e.g., LISA) and high-frequency ground-based detectors (e.g., LIGO, Virgo and GEO 600). As such, these sources are targets for proposed third-generation ground-based instruments, such as the Einstein Telescope which is currently in design study. Using galaxy merger trees and four different models of black hole accretion - which are meant to illustrate the potential of this new type of source rather than to yield precise event-rate predictions - we find that such detectors could observe a few to a few tens of seed black-hole merger events in three years and provide, possibly unique, information on the evolution of structure in the corresponding era. We show further that a network of detectors may be able to measure the luminosity distance to sources to a precision of ~40%, allowing us to be confident of the high-redshift nature of the sources.

"Hard" massive black hole (MBH) binaries embedded in steep stellar cusps can shrink via three-body slingshot interactions. We show that this process will inevitably be accompanied by a burst of stellar tidal disruptions, at a rate that can be several orders of magnitude larger than that appropriate for a single MBH. Our numerical scattering experiments reveal that: 1) a significant fraction of stars initially bound to the primary hole are scattered into its tidal disruption loss cone by gravitational interactions with the secondary hole, an enhancement effect that is more pronounced for very unequal-mass binaries; 2) about 25% (40%) of all strongly interacting stars are tidally disrupted by a MBH binary of mass ratio q=1/81 (q=1/243) and eccentricity 0.1; and 3) two mechanisms dominate the fueling of the tidal disruption loss cone, a Kozai non-resonant interaction that causes the secular evolution of the stellar angular momentum in the field of the binary, and the effect of close encounters with the secondary hole that change the stellar orbital parameters in a chaotic way. For a hard MBH binary of 10^7 solar masses and mass ratio 0.01, embedded in an isothermal stellar cusp of velocity dispersion sigma*=100 km/s, the tidal disruption rate can be as large as 1/yr. This is 4 orders of magnitude higher than estimated for a single MBH fed by two-body relaxation. When applied to the case of a putative intermediate-mass black hole inspiraling onto Sgr A*, our results predict tidal disruption rates ~0.05-0.1/yr.

Massive black holes are key components of the assembly and evolution of cosmic structures and a number of surveys are currently on-going or planned to probe the demographics of these objects and to gain insight into the relevant physical processes. Pulsar Timing Arrays (PTAs) currently provide the only means to observe gravitational radiation from massive black hole binary systems with masses >10^7 solar masses. The whole cosmic population produces a stochastic background that could be detectable with upcoming Pulsar Timing Arrays. Sources sufficiently close and/or massive generate gravitational radiation that significantly exceeds the level of the background and could be individually resolved. We consider a wide range of massive black hole binary assembly scenarios, we investigate the distribution of the main physical parameters of the sources, such as masses and redshift, and explore the consequences for Pulsar Timing Arrays observations. Depending on the specific massive black hole population model, we estimate that on average at least one resolvable source produces timing residuals in the range ~5-50 ns. Pulsar Timing Arrays, and in particular the future Square Kilometre Array (SKA), can plausibly detect these unique systems, although the events are likely to be rare. These observations would naturally complement on the high-mass end of the massive black hole distribution function future surveys carried out by the Laser Interferometer Space Antenna (LISA)

All the physical processes involved in the formation, merging, and accretion history of massive black holes along the hierarchical build--up of cosmic structures are likely to leave an imprint on the gravitational waves detectable by future space--borne missions, such as LISA. We report here the results of recent studies, carried out by means of dedicated simulations of black hole build--up, aiming at understanding the impact on LISA observations of two ingredients that are crucial in every massive black hole formation scenario, namely: (i) the nature and abundance of the first black hole seeds and (ii) the large gravitational recoils following the merger of highly spinning black holes. We predict LISA detection rates spanning two order of magnitude, in the range 3-300 events per year, depending on the detail of the assumed massive black hole seed model. On the other hand, large recoil velocities do not dramatically compromise the efficiency of LISA observations. The number of detections may drop substantially (by ~60%), in scenarios characterized byabundant light seeds, but if seeds are already massive and/or relatively rare, the detection rate is basically unaffected.

The discovery of hypervelocity binary stars (HVBs) in the Galactic halo would provide definite evidence of the existence of a massive black hole companion to Sgr A*. Here we use an hybrid approach to compute the rate of ejection and the total number of HVBs produced by a hypothetical intermediate-mass black hole (IMBH, M_2<10^5\msun) orbiting Sgr A*. Depending on the mass of M_2 and on the properties of binary stars in the central parsec of the Milky Way, we show that the number of undisrupted HVBs expected to be expelled from the Galactic Center before binary black hole coalescence ranges from zero to a few dozens at most. Therefore, the non-detection of stellar binaries in a complete survey of hypervelocity stars would not rule out the occurrence of an IMBH-Sgr A* in-spiralling event within the last few 10^8 years.

Massive black hole binary systems, with masses in the range ~10^4-10^10 \msun, are among the primary sources of gravitational waves in the frequency window ~10^-9 Hz - 0.1 Hz. Pulsar Timing Arrays (PTAs) and the Laser Interferometer Space Antenna (LISA) are the observational means by which we will be able to observe gravitational radiation from these systems. We carry out a systematic study of the generation of the stochastic gravitational-wave background from the cosmic population of massive black hole binaries. We consider a wide variety of assembly scenarios and we estimate the range of signal strength in the frequency band accessible to PTAs. We show that, taking into account the uncertainties surrounding the actual key model parameters, the amplitude lies in the interval h_c(f = 10^-8 Hz)~5x10^-16 - 8x10^-15. The most optimistic predictions place the signal level at a factor of ~3 below the current sensitivity of Pulsar Timing Arrays, but within the detection range of the complete Parkes PTA for a wide variety of models, and of the future Square-Kilometer-Array PTA for all the models considered here. We also show that at frequencies >10^-8 Hz the frequency dependency of the generated background follows a power-law significantly steeper than f^-2/3, that has been considered so far. Finally we show that LISA observations of individual resolvable massive black hole binaries are complementary and orthogonal to PTA observations of a stochastic background from the whole population in the Universe. In fact, the detection of gravitational radiation in both frequency windows will enable us to fully characterise the cosmic history of massive black holes.

We consider a potentially new class of gravitational wave sources consisting of a white dwarf coalescing into a massive black hole in the mass range ~10^4-10^5\msun. These sources are of particular interest because the gravitational wave signal produced during the inspiral phase can be detected by the Laser Interferometer Space Antenna (LISA) and is promptly followed, in an extended portion of the black hole and white dwarf mass parameter space, by an electro-magnetic signal generated by the tidal disruption of the star, detectable with X-ray, optical and UV telescopes. This class of sources could therefore yield a considerable number of scientific payoffs, that include precise cosmography at low redshift, demographics of black holes in the mass range ~10^4-10^5\Msun, insights into dynamical interactions and populations of white dwarfs in the cores of dwarf galaxies, as well as a new probe into the structure and equation of state of white dwarfs. By modelling the gravitational and electromagnetic radiation produced by these events, we find them detectable in both observational windows at a distance ~200 Mpc, and possibly beyond for selected regions of the parameter space. We also estimate the detection rate for a number of model assumptions about black hole and white dwarf mass functions and dynamical interactions: the rate is (not surprisingly) highly uncertain, ranging from ~0.01 yr^-1 to ~100 yr^-1. This is due to the current limited theoretical understanding and minimal observational constraints for these objects and processes. However, capture rate scaling arguments favor the high end of the above range, making likely the detection of several events during the LISA lifetime.

We develop a formalism for studying the dynamics of massive black hole binaries embedded in gravitationally-bound stellar cusps, and study the binary orbital decay by three-body interactions, the impact of stellar slingshots on the density profile of the inner cusp, and the properties of the ejected hypervelocity stars (HVSs). We find that the scattering of bound stars shrinks the binary orbit and increases its eccentricity more effectively than that of unbound ambient stars. Binaries with initial eccentricities e>0.3 and/or unequal-mass companions (M_2/M_1<0.1) can decay by three-body interactions to the gravitational wave emission regime in less than a Hubble time. The stellar cusp is significantly eroded, and cores as shallow as \rho\propto r^-0.7 may develop from a pre-existing singular isothermal density profile. A population of HVSs is ejected in the host galaxy halo, with a total mass ~M_2. We scale our results to the scattering of stars bound to Sgr A*, the massive black hole in the Galactic Center, by an inspiraling companion of intermediate mass. Depending on binary mass ratio, eccentricity, and initial slope of the stellar cusp, a core of radius ~0.1 pc typically forms in 1-10 Myr. On this timescale about 500-2500 HVSs are expelled with speeds sufficiently large to escape the gravitational potential of the Milky Way.

Recent numerical simulations of coalescences of highly spinning massive black hole binaries (MBHBs) suggest that the remnant can suffer a recoil velocity of the order of few thousands km/s. We study here, by means of dedicated simulations of black holes build--up, how such extreme recoils could affect the cosmological coalescence rate of MBHBs, placing a robust lower limit for the predicted number of gravitational wave (GW) sources detectable by future space--borne missions (such as LISA). We consider two main routes for black hole formation: one where seeds are light remnants of Population III stars (~10^2 \msun), and one where seeds are much heavier (>~10^4 \msun), formed via the direct gas collapse in primordial nuclear disks. We find that extreme recoil velocities do not compromise the efficient MBHB detection by LISA. If seeds are already massive and/or relatively rare, the detection rate is reduced by only ~15%. The number of detections drops substantially (by ~60%) if seeds are instead light and abundant, but in this case the number of predicted coalescences is so high that at least ~10 sources in a three year observation are guaranteed.

Hypervelocity stars (HVSs) are a natural consequence of the presence of a massive nuclear black hole (Sgr A*) in the Galactic Center. Here we use the Brown et al. sample of unbound and bound HVSs together with numerical simulations of the propagation of HVSs in the Milky Way halo to constrain three plausible ejection mechanisms: 1) the scattering of stars bound to Sgr A* by an inspiraling intermediate-mass black hole (IMBH); 2) the disruption of stellar binaries in the tidal field of Sgr A*; and 3) the two-body scattering of stars off a cluster of stellar-mass black holes orbiting Sgr A*. We compare the predicted radial and velocity distributions of HVSs with the limited-statistics dataset currently available, and show that the IMBH model appears to produce a spectrum of ejection velocities that is too flat. Future astrometric and deep wide-field surveys of HVSs should shed unambiguous light on the stellar ejection mechanism and probe the Milky Way potential on scales as large as 200 kpc.

The formation, merging, and accretion history of massive black holes along the hierarchical build--up of cosmic structures leaves a unique imprint on the background of gravitational waves at mHz frequencies. We study here, by means of dedicated simulations of black hole build--up, the possibility of constraining different models of black hole cosmic evolution using future gravitational wave space--borne missions, such as LISA. We consider two main scenarios for black hole formation, namely, one where seeds are light (~10^2 \msun, remnant of Population III stars), and one where seeds are heavy (>~10^4 \msun, direct collapse). In all the models we have investigated, massive black hole binary coalescences do not produce a stochastic GW background, but rather, a set of individual resolved events. Detection of several hundreds merging events in 3 year LISA mission will be the sign of a heavy seed scenario with efficient formation of black hole seeds in a large fraction of high redshift halos. On the other extreme, a low event rate, about few tens in 3 years, is peculiar of scenarios where either the seeds are light, and many coalescences do not fall into the LISA band, or seeds are massive, but rare. In this case a decisive diagnostic is provided by the shape of the mass distribution of detected events. Light binaries (m<10^4\msun) are predicted in a fairly large number in Population III remnants models, but are totally absent in direct collapse models. Finally, a further, helpful diagnostic of black hole formation models lies in the distribution of the mass ratios in binary coalescences. While heavy seed models predict that most of the detected events involve equal mass binaries, in the case of light seeds, mass ratios are equally distributed in the range 0.1-1.

We study the long-term evolution of massive black hole binaries (MBHBs) at the centers of galaxies using detailed scattering experiments to solve the full three-body problem. Ambient stars drawn from a isotropic Maxwellian distribution unbound to the binary are ejected by the gravitational slingshot. We construct a minimal, hybrid model for the depletion of the loss cone and the orbital decay of the binary, and show that secondary slingshots - stars returning on small impact parameter orbits to have a second super-elastic scattering with the MBHB - may considerably help the shrinking of the pair in the case of large binary mass ratios. In the absence of loss-cone refilling by two-body relaxation or other processes, the mass ejected before the stalling of a MBHB is half the binary reduced mass. About 50% of the ejected stars are expelled ejected in a "burst" lasting ~1E4 yrs M_6^1/4, where M_6 is the binary mass in units of 1E6 Msun. The loss cone is completely emptied in a few bulge crossing timescales, 1E7 yrs M_6^1/4. Even in the absence of two-body relaxation or gas dynamical processes, unequal mass and/or eccentric binaries with M_6 >0.1 can shrink to the gravitational wave emission regime in less than a Hubble time, and are therefore "safe" targets for the planned Laser Interferometer Space Antenna (LISA).

We explore the nature of possible electromagnetic counterparts of coalescences of massive black hole binaries at cosmological distances detectable by the Laser Interferometer Space Antenna (LISA). An electromagnetic precursor, during the last year of gravitational wave (GW)-driven inspiral, or an afterglow within few years after coalescence, may highlight the position in the sky of galaxies hosting LISA sources. We show that observations of precursors and afterglows are mutually exclusive, depending on the mass of the primary black hole. Precursors are expected to occur in binaries where the primary (more massive) black hole is heavier than ~10^7 \Msun. They may correspond to on-off states of accretion, i.e., to a bright X-ray source decaying into quiescence before black hole coalescence, and are likely associated to disturbed galaxies showing signs of ongoing starbursts. Coalescences of lighter binaries, with masses <5x10^6 \Msun, lack of any precursor, as gas is expected to be consumed long before the GW-driven orbital decay. Such events would not be hosted by (massive) galaxies with an associated starburst, given the slow binary inspiral time compared to the typical time scale of starbursts. By contrast, coalescence, for such light binaries, is followed by an electromagnetic afterglow, i.e., an off-on accretion state rising in <20 yrs. Using a cosmological merger tree algorithm, we show that future X-ray missions such as XEUS will be able to identify, in 20 yrs operation, almost all the massive BH binary detectable by LISA, and, in only 5 yrs, all the LISA sources at z>6.

We use full three-body scattering experiments to study the ejection of hypervelocity stars (HVSs) by massive black hole binaries (MBHBs) at the center of galaxies. Ambient stars drawn from a Maxwellian distribution unbound to the binary are expelled by the gravitational slingshot. Accurate measurements of thermally averaged hardening, mass ejection, and eccentricity growth rates (H, J, and K) for MBHBs in a fixed stellar background are obtained by numerical orbit integration from initial conditions determined by Monte Carlo techniques. Three-body interactions create a subpopulation of HVSs on nearly radial, corotating orbits, with a spatial distribution that is initially highly flattened in the inspiral plane of the MBHB, but becomes more isotropic with decreasing binary separation. The degree of anisotropy is smaller for unequal mass binaries and larger for stars with higher kick velocities. Eccentric MBHBs produce a more prominent tail of high-velocity stars and break planar symmetry, ejecting HVSs along a broad jet perpendicular to the semimajor axis. The jet two-sidedness decreases with increasing binary mass ratio, while the jet opening-angle increases with decreasing kick velocity and orbital separation. The detection of a numerous population of HVSs in the halo of the Milky Way by the next generation of large astrometric surveys like GAIA may provide a unique signature of the history, nature, and environment of the MBH at the Galactic center.

We study the long-term evolution of massive black hole binaries (MBHBs) at the centers of galaxies using detailed full three-body scattering experiments. Stars, drawn from a distribution unbound to the binary, are ejected by the gravitational slingshot. We quantify the effect of secondary slingshots -- stars returning on small impact parameter orbits to have a second super-elastic scattering with the MBHB -- on binary separation. Even in the absence of two-body relaxation or gas dynamical processes, very unequal mass binaries of mass M=10^7 solar masses can shrink to the gravitational wave emission regime in less than a Hubble time, and are therefore a target for the planned Laser Interferometer Space Antenna (LISA). Three-body interactions create a subpopulation of hypervelocity stars on nearly radial, corotating orbits, with a spatial distribution that is initially highly flattened in the inspiral plane of the MBHB, but becomes more isotropic with decreasing binary separation. The mass ejected is ~0.7 times the binary reduced mass, and most of the stars are ejected in an initial burst lasting much less than a bulge crossing time.

We compute the expected low-frequency gravitational wave signal from coalescing massive black hole (MBH) binaries at the center of galaxies. We follow the merging history of halos and associated holes via cosmological Monte Carlo realizations of the merger hierarchy from early times to the present in a LCDM cosmology. MBHs get incorporated through a series of mergers into larger and larger halos, sink to the centre owing to dynamical friction, accrete a fraction of the gas in the merger remnant to become more massive, and form a binary system. Stellar dynamical processes dominates the orbital evolution of the binary at large separations, while gravitational wave emission takes over at small radii, causing the final coalescence of the system. We discuss the observability of inspiraling MBH binaries by a low-frequency gravitational wave experiment such as the planned Laser Interferometer Space Antenna (LISA), discriminating between resolvable sources and unresolved confusion noise. Over a 3-year observing period LISA should resolve this GWB into discrete sources, detecting ~90 individual events above a S/N=5 confidence level, while expected confusion noise is well below planned LISA capabilities.

Massive black hole binaries, with masses in the range 1E3-1E8 Msun, are expected to be the most powerful sources of gravitational radiation at mHz frequencies, and hence are among the primary targets for the planned Laser Interferometer Space Antenna (LISA). We extend and refine our previous analysis (Sesana et al. 2004), detailing the gravitational wave signal expected from a cosmological population of massive black hole binaries. As done in our previous paper, we follow the merger history of dark matter halos, the dynamics of the massive black holes they host, and their growth via gas accretion and binary coalescences in a LCDM cosmology. Stellar dynamical processes dominates the orbital evolution of black hole binaries at large separations, while gravitational wave emission takes over at small radii, causing the final coalescence of the pairs. We show that the GW signal from this population, in a 3 year LISA observation, will be resolved into approx 90 discrete events with S/N>5, among which approx. 35 will be observed above threshold until coalescence. These "merging events" involve relatively massive binaries, M=10E5 Msun, in the redshift range 2<z<6. The remaining approx. 55 events come from higher redshift, less massive binaries (M=5E3 Msun at z>6) and, although their S/N integrated over the duration of the observation can be substantial, the final coalescence phase is at too high frequency to be directly observable by space-based interferometers such as LISA. LISA will be able to detect a fraction approx. 90% of all the coalescences of massive black hole binaries occurring at z=5. The residual confusion noise from unresolved massive black hole binaries is expected to be at least an order of magnitude below the estimated stochastic noise.

We compute the expected gravitational wave signal from coalescing massive black hole (MBH) binaries at the center of galaxies in a hierarchical structure formation scenario in which seed holes of intermediate mass form far up in the dark halo merger tree. The merger history of DM halos and MBHs is followed from z=20 to the present in a LCDM cosmology. MBHs get incorporated through halo mergers into larger and larger structures, sink to the center owing to dynamical friction against the DM background, accrete cold material in the merger remnant, and form MBH binary systems. Stellar dynamical interactions cause the hardening of the binary at large separations, while gravitational wave emission takes over at small radii and leads to the final coalescence of the pair. The integrated emission from inspiraling MBH binaries results in a gravitational wave background (GWB). The characteristic strain spectrum has the standard h_c(f)\propto f^{-2/3} behavior only in the range 1E-9<f<1E-6 Hz. At lower frequencies the orbital decay of MBH binaries is driven by the ejection of background stars, and h_c(f) \propto f. At higher frequencies, f>1E-6 Hz, the strain amplitude is shaped by the convolution of last stable circular orbit emission. We discuss the observability of inspiraling MBH binaries by the planned LISA. Over a 3-year observing period LISA should resolve this GWB into discrete sources, detecting ~60 (~250) individual events above a S/N=5 (S/N=1) confidence level. (Abridged)

We investigate a hierarchical structure formation scenario describing the evolution of a Super Massive Black Holes (SMBHs) population. The seeds of the local SMBHs are assumed to be 'pregalactic' black holes, remnants of the first POPIII stars. As these pregalactic holes become incorporated through a series of mergers into larger and larger halos, they sink to the center owing to dynamical friction, accrete a fraction of the gas in the merger remnant to become supermassive, form a binary system, and eventually coalesce. A simple model in which the damage done to a stellar cusps by decaying BH pairs is cumulative is able to reproduce the observed scaling relation between galaxy luminosity and core size. An accretion model connecting quasar activity with major mergers and the observed BH mass-velocity dispersion correlation reproduces remarkably well the observed luminosity function of optically-selected quasars in the redshift range 1<z<5. We finally asses the potential observability of the gravitational wave background generated by the cosmic evolution of SMBH binaries by the planned space-born interferometer LISA.