Black Hole Mass Gap Research Articles

Metal-enriched Pair-instability Supernovae: Effects of Rotation

In this paper, we revisit metal-enriched rotating pair-instability supernovae (PISNe) models for metallicities consistent with the Small Magellanic Cloud, the Large Magellanic Cloud (LMC), and 0.1Z ⊙. By calculating multiple models, we intend to clarify mass ranges and the ejected 56Ni masses from the PISNe, and mass-loss histories for progenitors. We find that the choice of the Wolf–Rayet (WR) mass-loss rates is important, and we adopt the recently proposed rate of Sander & Vink, which covers the mass ranges for PISNe progenitors. We show that slow rotation lowers the PISN range because the core mass increases by rotational mixing. On the other hand, when we assume a typical rotation speed for observed OB stars, the mass-loss increase becomes more significant, and the final stellar masses are lower than in nonrotating models. As a result, the typical mass range for bright supernovae (SNe), with a 56Ni mass higher than 10M ⊙ for these fast-rotating models is more than 400 and 350M ⊙ for LMC and 0.1Z ⊙ metallicities, respectively. It is interesting that unlike in previous works, we find oxygen-rich progenitors for most cases. This O-rich progenitor may be consistent with the recently identified PISN candidate SN2018ibb. He-rich progenitors are seen only for relatively dim and metal-poor (Z ≲ 0.1Z ⊙) PISNe. We also discuss the black hole mass gap for metal-enriched PISNe, and we show that the upper bound for the gap is lower than in the Population III case.

Detectability of Single Spinless Stellar-Mass Black Holes through Gravitational Lensing of Gravitational Waves with Advanced LIGO.

We investigate the detectability of gravitational waves that have been lensed by a spinless stellar-mass black hole, with respect to the advanced LIGO. By solving the full relativistic linear wave equations in the spacetime of a Schwarzschild black hole, we find that the strong gravity can create unique signals in the lensed waveform, particularly during the merger and ring-down stages. The differences in terms of fitting factor between the lensed waveform and best-fitted unlensed general relativity template with spin precession and higher-order multipoles are greater than 5% for the lens black hole mass within 70M_{⊙}<M_{lens}<133.33M_{⊙} under advanced LIGO's sensitivity. This is up to 5 times more detectable than the previous analysis based on the weak field approximation for a point mass and covers most part of the black hole mass gap predicted by stellar evolution theory. Based on Bayesian inference, the lensing feature can be distinguished with a signal-to-noise ratio of 12.5 for M_{lens}=70M_{⊙} and 19.2 for M_{lens}=250M_{⊙}, which is attainable for advanced LIGO.

Dark black holes in the mass gap

In the standard picture of stellar evolution, pair-instability — the energy loss in stellar cores due to electron-positron pair production — is predicted to prevent the collapse of massive stars into black holes with mass in the range between approximately 50 and 130 solar masses — a range known as the “black hole mass gap”. LIGO and Virgo detection of black hole binary mergers containing one or both black holes with masses in this mass gap thus challenges the standard picture, possibly pointing to an unexpected merger history, unanticipated or poorly understood astrophysical mechanisms, or new physics. Here, we entertain the possibility that a “dark sector” exists, consisting of dark electrons, dark protons, and electromagnetic-like interactions, but no nuclear forces. Dark stars would inevitably form given such dark sector constituents, possibly collapsing into black holes with masses within the mass gap. We study in detail the cooling processes necessary for successful stellar collapse in the dark sector and show that for suitable choices of the particle masses, we indeed predict populating the mass gap with dark sector black holes. In particular, we numerically find that the heavier of the two dark sector massive particles cannot be lighter than, approximately, the visible sector proton for the resulting dark sector black holes to have masses within the mass gap. We discuss constraints on this scenario and how to test it with future, larger black hole merger statistics.

A Study of Primordial Very Massive Star Evolution. II. Stellar Rotation and Gamma-Ray Burst Progenitors

