Towards a microscopic description of 12C+12C fusion at stellar energies

This paper presents a fully microscopic multichannel Resonating Group Method description of the 12C+12C fusion reaction at stellar energies, which successfully reproduces experimental scattering data, reveals the mixed nature of 24Mg resonances, predicts fusion hindrance at low energies, and lays the groundwork for reliable theoretical extrapolations to deep stellar burning temperatures.

The Gist

Imagine two massive, positively charged magnets trying to hug each other. Because they both have the same charge, they push each other away fiercely. This is the Coulomb barrier. In the universe, this happens when two Carbon-12 atoms (the building blocks of life and stars) try to smash together to form a heavier element, Magnesium-24. This process, called fusion, is the engine that powers massive stars and creates the heavy elements we are made of.

The problem? It's incredibly hard to make them touch. They repel each other so strongly that they only fuse at extremely high temperatures, deep inside stars. Trying to study this in a lab on Earth is like trying to hear a whisper in a hurricane; the signals are tiny, and the data is messy.

This paper by Pierre Descouvemont is like building a super-accurate, microscopic simulation of that hug to figure out exactly how it happens.

The Old Way vs. The New Way

The Old Way (The Single-Channel Approximation):
Imagine trying to understand a complex dance by watching only the lead dancer. Previous models looked at the two Carbon atoms just as two solid balls bouncing off each other. They ignored the fact that the atoms are actually made of smaller parts (protons and neutrons) that can rearrange themselves. It was like assuming the dancers never change their steps or partners. This led to predictions that didn't quite match reality.

The New Way (The Microscopic Multichannel Model):
Descouvemont's approach is like putting on X-ray glasses and watching the entire dance floor.

• The Microscopic View: Instead of treating the Carbon atoms as solid balls, he treats them as collections of 12 tiny particles (nucleons) each.
• The Multichannel View: He realizes that when two Carbon atoms get close, they don't just bounce or fuse directly. They can temporarily break apart and swap pieces. For example, one Carbon might spit out a tiny "alpha particle" (a helium nucleus), leaving behind a Neon-20 atom. The two atoms might dance as Carbon + Carbon or as Alpha + Neon.
• The Analogy: Imagine two couples dancing. In the old model, they just bumped into each other. In this new model, the dancers might suddenly swap partners mid-dance, forming two new couples for a split second before coming back together. The model accounts for all these possible partner-swaps.

What Did They Discover?

1. The "Molecular" Myth is Broken
For decades, scientists thought that when Carbon atoms fused, they formed "molecular states"—like two distinct molecules sticking together in a neat, predictable pattern.

• The Discovery: The simulation shows this isn't true. The resulting Magnesium-24 atom is a chaotic soup. The wave functions (the "dance moves" of the particles) are a messy mix of many different configurations. It's not a neat molecule; it's a highly mixed, complex state where the particles are constantly rearranging.

2. The Fusion "Traffic Jams" (Resonances)
The paper predicts that at certain energies, the fusion process speeds up dramatically because of resonances.

• The Analogy: Think of pushing a child on a swing. If you push at the exact right moment (the resonance), the swing goes huge with very little effort. If you push at the wrong time, nothing happens.
• The model finds that at specific energies near the "Coulomb barrier," the stars get a "perfect push." These are narrow and broad resonances that act like shortcuts, allowing fusion to happen more easily than previously thought.

3. The "Fusion Hindrance" Mystery
There is a debate in the scientific community: Does fusion get harder at very low energies (a phenomenon called "hindrance")?

• The Discovery: The model shows that the fusion rate (the "S-factor") does indeed drop off at very low energies. This supports the idea that fusion gets "stuck" or hindered when the energy is too low to overcome the repulsion, even with the help of quantum tunneling.

Why Does This Matter?

This work is a first step toward a perfect map of stellar evolution.

• For Stars: It helps us understand how massive stars burn their fuel and eventually explode as supernovae.
• For the Universe: It tells us how heavy elements (like the calcium in your bones or the iron in your blood) were forged in the fires of ancient stars.
• For the Future: While this model is a huge leap forward, the author admits it's still a "draft." The next step is to include even more complex "partner swaps" (involving protons and neutrons) to get a complete picture.

In a Nutshell:
Descouvemont built a high-definition, 3D simulation of two Carbon atoms trying to fuse. He found that they don't just bounce or stick neatly; they dance in a chaotic, multi-partner routine full of temporary breaks and swaps. This new understanding helps us finally decode the secret recipe of how stars cook the heavy elements of the universe.

Technical Summary

1. Problem Statement

The $^{12}\text{C} + ^{12}\text{C}$ fusion reaction is a critical process in stellar nucleosynthesis, driving carbon burning in massive stars and triggering the production of heavy elements. However, determining the fusion cross-section at stellar energies (typically $\approx 2.4$ MeV) is extremely challenging due to:

• Coulomb Repulsion: The small cross-sections at low energies make experimental measurement difficult.
• Resonance Uncertainty: The presence of resonances near the Gamow window is well-established but their origin and impact on reaction rates remain debated. Recent claims of low-energy resonances (e.g., via the Trojan Horse method) have been contested.
• Fusion Hindrance: There is ongoing debate regarding whether fusion hindrance (a suppression of the reaction rate at very low energies) occurs in this system.
• Theoretical Limitations: Existing models often rely on the optical model (phenomenological potentials) or time-dependent methods that cannot simultaneously describe fusion, elastic scattering, and the spectroscopy of the compound nucleus ($^{24}\text{Mg}$). Furthermore, many microscopic approaches use single-channel approximations or bound-state assumptions that fail to account for scattering conditions and channel interference.

