Progenitors Of Black Holes Research Articles

Investigating 39 Galactic Wolf-Rayet stars with VLTI/GRAVITY

Context. Wolf-Rayet stars (WRs) represent one of the final evolutionary stages of massive stars and are thought to be the immediate progenitors of stellar-mass black holes. Their multiplicity characteristics form an important anchor point in single and binary population models for predicting gravitational-wave progenitors. Recent spectroscopic campaigns have suggested incompatible multiplicity fractions and period distributions for N- and C-rich Galactic WRs (WNs and WCs) at both short and long orbital periods, in contradiction with evolutionary model predictions. Aims. In this work, we employed long-baseline infrared interferometry to investigate the multiplicity of WRs at long periods and explored the nature of their companions. We present a magnitude-limited (K < 9; V < 14) survey of 39 Galactic WRs, including 11 WN, 15 WC, and 13 H-rich WN (WNh) stars. Methods. We used the K-band instrument GRAVITY at the Very Large Telescope Interferometer (VLTI) in Chile. The sensitivity of GRAVITY at spatial scales of ∼1 to 200 milliarcseconds and flux contrast of 1% allowed an exploration of periods in the range 102 − 105 d and companions down to ∼5 M⊙. We carried out a companion search for all our targets, with the aim of either finding wide companions or calculating detection limits. We also explored the rich GRAVITY dataset beyond a multiplicity search to look for other interesting properties of the WR sample. Results. We detected wide companions with VLTI/GRAVITY for only four stars in our sample: WR 48, WR 89, WR 93, and WR 115. Combining our results with spectroscopic studies, we arrived at observed multiplicity fractions of fobsWN = 0.55 ± 0.15, fobsWC = 0.40 ± 0.13, and fobsWNh = 0.23 ± 0.12. The multiplicity fractions and period distributions of WNs and WCs are consistent in our sample. For single WRs, we placed upper limits on the mass of potential companions down to ∼5 M⊙ for WNs and WCs, and ∼7 M⊙ for WNh stars. In addition, we also found other features in the GRAVITY dataset, such as (i) a diffuse extended component contributing significantly to the K-band flux in over half the WR sample; (ii) five known spectroscopic binaries resolved in differential phase data, which constitutes an alternative detection method for close binaries; and (iii) spatially resolved winds in four stars: WR 16, WR 31a, WR 78, and WR 110. Conclusions. Our survey reveals a lack of intermediate- (a few hundred days) and long- (a few years to decades) period WR systems. The 200d peak in the period distributions of WR+OB and BH+OB binaries predicted by Case B mass-transfer binary evolution models is not seen in our data. The rich companionship of their O-type progenitors in this separation range suggests that the WR progenitor stars expand and interact with their companions, most likely through unstable mass transfer, resulting in either a short-period system or a merger.

The black hole low-mass X-ray binary V404 Cygni is part of a wide triple.

Evidence suggests that, when compact objects such as black holes and neutron stars form, they may receive a 'natal kick', during which the stellar remnant gains momentum. Observational evidence for neutron star kicks is substantial1,2, yet is limited for black hole natal kicks, and some proposed black hole formation scenarios result in very small kicks3-5. Here we report that the canonical black hole low-mass X-ray binary (LMXB) V404 Cygni is part of a wide hierarchical triple with a tertiary companion at least 3,500 astronomical units (AU) away from the inner binary. Given the orbital configuration, the black hole probably received a sub-5 km s-1 kick to have avoided unbinding the tertiary. This discovery lends support to the idea that at least some black holes form with nearly no natal kick. Furthermore, the tertiary in this system lends credence to evolutionary models of LMXBs involving a hierarchical triple structure6. Remarkably, the tertiary is evolved, indicating that the system formed 3-5 billion years ago and that the black hole has removed at least half a solar mass of matter from its evolved secondary companion. During the event in which the black hole formed, it is required that at least half of the mass of the black hole progenitor collapsed into the black hole; it may even have undergone a complete implosion, enabling the tertiary to remain loosely bound.

Sequential formation of supermassive stars and heavy seed BHs through the interplay of cosmological cold accretion and stellar radiative feedback

