Temperature-resolved sensitivities of $^{56}{\rm Ni}$ production to helium-burning reactions in pair-instability supernovae

This paper introduces a temperature-resolved Monte Carlo approach to demonstrate that pair-instability supernova nucleosynthesis, specifically $^{56}$Ni production, is most sensitive to variations in the triple-$\alpha$ and $^{12}$C($\alpha$,$\gamma$)$^{16}$O reaction rates at approximately $2.5 \times 10^8$ K, where these rates exert opposite influences on the pre-carbon-burning C/O composition.

The Gist

Imagine a massive star, hundreds of times heavier than our Sun, as a giant cosmic pressure cooker. As it burns through its fuel, it eventually faces a dramatic explosion called a Pair-Instability Supernova (PISN). When this happens, the star doesn't just fade away; it completely shatters, creating a brilliant flash of light powered by a massive amount of radioactive iron (specifically an isotope called Nickel-56).

Astronomers want to know exactly how bright these explosions will be, because that brightness tells us how much Nickel-56 was made. However, predicting this brightness is like trying to guess the outcome of a complex recipe when you aren't sure about the exact measurements of the ingredients.

The Problem: Uncertain Ingredients

In the life of a massive star, two specific nuclear reactions act as the main chefs during the "helium-burning" phase:

1. The Triple-Alpha Reaction: This is the "builder." It takes three helium particles and smashes them together to create Carbon.
2. The Carbon-Alpha Reaction: This is the "converter." It takes that newly made Carbon and turns it into Oxygen.

For decades, scientists have been unsure about the exact "speed" or "efficiency" of these two reactions. It's like knowing you need to bake a cake, but you don't know if your oven is set to 350°F or 400°F, or if your measuring cups are slightly off. Because these reactions compete with each other (one makes Carbon, the other eats it), even tiny uncertainties in their rates can change the final mix of Carbon and Oxygen inside the star. And that mix determines how violent the final explosion will be.

The Old Way vs. The New Way

Previously, scientists tried to solve this by saying, "Let's just assume these reactions could be twice as fast or half as fast everywhere in the star's life." They would run simulations with these extreme, uniform changes to see the best and worst-case scenarios.

But this is like saying, "Maybe my oven is broken at every temperature from 100°F to 500°F." In reality, the uncertainty might only matter at a specific temperature, like when the oven is preheating. The old method couldn't tell you when the uncertainty mattered most.

The New Approach: A Temperature-Specific Detective

The authors of this paper developed a new method, which they call a "Temperature-Resolved Monte Carlo approach."

Think of it like this: Instead of guessing the oven temperature for the whole day, they ran thousands of simulations where they randomly tweaked the reaction speeds at every single temperature step independently.

• At 100 million degrees, they might speed up the Carbon reaction.
• At 200 million degrees, they might slow down the Triple-Alpha reaction.
• At 300 million degrees, they might leave everything alone.

By running 10,000 different versions of the star's life with these random tweaks, they could look at the final result (the amount of Nickel-56) and ask: "Which specific temperature tweak caused the biggest change in the final explosion?"

The Big Discovery: The "Sweet Spot"

The study found a very specific "sweet spot" in the star's life. The reactions mattered most when the star's core was at a temperature of roughly 250 million degrees (2.5 × 10⁸ K).

Here is the interesting part:

• At this specific temperature, making the Carbon-Alpha reaction (the converter) faster led to more Nickel-56 in the explosion.
• Conversely, making the Triple-Alpha reaction (the builder) faster led to less Nickel-56.

Why? Because at this specific temperature, the balance between Carbon and Oxygen is being set. If you convert more Carbon into Oxygen early on, the star stays more compact and explodes more violently later, creating more Nickel. If you keep too much Carbon, it burns up too early, changing the star's structure and resulting in a weaker explosion.

The paper shows that the "recipe" for the star's final explosion is essentially printed onto the Carbon/Oxygen mix at this one specific temperature. If you get the rates right at 250 million degrees, you can predict the explosion's brightness much better.

A Real-World Test: SN 2018ibb

To show how this works, the authors looked at a real supernova candidate called SN 2018ibb. This star was observed to be extremely bright, suggesting it produced a huge amount of Nickel-56 (between 25 and 44 times the mass of our Sun).

When they applied their new method:

• If they assumed the star had a "normal" amount of heavy elements (metallicity), they couldn't reproduce that brightness, even with their best guesses.
• However, when they assumed the star was born in a very "clean" environment (very low metallicity), their model successfully matched the observed brightness.

This suggests that SN 2018ibb likely came from a very metal-poor star, and that the specific reaction rates at that 250-million-degree sweet spot were crucial in creating the massive explosion we saw.

Summary

In short, this paper is like finding the exact moment in a cooking process where a tiny change in the heat makes the difference between a burnt cake and a perfect one. The authors discovered that for massive stars, the "perfect moment" is when the core is at 250 million degrees. By focusing on the reaction rates at this specific temperature, we can finally understand why some of these cosmic explosions are so incredibly bright and use that knowledge to decode the history of the universe.

