🔭 astrophysics

Constraints on the $^{12}$ C $(\alpha, \gamma)^{16}$ O and $^{16}$ O+ $^{16}$ O Reaction Rates from Binary Black Holes Detected via Gravitational Wave Signals

By modeling low-metallicity helium stars with MESA, this study demonstrates that uncertainties in the $^{12}$ C $(\alpha, \gamma)^{16}$ O reaction rate significantly shift the predicted lower edge of the binary black hole mass gap, thereby allowing gravitational-wave observations to constrain the astrophysical S factor of this reaction to a range of 137.6–263.4 keV barn.

Original authors: Wenyu Xin, Xiaokun Hou, Xianfei Zhang, Shaolan Bi, Gang Zhao

Published 2026-03-23

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Wenyu Xin, Xiaokun Hou, Xianfei Zhang, Shaolan Bi, Gang Zhao

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant cosmic kitchen where massive stars are the chefs, cooking up elements in their fiery cores. The paper you're reading is like a team of astrophysicists trying to figure out exactly how these chefs work, specifically by looking at the "menu" of black holes they leave behind.

Here is the story of the paper, broken down into simple concepts and everyday analogies.

1. The Mystery of the "Missing Black Holes"

For a long time, astronomers knew about two types of black holes:

The Small Ones: Born from average-sized stars (like our Sun, but much bigger).
The Giant Ones: Born from the most massive stars in the universe.

But there was a weird gap in the middle. Theory said that stars in a certain weight range (roughly 50 to 130 times the mass of our Sun) shouldn't leave behind black holes at all. Instead, they should explode so violently that they blow themselves apart completely, leaving nothing behind. This is called the "Mass Gap."

However, recently, gravitational wave detectors (like LIGO) started hearing "thuds" from black hole mergers. Some of these black holes were sitting right inside that supposed "gap." It was like finding a cat in a room where you were 100% sure no cats could exist. This suggested our understanding of how stars cook their final meals was slightly off.

2. The Two Secret Ingredients

The authors of this paper decided to investigate the "recipe" inside these massive stars. They focused on two specific nuclear reactions (chemical cooking processes) that happen when stars are dying:

The Carbon-to-Oxygen Switch ( $^{12}\text{C} + \alpha \rightarrow ^{16}\text{O}$ ): Imagine a chef turning a block of carbon into oxygen. If this happens too fast or too slow, it changes how much fuel is left for the final explosion.
The Oxygen Fusion ( $^{16}\text{O} + ^{16}\text{O}$ ): This is like two blocks of oxygen smashing together.

The scientists asked: "If we tweak the speed of these reactions, does it change whether a star explodes completely or leaves a black hole behind?"

3. The Simulation: A Cosmic Test Kitchen

To find the answer, they didn't wait for real stars to die (which takes millions of years). Instead, they built a super-computer simulation using a program called MESA.

Think of MESA as a virtual flight simulator for stars.

They created thousands of "virtual stars" with different weights.
They tweaked the "recipe" by speeding up or slowing down the two reactions mentioned above.
They watched to see: Does this star explode and vanish? Or does it collapse into a black hole? If it becomes a black hole, how heavy is it?

4. The Big Discovery: One Ingredient Matters More

The results were fascinating, like finding out that one spice in a soup changes the whole flavor, while another barely makes a difference.

The Carbon-Oxygen Reaction (The Heavy Hitter): This was the main driver. When they changed the speed of this reaction, the "Mass Gap" moved around wildly.
- Analogy: Imagine the Mass Gap is a fence. Changing this reaction is like moving the fence posts. Depending on how you tweak it, the fence could move from being 100 meters away to being 45 meters away.
- Result: By adjusting this reaction, they could explain why we are seeing black holes in the "forbidden zone."
The Oxygen-Oxygen Reaction (The Minor Player): Changing this reaction was like adding a pinch of salt to a huge pot of soup. It changed the top of the black hole weight range slightly, but it didn't move the bottom edge (the start of the gap) very much.

5. Solving the Puzzle with Bayesian Inference

The authors didn't just guess; they used a statistical method called Bayesian Inference.