ABSTRACT Supermassive stars (SMSs) and heavy seed black holes, as their remnants, are promising candidates for supermassive black hole (SMBH) progenitors, especially for ones observed in the early universe $z\simeq 8.5-10$ by recent JWST observations. Expected cradles of SMSs are the atomic cooling haloes ($M_{\rm halo}\simeq 10^7\ \mathrm{M}_{\odot }$), where ‘cold accretion’ emerges and possibly forms SMSs. We perform a suit of cosmological radiation hydrodynamics simulations and investigate star formation after the emergence of cold accretion, solving radiative feedback from stars inside the halo. We follow the mass growth of the protostars for $\sim 3\ \mathrm{Myr}$, resolving the gas inflow down to $\sim 0.1\ \mathrm{pc}$ scales. We discover that, after cold accretion emerges, multiple SMSs of $m_{\star }\gtrsim 10^5\ \mathrm{M}_{\odot }$ form at the halo centre with the accretion rates maintained at $\dot{m}_{\star }\simeq 0.04\ \mathrm{M}_\odot \mathrm{yr}^{-1}$ for $\lesssim 3\ \mathrm{Myr}$. Cold accretion supplies gas at a rate of $\dot{M}_{\rm gas}\gtrsim 0.01-0.1\ \mathrm{M}_\odot \mathrm{yr}^{-1}$ from outside the halo virial radius to the central gas disc. Gravitational torques from spiral arms transport gas further inwards, which feeds the SMSs. Radiative feedback from stars suppresses $\mathrm{H}_2$ cooling and disc fragmentation, while photoevaporation is prevented by a dense envelope, which attenuates ionizing radiation. Our results suggest that cold accretion can bring efficient BH mass growth after seed formation in the later universe. Moreover, cold accretion and gas migration inside the central disc increase the mass concentration and provide a promising formation site for the extremely compact stellar clusters observed by JWST.

Binary black hole mergers from Population III star clusters

Binary black holes (BBHs) that are born from the evolution of Population III (Pop. III) stars are one of the main high-redshift targets for next-generation ground-based gravitational-wave (GW) detectors. Their predicted initial mass function and lack of metals make them the ideal progenitors of black holes above the upper edge of the pair-instability mass gap, that is, with a mass higher than ~134 (241) M⊙ for stars that become (or do not become) chemically homogeneous during their evolution. We investigate the effects of cluster dynamics on the mass function of BBHs that are born from Pop. III stars by considering the main uncertainties on the mass function of Pop. III stars, the orbital properties of the binary systems, the star cluster mass, and the disruption time. In our dynamical models, at least ~5% and up to 100% BBH mergers in Pop. III star clusters have a primary mass m1 above the upper edge of the pair-instability mass gap. In contrast, only ≲3% isolated BBH mergers have a primary mass above the gap, unless their progenitors evolved as chemically homogeneous stars. The lack of systems with a primary and/or secondary mass inside the gap defines a zone of avoidance with sharp boundaries in the plane of the primary mass-mass ratio. Finally, we estimate the merger rate density of BBHs. In the most optimistic case, we find a maximum of ℛ ~ 200 Gpc-3 yr-1 at ɀ ~ 15 for BBHs that formed via dynamical capture. For comparison, the merger rate density of isolated Pop. III BBHs is ℛ ≤ 10 Gpc−3 yr−1 for the same model of Pop. III star formation history.

Gaia BH1 and BH2 — Evolutionary Models with Overshooting of the Black Hole Progenitors within the Present-Day Binary Separation

Abstract Three black holes (BHs) in wide binaries — Gaia BH1, BH2 and BH3 — were recently discovered. The likely progenitors of the BHs were massive stars that experienced a supergiant phase, reaching radii of ∼1000R⊙, before collapsing to form the BH. Such radii are difficult to accommodate with the present-day orbits of BH1 and BH2 — with semi-major axes of 1.4 and 3.7 au, respectively. In this letter, we show that the maximal radii of the supergiants are not necessarily so large, and realistic stellar evolution models, with some assumed overshooting above the convective core into the radiative stellar envelope, produce substantially smaller maximal radii. The limited expansion of supergiants is consistent with the empirical Humphreys-Davidson limit — the absence of red supergiants above an upper luminosity limit, notably lower than the highest luminosity of main-sequence stars. We propose that the evolution that led to the formation of Gaia BH1 and BH2 simply did not involve an expansion to the cool supergiant phase.

