Abstract
Harmonic motion is one of the most important examples of motion in all of physics. Any vibration with a restoring force equal to Hooke’s law is generally caused by a simple harmonic oscillator. The potential for the harmonic ocillator is the natural solution every potential with small oscillations at the minimum. Almost all potentials in nature have small oscillations at the minimum, including many systems studied in quantum mechanics. Here, harmonic motion plays a fundamental role as a stepping stone in more rigorous applications. The Harmonic Oscillator is characterized by the its Schrödinger Equation. This equation is presented in section 1.1 of this manual. The harmonic oscillator has only discrete energy states as is true of the onedimensional particle in a box problem. The equation for these states is derived in section 1.2. An exact solution to the harmonic oscillator problem is not only possible, but also relatively easy to compute given the proper tools. It is one of the first applications of quantum mechanics taught at an introductory quantum level. Systems with nearly unsolvable equations are often broken down into smaller systems. The solution to this simple system can then be used on them. A firm understanding of the principles governing the harmonic oscillator is prerequisite to any substantial study of quantum mechanics. 
The classical potential for a harmonic oscillator is derivable from Hooke’s law. It is conventionally written:
 (1) 
Where is the natural frequency, k is the spring constant, and m is the mass of the body.
 (2) 
For convenience in this calculation, the potential for the harmonic oscillator is written
 (3) 
Placing this potential in the one dimensional, timeindependent Schrödinger equation, it yeilds:
 (4) 
which equals:
 (5) 
The Equation for the Quantum Harmonic Oscillator is a second order differential equation that can be solved using a power series. In following section, 2.2, the power series method is used to derive the wave function and the eigenenergies for the quantum harmonic oscillator.
The first step in the power series method is to perform a change of variables by introducing the dimensionless variable, y^{1}:
 (6) 
Substituting this new variable into Equation (5) above yeilds:
 (7) 
For very large values of y, the term is negligible in comparison to the y^{2} term. With this fact we can guess that (y) will be as e^{±}. The general solution to the differential equation is:
 (8) 
We must now calculate the derivatives of that will be substitued into the Schrödinger equation.
 (9) 
and
 (10) 
Putting the values from (8) and (10) in equation (7):
 (11) 
Canceling terms:
 (12) 
And dividing through by e^{}, we obtain:
 (13) 
The next step is to solve the second order differential equation (13) above for u(y) so that we can find an exact solution for (y). We begin by using a power series of y as the general solution to equation (13). u(y) then takes the form below:
 (14) 
In order to substitute this solution into equation (13) we have to solve the first and second derivatives of u(y):
 (15) 
and,
 (16) 
Replacing the u(y) values in equation (13) we have:
 (17) 
Simplifying the second term:
 (18) 
Replacing n with (n+2) in the first term. See section 4, Math Moves and Helpful Hints, for a discussion of this substitution.
 (19) 
and simplifying we have:
 (20) 
we can then solve for _{n+2}, because the coefficient for each power of y must equal zero (RHS of equation (20)):
 (21) 
Equation (21) is a series representation of all the expansion coefficients in terms of _{0} for the power series solution to equation (13). For large values of y, n is also very large. The ratio of _{n+1} and _{n} (from formula (21) for the coefficients of the power series expansion above) is very close to . Here we have a problem, because in the limit, grows faster than the exponential term in (y). The series must terminate in order for our solution to have any physical meaning. The best way to terminate the series is to equate the numerator in equation (21) with zero. We then have:
 (22) 
We can now solve for energy, E:
 (23) 
We have found an infinite number of energies for each energy level, n. The formula for the wave function is incomplete, however, because the power series solution is incomplete. The individual eigenstates of the Hamiltonian must be made orthogonal. Hermite polynomials need be incorporated into the final solution in order to do this. What results is a function of the form:
 (24) 
The following is a table for the first 5 eigenenergies and eigenstates for the harmonic oscillator:^{2}:

The Summation Substitution: Why is replacing n with (n+2) in the power series derivation mathematically legal?
Most textbooks do not expand on the rational for this substitution. The substitution makes perfect sense, however when each term of the summation is expanded and the derivatives for each term are taken. The steps below might help with the logic behind this part of the derivation of the solution to the harmonic oscillator equation: We are given:
 (25) 
Rewriting this summation in terms of its expansion:
 (26) 
Then taking the first and second derivatives of the expanded terms, we have:
 (27) 
 (28) 
This expansion shows that the first two terms in the second derivative (equation 27 above) are zero because the coefficients are zero. The summation for the second derivative actually begins with n+2. Hence the substitution in the power series derivation above.
Hermite Polynomials: What are Hermite polynomials?
The Hermite polynomial is defined as the solution to Hermite’s Differential equation. This polynomial is a direct result of solving the quantum harmonic oscillator differential equation. The Hermite’s Differential equation takes the familiar form:
 (29) 
Where n is a real, nonnegative number (n = 0, 1, 2, 3 )
Hermite polynomials form a complete orthogonal set on the interval  to + with respect to the function e^{y2 }. The orthogonality relationship can be shown as such^{4}:
 (30) 