The Analogy: Imagine you are trying to guess the exact temperature of a room, but you only have a thermometer that gives you a range (e.g., "It's between 68 and 72 degrees"). You also have a few different reports from other people saying what the temperature should be.
The scientists took the "forbidden zone" limits reported by other astronomers (the "reports") and ran them through their computer model.
They asked: "What specific speed for the Carbon-Oxygen reaction would make our model match the real-world observations?"

6. The Conclusion: A New Recipe Card

By combining their computer models with real gravitational wave data, they were able to narrow down the "recipe."

They concluded that the rate at which Carbon turns into Oxygen inside dying stars must be within a specific range (between roughly 137 and 263 "units" of reaction strength). If it's outside this range, the math doesn't add up with the black holes we are actually seeing.

In a nutshell:
The universe is telling us that the "Mass Gap" isn't a hard wall; it's more like a fuzzy zone. By figuring out exactly how stars cook their carbon and oxygen, this paper helps us understand why some stars leave behind heavy black holes while others vanish completely. It's a crucial step in understanding the life and death of the universe's most massive stars.

1. Problem Statement

Gravitational wave (GW) observations of binary black hole (BBH) mergers, particularly events like GW190521 and GW231123, have revealed black hole (BH) components with masses residing in or near the theoretical "pair-instability supernova (PISN) mass gap" (roughly 50–140 $M_\odot$ ). The existence of these objects challenges standard stellar evolution models.

The Core Issue: The location of the PISN mass gap is highly sensitive to nuclear reaction rates, specifically the $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ reaction (which determines the carbon-to-oxygen ratio in the core) and the $^{16}\text{O}+^{16}\text{O}$ reaction (which triggers oxygen burning).
The Conflict: Recent population analyses suggest the lower edge of the BH mass gap is between $\sim44$ and $68 M_\odot$ . However, standard stellar models with current reaction rate uncertainties often predict a gap starting at higher masses ( $\sim95-118 M_\odot$ ), leaving some observed BHs (e.g., the 103 $M_\odot$ component in GW231123) unexplained without invoking additional physics like rotation or convective overshooting.
Goal: The authors aim to quantify how uncertainties in the $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ and $^{16}\text{O}+^{16}\text{O}$ reaction rates affect the BH mass gap boundaries and use GW-derived mass gap constraints to infer the astrophysical S-factor of the $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ reaction.

2. Methodology

The study employs a multi-step approach combining stellar evolution simulations, surrogate modeling, and Bayesian inference.

Stellar Evolution Models:
- Code: MESA (Modules for Experiments in Stellar Astrophysics) version r12778.
- Progenitors: Low-metallicity ( $Z = 10^{-5}$ ) helium stars with initial masses $M(\text{He}) = 30 - 200 M_\odot$ .
- Physics: The models track evolution from the Zero-Age Horizontal Branch (ZAHB) through core collapse (CC), pulsational pair-instability supernovae (PPISN), and full pair-instability supernovae (PISN).
- Nuclear Networks: Uses the approx21 plus co56.net network. Key reaction rates include:
  - $3\alpha$ : Enhanced rate from Kibédi et al. (2020).
  - $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ : High-resolution table (2000 points) from deBoer et al. (2017), varied by $\pm 3\sigma$ (where $\sigma$ represents the uncertainty interval).
  - $^{16}\text{O}+^{16}\text{O}$ : Scaled by a global factor $f_{16\text{O}}$ (values: 0.1, 1, 10).
- Grid: Over 360 model configurations varying $\sigma_{^{12}\text{C}\alpha}$ (from -3 to +3) and $f_{16\text{O}}$ .
Surrogate Modeling (RBF):
- To avoid computationally expensive re-running of MESA for every parameter combination during inference, the authors constructed a Radial Basis Function (RBF) surrogate model.
- This model maps the input parameters ( $\sigma_{^{12}\text{C}\alpha}$ and $\sigma_{^{16}\text{O}} \equiv \log f_{16\text{O}}$ ) to the predicted lower edge of the BH mass gap ( $m$ ).
- A thin-plate spline kernel with a non-zero smoothing parameter was used to prevent overfitting the sparse MESA grid.
Bayesian Inference:
- A Bayesian framework was used to constrain the reaction rates.
- Prior: Bivariate normal distribution centered at 0 (standard rates).
- Likelihood: Defined by the discrepancy between the RBF-predicted mass gap lower edge ( $m(\theta)$ ) and literature values ( $m_0$ ) derived from GW population analyses (e.g., Antonini et al. 2025, Tong et al. 2025).
- Posterior: Sampled to derive constraints on $\sigma_{^{12}\text{C}\alpha}$ and $\sigma_{^{16}\text{O}}$ , which are then converted to the astrophysical S-factor at 300 keV ( $S_{300}$ ).