Intermediate-mass Black Hole Progenitors from Stellar Collisions in Dense Star Clusters

Very massive stars (VMSs) formed via a sequence of stellar collisions in dense star clusters have been proposed as the progenitors of massive black hole seeds. VMSs could indeed collapse to form intermediate-mass black holes, which would then grow by accretion to become the supermassive black holes observed at the centers of galaxies and powering high-redshift quasars. Previous studies have investigated how different cluster initial conditions affect the formation of a VMS, including mass segregation, stellar collisions, and binaries, among others. In this study, we investigate the growth of VMSs with a new grid of Cluster Monte Carlo star cluster simulations—the most expansive to date. The simulations span a wide range of initial conditions, varying the number of stars, cluster density, stellar initial mass function (IMF), and primordial binary fraction. We find a gradual shift in the mass of the most massive collision product across the parameter space; in particular, denser clusters born with top-heavy IMFs provide strong collisional regimes that form VMSs with masses easily exceeding 1000 M ⊙. Our results are used to derive a fitting formula that can predict the typical mass of a VMS formed as a function of the star cluster properties. Additionally, we study the stochasticity of this process and derive a statistical distribution for the mass of the VMS formed in one of our models, recomputing the model 50 times with different initial random seeds.

On the formation of a 33 solar-mass black hole in a low-metallicity binary

A 33M⊙ black hole (BH) was recently discovered in an 11.6-year binary only 590 pc from the Sun. The system, Gaia BH3, contains a 0.8M⊙ low-metallicity giant () that is a member of the ED-2 stellar stream. This paper investigates whether the system could have formed via isolated binary evolution. I construct evolutionary models for metal-poor massive stars with initial masses ranging from 35−55M⊙, which reach maximum radii of 1150−1800R⊙ as red supergiants. I then explore what combinations of initial orbit, mass loss, and natal kick can produce the period and eccentricity of Gaia BH3. Initial orbits wide enough to accommodate the BH progenitor as a red supergiant can match the observed period and eccentricity, but only if the BH formed with a significant natal kick (). These models are disfavored because such a kick would have ejected the binary from the ED-2 progenitor cluster. I conclude that Gaia BH3 likely formed through dynamical interactions, unless the BH progenitor did not expand to red supergiant dimensions. Only about 1 in 10,000 stars in the solar neighborhood have metallicities as low as Gaia BH3. This suggests that BH companions are dramatically over-represented at low-metallicity, though caveats related to small number statistics apply. The fact that the luminous star in Gaia BH3 has been a giant – greatly boosting its detectability – only for $$1% of the time since the system’s formation implies that additional massive BHs remain to be discovered with only moderately fainter companions. Both isolated and dynamically-formed BH binaries with orbits similar to Gaia BH3 are likely to be discovered in Gaia DR4.

Radio Emission from High-redshift Active Galactic Nuclei in the JADES and CEERS Surveys

Recent calculations indicate that radio emission from quasars at z ∼ 6–7 could be detected at much earlier stages of evolution, at z ∼ 14–15, by the Next-Generation Very Large Array (ngVLA) and the Square Kilometer Array (SKA). However, the James Webb Space Telescope has now discovered less luminous active galactic nuclei (AGNs) at z > 4 and a few massive black holes (BHs) at z > 10, which may be the progenitors of supermassive black holes (SMBHs) but at different stages of growth. Radio detections of these new AGNs would provide complementary measures of their properties and those of their host galaxies. Here we estimate radio flux densities for 19 new AGNs found by the JADES, CEERS, and UNCOVER surveys. We find that ngVLA should be able to detect most of these sources in targeted surveys with integration times of 10–100 hr (and in just 1 hr for a few of them) but most would require at least 100 hr of SKA time in spite of its greater sensitivities at low frequencies. In some cases, radio emission from the BH can be distinguished from that of H ii regions and supernovae in their host galaxies, which could be used to estimate their star formation rates. Such detections would be yet another example of the useful synergies between near-infrared and radio telescopes in SMBH science in the coming decade.

