The impact of new ($\alpha$, n) reaction rates on the weak s-process in metal-poor massive stars

This study demonstrates that updating specific ($\alpha$, n) reaction rates, particularly for $^{17}$O and $^{22}$Ne, significantly enhances the weak s-process yields in metal-poor massive stars by tens of times, with the $^{17}$O($\alpha$, n) reaction playing a dominant role across all evolutionary stages, especially in more massive models.

The Gist

Imagine the universe as a giant, cosmic kitchen. For billions of years, massive stars have been the master chefs, cooking up the heavy elements that make up everything around us—from the iron in your blood to the gold in your jewelry.

This paper is about a specific recipe these chefs use to make a particular batch of ingredients: the "Weak s-process." Think of this as the "lighter" version of heavy-element cooking, responsible for creating elements like Zinc, Gallium, and Germanium (the stuff between Copper and Zirconium on the periodic table).

Here is the story of what the authors discovered, explained simply:

1. The Ingredients: Neutrons are the "Magic Dust"

To cook these heavy elements, stars need neutrons. Think of neutrons as "magic dust" that gets sprinkled onto lighter atoms to make them heavier.

• The Main Chef (22Ne): Usually, the star gets its magic dust from a reaction involving Neon-22.
• The Spoiler (16O): But there's a problem. The star is full of Oxygen-16, which acts like a "neutron vacuum cleaner." It sucks up the magic dust before it can do its job.
• The Rescue (17O): However, if the star can turn some of that vacuum cleaner (Oxygen-16) into Oxygen-17, and then hit it with an alpha particle (a helium nucleus), it can actually release the magic dust back into the kitchen.

2. The Problem: Old Recipes vs. New Science

For a long time, astronomers used an old "cookbook" (called JINA REACLIB) to calculate how fast these chemical reactions happen. But recently, scientists in the lab have updated the measurements for two key reactions:

1. The Neon Reaction: How Neon-22 releases neutrons.
2. The Oxygen Reaction: How Oxygen-17 releases neutrons.

The authors of this paper asked: "What happens to our cosmic cooking if we swap the old cookbook for these new, more accurate recipes?"

3. The Experiment: Simulating Four Massive Stars

They built computer models of four massive stars (15, 20, 25, and 30 times the mass of our Sun) that are metal-poor.

• Why metal-poor? In the early universe, stars had very few heavy elements. In these stars, the "Neon Chef" is weak because there isn't much Neon to begin with. This makes the "Oxygen Rescue" even more critical.

4. The Big Discovery: The Oxygen Update is a Game-Changer

When they ran the simulation with the new numbers, they found two major things:

A. The Neon Update (The "Good" News, but limited)
The new recipe for Neon made the cooking slightly more efficient in the later stages of the star's life (when it's burning Carbon and Neon). It increased the production of heavy elements by about 20 to 25 times. That's a lot, but it's not the whole story.

B. The Oxygen Update (The "Huge" News)
This was the shocker. The new recipe for Oxygen-17 changed everything.

• In the old recipe, the Oxygen-17 reaction was too slow to matter much.
• In the new recipe, it becomes a neutron powerhouse.
• The Result: The production of heavy elements (like Gallium and Germanium) skyrocketed by over 100 times in some cases!

The Analogy: Imagine you are trying to fill a bucket with a leaky hose (the old recipe). You get a little water. Then, someone hands you a fire hose (the new Oxygen recipe). Suddenly, the bucket fills up instantly. The new Oxygen reaction is that fire hose.

5. Why Does This Matter?

• Massive Stars Matter More: The effect was even stronger in the heaviest stars (30 solar masses). This suggests that the biggest stars in the early universe were actually much better at making these specific elements than we thought.
• The "Missing" Elements: Astronomers have been looking at old, metal-poor stars and found they have more of these heavy elements than our old models predicted. This paper says, "Aha! We were using the wrong recipe. If we use the new one, the math finally matches the observations."
• The Solar System: When they averaged the results for all types of stars, they found that these new rates could explain why we have so much Zinc, Gallium, and Germanium in our solar system today.

6. The Catch: We Need Better Measurements

The paper ends with a call to action. The difference between the old and new recipes is massive (sometimes a factor of 100!). This means our current lab measurements of how Oxygen-17 behaves are still a bit shaky.

• The Takeaway: We need physicists to go back to the lab and measure these reactions even more precisely. If we get the numbers right, we can finally understand exactly how the universe cooked up the ingredients for our existence.

