MINI REVIEW article

Front. Astron. Space Sci., 09 July 2020

Sec. Cosmology

Volume 7 - 2020 | https://doi.org/10.3389/fspas.2020.00038

Binary Black Hole Mergers: Formation and Populations

  • 1. Physics and Astronomy Department Galileo Galilei, University of Padova, Padova, Italy

  • 2. INFN–Padova, Padova, Italy

  • 3. INAF–Osservatorio Astronomico di Padova, Padova, Italy

Abstract

We review the main physical processes that lead to the formation of stellar binary black holes (BBHs) and to their merger. BBHs can form from the isolated evolution of massive binary stars. The physics of core-collapse supernovae and the process of common envelope are two of the main sources of uncertainty about this formation channel. Alternatively, two black holes can form a binary by dynamical encounters in a dense star cluster. The dynamical formation channel leaves several imprints on the mass, spin and orbital properties of BBHs.

1. Black Hole Formation From Single Stars: Where We Stand Now

About 4 years ago, the LIGO detectors obtained the first direct detection of gravitational waves, GW150914 (Abbott et al., ; Abbott et al., ,), associated with the merger of two black holes (BHs). This event marks the dawn of gravitational wave astronomy: we now know that binary black holes (BBHs) exist, can reach coalescence by gravitational wave emission, and are composed of BHs with mass ranging from few solar masses to ~ 50 M. Here, we review the main physical processes that lead to the formation of BBHs and to their merger. We restrict our attention to stellar-born BHs. As to primordial BHs, which might form from gravitational instabilities in the early Universe, we refer the reader to Carr et al. () and Belotsky et al. (), and references therein. Before we start discussing binaries, we must briefly summarize the state-of-the-art knowledge about stellar-origin BHs: this is a necessary step to understand their pairing mechanisms.

Stellar-mass BHs are thought to be the final outcome of the evolution of a massive star (zero-age main sequence mass mZAMS ≳ 20 M). Hence the mass of the BH should be affected by the two main processes that influence the evolution of a single star: (i) mass loss by stellar winds and (ii) the final collapse.

1.1. Stellar Winds

Hot (>104 K) massive stars (mZAMS ≳ 30 M) lose a non-negligible fraction of their mass by line-driven winds. This process depends on metallicity (Z): the mass-loss rate by stellar winds can be described as ṁ ∝ Zβ, where Z is the absolute metallicity (see e.g., Vink et al., and references therein). The most recent models suggest that β is not constant, but depends at least on the luminosity of the star (Gräfener and Hamann, ; Vink et al., ; Chen et al., ): the closer the luminosity L is to the Eddington value LEdd, the higher the mass loss, basically canceling the dependence on metallicity when LLEdd.

In single stars, stellar winds uniquely determine the final mass of the star at the onset of collapse. If we consider a star with mZAMS = 90 M and metallicity Z = 0.02 (i.e., approximately solar), its final mass will be only mfin ~ 30 M; while the same star with Z < 0.0002 has mfin ≳ 0.8 mZAMS. The final mass of a star mfin is the strongest upper limit to the mass of the BH.

1.2. Core-Collapse Supernovae

For example, Fryer () and Fryer and Kalogera () suggest that if the final mass of the star is mfin ≳ 40 M, the fate of the star is to collapse to a BH directly, without supernova, because the binding energy of the outer stellar layers is too big to be overcome by the explosion. Fryer et al. () elaborate on these early results proposing that the mass of the compact object depends not only on mfin but also on the final mass of the carbon-oxygen core. Alternatively, O'Connor and Ott () proposed the role of the compactness parameter : if the compactness is small (e.g., ξ2.5 ≤ 0.2 − 0.4), the SN explosion is successful, otherwise we expect the star to collapse directly. All of these simplified models as well as more sophisticated ones (e.g., Ertl et al., ) point toward a similar direction: if the star ends its life with a large final mass, its carbon-oxygen core grows larger, its compactness is generally higher, and so on. Hence, we expect that metal-poor stars, which retain a larger fraction of their mass to the very end and develop larger cores, are more likely to collapse to BHs directly, producing larger BHs (e.g., Mapelli et al., , ; Zampieri and Roberts, ; Belczynski et al., ). This simplified picture seems to agree with observations, but must be taken with several grains of salt: we need a vigorous step forward in core-collapse SN simulations and theoretical models, before we can draw robust conclusions (e.g., Burrows et al., ).

1.3. Pair Instability

Core-collapse SNe are not the only mechanism that can end the life of a massive star. When the helium core of a star grows to ≥ 60 M and the central temperature reaches ~ 109 K, electron and positron pairs are produced at an efficient rate, leading to a softening of the equation of state. The star undergoes pair instability (PI, Ober et al., ; Bond et al., ; Heger et al., ; Woosley et al., ): oxygen, neon, and silicon are burned explosively and the entire star is disrupted leaving no remnant, unless its helium core is ≥ 130 M. In the latter case, the gravity of the outer layers is so big that the star collapses to a massive BH directly as an effect of PI (Heger et al., ). Smaller helium cores (~ 30−60 M) are associated with a less dramatic manifestation of PI: the softened equation of state drives oscillations of the core (pulsational PI, Barkat et al., ; Woosley et al., ; Chen et al., ; Yoshida et al., ); during each oscillation the star sheds some mass till it finds a new equilibrium to a lower core mass, but leaves a BH smaller than expected without pulsational PI (Belczynski et al., ; Spera and Mapelli, ; Woosley, , ; Marchant et al., ; Stevenson et al., ; Renzo et al., ).

