Massive Black Hole Mergers Research Articles

The gravitational wave background from massive black hole binaries in Illustris: spectral features and time to detection with pulsar timing arrays

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.

Reconstructing the dark sector interaction with LISA

We perform a forecast analysis of the ability of the LISA space-based interferometer to reconstruct the dark sector interaction using gravitational wave standard sirens at high redshift. We employ Gaussian process methods to reconstruct the distance-redshift relation in a model independent way. We adopt simulated catalogues of standard sirens given by merging massive black hole binaries visible by LISA, with an electromagnetic counterpart detectable by future telescopes. The catalogues are based on three different astrophysical scenarios for the evolution of massive black hole mergers based on the semi-analytic model of E. Barausse, Mon. Not. Roy. Astron. Soc. 423 (2012) 2533. We first use these standard siren datasets to assess the potential of LISA in reconstructing a possible interaction between vacuum dark energy and dark matter. Then we combine the LISA cosmological data with supernovae data simulated for the Dark Energy Survey. We consider two scenarios distinguished by the time duration of the LISA mission: 5 and 10 years. Using only LISA standard siren data, the dark sector interaction can be well reconstructed from redshift z∼1 to z∼3 (for a 5 years mission) and z∼1 up to z∼5 (for a 10 years mission), though the reconstruction is inefficient at lower redshift. When combined with the DES datasets, the interaction is well reconstructed in the whole redshift region from 0z∼ to z∼3 (5 yr) and z∼0 to z∼5 (10 yr), respectively. Massive black hole binary standard sirens can thus be used to constrain the dark sector interaction at redshift ranges not reachable by usual supernovae datasets which probe only the z≲ 1.5 range. Gravitational wave standard sirens will not only constitute a complementary and alternative way, with respect to familiar electromagnetic observations, to probe the cosmic expansion, but will also provide new tests to constrain possible deviations from the standard ΛCDM dynamics, especially at high redshift.

Constraining early and interacting dark energy with gravitational wave standard sirens: the potential of the eLISA mission

We perform a forecast analysis of the capability of the eLISA space-based interferometer to constrain models of early and interacting dark energy using gravitational wave standard sirens. We employ simulated catalogues of standard sirens given by merging massive black hole binaries visible by eLISA, with an electromagnetic counterpart detectable by future telescopes. We consider three-arms mission designs with arm length of 1, 2 and 5 million km, 5 years of mission duration and the best-level low frequency noise as recently tested by the LISA Pathfinder. Standard sirens with eLISA give access to an intermediate range of redshift 1 ≲ z ≲ 8, and can therefore provide competitive constraints on models where the onset of the deviation from ΛCDM (i.e. the epoch when early dark energy starts to be non-negligible, or when the interaction with dark matter begins) occurs relatively late, at z ≲ 6. If instead early or interacting dark energy is relevant already in the pre-recombination era, current cosmological probes (especially the cosmic microwave background) are more efficient than eLISA in constraining these models, except possibly in the interacting dark energy model if the energy exchange is proportional to the energy density of dark energy.

When and where did GW150914 form?

Abstract The recent LIGO detection of gravitational waves (GW150914), likely originating from the merger of two ∼ 30M⊙ black holes suggests progenitor stars of low metallicity ([Z/Z⊙] ≲ 0.3), constraining when and where the progenitor of GW150914 may have formed. We combine estimates of galaxy properties (star-forming gas metallicity, star formation rate and merger rate) across cosmic time to predict the low redshift black hole – black hole merger rate as a function of present day host galaxy mass, Mgal, the formation redshift of the progenitor system zf and different progenitor metallicities Zp. For $Z_\mathrm{p}\geqslant 0.1 \rm Z_{{\odot }}$, the signal is dominated by binaries in massive galaxies with zf ≃ 2 while below $Z_\mathrm{p}\leqslant 0.1\rm Z_{{\odot }}$ most mergers come from binaries formed around zf ≃ 0.5 in dwarf galaxies. Additional gravitational wave detections from merging massive black holes will provide constraints on the mass–metallicity relation and massive star formation at high redshifts.

A new route towards merging massive black holes