We calculate new evolutionary models of rotating primordial very massive stars, with initial mass from 100 M ⊙ to 200 M ⊙, for two values of the initial metallicity Z = 0 and Z = 0.0002. For the first time in this mass range, we consider stellar rotation and pulsation-driven mass loss, along with radiative winds. The models evolve from the zero-age main sequence until the onset of pair-instability. We discuss the main properties of the models during their evolution and then focus on the final fate and the possible progenitors of jet-driven events. All tracks that undergo pulsational-pair instability produce successful gamma-ray bursts (GRB) in the collapsar framework, while those that collapse directly to black holes (BH) produce jet-driven supernova events. In these latter cases, the expected black hole mass changes due to the jet propagation inside the progenitor, resulting in different models that should produce BH within the pair-instability black hole mass gap. Successful GRBs predicted here from zero metallicity, and very metal-poor progenitors, may be bright enough to be detected even up to redshift ∼20 using current telescopes such as the Swift-BAT X-ray detector and the JWST.

The Neutron Star to Black Hole Mass Gap in the Frame of the Jittering Jets Explosion Mechanism (JJEM)

I build a toy model in the frame of the jittering jets explosion mechanism (JJEM) of core collapse supernovae that incorporates both the stochastically varying angular momentum component of the material that the newly born neutron star (NS) accretes and the constant angular momentum component, and show that the JJEM can account for the ≃2.5–5M ⊙ mass gap between NSs and black holes (BHs). The random component of the angular momentum results from pre-collapse core convection fluctuations that are amplified by post-collapse instabilities. The fixed angular momentum component results from pre-collapse core rotation. For slowly rotating pre-collapse cores the stochastic angular momentum fluctuations form intermittent accretion disks (or belts) around the NS with varying angular momentum axes in all directions. The intermittent accretion disk/belt launches jets in all directions that expel the core material in all directions early on, hence leaving an NS remnant. Rapidly rotating pre-collapse cores form an accretion disk with angular momentum axis that is about the same as the pre-collapse core rotation. The NS launches jets along this axis and hence the jets avoid the equatorial plane region. Inflowing core material continues to feed the central object from the equatorial plane increasing the NS mass to form a BH. The narrow transition from slow to rapid pre-collapse core rotation, i.e., from an efficient to inefficient jet feedback mechanism, accounts for the sparsely populated mass gap.

New Determination of the 12C(α, γ)16O Reaction Rate and Its Impact on the Black-hole Mass Gap

We present a precise measurement of the asymptotic normalization coefficient (ANC) for the 16O ground state (GS) through the 12C(11B, 7Li)16O transfer reaction using the Quadrupole‐3‐Dipole (Q3D) magnetic spectrograph. The present work sheds light on the existing discrepancy of more than 2 orders of magnitude between the previously reported GS ANC values. This ANC is believed to have a strong effect on the 12C(α, γ)16O reaction rate by constraining the external capture to the 16O ground state, which can interfere with the high-energy tail of the 2+ subthreshold state. Based on the new ANC, we determine the astrophysical S-factor and the stellar rate of the 12C(α, γ)16O reaction. An increase of up to 21% in the total reaction rate is found within the temperature range of astrophysical relevance compared with the previous recommendation of a recent review. Finally, we evaluate the impact of our new rate on the pair-instability mass gap for black holes (BH) by evolving massive helium core stars using the MESA stellar evolution code. The updated 12C(α, γ)16O reaction rate decreases the lower and upper edges of the BH gap about 12% and 5%, respectively.

No Peaks without Valleys: The Stable Mass Transfer Channel for Gravitational-wave Sources in Light of the Neutron Star–Black Hole Mass Gap

Gravitational-wave (GW) detections are starting to reveal features in the mass distribution of double compact objects. The lower end of the black hole (BH) mass distribution is especially interesting as few formation channels contribute here and because it is more robust against variations in the cosmic star formation than the high-mass end. In this work we explore the stable mass transfer channel for the formation of GW sources with a focus on the low-mass end of the mass distribution. We conduct an extensive exploration of the uncertain physical processes that impact this channel. We note that, for fiducial assumptions, this channel reproduces the peak at ∼9 M ☉ in the GW-observed binary BH mass distribution remarkably well and predicts a cutoff mass that coincides with the upper edge of the purported neutron star–black hole (NS–BH) mass gap. The peak and cutoff mass are a consequence of the unique properties of this channel; namely (1) the requirement of stability during the mass transfer phases, and (2) the complex way in which the final compact object masses scale with the initial mass. We provide an analytical expression for the cutoff in the primary component mass and show that this adequately matches our numerical results. Our results imply that selection effects resulting from the formation channel alone can provide an explanation for the purported NS–BH mass gap in GW detections. This provides an alternative to the commonly adopted view that the gap emerges during BH formation.

