Search

Current observations of binary black-hole ({BBH}) merger events show support for a feature in the primary BH-mass distribution at $\sim\,35\,\mathrm{M}_{\odot}$, previously interpreted as a signature of pulsational pair-instability (PPISN) supernovae. Such supernovae are expected to map a wide range of pre-supernova carbon-oxygen (CO) core masses to a narrow range of BH masses, producing a peak in the BH mass distribution. However, recent numerical simulations place the mass location of this peak above $50\,\mathrm{M}_{\odot}$. Motivated by uncertainties in the progenitor's evolution and explosion mechanism, we explore how modifying the distribution of BH masses resulting from PPISN affects the populations of gravitational-wave (GW) and electromagnetic (EM) transients. To this end, we simulate populations of isolated {BBH} systems and combine them with cosmic star-formation rates. Our results are the first cosmological BBH-merger predictions made using the \textsc{binary\_c} rapid population synthesis framework. We find that our fiducial model does not match the observed GW peak. We can only explain the $35\,\mathrm{M}_{\odot}$ peak with PPISNe by shifting the expected CO core-mass range for PPISN downwards by $\sim{}15\,\mathrm{M}_{\odot}$. Apart from being in tension with state-of-the art stellar models, we also find that this is likely in tension with the observed rate of hydrogen-less super-luminous supernovae. Conversely, shifting the mass range upward, based on recent stellar models, leads to a predicted third peak in the BH mass function at $\sim{}64\,\mathrm{M}_{\odot}$. Thus we conclude that the $\sim{}35\,\mathrm{M}_{\odot}$ feature is unlikely to be related to PPISNe., Comment: Accepted for publication in MNRAS. 19 pages, 8 figures includings appendices

We present the software package binary_c-python which provides a convenient and easy-to-use interface to the binary_c framework, allowing the user to rapidly evolve individual systems and populations of stars. binary_c-python is available on Pip and on GitLab. binary_c-python contains many useful features to control and process the output of binary_c, like by providing binary_c-python with logging statements that are dynamically compiled and loaded into binary_c. Moreover, we have recently added standardised output of events like Roche-lobe overflow or double compact-object formation to binary_c, and automatic parsing and managing of that output in binary_c-python. binary_c-python uses multiprocessing to utilise all the cores on a particular machine, and can run populations with HPC cluster workload managers like HTCondor and Slurm, allowing the user to run simulations on large computing clusters. We provide documentation that is automatically generated based on docstrings and a suite of Jupyter notebooks. These notebooks consist of technical tutorials on how to use binary_c-python and use-case scenarios aimed at doing science. Much of binary_c-python is covered by unit tests to ensure reliability and correctness, and the test coverage is continually increased as the package is improved., Comment: Published in the JOSS

The first population of X-ray binaries (XRBs) is expected to affect the thermal and ionization states of the gas in the early Universe. Although these X-ray sources are predicted to have important implications for high-redshift observable signals, such as the hydrogen 21-cm signal from cosmic dawn and the cosmic X-ray background, their properties are poorly explored, leaving theoretical models largely uninformed. In this paper we model a population of X-ray binaries arising from zero metallicity stars. We explore how their properties depend on the adopted initial mass function (IMF) of primordial stars, finding a strong effect on their number and X-ray production efficiency. We also present scaling relations between XRBs and their X-ray emission with the local star formation rate, which can be used in sub-grid models in numerical simulations to improve the X-ray feedback prescriptions. Specifically, we find that the uniformity and strength of the X-ray feedback in the intergalactic medium is strongly dependant on the IMF. Bottom-heavy IMFs result in a smoother distribution of XRBs, but have a luminosity orders of magnitude lower than more top-heavy IMFs. Top-heavy IMFs lead to more spatially uneven, albeit strong, X-ray emission. An intermediate IMF has a strong X-ray feedback while sustaining an even emission across the intergalactic medium. These differences in X-ray feedback could be probed in the future with measurements of the cosmic dawn 21-cm line of neutral hydrogen, which offers us a new way of constraining population III IMF., Comment: Accepted for publication in MNRAS, 17 pages, 9 figures

Population synthesis relies on semi-analytic formulae to determine masses of compact objects from the (helium or carbon-oxygen) cores of collapsing stars. Such formulae are combined across mass ranges that span different explosion mechanisms, potentially introducing artificial features in the compact object mass distribution. Such artifacts impair the interpretation of gravitational-wave observations. We propose a "top-down" remnant mass prescription where we remove mass from the star for each possible mass-loss mechanism, instead of relying on the fallback onto a "proto-compact-object" to get the final mass. For one of these mass-loss mechanisms, we fit the metallicity-dependent mass lost to pulsational-pair instability supernovae from numerical simulations. By imposing no mass loss in the absence of pulses, our approach recovers the existing compact object masses prescription at the low mass end and ensures continuity across the core-collapse/pulsational-pair-instability regime. Our remnant mass prescription can be extended to include other mass-loss mechanisms at the final collapse., Comment: to appear in RNAAS, made with showyourwork, code available at https://github.com/mathren/top_down_compact_obj_mass

The theory for single stellar evolution predicts a gap in the mass distribution of black holes (BHs) between approximately 45-130M$_{\odot}$, the so-called "pair-instability mass gap". We examine whether BHs can pollute the gap after accreting from a stellar companion. To this end, we simulate the evolution of isolated binaries using a population synthesis code, where we allow for super-Eddington accretion. Under our most extreme assumptions, we find that at most about 2% of all merging binary BH systems contains a BH with a mass in the pair-instability mass gap, and we find that less than 0.5% of the merging systems has a total mass larger than 90M$_{\odot}$. We find no merging binary BH systems with a total mass exceeding 100M$_{\odot}$. We compare our results to predictions from several dynamical pathways to pair-instability mass gap events and discuss the distinguishable features. We conclude that the classical isolated binary formation scenario will not significantly contribute to the pollution of the pair-instability mass gap. The robustness of the predicted mass gap for the isolated binary channel is promising for the prospective of placing constraints on (i) the relative contribution of different formation channels, (ii) the physics of the progenitors including nuclear reaction rates, and (iii), tentatively, the Hubble parameter., Comment: 20 pages, 9 Figures, to be published in ApJ