Recent advances in gravitational-wave astronomy make the direct detection of gravitational waves from the merger of two stellar-mass compact objects a realistic prospect. Evolutionary scenarios towards mergers of double compact objects generally invoke common-envelope evolution which is poorly understood, leading to large uncertainties in merger rates. We explore the alternative scenario of massive overcontact binary (MOB) evolution, which involves two very massive stars in a very tight binary which remain fully mixed due to their tidally induced high spin. We use the public stellar-evolution code MESA to systematically study this channel by means of detailed simulations. We find that, at low metallicity, MOBs produce double-black-hole (BH+BH) systems that will merge within a Hubble time with mass ratios close to one, in two mass ranges, ~25...60msun and >~ 130msun, with pair instability supernovae (PISNe) being produced in-between. Our models are also able to reproduce counterparts of various stages in the MOB scenario in the local Universe, providing direct support for it. We map the initial parameter space that produces BH+BH mergers, determine the expected chirp mass distribution, merger times, Kerr parameters and predict event rates. We typically find that for Z~<Z_sun/10, there is one BH+BH merger for ~1000 core-collapse supernovae. The advanced LIGO (aLIGO) detection rate is more uncertain and depends on the metallicity evolution. Deriving upper and lower limits from a local and a global approximation for the metallicity distribution of massive stars, we estimate aLIGO detection rates (at design limit) of ~19-550 yr^(-1) for BH+BH mergers below the PISN gap and of ~2.1-370 yr^(-1) above the PISN gap. Even with conservative assumptions, we find that aLIGO should soon detect BH+BH mergers from the MOB scenario and that these could be the dominant source for aLIGO detections.

Science with the space-based interferometer eLISA: Supermassive black hole binaries

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.

Astrophysical constraints on massive black hole binary evolution from pulsar timing arrays

Abstract We consider the information that can be derived about massive black hole binary (MBHB) 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 MBHBs with chirp mass in the range 106–1011 M⊙ evolving solely due to radiation reaction. Our parametrized 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 for any particular mass.

The development of ground based gravitational wave astronomy and opportunities for Australia–China collaboration

This paper begins by reviewing the development of gravitational wave astronomy from the first predictions of gravitational waves to development of technologies across the entire gravitational wave spectrum, and then focuses on the current status of ground based gravitational wave detectors. With substantial improvements already demonstrated in early commissioning it is emphasised that Advanced detectors are on track for first detection of gravitational waves. The importance of a worldwide array of detectors is emphasised, and recent results are shown that demonstrate the continued advantage of a southern hemisphere detector. Finally it is shown that a north–south pair of 8 km arm length detectors would give rise to a dramatic improvement in event rate, enabling a pair of detectors to encompass a 64-times larger volume of the universe, to conduct a census on all stellar mass black hole mergers to [Formula: see text] and to observe neutron star mergers to a distance of [Formula: see text][Formula: see text]800 Mpc.

THE FASTEST UNBOUND STARS IN THE UNIVERSE

The discovery of hypervelocity stars (HVSs) leaving our galaxy with speeds of nearly 103 km s−1 has provided strong evidence of the existence of a massive compact object at the galaxy's center. HVSs ejected via the disruption of stellar binaries can occasionally yield a star with km s−1; here we show that this mechanism can be extended to massive black hole (MBH) mergers, where the secondary star is replaced by a MBH with mass . We find that stars that are originally bound to the secondary MBH are frequently ejected with km s−1, and occasionally with velocities ∼105 km s−1 (one third the speed of light). For this reason we refer to stars ejected from these systems as "semi-relativistic" hypervelocity stars (SHSs). Bound to no galaxy, the velocities of these stars are so great that they can cross a significant fraction of the observable universe in the time since their ejection (several Gpc). We demonstrate that if a significant fraction of MBH mergers undergoes a phase in which their orbital eccentricity is ≳0.5 and their periapse distance is tens of the primary's Schwarzschild radius, the space density of fast-moving ( km s−1) SHSs may be as large as 103 Mpc−3. Hundreds of SHSs will be giant stars that can be detected by future all-sky infrared surveys such as WFIRST or Euclid and proper motion surveys such as LSST, with spectroscopic follow-up being possible with the James Webb Space Telescope.

UNDERSTANDING BLACK HOLE MASS ASSEMBLY VIA ACCRETION AND MERGERS AT LATE TIMES IN COSMOLOGICAL SIMULATIONS