Investigating the Chemically Homogeneous Evolution Channel and Its Role in the Formation of the Enigmatic Binary Black Hole Progenitor Candidate HD 5980

Chemically homogeneous evolution (CHE) is a promising channel for forming massive binary black holes. The enigmatic, massive Wolf–Rayet binary HD 5980 A&B has been proposed to have formed through this channel. We investigate this claim by comparing its observed parameters with CHE models. Using MESA, we simulate grids of close massive binaries, then use a Bayesian approach to compare them with the stars’ observed orbital period, masses, luminosities, and hydrogen surface abundances. The most probable models, given the observational data, have initial periods ∼3 days, widening to the present-day ∼20 days orbit as a result of mass loss—correspondingly, they have very high initial stellar masses (≳150 M ⊙). We explore variations in stellar-wind mass loss and internal mixing efficiency, and find that models assuming enhanced mass loss are greatly favored to explain HD 5980, while enhanced mixing is only slightly favored over our fiducial assumptions. Our most probable models slightly underpredict the hydrogen surface abundances. Regardless of its prior history, this system is a likely binary black hole progenitor. We model its further evolution under our fiducial and enhanced wind assumptions, finding that both stars produce black holes with masses ∼19–37 M ⊙. The projected final orbit is too wide to merge within a Hubble time through gravitational waves alone. However, the system is thought to be part of a 2+2 hierarchical multiple. We speculate that secular effects with the (possible) third and fourth companions may drive the system to promptly become a gravitational-wave source.

Which Black Hole Is Spinning? Probing the Origin of Black Hole Spin with Gravitational Waves

Theoretical studies of angular momentum transport suggest that isolated stellar-mass black holes are born with negligible dimensionless spin magnitudes χ ≲ 0.01. However, recent gravitational-wave observations indicate ≳40% of binary black hole systems contain at least one black hole with a nonnegligible spin magnitude. One explanation is that the firstborn black hole spins up the stellar core of what will become the second-born black hole through tidal interactions. Typically, the second-born black hole is the “secondary” (less massive) black hole though it may become the “primary” (more massive) black hole through a process known as mass-ratio reversal. We investigate this hypothesis by analyzing data from the third gravitational-wave transient catalog using a “single-spin” framework in which only one black hole may spin in any given binary. Given this assumption, we show that at least 28% (90% credibility) of the LIGO–Virgo–KAGRA binaries contain a primary with significant spin, possibly indicative of mass-ratio reversal. We find no evidence for binaries that contain a secondary with significant spin. However, the single-spin framework is moderately disfavored (natural log Bayes factor lnB=3.1 ) when compared to a model that allows both black holes to spin. If future studies can firmly establish that most merging binaries contain two spinning black holes, it may call into question our understanding of formation mechanisms for binary black holes or the efficiency of angular momentum transport in black hole progenitors.

A Triple Scenario for the Formation of Wide Black Hole Binaries Such as Gaia BH1

Recently, several noninteracting black hole–stellar binaries have been identified in Gaia data—for example, Gaia BH1, where a Sun-like star is in a moderately eccentric (e = 0.44) 185 days orbit around a black hole. This orbit is difficult to explain through binary evolution. The present-day separation suggests the progenitor system would have undergone an episode of common-envelope evolution, but a common envelope should shrink the period below the observed one. Since the majority of massive stars form in higher-multiplicity systems, a triple evolution scenario is more likely for the progenitors of BH binaries. Here we show that such systems can indeed be more easily explained via evolution in hierarchical triple systems. von Zeipel–Lidov–Kozai oscillations or instabilities can delay the onset of the common-envelope phase in the inner binary of the triple, so that the black hole progenitor and low-mass star are more widely separated when it begins, leading to the formation of wider binaries. There are also systems with similar periods but larger eccentricities, where the BH progenitor is a merger product of the inner binary in the triple. Such mergers lead to a more top-heavy black hole mass function.

