Required math: calculus
Required physics: harmonic oscillator
Reference: Griffiths, David J. (2005), Introduction to Quantum Mechanics, 2nd Edition; Pearson Education – Sec 2.3.
The complete Schrödinger equation for the harmonic oscillator potential is
To solve this equation, we split the wave function into two factors: the first factor is the asymptotic behaviour for large , and the second is a function which we have yet to find. To simplify the notation we introduced two auxiliary variables.
The independent variable is related to the spatial variable by
and the parameter is related to the energy by
The parameter is the frequency of the oscillator.
After analyzing the asymptotic behaviour of the Schrödinger equation for the harmonic oscillator, we write the wave function in the form
Substituting this back into the Schrödinger equation gives us a differential equation for :
If we hurl this equation into mathematical software like Maple, it tells us that the solution involves two forms of Kummer functions, otherwise known as confluent hypergeometric functions of the first and second kinds. Apart from being able to impress your friends in the pub, these terms don’t really help us learn much about the physics. For that we need to solve the differential equation using a power series.
The idea is to propose a solution of the form
The theory behind Taylor series in elementary calculus assures us that for any ‘reasonable’ function (that is, pretty well any function found in physics), it is possible to write the function as a power series, so we should be able to find such a solution. At this stage, we can’t guarantee that such a solution will tell us much, but it’s worth a try.
To use the series, we need to calculate its first two derivatives:
The fancy footwork in the last line just relabels the summation index to make it more convenient for the next step, as we’ll see. To convince yourself it is the same series as the line above it, just write out the first 4 or 5 terms in the series and you’ll see it is the same. Notice also that in the first derivative series the first term is zero due to the factor of , so we don’t actually get a term with in it.
We now want to substitute these derivatives back into the differential equation 1 we want to solve. The reason we juggled the summation index in the second derivative is that we want the series in all three terms in the equation to contain terms rather than to some other power. This makes it easier to group together the terms with equal powers of .
Doing the substitution, we get:
From the mathematics of power series expansions, it is known that any given function’s expansion is unique (the proof takes us too far into pure mathematics so we’ll leave it for now). That means that, having decided on a value for , there is one and only one sequence of s that defines the function . So, since can be any value, the only way the above sum can be zero for all values of is if the coefficient of each power of vanishes separately. That is,
This, in turn, gives a recursion relation for the coefficients:
Since we are solving a second order differential equation we would expect to have two arbitrary constants that must be determined by initial conditions and normalization, and we see that since the recursion formula relates every second coefficient, we need to specify both and to be able to generate all the coefficients. If we start off with we get all the even coefficients:
There is a similar sequence of calculations for the odd coefficients starting with .
That’s about as far as we can go without using some external information to put some conditions on the series. As usual, we require the solution (the original solution, that is, ) to be normalizable. We now know that this solution has the form
So in order to be normalizable, the series will have to converge to some function that doesn’t expand to infinity as fast as . Otherwise the series term will kill off the negative exponential, and the overall wave function will not tend to zero as goes to infinity.
This seems like a difficult condition to check, but let’s have a look at the asymptotic behaviour (for large ) of the recursion formula 2.
Thus the ratio of two successive even terms (or two successive odd terms) in the series is
How does this compare with the series for an exponential function?
The Taylor series for is
The ratio of two successive terms from this series is
which for large is essentially the same as relation 3. And this is for only half (either even or odd terms) of the series; the other half will contribute another function of roughly equal size. Thus it looks like the series’ asymptotic behaviour is that of so the overall behaviour of the wave function is , which diverges and is therefore not normalizable.
This looks like a serious problem, but there is in fact a way out: if the series terminates after a finite number of terms, then the behaviour is that of a polynomial rather than an exponential, and multiplying any polynomial by will always give a normalizable function.
So if we can arrange things so that the recursion formula 2 gives for some , then clearly all further terms will be zero. The condition to be satisfied is therefore
where is some integer 0, 1, 2, 3, \ldots Note however, that each choice of , that is, each choice of where the series terminates, gives a different value for the energy. The lowest possible energy for the harmonic oscillator is when , and is and the energies increase at regular intervals of so the energy levels are all equally spaced.
It is more usual to give the energy formula as
with 0, 1, 2, 3, 4, \ldots
One note of caution here. Once we have chosen an energy level, this fixes the value of at which the series terminates. If is even, then the odd series must be zero right from the start, and vice versa. There is no way of getting both the even and odd series to terminate at some intermediate values in the same solution. So if we choose an even value of we must have to remove all the odd terms from the sum, and conversely if we choose to be odd, we must have to remove all the even terms.
The stationary states for the harmonic oscillator are therefore products of polynomials and the exponential factor. The polynomials turn out to be well-studied in mathematics and are known as Hermite polynomials. We will explore their properties in another post.
To summarize the behaviour of the quantum harmonic oscillator, we’ll list a few points.
- The harmonic oscillator potential is parabolic, and goes to infinity at infinite distance, so all states are bound states – there is no energy a particle can have that will allow it to be free.
- The energies are equally spaced, with spacing .
- The lowest energy is the ground state , so a particle always has positive, non-zero energy.
- The stationary states consist of either an even or an odd polynomial function multiplied by which is always even. Thus a stationary state is either an even or an odd function of (and hence of ).
- The polynomial functions in the stationary states are Hermite polynomials.