Last time, we found that the problem of the hydrogen atom could be split into a radial part and an angular part. Thanks to spherical symmetry, the angular part could be studied using angular momentum operators and spherical harmonics. We found that the 3D behavior of the electron could be reinterpreted as a 1D wavefunction of a particle in an effective potential which was the two-body interaction potential plus a “barrier” term which depended upon the angular momentum quantum number. Today, we’re going to solve the radial part of the problem and thereby find the eigenstates and eigenenergies of the hydrogen atom.

The technique we’ll employ has a certain charm, because we solved the first part, the angular dependence, using commutator relations, while as we shall see, the *radial* dependence can be solved with *anticommutator* relations.

**THE FAMILY OF COULOMB HAMILTONIANS**

We ended up with a family of Hamiltonians labeled by the angular momentum quantum number:

[tex]H_l^{\rm Coul} = -\frac{\hbar^2}{2\mu} \partial_r^2 -\frac{Ze^2}{r} + \frac{l(l + 1) \hbar^2}{2\mu r^2}.[/tex]

For each eigenstate of any [tex]H_l^{\rm Coul}[/tex], the hydrogen atom has a set of states labeled by the quantum number [tex]m[/tex], which takes the values

[tex]m = 0, \pm 1, \pm 2, \ldots, \pm l.[/tex]

However, because the original problem was spherically symmetric, the operator [tex]L_3[/tex] does not appear in the overall Hamiltonian (there’s no preferred axis), and so the energy will not depend upon [tex]m[/tex]. Geometrical symmetry implies that states of different [tex]m[/tex] but the same [tex]l[/tex] are *degenerate*; as we proceed, we’ll see an additional example of symmetry implying degeneracy unfold before us.

Carrying around all those constants in [tex]H_l^{\rm Coul}[/tex] would be a pain and not particularly useful, so let’s get rid of them. By appropriately rescaling our variables (craftily choosing “better rulers” for distance and energy) we can rewrite [tex]H_l^{\rm Coul}[/tex] in a cleaner form, preserving the essential features of the Hamiltonians without all that dimensionful baggage:

[tex]H_l^{\rm Coul} = -\partial_x^2 + \frac{l(l+1)}{x^2} – \frac{1}{x}.[/tex]

[I originally screwed this bit up; see the comments below.]

Now, we have a family of related Hamiltonians, labeled by a nonnegative integer. Why does this situation sound familiar? When we studied supersymmetry and shape invariance, we saw just such a hierarchy of Hamiltonians appear! So, let’s see if those ideas are applicable to this problem.

We should note one complication right away: when we explored SUSY QM and shape-invariant partner potentials, we worked with an original Hamiltonian whose ground-state energy was *zero.* Here, however, we’re looking at an electron *bound by an atom,* so the energies of the eigenstates will be *negative.* Shifts of the energy scale aren’t going to pose profound difficulties, but we should be on our guard.

**LIVING UP TO OUR SUPERPOTENTIAL**

Earlier, we took a *superpotential* and from it derived two *partner potentials.* Thus, if [tex]W[/tex] was the superpotential, then

[tex]A = \partial_x + W(x)[/tex]

and

[tex]A^\dag = -\partial_x + W(x),[/tex]

so our partner Hamiltonians would be

[tex]H^{(1)} = A^\dag A = -\partial_x^2 + W^2(x) – \partial_x W(x)[/tex]

and

[tex]H^{(2)} = AA^\dag = -\partial_x^2 + W^2(x) + \partial_x W(x).[/tex]

To study [tex]H_l^{\rm Coul}[/tex] in this framework, we’ll need a superpotential which, when put through these steps, yields both the [tex]1/x[/tex] Coulomb law and the angular momentum barrier. The superpotential will, therefore, depend upon [tex]l[/tex]. We choose the following:

[tex]W_l = -\frac{(l + 1)}{x} + \frac{1}{2(l + 1)}.[/tex]

This gives us a pair of operators for each value of [tex]l[/tex]:

[tex]A_l = \partial_x – \frac{l + 1}{x} + \frac{1}{2(l + 1)},[/tex]

and

[tex]A^\dag_l = -\partial_x – \frac{l + 1}{x} + \frac{1}{2(l + 1)}.[/tex]

The first Hamiltonian is given by one ordering,