2. Methodology

The author employs a fully microscopic, multichannel Resonating Group Method (RGM) combined with the Generator Coordinate Method (GCM). Key aspects of the methodology include:

• Microscopic Framework: The calculation uses a nucleon-nucleon interaction (Volkov V2 with Bartlett and Heisenberg terms) as the sole input, avoiding phenomenological imaginary potentials.
• Multichannel Approach: The wave function explicitly includes:
  • The entrance channel: $^{12}\text{C}(\text{g.s.}) + ^{12}\text{C}(\text{g.s.})$.
  • Inelastic channels: Excited states of $^{12}\text{C}$ ($2^+, 4^+$).
  • Transfer/Fusion exit channels: $\alpha + ^{20}\text{Ne}$ (including ground and excited states of $^{20}\text{Ne}$).
  • Note: Proton ($^{23}\text{Na} + p$) and neutron ($^{23}\text{Mg} + n$) channels are omitted due to computational constraints but are acknowledged as future work.
• Cluster Wave Functions:
  • $^{12}\text{C}$ and $^{20}\text{Ne}$ clusters are described using the shell model.
  • $^{12}\text{C}$ utilizes all $p$-shell configurations ($N=225$ Slater determinants).
  • $^{20}\text{Ne}$ utilizes $sd$-shell configurations for the valence nucleons ($N=4356$ Slater determinants).
  • The total basis size for the $^{12}\text{C}+^{12}\text{C}$ partition reaches $225^2 = 50,625$ Slater determinants.
• Scattering Formalism: The RGM wave functions are adapted for scattering problems using the microscopic R-matrix method to handle asymptotic boundary conditions and extract the scattering matrix ($U^J$).
• Fusion Cross-Section: Calculated via the transfer to $\alpha + ^{20}\text{Ne}$ channels, defined by the imaginary part of the scattering matrix elements.

3. Key Contributions

• First Fully Microscopic Multichannel Calculation: This is the first study to simultaneously treat $^{12}\text{C}+^{12}\text{C}$ fusion, elastic scattering, and $^{24}\text{Mg}$ spectroscopy using a fully microscopic approach that includes dominant open transfer channels ($\alpha + ^{20}\text{Ne}$) without phenomenological absorption potentials.
• Resolution of "Molecular State" Concept: The study challenges the traditional view of "pure molecular states" (dominated by a single $^{12}\text{C}+^{12}\text{C}$ configuration) in $^{24}\text{Mg}$.
• Consistent Description: It provides a unified theoretical framework where the same wave functions describe bound states, resonances, and scattering data.

4. Key Results

A. Elastic Scattering

• The multichannel model reproduces experimental elastic scattering data (ratios to Mott cross-sections) with excellent agreement, particularly at energies above the Coulomb barrier (6–10 MeV).
• Single-channel approximations fail to reproduce the oscillatory structure of the data, demonstrating that virtual couplings to excited and transfer channels are essential for accurate scattering descriptions.

B. $^{24}\text{Mg}$ Spectroscopy

• Wave Function Structure: Analysis reveals that $^{24}\text{Mg}$ states (both bound and resonant) are highly mixed configurations. The $^{12}\text{C}(0^+_1) + ^{12}\text{C}(0^+_1)$ component contributes at most ~10% to the wave function.
• Rejection of Pure Molecular States: The study contradicts the concept of pure "molecular states" found in single-channel calculations. Instead, the states are spread over multiple channels, including significant $\alpha + ^{20}\text{Ne}$ components.
• Resonances: The model predicts several resonances above the $^{12}\text{C}+^{12}\text{C}$ threshold, including $0^+$, $2^+$, and $4^+$ states.

C. Fusion Cross-Section and S-Factor

• Agreement with Data: The calculated fusion S-factor is consistent with available experimental data.
• Resonance Characteristics:
  • The model predicts both narrow ($2^+, 4^+$) and broad ($0^+, 2^+$) resonances near the Coulomb barrier.
  • Width Origin: The resonance widths are primarily determined by the $\alpha + ^{20}\text{Ne}$ exit channels, not the $^{12}\text{C}+^{12}\text{C}$ entrance channel (which has extremely small partial widths due to the Coulomb barrier).
• Fusion Hindrance: The calculated S-factor exhibits a decrease at low energies (below $\approx 2$ MeV). This provides microscopic support for the hypothesis of fusion hindrance, although the conclusion is currently limited to the $\alpha$-channel and requires verification with proton/neutron channels.

5. Significance and Future Outlook

• Astrophysical Impact: This work offers a reliable theoretical foundation for extrapolating the $^{12}\text{C}+^{12}\text{C}$ reaction rate to deep stellar burning temperatures, a critical input for stellar evolution models.
• Methodological Advancement: It demonstrates that modern computational capabilities allow for solving many-body scattering problems ($A=24$) with high precision, moving beyond phenomenological models.
• Future Directions: The author identifies the inclusion of proton ($^{23}\text{Na} + p$) and neutron ($^{23}\text{Mg} + n$) channels as the necessary next step to provide a complete description of the fusion process and to definitively confirm the fusion hindrance effect. Additionally, including negative-parity states of $^{20}\text{Ne}$ is suggested for further refinement.

In summary, Descouvemont presents a breakthrough in nuclear astrophysics by providing a consistent, microscopic, and multichannel description of the $^{12}\text{C}+^{12}\text{C}$ reaction, resolving long-standing discrepancies in scattering data and offering new insights into the structure of $^{24}\text{Mg}$ and the nature of fusion resonances.