Constraining the $N=16$ Shell Gap in $^{17}$C via Transfer to the Continuum in the $^{16}$C$(d,p)^{17}$C Reaction

This study extends a semi-microscopic reaction model to unbound states in the $^{16}$C(d,p)$^{17}$C reaction, demonstrating that the measured differential cross sections require a large $N=16$ shell gap exceeding 5 MeV to accurately describe the $1d_{3/2}$ single-particle strength.

The Gist

Imagine the atomic nucleus not as a solid ball, but as a bustling city. In this city, protons and neutrons are the citizens, and they live in specific "neighborhoods" called shells. Just like in a real city, there are rules about who can live where. Usually, these rules are set by "magic numbers" (like 2, 8, 20) that represent full, happy neighborhoods. When a shell is full, the nucleus is very stable.

However, in the world of "exotic" atoms (those with too many neutrons), the city layout changes. Some neighborhoods disappear, and new ones appear.

This paper is about investigating a specific neighborhood in a tiny, unstable city called Carbon-17 (17C). Scientists want to know if there is a massive "wall" (a shell gap) separating two specific neighborhoods: the 1d3/2 and the 2s1/2. If this wall is high enough, it means the city is very stable at a specific number of neutrons (N=16). If the wall is low, the city is chaotic.

Here is how the scientists figured this out, explained through a story:

1. The Setup: A Bumpy Ride

The nucleus of Carbon-17 is weird. It's made of a heavy, slightly squashed "core" (Carbon-16) and one lonely, weakly attached neutron. Because the core is squashed (deformed), it's like a city built on a hillside rather than flat ground.

To study this, the scientists performed a "transfer reaction." Imagine throwing a tennis ball (a deuteron, which is a proton and neutron stuck together) at the Carbon-16 city. The goal is to knock the proton out and leave the neutron behind, creating Carbon-17.

2. The Problem: The "Ghost" Neighbors

In the past, scientists could only study the "citizens" (neutrons) that stayed inside the city (bound states). But in Carbon-17, some neutrons are so energetic they don't stay; they immediately run away. These are the unbound states (the continuum).

Studying these "ghosts" that run away is incredibly hard. It's like trying to map a city by only looking at the people who are currently running out the front door. You can't just take a photo; you have to predict where they would have gone if they had stayed.

3. The Solution: The "Pixelated" Map

To solve this, the authors used a clever trick called Pseudo-States.

Imagine you are trying to draw a smooth, curved line on a computer screen, but your screen only has square pixels. You can't draw a perfect curve, but if you use enough tiny pixels, the curve looks smooth to the human eye.

The scientists did the same thing. They broke the infinite, smooth "runaway" energy of the neutron into a finite number of "pixels" (pseudo-states). This allowed them to use their computer models to simulate what happens when the neutron runs away, treating the "runaway" energy as a series of discrete steps.

They also had to deal with the Pauli Exclusion Principle. Think of this as a strict bouncer at a club. The core of the city is already full of neutrons. The new neutron wants to move in, but it can't take a spot that's already occupied. The scientists had to program their model to respect this "bouncer," ensuring the new neutron didn't try to live in a house that was already taken.

4. The Experiment: Testing the Wall Height

The team ran their simulation with different assumptions about the height of the "wall" (the shell gap) between the neighborhoods.

• Scenario A: They assumed the wall was short (small gap).
• Scenario B: They assumed the wall was tall (large gap).

They then compared their computer-generated "traffic patterns" (where the neutrons went and how much energy they had) against the actual data collected from a real experiment at a facility called GANIL in France.

5. The Verdict: The Wall Must Be Tall

The results were clear:

• When they assumed the wall was short (less than 5 MeV), the simulated "traffic" didn't match the real data. The neutrons ran away at the wrong speeds and in the wrong directions.
• When they assumed the wall was tall (greater than 5 MeV), the simulation matched the real-world data perfectly.

The Big Picture

This paper confirms that in the exotic world of Carbon-17, there is indeed a massive, high wall separating the neutron neighborhoods. This "N=16 shell gap" is a sign of stability in a place where we didn't expect it.

In simple terms:
The scientists built a sophisticated video game engine to simulate a tiny, unstable atom. They realized that to make the game look like reality, they had to assume there was a huge, unbreakable wall in the atom's structure. This proves that even in the wildest, most unstable parts of the atomic world, nature still likes to keep things organized and stable.

Technical Summary

1. Problem Statement

The study addresses the evolution of nuclear shell structure in exotic, weakly bound nuclei far from the stability line. Specifically, it focuses on the $N=16$ shell closure in neutron-rich carbon isotopes ($^{16,17}\text{C}$).

• Context: While the $N=16$ gap is well-established in Oxygen isotopes, its persistence in Carbon isotopes requires experimental verification.
• Challenge: Recent experimental data from GANIL on the $^{16}\text{C}(d, p)^{17}\text{C}$ reaction populated unbound (continuum) states of $^{17}\text{C}$. However, theoretical interpretation is difficult because:
  1. $^{17}\text{C}$ consists of a weakly bound neutron and a deformed core ($^{16}\text{C}$), requiring models that account for deformation.
  2. Standard reaction models often struggle with the continuum (unbound states) and Pauli-blocking effects (the exclusion of single-particle states already occupied by core nucleons).
  3. Previous analyses were limited to bound states; extending this to the continuum is necessary to accurately locate the $1d_{3/2}$ orbital and quantify the shell gap ($\Delta$).

