SUSY QM 2: Shape Invariance

Blogging on Peer-Reviewed ResearchAfter a while, you just get tired. An honest science blogger can only handle so much science jargon thrown around without meaning, only a limited amount of Choprawoo and quantum flapdoodle. How long can anyone with integrity, curiosity and a dose of genuine knowledge endure the trumpeting that, say, the brain’s limited ability to recover after injury is evidence for some quantum spirit? The brain is living flesh, made of living cells: by the same so-called logic, the scabbed knees of childhood are all evidence of quantum skin.

After a while, a deep reserve of psyche cries out, “Enough! If my freedom means aught, I must stop responding to these charlatans and move beyond. I must send a message of my own, a message which is not a reaction but an expression unto itself. I must sing the quantum genuine!”

In my case, this means another post on supersymmetric quantum mechanics.

RECAP

Last time, we deduced some interesting properties of Hamiltonians which can be factored into operators and adjoints:

H_1 = A^\dag A,\ H_2 = AA^\dag.

We observed that H_1 and H_2 are isospectral. That is, while the forms of their eigenfunctions may be different, the eigenvalues associated with those functions are the same; or, in physical terms, the wavefunctions have different shapes, but the energies match. The only exception is the ground state: if H_1 has a zero-energy ground state, then H_2 will not. Furthermore, the operator A maps eigenstates of H_1 into those of H_2, and the operator A^\dag maps eigenstates in the reverse direction:



We can sketch the full set of H_1 eigenstates as a ladder with, potentially, an infinite number of rungs. The eigenstates of H_2 sit right alongside them; the bottom “rung” is missing, as we mentioned earlier.

That “missing rung” has an important implication: instead of mapping the ground state of H_1 into a state of H_2, the operator A annihilates it.

A\ket{\psi_0^{(1)}} = 0.

We could, therefore, solve for the wavefunction \psi_0^{(1)}(x) by solving a first-order, linear differential equation — a simpler task than tackling the full, second-order Schrödinger Equation.

Now, we’ll make a qualitative statement which we’ll firm up in a moment. If the Hamiltonian H_2 were “roughly similar” to H_1, then we should be able to “pull the same trick again,” writing it as

H_2 = A_2^\dag A_2,

and inventing a third Hamiltonian,

H_3 = A_2 A_2^\dag,

which will be isospectral with H_2 except for a ground state satisfying

A_2\ket{\psi_0^{(2)}} = 0.

Then, why not, we do the same thing a whole bunch of times over again, building a whole hierarchy of Hamiltonia stretching off beyond the margin of our page.

Note that the “bottom rung” in each “ladder” satisfies a first-order differential equation, and that we can build any state of H_1 by stacking up a sequence of A^\dag operators. Start with a state annihilated by some A_n and then apply

A^\dag_{n-1}A^\dag_{n-2}\cdots A^\dag_{1}

to get an eigenstate of H_1 there on the left edge. This is reminiscent of the way we solve the harmonic oscillator using creation and annihilation operators, with the difference that instead of building up, we’re working in from the side.

DEFINITION OF SHAPE INVARIANCE

To see this in action, we need to solidify the notion of rough resemblance we invoked above. Suppose we have a potential V_1 associated with some system of interest. If we compute the superpotential W(x) and find that the partner potential V_2 has a similar shape to V_1, we can calculate the energy states of V_1 in a remarkably straightforward way. In this context, “similar shape” means that if V_1 depends upon a set of parameters, we can reach V_2 by substituting different values of those parameters into V_1. More mathematically, two partner potentials V_1 and V_2 are said to be shape invariant if they satisfy

V_2(x;a_1) = V_1(x;a_2) + R(a_1)

where a_1 and a_2 are sets of parameters, and a_2 is uniquely determined by a_1 (that is, a_2 = f(a_1) for some function f). The residual term R(a_1) does not depend upon the coordinate x. We shall see that if H_1 has a shape-invariant partner potential in H_2, the entire spectrum of H_1 can be computed through operator manipulations.

In the material which follows, we’ll work in “natural units,” to save ourselves from carrying around excess notational baggage. Our operators will therefore take the form

A = \partial_x + W(x)

and

A^\dag = -\partial_x + W(x).

HIERARCHY OF HAMILTONIANS

We are usually interested in examining the physics of a particular system, for which we have written an H_1. (We are restricted to work in one dimension; however, many higher-dimensional problems can be attacked through separation of variables. When studying the hydrogen atom, the radial dependence may be factored out of the Coulomb Hamiltonian to give a pseudo-one-dimensional potential; when studying the Landau effect, a convenient choice of gauge shows that along one axis, electrons behave as if in a harmonic oscillator potential.) As mentioned above, we can slide the energy scale up and down, so for convenience we fix the n = 0 eigenstate of H_1 to have zero energy:

E_0^{(1)}(a_1) = 0.

We can solve for the ground state knowing the superpotential to show that

\psi_0^{(1)}(x;a_1) = C\exp\left(-\int^x W(y;a_1)dy\right)

with C being a constant of normalization. (Having the freedom to choose this constant makes the lower limit of the integral arbitrary.)

To investigate the spectrum of H_1, we construct the series of Hamiltonia \{H_s\}, where s ranges over the positive integers:

H_s = -\partial_x^2 + V_1(x;a_s) + \sum_{k=1}^{s-1}R(a_k),

in which we have defined

a_s = f^{s-1}(a_1).

Here, f^{s-1}(a_1) denotes applying the transformation f to the parameter set a_1 a total of s-1 times. Each time, we get a new, uniquely determined set of parameters. With this definition, it is easy to see that H_s and H_{s+1} are SUSY partners:

H_{s+1} = -\partial_x^2 + V_1(x;a_{s+1}) + \sum_{k=1}^{s}R(a_k),

which is also equal to
 -\partial_x^2 + V_2(x;a_s) + \sum_{k=1}^{s-1}R(a_k).

Because H_s and H_{s+1} are SUSY partners, they share the same eigenenergies, save for the ground state, whose energy is just given by the sum of the residuals,

E_0^{(s)} = \sum_{k=1}^{s-1}R(a_k).

The SUSY between H_1 and H_2 means that they are isospectral; the n=1,2,3,\ldots eigenstates of H_1 have the same energies as the n = 0,1,2\ldots states of H_2. The same holds true for H_2 and H_3, H_3 and H_4, H_{17} with H_{18} and so on until we run out of bound states. (For a problem like the hydrogen atom, we never do: n ranges upward to infinity.) We “lose” one state, the lowest-energy ground state, every time we go from H_s to H_{s+1}. Following the degeneracies through and liberally applying the transitive law, we see that the n-th eigenstate of H_1 is degenerate with the ground state of H_n. Therefore, we can recover the entire spectrum of H_1 from the equation above.

We can get the eigenfunctions as well as the eigenvalues by this approach. Up to a normalization constant,

\psi_n^{(s+1)} \propto A_s \psi_{n+1}^{(s)},

so knowing the eigenfunctions of H_1 would allow us to compute those of H_2. The proportionality works in reverse, too:

\psi_n^{(s)} \propto A^\dag_{s+1}\psi_{n-1}^{(s+1)}.

Therefore, if we knew the ground-state wavefunction of H_n, we could apply a sequence of operators A^\dag_n,A^\dag_{n-1},\ldots until we reached a state in the spectrum of H_1. But, we do know the ground state of H_n: it is that state which is annihilated by the operator A_n. Remember that A_n involves a single first-order derivative \partial_x, so solving A_n \psi = 0 only involves solving a first-order differential equation. In this way, for any H_1 which has a shape-invariant potential (SIP), we can exactly solve for its eigenvalues and eigenfunctions.

As we noted earlier, this is reminiscent of the operator-based solution of the harmonic oscillator. With that in mind, the interested reader can for an exercise see what happens when the superpotential takes the simple form

W(x) = x.

READING

  • Fred Cooper, Avinash Khare, Uday Sukhatme, “Supersymmetry and Quantum Mechanics” Phys.Rept. 251 (1995): 267–385. Available as arXiv:hep-th/9405029.

SUSY QM SERIES

[Slashdot] [Digg] [Reddit] [del.icio.us] [Facebook] [Technorati] [Google] [StumbleUpon]

5 Responses to “SUSY QM 2: Shape Invariance”

  1. Great post. Thanks for the arXiv reference

  2. “Hamiltonia?” Really? That sounds like a pluralization of the pseudo-Latin “Hamiltonium.”
    As crappy as it sounds, I’d go with Hamiltonians. Though I could certainly be mistaken on the proper conjugation.

  3. Something in me had been wanting to say that for years, and I finally succumbed. With luck, it shall not happen again.

  4. You know, I actually followed that just fine without flinching until the first integral. The graphics help, what can I say? Anyway, I don’t know that I would ever need the math for a high school physics class–which I haven’t taught in five years, anyway: enrollment is down, sadly–but I enjoy the challenge of trying to follow it all.

  5. I saw this material in the first term of my junior year “at university,” so I doubt that it would be needed in high-school physics; we’d have to fix a whole host of broken things first, before we had to care about this! However, given how much time was wasted during my public-school education (pre-calculus? what in blazes was up with that?) it’s not impossibly hard to imagine speeding things up by a couple years.

    I’m sorry to hear that enrollment in physics classes is down, which I believe you mentioned before. We seem to be so adept at giving the next generation a poor start in life; if only there were prizes for it.