Inelastic nucleon-nucleus scattering from a microscopic point of view

This paper presents a parameter-free microscopic multiple scattering model based on the distorted-wave approximation and ab initio nuclear densities to successfully describe inelastic proton scattering off $^{12}$C across a range of energies, demonstrating the reliability of extending elastic scattering formalisms to inelastic transitions.

The Gist

Imagine you are trying to understand how a billiard ball (a proton) bounces off a complex, vibrating cluster of other balls (an atomic nucleus). Sometimes, the collision is simple: the ball bounces off, and the cluster stays the same. But often, the collision is "inelastic": the ball hits the cluster, and the cluster starts to wobble, spin, or vibrate in a new way, absorbing some of the energy.

This paper is about building a perfectly accurate map to predict exactly how that wobble happens, without using any "guesswork" or "fudge factors."

Here is the breakdown of their work using simple analogies:

1. The Goal: No More "Magic Numbers"

In the past, physicists trying to predict these collisions used "phenomenological" models. Think of this like a chef who knows a soup tastes good but doesn't know the recipe. They just keep adding salt and pepper (adjusting parameters) until the taste matches what they remember from a previous meal. It works for that specific meal, but if you change the ingredients, the chef has no idea what to do.

The authors of this paper wanted to be chefs who know the chemistry. They wanted to build a model from the ground up using only the fundamental laws of physics (the "ingredients") so that they could predict the result for any situation without ever having to tweak the recipe to fit the data.

2. The Method: The "Distorted Wave" Lens

To understand the collision, you can't just assume the proton flies in a straight line like a laser beam. The nucleus is heavy and has a strong gravitational-like pull (the nuclear force) that bends the proton's path.

• The Old Way (Plane Wave): Imagine throwing a ball in a vacuum. It goes straight. This is too simple for a nucleus.
• The New Way (Distorted Wave): Imagine throwing that ball through a thick, swirling fog. The fog bends the ball's path before it even hits the target. The authors use a mathematical tool called the Distorted-Wave Impulse Approximation (DWIA).
  • The Analogy: Think of the "fog" as the Optical Potential. It's a map of how the nucleus distorts the space around it. The authors calculated this map using only fundamental quantum mechanics, not by guessing.

3. The Three Ingredients (Potentials)

To make their prediction, they needed three specific "maps" (potentials):

1. The Incoming Map: How the proton's path is bent before it hits the nucleus (when the nucleus is calm).
2. The Outgoing Map: How the proton's path is bent after it hits the nucleus (when the nucleus is wobbling/excited).
3. The Transition Map: The specific "push" that causes the nucleus to go from calm to wobbling.

The Secret Sauce: Usually, scientists calculate these maps separately and hope they fit together. The authors did something clever: they realized all three maps are made of the exact same "dough," just with different "fillings."

• The "dough" is the fundamental force between two protons/neutrons (calculated using Chiral Effective Field Theory, a high-level quantum recipe).
• The "filling" is the density of the nucleus.
  • For the incoming map, they used the density of the calm nucleus.
  • For the outgoing map, they used the density of the wobbling nucleus.
  • For the transition map, they used the "difference" between the two.

By using the same high-quality "dough" for all three, they ensured the whole model was consistent.

4. The Test: The Carbon-12 Target

They tested their theory on a specific scenario: shooting protons at a Carbon-12 nucleus to make it vibrate into a specific excited state (the 2+ state at 4.44 MeV). They did this for protons moving at various speeds (energies between 65 and 300 MeV).

The Results:

• They compared their predictions to real-world experimental data.
• The Outcome: Their model, which had zero adjustable parameters (no "salt and pepper" added to fix the taste), matched the real data incredibly well, especially at higher speeds.
• They could predict not just how many protons bounced off, but the exact pattern of where they landed (the diffraction pattern), just like predicting the ripples in a pond after a stone is thrown.

5. Why This Matters

This paper is a big step forward because it proves we can understand complex nuclear reactions using only the fundamental rules of the universe, without needing to cheat by fitting numbers to past experiments.

• The Metaphor: It's like finally figuring out the laws of aerodynamics so well that you can design a plane that flies perfectly without ever having to test a model in a wind tunnel first.
• The Future: While their model works great for high speeds, it struggles a bit at very low speeds (where the "fog" gets too thick and complex). The authors admit this and suggest that future work will refine the model to handle those tricky low-speed collisions, eventually allowing us to predict reactions for all kinds of nuclei, not just Carbon.

In a nutshell: They built a "first-principles" simulator for nuclear collisions that is so accurate it matches reality without needing any manual tuning, proving that our understanding of the fundamental forces of nature is finally strong enough to predict complex atomic dances.

Technical Summary

1. Problem Statement

The paper addresses the challenge of describing inelastic nucleon-nucleus (NA) scattering from a fully microscopic, ab initio perspective. Traditional approaches to inelastic scattering often rely on phenomenological models that require empirical adjustments and fitting parameters (e.g., scaling transition densities to match experimental $B(E2)$ values), which limits their predictive power and universality.

The specific goal is to develop a theoretical framework that:

• Derives the reaction mechanism directly from fundamental nucleon-nucleon (NN) interactions.
• Avoids free adjustable parameters.
• Accurately describes the transition from the nuclear ground state to an excited state (specifically the $2^+$ state of $^{12}\text{C}$) across a wide range of projectile energies (65–300 MeV).

