Description of nucleon elastic scattering off $^6$Li with the four-body continuum-discretized coupled-channels method

This paper presents a semi-microscopic four-body continuum-discretized coupled-channels (CDCC) model using the JLM effective interaction that successfully describes neutron and proton elastic scattering off $^6$Li, along with relevant cross sections, across the 7 to 50 MeV energy range.

The Gist

The Big Picture: Why Do We Care?

Imagine a future where we have fusion power plants (like the sun, but on Earth) to generate clean energy. To make these plants work, we need a special "blanket" made of Lithium to catch fast-moving neutrons and turn them into fuel.

However, to build a safe and efficient blanket, we need to know exactly how neutrons behave when they hit Lithium atoms. Specifically, we need to know what happens when a neutron smashes into a Lithium-6 atom at high speeds (energies up to 50 MeV).

The problem? The current "instruction manuals" (nuclear databases) for this are a bit fuzzy or missing pieces for these high speeds. This paper is an attempt to write a better, more accurate manual.

The Problem with Old Maps

Scientists have tried to map these reactions before using two main approaches:

1. The "Ab Initio" Approach: Trying to calculate every single tiny interaction from scratch. It's like trying to predict the weather by tracking every single air molecule. It's incredibly accurate but too slow and complex for high-energy collisions.
2. The "Phenomenological" Approach: Using simple rules of thumb based on past experiments. It's fast, but it often breaks down when you go outside the range of data you've already seen.

This paper tries to find a "Goldilocks" solution: a Semi-Microscopic Model. It's detailed enough to be accurate but simplified enough to be computable.

The New Tool: The "Four-Body CDCC"

To understand the Lithium-6 atom, the authors treat it not as a solid ball, but as a fragile family of particles.

• The Old Way: They used to think of Lithium-6 as just an Alpha particle (a tight cluster of 2 protons and 2 neutrons) holding hands with a Deuteron (a proton and a neutron).
• The New Way (This Paper): They realized that at high speeds, that "Deuteron" hand-holding can break. So, they modeled Lithium-6 as a three-body family: Alpha + Proton + Neutron.

When a neutron comes in to hit this family, the whole system becomes a four-body dance (The incoming neutron + Alpha + Proton + Neutron).

They use a method called CDCC (Continuum-Discretized Coupled-Channels).

• The Analogy: Imagine a trampoline. If you jump on it gently, it bounces predictably. But if you jump hard, the trampoline fabric stretches, ripples, and might even tear (breakup).
• The CDCC method allows the scientists to calculate not just the bounce (elastic scattering), but also the ripples and tears (breakup channels) that happen when the neutron hits the Lithium.

The "Recipe" (The JLM Interaction)

To make the math work, they use a standard "recipe" for how particles interact, called the JLM interaction (named after Jeukenne, Lejeune, and Mahaux).

However, recipes aren't always perfect. Sometimes you need to adjust the amount of salt or sugar. In this paper, the "salt" and "sugar" are called Renormalization Factors.

• Real Part (The "Salt"): They found that the "salt" amount needed to be constant. They just needed to add a little bit more (1.1 times the standard amount) to get the physics right.
• Imaginary Part (The "Sugar"): This part is trickier. It represents how much energy is "lost" or absorbed during the collision. They found that the amount of "sugar" needed changes smoothly depending on how fast the neutron is moving. They created a simple formula to calculate exactly how much "sugar" to add at any speed.

The Results: A Better Map

The authors ran their new four-body simulation and compared it to real-world experimental data.

1. The "Breakup" Matters: When they ignored the possibility of the Lithium breaking apart (the old way), their predictions were wrong. When they included the breakup (the new way), the predictions matched the real data almost perfectly.
   • Analogy: It's like trying to predict how a car crash looks. If you assume the car stays in one piece, your prediction is wrong. You have to account for the crumpling and parts flying off to get the real picture.
2. The Sweet Spot: Their new model works beautifully for neutron energies between 7 MeV and 50 MeV.
   • Below 7 MeV: The model gets a bit wobbly (needs a more complex approach).
   • Above 50 MeV: The model starts to underestimate the reaction (needs a different high-energy approach).
3. Protons too: They tested the model with protons (positively charged particles) as well, and it worked just as well as it did for neutrons.

The Conclusion

This paper successfully built a reliable "instruction manual" for how neutrons and protons bounce off Lithium-6 at high speeds.

• Why it matters: This helps engineers design better fusion reactors (like the IFMIF facility mentioned) by giving them accurate data on how neutrons will behave.
• What's next: The authors say, "We've mastered the bounce and the total energy loss. Next, we want to use some fancy math tricks (Complex Scaling) to describe exactly how the Lithium breaks apart and what pieces fly off."

In a nutshell: They took a complex atomic collision, broke it down into a four-part family dance, adjusted their mathematical recipe slightly, and found a method that predicts the outcome with high accuracy for the energy range needed for future fusion power.

Technical Summary

1. Problem Statement and Motivation