Accretion is thought to primarily contribute to the mass accumulation history of supermassive black holes throughout cosmic time. While this may be true at high redshifts, at lower redshifts and for the most massive black holes mergers themselves might add significantly to the mass budget. We evolve SMBHs from $4 > z > 0$ using merger trees derived from hydrodynamical cosmological simulations of a cluster and void region, scaled to the observed value of the stellar mass fraction to account for overcooling. Mass gains from gas accretion proportional to bulge growth and BH-BH mergers are tracked, as are black holes that remain "orbiting" due to insufficient dynamical friction in a merger remnant, as well as those that are ejected due to gravitational recoil. We find that gas accretion remains the dominant source of mass accumulation in almost all SMBHs; mergers contribute $2.5\pm0.1\%$ for all SMBHs in the cluster and $1.0\pm0.1\%$ in the void since $z = 4$. However, mergers are significant for massive SMBHs. The fraction of mass accumulated from mergers for central BHs generally increases for larger values of the host bulge mass: in the void, the fraction is $2\%$ at $M_{*, bul} = 10^{10} M_{\odot}$, increasing to $4\%$ at $M_{*, bul} \gtrsim 10^{11} M_{\odot}$, and in the cluster it is $4\%$ at $M_{*, bul} = 10^{10} M_{\odot}$ and $23\%$ at $10^{12} M_{\odot}$. We find that $40\%$ of SMBHs and $\approx 8\%$ of the total SMBH mass is found orbiting in the cluster region at $z = 0$. The existence of orbiting and ejected SMBHs requires modification of the Soltan argument. We estimate this correction to the integrated accreted mass density of SMBHs to be in the range $6-21\%$, with a mean value of $11\pm3\%$. We also calculate the total energy output and strain from gravitational waves emitted by merging SMBHs, and obtain a signal potentially detectable by pulsar timing arrays.

THE FORMATION AND GRAVITATIONAL-WAVE DETECTION OF MASSIVE STELLAR BLACK HOLE BINARIES

If binaries consisting of two 100 Msun black holes exist they would serve as extraordinarily powerful gravitational-wave sources, detectable to redshifts of z=2 with the advanced LIGO/Virgo ground-based detectors. Large uncertainties about the evolution of massive stars preclude definitive rate predictions for mergers of these massive black holes. We show that rates as high as hundreds of detections per year, or as low as no detections whatsoever, are both possible. It was thought that the only way to produce these massive binaries was via dynamical interactions in dense stellar systems. This view has been challenged by the recent discovery of several stars with mass above 150 Msun in the R136 region of the Large Magellanic Cloud. Current models predict that when stars of this mass leave the main sequence, their expansion is insufficient to allow common envelope evolution to efficiently reduce the orbital separation. The resulting black-hole--black-hole binary remains too wide to be able to coalesce within a Hubble time. If this assessment is correct, isolated very massive binaries do not evolve to be gravitational-wave sources. However, other formation channels exist. For example, the high multiplicity of massive stars, and their common formation in relatively dense stellar associations, opens up dynamical channels for massive black hole mergers (e.g., via Kozai cycles or repeated binary-single interactions). We identify key physical factors that shape the population of very massive black-hole--black-hole binaries. Advanced gravitational-wave detectors will provide important constraints on the formation and evolution of very massive stars.

Ripple effects and oscillations in the broad Fe Kα line as a probe of massive black hole mergers

When a sufficiently massive satellite (or secondary) black hole is embedded in a gas disc around a (primary) supermassive black hole, it can open an empty gap in the disc. A gap-opening secondary close to the primary will leave an imprint in the broad component of the Fe Kα emission line, which varies in a unique and predictable manner. If the gap persists into the innermost disc, the effect consists of a pair of dips in the broad line which ripple blue-ward and red-ward from the line centroid energy, respectively, as the gap moves closer to the primary. This ripple effect could be unambiguously detectable and allow an electromagnetic monitoring of massive black hole mergers as they occur. As the mass ratio of the secondary to primary black hole increases to q ≳ 0.01, we expect the gap to widen, possibly clearing a central cavity in the inner disc, which shows up in the broad Fe Kα line component. If the secondary stalls at ≥ 102rg in its in-migration, due to low corotating gas mass, a detectable ripple effect occurs in the broad line component on the disc viscous time-scale as the inner disc drains and the outer disc is dammed. If the secondary maintains an accretion disc within a central cavity, due to dam bursting or leakage, a periodic ‘see-saw’ oscillation effect is exhibited in the observed line profile. Here, we demonstrate the range of ripple effect signatures potentially detectable with Astro-H and IXO/Athena, and oscillation effects potentially detectable with XMM–Newton or LOFT for a wide variety of merger and disc conditions, including gap width (or cavity size), disc inclination angle and emissivity profile, damming of the accretion flow by the secondary, and a minidisc around the satellite black hole. A systematic study of ripple effects would require a telescope effective area substantially larger than that planned for IXO/Athena. Future mission planning should take this into account. Observations of the ripple effect and periodic oscillations can be used to provide an early warning of gravitational radiation emission from the AGN. Once gravitational waves consistent with massive black hole mergers are detected, an archival search for the Fe Kα ripple effect or periodic oscillations will help in localizing their origin.