Optical Properties and Variability of the Be X-Ray Binary CPD-29 2176

Be X-ray binaries (Be XRBs) are high-mass X-ray binaries, with a neutron star or black hole orbiting and accreting material from a nonsupergiant B-star that is rotating at a near critical rate. These objects are prime targets to understand past binary interactions as the neutron star or black hole progenitor likely experienced Roche lobe overflow to spin up the Be star we observe now. The stellar variability can then allow us to explore the stellar structure of these objects. It was recently demonstrated that the high-mass X-ray binary CPD −29 2176 descended from an ultrastripped supernova and is a prime target to evolve into an eventual binary neutron star and kilonova. We present the photometric variability from both TESS and ASAS along with the spectral properties and disk variability of the system in this paper. All of the optical lines are contaminated with disk emission except for the He ii λ4686 absorption line. The disk variability timescales are not the same as the orbital timescale, but could be related to the X-ray outbursts that have been recorded by Swift. We end our study with a discussion comparing CPD −29 2176 to classical Be stars and other Be X-ray binaries, finding the stellar rotation to be near a frequency of 1.5 cycles day−1, and exhibiting incoherent variability in three frequency groups.

Dynamical formation of Gaia BH1 in a young star cluster

ABSTRACT Gaia BH1, the first quiescent black hole (BH) detected from Gaia data, poses a challenge to most binary evolution models: its current mass ratio is ≈0.1, and its orbital period seems to be too long for a post-common envelope system and too short for a non-interacting binary system. Here, we explore the hypothesis that Gaia BH1 formed through dynamical interactions in a young star cluster (YSC). We study the properties of BH-main sequence (MS) binaries formed in YSCs with initial mass 3 × 102–3 × 104 M⊙ at solar metallicity, by means of 3.5 × 104 direct N-body simulations coupled with binary population synthesis. For comparison, we also run a sample of isolated binary stars with the same binary population synthesis code and initial conditions used in the dynamical models. We find that BH-MS systems that form via dynamical exchanges populate the region corresponding to the main orbital properties of Gaia BH1 (period, eccentricity, and masses). In contrast, none of our isolated binary systems match the orbital period and MS mass of Gaia BH1. Our best-matching Gaia BH1-like system forms via repeated dynamical exchanges and collisions involving the BH progenitor star, before it undergoes core collapse. YSCs are at least two orders of magnitude more efficient in forming Gaia BH1-like systems than isolated binary evolution.

Massive binary black holes from Population II and III stars

ABSTRACT Population III stars, born from the primordial gas in the Universe, lose a negligible fraction of their mass via stellar winds and possibly follow a top-heavy mass function. Hence, they have often been regarded as the ideal progenitors of massive black holes (BHs), even above the pair instability mass gap. Here, we evolve a large set of Population III binary stars (metallicity Z = 10−11) with our population-synthesis code sevn, and compare them with Population II binary stars (Z = 10−4). In our models, the lower edge of the pair-instability mass gap corresponds to a BH mass of ≈86 (≈91) M⊙ for single Population III (II) stars. Overall, we find only mild differences between the properties of binary BHs (BBHs) born from Population III and II stars, especially if we adopt the same initial mass function and initial orbital properties. Most BBH mergers born from Population III and II stars have primary BH mass below the pair-instability gap, and the maximum secondary BH mass is <50 M⊙. Only up to ≈3.3 per cent (≈0.09 per cent) BBH mergers from Population III (II) progenitors have primary mass above the gap. Unlike metal-rich binary stars, the main formation channel of BBH mergers from Population III and II stars involves only stable mass transfer episodes in our fiducial model.

The role of stellar expansion on the formation of gravitational wave sources

