Microscopic Optical Potential from Brueckner-Hartree-Fock Theory

This paper develops a microscopic optical potential for nucleon-nucleus scattering based on Brueckner-Hartree-Fock theory and the local density approximation, demonstrating its quantitative agreement with phenomenological models and experimental data for neutron and proton scattering on calcium isotopes below 200 MeV.

The Gist

Imagine the atomic nucleus not as a solid ball, but as a bustling, crowded city made of tiny particles called protons and neutrons. When a new particle (a "projectile") tries to crash a party in this city, it doesn't just bounce off a wall; it interacts with the whole crowd. Sometimes it bounces off cleanly (elastic scattering), and sometimes it gets absorbed or causes chaos (reaction).

To predict how this particle will behave, physicists use a tool called an Optical Potential. Think of this potential as a "force field map" that tells the incoming particle how the city will push or pull on it.

For decades, scientists have used Phenomenological Potentials. These are like weather maps created by looking at past storms and drawing lines where the wind usually blows. They work great for cities we know well (stable atoms), but if you try to use them to predict the weather in a brand-new, strange city (exotic, unstable atoms), the predictions often fail because the map was never drawn for those conditions.

This paper introduces a new kind of map: The Microscopic Optical Potential (MOP).

Instead of guessing the shape of the map based on past data, the authors built it from the ground up using the fundamental laws of physics. Here is how they did it, explained simply:

1. Building the Blueprint (The Nuclear Matter Theory)

First, the authors looked at an infinite, perfect "city" of nuclear matter (a theoretical place where protons and neutrons are spread out evenly). They used a sophisticated mathematical engine called Brueckner-Hartree-Fock (BHF) theory.

• The Analogy: Imagine trying to understand how a single person moves through a crowd. Instead of watching one specific party, you simulate the physics of how every person interacts with every other person in a perfect, endless crowd. You calculate exactly how much they push, pull, and slow each other down.
• They included a tricky detail called Three-Body Forces. Usually, we think of interactions as pairs (Person A pushes Person B). But in a dense crowd, Person A might push Person B, but Person C standing right next to them changes how that push feels. The authors accounted for these "group dynamics."

2. Turning the Blueprint into a Real Map (The Local Density Approximation)

You can't use a map of an infinite city to navigate a real, finite city with edges and walls. The authors had to translate their infinite "perfect city" calculations into a map for real atoms (like Calcium-40 and Calcium-48).

• The Analogy: They used a technique called the Improved Local Density Approximation (ILDA). Imagine taking a high-resolution photo of the infinite city and then "smearing" it slightly to account for the fact that real cities have fuzzy edges. This allowed them to take their perfect physics calculations and apply them to the specific, lumpy shapes of real atoms.

3. Testing the New Map

To see if their new map worked, they sent "virtual particles" (neutrons and protons) crashing into Calcium atoms in a computer simulation. They measured three things:

1. Where the particles bounced: (Differential cross-sections)
2. How they spun: (Analyzing powers)
3. How many got stuck: (Total/Reaction cross-sections)

The Results:

• The Good News: Their new, physics-based map (MOP) matched the real experimental data almost as well as the old, guess-based maps (Koning-Delaroche potentials).
• The Big Win: Because their map is built on fundamental laws rather than past data, it is reliable for exotic nuclei. If scientists discover a new, unstable atom in a lab, they can use this map to predict how it will react, even if no one has ever seen that atom before.

Why This Matters

Think of the old method as having a dictionary of words you've already heard. If someone speaks a new language or a new word, you're stuck. This new method is like understanding the grammar and roots of the language. You can now understand and predict words you've never heard before.

