Exact Neutron-Proton Wavefunctions Using the Phase Function Method

This paper utilizes the Phase Function Method with optimized Morse potentials derived from the GRANADA analysis to calculate exact radial neutron-proton wavefunctions and phase shifts for various uncoupled channels, demonstrating excellent agreement with high-precision Nijmegen-II results across a wide range of laboratory energies.

The Gist

Imagine you are trying to understand how two tiny, invisible balls (a neutron and a proton) bounce off each other when they crash together. In the world of physics, scientists usually look at the "aftermath" of the crash—how much the balls scatter, how much energy they lose, or the angle they fly off at. They rarely get to see the actual "movie" of the collision happening in slow motion.

This paper by Anil Khachi is like a filmmaker who has figured out how to reconstruct that slow-motion movie, frame by frame, using a special mathematical camera called the Phase Function Method (PFM).

Here is a breakdown of what the paper does, using simple analogies:

1. The Goal: Reconstructing the Invisible Movie

Usually, physicists calculate a "phase shift." Think of this as a single number that tells you how much the collision "twisted" the path of the particles. It's like knowing a car took a sharp turn, but not seeing the road it drove on.

This paper goes a step further. Instead of just giving you the final turn number, it calculates the exact wavefunction.

• The Analogy: If the collision were a dance, the "phase shift" is just the final pose. The "wavefunction" is the entire choreography—the steps, the spins, and the movements at every single moment from the start of the dance to the end.
• The author calculates this dance for various "channels" (different ways the particles can spin and move relative to each other, labeled as S, P, and D waves).

2. The Tool: The "Morse" Trampoline

To calculate this dance, you need to know the rules of the interaction. What does the "floor" look like? Is it sticky? Is it bouncy? Is there a wall?

• The author uses a mathematical shape called the Morse Potential.
• The Analogy: Imagine the space between the two particles is a trampoline. Sometimes the trampoline dips down (attracting the particles together), and sometimes it has a stiff spring in the middle that pushes them apart (repulsion).
• The author didn't just guess the shape of this trampoline. He tuned it perfectly using a massive database of real-world experiment data (6,713 data points from 1950 to 2013). He adjusted the trampoline's springs until the math matched the real-world results perfectly.

3. The Method: The "Phase Function" Camera

The paper uses a technique called the Phase Function Method (PFM).

• The Analogy: Instead of trying to solve the whole dance at once (which is very hard), the PFM method builds the dance step-by-step as the particles get closer.
• It starts far away where the particles don't feel each other. As they move closer, the method calculates how the "dance steps" (the wave) change at every tiny fraction of a millimeter.
• It produces three things for every step of the way:
  1. Phase Shift (δ): How much the path has turned so far.
  2. Amplitude (A): How "loud" or strong the dance is at that point.
  3. Wavefunction (u): The actual shape of the dance at that specific distance.

4. The Results: Different Types of Dances

The author tested this method on different types of collisions (S, P, and D waves) and different speeds (energies).

• The S-Wave (The Straight-On Crash):

  • This is the simplest collision where the particles head straight for each other.
  • What happened: At low speeds, they are gently pulled together (like magnets). At high speeds, they hit a "hard core" in the middle that pushes them back. The paper shows exactly how the dance changes from a gentle pull to a hard bounce.
  • The Verdict: The author's "movie" matches the high-precision "movies" made by other famous physics teams (Nijmegen-II) almost perfectly.

• The P-Wave (The Glancing Blow):

  • Here, the particles have a bit of spin, so they don't hit head-on; they kind of graze each other.
  • What happened: Some of these collisions were purely "repulsive" (like two magnets with the same pole facing each other). The math showed the particles never really got close; they just bounced off an invisible wall. The author's method captured this "pushing away" perfectly.

• The D-Wave (The Complex Spin):

  • These are even more complex spins.
  • What happened: Because of the spin, there is a "centrifugal barrier" (like a spinning top that keeps things apart). The particles mostly feel the "middle" of the interaction, not the very center. The author's method worked very well here too, matching other experts' results.

5. The Conclusion: A Reliable New Camera

The paper claims that this "Phase Function Method" is a powerful, transparent, and accurate tool.