From the combination of PI, core-collapse SNe and stellar-wind mass loss prescriptions, we expect the mass spectrum of BHs to behave roughly as shown in Figure 1. In particular, PI is expected to carve a gap in the mass spectrum of BHs between ~ 50(−10, +20) M and ≈ 120 − 130 M. The uncertainty on this mass gap is mainly connected with uncertainties nuclear reaction rates (Farmer et al., ), on the collapse of the residual hydrogen envelope and on the role of stellar rotation (Mapelli et al., ). Within this framework, we predict a reasonable mass range for stellar-origin BHs to be ~ 3 − 65 M (assuming the most conservative value for the lower edge of pair-instability mass gap). If exotic metal-poor stars exist with mass mZAMS > 250 M, these might directly collapse to intermediate-mass BHs (IMBHs) with mass > 100 M.

Figure 1

), using PARSEC evolutionary tracks (Bressan et al., ) and the delayed model from Fryer et al. (). See Spera and Mapelli () for details.

2. Binary BH Formation in Isolation

The scenario highlighted in the previous section assumes that the progenitor star is single. But gravitational waves have shown the existence of BBHs with a very short orbital separation: the initial separation of a BBH must be of the order of few ten solar radii for the BBH to merge within a Hubble time by gravitational-wave emission. This challenges our understanding of binary star evolution. A close binary star undergoes several physical processes during its life, which can completely change its final fate (see e.g., Eggleton, ). The most important processes include mass transfer and common envelope, tides and natal kicks (Hurley et al., ).

Mass transfer and common envelope are crucial in this regard. After main sequence, a massive star can develop a stellar radius as large as several thousand solar radii. Hence, if this star is member of a binary system and its orbital separation is of few hundred to few thousand solar radii, the binary system undergoes Roche lobe overflow and possibly common envelope (Ivanova et al., ). If common envelope occurs between a BH and a giant companion star, the BH and the core of the giant star orbit about each other surrounded by the giant's envelope: they feel a strong gas drag from the envelope and lose kinetic energy, inspiralling about each other. This transfers thermal energy to the envelope, which might trigger the ejection of the envelope. If the envelope is not ejected, the binary system merges prematurely giving birth to a single BH. In contrast, if the envelope is ejected, the final binary system is composed of the BH and the core of the giant. Because of the spiral-in, the final semi-major axis of the binary is just few solar radii, much smaller than the initial one. If the naked core collapses to a BH without receiving a strong natal kick, the system becomes a BBH with a short orbital period, able to merge within a Hubble time. Unfortunately, our understanding of common envelope is still poor (see Fragos et al., for a recent simulation) and this uncertainty heavily affects our knowledge of BBH demography. The left-hand panel of Figure 2 is a schematic view of the isolated binary evolution channel through common envelope.

Figure 2

Several alternative scenarios to common envelope have been proposed (de Mink and Mandel, ; Mandel and de Mink, ; Marchant et al., ). For example, in the over-contact binary evolution, Marchant et al. () show that, when two massive stars in a tight binary are fast rotators, they remain fully mixed as a result of their tidally induced high spin; in this case, the binary avoids premature merger even if it is overfilling its Roche lobe and might evolve into a tight BBH.

The isolated binary evolution scenario has several characteristic signatures. In the common envelope isolated binary evolution scenario, the masses of the two BHs span from ~ 3 M up to ~ 45 M (see e.g., Giacobbo and Mapelli, ) and the mass ratios are preferentially close to 1 (although all mass ratios q = m2/m1 ≳ 0.1 are possible, see e.g., Giacobbo and Mapelli, ). Most processes in binary evolution tend to produce aligned spins (e.g., Rodriguez et al., ), while the magnitude of the spin is basically unconstrained (but see Qin et al., , ; Fuller and Ma, for some recent attempts to quantify spins). Mass transfer episodes and gravitational-wave decay are expected to efficiently damp eccentricity, so that almost all isolated binaries have near zero eccentricity in the LIGO-Virgo band. Finally, local merger rate densities span from a few to few thousand events Gpc−3 yr−1, depending on the details of common envelope and natal kicks (e.g., Dominik et al., ; Belczynski et al., ; Mapelli et al., ; Giacobbo and Mapelli, , ; Mapelli and Giacobbo, ; Neijssel et al., ; Santoliquido et al., ; Tang et al., ). The scenarios which include alternatives to common envelope predict an even stronger prevalence of systems with q ~ 1, a preferred mass range ~ 25 − 60 M (Marchant et al., ), high and aligned spins, zero eccentricity in the LIGO-Virgo band, and long delay times (≳ 3 Gyr, de Mink and Mandel, ). Local merger rate densities are expected to be ~ 10 Gpc−3 yr−1 (Mandel and de Mink, ), with large uncertainties.