Repeated Mergers of Black Hole Binaries: Implications for GW190521

The gravitational wave event GW190521 involves the merger of two black holes of ∼85 M ⊙ and ∼66 M ⊙ forming an intermediate-mass black hole (IMBH) of mass ∼142 M ⊙. Both progenitors are challenging to explain within standard stellar evolution as they are within the upper black hole mass gap. We propose a dynamical formation pathway for this IMBH based on multiple mergers in the core of a globular cluster. We identify such scenarios from analysis of a set of 58 N-body simulations using NBODY6-gpu. In one of our simulations, we observe a stellar black hole undergoing a chain of seven binary mergers within 6 Gyr, attaining a final mass of 97.8 M ⊙. We discuss the dynamical interactions that lead to the final IMBH product, as well as the evolution of the black hole population in that simulation. We explore statistically the effects of gravitational recoil on the viability of such hierarchical mergers. From the analysis of all 58 simulations we observe additional smaller chains, tentatively inferring that an IMBH formation through hierarchical mergers is expected in the lifetime of a median-mass globular cluster with probability 0.01 ≲ p ≲ 0.1 without gravitational merger recoil. Using this order-of-magnitude estimate we show that our results are broadly consistent with the rate implied by GW190521, assuming that gravitational recoil ejection of progenitors has a low probability. We discuss implications for future gravitational wave detections, emphasizing the importance of studying such formation pathways for black holes within the upper mass gap as a means to constrain such modeling.

Probing Extremal Gravitational-wave Events with Coarse-grained Likelihoods

Abstract As catalogs of gravitational-wave transients grow, new records are set for the most extreme systems observed to date. The most massive observed black holes probe the physics of pair-instability supernovae while providing clues about the environments in which binary black hole systems are assembled. The least massive black holes, meanwhile, allow us to investigate the purported neutron star–black hole mass gap, and binaries with unusually asymmetric mass ratios or large spins inform our understanding of binary and stellar evolution. Existing outlier tests generally implement leave-one-out analyses, but these do not account for the fact that the event being left out was by definition an extreme member of the population. This results in a bias in the evaluation of outliers. We correct for this bias by introducing a coarse-graining framework to investigate whether these extremal events are true outliers or whether they are consistent with the rest of the observed population. Our method enables us to study extremal events while testing for population model misspecification. We show that this ameliorates biases present in the leave-one-out analyses commonly used within the gravitational-wave community. Applying our method to results from the second LIGO–Virgo transient catalog, we find qualitative agreement with the conclusions of Abbott et al. GW190814 is an outlier because of its small secondary mass. We find that neither GW190412 nor GW190521 is an outlier.

Capability for detection of GW190521-like binary black holes with TianQin

The detection of GW190521 gains huge attention because it is the most massive binary that LIGO and Virgo ever confidently detected until the release of GWTC-3 (GW190426_190642 is more massive), and it is the first black hole merger whose remnant is believed to be an intermediate mass black hole. Furthermore, the primary black hole mass falls in the black hole mass gap, where the pair-instability supernova prevents the formation of astrophysical black holes in this range. In this paper, we systematically explore the prospect of TianQin on detecting GW190521-like sources. For sources with small orbital eccentricities, (i) TianQin could resolve up to a dozen of sources with signal-to-noise ratio (SNR) larger than 8. Even if the signal-to-noise ratio threshold increases to 12, TianQin could still detect GW190521-like binaries. (ii) The parameters of sources merging within several years would be precisely recovered. The precision of coalescence time and sky localization closes to $1\ {\rm s}$ and $1\ {\rm deg^{2}}$ respectively. This indicates that TianQin could provide early warnings for ground-based gravitational waves detectors and electromagnetic telescopes for these sources. Furthermore, TianQin could distinguish the formation channels of these sources by measuring the orbital eccentricities with a relative precision of $10^{-4}$. (iii) TianQin could constrain the Hubble constant with a $10\%$ precision with GW190521-like sources. Finally, for very eccentric GW190521-like sources, although their gravitational wave signal might be too weak for TianQin to detect, even the null detection of TianQin could still present a significant contribution to the understanding of the underlying science.

