(a) When I = 0, the radial wave function R(r) = x(r)/r satisfies the equation
the reduced mass being µ = M/2. By the change of variable
we can reduce the Schrodinger equation to Bessel's differential equation of order p
\(\frac{d^2J(x)}{dx^2} + \frac 1x \frac{dJ(c)}{dx} + \left(1 - \frac{p^2}{x^2}\right) J_p(x) = 0\)
Thus the (unnormalized) radial wave function is
(b) For bound states we require that for \(r \to \infty\), R(r) → 0, or Jp remains finite. This demands that p ≥ 0. R(r) must also be finite at r = 0, which means that x(0) = Jp(α) = 0.
This equation has an infinite number of real roots. For E = 2.2 MeV,
Figure shows the contours of Jp(x) for different values (indicated by right and top numbers) of the function in the x - p plane. The lowest zero of Jp(x) for p = 0.65 is 3.3, the next 6.6. Thus for \(\alpha \approx3.3\), the system has one l = 0 bound state, for which
which is a dimensionless constant.
(c) For \(\alpha \approx6.6\), there is an additional 1 = 0 bound state. Thus the minimum value of α for two l = 0 bound states is 6.6, for which