Gamow-Teller strength of $^{12,14,16}$C within deformed quasiparticle random-phase approximation

This study employs the deformed quasiparticle random-phase approximation (DQRPA) with Brückner $G$-matrix interactions to analyze Gamow-Teller transition strengths in carbon isotopes $^{12,14,16}$C, revealing that nuclear deformation significantly influences the strength distributions in $^{12}$C and $^{16}$C while a spherical limit successfully reproduces experimental data for $^{14}$C.

The Gist

Imagine the atomic nucleus not as a static, hard marble, but as a bustling, dynamic dance floor filled with tiny dancers (protons and neutrons). Sometimes, these dancers need to swap partners or change their rhythm to keep the whole system stable. This paper is about studying a very specific type of "dance move" called the Gamow-Teller (GT) transition.

Here is the story of what the researchers did, explained simply:

1. The Big Picture: Why Do We Care?

In the universe, stars are like giant nuclear furnaces. They burn fuel, create new elements, and sometimes explode as supernovas. To understand how stars live and die, we need to know how their nuclei react to neutrinos (ghostly particles that pass through everything).

The "Gamow-Teller transition" is the most common way a nucleus reacts to these neutrinos. If we can predict exactly how strong this reaction is, we can better understand:

• How stars explode.
• How heavy elements are made.
• How to date ancient artifacts (using Carbon-14).

2. The Problem: The "Shape" Matters

For a long time, physicists treated atomic nuclei like perfect spheres (like billiard balls). But the researchers in this paper realized that for some light atoms, like Carbon, this isn't true.

The Analogy: Imagine a group of people holding hands in a circle.

• Spherical Model: They stand in a perfect circle.
• Deformed Model: They stretch out into an oval or a rugby ball shape.

The paper focuses on three Carbon isotopes (versions of Carbon with different numbers of neutrons): Carbon-12, Carbon-14, and Carbon-16. The researchers asked: "Does the shape of the nucleus change how these dancers perform their moves?"

3. The Tools: A High-Tech Simulation

To answer this, they built a super-complex computer simulation called DQRPA (Deformed Quasiparticle Random-Phase Approximation).

• The Mean Field: This is the "floor" the dancers stand on. They used a realistic map of the forces inside the nucleus.
• The Residual Interaction: This is the "music" or the "rules" that tell the dancers how to react to each other when they move. They used a very precise set of rules derived from the CD-Bonn potential (a fancy way of saying they used a very accurate map of how protons and neutrons talk to each other).

4. The Findings: What Happened to the Three Carbons?

Carbon-12: The Shape-Shifter

• The Mystery: Experiments showed a specific "dance move" happening at a very low energy. Standard models (treating it as a sphere) predicted the move would happen at a higher energy, missing the mark by about 3 units.
• The Solution: The researchers realized that if they weakened the "spin-orbit" force (a force that keeps the dancers spinning in a specific way) and allowed the nucleus to become deformed (squashed or stretched), the simulation matched the experiment perfectly.
• The Takeaway: Carbon-12 is sensitive to its shape. If you ignore the deformation, you get the physics wrong.

Carbon-14: The Perfect Match

• The Mystery: Carbon-14 is famous for radiocarbon dating. It has two distinct "dance moves" (peaks in energy).
• The Solution: When the researchers treated Carbon-14 as a sphere (which it mostly is), their simulation perfectly matched the real-world data. They found that the interaction between the dancers (specifically the "particle-particle" and "particle-hole" forces) pushes the second move to the right energy level.
• The Takeaway: Sometimes, the simple spherical model works great, but you still need the right "music" (interaction forces) to get the timing right.

Carbon-16: The High-Flyer

