Predicting the Heaviest Black Holes below the Pair Instability Gap

Predicting the Heaviest Black Holes below the Pair Instability Gap

Traditionally, the pair instability (PI) mass gap is located between 50,and 130,$M_{odot}$, with stellar mass black holes (BHs) expected to “pile up” towards the lower PI edge. However, this lower PI boundary is based on the assumption that the star has already lost its hydrogen (H) envelope. With the announcement of an “impossibly” heavy BH of 85,$M_{odot}$ as part of GW,190521 located inside the traditional PI gap, we realised that blue supergiant (BSG) progenitors with small cores but large Hydrogen envelopes at low metallicity ($Z$) could directly collapse to heavier BHs than had hitherto been assumed. The question of whether a single star can produce such a heavy BH is important, independent of gravitational wave events. Here, we systematically investigate the masses of stars inside the traditional PI gap by way of a grid of 336 detailed MESA stellar evolution models calculated across a wide parameter space, varying stellar mass, overshooting, rotation, semi-convection, and $Z$. We evolve low $Z$ stars in the range $10^{-3} Traditionally, the pair instability (PI) mass gap is located between 50,and 130,$M_{odot}$, with stellar mass black holes (BHs) expected to “pile up” towards the lower PI edge. However, this lower PI boundary is based on the assumption that the star has already lost its hydrogen (H) envelope. With the announcement of an “impossibly” heavy BH of 85,$M_{odot}$ as part of GW,190521 located inside the traditional PI gap, we realised that blue supergiant (BSG) progenitors with small cores but large Hydrogen envelopes at low metallicity ($Z$) could directly collapse to heavier BHs than had hitherto been assumed. The question of whether a single star can produce such a heavy BH is important, independent of gravitational wave events. Here, we systematically investigate the masses of stars inside the traditional PI gap by way of a grid of 336 detailed MESA stellar evolution models calculated across a wide parameter space, varying stellar mass, overshooting, rotation, semi-convection, and $Z$. We evolve low $Z$ stars in the range $10^{-3}