The quantum harmonic oscillator is the quantum mechanical analogue of the classical harmonic oscillator. It is one of the most important model systems in quantum mechanics because, as in classical mechanics, a wide variety of physical situations can be reduced to it either exactly or approximately. In particular, a system near an equilibrium configuration can often be described in terms of one or more harmonic oscillators. Furthermore, it is one of the few quantum mechanical systems for which a simple exact solution is known.

The following discussion of the quantum harmonic oscillator relies on the article mathematical formulation of quantum mechanics.

One-dimensional harmonic oscillator

---

Diatomic molecules

---

In diatomic molecules, the natural frequency can be found by: ^*

$\backslash omega\; =\; \backslash sqrt\{\backslash frac\{k\}\{m\_r\}\}$

where

$\backslash omega\; =\; 2\; \backslash pi\; f$   is the angular frequency,

k is the bond force constant, and

$m\_r$ is the reduced mass.

Hamiltonian and energy eigenstates

In the one-dimensional harmonic oscillator problem, a particle of mass m is subject to a potential V(x) = (1/2)mω² x². The Hamiltonian of the particle is:

$H\; =\; \backslash frac\{p^2\}\{2m\}\; +\; \backslash frac\{1\}\{2\}\; m\; \backslash omega^2\; x^2$

where x is the position operator, and p is the momentum operator ($p\; =\; -i\; \backslash hbar\; \{\backslash partial\; \backslash over\; \backslash partial\; x\}$). The first term represents the kinetic energy of the particle, and the second term represents the potential energy in which it resides. In order to find the energy levels and the corresponding energy eigenstates, we must solve the time-independent Schrödinger equation,

$H\; \backslash left|\; \backslash psi\; \backslash right\backslash rangle\; =\; E\; \backslash left|\; \backslash psi\; \backslash right\backslash rangle$.

We can solve the differential equation in the coordinate basis, using a power series method. It turns out that there is a family of solutions,

$\backslash left\backslash langle\; x\; |\; \backslash psi\_n\; \backslash right\backslash rangle\; =\; \backslash frac\{1\}\{\backslash sqrt\{2^n\; n!\}\}\; \backslash cdot\; \backslash left(\backslash frac\{m\backslash omega\}\{\backslash pi\; \backslash hbar\}\backslash right)^\{1/4\}\; \backslash cdot\; \backslash exp$

\left(- \frac{m\omega x^2}{2 \hbar} \right) \cdot H_n\left(\sqrt{\frac{m\omega}{\hbar}} x \right)

$n\; =\; 0,\; 1,\; 2,\; \backslash ldots$

The first six solutions (n = 0 to 5) are shown on the right. The functions $H\_n$ are the Hermite polynomials:

$H\_n(x)=(-1)^n\; e^\{x^2\}\backslash frac\{d^n\}\{dx^n\}e^\{-x^2\}$

They should not be confused with the Hamiltonian, which is also denoted by H. The corresponding energy levels are

$E\_n\; =\; \backslash hbar\; \backslash omega\; \backslash left(n\; +\; \{1\backslash over\; 2\}\backslash right)$.

This energy spectrum is noteworthy for two reasons. Firstly, the energies are "quantized", and may only take the discrete values of $\backslash hbar\backslash omega$ times 1/2, 3/2, 5/2, and so forth. This is a feature of many quantum mechanical systems. In the following section on ladder operators, we will engage in a more detailed examination of this phenomenon. Secondly, the lowest achievable energy is not zero, but $\backslash hbar\backslash omega/2$, which is called the "ground state energy" or zero-point energy. It is not obvious that this is significant, because normally the zero of energy is not a physically meaningful quantity, only differences in energies. Nevertheless, the ground state energy has many implications, particularly in quantum gravity.

Note that the ground state probability density is concentrated at the origin. This means the particle spends most of its time at the bottom of the potential well, as we would expect for a state with little energy. As the energy increases, the probability density becomes concentrated at the "classical turning points", where the state's energy coincides with the potential energy. This is consistent with the classical harmonic oscillator, in which the particle spends most of its time (and is therefore most likely to be found) at the turning points, where it is the slowest. The correspondence principle is thus satisfied.