3. Key Contributions

Systematic Grid Expansion: The study provides a finer grid of simulations compared to previous works (e.g., Farmer et al. 2020), explicitly including the $^{16}\text{O}+^{16}\text{O}$ reaction rate variations alongside the $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ rate.
Differentiated Impact Analysis: The paper rigorously demonstrates that the two reaction rates have distinct impacts on the BH mass gap:
- $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ dominates the lower edge.
- $^{16}\text{O}+^{16}\text{O}$ primarily affects the upper edge.
Quantitative Constraints: It translates observational GW data into specific constraints on the nuclear physics of the $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ reaction, providing a range for $S_{300}$ that complements direct laboratory measurements.

4. Key Results

Impact of $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ :
- Varying the reaction rate from $-3\sigma$ to $+3\sigma$ shifts the predicted BH mass gap significantly.
- Lower Edge: Shifts from $\sim104 M_\odot$ (at $-3\sigma$ ) down to $\sim45 M_\odot$ (at $+3\sigma$ ).
- Upper Edge: Shifts from $\sim184 M_\odot$ to $\sim135 M_\odot$ .
- Mechanism: A higher rate converts more $^{12}\text{C}$ to $^{16}\text{O}$ , reducing the fuel for shell carbon burning, leading to earlier/explosive oxygen ignition and a lower mass threshold for PISN.
Impact of $^{16}\text{O}+^{16}\text{O}$ :
- Scaling $f_{16\text{O}}$ by factors of 0.1, 1, and 10 has a negligible effect on the lower edge of the mass gap (change $< 5 M_\odot$ ).
- It significantly shifts the upper edge (change $> 10 M_\odot$ ), with increasing $f_{16\text{O}}$ raising the upper limit.
Constraints on $S_{300}$ :
- Using the literature estimates for the BH mass gap lower edge ( $44 - 68 M_\odot$ ), the authors constrain the astrophysical S-factor at 300 keV ( $S_{300}$ ) for the $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ reaction.
- Resulting Range: $S_{300} \simeq 137.6 - 263.4$ keV barn.
- The specific value depends on which GW population analysis is adopted (e.g., Antonini et al. vs. Wang et al.), reflecting the current uncertainty in the observational mass gap boundary.
Comparison with Previous Studies:
- The default model in this work predicts a maximum BH mass of $\sim67.65 M_\odot$ (for $Z=10^{-5}$ ), which is higher than some previous studies (e.g., Farmer et al. 2020: $\sim47 M_\odot$ ).
- This difference is attributed to the use of an enhanced $3\alpha$ reaction rate and a high-resolution $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ table, which increase the core carbon abundance ( $X(^{12}\text{C})$ ), promoting stable shell burning and delaying the onset of PISN.

5. Significance

Bridging Nuclear Physics and Astrophysics: The paper successfully demonstrates how GW observations of black hole populations can serve as a probe for nuclear reaction rates that are difficult to measure directly in laboratories at stellar energies.
Refining Stellar Evolution: By isolating the effects of specific reaction rates, the study clarifies the physical mechanisms governing the PISN mass gap, highlighting the critical role of the $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ rate in determining the lower mass limit of the gap.
Implications for GW Events: The findings suggest that the existence of BHs in the "gap" (like the 103 $M_\odot$ component in GW231123) can be reconciled with standard stellar evolution if the $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ rate is on the higher end of current uncertainties, potentially reducing the need for exotic formation channels (like hierarchical mergers) for all such events, though the authors note that some components may still require additional physics like rotation.
Future Directions: The study identifies that while $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ is the dominant factor for the lower gap edge, future work must incorporate rotation and other reaction rates (like $3\alpha$ ) to fully resolve the remaining discrepancies in the BH mass spectrum.