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 12.2: The Quantum-mechanical Harmonic Oscillator

Restart | Show the energy eigenfunction and spectrum instead

Please wait for the animation to completely load.

The solution of the one-dimensional quantum harmonic oscillator problem begins with the time-independent Schrödinger equation

[−(*ħ*^{2}/2*m*)(*d*^{2}/*dx*^{2}) + (1/2)*m*ω^{2}*x*^{2}] ψ(*x*) = *E*ψ(*x*) , (12.4)

where *k*/2 = *m*ω^{2}/2. We can put this equation into a more standard form as

[*d*^{2}/*dx*^{2} − *m*^{2}ω^{2}*x*^{2}/*ħ*^{2 }+ 2*mE*/*ħ*^{2 }] ψ(*x*) = 0 . (12.5)

We can further simplify the differential equation by defining the following: β^{2 }= *m*ω/*ħ* and ξ^{2 }= β^{2}*x*^{2}, which when substituted into the differential equation, gives

[*d*^{2}/*d*ξ^{2} − ξ^{2 }+ 2*E*/(*ħ*ω)] ψ(ξ) = 0 . (12.6)

Even though this differential equation looks nothing like our original differential equation, the three terms are just the kinetic, potential, and total energies, respectively. Consider the two special limiting cases of Eq. (12.6):

**Case I:** −ξ^{2 }+ 2*E*/(*ħ*ω) ≈ −ξ^{2}, this situation results when the potential energy at a given position in the well is much greater than the total energy, *V * >>* E*. This is forbidden classically. We explicitly solve the differential equation

[*d*^{2}/*d*ξ^{2} − ξ^{2}] ψ(ξ) = 0 ,

which for large ξ yields the solutions ψ(ξ) = exp(±ξ^{2}/2).

**Case II:** −ξ^{2 }+ 2*E*/(*ħ*ω) ≈ 2*E*/(*ħ*ω), this corresponds to the case where the total energy is much greater than the potential energy, *E *>>* V*, which occurs near *x* = 0. We explicitly solve the differential equation

[*d*^{2}/*d*ξ^{2} +2*E*/(*ħ*ω)]ψ(ξ) = 0 ,

which yields the *real* solutions ψ(ξ) =* A*cos(*K*ξ) +* B*sin(*K*ξ) where *K*^{2 }= 2*E*/(*ħ*ω).

What does the total solution look like? In the region where E>>V we have an oscillating solution. However, in the classically forbidden region *V* >> *E*, we have some *leakage *of the energy eigenfunction into the classically-forbidden region. The well-behaved solution in this region is exp(−ξ^{2}/2). as it goes to zero for large ξ → ±∞.

Now that we have an idea of what the bound states should look like, we can find the entire solution to to Eq. (12.6). This solution can be written as Hermite polynomials weighted by a factor of exp(−ξ^{2}/2) = exp(−β^{2}*x*^{2}/2). These solutions can be written in terms of β and *x* as:

ψ_{n}(*x*) = *A*_{n}*H*_{n}(β*x*) exp(−β^{2}*x*^{2}/2) , (12.7)

where *n* = 0,1,2,3,4,…. The term *H*_{n}(β*x*) refers to the Hermite polynomials of order *n* and the *A*_{n} is the normalization factor which is equal to

*A*_{n }= ((*m*ω)^{1/2}/(2^{n}*n*!(*ħ*π)^{1/2})^{1/2} . (12.8)

In terms of the argument ξ, the first 4 Hermite polynomials are

*H*_{0}(ξ) = 1 *H*_{1}(ξ) = 2ξ *H*_{2}(ξ) = 4ξ^{2} − 2 *H*_{3}(ξ) = 8ξ^{3} − 12ξ .

In the animations, the energy eigenfunctions for a quantum harmonic oscillator are shown. **The animation uses ħ = 2m = 1 and ω = 2. **Note that in the well, the energy eigenfunction's amplitude and curviness change with position. These changes are due to the fact that the potential energy function itself changes with position. The energy eigenfunction is

*curvier*and has a smaller amplitude nearer the center of the well as compared to the positions closer to the classical turning point.

By brute force integration we can use these energy eigenfunctions to calculate the expectation values.^{2} We find that <*x*>_{n }=
0 and <*p*>_{n }= 0, as expected, and

<*x*^{2}>_{n }= *(ħ*/2*m*ω) (2*n *+ 1) (12.9)

and

<*p*^{2}>_{n }= *(mħ*ω/2) (2*n *+ 1) , (12.10)

and therefore

<*E*>_{n }= <*p*^{2}>_{n}/2*m* + *m*ω^{2}<*x*^{2}>_{n}/2 = (*n* + 1/2)*ħ*ω . (12.11)

This uniformity in the spacing of the harmonic oscillator energies is shown in the "spectrum" animation. To see the other bound states simply click-drag in the energy level diagram on the left to select a level. The selected level will turn red.** **Since we have chosen ω = 2 and *ħ* = 2*m* = 1, the energy spectrum is just E_{n }= (2*n* + 1).

^{2}There is an easier way. This method uses operators, called raising and lowering operators, to write the harmonic oscillator Hamiltonian. This method was pioneered by Dirac and is the easiest way to calculate expectation values.

« previous

next »