Filling the black hole mass gap: Avoiding pair instability in massive stars through addition of nonnuclear energy

In standard stellar evolution, stars with masses ranging from approximately $150\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ to $240\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ are expected to evolve to a pair-instability supernova with no black hole (BH) remnant. This evolutionary behavior leads to a predicted gap in the black hole mass function from approximately $50\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ to $140\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$. Yet the LIGO and Virgo collaborations [Phys. Rev. Lett. 125, 101102 (2020)] recently discovered black holes with masses of $66\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ and $85\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ in the gravitational-wave event GW190521. We propose a new method to populate the BH mass gap. If an energy source is added throughout the star in addition to nuclear fusion, it is possible for the altered evolution to avoid the complete destruction of a pair-instability supernova, and instead a BH remnant is left behind. An example of an extra energy source is dark matter annihilation within the star, but our results hold more generally. We show this phenomenon by exploring the effect of adding an energy source independent of temperature and density to a $180\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ star, using the mesa one-dimensional stellar evolution software. $\text{If}\ensuremath{\sim}50%$ of the star's energy is due to this new source, the star is capable of avoiding the pair instability entirely, reaching a presupernova mass of $119\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ before evolving into a BH remnant in the mass gap.

Find the Gap: Black Hole Population Analysis with an Astrophysically Motivated Mass Function

Abstract We introduce a novel black hole mass function that realistically models the physics of pair-instability supernovae with a minimal number of parameters. Applying this to all events in the LIGO-Virgo Gravitational-Wave Transient Catalog 2 (GWTC-2), we detect a peak at M BHMG = 46 − 6 + 17 M ⊙ . Repeating the analysis without the black holes from the event GW190521, we find this feature at M BHMG = 54 ± 6 M ⊙. These results establish the edge of the anticipated “black hole mass gap” at a value compatible with the expectation from standard stellar structure theory. The mass gap manifests itself as a discontinuity in the mass function and is populated by a distinct, less-abundant population of higher-mass black holes. We find that the primary black hole population scales with power-law index −1.95 ± 0.51 (−1.97 ± 0.44) with (without) GW190521, consistent with models of star formation. Using Bayesian techniques, we establish that our mass function fits a new catalog of black hole masses approximately as well as pre-existing phenomenological mass functions. We also remark on the implications of these results for constraining or discovering new phenomena in nuclear and particle physics.

Missing in axion: Where are XENON1T’s big black holes?

We pioneer the black hole mass gap as a powerful new tool for constraining new particles. A new particle that couples to the Standard Model – such as an axion – acts as an additional source of loss in the cores of population-III stars, suppressing mass lost due to winds and quenching the pair-instability. This results in heavier astrophysical black holes. As an example, using stellar simulations we show that the solar axion explanation of the recent XENON1T excess implies astrophysical black holes of ∼ 56 M ⊙ , squarely within the black hole mass gap predicted by the Standard Model.

Probing the Nuclear Equation of State from the Existence of a ∼2.6 M⊙ Neutron Star: The GW190814 Puzzle

