Showing 1–50 of 203 results for author: Belczynski, K

Unequal-mass, highly-spinning binary black hole mergers in the stable mass transfer formation channel

Authors: Aleksandra Olejak, Jakub Klencki, Xiao-Tian Xu, Chen Wang, Krzysztof Belczynski, Jean-Pierre Lasota

Abstract: The growing database of gravitational-wave (GW) detections with the binary black holes (BHs) merging in the distant Universe contains subtle insights into their formation scenarios. One of the puzzling properties of detected GW sources is the possible (anti)correlation between mass ratio q of BH-BH binaries and their effective spin. We use rapid binary evolution models to demonstrate that the isol… ▽ More The growing database of gravitational-wave (GW) detections with the binary black holes (BHs) merging in the distant Universe contains subtle insights into their formation scenarios. One of the puzzling properties of detected GW sources is the possible (anti)correlation between mass ratio q of BH-BH binaries and their effective spin. We use rapid binary evolution models to demonstrate that the isolated binary evolution followed by efficient tidal spin-up of stripped helium core produces a similar pattern in Xeff vs q distributions of BH-BH mergers. In our models, the progenitors of unequal BH-BH systems in the stable mass transfer formation scenario are more likely to efficiently shrink their orbits during the second Roche-lobe overflow than the binaries that evolve into nearly equal-mass component systems. This makes it easier for unequal-mass progenitors to enter the tidal spin-up regime and later merge due to GW emission. Our results are, however, sensitive to some input assumptions, especially, the stability of mass transfer and the angular momentum loss during non-conservative mass transfer. We note that mass transfer prescriptions widely adopted in rapid codes favor the formation of BH-BH merger progenitors with unequal masses and moderate separations. We compare our results with detailed stellar model grids and find reasonable agreement after appropriate calibration of the physics models. We anticipate that future detections of unequal-mass BH-BH mergers could provide valuable constraints on the role of the stable mass transfer formation channel. A significant fraction of BH-BH detections with mass ratio q in range (0.4 - 0.7) would be consistent with the mass ratio reversal scenario during the first, relatively conservative mass transfer, and a non-enhanced angular momentum loss during the second, highly non-conservative mass transfer phase. △ Less

Submitted 11 June, 2024; v1 submitted 18 April, 2024; originally announced April 2024.

Comments: Version accepted for publication in A&A

Evolution of rotating massive stars adopting a newer, self-consistent wind prescription at SMC metallicity

Authors: Alex Camilo Gormaz-Matamala, Jorge Cuadra, Sylvia Ekström, Georges Meynet, Michel Curé, Krzysztof Belczynski

Abstract: We use Geneva-evolution-code to run evolutionary tracks for stellar masses ranging from $20$ to $85$ $M_\odot$ at SMC metallicity ($Z=0.002$). We upgrade the recipe for stellar winds by adopting our self-consistent m-CAK prescription, which reduces the value of mass-loss rate by a factor between 2 and 6 depending on the mass range. The impact of our new winds is wide, and it can be divided betwe… ▽ More We use Geneva-evolution-code to run evolutionary tracks for stellar masses ranging from $20$ to $85$ $M_\odot$ at SMC metallicity ($Z=0.002$). We upgrade the recipe for stellar winds by adopting our self-consistent m-CAK prescription, which reduces the value of mass-loss rate by a factor between 2 and 6 depending on the mass range. The impact of our new winds is wide, and it can be divided between direct and indirect impact. For the most massive models ($60$ and $85$ $M_\odot$) with $\dot M\gtrsim2\times10^{-7}$ $M_\odot$ yr$^{-1}$, the impact is direct because lower mass loss make stars remove less envelope and therefore remain more massive and less chemically enriched at their surface at the end of their MS phase. For the less massive models ($20$ and $25$ $M_\odot$) with $\dot M\lesssim2\times10^{-8}$ $M_\odot$ yr$^{-1}$, the impact is indirect because lower mass loss make the stars keep high rotational velocities for a longer period of time, then extending the H-core burning lifetime and reaching the end of the MS with higher surface enrichment. Given that the conditions at the H-depletion change, the stars will lose more mass during their He-core burning stages anyways. For $M_\text{zams}=20$ to $40$ $M_\odot$, our models predict stars will evolve through the Hertzsprung gap, from O-type supergiants to BSG and finally RSG, with larger mass fractions of helium compared to old evolution models. New models also set down to $M_\text{zams}=85\,M_\odot$ the minimal initial mass required for a single star to become WR at metallicity $Z=0.002$. New values for $\dot M$ need to be complemented with upgrades in additional features such as convective core overshooting and distribution of rotational velocities, besides more detailed observations from projects such as XShootU, in order to provide a robust framework for the study of massive stars at low metallicity environments. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

The enigmatic origin of two dormant BH binaries: Gaia BH1 and Gaia BH2

Authors: Iwona Kotko, Sambaran Banerjee, Krzysztof Belczynski

Abstract: The two systems that have been confirmed as dormant (no X-ray emission detected) black hole (BH) - low-mass star binaries (Gaia BH1 and Gaia BH2) in the latest Gaia mission data release (DR3) are very intriguing in the context of their evolution. Both systems consist of $\sim 9 \mathrm{M_{\odot}}$ BH and $\sim 1 \mathrm{M_{\odot}}$ star orbiting each other on a wide, eccentric orbit ($e\sim 0.5$).… ▽ More The two systems that have been confirmed as dormant (no X-ray emission detected) black hole (BH) - low-mass star binaries (Gaia BH1 and Gaia BH2) in the latest Gaia mission data release (DR3) are very intriguing in the context of their evolution. Both systems consist of $\sim 9 \mathrm{M_{\odot}}$ BH and $\sim 1 \mathrm{M_{\odot}}$ star orbiting each other on a wide, eccentric orbit ($e\sim 0.5$). We argue that the formation of Gaia BH-like systems through the isolated binary evolution (IBE) channel under standard common envelope assumptions and from dynamical interactions in young massive and open clusters is equally probable, and the formation rate of those binaries is on order of $10^{-7} \mathrm{M_{\odot}}^{-1}$ in both channels. What plays an important role in formation of Gaia BH-like systems in the case of IBE channel is the mutual position of the natal kick velocity vector and the binary angular momentum vector. We find that the natal kicks of median magnitude $\sim 40$ km/s are preferred for the formation of Gaia BH1-like binaries. $\sim 94\%$ of those binaries is formed with BH spin misaligned to orbital axis by less than $40^{\circ}$. Gaia BH2-like binaries form if the low velocity natal kick (of median magnitude $\sim 20$ km/s) is directed within $15^{\circ}$ about the orbital plane. We follow the subsequent evolution of the binaries once formed as Gaia BH1 and Gaia BH2 systems to investigate their connection with the low-mass X-ray binary population. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: 15 pages, 11 figures, submitted to MNRAS

On the maximum black hole mass at solar metallicity

Authors: Amedeo Romagnolo, Alex C. Gormaz-Matamala, Krzysztof Belczynski

Abstract: In high metallicity environments the mass that black holes (BHs) can reach just after core-collapse widely depends on how much mass their progenitor stars lose via winds. On one hand new theoretical and observational insights suggest that early-stage winds should be weaker than what many canonical models prescribe. On the other hand the proximity to the Eddington limit should affect the formation… ▽ More In high metallicity environments the mass that black holes (BHs) can reach just after core-collapse widely depends on how much mass their progenitor stars lose via winds. On one hand new theoretical and observational insights suggest that early-stage winds should be weaker than what many canonical models prescribe. On the other hand the proximity to the Eddington limit should affect the formation of optically thick envelopes already during the earliest stages of stars with initial masses $M_{\rm ZAMS}\gtrsim 100$ M$_\odot$, hence resulting in higher mass-loss rates during the main sequence. We use the evolutionary codes MESA and Genec to calculate a suite of tracks for massive stars at solar metallicity Z$_\odot=0.014$ which incorporate these changes in our wind mass loss prescription. In our calculations we employ moderate rotation, high overshooting and magnetic angular momentum transport. We find a maximum BH mass $M_{\rm BH, max}=28.3$ M$_\odot$ at Z$_\odot$. The most massive BHs are predicted to form from stars with $M_{\rm ZAMS}\gtrsim 250$ M$_\odot$, with the BH mass directly proportional to its progenitor's $M_{\rm ZAMS}$. We also find in our models that at Z$_\odot$ almost any BH progenitor naturally evolves into a Wolf-Rayet star due to the combined effect of internal mixing and wind mass loss. These results are considerably different from most recent studies regarding the final mass of stars before their collapse into BHs. While we acknowledge the inherent uncertainties in stellar evolution modelling, our study underscores the importance of employing the most up-to-date physics in BH mass predictions. △ Less

Submitted 4 March, 2024; v1 submitted 30 November, 2023; originally announced November 2023.

Comments: Accepted for publication in ApJ Letters

The X-ray binaries in M83: will any of them form gravitational wave sources for LIGO/VIRGO/KAGRA?

Authors: Iwona Kotko, Krzysztof Belczynski

Abstract: There are 214 X-ray point-sources ($L_{\rm X}>10^{35} \mathrm{erg/s}$) identified as X-ray binaries (XRBs) in the nearby spiral galaxy M83. Since XRBs are powered by accretion onto a neutron star or a black hole from a companion/donor star these systems are promising progenitors of merging double compact objects (DCOs): black hole - black hole (BH-BH), black hole - neutron star (BH-NS), or neutron… ▽ More There are 214 X-ray point-sources ($L_{\rm X}>10^{35} \mathrm{erg/s}$) identified as X-ray binaries (XRBs) in the nearby spiral galaxy M83. Since XRBs are powered by accretion onto a neutron star or a black hole from a companion/donor star these systems are promising progenitors of merging double compact objects (DCOs): black hole - black hole (BH-BH), black hole - neutron star (BH-NS), or neutron star - neutron star (NS-NS) systems. The connection (i.e. XRBs evolving into DCOs) may provide some hints to the yet unanswered question: what is the origin of the LIGO/Virgo/KAGRA mergers? Available observations do not allow to determine what will be the final fate of the XRBs observed in M83. Yet, we can use evolutionary model of isolated binaries to reproduce the population of XRBs in M83 by matching model XRBs numbers/types/luminosities to observations. Knowing the detailed properties of M83 model XRBs (donor/accretor masses, their evolutionary ages and orbits) we follow their evolution to the death of donor stars to check whether any merging DCOs are formed. Although all merging DCOs in our isolated binary evolution model go through the XRB phase (defined as reaching X-ray luminosity from RLOF/wind accretion onto NS/BH above $10^{35}$ erg/s), only very few XRBs evolve to form merging (in Hubble time) DCOs. For M83 with its solar-like metallicity stars and continiuous star-formation we find that only $\sim 1-2\%$ of model XRBs evolve into merging DCOs depending on the adopted evolutionary physics. This is caused by (i) merger of donor star with compact object during common envelope phase, (ii) binary disruption at the supernova explosion of donor star, (iii) formation of a DCO on a wide orbit (merger time longer than Hubble time). △ Less

Submitted 25 October, 2023; v1 submitted 15 May, 2023; originally announced May 2023.

Comments: 17 pages, 8 figures, submitted to A&A; typos corrected, 7 figures added, sections Results and Discussion extended

Reconstructing the star formation rate for compact binary populations with the Einstein telescope

Authors: Neha Singh, Tomasz Bulik, Krzysztof Belczynski, Marek Cieslar, Francesca Calore