In a nutshell: The universe's "heavy element factory" is much more efficient than we thought, thanks to a specific reaction involving Oxygen that we finally got the math right on. It turns out, the "neutron vacuum cleaner" (Oxygen) can actually become a "neutron generator" if you hit it just right!

Technical Summary

1. Problem Statement

The weak s-process (ws-process) in massive stars ($M \gtrsim 12 M_\odot$) is the primary source of neutron-rich isotopes in the mass range $A \approx 60-90$. In metal-poor environments, the abundance of the primary neutron source, $^{22}\text{Ne}$, is limited by the initial metallicity (as it is synthesized from $^{14}\text{N}$). Consequently, the contribution of the $^{22}\text{Ne}(\alpha, n)^{25}\text{Mg}$ reaction is expected to be weak.

However, recent observations indicate that ws-process elements in metal-poor stars are more abundant than standard models predict. A potential explanation lies in the competition between neutron capture and release involving oxygen isotopes. In metal-poor stars, $^{16}\text{O}$ is abundant and acts as a "neutron poison" via $^{16}\text{O}(n, \gamma)^{17}\text{O}$. The fate of the resulting $^{17}\text{O}$ depends on competing reactions:

• $^{17}\text{O}(\alpha, n)^{20}\text{Ne}$ (releases neutrons, aiding the ws-process).
• $^{17}\text{O}(\alpha, \gamma)^{21}\text{Ne}$ (captures neutrons, hindering the process).

Uncertainties in the reaction rates for $^{17}\text{O}+\alpha$ and $^{22}\text{Ne}+\alpha$ significantly affect nucleosynthesis yields. This study investigates how newly proposed reaction rates for these channels alter ws-process yields in non-rotating, metal-poor massive stars compared to standard rates from the JINA REACLIB database.

2. Methodology

Stellar Evolution Models

• Code: Modules for Experiments in Stellar Astrophysics (MESA, v12778).
• Models: Four metal-poor stars with initial metallicity $Z = 10^{-3}$ ($0.1 Z_\odot$) and Zero-Age Main Sequence (ZAMS) masses of 15, 20, 25, and $30 M_\odot$.
• Evolution: Tracked from ZAMS until the onset of iron core collapse (infall velocity $\approx 10^3$ km s$^{-1}$).
• Physics:
  • Nuclear network: mesa 161.net (161 isotopes).
  • Resolution: High spatial resolution ($>3,500$ cells) and strict limits on density/temperature changes to ensure convergence.
  • Mixing: Convection, overshoot, and semiconvection included; global mixing neglected via a minimum diffusion coefficient.
  • Mass Cut: Defined at the location of the steepest pressure/density gradient ($M(V/U)_{max}$), separating the proto-neutron star from the ejecta.

Nucleosynthesis Post-Processing

• Code: WinNet (a one-zone nuclear reaction network code).
• Network: $\sim$2000 isotopes (from neutrons to Thorium, $Z=90$).
• Trajectories: Temperature and density profiles from MESA were mapped to WinNet to calculate detailed isotopic yields up to the onset of core collapse.
• Reaction Rate Scenarios: Four distinct "recipes" were tested for each model (Table 1):
  1. Default: JINA REACLIB (Cyburt et al. 2010) for both $^{17}\text{O}+\alpha$ and $^{22}\text{Ne}+\alpha$.
  2. New $^{22}\text{Ne}$: Wiescher et al. (2023) for $^{22}\text{Ne}+\alpha$; REACLIB for $^{17}\text{O}+\alpha$.
  3. New $^{17}\text{O}$: Best et al. (2013) for $^{17}\text{O}+\alpha$; REACLIB for $^{22}\text{Ne}+\alpha$.
  4. All New: Both Wiescher et al. (2023) and Best et al. (2013) rates applied.

3. Key Contributions

• Systematic Comparison: The first comprehensive study comparing the combined effects of updated $^{17}\text{O}+\alpha$ and $^{22}\text{Ne}+\alpha$ rates on the ws-process in non-rotating metal-poor massive stars.
• Rate Sensitivity Analysis: Quantified how specific updates in the $(\alpha, n)/(\alpha, \gamma)$ ratios at different temperatures (He, C, Ne, and O burning stages) impact neutron densities and final yields.
• Mass Dependence: Investigated how the enhancement of ws-process yields scales with stellar mass ($15-30 M_\odot$).
• Production Factors: Calculated integrated production factors (PFs) using a Salpeter Initial Mass Function (IMF) to assess the contribution of these stars to the chemical enrichment of the Galaxy.

4. Key Results

Impact of Reaction Rates on Neutron Density