On 14 August 2019, the LIGO/Virgo collaboration observed a compact object with mass ∼2.59−0.09+0.08M⊙, as a component of a system where the main companion was a black hole with mass ∼23M⊙. A scientific debate initiated concerning the identification of the low mass component, as it falls into the neutron star–black hole mass gap. The understanding of the nature of GW190814 event will offer rich information concerning open issues, the speed of sound and the possible phase transition into other degrees of freedom. In the present work, we made an effort to probe the nuclear equation of state along with the GW190814 event. Firstly, we examine possible constraints on the nuclear equation of state inferred from the consideration that the low mass companion is a slow or rapidly rotating neutron star. In this case, the role of the upper bounds on the speed of sound is revealed, in connection with the dense nuclear matter properties. Secondly, we systematically study the tidal deformability of a possible high mass candidate existing as an individual star or as a component one in a binary neutron star system. As the tidal deformability and radius are quantities very sensitive on the neutron star equation of state, they are excellent counters on dense matter properties. We conjecture that similar isolated neutron stars or systems may exist in the universe and their possible future observation will shed light on the maximum neutron star mass problem.

New physics and the black hole mass gap

In the Standard Model, the so-called pair instability from creation of electron-positron pairs reduces the mass of low-metallicity population-III stars as they collapse, unless they are heavy enough to make the process inefficient, leading to an expected mass gap in the black hole spectrum. In this work, the authors argue that future observations of the black hole population will allow one to test this prediction, and use it to put constraints on new particles, such as axions, that would cause an additional instability.

Formation of GW190521 from stellar evolution: the impact of the hydrogen-rich envelope, dredge-up, and 12C(α, γ)16O rate on the pair-instability black hole mass gap

ABSTRACT Pair-instability (PI) is expected to open a gap in the mass spectrum of black holes (BHs) between ≈40–65 and ≈120 M⊙. The existence of the mass gap is currently being challenged by the detection of GW190521, with a primary component mass of $85^{+21}_{-14}$ M⊙. Here, we investigate the main uncertainties on the PI mass gap: the 12C(α, γ)16O reaction rate and the H-rich envelope collapse. With the standard 12C(α, γ)16O rate, the lower edge of the mass gap can be 70 M⊙ if we allow for the collapse of the residual H-rich envelope at metallicity Z ≤ 0.0003. Adopting the uncertainties given by the starlib database, for models computed with the 12C(α, γ)16O rate $-1\, \sigma$, we find that the PI mass gap ranges between ≈80 and ≈150 M⊙. Stars with MZAMS &amp;gt; 110 M⊙ may experience a deep dredge-up episode during the core helium-burning phase, that extracts matter from the core enriching the envelope. As a consequence of the He-core mass reduction, a star with MZAMS = 160 M⊙ may avoid the PI and produce a BH of 150 M⊙. In the $-2\, {}\sigma {}$ case, the PI mass gap ranges from 92 to 110 M⊙. Finally, in models computed with 12C(α, γ)16O $-3\, {}\sigma {}$, the mass gap is completely removed by the dredge-up effect. The onset of this dredge-up is particularly sensitive to the assumed model for convection and mixing. The combined effect of H-rich envelope collapse and low 12C(α, γ)16O rate can lead to the formation of BHs with masses consistent with the primary component of GW190521.

Modified gravity and the black hole mass gap

We pioneer the black hole mass gap as a powerful new tool for constraining modified gravity theories. These theories predict fifth forces that alter the structure and evolution of population-III stars, exacerbating the pair-instability. This results in the formation of lighter astrophysical black holes and lowers both the upper and lower edges of the mass gap. These effects are explored using detailed numerical simulations to derive quantitative predictions that can be used as theoretical inputs for Bayesian data analysis. We discuss detection strategies in light of current and upcoming data as well as complications that may arise due to environmental screening. To demonstrate the constraining power of the mass gap, we present a novel test of the strong equivalence principle where we apply our results to an analysis of the first ten LIGO/Virgo binary black hole merger events to obtain a 7% bound on the relative difference between the gravitational constant experienced by baryonic matter, and that experienced by black holes, $\mathrm{\ensuremath{\Delta}}G/G$. The recent GW190521 event resulting from two black holes with masses in the canonical mass gap can be explained by modified gravity if the event originated from an unscreened galaxy where the strength of gravity is either enhanced or reduced by $\ensuremath{\sim}30%$ relative to its strength in the Solar System.

Populating the Upper Black Hole Mass Gap through Stellar Collisions in Young Star Clusters