ABSTRACT Massive stars are the progenitors of black holes and neutron stars, the mergers of which can be detected with gravitational waves (GW). The expansion of massive stars is one of the key factors affecting their evolution in close binary systems, but it remains subject to large uncertainties in stellar astrophysics. For population studies and predictions of GW sources, the stellar expansion is often simulated with the analytic formulae from Hurley et al. (2000). These formulae need to be extrapolated and are often considered outdated. In this work, we present five different prescriptions developed from 1D stellar models to constrain the maximum expansion of massive stars. We adopt these prescriptions to investigate how stellar expansion affects mass transfer interactions and in turn the formation of GW sources. We show that limiting radial expansion with updated 1D stellar models, when compared to the use of Hurley et al. (2000) radial expansion formulae, does not significantly affect GW source properties (rates and masses). This is because most mass transfer events leading to GW sources are initialized in our models before the donor star reaches its maximum expansion. The only significant difference was found for the mass distribution of massive binary black hole mergers (Mtot > 50 M⊙) formed from stars that may evolve beyond the Humphreys–Davidson limit, whose radial expansion is the most uncertain. We conclude that understanding the expansion of massive stars and the origin of the Humphrey–Davidson limit is a key factor for the study of GW sources.

Understanding the progenitor formation galaxies of merging binary black holes

ABSTRACT With nearly a hundred gravitational wave detections, the origin of black hole mergers has become a key question. Here, we focus on understanding the typical galactic environment in which binary black hole (BBH) mergers arise. To this end, we synthesize progenitors of BBH mergers as a function of the redshift of progenitor formation, present-day formation galaxy mass, and progenitor stellar metallicity for 240 star formation and binary evolution models. We provide guidelines to infer the formation galaxy properties and time of formation, highlighting the interplay between the star formation rate and the efficiency of forming merging BBHs from binary stars, both of which strongly depend on metallicity. We find that across models, over 50 per cent of BBH mergers have a progenitor metallicity of a few tenths of Solar metallicity, however, inferring formation galaxy properties strongly depends on both the binary evolution model and global metallicity evolution. The numerous, low-mass black holes (≲ 15 M⊙) trace the bulk of the star formation in galaxies heavier than the Milky Way (MGal ≳ 1010.5 M⊙). In contrast, heavier BBH mergers typically stem from larger black holes forming in lower metallicity dwarf galaxies (MGal ≲ 109 M⊙). We find that the progenitors of detectable BBHs tend to arise from dwarf galaxies at a lower formation redshift (≲1). We also produce a posterior probability of the progenitor environment for any detected gravitational wave signal. For the massive GW150914 merger, we show that it likely came from a very low-metallicity (Z ≲ 0.025 Z⊙) environment.

Constraining the merger history of primordial-black-hole binaries from GWTC-3

Primordial black holes (PBHs) can be not only cold dark matter candidates but also progenitors of binary black holes observed by LIGO-Virgo-KAGRA (LVK) Collaboration. The PBH mass can be shifted to the heavy distribution if multi-merger processes occur. In this work, we constrain the merger history of PBH binaries using the gravitational wave events from the third Gravitational-Wave Transient Catalog (GWTC-3). Considering four commonly used PBH mass functions, namely the log-normal, power-law, broken power-law, and critical collapse forms, we find that the multi-merger processes make a subdominant contribution to the total merger rate. Therefore, the effect of merger history can be safely ignored when estimating the merger rate of PBH binaries. We also find that GWTC-3 is best fitted by the log-normal form among the four PBH mass functions and confirm that the stellar-mass PBHs cannot dominate cold dark matter.

3D hydrodynamics simulations of core convection in supermassive main-sequence stars

ABSTRACT Supermassive stars are Population III stars with masses exceeding $10^4\, {\rm M}_{\odot }$ that could be the progenitors of the first supermassive black holes. Their interiors are in a regime where radiation pressure dominates the equation of state. In this work, we use the explicit gas dynamics code ppmstar to simulate the hydrogen-burning core of a $10^4\, {\rm M}_{\odot }$ supermassive main-sequence star. These are the first three-dimensional hydrodynamics simulations of core convection in supermassive stars. We perform a series of 10 simulations at different heating rates and on Cartesian grids with resolutions of 7683, 11523, and 17283. We examine different properties of the convective flow, including its large-scale morphology, its velocity spectrum, and its mixing properties. We conclude that the radiation pressure-dominated nature of the interior does not noticeably affect the behaviour of convection compared to the case of core convection in a massive main-sequence star where gas pressure dominates. Our simulations also offer support for the use of mixing-length theory in one-dimensional models of supermassive stars.