• The Mystery: Carbon-16 has extra neutrons, making it "open-shell" (like a dance floor with extra space).
• The Solution: Because Carbon-16 is deformed (stretched out), the extra neutrons can jump to much higher energy levels. The simulation predicted a whole new set of "high-flying" dance moves above 15 MeV that wouldn't exist if the nucleus were a perfect sphere.
• The Takeaway: Deformation creates new possibilities. It mixes the dancers' configurations, leading to high-energy reactions that we haven't seen in experiments yet (because we haven't looked hard enough), but the theory predicts they are there.

5. The "Spin-Orbit" Secret Sauce

A key discovery in this paper is the role of the spin-orbit force.

• Analogy: Imagine the dancers are spinning tops. The spin-orbit force is the friction that keeps them spinning upright.
• Discovery: In light nuclei like Carbon, if you reduce this "friction" (spin-orbit strength), the nucleus changes shape, and the energy levels of the dance moves shift. This explains why Carbon-12 behaves the way it does.

Summary

This paper is like a detective story where the scientists used a high-tech simulation to figure out why Carbon atoms behave the way they do. They found that:

1. Shape matters: Carbon-12 isn't a perfect ball; it's squashed, and that changes its behavior.
2. Forces matter: The specific rules of how neutrons and protons interact are crucial for getting the energy levels right.
3. Prediction: They predicted new, high-energy reactions in Carbon-16 that future experiments might find.

By understanding these tiny "dance moves," we get a clearer picture of how stars work and how the elements in our universe were created.

Technical Summary

1. Problem Statement

The study addresses the theoretical description of Gamow-Teller (GT) transition strength distributions in light carbon isotopes ($^{12,14,16}\text{C}$). While GT transitions are crucial for understanding weak interaction processes (such as neutrino-nucleus reactions relevant to supernovae) and nuclear structure, several challenges remain:

• Deformation Effects: Many light nuclei exhibit deformation, which alters shell evolution and particle-hole ($p-h$) and particle-particle ($p-p$) interactions. Standard spherical models often fail to capture these effects accurately.
• Spin-Orbit Interaction: In light $p$-shell nuclei, the spin-orbit coupling is moderate, leading to intermediate coupling schemes where deformation may arise from reduced spin-orbit strength.
• Residual Interactions: Accurate modeling requires realistic residual interactions derived from nucleon-nucleon potentials, rather than phenomenological adjustments, to explain the fragmentation and energy shifts of GT peaks.
• Specific Challenges:
  • $^{12}\text{C}$: Theoretical models often overestimate the excitation energy of the low-lying GT peak compared to experimental $(p,n)$ data.
  • $^{14}\text{C}$: The distribution features two distinct peaks, requiring precise tuning of residual interactions.
  • $^{16}\text{C}$: As an open-shell nucleus with extra neutrons, it may exhibit high-lying GT strength and deformation-induced configuration mixing not seen in lighter isotopes.

2. Methodology

The authors employ a Deformed Quasiparticle Random Phase Approximation (DQRPA) framework with the following specific components:

• Mean Field: The single-particle states and energies are derived from Deformed Skyrme Hartree-Fock-Bogoliubov (DSHFB) calculations using a harmonic oscillator basis. Three Skyrme parameter sets were tested: SLy4, SkP, and SGII.
• Residual Interaction: Unlike standard Skyrme-based QRPA where the residual force is derived from the same energy density functional, this study uses a realistic residual interaction derived from the Brueckner $G$-matrix based on the CD-Bonn potential. This interaction is calculated in the spherical basis and transformed into the deformed basis.
• Pairing: Pairing correlations are treated within the BCS approximation using a density-dependent pairing force. The pairing strength ($V_0$) was adjusted to reproduce the GT peak position in $^{18}\text{O}$ and to ensure finite pairing gaps in the open-shell $^{16}\text{C}$.
• Key Variations:
  • Spin-Orbit Strength: For $^{12}\text{C}$, the spin-orbit strength ($W_0$) was reduced to 60% of its standard value to investigate its impact on deformation and shell structure.
  • Interaction Strengths: The particle-hole ($g_{ph}$) and particle-particle ($g_{pp}$) renormalization factors were systematically varied to examine their impact on GT strength distributions.
• Formalism: The calculation utilizes the equation of motion for phonon operators acting on a BCS ground state, explicitly treating axial deformation ($K$ quantum number) and parity.

3. Key Contributions

• Hybrid Framework: The paper successfully implements a hybrid approach combining a Skyrme mean field with a realistic $G$-matrix residual interaction, allowing for the inclusion of finite-range correlations beyond standard Skyrme QRPA.
• Deformation in Light Nuclei: It explicitly demonstrates how nuclear deformation, driven by reduced spin-orbit strength in $^{12}\text{C}$ and open-shell effects in $^{16}\text{C}$, significantly alters GT strength distributions.
• Systematic Analysis of Residual Interactions: The study isolates the effects of $p-p$ and $p-h$ interactions, showing that the repulsive $p-p$ interaction is critical for shifting GT peaks to match experimental data, particularly in $^{12}\text{C}$ and $^{14}\text{C}$.
• Prediction for $^{16}\text{C}$: It provides the first theoretical prediction of GT strength for $^{16}\text{C}$, identifying high-lying strength associated with deformation-induced configuration mixing.

4. Results

$^{12}\text{C}$

• Deformation: Standard Skyrme forces predict a spherical ground state. However, reducing the spin-orbit strength to 60% induces an oblate deformation ($\beta_2 \approx -0.34$).
• GT Strength:
  • In the spherical limit, the calculated GT peak is shifted $\sim 3$ MeV higher than experimental data.
  • The reduced spin-orbit strength enhances configuration mixing between $0p_{3/2}$ and $0p_{1/2}$ orbitals. In the BCS limit, this places the GT peak near $E_x = 0$ MeV, consistent with data.
  • However, the inclusion of the residual interaction (QRPA) shifts the peak upward again by a few MeV.
  • The repulsive $p-p$ interaction ($g_{pp} > 1.0$) is found to shift the low-lying peaks to higher energies, while the $p-h$ interaction has a minimal shifting effect in this light deformed system.

$^{14}\text{C}$

• Structure: The nucleus remains spherical ($\beta_2 \approx 0$) across all parameter sets.
• GT Strength: The experiment shows two peaks: a weak one near the ground state and a strong one at $\sim 4$ MeV.
  • Theoretical calculations reproduce the two-peak structure.
  • The second peak (dominated by the $\nu 0p_{3/2} \to \pi 0p_{1/2}$ configuration) is sensitive to residual interactions. Increasing both $g_{pp}$ and $g_{ph}$ by $\sim 30\%$ shifts this peak to higher excitation energies, achieving excellent agreement with experimental $(p,n)$ data.
  • The main configuration for the low-lying state is $\nu 0p_{1/2} \to \pi 0p_{1/2}$, while the high-lying state is dominated by $\nu 0p_{3/2} \to \pi 0p_{1/2}$.

$^{16}\text{C}$

• Deformation: The SGII parameter set predicts a prolate deformation ($\beta_2 \approx 0.14$), whereas SLy4 and SkP predict spherical shapes.
• GT Strength Distribution:
  • Low-lying peaks: Similar to $^{14}\text{C}$, arising from $p$-shell configurations.
  • High-lying strength: A distinct feature is the appearance of significant GT strength above 15 MeV. This is attributed to deformation-induced configuration mixing involving the $0d_{5/2}$ and $0d_{3/2}$ orbitals ($\nu 0d_{5/2} \to \pi 0d_{5/2}$, etc.).
  • Ikeda Sum Rule (ISR): The deformed SGII case shows a reduced ISR exhaustion ($\sim 88\%$) compared to spherical cases ($\sim 96-99\%$), indicating strength fragmentation into high-lying states and continuum configurations.

5. Significance

• Astrophysical Relevance: The accurate modeling of GT strengths in carbon isotopes is vital for calculating neutrino-nucleus cross-sections, which govern energy transport and nucleosynthesis in supernovae and stellar evolution.
• Nuclear Structure Insight: The results highlight that deformation and spin-orbit coupling are not negligible even in light $p$-shell nuclei. The study validates the necessity of treating deformation consistently within the QRPA framework.
• Model Validation: The success of the DQRPA with realistic $G$-matrix interactions in reproducing experimental data (especially the peak positions in $^{14}\text{C}$ and the deformation effects in $^{12}\text{C}$) suggests this hybrid approach is a robust tool for studying weak interaction processes in nuclei where full self-consistency is computationally challenging.
• Future Directions: The authors note that the current formalism excludes $T=0$ (isoscalar) pairing, which is expected to be significant in $N \approx Z$ nuclei like $^{12}\text{C}$. Including $T=0$ pairing is identified as a critical next step for a complete description of low-lying GT states.