Abstract: The Einstein Telescope (ET) is a proposed third-generation, wide-band gravitational wave (GW) detector. Given its improved detection sensitivity in comparison to the second-generation detectors, it will be capable of exploring the Universe with GWs up to very high redshifts. In this paper, we present a population-independent method to infer the functional form of star formation rate density (SFR)… ▽ More The Einstein Telescope (ET) is a proposed third-generation, wide-band gravitational wave (GW) detector. Given its improved detection sensitivity in comparison to the second-generation detectors, it will be capable of exploring the Universe with GWs up to very high redshifts. In this paper, we present a population-independent method to infer the functional form of star formation rate density (SFR) for different populations of compact binaries originating in stars from Population (Pop) I+II and Pop III using ET as a single instrument. We use an algorithm to answer three major questions regarding the SFR of different populations of compact binaries. Specifically, these questions refer to the termination redshift of the formation of Pop III stars, the redshift at peak SFR, and the functional form of SFR at high redshift, all of which remain to be elucidated. We show that the reconstruction of SFR as a function of redshift for the different populations of compact binaries is independent of the time-delay distributions up to $z \sim 14,$ and that the accuracy of the reconstruction only strongly depends on this distribution at higher redshifts of $z\gtrsim 14$. We define the termination redshift for Pop III stars as the redshift where the SFR drops to 1\% of its peak value. In this analysis, we constrain the peak of the SFR as a function of redshift and show that ET as a single instrument can distinguish the termination redshifts of different SFRs for Pop III stars, which have a true separation of at least $Δz \sim 2$. The accurate estimation of the termination redshift depends on correctly modelling the tail of the time-delay distribution, which constitutes delay times of $\gtrsim 8$ Gyr. △ Less

Submitted 5 November, 2023; v1 submitted 3 April, 2023; originally announced April 2023.

Comments: 19 pages, 9 figures, Accepted for publication in A&A

Report number: LAPTH-016/23

Journal ref: A&A 681, A56 (2024)

Synthetic Population of Binary Cepheids. II. The effect of companion light on the extragalactic distance scale

Authors: Paulina Karczmarek, Gergely Hajdu, Grzegorz Pietrzyński, Wolfgang Gieren, Weronika Narloch, Radosław Smolec, Grzegorz Wiktorowicz, Krzysztof Belczynski

Abstract: Because of their period-luminosity relation (PLR), classical Cepheids play a key role in the calibration of the extragalactic distance scale and the determination of the Hubble-Lemaître constant $H_0$. Recent findings show that the majority of classical Cepheids should be in binary or multiple systems, which might undermine their accuracy, as the extra -- and unaccounted for -- light from the comp… ▽ More Because of their period-luminosity relation (PLR), classical Cepheids play a key role in the calibration of the extragalactic distance scale and the determination of the Hubble-Lemaître constant $H_0$. Recent findings show that the majority of classical Cepheids should be in binary or multiple systems, which might undermine their accuracy, as the extra -- and unaccounted for -- light from the companions of Cepheids causes a shift in the PLR. We quantify this shift using synthetic populations of binary Cepheids that we developed for this purpose, as described in Paper I of this series. We find that while all PLRs are shifted toward brighter values due to the excess light from the companions, the bias in the relative distance modulus between two galaxies hosting binary Cepheids can be either positive or negative, depending on the percentage of binary Cepheids in them. If the binarity percentage in the two galaxies is similar, the effect of binarity is canceled. Otherwise, it introduces a shift in the distance modulus of the order of millimags in the near-infrared passbands and Wesenheit indices, and tens of millimags in the visual domain; its exact value depends on the variant of the synthetic population (a unique combination of metallicity, star formation history, shape and location of the instability strip, and initial parameter distributions). Such shifts in distance moduli to type Ia supernova host galaxies introduce an additional statistical error on $H_0$, which however does not prevent measuring $H_0$ with a precision of 1%. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 16 pages, 11 figures, accepted for publication in The Astrophysical Journal

Symmetry breaking in merging binary black holes from young massive clusters and isolated binaries

Authors: Sambaran Banerjee, Aleksandra Olejak, Krzysztof Belczynski