Abstract Theoretical modeling of massive stars predicts a gap in the black hole (BH) mass function above ∼40–50 M ⊙ for BHs formed through single star evolution, arising from (pulsational) pair-instability supernovae (PISNe). However, in dense star clusters, dynamical channels may exist that allow construction of BHs with masses in excess of those allowed from single star evolution. The detection of BHs in this so-called “upper-mass gap” would provide strong evidence for the dynamical processing of BHs prior to their eventual merger. Here, we explore in detail the formation of BHs with masses within or above the pair-instability gap through collisions of young massive stars in dense star clusters. We run a suite of 68 independent cluster simulations, exploring a variety of physical assumptions pertaining to growth through stellar collisions, including primordial cluster mass segregation and the efficiency of envelope stripping during collisions. We find that as many as ∼20% of all BH progenitors undergo one or more collisions prior to stellar collapse and up to ∼1% of all BHs reside within or above the pair-instability gap through the effects of these collisions. We show that these BHs readily go on to merge with other BHs in the cluster, creating a population of massive BH mergers at a rate that may compete with the “multiple-generation” merger channel described in other analyses. This has clear relevance for the formation of very massive BH binaries as recently detected by the Laser Interferometer Gravitational-Wave Observatory/Virgo in GW190521. Finally, we describe how stellar collisions in clusters may provide a unique pathway to PISNe and briefly discuss the expected rate of these events and other electromagnetic transients.

The impact of stellar rotation on the black hole mass-gap from pair-instability supernovae

Models of pair-instability supernovae (PISNe) predict a gap in black hole (BH) masses between ∼45 M⊙ and 120 M⊙, which is referred to as the upper BH mass-gap. With the advent of gravitational-wave astrophysics, it has become possible to test this prediction, and there is an important associated effort to understand which theoretical uncertainties modify the boundaries of this gap. In this work we study the impact of rotation on the hydrodynamics of PISNe, which leave no compact remnant, as well as the evolution of pulsational-PISNe (PPISNe), which undergo thermonuclear eruptions before forming a compact object. We perform simulations of nonrotating and rapidly rotating stripped helium stars in a metal-poor environment (Z⊙/50) in order to resolve the lower edge of the upper mass-gap. We find that the outcome of our simulations is dependent on the efficiency of angular momentum transport: models that include efficient coupling through the Spruit-Tayler dynamo shift the lower edge of the mass-gap upward by ∼4%, while simulations that do not include this effect shift it upward by ∼15%. From this, we expect that the lower edge of the upper mass-gap is dependent on BH spin, which can be tested as the number of observed BH mergers increases. Moreover, we show that stars undergoing PPISNe have extended envelopes (R ∼ 10 − 1000 R⊙) at iron-core collapse, making them promising progenitors for ultra-long gamma-ray bursts.

Sensitivity of the lower edge of the pair-instability black hole mass gap to the treatment of time-dependent convection

ABSTRACT Gravitational-wave detections are now probing the black hole (BH) mass distribution, including the predicted pair-instability mass gap. These data require robust quantitative predictions, which are challenging to obtain. The most massive BH progenitors experience episodic mass ejections on time-scales shorter than the convective turnover time-scale. This invalidates the steady-state assumption on which the classic mixing length theory relies. We compare the final BH masses computed with two different versions of the stellar evolutionary code $\tt{MESA}$: (i) using the default implementation of Paxton et al. (2018) and (ii) solving an additional equation accounting for the time-scale for convective deceleration. In the second grid, where stronger convection develops during the pulses and carries part of the energy, we find weaker pulses. This leads to lower amounts of mass being ejected and thus higher final BH masses of up to ∼$5\, \mathrm{M}_\odot$. The differences are much smaller for the progenitors that determine the maximum mass of BHs below the gap. This prediction is robust at $M_{\rm BH, max}\simeq 48\, \mathrm{M}_\odot$, at least within the idealized context of this study. This is an encouraging indication that current models are robust enough for comparison with the present-day gravitational-wave detections. However, the large differences between individual models emphasize the importance of improving the treatment of convection in stellar models, especially in the light of the data anticipated from the third generation of gravitational-wave detectors.