• Why it matters: It proves you can take a simple, well-tuned mathematical model (the Morse potential) and use this specific method to generate the exact wavefunctions of the collision.
• The Limitation: The paper admits it only looked at "uncoupled" states (simple dances where the spin doesn't get tangled with the orbit). It notes that "coupled" states (where the spin and orbit get tangled up, like a complex tango) are too complicated for this specific version of the math and will need to be studied in a future paper.

In summary: The author built a mathematical camera that films the invisible dance of neutrons and protons. By tuning the camera with real-world data, he produced a movie that looks exactly like the ones made by the most expensive, high-tech physics labs, proving that his simpler, step-by-step method works beautifully.

Technical Summary

Technical Summary: Exact Neutron-Proton Wavefunctions Using the Phase Function Method

Problem Statement
The primary objective of theoretical nuclear physics in scattering problems is often to derive the wavefunction, from which other observables like phase shifts, cross-sections, and low-energy parameters can be extracted. However, experimentalists do not directly measure wavefunctions; they measure differential and total cross-sections, which are related to phase shifts. While various approximation techniques (e.g., Born approximation, Brysk's approximation) exist to solve the Schrödinger equation for scattering phases, this work focuses on obtaining exact radial wavefunctions, amplitude functions, and distance-dependent phase shifts for neutron-proton ($n-p$) scattering. The study specifically addresses uncoupled S, P, and D channels, aiming to provide a transparent and accurate description of these states using the Phase Function Method (PFM) combined with inverse potentials derived from the Morse function.

Methodology
The study employs the Phase Function Method (PFM), a technique that transforms the second-order Schrödinger differential equation into a first-order non-homogeneous Riccati-type differential equation for the radial phase shift, $\delta_\ell(k, r)$.

1. Potential Model: The interaction is modeled using the Morse potential, $V(r) = V_0(e^{-2(r-r_m)/a_m} - 2e^{-(r-r_m)/a_m})$. This potential was chosen for its exact solvability, shape invariance, and the existence of an exact analytical expression for the $^1S_0$ state.
2. Parameter Optimization: The Morse potential parameters ($V_0$, $r_m$, $a_m$) were not fitted arbitrarily. Instead, they were optimized using the comprehensive GRANADA partial wave analysis, which comprises 6,713 experimental $n-p$ phase shift data points from 1950 to 2013. The optimization minimized the mean square error (MSE) between the model and experimental data.
3. Numerical Solution:
   • The Riccati equation for the phase shift $\delta_\ell(k, r)$ is solved numerically using the 5th-order Runge-Kutta (RK-5) method with the initial condition $\delta_\ell(0) = 0$.
   • Specific recurrence relations are used to adapt the Riccati-Bessel and Riccati-Neumann functions for different angular momentum states ($\ell = 0, 1, 2$), yielding specific forms of the PFM equation for S, P, and D waves.
   • Once the phase function is determined, the amplitude function $A_\ell(r)$ and the exact radial wavefunction $u_\ell(r)$ are calculated using coupled first-order differential equations derived within the PFM framework.
4. Scope: Calculations were performed for laboratory energies $E_{lab} = [1, 10, 50, 100, 150, 250, 350]$ MeV, with detailed radial dependence presented up to 5 fm. The study is restricted to uncoupled channels ($^1S_0, ^1P_1, ^3P_0, ^3P_1, ^1D_2, ^3D_2$), explicitly excluding coupled-channel effects (such as the $^3S_1-^3D_1$ mixing due to tensor forces).

Key Contributions

• Exact Wavefunctions: The paper provides explicit, distance-dependent radial wavefunctions $u(r)$, amplitude functions $A(r)$, and phase functions $\delta(r)$ for various uncoupled $n-p$ scattering channels, rather than just asymptotic phase shifts.
• Inverse Potential Application: It demonstrates the utility of using inverse Morse potentials, optimized against a massive experimental dataset (GRANADA), as a zeroth-order reference for high-precision scattering calculations.
• Validation against Nijmegen-II: The calculated wavefunctions are directly compared with the high-precision Nijmegen-II interaction results, serving as a benchmark for the accuracy of the PFM approach.

Results
The results are presented for laboratory energies of 10, 50, 150, and 250 MeV across six uncoupled channels:

• $^1S_0$ State: The radial phase shift $\delta_0(r)$ shows energy-dependent behavior, accumulating gradually at low energies but exhibiting reduced accumulation and partial reversal at small radii as energy increases, reflecting the short-range repulsive core. The wavefunctions show excellent agreement with Nijmegen-II, particularly at higher energies, with smooth asymptotic matching.
• $^1P_1$ State: The inverse potential is found to be purely repulsive. Consequently, the phase shift remains negative across the entire interaction range. The wavefunctions exhibit suppressed amplitude at short distances, characteristic of repulsive P-wave interactions, and align well with Nijmegen-II results.
• $^3P_0$ and $^3P_1$ States: These channels display distinct behaviors. The $^3P_0$ state shows a transition from negative to positive phase shifts, indicating an interplay between short-range repulsion and intermediate-range attraction. The $^3P_1$ state is dominated by repulsion, maintaining a negative phase shift. Both show good agreement with Nijmegen-II, with accuracy improving at higher energies.
• $^1D_2$ and $^3D_2$ States: Due to the centrifugal barrier, D-waves are less sensitive to the short-range core. Both states are predominantly attractive, resulting in positive phase shifts. The wavefunctions show very good agreement with Nijmegen-II, with minor discrepancies at short distances attributed to the simplified form of the Morse potential rather than the PFM method itself.

Significance and Claims
The paper claims that the Phase Function Method, when combined with inverse Morse potentials optimized against experimental data, provides a reliable, accurate, and physically transparent framework for describing uncoupled neutron-proton scattering states.

• Accuracy: The obtained wavefunctions show "excellent agreement" with the high-precision Nijmegen-II results, validating the PFM approach for generating realistic scattering wavefunctions.
• Transparency: The method offers a clear physical interpretation of how attractive and repulsive regions of the potential contribute to the total phase shift (positive vs. negative accumulation).
• Robustness: The absence of numerical instabilities or unphysical oscillations in the amplitude functions confirms the stability of the first-order PFM formulation.
• Limitations and Future Work: The authors modestly note that the current work is restricted to uncoupled channels. They acknowledge that coupled-channel effects (e.g., tensor interactions in the $^3S_1-^3D_1$ channel) are not included and require a more complex matrix-valued formulation. The extension to coupled channels and other two-body systems (e.g., $n\alpha$, $p\alpha$) is identified as a natural continuation of this study.