Microscopic study of nuclei synthesis in pycnonuclear reaction $^{12}$C + $^{12}$C in neutron stars

This paper employs a microscopic cluster model with folding potentials and the Multiple Internal Reflections method to demonstrate that the synthesis of $^{24}$Mg via $^{12}$C + $^{12}$C pycnonuclear reactions in neutron stars is most probable through the formation of a new excited compound nucleus in quasibound states, offering a more precise description than traditional Woods-Saxon potentials.

The Gist

The Setting: A Cosmic Pressure Cooker

Imagine a neutron star. It's not just a star; it's a cosmic monster. It's so heavy that a teaspoon of its material would weigh as much as a mountain. Inside, the atoms are crushed together so tightly that they are practically touching, like sardines in a can that has been squeezed until the cans themselves are deforming.

In this environment, a special kind of cooking happens called a pycnonuclear reaction. "Pycno" means dense. Unlike normal stars, which need heat to cook atoms together, neutron stars are cold. Instead, they use pressure. The atoms are squeezed so close that they are forced to collide and fuse, creating new, heavier elements.

The Problem: The Old Map vs. The New GPS

For decades, scientists tried to predict what happens when two Carbon-12 atoms (the building blocks of life) crash into each other in this pressure cooker. They used an old map called the Woods-Saxon potential.

Think of this old map like a blurry, low-resolution photo of a mountain. It tells you roughly where the peak is and where the valley is, but it misses the tiny details. It assumes the atoms are smooth, featureless balls.

This paper says: "That map is wrong for this terrain."

The authors, led by S. P. Maydanyuk and his team, decided to build a high-definition 3D GPS. They didn't treat the Carbon atoms as smooth balls. Instead, they looked inside the atoms, seeing them as clusters of smaller particles (protons and neutrons) arranged in a specific, wobbly structure. They used a method called the Folding Approximation, which is like taking a detailed blueprint of the atom's internal structure and "folding" it into the calculation to see exactly how the atoms interact when they get close.

The Discovery: The "Hidden Room" in the Castle

When two Carbon atoms approach each other in a neutron star, they have to climb a giant energy hill (the Coulomb barrier) to fuse. Usually, they bounce off or tunnel through slowly.

The authors used a new mathematical technique called the Method of Multiple Internal Reflections. Imagine a hallway with mirrors at both ends. If you shout, the sound bounces back and forth. The authors tracked these "echoes" of the quantum waves inside the atoms.

Here is the big surprise:
They found that at certain specific energy levels, the atoms don't just bounce or tunnel; they get "stuck" in a sweet spot. They call these Quasi-Bound States.

• The Old View: The atoms vibrate a little (zero-point vibration) and might occasionally fuse. It's a slow, rare event.
• The New View: There are specific "parking spots" (resonance states) where the atoms lock together with incredible stability.

The Analogy:
Imagine trying to push two magnets together.

• Old Theory: You push them, they repel, and sometimes, if you push hard enough, they snap together.
• New Theory: The authors found that if you push them at exactly the right rhythm, they snap into a "magnetic embrace" that is 10,000,000,000,000,000,000,000,000,000,000 times more likely to happen than the old theory predicted.

The Result: A New Kind of Magnesium

When two Carbon-12 atoms fuse, they become Magnesium-24.

The paper concludes that in the dense heart of a neutron star, this fusion doesn't happen randomly. It happens most efficiently when the atoms are in these special Quasi-Bound States.

1. It's much more likely: The probability of creating Magnesium-24 in these states is astronomically higher than previously thought.
2. It's safer: At the first energy level, the new Magnesium nucleus has a "shield" (a barrier) that prevents it from falling apart immediately. It's a stable, new baby nucleus.
3. The Map Changed: The "Folding" method (the high-res GPS) showed that the energy levels where this happens are different from the old "Woods-Saxon" map. If you use the old map, you might miss the party entirely.

Why Does This Matter?

Neutron stars are the universe's ultimate laboratories. They create the heavy elements that make up planets and people.

• Before: We thought the creation of elements in these stars was a slow, grinding process based on rough estimates.
• Now: We know there are specific "sweet spots" where the universe's chemistry accelerates dramatically.

The authors have essentially told us that the universe isn't just a chaotic smash-up; it's a highly tuned machine. If you know the exact rhythm (the quasi-bound energy), you can predict exactly how the stars cook up the ingredients for the next generation of planets.

In a Nutshell

The scientists took a microscope to the heart of a neutron star. They realized that the old way of calculating how atoms fuse was like looking at a cloud and guessing its shape. They built a new model that sees the individual water droplets. They discovered that the atoms have "secret handshake" energies where they fuse together almost instantly, creating new elements much faster and more efficiently than we ever imagined.

Technical Summary

1. Problem Statement

The paper addresses the mechanism of pycnonuclear reactions (nuclear reactions driven by high density rather than temperature) occurring in the dense crusts of neutron stars. Specifically, it investigates the synthesis of heavier nuclei, focusing on the $^{12}\text{C} + ^{12}\text{C} \to ^{24}\text{Mg}$ reaction.

Key Challenges Identified:

• Limitations of Previous Models: Traditional studies often rely on the Woods-Saxon potential and semiclassical approximations. These models assume nuclei collide from asymptotic distances (relevant for terrestrial beam experiments) and often ignore the complex quantum interference effects occurring when nuclei are already in close proximity within a stellar lattice.
• Zero-Point Energy Assumption: Previous estimates of reaction rates often relied on the energy of zero-point vibrations (ground state of a harmonic oscillator near the lattice saddle point).
• Lack of Microscopic Precision: There is a need for a fully microscopic description of the interaction between nuclei at close distances (tens of femtometers) without relying on phenomenological parameters fitted to experimental scattering cross-sections, which are unavailable for dense stellar matter.

