Comparison of Pauli projection and supersymetric transformation methods for three-body nuclear structure and reactions

This paper compares two methods for eliminating Pauli-forbidden states in three-body nuclear systems—Pauli projection and supersymmetric transformation—finding that while Pauli projection better matches experimental deuteron-${}^4$He scattering data, neither method demonstrates clear superiority for bound and resonant states, though systematic differences between them are observed.

The Gist

Imagine you are trying to build a house out of Lego bricks. In the world of nuclear physics, the "house" is an atom's nucleus, and the "bricks" are protons and neutrons.

This paper is about a specific problem that happens when you try to build a house with a core (a tight cluster of bricks already stuck together) and then try to add two extra bricks (nucleons) to the outside.

The Problem: The "No-Go" Zones

Inside the core, the bricks are already packed very tightly. Quantum physics (the rules of the very small) says that the outside bricks cannot just squeeze into the exact same spots the inside bricks are already occupying. If they did, it would be like trying to put two people in the same chair at the same time—it's forbidden!

In physics, we call these forbidden spots "Pauli-forbidden states." They are like deep, dark basements in the house that the new bricks are not allowed to enter.

The Two Solutions

The author, A. Deltuva, is testing two different ways to make sure the new bricks respect these "No-Go" zones while still building the house correctly. Think of these as two different construction strategies:

1. The "Bouncer" Method (Pauli Projection - PP)
Imagine hiring a strict bouncer at the door of the basement. If a brick tries to enter the forbidden zone, the bouncer kicks it out immediately with a massive, invisible force.

• How it works: The math adds a "repulsive wall" that pushes the forbidden states away to infinity.
• The vibe: It's a bit aggressive and non-local (the bouncer can reach across the room to kick you out).

2. The "Shape-Shifter" Method (Supersymmetric Transformation - SS)
Imagine you don't hire a bouncer. Instead, you magically reshape the basement itself. You turn the floor of the forbidden zone into a steep, slippery slide that leads nowhere, while keeping the rest of the house exactly the same.

• How it works: The math changes the potential energy landscape so the forbidden state disappears naturally, but the "allowed" states remain untouched.
• The vibe: It's a smoother, more elegant transformation that keeps the local geometry intact.

What Did They Find?

The author ran complex computer simulations (using "Faddeev equations," which are like super-advanced blueprints) to see which method builds a better house.

1. When things are flying around (Scattering):
They looked at what happens when a deuteron (a tiny two-brick cluster) smashes into a Helium-4 core.

• The Verdict: The Bouncer (PP) method won. The experimental data (what actually happens in real life) matched the Bouncer's predictions much better. The Shape-Shifter (SS) method and a third "Hard Core" method were too similar to each other and didn't match reality as well.
• Analogy: If you throw a ball at a wall, the Bouncer method predicts the bounce perfectly. The Shape-Shifter method predicts a bounce that's slightly off.

2. When things are stuck together (Bound States):
They looked at stable nuclei like Lithium-11 or Oxygen-18.

• The Verdict: It's a tie, but with a twist. Neither method is clearly "better" at predicting the exact weight (binding energy) of the nucleus.
• The Twist: The Shape-Shifter (SS) method tends to build tighter, heavier houses (more binding energy) than the Bouncer (PP) method.
• Why? The Bouncer method forces the bricks to be very "jittery" (high kinetic energy) to avoid the forbidden zones, which makes the house slightly looser. The Shape-Shifter method lets the bricks sit a bit more comfortably, making the house tighter.

3. The Resonance (The Unstable House):
They also looked at unstable nuclei (like 16-Be) that fall apart quickly.

• The Verdict: Again, the Shape-Shifter method predicted these unstable houses would fall apart at slightly lower energies (they are "tighter" before they break) compared to the Bouncer method.

The Big Picture

The paper concludes that while the Bouncer (PP) method is the champion for predicting how particles scatter and bounce off each other (matching real-world experiments), the Shape-Shifter (SS) method is a very close second for calculating how stable nuclei hold together.

However, they aren't identical. They produce systematic differences:

• Bouncer (PP): Higher energy, slightly looser binding, matches scattering data best.
• Shape-Shifter (SS): Lower energy, tighter binding, slightly different internal structure.

The Takeaway:
If you are a nuclear physicist trying to understand how a nucleus reacts to a collision, use the Bouncer. If you are studying the internal structure of a stable nucleus, you can use either, but you must remember that the Shape-Shifter will give you a slightly "heavier" result. The choice of method changes the details of the house, even if the overall shape looks the same.

Technical Summary

1. Problem Statement

In few-body nuclear physics, simplified models often treat a nucleus as a "core" plus external nucleons (e.g., a core plus two neutrons). However, the fermionic nature of nucleons imposes antisymmetry constraints. The effective potentials used in these models often support deeply bound states that correspond to internal core nucleons. These states are Pauli-forbidden for the external nucleons.

If not handled correctly, these forbidden states contaminate the physical spectrum. The paper addresses the comparison of two primary methods for eliminating these Pauli-forbidden states in three-body calculations:

1. Pauli Projection (PP): Adding a strong, nonlocal repulsive term to the potential to project out forbidden states.
2. Supersymmetric (SS) Transformation: Applying a transformation to the potential to create a local repulsive $r^{-2}$ term that removes the deepest states while preserving phase equivalence for higher states.