Merging binary black holes formed through double-core evolution

Context. To date, various formation channels of merging events have been heavily explored with the detection of nearly 100 double black hole (BH) merger events reported by the LIGO-Virgo-KAGRA (LVK) Collaboration. In this paper, we systematically investigate an alternative formation scenario: binary BHs (BBHs) formed through double helium stars (hereafter, “double-core evolution channel”). In this scenario, two helium stars (He-rich stars) could be the outcome of the classical isolated binary evolution scenario with and without the common envelope (CE) phase (i.e., CE channel and stable mass transfer channel) or, alternatively, of massive close binaries evolving chemically homogeneously (i.e., CHE channel). Aims. We study the properties (i.e., the chirp masses and the effective spins) of BBHs formed through the double-core evolution and investigate the impact of different efficiencies of angular momentum transport within massive He-rich stars on double-core evolution. Methods. We performed detailed stellar structure and binary evolution calculations that take into account internal rotation and mass loss of He-rich stars as well as tidal interactions in binaries. We systematically studied the parameter space of initial binary He-rich stars, including the initial mass and metallicity of He-rich stars as well as initial orbital periods. Apart from direct core collapse with mass and angular momentum conserved, we also follow the framework in Batta & Ramirez-Ruiz (2019, ArXiv e-prints [arXiv:1904.04835]) to estimate the mass and spin of the resulting BHs. Results. We show that the radii of massive He-rich stars decrease as a function of time, which comes mainly from mass loss and mixing in high metallicity and from mixing in low metallicity. For double He-rich stars with equal masses in binaries, we find that tides start to be at work on the zero age helium main sequence (i.e., the time when a He-rich star starts to burn helium in the core, which is analogous to zero age main sequence for core hydrogen burning) for initial orbital periods not longer than 1.0 day, depending on the initial metallicities. In addition to the stellar mass-loss rate and tidal interactions in binaries, we find that the role of the angular momentum transport efficiency in determining the resulting BH spins becomes stronger when considering BH progenitors originated from a higher metal-metallicity environment. We highlight that the double-core evolution scenario does not always produce fast-spinning BBHs and compare the properties of the BBHs reported from the LVK with our modeling. Conclusions. After detailed binary calculations of double-core evolution, we have confirmed that the spin of the BH is not only determined by the interplay of the binary’s different initial conditions (metallicity, mass, and orbital period) but is also dependent on the angular momentum transport efficiency within its progenitor. We predict that with the sensitivity improvements to the LVK’s next observing run (O4), the sample of merging BBHs will contain more sources with positive but moderate (even high) χeff and part of the events will likely show to have been formed through the double-core evolution channel.

POSYDON: A General-purpose Population Synthesis Code with Detailed Binary-evolution Simulations

Most massive stars are members of a binary or a higher-order stellar system, where the presence of a binary companion can decisively alter their evolution via binary interactions. Interacting binaries are also important astrophysical laboratories for the study of compact objects. Binary population synthesis studies have been used extensively over the last two decades to interpret observations of compact-object binaries and to decipher the physical processes that lead to their formation. Here, we present POSYDON, a novel, publicly available, binary population synthesis code that incorporates full stellar structure and binary-evolution modeling, using the MESA code, throughout the whole evolution of the binaries. The use of POSYDON enables the self-consistent treatment of physical processes in stellar and binary evolution, including: realistic mass-transfer calculations and assessment of stability, internal angular-momentum transport and tides, stellar core sizes, mass-transfer rates, and orbital periods. This paper describes the detailed methodology and implementation of POSYDON, including the assumed physics of stellar and binary evolution, the extensive grids of detailed single- and binary-star models, the postprocessing, classification, and interpolation methods we developed for use with the grids, and the treatment of evolutionary phases that are not based on precalculated grids. The first version of POSYDON targets binaries with massive primary stars (potential progenitors of neutron stars or black holes) at solar metallicity.