This is crucial for the future of nuclear physics. As we build better accelerators to create "exotic" atoms (atoms that don't exist naturally on Earth), we need tools that can predict their behavior without needing to test them first. This paper provides a robust, "from-scratch" tool that could help us understand everything from how stars explode to how we might create new materials.

In a nutshell: The authors stopped guessing the rules of the nuclear game and started calculating them from first principles. They built a universal translator that works for both familiar atoms and the strange, exotic ones we are just beginning to discover.

Technical Summary

1. Problem Statement

The study of nuclear reactions, particularly with exotic nuclei far from the stability line, relies heavily on the Optical Model, which uses a complex, energy-dependent potential to describe projectile-target interactions.

• Limitation of Phenomenology: Traditional optical potentials (e.g., the global Koning-Delaroche or KD potential) are phenomenological. They rely on adjustable parameters fitted to data from stable nuclei. When applied to exotic nuclei (far from stability), these models suffer from uncontrolled uncertainties and lack a rigorous theoretical foundation.
• Need for Microscopy: There is a critical need for a Microscopic Optical Potential (MOP) derived directly from realistic nucleon-nucleon (NN) interactions. Such a potential would guarantee reliable predictions for both stable and exotic nuclei without free parameters.
• Gap in Current Theory: While Brueckner-Hartree-Fock (BHF) theory has been successful in describing infinite nuclear matter properties, extending these results to finite nuclei (where density varies and surface effects dominate) remains a challenge. Previous BHF-based MOPs often addressed only the volume part or lacked a consistent treatment for finite nuclei.

2. Methodology

The authors develop a microscopic optical potential for finite nuclei by bridging the gap between infinite nuclear matter calculations and finite nuclear systems.

• Theoretical Framework (BHF):

  • Calculations are performed using the Brueckner-Hartree-Fock (BHF) approach.
  • Interactions: The realistic Argonne V18 nucleon-nucleon interaction is used, augmented by a consistent microscopic Three-Body Force (TBF) approximated as an effective two-body force.
  • Self-Energy: The nucleon self-energy ($\Sigma$) is calculated, including the BHF mean field, core polarization, and the TBF rearrangement term.
  • Nuclear Matter: The self-energy is computed for symmetric and asymmetric nuclear matter across a range of densities ($\rho = 0.06 - 0.20 \text{ fm}^{-3}$) and energies ($10 - 200 \text{ MeV}$).

• Parameterization:

  • The BHF results for the self-energy are parameterized into analytic forms following the Jeukenne-Lejeune-Mahaux (JLM) ansatz.
  • The potential is split into isoscalar ($V_0, W_0$) and isovector ($V_1, W_1$) components, dependent on density ($\rho$), energy ($E$), and isospin asymmetry ($\beta$).
  • Analytic expressions are fitted to the BHF data points, ensuring the potential vanishes naturally at zero density.

• Construction for Finite Nuclei (ILDA):

  • To apply the nuclear matter potential to finite nuclei, the authors use an Improved Local Density Approximation (ILDA).
  • Finite-Range Correction: Unlike the standard LDA (a zero-range approximation), ILDA includes finite-range effects by convolving the local potential with a Gaussian form factor (width $t = 1.4 \text{ fm}$). This corrects for the underestimation of surface diffuseness in standard LDA.
  • Density Distribution: Nuclear densities for $^{40}\text{Ca}$ and $^{48}\text{Ca}$ are calculated using the Hartree-Fock (HF) approximation with the LNS5 Skyrme interaction. The LNS5 interaction was specifically constrained by BHF results to ensure consistency.
  • Spin-Orbit Potential: The spin-orbit term is derived from the HF calculation using the LNS5 interaction, as the spin-orbit contribution vanishes in infinite nuclear matter.

3. Key Contributions

1. First Consistent Finite-Nucleus MOP from BHF: The paper presents a complete MOP for finite nuclei derived directly from BHF theory with realistic NN and TBF interactions, moving beyond previous studies that were limited to nuclear matter or the volume part of the potential.
2. Analytic Parameterization: The authors provide explicit analytic formulas (with coefficients listed in tables) for the real and imaginary parts of the optical potential. This makes the MOP readily usable for analyzing experimental data on exotic nuclei without re-running complex many-body calculations.
3. Improved Local Density Approximation: The implementation of ILDA with a fixed Gaussian width ($t=1.4 \text{ fm}$) successfully addresses the surface diffuseness issues inherent in standard LDA, improving the agreement with phenomenological potentials.
4. Systematic Validation: The potential is rigorously tested against neutron and proton scattering on $^{40}\text{Ca}$ and $^{48}\text{Ca}$, covering elastic differential cross-sections, analyzing powers, and total/reaction cross-sections up to 200 MeV.

4. Results

• Potential Structure:

  • The central real and imaginary potentials derived from BHF show quantitative consistency with the global phenomenological KD potentials.
  • The isovector components (crucial for exotic nuclei) are explicitly calculated and parameterized.
  • The imaginary part exhibits a surface peak at low energies, though the microscopic potential shows slightly different behavior compared to KD at very low energies and for protons at high energies.

• Scattering Observables ($n/p + ^{40,48}\text{Ca}$):

  • Differential Cross Sections: The MOP reproduces experimental data well across the 10–200 MeV range. For $^{48}\text{Ca}$, the MOP outperforms the global KD potential, likely because KD was not fitted to $^{48}\text{Ca}$ data.
  • Analyzing Powers ($A_y$): The MOP correctly predicts the positions of peaks and valleys for $^{40}\text{Ca}$. For $^{48}\text{Ca}$, it captures the general trend, though some phase shifts are observed at higher energies.
  • Total/Reaction Cross Sections:
    • The MOP generally agrees with data but tends to underestimate reaction cross-sections in the 20–80 MeV range and overestimate them above 80 MeV.
    • The authors note that using a different Skyrme interaction (SKP) for the density distribution improves the total cross-section results, suggesting sensitivity to the density profile.
    • Replacing the microscopic imaginary part with the phenomenological KD imaginary part significantly improves the reaction cross-section description, indicating that the microscopic imaginary part (currently based on first-order mass operators) needs refinement (e.g., higher-order contributions).

5. Significance and Future Outlook

• Predictive Power for Exotic Nuclei: By grounding the optical potential in microscopic theory (BHF + TBF), this work provides a reliable tool for analyzing scattering data on rare isotopes where phenomenological fits are unreliable.
• Analytic Utility: The provision of analytic forms allows for immediate integration into reaction codes for global nuclear data evaluations.
• Future Improvements: The authors identify that the imaginary part of the potential requires higher-order many-body contributions to better describe reaction cross-sections at low energies. Future work will extend these calculations to a wider mass range ($56 \le A \le 208$) and explore different realistic interactions to quantify theoretical uncertainties. Additionally, the framework aims to extract information on the symmetry energy near saturation density from new scattering data.

In conclusion, this paper successfully bridges the gap between microscopic many-body theory and nuclear reaction phenomenology, offering a robust, parameter-free optical potential that performs comparably to global phenomenological models while offering superior predictive capabilities for exotic nuclear systems.