Technical Summary

Technical Summary: Temperature-Resolved Sensitivities of 56Ni Production in Pair-Instability Supernovae

Problem Statement
The amount of radioactive $^{56}\text{Ni}$ synthesized in Pair-Instability Supernovae (PISNe) is a primary determinant of the supernova light curve brightness and a key observable for identifying these events. However, theoretical predictions of the $^{56}\text{Ni}$ yield suffer from significant uncertainties, particularly regarding nuclear reaction rates during the helium-burning phase. Previous studies have addressed these uncertainties by uniformly scaling reaction rates (e.g., the triple-$\alpha$ and $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ reactions) across the entire temperature range. This approach fails to identify which specific temperature regimes drive the sensitivity of the final yield, as reaction cross-sections and dominant physical processes vary significantly with energy. Furthermore, treating rates as uniformly uncertain ignores the competitive nature of these reactions in establishing the pre-explosive Carbon/Oxygen (C/O) core composition.

Methodology
The authors propose a temperature-resolved Monte Carlo (MC) approach to isolate the specific temperature windows where variations in helium-burning reaction rates most strongly influence $^{56}\text{Ni}$ production.

• Stellar Evolution: The study utilizes the MESA code (version r15140) to simulate the evolution of Very Massive Stars (VMSs), specifically a hydrogen-envelope-stripped model with an initial mass of $100\,M_\odot$ and metallicity $Z=10^{-5}$.
• Reaction Rate Sampling: Instead of applying a single global multiplier to reaction rates, the authors generate thousands of stellar evolution models where the multipliers for the triple-$\alpha$ and $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ reactions are sampled independently at each temperature point from a uniform distribution $[0.5, 2.0]$.
• Statistical Analysis: For each model, the final $^{56}\text{Ni}$ mass ($M_{56\text{Ni}}$) is recorded. The authors calculate the Pearson correlation coefficient ($r$) between the sampled rate multiplier at a specific temperature $T$ and the resulting $M_{56\text{Ni}}$. This is performed for both reactions across the temperature range.
• Compositional Imprint: To explain the correlation results, the authors analyze the "effective produced C/O number ratio" at the onset of carbon burning ($T_c = 10^8.7$ K) as a function of the ratio of the sampled multipliers ($F_i(T)$) at various temperatures.

Key Results

1. Peak Sensitivity Temperature: The correlation analysis reveals that both the triple-$\alpha$ and $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ reactions exhibit their strongest sensitivity to the final $^{56}\text{Ni}$ yield at a specific temperature of $T \simeq 2.5 \times 10^8$ K.
2. Opposite Correlation Signs: At this temperature, the two reactions show opposite correlation signs:
   • $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$: Shows a strong positive correlation ($r \simeq 0.64$). An increase in this rate leads to higher $^{56}\text{Ni}$ yields.
   • Triple-$\alpha$: Shows a negative correlation ($r \simeq -0.55$). An increase in this rate leads to lower $^{56}\text{Ni}$ yields.
3. Physical Mechanism: The sensitivity at $2.5 \times 10^8$ K corresponds to the regime where the ratio of the reaction rates is most clearly imprinted on the pre-carbon-burning C/O composition.
   • A higher $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ rate converts more $^{12}\text{C}$ to $^{16}\text{O}$, resulting in an oxygen-rich core. This leads to a more compact progenitor structure, triggering more violent oxygen burning and higher $^{56}\text{Ni}$ production.
   • Conversely, a higher triple-$\alpha$ rate increases the carbon abundance. This enhances pre-explosive carbon burning, which injects energy earlier, alters the stellar structure, and suppresses subsequent $^{56}\text{Ni}$ production.
4. Convergence: The study demonstrates that candidate temperature windows for dominant sensitivities emerge after approximately $O(10^2)$ successful stellar evolution models, although larger samples are required for robust amplitude estimates.

Significance and Claims
The paper claims that PISN nucleosynthesis can serve as a probe for helium-burning reaction rates specifically in the low-temperature regime around $2.5 \times 10^8$ K. By identifying that the final yield is most sensitive to rate variations in this specific window rather than the entire temperature range, the authors provide a framework for:

• Interpreting Observations: Future PISN observations can be used to constrain reaction rates in this specific temperature regime.
• Guiding Nuclear Physics: The results highlight specific low-energy regimes where improved nuclear cross-section measurements would have the most significant impact on PISN predictions.
• Methodological Advancement: The temperature-resolved MC approach offers a more efficient and physically grounded alternative to uniform scaling, capable of disentangling competing nuclear effects (production vs. destruction channels) that uniform scaling obscures.

The authors apply this framework to the candidate SN 2018ibb, suggesting that its observed $^{56}\text{Ni}$ mass ($\sim 25\text{--}44\,M_\odot$) is consistent with models only if the progenitor had a lower metallicity ($Z \approx 10^{-2} Z_\odot$) than typically expected for the bulk of the PISN population, illustrating the utility of their sensitivity analysis in interpreting specific astrophysical candidates.