• $^{17}\text{O}+\alpha$ (Best et al. 2013):
  • Below $0.7$ GK (He burning), rates are similar to REACLIB.
  • Above $0.7$ GK (C, Ne, O burning), the $(\alpha, n)/(\alpha, \gamma)$ ratio increases dramatically (by tens of times) compared to REACLIB.
  • Effect: Significantly increases neutron release throughout all burning stages (He, C, Ne), reducing neutron consumption by $^{16}\text{O}$.
• $^{22}\text{Ne}+\alpha$ (Wiescher et al. 2023):
  • The $(\alpha, n)/(\alpha, \gamma)$ ratio is 1.2–2.0 times higher than REACLIB below $1.5$ GK, but increases by tens of times above $1.5$ GK.
  • Effect: Primarily enhances neutron production during C and Ne burning stages. It suppresses the $^{22}\text{Ne}(\alpha, \gamma)^{26}\text{Mg}$ channel, preserving $^{22}\text{Ne}$ for neutron release later in the evolution.

Isotopic Yields ($Z=31-40$, "Ga-Zr")

• Magnitude of Enhancement:
  • Default: Baseline yields.
  • New $^{17}\text{O}$: Increases total $X(\text{Ga-Zr})$ by a factor of ~4.8.
  • New $^{22}\text{Ne}$: Increases total $X(\text{Ga-Zr})$ by a factor of ~23.8.
  • Combined (All New): Increases total $X(\text{Ga-Zr})$ by a factor of ~113.6 (over two orders of magnitude).
• Burning Stage Contribution:
  • With default rates, ~51% of production occurs in He burning and 41% in C burning.
  • With new $^{22}\text{Ne}$ rates, the contribution shifts heavily to C and Ne burning (89% combined) because $^{22}\text{Ne}$ is not exhausted until late stages.
  • With new $^{17}\text{O}$ rates, production increases across all stages, but the relative contribution percentages remain similar to the default case.

Mass Distribution and Stellar Mass Dependence

• Spatial Distribution: Without mixing, the abundance profile shows a "two-bump" structure: one peak in the C/Ne/O burning shells ($M_r \approx 2.2-3.6 M_\odot$) and another in the He burning shell ($M_r \approx 5.0-5.9 M_\odot$).
• Mass Scaling: The enhancement due to new $^{17}\text{O}$ rates becomes more significant in more massive stars ($30 M_\odot$).
• Heavier Elements ($Z > 40$): The new rates have minimal impact on elements heavier than Zr (Nb-Th), which are produced primarily during He burning where the new rates have less effect.

Production Factors (Galactic Enrichment)

• Using a Salpeter IMF, the study found that the new $^{17}\text{O}$ rates significantly increase the production factors for elements from Zn to Rb (by $>0.5$ dex).
• The combined new rates lead to an overproduction of Ga, Ge, As, and Se relative to solar abundances in metal-poor models, suggesting that either the rates need further refinement or metallicity dependence must be considered (e.g., models with $Z=0.5 Z_\odot$).

5. Significance and Conclusions

1. Resolution of Discrepancy: The study suggests that updated nuclear reaction rates, particularly for $^{17}\text{O}+\alpha$, can resolve the discrepancy between predicted low ws-process yields in metal-poor stars and observational data.
2. Dominance of $^{17}\text{O}$: While $^{22}\text{Ne}$ is the primary neutron source, the efficiency of the ws-process in metal-poor stars is critically dependent on the $^{17}\text{O}+\alpha$ rates, which determine whether neutrons are recycled or lost to $^{16}\text{O}$.
3. Experimental Urgency: The results show that different literature values for $^{17}\text{O}+\alpha$ rates can cause order-of-magnitude differences in yields, whereas $^{22}\text{Ne}+\alpha$ variations cause factors of 3–5. This highlights an urgent need for precise experimental measurements of $^{17}\text{O}(\alpha, n)^{20}\text{Ne}$ and $^{17}\text{O}(\alpha, \gamma)^{21}\text{Ne}$ across the temperature range of 0.1–10 GK.
4. Modeling Limitations: The authors note that their post-processing approach neglects convective mixing, which could further enhance yields by transporting fresh fuel ($^{22}\text{Ne}$, $^{17}\text{O}$) into burning shells. Coupled evolution and nucleosynthesis calculations are recommended for future work.

In summary, the paper demonstrates that nuclear physics uncertainties in light-element reactions are the dominant factor governing the weak s-process in metal-poor massive stars, with the potential to increase heavy element production by tens to over a hundred times depending on the adopted rates.