We study the impact of mass-transfer physics on the observable properties of binary black hole populations that formed through isolated binary evolution. We used the POSYDON framework to combine detailed MESA binary simulations with the COSMIC population synthesis tool to obtain an accurate estimate of merging binary black hole observables with a specific focus on the spins of the black holes. We investigate the impact of mass-accretion efficiency onto compact objects and common-envelope efficiency on the observed distributions of the effective inspiral spin parameter χeff, chirp mass Mchirp, and binary mass ratio q. We find that low common envelope efficiency translates to tighter orbits following the common envelope and therefore more tidally spun up second-born black holes. However, these systems have short merger timescales and are only marginally detectable by current gravitational-wave detectors as they form and merge at high redshifts (z ∼ 2), outside current detector horizons. Assuming Eddington-limited accretion efficiency and that the first-born black hole is formed with a negligible spin, we find that all non-zero χeff systems in the detectable population can come only from the common envelope channel as the stable mass-transfer channel cannot shrink the orbits enough for efficient tidal spin-up to take place. We find that the local rate density (z ≃ 0.01) for the common envelope channel is in the range of ∼17–113 Gpc−3 yr−1, considering a range of αCE ∈ [0.2, 5.0], while for the stable mass transfer channel the rate density is ∼25 Gpc−3 yr−1. The latter drops by two orders of magnitude if the mass accretion onto the black hole is not Eddington limited because conservative mass transfer does not shrink the orbit as efficiently as non-conservative mass transfer does. Finally, using GWTC-2 events, we constrained the lower bound of branching fraction from other formation channels in the detected population to be ∼0.2. Assuming all remaining events to be formed through either stable mass transfer or common envelope channels, we find moderate to strong evidence in favour of models with inefficient common envelopes.

Context. Ultra-luminous X-ray sources (ULXs) are those X-ray sources located away from the centre of their host galaxy with luminosities exceeding the Eddington limit of a stellar-mass black hole (LX >  1039 erg s−1). Observed X-ray variability suggests that ULXs are X-ray binary systems. The discovery of X-ray pulsations in some of these objects (e.g. M82 X-2) suggests that a certain fraction of the ULX population may have a neutron star as the accretor. Aims. We present systematic modelling of low- and intermediate-mass X-ray binaries (LMXBs and IMXBs; donor-star mass range 0.92–8.0 M⊙ and neutron-star accretors) to explain the formation of this sub-population of ULXs. Methods. Using MESA, we explored the allowed initial parameter space of binary systems consisting of a neutron star and a low- or intermediate-mass donor star that could explain the observed properties of ULXs. These donors are transferring mass at super-Eddington rates while the accretion is limited locally in the accretion disc by the Eddington limit. Thus, our simulations take into account beaming effects and also include stellar rotation, tides, general angular momentum losses, and a detailed and self-consistent calculation of the mass-transfer rate. Results. Exploring the initial parameters that lead to the formation of neutron-star ULXs, we study the conditions that lead to dynamical stability of these systems, which depends strongly on the response of the donor star to mass loss. Using two values for the initial neutron star mass (1.3 M⊙ and 2.0 M⊙), we present two sets of mass-transfer calculation grids for comparison with observations of NS ULXs. We find that LMXBs/IMXBs can produce NS-ULXs with typical time-averaged isotropic-equivalent X-ray luminosities of between 1039 and 1041 erg s−1 on a timescale of up to ∼1.0 Myr for the lower luminosities. Finally, we estimate their likelihood of detection, the types of white-dwarf remnants left behind by the donors, and the total amount of mass accreted by the neutron stars. Conclusions. We show that observed super-Eddington luminosities can be achieved in LMXBs/IMXBs undergoing non-conservative mass transfer while assuming geometrical beaming. We also compare our results to the observed pulsating ULXs and infer their initial parameters. Our results suggest that a large subset of the observed pulsating ULX population can be explained by LMXBs/IMXBs in a super-Eddington mass-transfer phase.

Long-duration gamma-ray bursts are thought to be associated with the core-collapse of massive, rapidly spinning stars and the formation of black holes. However, efficient angular momentum transport in stellar interiors, currently supported by asteroseismic and gravitational-wave constraints, leads to predominantly slowly-spinning stellar cores. Here, we report on binary stellar evolution and population synthesis calculations, showing that tidal interactions in close binaries not only can explain the observed subpopulation of spinning, merging binary black holes but also lead to long gamma-ray bursts at the time of black-hole formation. Given our model calibration against the distribution of isotropic-equivalent energies of luminous long gamma-ray bursts, we find that ≈10% of the GWTC-2 reported binary black holes had a luminous long gamma-ray burst associated with their formation, with GW190517 and GW190719 having a probability of ≈85% and ≈60%, respectively, being among them. Moreover, given an assumption about their average beaming fraction, our model predicts the rate density of long gamma-ray bursts, as a function of redshift, originating from this channel. For a constant beaming fraction fB ∼ 0.05 our model predicts a rate density comparable to the observed one, throughout the redshift range, while, at redshift z ∈ [0, 2.5], a tentative comparison with the metallicity distribution of observed LGRB host galaxies implies that between 20% to 85% of the observed long gamma-ray bursts may originate from progenitors of merging binary black holes. The proposed link between a potentially significant fraction of observed, luminous long gamma-ray bursts and the progenitors of spinning binary black-hole mergers allows us to probe the latter well outside the horizon of current-generation gravitational wave observatories, and out to cosmological distances.

No abstract

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.