With the orthogonality condition met, piecewise continuous function such as the solution to the quantum harmonic oscillator can be expressed in terms of the equation for Hermite polynomials:
 (31) 
Where
 (32) 
The Hermite Polynomial is graphed below^{5}:
1. The Schrdinger equation for the one dimensional harmonic oscillator is reduced to the following equation for the polynomial u(y):
 (33) 
a. Determine the recursion relation separately for the even and the odd energy states.
b. Derive an equation for energy for both the even and the odd states from the recursion relations above.
c. Find the first three energy values from each of the equations obtained in (b).
Solution
Part a. We begin by considering only the even states. The general solution given in the power series derivation above may be modified so that we include only these states:
 (34) 
Taking the first and second derivative of u(y) we obtain:
 (35) 
and
 (36) 
Putting these values for the derivative back into the differential equation:
 (37) 
simplifying the second term we have:
 (38) 
Replacing n with (n+1) in the first summation, the equation can be rewritten:
 (39) 
simplifying this:
 (40) 
Again, since the coefficient must equal zero, (see RHS of equation), we can solve for _{(n+1)} and obtain the recursion relation for the even states:
 (41) 
A similar method may be used to find the recursion relation for the odd states. This time we use a general solution of the form below so as to attain only odd states:
 (42) 
After filling this solution into the differential equation and following the same steps as shown for the even states, one obtains the recursion relation:
 (43) 
Part b. As in the power series derivation above, we know that the numerator must go to zero. We can therefore, solve for the energy and obtain the equation for even states:
 (44) 
and the odd states:
 (45) 
Part c.
The first three energy values for the even and odd states are listed below.
Even States Odd States
n = 0 E = E =
n = 1 E = E =
n = 2 E = E =
2. Find the energy levels of a particle moving in a potential field of the shape^{3}
Solution
The first step in solving this equation is to look at the boundary conditions. As x 0, the wave function should fall to zero. For x > 0, the wave function satisfies the differential equation for the harmonic oscillator. Since the odd wave functions for the harmonic oscillator tend toward zero as x 0, we can conclude that the equation for the odd states in Problem 1 above is the solution to the problem:
 (48) 
3. A Harmonic Oscillator is in the initial state: (x,0) = _{n}(x), that is, an eigenstate of . What is (x,t)?^{2}
Solution
From Postulate IV of Quantum Mechanics we can show that:
 (49) 
Therefore, the timedependent wave function for the problem here is:
 (50) 
4. Show that in the nth state of the harmonic oscillator:^{2}
 (51) 
 (52) 
Solution
 (53) 
Therefore,
 (54) 
It can be shown that because of orthogonality:
 (55) 
Therefore,
 (56) 
and
 (57) 
A similar case may be shown for the momentum operator
5. Using the uncertainty principles between x and p derived in Problem 3, derive the zeropoint energy:
 (58) 
for a harmonic oscillator with natural frequency _{0}.^{2}
Solution At E_{0}, the kinetic energy of the system equals the potential energy of the system:
 (59) 
Where
 (60) 
From the relationship developed in Problem 3,
 (61) 
and
 (62) 
So
 (63) 
and
 (64) 
Therefore
 (65) 
Or
 (66) 
Then
 (67) 
And
 (68) 
With the relationship between k and in equation (56),
 (69) 
The following list of references proved excellent resources on the topic of the quantum harmonic oscillator.
[1] Dicke, Robert H. & Wittke, James P., Introduction to Quantum Mechanics, 1960, Addison Wesley, San Francisco, CA.
[2] Gasiorowicz, Stephen, Quantum Physics, 1974, John Wiley & Sons, New York.
[3] Liboff, Richard L., Introductory Quantum Mechanics, 2003, Addison Wesley, San Francisco, CA.
[4] Saxon David S., Elementary Quantum Mechanics, 1968, HoldenDay, San Francisco, CA.
[5] Schiff, Leonard I., Quantum Mechanics, 1955, McGrawHill Book Company, New York.
[6] Flügge, S., Practical Quantum Mechanics, 1974, SpringerVerlag, New York.
[7] Gol’dman, I. I. & Krivchenkov, V. D., Problems in Quantum Mechanics, 1961, Pergamon Press, AddisonWesley Publishing Co., Reading, MA.
[8] http://www.efunda.com/math/Hermite/index.cfm
[9] http://mathworld.wolfram.com/HermitePolynomial.html
[10] http://www.sciencedaily.com/encyclopedia/hermite polynomials
[11] http://sps.physics.gatech.edu/Quantum Harmonic Oscillator Lecture.pdf
[12] http://www.sciencedaily.com/encyclopedia/quantum harmonic oscillator
[13] http://scienceworld.wolfram.com/physics/SimpleHarmonicOscillatorQuantumMechanical.html
[14] http://vergil.chemistry.gatech.edu/notes/ho/node3.html
[15] http://www.physics.orst.edu/~rubin/CPbook/manuel.97/Harmos.html
[16] http://www3.adnc.com/~topquark/quantum/qhomain.html
[17] http://www.physics.orst.edu/ rubin/nacphy/CPapplets/QMWAVE/Harmos2/Harmos2.html