Physlets run in a Java-enabled browser on the latest Windows & Mac operating systems.

If Physlets do not run, click here for help on updating Java and setting Java security.

# Section 13.7: The Coulomb Potential for the Idealized Hydrogen Atom

Please wait for the animation to completely load.

Now consider the radial part of the Schrödinger equation in Eq. (13.21) written as

[−(*ħ*^{2}/2μ)(1/*r*^{2}) *d*/*dr* (*r ^{2}*

*d*/

*dr*) +

*l*(

*l*+ 1)

*ħ*

^{2}/(

*2μr*) +

^{2}*V*(

*r*)]

*R*(

*r*) =

*E*

*R*(

*r*) . (13.26)

We want to rewrite these terms, especially the derivative, into a more standard form. The substitution, *R*(*r*) = *u*(*r*)/*r*, simplifies the derivative term, in the above equation to yield

[−(*ħ*^{2}/2μ) (*d*^{2}/*dr*^{2}) + *l*(*l* + 1)*ħ*^{2}/(2μ*r*^{2}) +* V*(*r*)] *u*(*r*) = E*u*(*r*) .

When we find solutions to this equation, we must keep in mind that we are solving for *u*(*r*), not *R*(*r*), and that we must divide *u*(*r*) by *r* to give the true radial wave function, *R*(*r*). Before looking at a particular *V*(*r*), we look at the general equation and interpret terms. We see what we can interpret as an *effective potential*:

*V*_{eff }= *l*(*l* + 1)*ħ*^{2}/(2μ*r*^{2}) + *V*(*r*) ,

where *l*(*l* + 1)*ħ*^{2}/(2μ*r*^{2}) is the *potential* associated with the so-called the centrifugal barrier. Now consider the following Coulomb potential, *V* = −*e*^{2}/*r*, which describes the potential energy function for an electron in the proximity of a proton: the potential responsible for the basic structure of the hydrogen atom.^{4} When we insert this Coulomb potential in the radial differential equation, we have a differential equation that describes the electron:

[−(*ħ*^{2}/2μ_{e}) (*d*^{2}/*dr*^{2}) + *l*(*l* + 1)*ħ*^{2}/(2μ_{e}*r*^{2}) − *e*^{2}/*r*] *u*(*r*) = *Eu*(*r*) . (13.27)

where μ_{e} is the electron's mass and *e* is the charge of the electron. In Animation 1 the effective potential for the Coulomb problem *V*_{eff }= *l*(*l* + 1)*ħ*^{2}/(2μ_{e}*r*^{2}) + *V*(*r*) is shown for *l* = 0, 1, 2. Notice that as *l* gets bigger, the centrifugal barrier increases as well. To get this equation into standard form, divide by −*ħ*^{2}/2μ_{e}, which yields

[ (*d*^{2}/*dr*^{2}) − *l*(*l* + 1)/*r*^{2 }+ 2μ_{e}*e*^{2}/(*ħ*^{2}*r*)] *u*(*r*) = (−2μ_{e}*E*/*ħ*^{2}) *u*(*r*) .

We now define κ^{2 }= (−2μ_{e}*E*/*ħ*^{2}) (which is real since *E* < 0) and the dimensionless quantities ρ = κ*r* and ρ_{0 }= 2μ_{e}*e*^{2}/(κ *ħ*^{2}). Making these substitutions yields

[ (*d*^{2}/*d*ρ^{2}) − *l*(*l* + 1)/ρ^{2 }+ρ_{0}/ρ − 1] *u* = 0 .

We begin our analysis of the solutions of this differential equation by considering the two special limiting cases:

**Case I:** ρ→ 0 (r→ 0). In this case the centrifugal barrier dominates in Eq. (13.46)

[ (*d*^{2}/*d*ρ^{2}) − *l*(*l* + 1)/ρ^{2}]
*u* = 0 . (13.28)

We find that the general solution to this equation is *u* = *A*ρ^{l+1}+*B*ρ^{−l }and therefore the normalizable piece is just *u* ∝ ρ^{l+1 }∝ *r*^{l+1}.

**Case II:** ρ → ∞ (*r* → ∞). In this case the centrifugal term, *l*(*l* + 1)/ρ^{2}, and the potential energy, ρ_{0}/ρ, vanish from Eq. (13.28) at large ρ. This leaves

[ (*d*^{2}/*d*ρ^{2}) − 1] *u* = 0 . (13.48)

which for *E* < 0 gives the normalizable solution *u* = exp(−ρ) = exp(−κ*r*).

Now that we have an idea of what the bound states should look like asymptotically, we can find the entire solution. After much algebra we first find that

*E* = μ_{e}e^{4}/(2*n*^{2}*ħ*^{2}) = −*R* (1/*n*^{2}) ,

where *R *= μ_{e}*e*^{4}/(2*ħ*^{2}) is the Rydberg and is 13.6 eV. This result describes the energy levels for the Coulomb problem, and hence, the basic energy level structure for the hydrogen atom. We now simplify ρ. We use ρ = κ*r* and the definition of κ to find

ρ = μ_{e}*e*^{2}/(*nħ*^{2}) *r* = (*r*/*na*_{0}) ,

where *a*_{0 }= *ħ*^{2}/μ*e*^{2} is the Bohr radius. We can make use of further substitutions, this time yielding the radial wave functions

*R*_{nl}(*r*) = *A*_{nl }e^{−r/na}0 [(r/*na*_{0})^{l+1}/*r*] *v*_{n}(*r*/*na*_{0}) , (13.29)

where *A*_{nl }= [(2/(*na*_{0})^{3}(*n* − *l* − 1)!/(2*n*[(*n* + *l*)!]^{3})] is the normalization for the radial energy eigenfunction. In addition, *v*_{n}( *r*/*na*_{0}) = *L*^{2l+1}_{n−l−1}(2*r*/*na*_{0}) are the associated Laguerre polynomials. The *unnormalized* radial wave functions, above without *A*_{nl}, are shown in Animation 2. In the animation, distances are given in terms of Bohr radii, *a*_{0}. You may enter values of *n* and *l* and see the radial energy eigenfunction that results.

We find that the entire energy eigenfunction, properly normalized, is simply the product of the radial and angular solutions:

ψ_{nlm }= *R*_{nl}(r) *Y*_{l}^{m}(θ,φ) .

Note that given that there are *n*^{2} states per *n* value, and that the energy just depends on *n*, the solutions have an *n*^{2} energy degeneracy.

*When you get a good-looking graph, right-click on it to clone the graph and resize it for a better view.*

^{4}In Chapter 14 we discuss corrections to the Coulomb potential which are responsible for the remaining structure in the hydrogen spectral lines. We also generalize the Coulomb potential to include *hydrogenic atoms*, those with one electron and *Z* protons.

« previous

next »