Part 13, 14, 15
Outline, Schedule

PHYS 424 Notes

Part 14

  1. Chapter 4
    1. Spherical Coordinates
      1. Three dimensions

        To go to three dimensions in the Schroedinger equation, simply repeat what we have done before using vectors:

        Hy = i hbar dy/dt

        H = px2/(2m) + V(x)         H = p2/(2m) + V(r)

        In general, the potential is V(r), but in most interesting cases the dependence will be on the magnitude of r only. With

        px = - i hbar d/dx         p = - i hbar Ñ

        so that the Schroedinger equation is

        - [hbar2/(2m)] Ñ2 y + V(r)y = E y

        From calculus, Ñ = 1r /r + 1q (1/r) /q + 1f(1/r sinq) /f, where each "denominator" is ds along the indicated direction. Unfortunately, Ñ2 is not so simply interpreted, although the only really tricky term is /q:

        Ñ2 = (1/r2)/r (r2/r) + 1/(r2sinq) /q (sinq /q) + 1/(r2 sin2q) 2/f2

        There is an overall 1/r2, with the placement of factors of r inside derivatives given by units, and the /f term is completely natural.

      2. Separation of variables

        Take y = R(r) Y(q, f), substitute into the Schroedinger equation, and multiply by
        -2 m r2 / (hbar2 R Y) to get

        (1/R) d/dr (r2 dR/dr) - (2 m r2 / hbar2)[V(r) - E] = - (1/Y){(1/sinq)d/dq (sinq dY/dq) + (1/sin2q) d2Y/df2} = l(l+1)

        where the last part of the equation is due to the fact that we have a function of r alone equal to a function of q, f alone. The form of the constant is chosen for later convenience. We can't do anything about the radial equation until we make some choices for V(r), but the angular equation can be further simplified by taking the form

        Y(q,f) = Q(q) F(f),

        substituting into the equation and dividing by QF.

        This time we get a function of q equal to a function of f so that each must be a constant. We choose to call the constant m2 [watch out, this is the same as the symbol for the mass], so that we have

        (1/Q) [sinq (d/dq) (sinq dQ/dq)] + l(l+1)sin2q = - (1/F) d2F/df2 = m2

        Altogether we have

        d/dr (r2 dR/dr) - (2 m r2 / hbar2) [V(r) - E + hbar2 l(l+1)/(2 m r2] R = 0
        sinq (d/dq) (sinq dQ/dq) + [l(l+1)sin2q - m2] Q = 0
        d2F/df2 + m2 F = 0

        The f equation gives     F(f) = eimf,     with m zero or a positive or negative integer.

        Because the equation depends on m2 instead of m, the solutions with positive and negative m are the required two solutions. The overall constant may be lumped into Q. The quantity m must be an integer so that F will be the same function at f and f + 2p.

        The q equation depends on both l and m, so its solution must also. For physical wavefunctions, l must be a nonnegative integer. The properties of the solutions are given in the text. The notation is

        Q(q) = A Plm(cosq)     Q(q) F(f) = Ylm (q, f)

        and the immediately important properties are

        Ylm (q, f) = sqrt{[(2l+1)/4p] [(l-m)! / [(l+m)!]} eimf Plm(cosq),     m>0
        Yl-m (q, f) = (-1)m Ylm (q, f)
        ò 0 2p  df ei(m-m')f = (2p)dmm'
        ò 0 p  sinq dq Plm(cosq) Pl 'm(cosq) = [2/(2l+1)][(l+m)!/(l-m)!]
        Note m's must be the same, and m³0.
        ò 0 p  ò 0 2p  [Yl 'm' (q, f)]* [Ylm (q, f)] sinq dq df = dll ' dmm'

        Ylm is more often written Ylm. The Ylm are a complete set of square-integrable angular functions.

Last Revised 01/10/18

Part 13, 14, 15
Outline, Schedule