Gravitational Waves and Time-Domain Astronomy

AbstractThe gravitational-wave window onto the universe will open in roughly five years, when Advanced LIGO and Virgo achieve the first detections of high-frequency gravitational waves, most likely coming from compact binary mergers. Electromagnetic follow-up of these triggers, using radio, optical, and high energy telescopes, promises exciting opportunities in multi-messenger time-domain astronomy. In the decade, space-based observations of low-frequency gravitational waves from massive black hole mergers, and their electromagnetic counterparts, will open up further vistas for discovery. This two-part workshop featured brief presentations and stimulating discussions on the challenges and opportunities presented by gravitational-wave astronomy. Highlights from the workshop, with the emphasis on strategies for electromagnetic follow-up, are presented in this report.

Sky localization of complete inspiral-merger-ringdown signals for nonspinning massive black hole binaries

We investigate the capability of LISA to measure the sky position of equal-mass, nonspinning black hole binaries, combining for the first time the entire inspiral-merger-ringdown signal, the effect of the LISA orbits, and the complete three-channel LISA response. We consider an ensemble of systems near the peak of LISA's sensitivity band, with total rest mass of $2\ifmmode\times\else\texttimes\fi{}{10}^{6}\text{ }\text{ }{\mathrm{M}}_{\ensuremath{\bigodot}}$, a redshift of $z=1$, and randomly chosen orientations and sky positions. We find median sky localization errors of approximately $\ensuremath{\sim}3$ arcminutes. This is comparable to the field of view of powerful electromagnetic telescopes, such as the James Webb Space Telescope, that could be used to search for electromagnetic signals associated with merging massive black holes. We investigate the way in which parameter errors decrease with measurement time, focusing specifically on the additional information provided during the merger-ringdown segment of the signal. We find that this information improves all parameter estimates directly, rather than through diminishing correlations with any subset of well-determined parameters. Although we have employed the baseline LISA design for this study, many of our conclusions regarding the information provided by mergers will be applicable to alternative mission designs as well.

Resonant trapping of stars by merging massive black hole binaries

A massive black hole binary might resonantly trap a star (e.g. a white dwarf) and the gas released by its tidal disruption might emit electromagnetic wave signals around the coalescence of the binary. With post-Newtonian equations of motion including gravitational radiation reaction, we numerically studied resonant trappings by black hole binaries with mass ratio 1/100. It is found that 2:1 (and simultaneously 4:2) mean motion resonances of the binaries would be strong and could, in principle, draw small third objects deep into relativistic regimes (e.g. ~10 Schwarzschild radii). The inclinations of the trapped objects could increase significantly and, in some cases, retrograde orbits could be realized eventually.

The stellar kinematic signature of massive black hole binaries

The stalling radius of a merging massive binary black hole (BBH) is expected to be below 0".1 even in nearby galaxies (Yu 2002), and thus BBHs are not expected to be spatially resolved in the near future. However, as we show below, a BBH may be detectable through the significantly anisotropic stellar velocity distribution it produces on scales 5-10 times larger than the binary separation. We calculate the velocity distribution of stable orbits near a BBH by solving the restricted three body problem for a BBH embedded in a bulge potential. We present high resolution maps of the projected velocity distribution moments, based on snapshots of ~ 10^8 stable orbits. The kinematic signature of a BBH in the average velocity maps is a counter rotating torus of stars outside the BBH Hill spheres. The velocity dispersion maps reveal a dip in the inner region, and an excess of 20-40% further out, compared to a single BH of the same total mass. More pronounced signatures are seen in the third and fourth Gauss-Hermite velocity moments maps. The detection of these signatures may indicate the presence of a BBH currently, or at some earlier time, which depends on the rate of velocity phase space mixing following the BBH merger.

TIME-DEPENDENT MODELS FOR THE AFTERGLOWS OF MASSIVE BLACK HOLE MERGERS