2. Methodology

The authors employ a rigorous, fully microscopic approach combining the Cluster Model with the Method of Multiple Internal Reflections (MIR).

A. Theoretical Framework

1. Cluster Model & Folding Approximation:

   • The system is treated as two interacting $^{12}\text{C}$ clusters.
   • Wave Functions: The internal structure of $^{12}\text{C}$ is described using a many-particle shell model with an oscillator potential (Elliott quantum numbers $(04)$), representing an oblate shape.
   • Interaction Potential: Instead of a phenomenological Woods-Saxon potential, the authors derive the folding potential ($V^{(F)}$). This is calculated by folding a semi-realistic nucleon-nucleon (NN) interaction (Gaussian form) with the density distributions of the two $^{12}\text{C}$ nuclei.
   • Two Forms of Folding Potential:
     • S-form: A simplified scalar approximation.
     • F-form: A more rigorous calculation including tensor components and the full p-shell structure, which reveals a unique double-barrier structure.

2. Method of Multiple Internal Reflections (MIR):

   • This quantum mechanical method solves the Schrödinger equation for nuclei starting from close distances (typical of lattice sites in neutron stars, $R_0 \approx 92.5$ fm) rather than infinity.
   • The potential barrier is discretized into rectangular steps.
   • The method accounts for interference between incoming and outgoing fluxes within the nuclear region, a feature often neglected in standard fusion models.
   • It calculates the probability of compound nucleus formation ($P_{cn}$) by integrating the probability density between internal turning points.

3. Resonance State Detection:

   • The authors use the stabilization method (analyzing eigenenergies as a function of basis size) and analyze phase shifts, oscillation coefficients, and internal wave function weights to identify quasibound states.

3. Key Contributions

1. Development of the F-form Folding Potential: The paper introduces a new, highly accurate folding potential (F-form) that incorporates the p-shell structure of $^{12}\text{C}$. This potential reveals a second, smaller barrier ($V \approx -56.6$ MeV at $r=3.12$ fm) inside the interaction region, a phenomenon previously unobserved in nuclear reaction potentials.
2. Microscopic Derivation without Phenomenology: The interaction parameters are derived strictly from two-nucleon interactions and shell-model wave functions, eliminating the need to fit unknown parameters to experimental data (which is impossible for neutron star conditions).
3. Quantum Flux Interference: The study demonstrates that interference between ingoing and outgoing fluxes in the dense medium plays a dominant role in determining fusion probabilities, a factor missed by semiclassical approaches.
4. Identification of Quasibound States: The authors establish that the most probable synthesis of $^{24}\text{Mg}$ occurs not at zero-point vibration energies, but at specific quasibound state energies where the compound nucleus formation probability is maximized.

4. Key Results

A. Potential Structure and Resonances

• Potential Comparison: The folding potentials (S-form and F-form) differ significantly from the Woods-Saxon potential, particularly in the depth of the potential well and the barrier structure.
• Resonance States: A spectrum of resonance states for $^{24}\text{Mg}$ was calculated (up to $L=20$).
  • States split into two rotational bands.
  • Narrow states: Widths $< 20$ keV, average cluster distance $\approx 1.9$ fm.
  • Wide states: Widths $\approx 2.0$ MeV, average cluster distance $\approx 5.0$ fm.
  • The interacting $^{12}\text{C}$ nuclei behave as a rigid rotating body.

B. Quasibound States and Synthesis Probability

• New Maxima: The probability of compound nucleus formation ($P_{cn}$) exhibits clear maxima at specific energies (quasibound states).
• Energy Shifts: The energies of these quasibound states differ significantly depending on the potential used:
  • Woods-Saxon: $E_1 \approx 4.88$ MeV.
  • Folding S-form: $E_1 \approx 2.49$ MeV.
  • Folding F-form: $E_1 \approx 3.49$ MeV.
• Barrier Penetration: For the first quasibound energy, the energy is lower than the barrier maximum for all potentials. This implies the compound nucleus is formed in a state where a barrier prevents immediate decay via tunneling, effectively creating a stable (or long-lived) new nucleus $^{24}\text{Mg}$.

C. Comparison with Zero-Point Vibrations

• The probability of forming a compound nucleus in these quasibound states is orders of magnitude higher (e.g., ratio $\sim 10^{30}$) than in states of zero-point vibrations.
• This suggests that previous estimates of pycnonuclear reaction rates based solely on zero-point energies may be significantly underestimated or physically inaccurate regarding the dominant reaction channel.

5. Significance and Conclusions

• Revised Picture of Stellar Nucleosynthesis: The study fundamentally alters the understanding of how heavy nuclei are synthesized in neutron stars. It posits that synthesis is driven by quasibound resonant states rather than simple zero-point tunneling.
• Methodological Advancement: The combination of the folding approximation (for microscopic interaction) and the MIR method (for close-distance scattering) provides a high-precision ($10^{-14}$ accuracy in flux conservation) tool for studying nuclear processes in dense matter.
• Implications for Neutron Star Evolution: Since the $^{12}\text{C} + ^{12}\text{C}$ reaction is critical for the ignition of Type Ia supernovae and the evolution of white dwarfs, the identification of these specific quasibound states and their higher formation probabilities could lead to revised models of stellar heating, cooling, and explosion mechanisms.
• Future Outlook: The paper advocates for the development of similar microscopic quantum mechanical methods to study other nuclear reactions in dense media, moving away from phenomenological potentials toward first-principles calculations based on nucleon-nucleon interactions.