Quantum Mechanics — Lecture notes for PHYS223

XIV Central potentials

The isotropic harmonic oscillator is an example of a spherical symmetric potential, where the potential energy only depends on the radial distance from the origin, r=|𝐫|. These potentials are called central potentials. Another example is the Coulomb potential

Vc=q1q24πε0r (235)

of two charges q1 and q2. For central potentials the Schrödinger equation

Eψ(𝐫)=-22μΔψ(𝐫)+V(r)ψ(𝐫) (236)

can also be separated in three one-dimensional problems, but only when one works in spherical polar coordinates. Note: Here mass is denoted by the symbol μ since we need m later for a different purpose.

XIV.1 Spherical polar coordinates

The three spherical coordinates are the radial coordinate r=x2+y2+z2, the azimuthal angle ϕ=arctan(y/x), and the polar angle θ=arccos(z/r).

Figure 1: Spherical polar coordinates.

The cartesian coordinates can be written as

x = rsinθcosϕ, (237)
y = rsinθsinϕ, (238)
z = rcosθ. (239)

The Laplace operator is given by

Δ=1r2r(r2r)+1r2sinθθ(sinθθ)+1r2sin2θ2ϕ2. (240)

We abbreviate the angular part as

Δθ,ϕ=1sinθθ(sinθθ)+1sin2θ2ϕ2. (241)

XIV.2 Separation of variables in spherical polar coordinates

We separate the radial dependence from the angular dependence in the wavefunction, ψ(x,y,z)=R(r)Y(θ,ϕ). From the Schrödinger equation (236) one finds that the angular part solves the equation Δθ,ϕY=-l(l+1)Y, i.e.,

1sinθθ(sinθYθ)+1sin2θ2Yϕ2=-l(l+1)Y, (242)

while the radial part solves the equation

-22μr2ddr(r2dRdr)+(V(r)+2l(l+1)2μr2)R=ER, (243)

where l is a constant.

XIV.3 Angular part

Since Eq. (242) for the angular part does not depend on the potential V(r) we can solve this equation once and for all. The solutions are the spherical harmonics

Ylm(θ,ϕ)=(-1)m(2l+1)4π(l-m)!(l+m)!Plm(cosθ)eimϕ, (244)

for m0, Yl,-m(θ,ϕ)=(-1)mYl,m*(θ,ϕ) for m<0.

Here l=0,1,2,3, is called the azimuthal quantum number and m=0,±1,±2,,±l is called the magnetic quantum number. Note that for each l there are 2l+1 values of m.

For m=0 the functions Plm(s) are the Legendre polynomials

Pl0(s)=12ll!dldsl(s2-1)l,l=0,1,2,3,. (245)

Hence

P00(s) = 1, (246)
P10(s) = s, (247)
P20(s) = (3s2-1)/2, (248)
P30(s) = (5s3-3s)/2. (249)

For m0 the functions Plm(s) are the associated Legendre polynomials

Plm(s)=Pl-m(s)=(1-s2)|m|/2d|m|ds|m|Pl0(s). (250)

The first spherical harmonics are

Y00(θ,ϕ) = 14πfor l=0m=0, (251)
Y10(θ,ϕ) = 34πcosθfor l=1m=0, (252)
Y1,±1(θ,ϕ) = 38πsinθe±iϕfor l=1m=±1, (253)
Y20(θ,ϕ) = 516π(3cos2θ-1), (254)
Y2,±1(θ,ϕ) = 158πsinθcosθe±iϕ, (255)
Y2,±2(θ,ϕ) = 1532πsin2θe±2iϕ, (256)

XIV.4 Quantisation of angular momentum

The angular momentum operator is given by 𝐋^=𝐫^×𝐩^, hence L^x=y^p^z-z^p^y etc. The squared length is given by 𝐋^2=L^x2+L^y2+L^z2. One can show that

𝐋^2Ylm(θ,ϕ)=2l(l+1)Ylm, (257)
L^zYlm(θ,ϕ)=mYlm. (258)

Hence the spherical harmonics are eigenfunctions of 𝐋^2 and L^z, with eigenvalues 2l(l+1) and m, respectively.

XIV.5 Degeneracy of the radial part

Equation (243) for the radial part R(r) depends on the potential and hence has to be solved separately for every problem in order to find the energies E. However, since this equation only depends on the azimuthal quantum number l but not on the magnetic quantum number m, each energy level is at least 2l+1 fold degenerate, since for each l there are 2l+1 solutions Ylm of the angular equation with different values of m.