[tex]H_l^{(1)} = A_l^\dag A_l,[/tex]

which works out to be

[tex]H_l^{(1)} = -\partial_x^2 + \frac{l(l+1)}{x^2} – \frac{1}{x} + \frac{1}{4(l+1)^2},[/tex]

and when we define a Hamiltonian with the other ordering,

[tex]H_l^{(2)} = A_l A^\dag_l,[/tex]

we get as the partner Hamiltonian

[tex]H_l^{(2)} = -\partial_x^2 + \frac{(l+1)(l+2)}{x^2} – \frac{1}{x} + \frac{1}{4(l+1)^2}.[/tex]

By comparing these two expressions, we get that

[tex]H_l^{(2)} = H_{l + 1}^{(1)} + \frac{1}{4(l + 1)^2} – \frac{1}{4(l + 2)^2}.[/tex]

Here, we see the shape-invariance criterion: [tex]H_l^{(2)}[/tex] is a Hamiltonian which depends upon the parameter [tex]l[/tex], and we can reach [tex]H_l^{(2)}[/tex] by taking its partner [tex]H_l^{(1)}[/tex], tweaking a parameter and adding an offset which does not depend upon [tex]x[/tex]. We can see the energy-axis shift mentioned above if we rewrite these operators in terms of the Coulomb Hamiltonian. First, we see that

[tex]H_l^{(1)} = H_l^{\rm Coul} + \frac{1}{4(l + 1)^2},[/tex]

so the operators [tex]H_l^{(1)}[/tex] and [tex]H_l^{\rm Coul}[/tex] will have the same eigenstates, although their eigenvalues will be separated by a term dependent upon [tex]l[/tex]. We also have

[tex]H_l^{(2)} = H_{l + 1}^{\rm Coul} + \frac{1}{4(l + 1)^2}.[/tex]

**BUILDING UP THE STATES**

We are now in a position to employ the diagrammatic approach we built up earlier. Let’s label the eigenstates of a particular Hamiltonian in the hierarchy by [tex]\nu[/tex], such that [tex]\nu = 0[/tex] for the ground state. Each state in our diagram will then be labeled by [tex]l[/tex] and [tex]\nu[/tex], which we write in Dirac notation as [tex]|\nu, l\rangle[/tex]. The isospectrality between partner potentials implies that state [tex]\nu[/tex] of Hamiltonian [tex]H_l^{(1)}[/tex] has the same energy as state [tex]\nu – 1[/tex] of [tex]H_l^{(2)}[/tex]. The state annihilated by [tex]A_l[/tex] is defined by

[tex]A_l |0, l\rangle = 0,[/tex]

which tells us that

[tex]H_l^{(1)} |0, l\rangle = 0.[/tex]

In turn, because of the energy shift, the Coulomb Hamiltonian satisfies

[tex]H_l^{\rm Coul} |0, l\rangle = -\frac{1}{4(l + 1)^2} |0, l\rangle.[/tex]

Likewise, for [tex]l + 1[/tex],

[tex]H_{l+1}^{(1)} |0, l + 1\rangle = 0.[/tex]

Using the shape-invariance condition yields

[tex]H_l^{(2)} |0, l + 1\rangle = \left[\frac{1}{4(l + 1)^2} - \frac{1}{4(l + 2)^2}\right]|0, l + 1\rangle.[/tex]

By the other key property, isospectrality, we know that [tex]H_l^{(1)}[/tex] has a state of the same energy, one notch higher in [tex]\nu[/tex]:

[tex]H_l^{(1}} |1, l\rangle = \left[\frac{1}{4(l + 1)^2} - \frac{1}{4(l + 2)^2}\right]|1, l\rangle.[/tex]

Again, the physical energy of the hydrogen atom can be found by adding a constant:

[tex]H_l^{\rm Coul} |1, l\rangle = -\frac{1}{4(l + 2)^2} |1, l\rangle.[/tex]

Notice how as we moved from [tex]\nu = 0[/tex] to [tex]\nu = 1[/tex], the value in the denominator went from [tex]l + 1[/tex] to [tex]l + 2[/tex]? The rest of the states can be built by working “leftwards” in the state diagram:

[tex]|\nu, l\rangle = A^\dag_l A^\dag_{l + 1} \cdots A^\dag_{l + \nu – 1} |0, l + \nu\rangle.[/tex]

