Mass transfer in binaries with massive donors and compact companions, when the donors rapidly evolve after their main sequence, determines the formation rates of merging double stellar-mass black ...hole (BH) binaries formed outside clusters. This mass transfer was previously postulated to be unstable and was expected to lead to a common envelope event. The common envelope event then ends with either the merger of the two stars or formation of a binary that eventually may become a merging double BH. We revisit the stability of this mass transfer and find an unanticipated third outcome: for a large range of binary orbital separations, this mass transfer is stable. This newly found stability allows us to reconcile the empirical rate obtained by LIGO, 9-240 Gpc... yr..., with the theoretical rate for double BH binary mergers predicted by population synthesis studies by excluding a channel that predicts a merger rate above 1000 Gpc... yr... Furthermore, the stability of the mass transfer leads to the formation of ultraluminous X-ray sources. The theoretically predicted formation rates of bright ultraluminous X-ray sources powered by a stellar-mass BH are high enough to explain the number of observed bright ultraluminous X-ray sources. (ProQuest: ... denotes formulae/symbols omitted.)
ABSTRACT The initial mass function (IMF), binary fraction, and distributions of binary parameters (mass ratios, separations, and eccentricities) are indispensable inputs for simulations of stellar ...populations. It is often claimed that these are poorly constrained, significantly affecting evolutionary predictions. Recently, dedicated observing campaigns have provided new constraints on the initial conditions for massive stars. Findings include a larger close binary fraction and a stronger preference for very tight systems. We investigate the impact on the predicted merger rates of neutron stars and black holes. Despite the changes with previous assumptions, we only find an increase of less than a factor of 2 (insignificant compared with evolutionary uncertainties of typically a factor of 10-100). We further show that the uncertainties in the new initial binary properties do not significantly affect (within a factor of 2) our predictions of double compact object merger rates. An exception is the uncertainty in IMF (variations by a factor of 6 up and down). No significant changes in the distributions of final component masses, mass ratios, chirp masses, and delay times are found. We conclude that the predictions are, for practical purposes, robust against uncertainties in the initial conditions concerning binary parameters, with the exception of the IMF. This eliminates an important layer of the many uncertain assumptions affecting the predictions of merger detection rates with the gravitational wave detectors aLIGO/aVirgo.
Aims.
We present an open-access database that includes a synthetic catalog of black holes (BHs) in the Milky Way, divided by the components disk, bulge, and halo.
Methods.
To calculate the evolution ...of single and binary stars, we used the updated population synthesis code StarTrack. We applied a new model of the star formation history and chemical evolution of Galactic disk, bulge, and halo that was synthesized from observational and theoretical data. This model can be easily employed for other studies of population evolution.
Results.
We find that at the current Milky Way (disk+bulge+halo) contains about 1.2 × 10
8
single BHs with an average mass of about 14
M
⊙
, and 9.3 × 10
6
BHs in binary systems with an average mass of 19
M
⊙
. We present basic statistical properties of the BH population in three Galactic components such as the distributions of BH masses, velocities, or the numbers of BH binary systems in different evolutionary configurations.
Conclusions.
The metallicity of a stellar population has a significant effect on the final BH mass through the stellar winds. The most massive single BH in our simulation of 113
M
⊙
originates from a merger of a BH and a helium star in a low-metallicity stellar environment in the Galactic halo. We constrain that only ∼0.006% of the total Galactic halo mass (including dark matter) can be hidden in the form of stellar origin BHs. These BHs cannot be detected by current observational surveys. We calculated the merger rates for current Galactic double compact objects (DCOs) for two considered common-envelope models: ∼3–81 Myr
−1
for BH-BH, ∼1–9 Myr
−1
for BH-neutron star (NS), and ∼14–59 Myr
−1
for NS-NS systems. We show the evolution of the merger rates of DCOs since the formation of the Milky Way until the current moment with the new star formation model of the Galaxy.
Aims.
It is speculated that a merger of two massive stellar-origin black holes in a dense stellar environment may lead to the formation of a massive black hole in the pair-instability mass gap ...(∼50−135
M
⊙
). Such a merger-formed black hole is expected to typically have a high spin (
a
∼ 0.7). If such a massive black hole acquires another black hole it may lead to another merger detectable by LIGO/Virgo in gravitational waves. Acquiring a companion may be hindered by gravitational-wave kick/recoil, which accompanies the first merger and may quickly remove the massive black hole from its parent globular or nuclear cluster. We test whether it is possible for a massive merger-formed black hole in the pair-instability gap to be retained in its parent cluster and have low spin. Such a black hole would be indistinguishable from a primordial black hole.
Methods.
We employed results from numerical relativity calculations of black hole mergers to explore the range of gravitational-wave recoil velocities for various combinations of merging black hole masses and spins. We compared merger-formed massive black hole speeds with typical escape velocities from globular and nuclear clusters.
Results.
We show that a globular cluster is highly unlikely to form and retain a ∼100
M
⊙
black hole if the spin of the black hole is low (
a
≲ 0.3). Massive merger-formed black holes with low spins acquire high recoil speeds (≳ 200 km s
−1
) from gravitational-wave kick during formation that exceed typical escape speeds from globular clusters (∼ 50 km s
−1
). However, a very low-spinning (
a
∼ 0.1) and massive (∼100
M
⊙
) black hole could be formed and retained in a galactic nuclear star cluster. Even though such massive merger-formed black holes with such low spins acquire high speeds during formation (∼ 400 km s
−1
), they may avoid ejection since massive nuclear clusters have high escape velocities (∼ 300−500 km s
−1
). A future detection of a massive black hole in the pair-instability mass gap with low spin would therefore not be proof of the existence of primordial black holes, which are sometimes claimed to have low spins and arbitrarily high masses.
The treatment and criteria for development of unstable Roche lobe overflow (RLOF) that leads to the common envelope (CE) phase have hindered the area of evolutionary predictions for decades. In ...particular, the formation of black hole-black hole (BH-BH), black hole-neutron star (BH-NS), and neutron star-neutron star (NS-NS) merging binaries depends sensitively on the CE phase in classical isolated binary evolution model. All these mergers are now reported as LIGO/Virgo sources or source candidates. CE is even considered by some as a mandatory phase in the formation of BH-BH, BH-NS, or NS-NS mergers in binary evolution models. At the moment, there is no full first-principles model for the development of the CE. We employed the
StarTrack
population synthesis code to test the current advancements in studies on the stability of RLOF for massive donors to assess their effect on the LIGO/Virgo source population. In particular, we allowed for more restrictive CE development criteria for massive donors (
M
> 18
M
⊙
). We also tested a modified condition for switching between different types of stable mass transfer and between the thermal or nuclear timescale. The implemented modifications significantly influence the basic properties of merging double compact objects, sometimes in non-intuitive ways. For one of the tested models, with restricted CE development criteria, the local merger rate density for BH-BH systems increased by a factor of 2–3 due to the emergence of a new dominant formation scenario without any CE phase. We find that the changes in highly uncertain assumptions on RLOF physics may significantly affect: (i) the local merger rate density; (ii) shape of the mass and mass ratio distributions; and (iii) dominant evolutionary formation (with and without CE) scenarios of LIGO/Virgo sources. Our results demonstrate that without sufficiently strong constraints on RLOF physics, it is not possible to draw fully reliable conclusions about the population of double compact object systems based on population synthesis studies.
The distributions of the initial main-sequence binary parameters are one of the key ingredients in obtaining evolutionary predictions for compact binary (BH–BH/BH–NS/NS–NS) merger rates. Until now, ...such calculations were done under the assumption that initial binary parameter distributions were independent. For the first time, we implement empirically derived inter-correlated distributions of initial binary parameters primary mass (M1), mass ratio (q), orbital period (P), and eccentricity (e). Unexpectedly, the introduction of inter-correlated initial binary parameters leads to only a small decrease in the predicted merger rates by a factor of ≲2–3 relative to the previously used non-correlated initial distributions. The formation of compact object mergers in the isolated classical binary evolution favours initial binaries with stars of comparable masses (q ≈ 0.5–1) at intermediate orbital periods (log P (days) = 2–4). New distributions slightly shift the mass ratios towards lower values with respect to the previously used flat q distribution, which is the dominant effect decreasing the rates. New orbital periods (∼1.3 more initial systems within log P (days) = 2–4), together with new eccentricities (higher), only negligibly increase the number of progenitors of compact binary mergers. Additionally, we discuss the uncertainty of merger rate predictions associated with possible variations of the massive-star initial mass function (IMF). We argue that evolutionary calculations should be normalized to a star formation rate (SFR) that is obtained from the observed amount of UV light at wavelength 1500 Å (an SFR indicator). In this case, contrary to recent reports, the uncertainty of the IMF does not affect the rates by more than a factor of ∼2. Any change to the IMF slope for massive stars requires a change of SFR in a way that counteracts the impact of IMF variations on compact object merger rates. In contrast, we suggest that the uncertainty in cosmic SFR at low metallicity can be a significant factor at play.
Anisotropy of gravitational wave (GW) emission results in a net momentum gained by the black hole (BH) merger product, leading to a recoil velocity up to ∼103 km s−1, which may kick it out of a ...globular cluster (GC). We estimate GW kick retention fractions of merger products assuming different models for BH spin magnitude and orientation. We check how they depend on BH–BH merger time and properties of the cluster. We analyse the implications of GW kick retention fractions on intermediate massive BH formation by repeated mergers in a GC. We also calculate final spin of the merger product, and investigate how it correlates with other parameters: effective spin of the binary and gravitational kick velocity. We used data from MOCCA (MOnte Carlo Cluster simulAtor) GC simulations to get a realistic sample of BH–BH mergers, assigned each BH spin value according to a studied model, and calculated recoil velocity and final spin based on most recent theoretical formulas. We discovered that for physically motivated models, GW kick retention fractions are about 30 per cent 30 per cent and display small dependence on assumptions about spin, but are much more prone to cluster properties. In particular, we discovered a strong dependence of GW kick retention fractions on cluster density. We also show that GW kick retention fractions are high in final life stages of the cluster, but low at the beginning. Finally, we derive formulas connecting final spin with effective spin for primordial binaries, and with maximal effective spin for dynamical binaries.
The first neutron star-neutron star (NS-NS) merger was discovered on August 17, 2017 through gravitational waves (GW170817) and followed with electromagnetic observations. This merger was detected in ...an old elliptical galaxy with no recent star formation. We perform a suite of numerical calculations to understand the formation mechanism of this merger. We probe three leading formation mechanisms of double compact objects: classical isolated binary star evolution, dynamical evolution in globular clusters, and nuclear cluster formation to test whether they are likely to produce NS-NS mergers in old host galaxies. Our simulations with optimistic assumptions show current NS-NS merger rates at the level of 10−2 yr−1 from binary stars, 5 × 10−5 yr−1 from globular clusters, and 10−5 yr−1 from nuclear clusters for all local elliptical galaxies (within 100 Mpc3). These models are thus in tension with the detection of GW170817 with an observed rate of 1.5−1.2+3.2$1.5^{+3.2}_{-1.2}$1.5−1.2+3.2 yr−1 (per 100 Mpc3; LIGO/Virgo 90% credible limits). Our results imply that either the detection of GW170817 by LIGO/Virgo at their current sensitivity in an elliptical galaxy is a statistical coincidence; that physics in at least one of our three models is incomplete in the context of the evolution of stars that can form NS-NS mergers; or that another very efficient (unknown) formation channel with a long delay time between star formation and merger is at play.
Abstract
The LIGO/Virgo gravitational-wave observatories have detected at least 50 double black hole (BH) coalescences. This sample is large enough to have allowed several recent studies to draw ...conclusions about the implied branching ratios between isolated binaries versus dense stellar clusters as the origin of double BHs. It has also led to the exciting suggestion that the population is highly likely to contain primordial BHs. Here we demonstrate that such conclusions cannot yet be robust because of the large current uncertainties in several key aspects of binary stellar evolution. These include the development and survival of a common envelope, the mass and angular-momentum loss during binary interactions, mixing in stellar interiors, pair-instability mass loss, and supernova outbursts. Using standard tools such as the rapid population synthesis codes
StarTrack
and
COMPAS
and the detailed stellar evolution code
MESA
, we examine as a case study the possible future evolution of Melnick 34, the most massive known binary star system (with initial component masses of 144
M
⊙
and 131
M
⊙
). We show that, despite its fairly well-known orbital architecture, various assumptions regarding stellar and binary physics predict a wide variety of outcomes: from a close BH–BH binary (which would lead to a potentially detectable coalescence), through a wide BH–BH binary (which might be seen in microlensing observations), or a Thorne–Żytkow object, to a complete disruption of both objects by a pair-instability supernova. Thus, because the future of massive binaries is inherently uncertain, sound predictions about the properties of BH–BH systems formed in the isolated binary evolution scenario are highly challenging at this time. Consequently, it is premature to draw conclusions about the formation channel branching ratios that involve isolated binary evolution for the LIGO/Virgo BH–BH merger population.