Microscopic optical potential framework applied to neutron scattering on deformed $^{48,50}$Cr

This paper presents a unified microscopic framework based on a symmetry-restored multi-excitation generator coordinate method to derive non-local optical potentials from an effective Hamiltonian, successfully applying it to calculate neutron scattering cross sections for deformed chromium isotopes while consistently treating nuclear structure and reaction observables.

The Gist

Imagine you are trying to predict how a billiard ball (a neutron) will bounce off a complex, squishy, spinning target (an atomic nucleus).

For decades, physicists have used a "best guess" map to predict these bounces. This map, called an Optical Potential, is like a weather forecast: it's built by looking at past storms (experimental data) and adjusting the model until it fits. It works great for places we've visited before, but if you try to use it for a new, strange, or radioactive nucleus (a place we've never been), the forecast often fails.

This paper presents a new way to build that map. Instead of guessing based on past weather, the authors build the map from the ground up, using the fundamental laws of physics that govern the nucleus itself.

Here is the story of how they did it, broken down into simple concepts:

1. The Problem: The "Squishy" Target

Atomic nuclei aren't perfect, hard spheres. They are more like spinning, squishy jellybeans. Some are round, but many (like the Chromium isotopes studied here) are deformed—they look like rugby balls or even pretzels.

When a neutron hits this jellybean, it doesn't just bounce off the surface. It can:

• Make the jellybean wobble (vibration).
• Spin it faster (rotation).
• Get temporarily stuck inside before popping out.

Traditional models often treat the nucleus as a static, smooth ball. This misses all the "wiggles" and internal chaos, leading to inaccurate predictions.

2. The Solution: The "Symmetry-Restored" Blueprint

The authors used a sophisticated mathematical tool called the Generator Coordinate Method (GCM). Think of this as a way to take a 3D scan of the nucleus from every possible angle and shape.

• Breaking Symmetry: First, they let the nucleus be messy. They allowed it to be deformed, spinning, and vibrating, ignoring the fact that nature usually prefers things to look the same from every angle (symmetry).
• Restoring Symmetry: Then, they used a mathematical "projector" to force the messy, squishy shapes back into a clean, symmetrical state. This is like taking a crumpled piece of paper and ironing it out, but keeping the memory of the wrinkles in the final result.

This process creates a Microscopic Optical Potential. Instead of a smooth, featureless wall, the potential is a detailed, bumpy landscape that reflects the actual internal structure of the nucleus.

3. The Missing Pieces: The "Fill-in-the-Blanks" Trick

Here is the tricky part: To build this perfect map, you need to know every possible state the nucleus can be in. But even with supercomputers, you can only calculate a fraction of them. It's like trying to draw a map of a city but only being able to visit 10% of the streets.

If you stop there, your map is incomplete, and the billiard ball (neutron) will get lost.

The authors used a clever trick called Sum Rules.

• Imagine you have a bucket of sand (the total strength of the nucleus). You've scooped out 80% of the sand into a jar (the states you calculated).
• You know exactly how much sand should be in the bucket based on physics laws.
• The "Sum Rule" tells you exactly how much sand is missing (the 20%).
• Instead of trying to find every single grain of missing sand, they created a "Mean Field"—a smooth, average representation of that missing 20%.

This allows them to fill in the gaps of their map without needing to calculate every single impossible state, ensuring the neutron has a complete path to follow.

4. The Result: A Better Map for the Future

They tested this new method on Chromium-48 and Chromium-50. These are important because Chromium is used in stainless steel for nuclear reactors. Knowing exactly how neutrons bounce off them is crucial for reactor safety and efficiency.

• The Old Way: Used a smooth, averaged map. It missed the peaks and valleys of the actual data.
• The New Way: Used the detailed, "squishy" map. It matched the real-world experimental data almost perfectly, especially in the energy range where nuclear reactors operate (1 to 15 MeV).

Why Does This Matter?

Think of this as moving from drawing a map by hand (phenomenological models) to using a satellite GPS (microscopic models).

1. Predicting the Unknown: Because this method is built from fundamental laws, not just past data, it can predict how neutrons will behave on nuclei we haven't even discovered yet, or on exotic stars in the universe.
2. Safer Reactors: Better predictions mean we can design nuclear reactors that are safer and more efficient.
3. Unified Theory: It finally connects the "structure" of the nucleus (what it looks like inside) with the "reaction" (how it behaves when hit). Before, these were two separate worlds; now, they are one.

In a nutshell: The authors built a high-definition, 3D map of how neutrons interact with squishy, spinning nuclei. They did it by calculating the nucleus's internal "wiggles," using a smart trick to fill in the missing data, and proving that this physics-first approach works better than the old "guess-and-check" methods.

Technical Summary

1. Problem Statement

Nuclear reaction modeling, particularly for neutron scattering, relies heavily on the optical potential to describe particle scattering while accounting for all other reaction channels.

• Limitations of Current Methods: The most widely used optical potentials are phenomenological, meaning they are fitted to experimental data. While effective for stable nuclei near the fitting region, they lack predictive power for exotic nuclei, radioactive beams, and astrophysical environments where experimental data is scarce or non-existent.
• The Challenge: Deriving optical potentials directly from microscopic nuclear structure (wavefunctions and Hamiltonians) is theoretically exact but computationally difficult. It requires a consistent treatment of nuclear structure (bound states) and reaction dynamics (continuum), especially for deformed nuclei where collective correlations are strong.
• Specific Context: Chromium isotopes (specifically $^{48}\text{Cr}$ and $^{50}\text{Cr}$) are critical components in stainless steel used in nuclear reactors. Uncertainties in their neutron cross-sections impact reactor safety and critical mass calculations, yet high-precision experimental data is limited.

2. Methodology

The authors propose and implement a unified framework that derives non-local optical potentials directly from a microscopic many-body solution. The core components are:

A. Symmetry-Restored Generator Coordinate Method (GCM)

• Effective Hamiltonian: The authors use an effective Hamiltonian ($\hat{H} = \hat{H}_0 + \hat{H}_P + \hat{H}_Q$) designed to reproduce Density Functional Theory (DFT) results. It includes a spherical single-particle part, a pairing interaction (seniority pairing), and a quadrupole-quadrupole interaction.
• Generator Coordinates: To capture collective behavior (vibrational, pairing, rotational), the GCM uses a 5-dimensional landscape of generator coordinates:
  • Triaxial quadrupole deformation ($\beta, \gamma$).
  • Cranking frequency (angular momentum around the x-axis).
  • Proton and neutron pairing strengths.
• Symmetry Restoration: The method explicitly breaks symmetries (particle number, angular momentum) in the intrinsic Hartree-Fock-Bogoliubov (HFB) states and then restores them via projection operators ($\hat{P}^I_{MK}, \hat{P}^N, \hat{P}^Z$).
• Basis Construction: The basis includes HFB vacua and excited states constructed using a Bogoliubov coupled-cluster (BCC) singles operator to capture particle-hole excitations.

B. Green's Function Formalism

• The optical potential $V$ is derived from the self-energy $\Sigma$ using the Dyson equation: $G = G_0 + G_0 \Sigma G$.
• The self-energy is calculated from the Green's function matrix elements, which depend on the spectroscopic amplitudes ($\sigma$) and energy levels of the $A \pm 1$ systems (target nucleus $\pm$ one neutron).
• Spectroscopic Amplitudes: Calculated efficiently between projected basis states using Pfaffians of overlap matrices, handling the non-orthogonality of the GCM basis.

C. Sum Rules and Completion

• Since the GCM basis is finite and incomplete, the authors employ sum rules (zeroth and energy-weighted) to estimate the contribution of missing states.
• Completion States: A set of "completion states" is constructed to satisfy the sum rules, effectively representing an average generalized mean-field for the missing strength. This ensures the Green's function has the correct asymptotic behavior and total strength.

D. Regularization and Scattering

• Regularization: To prevent unphysical oscillations (Gibbs phenomenon) arising from truncating the basis in momentum space, a smoothing regulator is applied to the potential matrix elements.
• Scattering Calculation: The resulting non-local, energy-dependent optical potential is used in the Lippmann-Schwinger equation to calculate phase shifts, differential cross-sections, and total reaction cross-sections.

3. Key Contributions

1. Unified Framework: The paper presents a fully microscopic framework that links nuclear structure (Hamiltonian solutions) directly to reaction observables without phenomenological tuning of the potential.
2. Deformed Nuclei Application: It successfully applies this method to deformed, triaxial nuclei ($^{48}\text{Cr}$ and $^{50}\text{Cr}$), a regime where standard spherical mean-field approaches often fail.
3. Sum Rule Completion: The rigorous implementation of sum rules to complete the basis allows the method to be applied even when the many-body solution is not exhaustive, estimating the contribution of high-energy states via a mean-field average.
4. Non-Local Potential Analysis: The study provides detailed insights into the non-locality of the optical potential, showing that absorption (imaginary part) is more significant at the nuclear surface, a feature consistent with other microscopic approaches but derived here from a different many-body method.

4. Results

• Cross-Section Accuracy:
  • For $^{48}\text{Cr}$ and $^{50}\text{Cr}$, the calculated total neutron scattering cross-sections show excellent agreement with evaluated nuclear data (ENDF/B-VIII.0 and JENDL-5) in the 2–10 MeV energy range.
  • The results significantly outperform calculations using only the spherical mean-field potential ($H_A$), particularly in smoothing out unphysical peaks associated with single-particle shells (e.g., $g_{7/2}$ and $h_{11/2}$).
• Differential Scattering:
  • Calculated differential elastic scattering cross-sections for $^{50}\text{Cr}$ at 5 energies match experimental data and evaluated libraries better than phenomenological models (like JLM) which tend to overestimate the elastic channel due to insufficient absorption.
• Convergence and Sensitivity:
  • The method is sensitive to the number of basis states ($n_{grid}$) and the treatment of low-energy resonances.
  • Low-energy convergence ($<2$ MeV) remains challenging due to the dominance of individual s-wave resonances, which require explicit coupling to the continuum (a future extension).
  • The imaginary part of the potential (controlled by the parameter $\eta$) has a minor impact on the total cross-section in the 2–10 MeV range but is crucial for damping unphysical resonances.

5. Significance and Outlook

• Predictive Power: This framework enables the prediction of neutron cross-sections for nuclei far from stability (exotic nuclei) where no experimental data exists, by relying solely on the underlying nuclear Hamiltonian.
• Nuclear Technology: By reducing uncertainties in Chromium cross-sections, this work directly supports the design and safety analysis of nuclear reactors and fusion devices.
• Theoretical Advancement: It demonstrates that reaction observables can serve as constraints on nuclear structure models, creating a feedback loop that refines our understanding of the nuclear force and many-body correlations.
• Future Work: The authors plan to extend the framework to include inelastic scattering (excited target states) and improve the treatment of the low-energy resonant region by explicitly coupling the many-body states to the continuum.

In summary, this work represents a significant step toward a parameter-free, unified description of nuclear structure and reactions, proving that microscopic many-body methods can yield high-precision optical potentials for complex, deformed nuclei.