• Context: Lithium isotopes (specifically ⁶Li and ⁷Li) are critical components in fusion reactor "blankets" (e.g., the International Fusion Material Irradiation Facility, IFMIF) for decelerating neutrons and breeding tritium.
• The Gap: While nuclear databases (like FENDL) cover data up to 150 MeV, there is a lack of reliable high-energy neutron data for lithium in the 7 MeV to 50 MeV range. This energy range is crucial for understanding particle production in fusion environments.
• Limitations of Existing Models:
  • Ab initio approaches are computationally prohibitive at energies above 20 MeV due to the vast number of reaction channels.
  • Previous semi-microscopic models (e.g., Matsumoto et al., 2011) using a three-body α + d model for ⁶Li and CDCC (Continuum-Discretized Coupled-Channels) successfully described data up to ~25 MeV but failed to reproduce total cross-sections above 25 MeV.
  • Previous models often neglected "closed channels" (channels with negative relative energy) and used a simplified two-body structure for ⁶Li.

2. Methodology

The authors developed a semi-microscopic reaction model based on the Four-Body Continuum-Discretized Coupled-Channels (CDCC) method.

• Structure of ⁶Li: Instead of the previous α + d (two-body) model, the authors adopted an α + p + n three-body model. This allows for a more accurate description of the ⁶Li breakup continuum and includes the α + p + n configuration explicitly.
• Formalism:
  • The total wave function is expanded in terms of eigenstates of the ⁶Li Hamiltonian, calculated using the Gaussian Expansion Method (GEM).
  • The calculation includes both open channels (positive energy, leading to scattering) and closed channels (negative energy, bound or virtual states), the latter of which were previously omitted in similar studies.
  • Coupling potentials are generated via a single-folding model using the JLM (Jeukenne-Lejeune-Mahaux) g-matrix effective interaction.
• Interaction Parameters:
  • The JLM interaction is renormalized by two factors: $N_V$ (real part) and $N_W$ (imaginary part).
  • These factors were treated as free parameters to be fitted against experimental data, rather than fixed to standard optical model values.
• Energy Range: Calculations were performed for nucleon (neutron and proton) scattering off ⁶Li from 4 MeV to 72 MeV, with a focus on the 7–50 MeV range.

3. Key Contributions

1. Four-Body CDCC Implementation: Successfully applied a four-body CDCC framework (incident nucleon + α + p + n) to nucleon-⁶Li scattering, significantly expanding the model space compared to previous three-body (α + d) approaches.
2. Inclusion of Closed Channels: Explicitly included closed channels in the coupled-channels equations. The study demonstrates that these channels are essential for accurate descriptions at lower energies (below ~24 MeV).
3. Systematic Parameter Determination: Determined the energy dependence of the JLM renormalization factors ($N_V$ and $N_W$) up to 50 MeV, moving beyond the assumption of constant parameters used in earlier works.
4. Validation of Applicability: Established a reliable energy window (7–50 MeV) where this semi-microscopic approach is valid, bridging the gap between low-energy ab initio methods and high-energy multiple scattering theories.

4. Results

• Renormalization Factors:
  • Real Part ($N_V$): Found to be a constant value of 1.1 across the energy range.
  • Imaginary Part ($N_W$): Found to have a smooth energy dependence, fitted by the function $N_W = 0.213 \ln(E) - 0.247$.
• Elastic Scattering:
  • The four-body CDCC model with the fitted parameters successfully reproduced the angular distributions for both neutron and proton elastic scattering across 14 energies (up to 50 MeV).
  • Breakup Importance: Calculations without breakup channels (single-channel) failed to reproduce experimental data, confirming the necessity of the CDCC approach.
  • Closed Channel Importance: At low energies (<24 MeV), calculations excluding closed channels showed significant deviations, proving their role in the reaction dynamics.
• Total Cross Sections:
  • The model accurately reproduced the neutron total cross section ($\sigma_{tot}$) and total reaction cross section ($\sigma_R$) for energies above 7 MeV.
  • Below 7 MeV, the model undershoots the data, indicating the need for more sophisticated many-body approaches at very low energies.
  • Above 50 MeV, the model tends to undershoot the total cross section, suggesting limitations of the JLM interaction at higher energies.
• Proton Scattering: The model, calibrated on neutron data, successfully predicted proton elastic scattering and reaction cross sections up to 50 MeV and even at 72 MeV (though the latter is an extrapolation).

5. Significance and Future Outlook

• Nuclear Data Science: The constructed model provides a reliable theoretical basis for generating neutron reaction data for ⁶Li in the 7–50 MeV range, which is vital for the design and safety analysis of fusion reactors like IFMIF.
• Methodological Advancement: The study validates that semi-microscopic CDCC with a three-body target structure and closed-channel inclusion is a robust alternative to fully microscopic approaches in the intermediate energy regime.
• Future Work:
  • The authors note that while the scattering matrices for discretized breakup channels are calculated, generating continuous observables requires complex scaling methods and channel selection techniques (e.g., CSLS) to disentangle $\alpha+d$ and $\alpha+p+n$ components.
  • Future papers will apply these techniques to describe inelastic scattering and breakup reactions of ⁶Li.
  • Plans are underway to extend this model to ⁷Li using an $\alpha + p + n + n$ four-body model.

In conclusion, this paper presents a significant step forward in modeling nucleon-lithium interactions, offering a validated, semi-microscopic tool that fills a critical gap in nuclear data for fusion applications.