Search

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., Comment: 12 pages, 3 figures, 5 tables. These proceedings are a summary of the results presented in arXiv:1703.09722

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., 12 pages, 9 figures, published in MNRAS

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., 5 figures, 1 table

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., Accepted for publication by the Monthly Notices of the Royal Astronomical Society

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

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

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

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., Comment: 14 pages, 11 figures, 3 tables. Resubmitted to the Astrophysical Journal after some major revision of the results including a proper estimate of the intrinsic mass of the binary candidates

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), 12 pages, 10 figures, accepted for publication in MNRAS. Results revised (differences within a factor of two) after a bug in the code for generating the timing residuals has been fixed

The black hole binary (BHB) coalescence rates inferred from the Advanced LIGO detection of GW150914 imply an unexpectedly loud gravitational-wave (GW) sky at millihertz frequencies accessible to the Evolved Laser Interferometer Space Antenna (eLISA), with several outstanding consequences. First, up to thousands of BHBs will be individually resolvable by eLISA; second, millions of nonresolvable BHBs will build a confusion noise detectable with a signal-to-noise ratio of a few to hundreds; third - and perhaps most importantly - up to hundreds of BHBs individually resolvable by eLISA will coalesce in the Advanced LIGO band within 10 y. eLISA observations will tell Advanced LIGO and all electromagnetic probes weeks in advance when and where these BHB coalescences will occur, with uncertainties of

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., Comment: 14 pages, 8 figures, 1 table. Updated to match published version

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., Comment: 17 pages, 25 figures, accepted for publication in MNRAS

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., Comment: 3 figures; small changes to match version published in Physical Review Letters

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., Comment: JCAP style, 48 pages, 16 figures, 10 tables; v2: minor addition to match published version; v3: acknowledgements added

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., 28 pages, 13 figures, 7 tables

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., 6 pages 4 figures, submitted to MNRAS letters

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., 5 pages, 5 figures

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., Accepted by MNRAS Letters, 6 pages, 3 figures

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, Comment: 19 pages, 11 figures, accepted for publication in MNRAS

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., 6 pages, 2 figures, 1 table. Accepted for publication in Physical Review Letters

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.

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., Comment: 22 pages, 3 figures. Minor change (added two references) to match version accepted in the Proceedings of LISA Symposium X (Journal of Physics: Conference Series)

The parallel session C5 on Space-Based Detectors gave a broad overview over the planned space missions related to gravitational wave detection. Overviews of the revolutionary science to be expected from LISA was given by Alberto Sesana and Sasha Buchman. The launch of LISA Pathfinder (LPF) is planned for 2015. This mission and its payload “LISA Technology Package” will demonstrate key technologies for LISA. In this context, reference masses in free fall for LISA, and gravitational physics in general, was described by William Weber, laser interferometry at the pico-metre level and the optical bench of LPF was presented by Christian Killow and the performance of the LPF optical metrology system by Paul McNamara. While LPF will not yet be sensitive to gravitational waves, it may nevertheless be used to explore fundamental physics questions, which was discussed by Michele Armano. Some parts of the LISA technology that are not going to be demonstrated by LPF, but under intensive development at the moment, were presented by Oliver Jennrich and Oliver Gerberding. Looking into the future, Japan is studying the design of a mid-frequency detector called DECIGO, which was discussed by Tomotada Akutsu. Using atom interferometry for gravitational wave detection has also been recently proposed, and it was critically reviewed by Peter Bender. In the nearer future, the launch of GRACE Follow-On (for Earth gravity observation) is scheduled for 2017, and it will include a Laser Ranging Interferometer as technology demonstrator. This will be the first inter-spacecraft laser interferometer and has many aspects in common with the LISA long arm, as discussed by Andrew Sutton.

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., Comment: 26 pages, 18 figures, 6 tables. Replaced to match the final version accepted for publication in The Astrophysical Journal, after major revisions following the referee's suggestions

International audience; This paper summarizes parallel session C1: Pulsar Timing Arrays of the Amaldi10/GR20 Meeting held in Warsaw, Poland in July 2013. The session showcased recent results from pulsar timing array collaborations, advances in modelling the gravitational-wave signal, and new methods to search for and characterize gravitational waves in pulsar timing array observations.

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., 16 pages, 8 figures, 1 table, accepted for publication in MNRAS

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., 10 pages, 1 figure, 1 table, prepared for the Proceedings of the 26th Texas Symposium on Relativistic Astrophysics, Sao Paulo, December 16-20 2012

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., 12 pages, 4 figures, 1 table; submitted to Physical Review D

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., Comment: 14 pages, 2 Tabes, 11 Figures, submitted to MNRAS, comments welcome

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., Comment: 26 pages, 9 figures, invited article for the focus issue on astrophysical black holes in Classical and Quantum Gravity, guest editors: D. Merritt and L. Rezzolla. Submitted

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., 13 pages, 12 Figures, Astronomy and Astrophysics accepted

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., Comment: 20 pages, 8 figures, proceedings of the 9th Amaldi Conference on Gravitational Waves. Final journal version. For a longer exposition of the eLISA science case, see http://arxiv.org/abs/1201.3621

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., Comment: 5 pages, 3 figures, 2 tables. Minor changes to address comments by the referees

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., Based on the talk given at GR19. 10 pages, 2 figures. Revised version, accepted in CQG

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., Comment: Submitted to ApJ. v1 original submission; v2 minor changes in response to referee

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., Comment: 9 pages, 1 figure, accepted by MNRAS. Comments are welcome

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., Comment: 14 pages, 11 figures, submitted to Phys. Rev. D

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., 12 pages, 3 figures, submitted to Class. Quantum Grav. for proceedings of 8th LISA Symposium; v2 minor changes for consistency with accepted version

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), Comment: 20 pages, 12 figures, 2 color figures, submitted to Phys. Rev. D

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., accepted by Classical and Quantum Gravity. Review talk for the Amaldi8 conference series

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., Comment: 24 pages, 16 figures, submitted to Phys. Rev. D

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., Comment: 15 pages, 9 figures, accepted for publication in the Astrophysical Journal

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., Comment: 14 pages, 6 figures, 2 tables, accepted for proceedings of 13th GWDAW meeting

"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., Comment: 5 pages, 3 figures, accepted for publication in the Astrophysical Journal Letters

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., Accepted for publication by MNRAS, minor changes

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., Comment: 5 pages, 2 figures, 1 table, accepted to ApJ letters; v2 contains more technical details in response to referee's comments, 1 new figure, table removed

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, 5 pages, 3 figures, accepted for publication in MNRAS letters