2. Methodology

The authors developed a theoretical framework combining a semi-microscopic structure model with a transfer-to-the-continuum reaction formalism.

A. Structure Model: NAMD

• Model: The Nilsson+AMD (NAMD) model is used. It treats $^{17}\text{C}$ as a two-body system: a valence neutron and a deformed $^{16}\text{C}$ core.
• Hamiltonian: Includes kinetic energy, a spin-orbit term, and a valence-core interaction derived by convoluting effective nucleon-nucleon interactions with microscopic core densities calculated via Antisymmetrized Molecular Dynamics (AMD).
• Pauli Blocking: Two methods are tested to handle the Pauli exclusion principle:
  1. No-Blocking (NoB): States dominated by occupied Nilsson orbitals are simply discarded.
  2. Total Blocking (TB): The $N$ lowest energy Nilsson levels are explicitly excluded from the calculation.
• Basis: The Hamiltonian is diagonalized using a Pseudo-States (THO) basis (Transformed Harmonic Oscillator). This discretizes the continuum into a set of discrete states, allowing for the calculation of resonances and continuum properties.

B. Reaction Formalism: ADWA + Pseudo-States

• Framework: The Adiabatic Distorted Wave Approximation (ADWA) is employed for the $^{16}\text{C}(d, p)^{17}\text{C}$ reaction.
• Continuum Treatment:
  • The standard ADWA is adapted for unbound states using the pseudo-states discretization method.
  • Scattering Amplitudes: Calculated using the FRESCO code in the prior form of the transition potential ($\Delta U = V_{nC} + U_{pC} - U_{dC}$). The prior form is chosen to avoid convergence issues associated with the slow decay of unbound state wavefunctions in the post form.
  • Vertex Functions: The scattering amplitude relies on vertex functions extracted from the NAMD structure model.
  • Reconstruction: The discrete scattering amplitudes are converted back into a continuous energy distribution using the energy density of pseudo-states ($G_{\alpha i}^J$), which represents the overlap between pseudo-states and exact continuum wavefunctions.

3. Key Contributions

1. Extension to Continuum: The work successfully extends the NAMD model, previously validated for bound states, to transfer-to-the-continuum reactions.
2. Methodological Integration: It demonstrates a practical tool for combining deformed two-body structure models with ADWA reaction theory, specifically handling the discretization of the continuum via pseudo-states.
3. Pauli Blocking Validation: It provides a comparative analysis of NoB vs. TB methods in the context of continuum transfer, highlighting the necessity of proper Pauli blocking for bound states.
4. Constraint on Shell Gap: It uses the reaction data to place a quantitative lower bound on the $N=16$ shell gap in $^{17}\text{C}$.

4. Results

• Comparison with GANIL Data:
  • The calculated energy differential cross sections ($d\sigma/dE$) were compared with experimental data at 17.2 MeV/nucleon.
  • Bound States: The TB method significantly outperforms the NoB method in reproducing the cross-section for bound states, confirming the importance of Pauli blocking.
  • Continuum States: Both methods yield similar results for the continuum, but the TB method is preferred for overall consistency.
• Shell Gap Analysis:
  • The population of $3/2^+$ states (associated with the $1d_{3/2}$ orbital) is the dominant feature in the continuum region ($E > 2$ MeV).
  • The authors tested variations of the NAMD model by reducing the spin-orbit splitting to lower the energy of the $1d_{3/2}$ orbital, effectively reducing the shell gap ($\Delta$).
  • Findings:
    • When $\Delta$ was reduced to 5.0 MeV and 4.0 MeV, the calculated cross-section peak shifted to lower energies (around 2 MeV) and became too concentrated, disagreeing with the experimental data.
    • The experimental data shows a distribution concentrated between 3 and 6 MeV.
  • Conclusion: To reproduce the data, the $1d_{3/2}$ orbital must lie more than 3 MeV above the neutron separation threshold. This implies a shell gap $\Delta > 5$ MeV.
  • The standard NAMD (TB) model predicts $\Delta = 5.8$ MeV, which is consistent with the experimental lower limit of $5.08^{+0.43}_{-0.33}$ MeV reported in previous studies.

5. Significance

• Validation of Shell Evolution: The results strongly support the existence of a robust $N=16$ shell gap in neutron-rich carbon isotopes, confirming that the magic number persists even in this region of the nuclear chart.
• Theoretical Tool: The developed formalism provides a robust framework for studying other exotic nuclei where deformation and continuum effects are critical (e.g., $^{18}\text{C}$, $^{30}\text{Ne}$).
• Future Applications: The method is applicable to other transfer reactions populating unbound states, aiding in the mapping of "islands of inversion" and the evolution of shell structure in unstable nuclei.

In summary, the paper successfully combines advanced nuclear structure modeling with reaction theory to resolve a specific nuclear physics question, demonstrating that the $N=16$ shell gap in $^{17}\text{C}$ is large ($>5$ MeV) and that proper treatment of Pauli blocking and the continuum is essential for accurate theoretical predictions.