Ladder operator method

The power series solution, though straightforward, is rather tedious. The "ladder operator" method, due to Paul Dirac, allows us to extract the energy eigenvalues without directly solving the differential equation. Furthermore, it is readily generalizable to more complicated problems, notably in quantum field theory. Following this approach, we define the operators a and its adjoint a^†

$\backslash begin\{matrix\}$

a &=& \sqrt{m\omega \over 2\hbar} \left(x + {i \over m \omega} p \right) \\ a^{\dagger} &=& \sqrt{m \omega \over 2\hbar} \left( x - {i \over m \omega} p \right) \end{matrix}

The operator a is not Hermitian since it and its adjoint a^† are not equal.

In deriving the form of a^†, we have used the fact that the operators x and p, which represent observables, are Hermitian. These observable operators can be expressed as a linear combination of the ladder operators as

$\backslash begin\{matrix\}$

x &=& \sqrt{\hbar \over 2m\omega} \left( a^{\dagger} + a \right) \\ p &=& i \sqrt \pi^{-1/4} \hbox{exp} (-u^2 / 2) H_n(u)

$E\_n\; =\; n\; +\; \{1\backslash over\; 2\}$.

To avoid confusion, we will not adopt these natural units in this article. However, they frequently come in handy when performing calculations.

N-dimensional harmonic oscillator

---

The one-dimensional harmonic oscillator is readily generalizable to N dimensions, where N = 1, 2, 3, ... . In one dimension, the position of the particle was specified by a single coordinate, x. In N dimensions, this is replaced by N position coordinates, which we label x₁, ..., x_N. Corresponding to each position coordinate is a momentum; we label these p₁, ..., p_N. The canonical commutation relations between these operators are

$\backslash begin\{matrix\}$

\left, p_j \right &=& i\hbar\delta_{i,j} \\ \left, x_j \right &=& 0 \\ \left, p_j \right &=& 0 \end{matrix}.

The Hamiltonian for this system is

$H\; =\; \backslash sum\_\{i=1\}^N\; \backslash left(\; \{p\_i^2\; \backslash over\; 2m\}\; +\; \{1\backslash over\; 2\}\; m\; \backslash omega^2\; x\_i^2\; \backslash right)$.

As the form of this Hamiltonian makes clear, the N-dimensional harmonic oscillator is exactly analogous to N independent one-dimensional harmonic oscillators with the same mass and spring constant. In this case, the quantities x₁, ..., x_N would refer to the positions of each of the N particles. This is a happy property of the r² potential, which allows the potential energy to be separated into terms depending on one coordinate each.

This observation makes the solution straightforward. For a particular set of quantum numbers {n} the energy eigenfunctions for the N-dimensional oscillator are expressed in terms of the 1-dimensional eigenfunctions as:

$$

\langle \mathbf{x}|\psi_{\{n\}}\rangle =\prod_{i=1}^N\langle x_i|\psi_{n_i}\rangle

In the ladder operator method, we define N sets of ladder operators,

$\backslash begin\{matrix\}$

a_i &=& \sqrt{m\omega \over 2\hbar} \left(x_i + {i \over m \omega} p_i \right) \\ a^{\dagger}_i &=& \sqrt{m \omega \over 2\hbar} \left( x_i - {i \over m \omega} p_i \right) \end{matrix}.

By a procedure analogous to the one-dimensional case, we can then show that each of the a_i and a^†_i operators lower and raise the energy by ℏω respectively. The energy levels of the system are

$E\; =\; \backslash hbar\; \backslash omega\; \backslash left+\; \backslash cdots\; +\; n\_N)\; +\; \{N\backslash over\; 2\}\backslash right$.

$\backslash ,\; n\_i\; =\; 0,\; 1,\; 2,\; \backslash dots$

As in the one-dimensional case, the energy is quantized. The ground state energy is N times the one-dimensional energy, as we would expect using the analogy to N independent one-dimensional oscillators. There is one further difference: in the one-dimensional case, each energy level corresponds to a unique quantum state. In N-dimensions, except for the ground state, the energy levels are degenerate, meaning there are several states with the same energy.