The Laser Interferometer Space Antenna (LISA) will detect gravitational wave signals from coalescing pairs of massive black holes in the total mass range (10^5 - 10^7)/Msol out to cosmological distances. Identifying and monitoring the electromagnetic counterparts of these events would enable cosmological studies and offer new probes of gas physics around well-characterized massive black holes. Milosavljevic & Phinney (2005) proposed that a circumbinary disk around a binary of mass ~10^6 Msol will emit an accretion-powered X-ray afterglow approximately one decade after the gravitational wave event. We revisit this scenario by using Green's function solutions to calculate the temporal viscous evolution and the corresponding electromagnetic signature of the circumbinary disk. Our calculations suggest that an electromagnetic counterpart may become observable as a rapidly brightening source soon after the merger, i.e. several years earlier than previously thought. The afterglow can reach super-Eddington luminosities without violating the local Eddington flux limit. It is emitted in the soft X-ray by the innermost circumbinary disk, but it may be partially reprocessed at optical and infrared frequencies. We also find that the spreading disk becomes increasingly geometrically thick close to the central object as it evolves, indicating that the innermost flow could become advective and radiatively inefficient, and generate a powerful outflow. We conclude that the mergers of massive black holes detected by LISA offer unique opportunities for monitoring on humanly tractable timescales the viscous evolution of accretion flows and the emergence of outflows around massive black holes with precisely known masses, spins and orientations.

Ring formation from an oscillating black hole

Massive black hole (BH) mergers can result in the merger remnant receiving a "kick", of order 200 km s$^{-1}$ or more, which will cause the remnant to oscillate about the galaxy centre. Here we analyze the case where the BH oscillates through the galaxy centre perpendicular or parallel to the plane of the galaxy for a model galaxy consisting of an exponential disk, a Plummer model bulge, and an isothermal dark matter halo. For the perpendicular motion we find that there is a strong resonant forcing of the disk radial motion near but somewhat less than the "resonant radii" $r_R$ where the BH oscillation frequency is equal one-half, one-fourth, (1/6, etc.) of the radial epicyclic frequency in the plane of the disk. Near the resonant radii there can be a strong enhancement of the radial flow and disk density which can lead to shock formation. In turn the shock may trigger the formation of a ring of stars near $r_R$. As an example, for a BH mass of $10^8 M_\odot$ and a kick velocity of 150 km s$^{-1}$, we find that the resonant radii lie between 0.2 and 1 kpc. For BH motion parallel to the plane of the galaxy we find that the BH leaves behind it a supersonic wake where star formation may be triggered. The shape of the wake is calculated as well as the slow-down time of the BH. The differential rotation of the disk stretches the wake into ring-like segments.

Gravitational recoil: effects on massive black hole occupation fraction over cosmic time

We assess the influence of massive black hole (MBH) ejections from galaxy centres due to gravitational radiation recoil, along the cosmic merger history of the MBH population. We discuss the ‘danger’ of recoil for MBHs as a function of different MBH spin-orbit configurations and of the host halo cosmic bias, and on how that reflects on the occupation fraction of MBHs. We assess ejection probabilities for mergers occurring in a gas-poor environment, in which the MBH binary coalescence is driven by stellar dynamical processes and the spin-orbit configuration is expected to be isotropically distributed. We contrast this case with the ‘aligned’ case. The latter is the more realistic situation for gas-rich, i.e. ‘wet’, mergers, which are expected for high-redshift galaxies. We find that if all haloes at z > 5–7 host an MBH, the probability of the Milky Way (or similar size galaxy) to host an MBH today is less than 50 per cent, unless MBHs form continuously in galaxies. The occupation fraction of MBHs, intimately related to halo bias and MBH formation efficiency, plays a crucial role in increasing the retention fraction. Small haloes, with shallow potential wells and low escape velocities, have a high ejection probability, but the MBH merger rate is very low along their galaxy formation merger hierarchy: MBH formation processes are likely inefficient in such shallow potential wells. Recoils can decrease the overall frequency of MBHs in small galaxies to ∼60 per cent, while they have little effect on the frequency of MBHs in large galaxies (at most a 20 per cent effect).

Gravitational wave bursts from collisions of primordial black holes in clusters

The rate of gravitational wave bursts from the mergers of massive primordial black holes in clusters is calculated. Such clusters of black holes can be formed through phase transitions in the early Universe. The central black holes in clusters can serve as the seeds of supermassive black holes in galactic nuclei. The expected burst detection rate by the LISA gravitational wave detector is estimated.