Abstract: Properties of the to-date-observed binary black hole (BBH) merger events suggest a preference towards spin-orbit aligned mergers. Naturally, this has caused widespread interest and speculations regrading implications on various merger formation channels. Here we show that (i) not only the BBH-merger population from isolated binaries, but also (ii) BBH population formed in young massive clusters (Y… ▽ More Properties of the to-date-observed binary black hole (BBH) merger events suggest a preference towards spin-orbit aligned mergers. Naturally, this has caused widespread interest and speculations regrading implications on various merger formation channels. Here we show that (i) not only the BBH-merger population from isolated binaries, but also (ii) BBH population formed in young massive clusters (YMC) would possess an asymmetry in favour of aligned mergers, in the distribution of the events' effective spin parameter ($χ_{\rm eff}$). In our analysis, we utilize BBH-merger outcomes from state-of-the-art N-body evolutionary models of YMCs and isolated binary population synthesis. We incorporate, for the first time in such an analysis, misalignments due to both natal kicks and dynamical encounters. The YMC $χ_{\rm eff}$ distribution has a mean (an anti-aligned merger fraction) of $\langleχ_{\rm eff}\rangle\leq0.04$ ($f_X-\approx40\%$), which is smaller (larger) than but consistent with the observed asymmetry of $\langleχ_{\rm eff}\rangle\approx0.06$ ($f_X-\approx28\%$) as obtained from the population analysis by the LIGO-Virgo-KAGRA collaboration. In contrast, isolated binaries alone tend to produce a much stronger asymmetry; for the tested physical models, $\langleχ_{\rm eff}\rangle\approx0.25$ and $f_X-\lesssim7\%$. Although the YMC $χ_{\rm eff}$ distribution is more similar to the observed counterpart, none of the channels correctly reproduce the observed distribution. Our results suggest that further extensive model explorations for both isolated-binary and dynamical channels as well as better observational constraints are necessary to understand the physics of 'the symmetry breaking' of the BBH-merger population. △ Less

Submitted 13 June, 2023; v1 submitted 21 February, 2023; originally announced February 2023.

Comments: 20 pages including Appendix, 7 figures, 3 tables. Extended results and discussions; main conclusions are unchanged. Accepted for publication in ApJ

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

Authors: A. Romagnolo, K. Belczynski, J. Klencki, P. Agrawal, T. Shenar, D. Szécsi

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… ▽ More 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 for stars beyond 50 solar masses 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 initialised before the donor star reaches its maximum expansion. The only significant difference was found for the mass distribution of massive binary black hole mergers (total mass > 50 solar masses) 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. △ Less

Submitted 1 August, 2023; v1 submitted 28 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS

Can neutron star mergers alone explain the r-process enrichment of the Milky Way?

Authors: Chiaki Kobayashi, Ilya Mandel, Krzysztof Belczynski, Stephane Goriely, Thomas H. Janka, Oliver Just, Ashley J. Ruiter, Dany Van Beveren, Matthias U. Kruckow, Max M. Briel, Jan J. Eldridge, Elizabeth Stanway

Abstract: Comparing Galactic chemical evolution models to the observed elemental abundances in the Milky Way, we show that neutron star mergers can be a leading r-process site only if at low metallicities such mergers have very short delay times and significant ejecta masses that are facilitated by the masses of the compact objects. Namely, black hole-neutron star mergers, depending on the black-hole spins,… ▽ More Comparing Galactic chemical evolution models to the observed elemental abundances in the Milky Way, we show that neutron star mergers can be a leading r-process site only if at low metallicities such mergers have very short delay times and significant ejecta masses that are facilitated by the masses of the compact objects. Namely, black hole-neutron star mergers, depending on the black-hole spins, can play an important role in the early chemical enrichment of the Milky Way. We also show that none of the binary population synthesis models used in this paper, i.e., COMPAS, StarTrack, Brussels, ComBinE, and BPASS, can currently reproduce the elemental abundance observations. The predictions are problematic not only for neutron star mergers, but also for Type Ia supernovae, which may point to shortcomings in binary evolution models. △ Less

Submitted 18 December, 2022; v1 submitted 8 November, 2022; originally announced November 2022.

Comments: Accepted for publication in the Astrophysical Journal Letters, 15 pages, 4 figures, 1 table

Abundances and Transients from Neutron Star-White Dwarf Mergers

Authors: Mark Alexander Kaltenborn, Chris L. Fryer, Ryan T. Wollaeger, Krzysztof Belczynski, Wesley Even, Chryssa Kouveliotou

Abstract: We systematically investigate the mergers of neutron star-white dwarf binaries from beginning to end, with focus on the properties of the inflows and outflows in accretion disks and their electromagnetic emissions. Using population synthesis models, we determine a subset of these binaries in which the white dwarf companion undergoes unstable mass transfer and complete tidal disruption, forming a l… ▽ More We systematically investigate the mergers of neutron star-white dwarf binaries from beginning to end, with focus on the properties of the inflows and outflows in accretion disks and their electromagnetic emissions. Using population synthesis models, we determine a subset of these binaries in which the white dwarf companion undergoes unstable mass transfer and complete tidal disruption, forming a large accretion disk around the neutron star. The material evolves according to an advection-dominated accretion disk model with nuclear burning, neutrino-emissions, and disk-surface wind ejection. The extreme dynamics of the entire process has proven difficult for analytic analysis, and thus currently the properties are poorly understood. The outflows from the mergers are iron- and nickel-rich, giving rise to optical and infrared emissions powered from the decay of the radioactive iron-type isotopes, calculated via the SuperNu light-curve code. We find these systems capable of powering bright, yet short-lived, optical transients with the potential to power gamma-ray bursts. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: LA-UR-22-29879

The Effect of Supernova Convection On Neutron Star and Black Hole Masses

Authors: Chris L. Fryer, Aleksandra Olejak, Krzysztof Belczynski

Abstract: Our understanding of the convective-engine paradigm driving core-collapse supernovae has been used for 2 decades to predict the remnant mass distribution from stellar collapse. These predictions improve as our understanding of this engine increases. In this paper, we review our current understanding of convection (in particular, the growth rate of convection) in stellar collapse and study its effe… ▽ More Our understanding of the convective-engine paradigm driving core-collapse supernovae has been used for 2 decades to predict the remnant mass distribution from stellar collapse. These predictions improve as our understanding of this engine increases. In this paper, we review our current understanding of convection (in particular, the growth rate of convection) in stellar collapse and study its effect on the remnant mass distribution. We show how the depth of the mass gap between neutron stars and black holes can help probe this convective growth. We include a study of the effects of stochasticity in both the stellar structure and the convective seeds caused by stellar burning. We study the role of rotation and its effect on the pair-instability mass gap. Under the paradigm limiting stellar rotation to those stars in tight binaries, we determine the effect of rotation on the remnant mass distribution. △ Less

Submitted 27 April, 2022; originally announced April 2022.

Comments: 17 pages, accepted by ApJ

Report number: LA-UR-21-32205

Black hole - black hole total merger mass and the origin of LIGO/Virgo sources

Authors: K. Belczynski, Z. Doctor, M. Zevin, A. Olejak, S. Banerjee, D. Chattopadhyay

Abstract: The LIGO-Virgo-KAGRA (LVK) Collaboration has reported nearly 100 BH-BH mergers. The LVK provides estimates of rates, masses, effective spins, and redshifts for these mergers. Yet, the formation channel(s) of the mergers remains uncertain. One way to search for a formation site is to contrast properties of detected BH-BH mergers with different models of BH-BH merger formation. Our study is designed… ▽ More The LIGO-Virgo-KAGRA (LVK) Collaboration has reported nearly 100 BH-BH mergers. The LVK provides estimates of rates, masses, effective spins, and redshifts for these mergers. Yet, the formation channel(s) of the mergers remains uncertain. One way to search for a formation site is to contrast properties of detected BH-BH mergers with different models of BH-BH merger formation. Our study is designed to investigate the usefulness of the total BH-BH merger mass and its evolution with redshift in establishing the origin of gravitational-wave sources. We find that the average intrinsic BH-BH total merger mass shows exceptionally different behavior for the models that we adopt for our analysis. In the local universe (z=0) the average merger mass changes from Mtot,int=25Msun for CE binary evolution and open clusters formation channels, to Mtot,int=30Msun for the stable-RLOF binary channel, to Mtot,int=45Msun for the globular cluster channel. These differences are even more pronounced at larger redshifts. However, these differences are diminished when considering LVK O3 detector sensitivity. Comparison with LVK O3 data shows that none of our adopted models can match the data despite large errors on BH-BH masses and redshifts. We emphasize that our conclusions are derived from a small set of 6 models that are subject to numerous known uncertainties. We also note that BH-BH mergers may originate from a mix of several channels and that other (than those adopted here) BH-BH formation channels may exist. △ Less

Submitted 16 August, 2022; v1 submitted 25 April, 2022; originally announced April 2022.

Comments: ApJ: accepted. Revised and updated version

The role of supernova convection for the lower mass gap in the isolated binary formation of gravitational wave sources

Authors: Aleksandra Olejak, Chris L. Fryer, Krzysztof Belczynski, Vishal Baibhav

Abstract: Understanding astrophysical phenomena involving compact objects requires an insight about the engine behind core-collapse supernovae (SNe) and the fate of the stellar collapse of massive stars. In particular, this insight is crucial in developing an understanding of the origin and formation channels of detected population of BH-BH, BH-NS and NS-NS mergers. To gain this understanding, we must tie o… ▽ More Understanding astrophysical phenomena involving compact objects requires an insight about the engine behind core-collapse supernovae (SNe) and the fate of the stellar collapse of massive stars. In particular, this insight is crucial in developing an understanding of the origin and formation channels of detected population of BH-BH, BH-NS and NS-NS mergers. To gain this understanding, we must tie our current knowledge of pre-SN stars properties and their potential explosions to the final NS or BH mass distribution. The timescale of convection growth may have a large effect on the strength of SN explosion and therefore also on the mass distribution of stellar remnants. In this study we adopt the new formulas for the relation between the pre-SN star properties and its remnant from Fryer et al. 2022 into StarTrack population synthesis code and check how they impact double compact object (DCO) mergers formed via isolated binary evolution. The new formulas give one ability to test a wide spectrum of assumptions on the convection growth time. In particular, different variants allow for a smooth transition between having a deep lower mass gap and a remnant mass distribution filled by massive NSs and low mass BHs. In this paper we present distribution of masses, mass ratios and the local merger rate densities of DCO mergers for different variants of new remnant mass formulas. We test them together with different approaches to other highly uncertain processes. We find that mass distribution of DCO mergers up to m_1+m_2 < 35 Msun is sensitive to adopted assumption on SN convection growth timescale. Between the two extreme tested variants the probability of compact object formation within the lower mass gap may differ up to 2 orders of magnitude. The mass ratio distribution of DCO mergers is significantly influenced by SN model only for our standard mass transfer stability criteria. △ Less

Submitted 19 August, 2022; v1 submitted 19 April, 2022; originally announced April 2022.

Comments: 22 pages, version accepted for publication in MNRAS

Synthetic population of binary Cepheids. I. The effect of metallicity and initial parameter distribution on characteristics of Cepheids' companions

Authors: Paulina Karczmarek, Radosław Smolec, Gergely Hajdu, Grzegorz Pietrzyński, Wolfgang Gieren, Weronika Narloch, Grzegorz Wiktorowicz, Krzysztof Belczynski

Abstract: The majority of classical Cepheids are binary stars, yet the contribution of companions' light to the total brightness of the system has been assumed negligible and lacked a thorough, quantitative evaluation. We present an extensive study of synthetic populations of binary Cepheids, that aims to characterize Cepheids' companions (e.g. masses, evolutionary and spectral types), quantify their contri… ▽ More The majority of classical Cepheids are binary stars, yet the contribution of companions' light to the total brightness of the system has been assumed negligible and lacked a thorough, quantitative evaluation. We present an extensive study of synthetic populations of binary Cepheids, that aims to characterize Cepheids' companions (e.g. masses, evolutionary and spectral types), quantify their contribution to the brightness and color of Cepheid binaries, and assess the relevance of input parameters on the results. We introduce a collection of synthetic populations, which vary in metal content, initial parameter distribution, location of the instability strip edges, and star formation history. Our synthetic populations are free from the selection bias, while the percentage of Cepheid binaries is controlled by the binarity parameter. We successfully reproduce recent theoretical and empirical results: the percentage of binary Cepheids with main sequence (MS) companions, the contrast-mass ratio relation for binary Cepheids with MS companions, the appearance of binary Cepheids with giant evolved companions as outlier data points above the period-luminosity relation. Moreover, we present the first estimation of the percentage of binary Cepheids in the Large Magellanic Cloud and announce the quantification of the effect of binarity on the slope and zero-point of multiband period-luminosity relations, which will be reported in the next paper of this series. △ Less

Submitted 27 March, 2023; v1 submitted 1 April, 2022; originally announced April 2022.

Comments: 29 pages, 19 figures (5 in the Appendix), accepted for ApJ, replaced with a version after editor's corrections

Journal ref: The Astrophysical Journal, Volume 930, Issue 1, id.65 (2022)

Modelling the formation of the first two neutron star-black hole mergers, GW200105 and GW200115: metallicity, chirp masses and merger remnant spins

Authors: Debatri Chattopadhyay, Simon Stevenson, Floor Broekgaarden, Fabio Antonini, Krzysztof Belczynski

Abstract: The two neutron star-black hole mergers (GW200105 and GW200115) observed in gravitational waves by advanced LIGO and Virgo, mark the first ever discovery of such binaries in nature. We study these two neutron star-black hole systems through isolated binary evolution, using a grid of population synthesis models. Using both mass and spin observations (chirp mass, effective spin and remnant spin) of… ▽ More The two neutron star-black hole mergers (GW200105 and GW200115) observed in gravitational waves by advanced LIGO and Virgo, mark the first ever discovery of such binaries in nature. We study these two neutron star-black hole systems through isolated binary evolution, using a grid of population synthesis models. Using both mass and spin observations (chirp mass, effective spin and remnant spin) of the binaries, we probe their different possible formation channels in different metallicity environments. Our models only support LIGO data when assuming the black hole is non spinning. Our results show a strong preference that GW200105 and GW200115 formed from stars with sub-solar metallicities $Z\lesssim 0.005$. Only two metal-rich ($Z=0.02$) models are in agreement with GW200115. We also find that chirp mass and remnant spins jointly aid in constraining the models, whilst the effective spin parameter does not add any further information. We also present the observable (i.e. post selection effects) median values of spin and mass distribution from all our models, which maybe used as a reference for future mergers. Further, we show that the remnant spin parameter distribution exhibits distinguishable features in different neutron star-black hole sub-populations. We find that non-spinning, first born black holes dominate significantly the merging neutron star-black hole population, ensuring electromagnetic counterparts to such mergers a rare affair. △ Less

Submitted 5 May, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: accepted (MNRAS), comments welcome

Report number: Volume 513, Issue 4, Pages 5780--5789

Journal ref: Monthly Notices of the Royal Astronomical Society, July 2022

VLT/MUSE and ATCA observations of the host galaxy of the short GRB 080905A at z=0.122

Authors: A. M. Nicuesa Guelbenzu, S. Klose, P. Schady, K. Belczynski, D. H. Hartmann, L. K. Hunt, M. J. Michałowski

Abstract: Short-GRB progenitors could come in various flavors, depending on the nature of the merging compact stellar objects (including a stellar-mass black hole or not) or depending on their ages (millions or billions of years). At a redshift of z=0.122, the nearly face-on spiral host of the short GRB 080905A is one of the closest short-GRB host galaxies identified so far. This made it a preferred target… ▽ More Short-GRB progenitors could come in various flavors, depending on the nature of the merging compact stellar objects (including a stellar-mass black hole or not) or depending on their ages (millions or billions of years). At a redshift of z=0.122, the nearly face-on spiral host of the short GRB 080905A is one of the closest short-GRB host galaxies identified so far. This made it a preferred target to explore spatially resolved star-formation and to investigate the afterglow position in the context of its star formation structures. We used VLT/MUSE integral-field unit observations, supplemented by ATCA 5.5/9.0 GHz radio-continuum measurements and publicly available HST data, to study the star-formation activity in the GRB 080905A host galaxy. The MUSE observations reveal that the entire host is characterized by strong line emission. Using the Halpha line flux, we measure for the entire galaxy an SFR of about 1.6 Msun/yr, consistent with its non-detection by ATCA. Several individual star-forming regions are scattered across the host. The most luminous region has a Halpha luminosity that is nearly four times as high as the luminosity of the Tarantula nebula in the Large Magellanic Cloud. Even though star-forming activity can be traced as close to about 3 kpc (in projection) distance to the GRB explosion site, stellar population synthesis calculations show that none of the Halpha-bright star-forming regions is a likely birthplace of the short-GRB progenitor. △ Less

Submitted 9 December, 2021; originally announced December 2021.

Comments: The Astrophysical Journal Supplement Series, in press

Exploring compact binary populations with the Einstein Telescope

Authors: Neha Singh, Tomasz Bulik, Krzysztof Belczynski, Abbas Askar

Abstract: The Einstein Telescope (ET), a wide-band, future third generation gravitational wave detector, is expected to have detection rates of $\sim 10^5 - 10^6$ binary black hole (BBH) detections and $\sim 7 \times 10^4$ binary neutron star (BNS) detections in one year. The coalescence of compact binaries with a total mass of 20 - 100 $M_{\odot}$, typical of BH-BH or BH-NS binaries, will be visible up to… ▽ More The Einstein Telescope (ET), a wide-band, future third generation gravitational wave detector, is expected to have detection rates of $\sim 10^5 - 10^6$ binary black hole (BBH) detections and $\sim 7 \times 10^4$ binary neutron star (BNS) detections in one year. The coalescence of compact binaries with a total mass of 20 - 100 $M_{\odot}$, typical of BH-BH or BH-NS binaries, will be visible up to redshift $z\approx 20$ and even higher, thus facilitating the understanding of the dark era of the Universe preceding the birth of the first stars. The ET will therefore be a crucial instrument for population studies. We analysed the compact binaries originating in stars from (i) Population (Pop) I+II, (ii) Pop III, and (iii) globular clusters (GCs), with the single ET instrument, using the ET-D design sensitivity for the analysis. We estimated the constraints on the chirp mass, redshift, and merger rate as function of redshift for these classes of compact object binaries. We conclude that the ET as a single instrument is capable of detecting and distinguishing different compact binary populations separated in chirp mass - redshift space. While compact binaries originating in stars from Pop III are clearly distinguishable, owing to the separation in chirp mass - redshift space, the other two populations, Pop I+II, and GCs, can be distinguished with just 500 detections, corresponding to an observation time of $\sim 1$ hr. The mass distribution characteristics of such different compact binary populations can also be estimated with the single ET instrument. △ Less

Submitted 16 July, 2022; v1 submitted 7 December, 2021; originally announced December 2021.

Comments: 16 pages, 10 figures, Accepted in Astronomy & Astrophysics

Report number: Virgo document number VIR-0785A-21 and ET document number ET-0010A-22

Journal ref: A&A 667, A2 (2022)

Common origin of black holes in high mass X-ray binaries and in gravitation-wave sources

Authors: K. Belczynski, C. Done, S. Hagen, J. -P. Lasota, K. Sen

Abstract: BH HMXB systems are likely the progenitors of the BH-BH mergers detected in gravitational waves by LIGO/Virgo/KAGRA (LVK). Yet merging BHs reach higher masses (~ $100 M_{\odot}$) than BHs in HMXBs ($\sim 20M_{\odot}$) and have typically lower spins ($a_{\rm BH}\lesssim 0.25$ with a larger values tail) than often claimed for BH in HMXBs ($a_{\rm BH}\gtrsim 0.9$). This could suggest that these two c… ▽ More BH HMXB systems are likely the progenitors of the BH-BH mergers detected in gravitational waves by LIGO/Virgo/KAGRA (LVK). Yet merging BHs reach higher masses (~ $100 M_{\odot}$) than BHs in HMXBs ($\sim 20M_{\odot}$) and have typically lower spins ($a_{\rm BH}\lesssim 0.25$ with a larger values tail) than often claimed for BH in HMXBs ($a_{\rm BH}\gtrsim 0.9$). This could suggest that these two classes of systems belong to different populations, but we show here that this does not have to be the case. The difference in masses is easily explained as the known HMXB-BHs are in galaxies with relatively high metallicity, so their progenitor stars are subject to strong mass loss from winds, leading to relatively low mass BH at core collapse. But LVK can also detect BHs from low-metallicity galaxies that naturally produce more massive stellar-origin BHs. However, a difference in spin is more difficult to explain. Models with efficient angular momentum transport in stellar interiors produces slowly spinning progenitors of both LVK and HMXB BHs. Known HMXBs have orbital periods too long for efficient tidal spin-up, and are also unlikely to have undergone significant accretion spin-up. We show that the derived value of the BH spin depends strongly on how the HMXB accretion disc emission is modelled. We argue that since Cyg X-1 is never observed to be in a soft spectral state, the appropriate spectral models must take into account the Comptonization of the disc photosphere. We show that such models {are consistent with low spin values: $a_{\rm BH}\sim 0.1$.} This has been confirmed recently by other teams for both Cyg X-1 and LMC X-1, and we show this is also the case for M33 X-7. We conclude that all known HMXB BHs can have low spin, in accordance with the results of stellar evolution models. Hence, current observations are consistent with LVK BHs and HMXB BHs belonging to the same population. △ Less

Submitted 7 July, 2024; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: Accepted for publication in Astronomy and Astrophysics

A Nuclear Equation of State Inferred from Stellar r-process Abundances

Authors: Erika M. Holmbeck, Richard O'Shaughnessy, Vera Delfavero, Krzysztof Belczynski

Abstract: Binary neutron star mergers (NSMs) have been confirmed as one source of the heaviest observable elements made by the rapid neutron-capture (r-) process. However, modeling NSM outflows -- from the total ejecta masses to their elemental yields -- depends on the unknown nuclear equation of state (EOS) that governs neutron-star structure. In this work, we derive a phenomenological EOS by assuming that… ▽ More Binary neutron star mergers (NSMs) have been confirmed as one source of the heaviest observable elements made by the rapid neutron-capture (r-) process. However, modeling NSM outflows -- from the total ejecta masses to their elemental yields -- depends on the unknown nuclear equation of state (EOS) that governs neutron-star structure. In this work, we derive a phenomenological EOS by assuming that NSMs are the dominant sources of the heavy element material in metal-poor stars with r-process abundance patterns. We start with a population synthesis model to obtain a population of merging neutron star binaries and calculate their EOS-dependent elemental yields. Under the assumption that these mergers were responsible for the majority of r-process elements in the metal-poor stars, we find parameters representing the EOS for which the theoretical NSM yields reproduce the derived abundances from observations of metal-poor stars. For our proof-of-concept assumptions, we find an EOS that is slightly softer than, but still in agreement with, current constraints, e.g., by the Neutron Star Interior Composition Explorer, with $R_{1.4}=12.25\pm 0.03$~km and $M_{\textrm TOV}$ of $2.17\pm 0.03$~M$_\odot$(statistical uncertainties, neglecting modeling systematics). △ Less

Submitted 28 February, 2022; v1 submitted 12 October, 2021; originally announced October 2021.

Comments: 11 pages, 8 figures

The implications of high BH spins on the origin of BH-BH mergers

Authors: Aleksandra Olejak, Krzysztof Belczynski

Abstract: The LIGO/Virgo collaboration has reported 50 black hole-black hole (BH-BH) mergers and 8 candidates recovered from digging deeper into the detectors noise. The majority of these mergers have low effective spins pointing toward low BH spins and efficient angular momentum transport in massive stars as proposed by several models (e.g., the Tayler-Spruit dynamo). However, out of these 58 mergers, 7 ar… ▽ More The LIGO/Virgo collaboration has reported 50 black hole-black hole (BH-BH) mergers and 8 candidates recovered from digging deeper into the detectors noise. The majority of these mergers have low effective spins pointing toward low BH spins and efficient angular momentum transport in massive stars as proposed by several models (e.g., the Tayler-Spruit dynamo). However, out of these 58 mergers, 7 are consistent with having high effective spin parameter (chi_eff>0.3). Additionally, 2 events seem to have high effective spins sourced from the spin of the primary (more massive) BH. These particular observations could be used to discriminate between the isolated binary and dynamical formation channels. It might seem that high BH spins point to a dynamical origin if angular momentum transport in stars is efficient and forms low-spinning BHs. In such a case dynamical formation is required to produce second and third generations of BH-BH mergers with typically high-spinning BHs. Here we show, however, that isolated binary BH-BH formation also naturally reproduces such highly spinning BHs. Our models start with efficient angular momentum transport in massive stars that is needed to reproduce the majority of BH-BH mergers with low effective spins. Later, some of the binaries are subject to a tidal spin-up allowing the formation of a moderate fraction (~10%) of BH-BH mergers with high effective spins (chi_eff>=0.4-0.5). In addition, isolated binary evolution can produce a small fraction of BH-BH mergers with almost maximally spinning primary BHs. Therefore, the formation scenario of these atypical BH-BH mergers remains to be found. △ Less

Submitted 11 October, 2021; v1 submitted 14 September, 2021; originally announced September 2021.

Comments: The manuscript version after the first review in ApJL. Comments welcome

The Uncertain Future of Massive Binaries Obscures the Origin of LIGO/Virgo Sources

Authors: K. Belczynski, A. Romagnolo, A. Olejak, J. Klencki, D. Chattopadhyay, S. Stevenson, M. Coleman Miller, J. -P. Lasota, P. A. Crowther

Abstract: The LIGO/Virgo gravitational--wave observatories have detected 50 BH-BH coalescences. This sample is large enough to have allowed several recent studies to draw conclusions about the 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 black holes.… ▽ More The LIGO/Virgo gravitational--wave observatories have detected 50 BH-BH coalescences. This sample is large enough to have allowed several recent studies to draw conclusions about the 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 black holes. 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 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. We show that, despite its 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-Zytkow object, to a complete disruption of both objects by pair-instability supernovae. Thus since the future of massive binaries is inherently uncertain, sound predictions about the properties of BH-BH systems are highly challenging at this time. Consequently, drawing conclusions about the formation channels for the LIGO/Virgo BH-BH merger population is premature. △ Less

Submitted 24 August, 2021; originally announced August 2021.

Comments: comments welcomed

Testing the presence of a dormant black hole inside HR 6819

Authors: A. Romagnolo, A. Olejak, A. Hypki, G. Wiktorowicz, K. Belczynski

Abstract: HR 6819 was recently reported to be a triple system with a dormant black hole (BH). The inner binary system was defined as a star estimated to be at the end of its main sequence and a dormant BH. As the inner binary is not resolved, the third component may actually just be spatially coinciding with the inner binary. In this study we test whether the system inner binary can be reconstructed using t… ▽ More HR 6819 was recently reported to be a triple system with a dormant black hole (BH). The inner binary system was defined as a star estimated to be at the end of its main sequence and a dormant BH. As the inner binary is not resolved, the third component may actually just be spatially coinciding with the inner binary. In this study we test whether the system inner binary can be reconstructed using the isolated binary evolution in the Galactic field or through the dynamical evolution within globular star clusters. Our goal is to understand the formation of the HR 6819 inner binary. To simulate the inner binary evolution we assumed that the influence of the third body on the inner binary is negligible. We created synthetic populations of BH main-sequence binaries for the Galactic disc and the Galactic globular clusters to compare to the reported parameters of the HR 6819 inner binary. We have adopted very optimistic input physics, in terms of common envelope evolution and BH formation, for the formation of binaries similar to the reported inner HR 6819 binary. Despite our optimistic assumptions we cannot form systems like the inner HR 6819 binary in globular clusters. Even with our extreme assumptions, the formation of an HR 6819-like binary in the Galactic field population is not expected. We argue that if a dormant BH actually exists in the reported configuration inside HR 6819, its presence cannot easily be explained by our models based on isolated and dynamical binary evolution. △ Less

Submitted 9 September, 2022; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 11 pages

Journal ref: A&A 667, A55 (2022)

Wind-powered ultraluminous X-ray sources

Authors: Grzegorz Wiktorowicz, Jean-Pierre Lasota, Krzysztof Belczynski, Youjun Lu, Jifeng Liu, Krystian Iłkiewicz

Abstract: Although ultra-luminous X-ray sources (ULX) are important for astrophysics due to their extreme apparent super-Eddington luminosities, their nature is still poorly known. Theoretical and observational studies suggest that ULXs could be a diversified group of objects composed of low-mass X-ray binaries, high-mass X-ray binaries and marginally also systems containing intermediate-mass black holes, w… ▽ More Although ultra-luminous X-ray sources (ULX) are important for astrophysics due to their extreme apparent super-Eddington luminosities, their nature is still poorly known. Theoretical and observational studies suggest that ULXs could be a diversified group of objects composed of low-mass X-ray binaries, high-mass X-ray binaries and marginally also systems containing intermediate-mass black holes, which is supported by their presence in a variety of environments. Observational data on the ULX donors could significantly boost our understanding of these systems, but only a few were detected. There are several candidates, mostly red supergiants (RSGs), but surveys are typically biased toward luminous near-infrared objects. Nevertheless, it is worth exploring if RSGs can be members of ULX binaries. In such systems matter accreted onto the compact body would have to be provided by the stellar wind of the companion, since a Roche-lobe overflow could be unstable for relevant mass-ratios. Here we present a comprehensive study of the evolution and population of wind-fed ULXs and provide a theoretical support for the link between RSGs and ULXs. Our estimated upper limit on contribution of wind-fed ULX to the overall ULX population is $\sim75$--$96\%$ for young ($<100$ Myr) star forming environments, $\sim 49$--$87\%$ for prolonged constant star formation (e.g., disk of Milky Way), and $\lesssim1\%$ for environments in which star formation ceased long time ($>2$ Gyr) ago. We show also that some wind-fed ULXs (up to $6\%$) may evolve into merging double compact objects (DCOs), but typical systems are not viable progenitors of such binaries because of their large separations. We demonstrate that, the exclusion of wind-fed ULXs from population studies of ULXs, might have lead to systematical errors in their conclusions. △ Less

Submitted 17 August, 2021; v1 submitted 2 March, 2021; originally announced March 2021.

Comments: 14 pages, 9 figures, 2 tables, accepted to ApJ

The impact of common envelope development criteria on the formation of LIGO/Virgo sources

Authors: Aleksandra Olejak, Krzysztof Belczynski, Natalia Ivanova

Abstract: The treatment and criteria for development of unstable Roche lobe overflow (RLOF) that leads to the common envelope (CE) phase have hindered the 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 e… ▽ More The treatment and criteria for development of unstable Roche lobe overflow (RLOF) that leads to the common envelope (CE) phase have hindered the 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. At the moment, there is no full first-principles model for development of CE. We employ the Startrack population synthesis code to test the current advancements in studies on stability of RLOF for massive donors to assess their effect on LIGO/Virgo source population. In particular, we allow for more restrictive CE development criteria for massive donors. We also test a modified condition for switching between different types of stable mass transfer. Implemented modifications significantly influence basic properties of merging double compact objects, sometimes in non-intuitive way. For one of tested models with restricted CE development criteria local merger rate density for BH-BH systems increased due to 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) 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, one is not able to draw fully reliable conclusions about the population of double compact object systems based on population synthesis studies. △ Less