The degeneracy can be calculated relatively easily, as an example, consider the 3-dimensional case: Define n = n₁ + n₂ + n₃. All states with the same n will have the same energy. For a given n, we choose a particular n₁. Then n₂ + n₃ = n − n₁. There are n − n₁ + 1 possible groups {n₂, n₃}. n₂ can take on the values 0 to n − n₁, and for each n₂ the value of n₃ is fixed. The degree of degeneracy therefore is:

$$

g_n = \sum_{n_1=0}^n n - n_1 + 1 = \sum_{n_1=0}^n n + 1 - \sum_{n_1=0}^n n_1 = (n+1)(n+1) - \frac{n(n+1)}{2} = \frac{(n+1)(n+2)}{2}

Related problems

---

The quantum harmonic oscillator can be extended in many interesting ways. We will briefly discuss two of the more important extensions, the anharmonic oscillator and coupled harmonic oscillators.

Anharmonic oscillator

As mentioned in the introduction, a system residing "near" the minimum of some potential may be treated as a harmonic oscillator. In this approximation, we Taylor-expand the potential energy around the minimum and discard terms of third or higher order, resulting in an approximate quadratic potential. Once we have studied the system in this approximation, we may wish to investigate the corrections due to the discarded higher-order terms, particularly the third-order term.

The anharmonic oscillator Hamiltonian is the harmonic oscillator Hamiltonian with an additional x³ potential:

$H\; =\; \{p^2\; \backslash over\; 2m\}\; +\; \{1\backslash over\; 2\}\; m\; \backslash omega^2\; x^2\; +\; \backslash lambda\; x^3$

If the harmonic approximation is valid, the coefficient λ is small compared to the quadratic term. We may therefore use perturbation theory to determine the corrections to the states and energy levels imposed by the anharmonic term. This task may be simplified by using the ladder operators to rewrite the anharmonic term as

$\backslash lambda\; \backslash left(\{\backslash hbar\; \backslash over\; 2m\backslash omega\}\backslash right)^\{3\backslash over2\}\; (a\; +\; a^\backslash dagger)^3.$

It turns out that the correction to the energies vanish to first-order in λ. The second-order corrections are given by the usual formula in perturbation theory:

$\backslash Delta\; E^\{(2)\}\; =\; \backslash lambda^2\; \backslash left\backslash langle\; \backslash psi\_E\; \backslash right|\; x^3\; \{1\; \backslash over\; E\; -\; H\_0\}\; x^3\; \backslash left|\; \backslash psi\_E\; \backslash right\backslash rangle.$

This is straightforward, though tedious, to evaluate. One failing of this method, however, is that it does not take into account the possibility of the particle tunnelling out, since it is no longer bound on both sides.

Coupled harmonic oscillators

In this problem, we consider N equal masses which are connected to their neighbors by springs, in the limit of large N. The masses form a linear chain in one dimension, or a regular lattice in two or three dimensions.

As in the previous section, we denote the positions of the masses by x₁, x₂, ..., as measured from their equilibrium positions (i.e. x_k = 0 if particle k is at its equilibrium position.) In two or more dimensions, the xs are vector quantities. The Hamiltonian of the total system is

$H\; =\; \backslash sum\_\{i=1\}^N\; \{p\_i^2\; \backslash over\; 2m\}\; +\; \{1\backslash over\; 2\}\; m\; \backslash omega^2\; \backslash sum\_\{\backslash \{ij\backslash \}\; (nn)\}\; (x\_i\; -\; x\_j)^2$

The potential energy is summed over "nearest-neighbor" pairs, so there is one term for each spring.

Remarkably, there exists a coordinate transformation to turn this problem into a set of independent harmonic oscillators, each of which corresponds to a particular collective distortion of the lattice. These distortions display some particle-like properties, and are called phonons. Phonons occur in the ionic lattices of many solids, and are extremely important for understanding many of the phenomena studied in solid state physics.

See also

---

• Gas in a harmonic trap
• Creation and annihilation operators
• Coherent state
• Morse potential

References

---

External links

---

• Quantum Harmonic Oscillator

Quantum models

Oscillateur harmonique quantique | אוסצילטור הרמוני קוונטי