7.3: Central Forces
( \newcommand{\kernel}{\mathrm{null}\,}\)
Separation of Variables
Working in spherical coordinates is significantly more difficult than working in Cartesian coordinates. So why do it? Because point-like particles are sources for spherically-symmetric potentials that affect other particles. Understanding how to work in spherical coordinates is essential for solving the hydrogen atom in particular. First a reminder of the coordinates themselves:
Figure 7.3.1 – Spherical Coordinates
We begin with the stationary-state Schrödinger equation in three dimensions. We need to write the Laplacian operator found in the Hamiltonian (Equation 7.1.7) in spherical coordinates. This we can simply look up (deriving it from the transformation between coordinate systems is no fun):
∇2=1r2[∂∂r(r2∂∂r)+(1sinθ)∂∂θ(sinθ∂∂θ)+(1sin2θ)∂2∂ϕ2]
The stationary-state Schrodinger equation in spherical coordinates is therefore:
−ℏ22mr2[∂∂r(r2∂∂r)+(1sinθ)∂∂θ(sinθ∂∂θ)+(1sin2θ)∂2∂ϕ2]ψ(r,θ,ϕ)+V(r,θ,ϕ)ψ(r,θ,ϕ)=Eψ(r,θ,ϕ)
Our goal is to collect similar variables, and before we use our usual trick, we can get a head start on that by removing the r−2 factor from the second and third terms. Collecting all the terms with factors of r on one side of the equation then gives:
(1sinθ)∂∂θ(sinθ∂∂θ)ψ(r,θ,ϕ)+(1sin2θ)∂2∂ϕ2ψ(r,θ,ϕ)=2mr2ℏ2[V(r,θ,ϕ)−E]ψ(r,θ,ϕ)−∂∂r(r2∂∂r)ψ(r,θ,ϕ)
Now follow our two-step separation trick: Write the wave function as a product of single-variable functions...
ψ(r,θ,ϕ)=R(r)Θ(θ)Φ(ϕ)
...and divide the whole equation by the wave function:
1Θ(1sinθ)ddθ(sinθddθ)Θ+1Φ(1sin2θ)d2dϕ2Φ=2mr2ℏ2[V(r,θ,ϕ)−E]−1Rddr(r2ddr)R
We haven't quite separated the angular variables yet, but let's come back to that in a bit. First, we note that the only thing keeping the right-hand side of the equation from containing only functions of r is the potential. This method therefore only helps us when the potential is a function of r, which is the case for point-source potentials, as discussed above. With the left-hand side only a function of the angular coordinates, and the right hand side only functions of r, they must both equal a common constant (which we will for the time being call "C1"), giving us (after a bit of rearranging) what is known as the radial equation:
−ℏ22mr2ddr(r2ddr)R(r)+ℏ2C12mr2R(r)+V(r)R(r)=ER(r)
Plugging C1 into the right-hand side of Equation 7.3.5, we now set out to separate the angular functions:
1Θ(1sinθ)ddθ(sinθddθ)Θ+1Φ(1sin2θ)d2dϕ2Φ=C1
Multiply the equation by sin2θ and collect the functions of each variable to get:
1Θ[sinθddθ(sinθddθ)Θ−C1sin2θ]=−1Φd2dϕ2Φ
Separated at last, we can state that both sides of this equation equal a constant (which will will call "C2"), giving us two more differential equations in single variables, the polar equation (in θ) and the asimuthal equation (in ϕ):
sinθddθ(sinθddθ)Θ(θ)−(C1sin2θ+C2)Θ(θ)=0
d2dϕ2Φ(ϕ)+C2Φ(ϕ)=0
While the solution to the radial equation depends upon the details of the potential, the solutions to the polar and azimuthal equations do not (other than the fact that the potential is strictly a function of r). We have seen the azimuthal differential equation before, and we can write the solution this way:
Φ(ϕ)=Ae±i√C2ϕ
We cannot go any further to determine constants in this solution without boundary conditions. In our previous encounters with this differential equation, we found our boundary conditions from the potential (e.g. the wave function must vanish at the walls of the infinite square well). That is not the case here, but we do have a condition we can impose: If r and θ are held fixed and we change ϕ by some multiple of 2π, then we come back to the same point in space. We therefore insist that the function Φ returns to the same value every time ϕ changes by some multiple of 2π. This is ensured by insisting that the constant √C2 is an integer. For reasons that will become clear later (and with the ancillary benefit of not confusing it with the mass of the particle), we name this integer "ml," giving us:
Φml(ϕ)=Aeimlϕ,ml=0,±1,±2,…
The solution to the polar equation, while it doesn't depend upon the potential, is no picnic, and we will not delve into that bit of mathematics. Its solution is called a Legendre polynomial, and like the azimuthal case, there is a periodic boundary condition involved, which once again introduces an integer-valued variable, this time related to the constant C1, which we call "l." Specifically, it turns out:
C1=−l(l+1),l=0,1,2,…
Notice that the polar equation also includes the constant C2, which means its solution also depends upon ml. Indeed, this links the integers l and ml intimately, which explains the subscript on ml. For this reason, the functions Θlml(θ) and Φml(ϕ) are usually thrust together to make a single function (called spherical harmonic functions, which have been solved, and can simply be looked-up) of both variables, and this new function includes both quantum numbers (typically the "l" subscript in ml is suppressed in the variable representing this function, as the relationship between the two quantum numbers m and l is understood, and in this context m is not going to be confused with the mass):
Ylm(θ,ϕ)=Θlml(θ)Φml(ϕ)
[Note: In the literature, the indices defining spherical harmonics are frequently written with the ml index raised: Yml.]
For a given value of l, it so happens that ml can take on only a limited number of values:
That is, the absolute value of ml can never exceed l.
Plugging in for the value of C1 removes the unknown constant from the radial equation:
−ℏ22mr2ddr(r2ddr)R(r)−ℏ2l(l+1)2mr2R(r)+V(r)R(r)=ER(r)
If the particle is bound, then this equation will give quantized partial wave functions Rnl(r) (from the presence of the l(l+1) in the radial equation, it seems clear that the radial wave function will have dependence on the quantum number l as well as the quantum number n that comes from its own differential equation), and we will once again have three quantum numbers: n, l, and ml.
Just for completeness, we can put the wave function of the stationary state back together:
ψnlml(r,θ,ϕ)=ARnl(r)Ylm(θ,ϕ)
All that remains is to plug in the radial potential and use the boundary conditions to solve the differential equation for Rnl(r).
While these names may not make sense just yet, these quantum numbers are referred to as the principle quantum number (n), the orbital quantum number (l), and the magnetic quantum number (ml).
The Free Particle
From our work with the free particle in Cartesian coordinates, it is clear that it is the easiest coordinate system to use when one wishes to discuss plane waves. But as we have mentioned many times before, plane waves are not the only energy eigenstate wave functions. There is in fact a common energy eigenstate for which spherical coordinates are ideal. It is called a spherical wave, as it emanates radially outward from (or inward toward) a single point, which of course is the origin of our spherical coordinates. This is more often encountered when discussing light rather than massive particles (as light is frequently emitted from a point source), but as we have seen before, quantum mechanics treats both as "quanta."
For reasons we will see when we discuss the hydrogen atom next, a radially spreading (or contracting) wave function has an l-value equal to zero (and therefore a zero ml-value as well), and it is left as an exercise to show that this wave function satisfies the free-particle (V(r)=0), radially-outward (l=0) radial equation:
ψ(r)=Ae±ikrr,E=ℏ2k22m
The positive-valued exponent corresponds to the wave that moves radially-outward. Notice that the probability density is simply 1r2, which reflects an inverse-square reduction in intensity when one converts this single particle wave function to a steady-state stream of particles. Also notice that like the plane wave, this wave function cannot be normalized. And finally, as is always the case with energy eigenstates, we can attach the time-dependence by multiplying the stationary-state solution by the usual e−iωt.