Submitted 10 May, 2021; v1 submitted 10 February, 2021; originally announced February 2021.

Comments: 20 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 651, A100 (2021)

A supernova remnant associated with a nascent black hole low-mass X-ray binary

Authors: N. I. Maxted, A. J. Ruiter, K. Belczynski, I. R. Seitenzahl, R. M. Crocker

Abstract: Energy released when the core of a high-mass star collapses into a black hole often powers an explosion that creates a supernova remnant. Black holes have limited windows of observability, and consequently are rarely identified in association with supernova remnants. Analysing multi-messenger data, we show that MAXI J1535-571 is the black hole produced in the stellar explosion that gave rise to th… ▽ More Energy released when the core of a high-mass star collapses into a black hole often powers an explosion that creates a supernova remnant. Black holes have limited windows of observability, and consequently are rarely identified in association with supernova remnants. Analysing multi-messenger data, we show that MAXI J1535-571 is the black hole produced in the stellar explosion that gave rise to the supernova remnant G323.7-1.0, making it the first case of an association between a black hole low-mass X-ray binary and a supernova remnant. Given this connection, we can infer from our modelling that the progenitor system was a close binary whose primary star had an initial mass of approx. 23-35 solar masses with a companion star about 10 times less massive. △ Less

Submitted 28 October, 2020; originally announced October 2020.

