**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, … 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, …

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.

## Comments

The harmonic series has two sums either log infinity or -log0.

## Trackbacks

[…] that is left over. This latter differential equation can be solved by proposing a solution as an infinite series, in which a sum of powers of is plugged into the differential equation and conditions are imposed […]

[…] No Comments Post a comment or leave a trackback: Trackback URL. « Taylor series Harmonic oscillator – series solution » LikeBe the first to like this […]

[…] differential equation shows up during the solution of the Schrödinger equation for the harmonic oscillator. The differential equation can be written […]

[…] can solve this equation using a power series, and we find that in order for the solution to be normalizable, the power series must contain only […]

[…] solution technique to try is that of proposing as a power series (much as we did in solving the harmonic oscillator system in quantum mechanics). That is, we try a solution of the […]

[…] Harmonic oscillator – series solution […]

[…] task here is to solve 1 by using the same method as for the harmonic oscillator. We propose a solution of the […]

[…] asymptotic behaviour of the equation for large and small , factor out this behaviour and then use a series to try to find the solution of what’s […]

[…] In the solution of the Schrödinger equation for the harmonic oscillator, we found that the wave function can be expressed as a power series: […]

[…] last equation is the same as that obtained from the Schrödinger equation earlier (with different variable […]

[…] the operators have the form of the harmonic oscillator hamiltonian, we know that their eigenvalues are . Since , the eigenvalues of must be of form , so since , which must be an […]

[…] has the same form as the harmonic oscillator equation, so it will have the same solution if we define and follow through the same steps. This will give […]

[…] has the same form as the harmonic oscillator equation, so it will have the same solution if we define and follow through the same steps. This will give […]