3. Binary BH Formation in Star Clusters

Star clusters are among the densest places in the Universe. There is a plethora of star clusters, with their distinguishing features: (i) globular clusters (Gratton et al., ) are old (~ 12 Gyr) and massive systems (~ 104−6 M), (ii) nuclear star clusters can be even more massive (~ 107 M) and lie at the center of many galaxies, in some cases coexisting with the supersessive BH (Neumayer et al., ), (iii) open clusters and young star clusters (Portegies Zwart et al., ) are generally less massive (up to ~ 105 M) and short lived (less than a few Gyr), but are the main birthplace of massive stars in the local Universe (Lada and Lada, ).

The central density of star clusters is sufficiently high (≳ 103 stars pc−3) and their typical velocity dispersion sufficiently low (from a few to a few tens of km s−1, possibly with the exception of nuclear star clusters) that their central two-body relaxation time (Spitzer, ) is shorter than their lifetime. This has one fascinating implication: the orbits of stars and binary stars in a star cluster are constantly perturbed by dynamical encounters with other cluster members. This process affects the formation and the evolution of binary BHs in multiple ways (e.g., Portegies Zwart and McMillan, ).

Dynamical Exchanges

Dynamical exchanges occur when a binary system interacts with a single stellar object and the latter replaces one of the members of the binary. We have known for a long time that massive objects are more likely to acquire companions by dynamical exchanges (Hills and Fullerton, ). Since BHs are among the most massive objects in a star cluster, they are very efficient in forming new binaries through exchanges (e.g., Ziosi et al., ).

During a three-body encounter, a binary star exchanges a fraction of its internal energy with the third body. If the binary is particularly tight (hard binary), such encounters tend to harden the binary star, i.e., to increase its binding energy by reducing its semi-major axis (dynamical hardening). In the case of a BBH, this hardening might speed up the merger, because it drives the semi-major axis of the BBH in the regime where orbital decay by gravitational waves becomes efficient (see e.g., Figure 10 of Mapelli, ). On the other hand, the least massive BBHs can even be ionized, i.e., split by strong dynamical encounters with massive intruders.

Mergers of Massive Stars

Mergers of massive stars are common in dense young star clusters, because of the short dynamical friction timescale (Portegies Zwart et al., ). Under some assumptions, these mergers can lead to the formation of massive BHs (mBH > 60 M), with mass in the pair-instability gap (Di Carlo et al., ). In star clusters, such massive BHs can acquire a companion by dynamical exchanges, leading to the formation of BBHs in the mass gap. A fast sequence of stellar mergers in the dense core of a young star cluster (also known as runaway collision, Portegies Zwart et al., ; Giersz et al., ) might even lead to the formation of intermediate-mass BHs (IMBHs), i.e., BHs with mass mBH > 100 M, especially at low metallicity (Mapelli, ).

The dynamical processes we briefly summarized above (and in the right-hand panel of Figure 2) leave a clear imprint on BBHs. First, dynamically formed BBHs extend to higher masses than isolated BBHs: they might even be in the pair-instability mass gap or in the IMBH regime (Di Carlo et al., ; Rodriguez et al., ). Secondly, dynamical exchanges randomize the spin direction, leading to an isotropic distribution of BH spins. In contrast, isolated BBHs have a preference for aligned spins (Rodriguez et al., ; Gerosa et al., ). Third, dynamics can trigger the merger of BBHs with non-zero eccentricity even in the LIGO-Virgo band (Samsing et al., ; Rodriguez et al., ; Samsing, ; Zevin et al., ). These signatures provide an unique opportunity to differentiate among the isolated and the dynamical formation channel when the number of gravitational-wave detections will be of the order of a few hundreds (e.g., Zevin et al., ; Bouffanais et al., ).

Statements

Author contributions

MM collected the material and wrote the review.

Funding

MM acknowledges financial support from the European Research Council for the ERC Consolidator grant DEMOBLACK, under contract no. 770017.

Conflict of interest

The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Summary

Keywords

stars: black holes, black hole physics, Galaxy: open clusters and associations: general, stars: kinematics and dynamics, gravitational waves

Citation

Mapelli M (2020) Binary Black Hole Mergers: Formation and Populations. Front. Astron. Space Sci. 7:38. doi: 10.3389/fspas.2020.00038

Received

25 February 2020

Accepted

03 June 2020

Published

09 July 2020

Volume

7 - 2020

Edited by

Rosalba Perna, Stony Brook University, United States

Reviewed by

Maxim Yurievich Khlopov, UMR7164 Astroparticule et Cosmologie (APC), France; Ataru Tanikawa, The University of Tokyo, Japan

Updates

Copyright

*Correspondence: Michela Mapelli

This article was submitted to Cosmology, a section of the journal Frontiers in Astronomy and Space Sciences

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics