Bound states of the $D$-dimensional Schrödinger equation for the generalized Woods-Saxon potential

This paper derives approximate analytical solutions for the bound states of the $D$-dimensional Schrödinger equation with a generalized Woods-Saxon potential using the Pekeris approximation, Nikiforov-Uvarov, and Supersymmetric quantum mechanics methods to obtain energy eigenvalues and wave functions, which are then applied to calculate the neutron system of a $^{56}\text{Fe}$ nucleus in two and three dimensions.

The Gist

Imagine the inside of an atomic nucleus as a crowded, bustling city. The particles (like neutrons) moving inside aren't just floating freely; they are trapped in a specific "neighborhood" created by the collective pull of all the other particles. Physicists call this neighborhood a potential well.

This paper is like a team of architects and mathematicians trying to draw a perfect map of how a single neutron moves within this neighborhood, specifically for a city made of Iron-56 (a common type of iron atom).

Here is the breakdown of their work in simple terms:

1. The Map They Are Using: The "Generalized Woods-Saxon"

For decades, physicists have used a specific shape to describe this nuclear neighborhood, called the Woods-Saxon potential. Think of it like a smooth, round bowl where the bottom is flat and the sides slope gently up.

However, the authors in this paper are using a "Generalized" version of this bowl. Imagine taking that smooth bowl and adding a little extra "lip" or a small dip right near the rim. This extra feature (called the surface term) is crucial because it helps explain how particles bounce off the edge of the nucleus or get temporarily stuck there.

2. The Problem: The "Centrifugal Force" Spin

The math gets tricky when the neutron isn't just sitting still; it's spinning or orbiting. In physics, this spinning creates a "centrifugal force" that tries to fling the particle outward, like water flying off a spinning wet dog.

In the real world, this force depends on the distance from the center ($1/r^2$). But the "Generalized Woods-Saxon" bowl has a very specific, curved shape that doesn't play nicely with that simple $1/r^2$ rule. It's like trying to fit a square peg into a round hole. The math becomes so complex that it's impossible to solve exactly for a spinning particle using standard methods.

3. The Solution: The "Pekeris Approximation"

To fix this, the authors use a clever trick called the Pekeris approximation.

Imagine you are trying to draw a perfect circle, but your hand is shaking. Instead of trying to draw the whole circle at once, you zoom in on a tiny spot and draw a straight line that looks almost like the curve. Then you move to the next spot and draw another line. If you do this enough times, you get a shape that is close enough to the circle for your purposes.

The authors do this with the spinning force. They replace the difficult "spinning rule" with a simpler formula that looks very similar to the shape of their nuclear bowl. This allows them to solve the math puzzle.

4. The Tools: Two Different Keys to the Same Lock

To solve the resulting equations, the authors use two different mathematical "keys":

• The Nikiforov-Uvarov (NU) Method: Think of this as a systematic, step-by-step recipe for solving complex equations.
• Supersymmetric Quantum Mechanics (SUSY QM): Think of this as a shortcut that uses the relationship between different types of energy states to find the answer without doing all the heavy lifting.

The Big Reveal: When they used both keys, they got the exact same answer. This is a huge deal in science. It means their map is reliable. If two different methods lead to the same destination, you know you aren't lost.

5. What They Found: The "Energy Menu"

Using their new map and tools, they calculated the specific energy levels (how much energy the neutron has) for a neutron inside an Iron-56 nucleus.

• Finite Menu: They found that the neutron can only exist at certain specific energy levels, like rungs on a ladder. It can't be anywhere in between.
• The Rules of the Ladder: They discovered that the number of rungs (bound states) depends on how deep the bowl is and how wide the "lip" is. If the bowl isn't deep enough, the neutron escapes (it's not a "bound state").
• Dimensions: They did this math for a 2-dimensional world (flat) and a 3-dimensional world (our reality), and then showed how to predict what would happen in 4, 5, or even higher dimensions just by shifting the numbers slightly.

6. The Results: What the Data Says

They ran the numbers for a neutron in an Iron-56 nucleus:

• Low Spinning (Low "l"): The neutron can stay trapped in the nucleus at several different energy levels.
• High Spinning (High "l"): As the neutron spins faster, the "centrifugal force" pushes it out. Eventually, the force is so strong that the neutron can no longer stay trapped, no matter how deep the bowl is. The "ladder" runs out of rungs.
• The "Unbound" Zone: For certain combinations of spinning and energy, the math says the neutron simply cannot exist in a stable state there. It's like trying to park a car on a cliff edge; it just falls off.

Summary

In short, this paper is a mathematical success story. The authors took a very difficult problem (describing a spinning particle in a complex nuclear bowl), used a smart approximation to simplify the geometry, and solved it using two different advanced methods that agreed with each other. They provided a precise list of energy levels for neutrons in an Iron nucleus, showing exactly which states are stable and which ones are impossible.

Technical Summary

1. Problem Statement

The paper addresses the challenge of finding analytical solutions to the hyper-radial Schrödinger equation for a particle moving in a Generalized Woods-Saxon (GWS) potential in $D$-dimensional space ($D \ge 2$).

• The Potential: The GWS potential includes a standard volume term and a surface term (derivative of the Woods-Saxon), defined as:
  $$V(r) = -\frac{V_0}{1 + e^{(r-R_0)/a}} + \frac{W e^{(r-R_0)/a}}{\left(1 + e^{(r-R_0)/a}\right)^2}$$
  where $V_0$ and $W$ are potential depths, $R_0$ is the radius, and $a$ is the surface diffuseness.
• The Difficulty: While the potential is analytically solvable for zero angular momentum ($l=0$), the inclusion of the centrifugal term ($\propto 1/r^2$) for $l \neq 0$ makes the Schrödinger equation non-exactly solvable using standard methods (Nikiforov-Uvarov, Supersymmetric Quantum Mechanics, etc.).
• The Goal: To obtain approximate analytical solutions for energy eigenvalues and hyper-radial wave functions for arbitrary angular momentum $l$ and arbitrary dimension $D$, and to analyze the bound state spectrum for a neutron in a $^{56}\text{Fe}$ nucleus.

2. Methodology

The authors employ a two-pronged approach to solve the problem, utilizing two distinct analytical methods to verify consistency.

A. The Pekeris Approximation

To handle the centrifugal barrier term $V_l(r) = \frac{\hbar^2 \tilde{l}(\tilde{l}+1)}{2\mu r^2}$ (where $\tilde{l} = l + \frac{D-3}{2}$), the authors apply the Pekeris approximation.

• They expand the centrifugal term in a Taylor series around the minimum of the effective potential ($r_{min}$).
• The $1/r^2$ term is approximated by a series of exponentials matching the form of the Woods-Saxon potential:
  $$\frac{1}{r^2} \approx \frac{1}{R_0^2} \left( C_0 + \frac{C_1}{1+e^{\alpha x}} + \frac{C_2}{(1+e^{\alpha x})^2} \right)$$
• The coefficients $C_0, C_1, C_2$ are determined by comparing the Taylor expansion of the centrifugal term with the expansion of the approximated form around the equilibrium point. This transforms the effective potential into a form solvable by standard analytical techniques.

B. Analytical Solution Methods

Once the effective potential is approximated, the authors solve the resulting differential equation using two independent methods:

1. Nikiforov-Uvarov (NU) Method:

   • The Schrödinger equation is transformed into a generalized hypergeometric-type differential equation via a coordinate transformation $z = (1 + e^{(r-R_0)/a})^{-1}$.
   • The method utilizes specific polynomials ($\sigma, \tilde{\sigma}, \tau$) to derive the energy eigenvalues and wave functions expressed in terms of Jacobi polynomials.

2. Supersymmetric Quantum Mechanics (SUSY QM):

   • The authors construct a superpotential $W(r)$ to generate partner potentials.
   • They utilize the concept of shape invariance to determine the ground state energy and wave function.
   • The excited state energies are derived by summing the remainder terms of the shape-invariant potentials.

Verification: The authors demonstrate that both methods yield identical expressions for the energy eigenvalues and that the wave functions derived from both methods are mathematically equivalent (transformable into one another).

3. Key Contributions and Results

A. General Energy Spectrum

The paper derives a closed-form expression for the energy eigenvalues $E_{nrl}^{(D)}$ for the GWS potential in $D$ dimensions. The energy depends on:

• Potential parameters: $V_0, W, R_0, a$.
• Quantum numbers: Radial $n_r$, orbital $l$, and dimension $D$.
• The derived formula shows that the energy spectrum is finite; bound states only exist if specific inequalities regarding the potential depth and quantum numbers are satisfied.

B. Dimensional Scaling Relations

A significant theoretical contribution is the identification of a scaling law connecting energy spectra across different dimensions. The authors show that the term $\tilde{l}(\tilde{l}+1)$ remains invariant under the transformation:
$$l \to l \pm N, \quad D \to D \mp 2N$$
This allows the calculation of energy spectra for high dimensions ($D > 3$) based on the results of $D=2$ and $D=3$:

• $E^{(4)}_{n_r, l} = E^{(2)}_{n_r, l+1}$
• $E^{(5)}_{n_r, l} = E^{(3)}_{n_r, l+1}$
• And so on.

C. Numerical Application: The $^{56}\text{Fe}$ Nucleus

The authors applied their formulas to calculate bound states for a neutron moving in the potential field of a $^{56}\text{Fe}$ nucleus ($A=56$).

• Parameters: $V_0 \approx 47.78$ MeV, $R_0 \approx 4.916$ fm, $a = 0.65$ fm.
• Cases Analyzed: Calculations were performed for $D=2$ and $D=3$ with surface potential depths $W = 100$ MeV and $W = 200$ MeV.
• Findings:
  • Existence of Bound States: Bound states exist only for specific combinations of $n_r$ and $l$. For $n_r=0$, bound states exist up to $l=5$ (for $W=100$) and $l=6$ (for $W=200$) in $D=2$.
  • Non-Physical Solutions: For $n_r \ge 1$, the authors found that the calculated energies often violate the condition $V_{eff, min} < E < V_l$, resulting in "unbound" or non-physical mathematical solutions. This suggests that for higher radial excitations, the standard GWS potential may require spin-orbit or pseudospin symmetry corrections to describe bound states accurately.
  • Dimensional Effect: Increasing the dimension $D$ or the orbital quantum number $l$ increases the energy of the bound states (pushing them closer to zero), effectively reducing the number of available bound states due to the repulsive centrifugal barrier.

4. Significance and Conclusion

• Methodological Validation: The paper successfully validates the consistency between the Nikiforov-Uvarov and SUSY QM methods for the GWS potential, correcting previous inconsistencies found in literature regarding the $l=0$ case.
• Physical Insight: It highlights the limitations of the standard GWS potential in describing high-excitation states ($n_r \ge 1$) without additional symmetry considerations (spin/pseudospin).
• Generalizability: The derived scaling relations for $D$-dimensions provide a powerful tool for studying quantum systems in higher-dimensional spaces without re-solving the differential equation from scratch.
• Nuclear Physics Application: The specific calculations for $^{56}\text{Fe}$ provide a benchmark for understanding single-particle energy levels and the behavior of neutrons in nuclear potentials under varying dimensional constraints.

In summary, this work provides a rigorous analytical framework for solving the D-dimensional Schrödinger equation with a generalized Woods-Saxon potential, offering both theoretical scaling laws and practical numerical results for nuclear systems.