2. Methodology

The authors employ a Distorted-Wave Impulse Approximation (DWIA) framework grounded in Watson multiple scattering theory. The methodology integrates several advanced ab initio components:

A. Theoretical Framework

• Distorted-Wave Approximation: Unlike the Born approximation (which assumes plane waves), this approach accounts for the distortion of incoming and outgoing nucleon trajectories due to the nuclear potential.
• Three Potentials: The formalism requires the calculation of three distinct potentials:
  1. Initial Distorting Potential ($U_{el}$): Describes the interaction in the ground state (elastic optical potential).
  2. Final Distorting Potential ($U_{ex}$): Describes the interaction in the excited state.
  3. Transition Potential ($U_{tr}$): Drives the transition from the ground state to the excited state.
• Impulse Approximation: The interaction is approximated as a sum of single-nucleon interactions, where the projectile scatters off individual nucleons in the target using a free NN $t$-matrix.

B. Microscopic Inputs

• Nuclear Interactions: The NN interaction is derived from Chiral Effective Field Theory (Chiral EFT) at fifth order (N4LO) with a 500 MeV cutoff. Three-nucleon (3N) forces are included at N2LO.
• Nuclear Structure: The target nucleus ($^{12}\text{C}$) is described using the No-Core Shell Model (NCSM).
  • Wavefunctions are expanded in a harmonic oscillator basis ($N_{max}=8$, $\hbar\omega=16$ MeV).
  • Interactions are softened using the Similarity Renormalization Group (SRG) with a cutoff $\lambda_{SRG} = 2.0 \text{ fm}^{-1}$.
• Density Folding: The potentials are obtained by folding the nonlocal ab initio one-body nuclear densities (ground state, excited state, and transition density) with the NN $t$-matrix.
  • Consistency: The same Chiral EFT interaction is used to generate both the nuclear densities and the NN $t$-matrix, ensuring theoretical consistency.

C. Computational Algorithm

1. Calculate the optical potentials ($U_{el}, U_{ex}, U_{tr}$) by folding the NCSM densities with the Chiral EFT $t$-matrix.
2. Solve the Schrödinger equation with these nonlocal potentials to generate the initial ($|\psi^{(+)}\rangle$) and final ($\langle\psi^{(-)}_{*}|$) distorted waves.
3. Compute the inelastic transition amplitude $T_{fi}^{inel}$ using the matrix element $\langle\psi^{(-)}_{*}| U_{tr} |\psi^{(+)}\rangle$.
4. Derive the differential cross-sections without any empirical normalization.

3. Key Contributions

• Fully Microscopic Formalism: This is one of the first applications of a fully coherent, parameter-free microscopic multiple scattering approach to inelastic NA scattering. It extends previous work on elastic scattering to inelastic processes.
• Unified Treatment of Potentials: The authors demonstrate that the initial, final, and transition potentials share the same formal expression, differing only in the specific nuclear density matrix used (ground, excited, or transition density).
• Consistency: By using the same Chiral EFT interaction for structure (densities) and reaction ($t$-matrix) calculations, the model eliminates ambiguities often introduced by mixing different interaction models.
• No Free Parameters: The model contains no adjustable parameters or phenomenological scaling factors (e.g., no scaling of transition densities to match $B(E2)$ values).

4. Results

The model was benchmarked against experimental data for inelastic proton scattering off $^{12}\text{C}$, specifically the transition to the first $2^+$ excited state at 4.44 MeV.

• Energy Range: Calculations were performed for projectile energies between 65 MeV and 300 MeV.
• High Energy Agreement (100–300 MeV): The model shows excellent agreement with experimental differential cross-sections. It accurately reproduces:
  • The absolute magnitude of the cross-sections.
  • The angular distributions, including the positions and depths of diffraction minima and maxima.
  • The results are comparable to or better than other sophisticated models (e.g., coupled-channel calculations using phenomenological interactions or g-folding models).
• Low Energy Discrepancy (65 MeV): At lower energies, the model slightly underestimates the magnitude of the cross-section. The authors attribute this to the breakdown of the Impulse Approximation at lower energies, where multi-step processes and coupled-channel effects become more significant.
• Angular Distribution: The model captures the diffraction patterns well but shows some overestimation at very small and very large scattering angles, suggesting that neglected terms (like tensor forces) or higher-order effects might be necessary for a complete description in those regions.

5. Significance

• Predictive Power: The study validates the predictive capability of ab initio nuclear theory in reaction physics. The ability to reproduce experimental data without fitting parameters demonstrates that the underlying Chiral EFT interactions and NCSM structure are sufficient to describe complex inelastic processes.
• Robustness: The consistency of results across a wide energy range (up to 300 MeV) highlights the robustness of the distorted-wave approach when applied with high-order Chiral EFT interactions.
• Future Directions: The work establishes a foundation for extending ab initio methods to more complex reactions, such as transfer and knockout processes, and to nuclei further from stability. It also identifies the need for future refinements, such as including coupled-channel effects and tensor terms in the potential, to improve accuracy at lower energies and extreme scattering angles.

In conclusion, the paper successfully bridges the gap between nuclear structure and reaction theory, providing a rigorous, parameter-free framework for understanding inelastic nucleon-nucleus scattering.