Comments: 4 pages text + 2 figures (8 pages total). To be submitted. Comments and feedback welcome!

The most ordinary formation of the most unusual double black hole merger

Authors: Krzysztof Belczynski

Abstract: LIGO/Virgo Collaboration reported the detection of the most massive black hole - black hole (BH-BH) merger up to date with component masses of 85 Msun and 66 Msun (GW190521). Motivated by recent observations of massive stars in the 30 Doradus cluster in the Large Magellanic Cloud (>200 Msun; e.g. R136a) and employing newly estimated uncertainties on pulsational pair-instability mass-loss (that all… ▽ More LIGO/Virgo Collaboration reported the detection of the most massive black hole - black hole (BH-BH) merger up to date with component masses of 85 Msun and 66 Msun (GW190521). Motivated by recent observations of massive stars in the 30 Doradus cluster in the Large Magellanic Cloud (>200 Msun; e.g. R136a) and employing newly estimated uncertainties on pulsational pair-instability mass-loss (that allow for possibility of forming BHs with mass up to 90Msun) we show that it is trivial to form such massive BH-BH mergers through the classical isolated binary evolution (with no assistance from either dynamical interactions or exotica). A binary consisting of two massive (180+150 Msun) Population II stars (Z=0.0001) evolves through a stable Roche lobe overflow and common envelope episode. Both exposed stellar cores undergo direct core-collapse and form massive BHs while avoiding pair-instability pulsation mass-loss or total disruption. LIGO/Virgo observations show that the merger rate density of light BH-BH mergers (both components: <50 Msun) is of the order of 10-100 Gpc^-3 yr^-1, while GW190521 indicates that the rate of heavier mergers is 0.02-0.43 Gpc^-3 yr^-1. Our model (with standard assumptions about input physics) but extended to include 200 Msun stars and allowing for the possibility of stellar cores collapsing to 90 Msun BHs produces the following rates: 63 Gpc^-3 yr^-1 for light BH-BH mergers and 0.04 Gpc^-3 yr^-1 for heavy BH-BH mergers. We do not claim that GW190521 was formed by an isolated binary, but it appears that such a possibility can not be excluded. △ Less

Submitted 13 November, 2020; v1 submitted 28 September, 2020; originally announced September 2020.

Comments: after revision: now only 2 sigma deviation in C12 fusion rate required

Black hole - neutron star mergers: The first mass gap and kilonovae

Authors: P. Drozda, K. Belczynski, R. O'Shaughnessy, T. Bulik, C. L. Fryer

Abstract: Observations of X-ray binaries indicate a dearth of compact objects in the mass range from $\sim 2-5 M_{\odot}$. The existence of this (first mass) gap has been used to discriminate between proposed engines behind core-collapse supernovae. From LIGO/Virgo observations of binary compact remnant masses, several candidate first mass gap objects (either neutron stars (NSs) or black holes (BHs)) were i… ▽ More Observations of X-ray binaries indicate a dearth of compact objects in the mass range from $\sim 2-5 M_{\odot}$. The existence of this (first mass) gap has been used to discriminate between proposed engines behind core-collapse supernovae. From LIGO/Virgo observations of binary compact remnant masses, several candidate first mass gap objects (either neutron stars (NSs) or black holes (BHs)) were identified during the O3 science run. Motivated by these new observations, we study the formation of BH-NS mergers in the framework of isolated classical binary evolution, using population synthesis methods to evolve large populations of binary stars (Population I and II) across cosmic time. We present results on the NS to BH mass ratios ($q=M_{\rm NS}/M_{\rm BH}$) in merging systems, showing that although systems with a mass ratio as low as $q=0.02$ can exist, typically BH-NS systems form with moderate mass ratios $q=0.1-0.2$. If we adopt a delayed supernova engine, we conclude that $\sim 30\%$ of BH-NS mergers may host at least one compact object in the first mass gap (FMG$^\circ$). Even allowing for uncertainties in the processes behind compact object formation, we expect the fraction of BH-NS systems ejecting mass during the merger to be small (from $\sim 0.6-9\%$). In our reference model, we assume: (i) the formation of compact objects within the FMG, (ii) natal NS/BH kicks decreased by fallback, (iii) low BH spins due to Tayler-Spruit angular momentum transport in massive stars. We find that $\lesssim 1\%$ of BH-NS mergers will have any mass ejection and about the same percentage will produce kilonova bright enough to have a chance of being detected with a large (Subaru-class) $8$m telescope. Interestingly, all these mergers will have both a BH and an NS in the FMG. △ Less

Submitted 29 December, 2022; v1 submitted 14 September, 2020; originally announced September 2020.

Comments: 15 pages, 12 figures

Journal ref: Astronomy & Astrophysics, Volume 667, (2022), id.A126, 15 pp

Non-interacting black hole binaries with Gaia and LAMOST

Authors: Grzegorz Wiktorowicz, Youjun Lu, Łukasz Wyrzykowski, Haotong Zhang, Jifeng Liu, Stephen Justham, Krzysztof Belczynski

Abstract: Until recently, black holes (BHs) could be discovered only through accretion from other stars in X-ray binaries, or in merging double compact objects. Improvements in astrometric and spectroscopic measurements have made it possible to detect BHs also in non-interacting BH binaries (nBHB) through a precise analysis of the companion's motion. In this study, using an updated version of the Startrack… ▽ More Until recently, black holes (BHs) could be discovered only through accretion from other stars in X-ray binaries, or in merging double compact objects. Improvements in astrometric and spectroscopic measurements have made it possible to detect BHs also in non-interacting BH binaries (nBHB) through a precise analysis of the companion's motion. In this study, using an updated version of the Startrack binary-star population modelling code and a detailed model of the Milky Way (MW) galaxy we calculate the expected number of detections for Gaia and LAMOST surveys. We develop a formalism to convolve the binary population synthesis output with a realistic stellar density distribution, star-formation history (SFH), and chemical evolution for the MW, which produces a probability distribution function of the predicted compact-binary population over the MW. This avoids the additional statistical uncertainty which is introduced by methods which Monte Carlo sample from binary population synthesis output to produce one potential specific realisation of the MW compact-binary distribution, and our method is also comparatively fast to such Monte Carlo realisations. Specifically, we predict $\sim41$-$340$ nBHBs to be observed by Gaia, although the numbers may drop to $\sim10$-$70$ if the recent ($\lesssim100\;$ Myr) star formation is low ($\sim1\;M_\odot$/yr ). For LAMOST we predict $\lesssim14$ detectable nBHBs, which is lower partially because its field-of-view covers just $\sim6\%$ of the Galaxy. △ Less

Submitted 29 December, 2020; v1 submitted 15 June, 2020; originally announced June 2020.

Comments: 23 pages, 15 figures, 8 tables

Journal ref: 2020 ApJ 905 134

Phase-dependent study of near-infrared disk emission lines in LB-1

Authors: Jifeng Liu, Zheng Zheng, Roberto Soria, Jesus Aceituno, Haotong Zhang, Youjun Lu, Song Wang, Wolf-Rainer Hamann, Lida M. Oskinova, Varsha Ramachandran, Hailong Yuan, Zhongrui Bai, Shu Wang, Brendan J. McKee, Jianfeng Wu, Junfeng Wang, Mario Lattanzi, Krzysztof Belczynski, Jorge Casares, Sergio Simon-Diaz, Jonay I. González Hernández, Rafael Rebolo

Abstract: The mass, origin and evolutionary stage of the binary system LB-1 has been the subject of intense debate, following the claim that it hosts an $\sim$70$M_{\odot}$ black hole, in stark contrast with the expectations for stellar remnants in the Milky Way. We conducted a high-resolution, phase-resolved spectroscopic study of the near-infrared Paschen lines in this system, using the 3.5-m telescope at… ▽ More The mass, origin and evolutionary stage of the binary system LB-1 has been the subject of intense debate, following the claim that it hosts an $\sim$70$M_{\odot}$ black hole, in stark contrast with the expectations for stellar remnants in the Milky Way. We conducted a high-resolution, phase-resolved spectroscopic study of the near-infrared Paschen lines in this system, using the 3.5-m telescope at Calar Alto Observatory. We find that Pa$β$ and Pa$γ$ (after proper subtraction of the stellar absorption component) are well fitted with a standard double-peaked model, typical of disk emission. We measured the velocity shifts of the red and blue peaks at 28 orbital phases: the line center has an orbital motion in perfect antiphase with the stellar motion, and the radial velocity amplitude ranges from 8 to 13 km/s for different choices of lines and profile modelling. We interpret this curve as proof that the disk is tracing the orbital motion of the primary, ruling out the circumbinary disk and the hierarchical triple scenarios. The phase-averaged peak-to-peak half-separation (proxy for the projected rotational velocity of the outer disk) is $\sim$70 km s$^{-1}$, larger than the stellar orbital velocity and also inconsistent with a circumbinary disk. From those results, we infer a primary mass 4--8 times higher than the secondary mass. Moreover, we show that the ratio of the blue and red peaks (V/R intensity ratio) has a sinusoidal behaviour in phase with the secondary star, which can be interpreted as the effect of external irradiation by the secondary star on the outer disk. Finally, we briefly discuss our findings in the context of alternative scenarios recently proposed for LB-1. Definitive tests between alternative solutions will require further astrometric data from $Gaia$. △ Less

Submitted 26 May, 2020; originally announced May 2020.

Comments: To be submitted to ApJ. Comments are welcome

The Origin of inequality: isolated formation of a 30+10Msun binary black-hole merger

Authors: A. Olejak, M. Fishbach, K. Belczynski, D. E. Holz, J. -P. Lasota, M. C. Miller, T. Bulik

Abstract: The LIGO/Virgo collaboration has reported the detection of GW190412, a BH-BH merger with the most unequal masses to date: 24.4-34.7 Msun and 7.4-10.1 Msun (a mass ratio of q=0.21-0.41). Additionally, GW190412's effective spin was estimated to be Xeff=0.14-0.34, with the spin of the primary BH in the range a=0.17-0.59. Based on this and prior detections, about 10 percent of BH-BH mergers have q<0.4… ▽ More The LIGO/Virgo collaboration has reported the detection of GW190412, a BH-BH merger with the most unequal masses to date: 24.4-34.7 Msun and 7.4-10.1 Msun (a mass ratio of q=0.21-0.41). Additionally, GW190412's effective spin was estimated to be Xeff=0.14-0.34, with the spin of the primary BH in the range a=0.17-0.59. Based on this and prior detections, about 10 percent of BH-BH mergers have q<0.4. Major BH-BH formation channels tend to produce BH-BH mergers with comparable masses (typically with q>0.5). Here we test whether the classical isolated binary evolution channel can produce mergers resembling GW190412. We show that our standard binary evolution scenario, with the typical assumptions on input physics we have used in the past, produces such mergers (masses and spins). For this particular model of the input physics the overall BH-BH merger rate density in the local Universe (z=0) is: 73.5 Gpc^-3 yr^-1, while for systems with q<0.41 the rate density is: 6.8 Gpc^-3 yr^-1. As GW190412 shows some weak evidence for misaligned spins, we provide distribution of precession parameter in our models and conclude that if among the new LIGO/Virgo detections the evidence of system precession is strong and more than 10 percent of BH-BH mergers have large in-plane spin components (Xp>0.5) then common envelope isolated binary BH-BH formation channel can be excluded as their origin. △ Less

Submitted 5 September, 2020; v1 submitted 24 April, 2020; originally announced April 2020.

Comments: ApJ Letters: accepted

Formation of low-spinning 100 Msun black holes

Authors: Krzysztof Belczynski, Sambaran Banerjee

Abstract: It is speculated that a merger of two massive stellar-origin BHs in a dense stellar environment may lead to the formation of a massive BH in the pair-instability mass gap (50-135 Msun). Such a merger-formed BH is expected to typically have a high spin (a=0.7). If such a massive BH acquires another BH it may lead to another merger detectable by LIGO/Virgo in gravitational waves. Acquiring a compani… ▽ More It is speculated that a merger of two massive stellar-origin BHs in a dense stellar environment may lead to the formation of a massive BH in the pair-instability mass gap (50-135 Msun). Such a merger-formed BH is expected to typically have a high spin (a=0.7). If such a massive BH acquires another BH 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 BH from its parent globular or nuclear cluster. We test whether it is possible for a massive merger-formed BH in the pair-instability gap to be retained in its parent cluster and have low spin. Such a BH would be indistinguishable from a primordial BH. 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 BH masses and spins. We compared merger-formed massive BH speeds with typical escape velocities from globular and nuclear clusters. We show that a globular cluster is highly unlikely to form and retain a 100 Msun BH if the spin of the BH is low (a<0.3) as such BHs acquire high recoil speeds (>200 km/s) that exceed typical escape speeds from globular clusters (50 km/s). However, a very low-spinning (a=0.1) and massive (100 Msun) BH could be formed and retained in a galactic nuclear star cluster. Even though such massive merger-formed BHs with such low spins acquire high speeds during formation (400 km/s), they may avoid ejection since massive nuclear clusters have high escape velocities (300-500 km/s). A future detection of a massive BH in the pair-instability mass gap with low spin would therefore not be proof of the existence of primordial BHs, which are sometimes claimed to have low spins and arbitrarily high masses. △ Less

Submitted 23 August, 2020; v1 submitted 19 February, 2020; originally announced February 2020.

Comments: published: A&A 640, L20 (2020)

Orbital period dependent spatial scatter of Galactic Black Hole Transients

Authors: P. Gandhi, A. Rao, P. A. Charles, K. Belczynski, T. J. Maccarone, K. Arur, J. M. Corral-Santana

Abstract: There remain significant uncertainties in the origin and evolution of black holes in binary systems, in particular regarding their birth sites and the influence of natal kicks. These are long-standing issues, but their debate has been reinvigorated in the era of gravitational wave detections and the improving precision of astrometric measurements. Using recent and archival characterisation of Gala… ▽ More There remain significant uncertainties in the origin and evolution of black holes in binary systems, in particular regarding their birth sites and the influence of natal kicks. These are long-standing issues, but their debate has been reinvigorated in the era of gravitational wave detections and the improving precision of astrometric measurements. Using recent and archival characterisation of Galactic black hole X-ray binaries (BHXBs), we report here an apparent anticorrelation between P{orb} (system orbital periods) and scatter in z (elevation above the Galactic plane). The absence of long period sources at high z is not an obvious observational bias, and two possible explanatory scenarios are qualitatively explored: (1) a disc origin for BHXBs followed by natal kicks producing the scatter in z, with only the tightest binaries preferentially surviving strong kicks; (2) a halo origin, with P{orb} shortening through dynamical interactions in globular clusters (GCs). For the latter case, we show a correspondence in z-scatter between BHXBs and the GCs with most compact core radii of <0.1pc. However, the known absence of outbursting BHXB transients within Galactic GCs remains puzzling in this case, in contrast to the multitude of known GC neutron star XRBs. These results provide an interesting observational constraint for any black hole binary evolutionary model to satisfy. △ Less

Submitted 30 April, 2020; v1 submitted 3 February, 2020; originally announced February 2020.

Comments: MNRAS Letters in press following minor revision. Includes Online Appendix table

The Formation of a 70 Msun Black Hole at High Metallicity

Authors: K. Belczynski, R. Hirschi, E. A. Kaiser, Jifeng Liu, J. Casares, Youjun Lu, R. O'Shaughnessy, A. Heger, S. Justham, R. Soria

Abstract: A 70Msun BH was discovered in Milky Way disk in a long period and almost circular detached binary system (LB-1) with a high metallicity 8Msun B star companion. Current consensus on the formation of BHs from high metallicity stars limits the black hole mass to be below 20Msun. Using simple evolutionary model, we show that the formation of a 70Msun BH in high metallicity environment is possible if s… ▽ More A 70Msun BH was discovered in Milky Way disk in a long period and almost circular detached binary system (LB-1) with a high metallicity 8Msun B star companion. Current consensus on the formation of BHs from high metallicity stars limits the black hole mass to be below 20Msun. Using simple evolutionary model, we show that the formation of a 70Msun BH in high metallicity environment is possible if stellar wind mass loss rates are reduced by factor of 5. As observations indicate, a fraction of massive stars have surface magnetic fields which may quench the wind mass-loss, independently of stellar mass and metallicity. We also computed detailed stellar evolution models and we confirm such a scenario. A non-rotating 85Msun model at Z=0.014 with decreased winds ends up as a 71Msun star prior core-collapse with a 32Msun helium core and a 28Msun CO core. Such star avoids pair-instability pulsation supernova mass loss and may form a 70Msun BH in the direct collapse. Stars that can form such BHs expand to significant size with radius of R>600Rsun, exceeding the size of LB-1 orbit. Therefore, we can explain the formation of BHs upto 70Msun at high metallicity and this result is independent from LB-1. However, if LB-1 hosts a massive BH we are unable to explain how such a binary star system could have formed without invoking some exotic scenarios. △ Less

Submitted 11 January, 2020; v1 submitted 27 November, 2019; originally announced November 2019.

Comments: ApJ: after 2nd review

A wide star-black-hole binary system from radial-velocity measurements

Authors: Jifeng Liu, Haotong Zhang, Andrew W. Howard, Zhongrui Bai, Youjun Lu, Roberto Soria, Stephen Justham, Xiangdong Li, Zheng Zheng, Tinggui Wang, Krzysztof Belczynski, Jorge Casares, Wei Zhang, Hailong Yuan, Yiqiao Dong, Yajuan Lei, Howard Isaacson, Song Wang, Yu Bai, Yong Shao, Qing Gao, Yilun Wang, Zexi Niu, Kaiming Cui, Chuanjie Zheng , et al. (30 additional authors not shown)

Abstract: All stellar mass black holes have hitherto been identified by X-rays emitted by gas that is accreting onto the black hole from a companion star. These systems are all binaries with black holes below 30 M$_{\odot}$$^{1-4}$. Theory predicts, however, that X-ray emitting systems form a minority of the total population of star-black hole binaries$^{5,6}$. When the black hole is not accreting gas, it c… ▽ More All stellar mass black holes have hitherto been identified by X-rays emitted by gas that is accreting onto the black hole from a companion star. These systems are all binaries with black holes below 30 M$_{\odot}$$^{1-4}$. Theory predicts, however, that X-ray emitting systems form a minority of the total population of star-black hole binaries$^{5,6}$. When the black hole is not accreting gas, it can be found through radial velocity measurements of the motion of the companion star. Here we report radial velocity measurements of a Galactic star, LB-1, which is a B-type star, taken over two years. We find that the motion of the B-star and an accompanying H$α$ emission line require the presence of a dark companion with a mass of $68^{+11}_{-13}$ M$_{\odot}$, which can only be a black hole. The long orbital period of 78.9 days shows that this is a wide binary system. The gravitational wave experiments have detected similarly massive black holes$^{7,8}$, but forming such massive ones in a high-metallicity environment would be extremely challenging to current stellar evolution theories$^{9-11}$. △ Less

Submitted 27 November, 2019; originally announced November 2019.

Comments: Published in Nature on Nov 28, 2019

Journal ref: Nature, 575, 618-621 (2019)

The Impact of Nuclear Physics Uncertainties on Galactic Chemical Evolution Predictions

Authors: Benoit Côté, Pavel Denissenkov, Falk Herwig, Chris L. Fryer, Krzysztof Belczynski, Nicole Vassh, Matthew R. Mumpower, Jonas Lippuner, Marco Pignatari, Ashley J. Ruiter

Abstract: Modeling the evolution of the elements in the Milky Way is a multidisciplinary and challenging task. In addition to simulating the 13 billion years evolution of our Galaxy, chemical evolution simulations must keep track of the elements synthesized and ejected from every astrophysical site of interest (e.g., supernova, compact binary merger). The elemental abundances of such ejecta, which are a fun… ▽ More Modeling the evolution of the elements in the Milky Way is a multidisciplinary and challenging task. In addition to simulating the 13 billion years evolution of our Galaxy, chemical evolution simulations must keep track of the elements synthesized and ejected from every astrophysical site of interest (e.g., supernova, compact binary merger). The elemental abundances of such ejecta, which are a fundamental input for chemical evolution codes, are usually taken from theoretical nucleosynthesis calculations performed by the nuclear astrophysics community. Therefore, almost all chemical evolution predictions rely on the nuclear physics behind those calculations. In this proceedings, we highlight the impact of nuclear physics uncertainties on galactic chemical evolution predictions. We demonstrate that nuclear physics and galactic evolution uncertainties both have a significant impact on interpreting the origin of neutron-capture elements in our Solar System. Those results serve as a motivation to create and maintain collaborations between the fields of nuclear astrophysics and galaxy evolution. △ Less

Submitted 31 October, 2019; originally announced November 2019.

Comments: 8 pages, 2 figures, NPA-IX proceedings

The connection between merging double compact objects and the Ultraluminous X-ray Sources

Authors: Samaresh Mondal, Krzysztof Belczyński, Grzegorz Wiktorowicz, Jean-Pierre Lasota, Andrew R. King

Abstract: We explore the different formation channels of merging double compact objects (DCOs: BH-BH/BH-NS/NS-NS) that went through a ultraluminous X-ray phase (ULX: X-ray sources with apparent luminosity exceeding $10^{39}$ erg s$^{-1}$). There are many evolutionary scenarios which can naturally explain the formation of merging DCO systems: isolated binary evolution, dynamical evolution inside dense cluste… ▽ More We explore the different formation channels of merging double compact objects (DCOs: BH-BH/BH-NS/NS-NS) that went through a ultraluminous X-ray phase (ULX: X-ray sources with apparent luminosity exceeding $10^{39}$ erg s$^{-1}$). There are many evolutionary scenarios which can naturally explain the formation of merging DCO systems: isolated binary evolution, dynamical evolution inside dense clusters and chemically homogeneous evolution of field binaries. It is not clear which scenario is responsible for the majority of LIGO/Virgo sources. Finding connections between ULXs and DCOs can potentially point to the origin of merging DCOs as more and more ULXs are discovered. We use the StarTrack population synthesis code to show how many ULXs will form merging DCOs in the framework of isolated binary evolution. Our merger rate calculation shows that in the local Universe typically 50% of merging BH-BH progenitor binaries have evolved through a ULX phase. This indicates that ULXs can be used to study the origin of LIGO/Virgo sources. We have also estimated that the fraction of observed ULXs that will form merging DCOs in future varies between 5% to 40% depending on common envelope model and metallicity. △ Less

Submitted 28 November, 2019; v1 submitted 10 September, 2019; originally announced September 2019.

Comments: 15 pages, 7 figures, Accepted in MNRAS

Synthetic catalog of black holes in the Milky Way

Authors: A. Olejak, K. Belczynski, T. Bulik, M. Sobolewska

Abstract: We present an open-access database which includes a synthetic catalog of black holes in the Milky Way. To calculate evolution of single and binary stars we used updated population synthesis code StarTrack. We applied a new model of star formation history and chemical evolution of Galactic disk, bulge and halo synthesized from observational and theoretical data. We find that at the current moment M… ▽ More We present an open-access database which includes a synthetic catalog of black holes in the Milky Way. To calculate evolution of single and binary stars we used updated population synthesis code StarTrack. We applied a new model of star formation history and chemical evolution of Galactic disk, bulge and halo synthesized from observational and theoretical data. We find that at the current moment Milky Way (disk+bulge+halo) contains about 1.2 x 10^8 single black holes with average mass of about 14 Msun and 9.3 x 10^6 BHs in binary systems with average mass of 19 Msun. We present basic statistical properties of BH populations such as distributions of single and binary BH masses, velocities, orbital parameters or numbers of BH binary systems in different evolutionary configurations. We find that the most massive BHs are formed in mergers of binary systems, such as BH-MS, BH+He, BH-BH. The metallicity of stellar population has a significant impact on the final BH mass due to the stellar winds. Therefore the most massive single BH in our simulation, 113 Msun, originates from a merger of a helium star and a black hole in a low metallicity stellar environment in Galactic halo. The most massive BH in binary system is 60 Msun and was also formed in Galactic halo. We constrain that only 0.006% of total Galactic halo mass (including dark matter) could be hidden in the form of stellar origin BHs which are not detectable by current observational surveys. Galactic binary BHs are minority (10% of all Galactic BHs) and most of them are in BH-BH systems. The current Galactic merger rates for two considered common envelope models which are: 3-81 Myr^-1 for BH-BH, 1-9 Myr^-1, for BH-NS and 14-59 Myr^-1 for NS-NS systems. Data files are available at https://bhc.syntheticuniverse.org/. △ Less

Submitted 7 April, 2020; v1 submitted 23 August, 2019; originally announced August 2019.

Comments: 21 pages, A&A accepted, data from catalog available online

Journal ref: A&A 638, A94 (2020)

Populations of stellar mass Black holes from binary systems

Authors: Grzegorz Wiktorowicz, Łukasz Wyrzykowski, Martyna Chruslinska, Jakub Klencki, Krzysztof A. Rybicki, Krzysztof Belczynski

Abstract: In large and complicated stellar systems like galaxies it is difficult to predict the number and characteristics of a black hole population. Such populations may be modelled as an aggregation of homogeneous (i.e. having uniform star formation history and the same initial chemical composition) stellar populations. Using realistic evolutionary models we predict the abundances and properties of black… ▽ More In large and complicated stellar systems like galaxies it is difficult to predict the number and characteristics of a black hole population. Such populations may be modelled as an aggregation of homogeneous (i.e. having uniform star formation history and the same initial chemical composition) stellar populations. Using realistic evolutionary models we predict the abundances and properties of black holes formed from binaries in these environments. We show that the black hole population will be dominated by single black holes originating from binary disruptions and stellar mergers. Furthermore, we discuss how black hole populations are influenced by such factors as initial parameters, metallicity, initial mass function, and natal kick models. As an example application of our results, we estimate that about 26 microlensing events to happen every year in the direction of the Galactic Bulge due to black holes in a survey like OGLE-IV. Our results may be used to perform in-depth studies related to realistic black hole populations, e.g. observational predictions for space survey missions like Gaia, or Einstein Probe. We prepared a publicly available database with the raw data from our simulations to be used for more in-depth studies. △ Less

Submitted 13 September, 2019; v1 submitted 26 July, 2019; originally announced July 2019.

Comments: 20 pages, 16 figures, revised version

BSE versus StarTrack: implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences

Authors: Sambaran Banerjee, Krzysztof Belczynski, Christopher L. Fryer, Peter Berczik, Jarrod R. Hurley, Rainer Spurzem, Long Wang

Abstract: The masses of stellar-remnant black holes (BH), as a result of their formation via massive single- and binary-stellar evolution, is of high interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star (BNS) mergers. Here we present new developments in the N-body evolution program NBODY7 in regards to its stellar-remnant formation and related schemes. We… ▽ More The masses of stellar-remnant black holes (BH), as a result of their formation via massive single- and binary-stellar evolution, is of high interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star (BNS) mergers. Here we present new developments in the N-body evolution program NBODY7 in regards to its stellar-remnant formation and related schemes. We demonstrate that the newly-implemented stellar-wind and remnant-formation schemes in the NBODY7 code's BSE sector, such as the 'rapid' and the 'delayed' supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produces neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main sequence (ZAMS) mass and metallicity, with those from the StarTrack population-synthesis program. We also demonstrate the new implementations of various natal-kick mechanisms on NSs and BHs such as the 'convection-asymmetry-driven', 'collapse-asymmetry-driven', and 'neutrino-emission-driven' kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. We find that the SN material fallback causes the convection-asymmetry kick to effectively retain similar number and mass of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry kick would cause nearly all BHs to retain in clusters irrespective of remnant formation model and metallicity, whereas the inference of a large number of BHs in GCs would potentially rule out the neutrino-driven kick mechanism. Pre-SN mergers of massive primordial binaries would cause BH masses to deviate from the single-star ZAMS mass-remnant mass relation. Such mergers, at low metallicities, can produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster. △ Less

Submitted 9 April, 2020; v1 submitted 20 February, 2019; originally announced February 2019.

Comments: 27 pages, 15 figures, 2 tables. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 639, A41 (2020)

Nucleosynthetic Yields from Neutron Stars Accreting in Binary Common Envelopes

Authors: James D. Keegans, Chris L. Fryer, Samuel W. Jones, Benoit Cote, Krzysztof Belczynski, Falk Herwig, Marco Pignatari, Alison M. Laird, Christian Aa. Diget

Abstract: Massive-star binaries can undergo a phase where one of the two stars expands during its advanced evolutionary stage as a giant and envelops its companion, ejecting the hydrogen envelope and tightening its orbit. Such a common envelope phase is required to tighten the binary orbit in the formation of many of the observed X-ray binaries and merging compact binary systems. In the formation scenario f… ▽ More Massive-star binaries can undergo a phase where one of the two stars expands during its advanced evolutionary stage as a giant and envelops its companion, ejecting the hydrogen envelope and tightening its orbit. Such a common envelope phase is required to tighten the binary orbit in the formation of many of the observed X-ray binaries and merging compact binary systems. In the formation scenario for neutron star binaries, the system might pass through a phase where a neutron star spirals into the envelope of its giant star companion. These phases lead to mass accretion onto the neutron star. Accretion onto these common-envelope-phase neutron stars can eject matter that has undergone burning near to the neutron star surface. This paper presents nucleosynthetic yields of this ejected matter, using population synthesis models to study the importance of these nucleosynthetic yields in a galactic chemical evolution context. Depending on the extreme conditions in temperature and density found in the accreted material, both proton-rich and neutron-rich nucleosynthesis can be obtained, with efficient production of neutron rich isotopes of low Z material at the most extreme conditions, and proton rich isotopes, again at low Z, in lower density models. Final yields are found to be extremely sensitive to the physical modeling of the accretion phase. We show that neutron stars accreting in binary common envelopes might be a new relevant site for galactic chemical evolution, and therefore more comprehensive studies are needed to better constrain nucleosynthesis in these objects. △ Less

Submitted 5 February, 2019; originally announced February 2019.

Comments: 21 pages, 14 figures

Report number: LA-UR-18-29776

Multiband gravitational-wave event rates and stellar physics

Authors: Davide Gerosa, Sizheng Ma, Kaze W. K. Wong, Emanuele Berti, Richard O'Shaughnessy, Yanbei Chen, Krzysztof Belczynski

Abstract: Joint gravitational-wave detections of stellar-mass black-hole binaries by ground- and space-based observatories will provide unprecedented opportunities for fundamental physics and astronomy. We present a semianalytic method to estimate multiband event rates by combining selection effects of ground-based interferometers (like LIGO/Virgo) and space missions (like LISA). We forecast the expected nu… ▽ More Joint gravitational-wave detections of stellar-mass black-hole binaries by ground- and space-based observatories will provide unprecedented opportunities for fundamental physics and astronomy. We present a semianalytic method to estimate multiband event rates by combining selection effects of ground-based interferometers (like LIGO/Virgo) and space missions (like LISA). We forecast the expected number of multiband detections first by using information from current LIGO/Virgo data, and then through population synthesis simulations of binary stars. We estimate that few to tens of LISA detections can be used to predict mergers detectable on the ground. Conversely, hundreds of events could potentially be extracted from the LISA data stream using prior information from ground detections. In general, the merger signal of binaries observable by LISA is strong enough to be unambiguously identified by both current and future ground-based detectors. Therefore third-generation detectors will not increase the number of multiband detections compared to LIGO/Virgo. We use population synthesis simulations of isolated binary stars to explore some of the stellar physics that could be constrained with multiband events, and we show that specific formation pathways might be overrepresented in multiband events compared to ground-only detections. △ Less

Submitted 9 May, 2019; v1 submitted 31 January, 2019; originally announced February 2019.

Comments: 17 pages, 11 figures. Database and python code available at https://github.com/dgerosa/spops - Published in PRD

Journal ref: Phys. Rev. D 99, 103004 (2019)

Binary neutron star formation and the origin of GW170817

Authors: K. Belczynski, T. Bulik, A. Olejak, M. Chruslinska, N. Singh, N. Pol, L. Zdunik, R. O'Shaughnessy, M. McLaughlin, D. Lorimer, O. Korobkin, E. P. J. van den Heuvel, M. B. Davies, D. E. Holz

Abstract: The first neutron star-neutron star merger was detected in by LIGO/Virgo in a galaxy in which the majority of star formation was taking place a long time ago (11 Gyr). LIGO/Virgo estimated that local cosmic NS-NS merger rate is 110-3840 Gpc^-3 yr^-1. Only some extreme evolutionary models can generate NS-NS merger rates in old host galaxies consistent with the LIGO/Virgo estimate. However, we show… ▽ More The first neutron star-neutron star merger was detected in by LIGO/Virgo in a galaxy in which the majority of star formation was taking place a long time ago (11 Gyr). LIGO/Virgo estimated that local cosmic NS-NS merger rate is 110-3840 Gpc^-3 yr^-1. Only some extreme evolutionary models can generate NS-NS merger rates in old host galaxies consistent with the LIGO/Virgo estimate. However, we show that these models generate rates exceeding empirical Galactic NS-NS merger rates based on the large population of Milky Way NS-NS binaries. Typically, current evolutionary models produce NS-NS merger rates that are consistent with the Milky Way empirical rates. However, these models generate local cosmic NS-NS merger rate in old host galaxies that are below the LIGO/Virgo estimate. The reason behind this tension is the predicted delay time distribution between star formation and NS-NS mergers that favors short delays. Evolutionary models produce a generic steep power-law NS-NS delay time distribution. This limits NS-NS merger rates in old host galaxies. However, we show that such distribution is consistent with observations of Galactic NS-NS binaries; 50% of which show very long merger times (much longer than Hubble time). Once model distributions are convolved with continuous prolonged (10 Gyr) star formation in the Galactic disk, then 20-70% (depending on a model) of the predicted NS-NS population has very long current Galactic merger times. Although NS-NS binaries are formed predominantly with short delay times, many of short delay time systems merge and do not make it to the present, while long delay time systems survive and contribute to the current Galactic NS-NS population. This study highlights the tension between the current evolutionary predictions and the observation of the first NS-NS merger in an old host galaxy. △ Less

Submitted 25 December, 2018; originally announced December 2018.

Comments: comments welcome

Remnants and ejecta of thermonuclear electron-capture supernovae: Constraining oxygen-neon deflagrations in high-density white dwarfs

Authors: S. Jones, F. K. Roepke, C. Fryer, A. J. Ruiter, I. R. Seitenzahl, L. R. Nittler, S. T. Ohlmann, R. Reifarth, M. Pignatari, K. Belczynski

Abstract: (Abridged) The explosion mechanism of electron-capture supernovae (ECSNe) remains equivocal. We attempt to constrain the explosion mechanism (neutron-star-forming implosion or thermonuclear explosion) and the frequency of occurrence of ECSNe using nucleosynthesis simulations of the latter scenario, population synthesis, the solar abundance distribution, pre-solar meteoritic oxide grain isotopic ra… ▽ More (Abridged) The explosion mechanism of electron-capture supernovae (ECSNe) remains equivocal. We attempt to constrain the explosion mechanism (neutron-star-forming implosion or thermonuclear explosion) and the frequency of occurrence of ECSNe using nucleosynthesis simulations of the latter scenario, population synthesis, the solar abundance distribution, pre-solar meteoritic oxide grain isotopic ratio measurements and the white dwarf mass-radius relation. Tracer particles from 3d hydrodynamic simulations were post-processed with a large nuclear reaction network in order to determine the complete compositional state of the bound ONeFe remnant and the ejecta, and population synthesis simulations were performed in order to estimate the ECSN rate with respect to the CCSN rate. The 3d deflagration simulations drastically overproduce the neutron-rich isotopes $^{48}$Ca, $^{50}$Ti, $^{54}$Cr, $^{60}$Fe and several of the Zn isotopes relative to their solar abundances. Using the solar abundance distribution as our constraint, we place an upper limit on the frequency of thermonuclear ECSNe as 1$-$3~\% the frequency at which core-collapse supernovae (FeCCSNe) occur. This is on par with or 1~dex lower than the estimates for ECSNe from single stars. The upper limit from the yields is also in relatively good agreement with the predictions from our population synthesis simulations. The $^{54}$Cr/$^{52}$Cr and $^{50}$Ti/$^{48}$Ti isotopic ratios in the ejecta are a near-perfect match with recent measurements of extreme pre-solar meteoritc oxide grains, and $^{53}$Cr/$^{52}$Cr can also be matched if the ejecta condenses before mixing with the interstellar medium. Theoretical mass-radius relations for the bound ONeFe WD remnants of these explosions are apparently consistent with several observational WD candidates. △ Less

Submitted 4 January, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: 23 pages, 13 figures, 7 tables; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 622, A74 (2019)

Wind Roche Lobe Overflow as a way to make type Ia supernova from the widest symbiotic systems

Authors: Krystian Ilkiewicz, Joanna Mikolajewska, Krzysztof Belczynski, Grzegorz Wiktorowicz, Paulina Karczmarek

Abstract: Symbiotic stars are interacting binaries with one of the longest orbital periods. Since they can contain a massive white dwarf with a high accretion rate they are considered a possible type Ia supernovae (SNe Ia) progenitors. Among symbiotic binaries there are systems containing a Mira donor, which can have orbital periods of a few tens of years and more. This subclass of symbiotic stars due to th… ▽ More Symbiotic stars are interacting binaries with one of the longest orbital periods. Since they can contain a massive white dwarf with a high accretion rate they are considered a possible type Ia supernovae (SNe Ia) progenitors. Among symbiotic binaries there are systems containing a Mira donor, which can have orbital periods of a few tens of years and more. This subclass of symbiotic stars due to their very large separation usually was not considered promising SNe Ia progenitors. We analysed evolution of one of the well studied symbiotic star with a Mira donor, V407 Cyg. We showed that the standard evolution model predicts that the system will not become a SN Ia. However, by simply adding a Wind Roche Lobe Overflow as one of the mass transfer modes we predict that the white dwarf in V407 Cyg will reach the Chandrasekhar limit in 40-200 Myr. △ Less

Submitted 6 December, 2018; originally announced December 2018.

Comments: 4 pages, 3 figures, submitted to MNRAS, comments welcome

The observed vs total population of ULXs

Authors: Grzegorz Wiktorowicz, Jean-Pierre Lasota, Matthew Middleton, Krzysztof Belczynski

Abstract: We have analyzed how anisotropic emission of radiation affects the observed sample of ultraluminous X-ray sources (ULXs) by performing simulations of the evolution of stellar populations, employing recent developments in stellar and binary physics, and by utilizing a geometrical beaming model motivated by theory and observation. Whilst ULXs harboring black hole accretors (BH ULXs) are typically em… ▽ More We have analyzed how anisotropic emission of radiation affects the observed sample of ultraluminous X-ray sources (ULXs) by performing simulations of the evolution of stellar populations, employing recent developments in stellar and binary physics, and by utilizing a geometrical beaming model motivated by theory and observation. Whilst ULXs harboring black hole accretors (BH ULXs) are typically emitting isotropically, the majority of ULXs with neutron star accretors (NS ULXs) are found to be beamed. These findings confirm previous assertions that a significant fraction of ULXs are hidden from view due to a substantial misalignment of the emission beam and the line-of-sight. We find the total number of NS ULXs in regions with constant star formation, solar metallicity, and ages above ~1 Gyr to be higher than the BH ULXs, although observationally both populations are comparable. For lower metallicities BH ULX dominate both the total and observed ULX populations. As far as burst star-formation is concerned, young ULX populations are dominated by BH ULXs, but this changes as the population ages and, post star-formation, NS ULXs dominate both the observed and total population of ULXs. We also compare our simulation output to a previous analytical prediction for the relative ratio of BH to NS ULXs in idealized flux-limited observations and find broad agreement for all but the lowest metallicities. In so doing we find that in such surveys the observed ULX population should be heavily dominated by black-hole systems rather than by systems containing neutron stars. △ Less

Submitted 1 May, 2019; v1 submitted 21 November, 2018; originally announced November 2018.

Comments: 11 pages, 8 figures, minor corrections, accepted to ApJ

Journal ref: The Astrophysical Journal, Volume 875, Issue 1, article id. 53, 12 pp. (2019)

The influence of the distribution of cosmic star formation at different metallicities on the properties of merging double compact objects

Authors: Martyna Chruslinska, Gijs Nelemans, Krzysztof Belczynski

Abstract: Binaries that merge within the local Universe originate from progenitor systems that formed at different times and in various environments. The efficiency of formation of double compact objects is highly sensitive to metallicity of the star formation. Therefore, to confront the theoretical predictions with observational limits resulting from gravitational waves observations one has to account for… ▽ More Binaries that merge within the local Universe originate from progenitor systems that formed at different times and in various environments. The efficiency of formation of double compact objects is highly sensitive to metallicity of the star formation. Therefore, to confront the theoretical predictions with observational limits resulting from gravitational waves observations one has to account for the formation and evolution of progenitor stars in the chemically evolving Universe. In particular, this requires knowledge of the distribution of cosmic star formation rate at different metallicities and times, probed by redshift (SFR(Z,z)). We investigate the effect of the assumed SFR(Z,z) on the properties of merging double compact objects, in particular on their merger rate densities. Using a set of binary evolution models from Chruslinska et al. (2018) we demonstrate that the reported tension between the merger rates of different types of double compact objects and current observational limits in some cases can be resolved if a SFR(Z,z) closer to that expected based on observations of local star-forming galaxies is used, without the need for changing the assumptions about the evolution of progenitor stars of different masses. This highlights the importance of finding tighter constraints on SFR(Z,z) and understanding the associated uncertainties. △ Less

Submitted 8 November, 2018; originally announced November 2018.

Comments: 6 pages, 4 figures, resubmitted to MNRAS after minor revision

Neutron Star Mergers Might not be the Only Source of r-Process Elements in the Milky Way

Authors: Benoit Côté, Marius Eichler, Almudena Arcones, Camilla J. Hansen, Paolo Simonetti, Anna Frebel, Chris L. Fryer, Marco Pignatari, Moritz Reichert, Krzysztof Belczynski, Francesca Matteucci

Abstract: Probing the origin of r-process elements in the universe represents a multi-disciplinary challenge. We review the observational evidence that probe the properties of r-process sites, and address them using galactic chemical evolution simulations, binary population synthesis models, and nucleosynthesis calculations. Our motivation is to define which astrophysical sites have significantly contribute… ▽ More Probing the origin of r-process elements in the universe represents a multi-disciplinary challenge. We review the observational evidence that probe the properties of r-process sites, and address them using galactic chemical evolution simulations, binary population synthesis models, and nucleosynthesis calculations. Our motivation is to define which astrophysical sites have significantly contributed to the total mass of r-process elements present in our Galaxy. We found discrepancies with the neutron star (NS-NS) merger scenario. Assuming they are the only site, the decreasing trend of [Eu/Fe] at [Fe/H]\,$>-1$ in the disk of the Milky Way cannot be reproduced while accounting for the delay-time distribution (DTD) of coalescence times ($\propto~t^{-1}$) derived from short gamma-ray bursts and population synthesis models. Steeper DTD functions ($\propto~t^{-1.5}$) or power laws combined with a strong burst of mergers before the onset of Type~Ia supernovae can reproduce the [Eu/Fe] trend, but this scenario is inconsistent with the similar fraction of short gamma-ray bursts and Type~Ia supernovae occurring in early-type galaxies, and reduces the probability of detecting GW170817 in an early-type galaxy. One solution is to assume an extra production site of Eu that would be active in the early universe, but would fade away with increasing metallicity. If this is correct, this extra site could be responsible for roughly 50% of the Eu production in the early universe, before the onset of Type~Ia supernovae. Rare classes of supernovae could be this additional r-process source, but hydrodynamic simulations still need to ensure the conditions for a robust r-process pattern. △ Less

Submitted 4 March, 2019; v1 submitted 10 September, 2018; originally announced September 2018.

Comments: 29 pages, 10 figures, 1 table, re-submitted to ApJ (revised version)

Impact of inter-correlated initial binary parameters on double black hole and neutron star mergers

Authors: J. Klencki, M. Moe, W. Gladysz, M. Chruslinska, D. E. Holz, K. Belczynski

Abstract: 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. Here, we implement empirically derived inter-correlated distributions of initial binary… ▽ More 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. Here, 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 favors 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 smaller values with respect to the previously used flat q distribution, which is the dominant effect decreasing the rates. New orbital periods only negligibly increase the number of progenitors. 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Å (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 the merger rates. In contrast, we suggest that the uncertainty in cosmic SFR at low metallicity can be a significant factor at play. △ Less

Submitted 13 November, 2018; v1 submitted 23 August, 2018; originally announced August 2018.

Comments: accepted for publication in A&A

Journal ref: A&A 619, A77 (2018)

Spin orientations of merging black holes formed from the evolution of stellar binaries

Authors: Davide Gerosa, Emanuele Berti, Richard O'Shaughnessy, Krzysztof Belczynski, Michael Kesden, Daniel Wysocki, Wojciech Gladysz

Abstract: We study the expected spin misalignments of merging binary black holes (BHs) formed in isolation by combining state-of-the-art population-synthesis models with efficient post-Newtonian evolutions, thus tracking sources from stellar formation to gravitational-wave detection. We present extensive predictions of the properties of sources detectable by both current and future interferometers. We accou… ▽ More We study the expected spin misalignments of merging binary black holes (BHs) formed in isolation by combining state-of-the-art population-synthesis models with efficient post-Newtonian evolutions, thus tracking sources from stellar formation to gravitational-wave detection. We present extensive predictions of the properties of sources detectable by both current and future interferometers. We account for the fact that detectors are more sensitive to spinning BH binaries with suitable spin orientations and find that this significantly impacts the population of sources detectable by LIGO, while this is not the case for 3rd-generation detectors. We find that three formation pathways, differentiated by the order of core collapse and common-envelope phases, dominate the observed population, and that their relative importance critically depends on the recoils imparted to BHs at birth. Our models suggest that measurements of the "effective-spin" parameter $χ_{\rm eff}$ will allow for powerful constraints. For instance, we find that the role of spin magnitudes and spin directions in $χ_{\rm eff}$ can be largely disentangled, and that the symmetry of the effective-spin distribution is a robust indicator of the binary's formation history. Our predictions for individual spin directions and their precessional morphologies confirm and extend early toy models, while exploring substantially more realistic and broader sets of initial conditions. Our main conclusion is that specific subpopulations of BH binaries will exhibit distinctive precessional dynamics: these classes include (but are not limited to) sources where stellar tidal interactions act on sufficiently short timescales, and massive binaries produced in pulsational pair-instability supernovae. Measurements of BH spin orientations have enormous potential to constrain specific evolutionary processes in the lives of massive binary stars. △ Less

Submitted 20 October, 2018; v1 submitted 7 August, 2018; originally announced August 2018.

Comments: 21 pages, 16 figures. Database and python code available at https://github.com/dgerosa/spops - Published in PRD

Journal ref: Phys. Rev. D 98, 084036 (2018)