While previous studies suggested differences, it remained unclear whether the SS method (which preserves phase equivalence) yields results distinct from the PP method, particularly in scattering observables and resonance properties.

2. Methodology

The author employs a rigorous momentum-space framework using Faddeev-type equations to solve for bound, resonant, and scattering states in systems consisting of a nuclear core and two nucleons.

• Equations: The study utilizes the Alt-Grassberger-Sandhas (AGS) version of the Faddeev equations for transition operators ($U_{\beta\alpha}$). These are solved in the partial-wave representation.
• Potentials:
  • Two-nucleon interaction: High-precision charge-dependent Bonn (CD Bonn) potential.
  • Nucleon-core interaction: Standard Woods-Saxon potential with spin-orbit terms.
  • Pauli Elimination:
    • PP Model: Modifies the potential with a nonlocal projection term $v_P = v + |\phi_P\rangle \Gamma_P \langle \phi_P|$, where $\Gamma_P$ is a large repulsive strength.
    • SS Model: Transforms the potential to generate a local repulsive core that behaves like $r^{-2}$ at short distances, eliminating the forbidden state while keeping the spectrum of allowed states identical (phase equivalent).
• Systems Studied:
  • Scattering: Deuteron-$^4$He elastic scattering and breakup.
  • Bound States: Weakly bound halo nuclei ($^{11}$Li, $^{19}$B) and tightly bound systems with excited states ($^{16}$C, $^{18}$O, $^{20}$C).
  • Resonances: The unbound nucleus $^{16}$Be.
• Coulomb Handling: For charged systems (deuteron-$^4$He), the screening and renormalization method is used to include Coulomb forces.

3. Key Contributions

• Systematic Comparison: Unlike previous works that focused on specific halo nuclei, this paper provides a comprehensive comparison of PP and SS methods across a wide range of nuclear systems (weakly bound, tightly bound, resonant, and scattering).
• Momentum-Space Implementation: The study rigorously implements the nonlocal PP potentials and local SS potentials within the same momentum-space integral equation framework, ensuring a fair comparison without coordinate-space discretization artifacts.
• Resolution of Phase Equivalence Debate: The paper clarifies whether the phase-equivalent nature of the SS method leads to different physical predictions compared to the nonlocal PP method, specifically in the context of continuum and resonance states.

4. Key Results

A. Scattering Observables (Deuteron-$^4$He)

• Experimental Preference: Experimental data for deuteron-$^4$He elastic scattering and breakup strongly favor the Pauli Projection (PP) model.
• SS vs. RC: The SS model results are very close to the Repulsive Core (RC) model (which violates phase equivalence) and both differ significantly from the PP model.
• Conclusion: The PP method is superior for reproducing scattering data, suggesting that the specific nonlocal structure of the projection term is physically more relevant than the local $r^{-2}$ singularity of the SS transformation in this continuum regime.

B. Bound States (Ground and Excited)

• Binding Energies: The SS method systematically predicts larger binding energies (more bound) for ground states and $2^+$ excited states compared to the PP method. Conversely, for $0^+_2$ excited states, the PP method predicts larger binding energies.
• Partial Wave Dependence:
  • When only low partial waves (containing Pauli-forbidden states) are included, the results are similar (due to Efimov physics/universality).
  • The divergence arises when higher partial waves ($L \ge 2$) are included. The SS model shows a significantly larger increase in binding energy as higher waves are added.
• Kinetic Energy & Momentum:
  • The PP model yields a significantly higher expectation value for kinetic energy ($E_K$) than the SS model.
  • This is attributed to the PP method explicitly enforcing orthogonality to the forbidden state, inducing higher momentum components.
  • The SS model allows for a larger weight of large-distance components in low partial waves, making the coupling to higher waves more efficient.
• Momentum Distributions: The core-neutron relative momentum distributions differ drastically. The forbidden state in the PP model shows a vanishing overlap (zero at the peak), whereas the SS model retains a non-zero overlap weight (0.166) and lacks a deep minimum in the momentum distribution.

C. Resonances ($^{16}$Be)

• The SS model predicts resonance positions at lower energies (more bound) than the PP model, consistent with the trend observed in bound states.
• Both models predict a very weakly bound excited state in $^{20}$C ($0^+_2$) that is not yet observed experimentally.
• Current potentials overestimate resonance energies compared to experimental data, suggesting a need for attractive three-body forces.

5. Significance and Conclusions

• No Single "Best" Method: There is no universal superiority of one method over the other for all observables.
  • Scattering: The PP method is preferred by experimental data.
  • Bound/Resonant States: The methods yield systematic differences rather than clear superiority. The choice of method affects the predicted binding energies and momentum distributions.
• Physical Implications: The differences stem from how the two methods handle the wave function's momentum components and orthogonality. The PP method's explicit orthogonality constraint leads to higher kinetic energy and different momentum distributions, which could significantly impact breakup reaction amplitudes.
• Future Directions: To draw definitive conclusions for exotic nuclei (like $^{16}$Be or $^{20}$C), better-constrained nucleon-nucleus potentials are required. The study highlights that the interplay between low and high partial waves is critical in determining the differences between PP and SS approaches.

In summary, while the SS method offers mathematical elegance through phase equivalence, the PP method appears more robust for describing scattering data. For bound states, the choice of method introduces systematic shifts in binding energies and wave function structures that must be accounted for in precision nuclear structure calculations.