This just means repeating the process we went through earlier, meaning that the denominator of the eigenvalue will depend on [tex]l + \nu + 1[/tex]. Note that the energy of the state labeled by [tex]\nu[/tex] and [tex]l[/tex] depends only upon their sum. If we restore the constants we scaled away earlier, we get that

[tex]E_{\nu l} = -\frac{\mu e^4 Z^2}{2\hbar^2} \frac{1}{(l + \nu + 1)^2},[/tex]

so we can call [tex]l + \nu + 1[/tex] by the new name [tex]n[/tex] and write

[tex]E_n = -\frac{\mu e^4 Z^2}{2\hbar^2} \frac{1}{n^2}.[/tex]

This is the energy of the [tex]n{\rm th}[/tex] eigenstate of the hydrogenic Hamiltonian. We could work out explicit expressions for states of our choice, but before we plug-and-chug any specific examples, let’s take a breather and look at a picture.

**THE SLIGHTLY BIGGER PICTURE**

The diagram we drew before now looks like this:

The eigenstate diagram shows *more symmetry* than we had expected on geometrical considerations alone: the only quantum number necessary for determining the energy is [tex]n[/tex], which corresponds to the energy-level index in Bohr’s model of the atom. If our interaction potential were not shape invariant, this degeneracy would be broken.

This is the way I learned to solve the hydrogenic atom in the misty days of my undergraduacy. The only textbook I know of which takes an approach like this is Ohanian’s *Principles of Quantum Mechanics*; other than a handful of universities, most schools attack the problem by plowing into Schrödinger’s second-order differential equation and eventually finding a recursion relation for the Laguerre polynomials. Prof. Rajagopal‘s lecture notes call the standard method “much more painful,” and as for why most textbooks follow that route, “Go figure.” I suspect that too many teachers of quantum mechanics have been bitten by the Matrix Zombie and think that mathematics beyond differential equations is just too hard for introductory classes. Rather than making the time investment necessary to use “more advanced” techniques, they solve problems in laborious and rather unilluminating ways.

Unfortunately, MIT’s OpenCourseWare project doesn’t provide the lecture notes we used, or any later editions thereof; the site for 8.05 Quantum Physics II just lists the sections of textbooks which should be read, instead of providing actual juicy PDFs. This post, in particular, was based on the 8.05 material, while my earlier overview of the general superalgebra machinery mostly follows Fred Cooper, Avinash Khare and Uday Sukhatme’s review article, “Supersymmetry and Quantum Mechanics” (1994). As that review explains, Schrödinger himself solved for the hydrogen atom eigenstates with a method rather like this, in 1940; many years later, the supersymmetric context of that “factorization” method was discovered.

From here, we can go in several directions. After perhaps working a few examples, we can head towards the relativistic regime and find SUSY-based solutions to the Dirac Equation. Also, we can look back at classical mechanics and relate these ideas to the Laplace-Runge-Lenz vector, an avenue which will eventually lead us to *superalgebras with central charge* and BPS bounds. I’m also strongly tempted to look at the application of SUSY to diffusion problems via the Fokker-Planck Equation.

**SUSY QM SERIES:**

- Part 1: Superalgebra
- Part 2: Shape Invariance
- Part 3: Two-Body Problem
- Part 4: Separation of Variables
- Part 5: The Hydrogen Atom
- Part 6: Intermezzo: The Dirac Equation
- Part 7: The 1D Dirac Hamiltonian

## 3 Comments

Okay, you chose [tex]x^2=\frac{2\mu}{\hbar^2}r^2[/tex] and that rescales the differential term and the abgular momentum barrier, but why does it handle the potential term, getting rid of the [tex]Ze^2[/tex]?

(did I get the TeX right this time?)

Argh. This is what I get for writing too fast and not proofreading the whole thing. The scaling which

shouldbe done is, writing [tex]a[/tex] for the Bohr radius,[tex]x = r\left(\frac{2Z}{a}\right),[/tex]

to rescale the distance,

and[tex]\tilde{E}_{\nu l} = \frac{\hbar^2}{2\mu e^4 Z^2} E_{\nu l},[/tex]

to rescale the energy.

Uh-oh, I just got a referral to this page from a Google search, {dirac equation and susy qm}. Looks like I better get my act together and write that next post! In the meantime, an introduction can be found in the review article cited above.