# Quantum harmonic oscillator

Quantum harmonic oscillator

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 an arbitrary potential can be approximated as a harmonic potential at the vicinity of a stable equilibrium point. Furthermore, it is one of the few quantum mechanical systems for which a simple exact solution is known.

One-dimensional harmonic oscillator

Hamiltonian and energy eigenstates

In the one-dimensional harmonic oscillator problem, a particle of mass $displaystyle m$ is subject to a potential $displaystyle V\left(x\right) = frac\left\{1\right\}\left\{2\right\} m omega^2 x^2$. In classical mechanics, $displaystyle m omega^2 = k$is called the spring stiffness coefficient, force constant or spring constant, and $displaystyle omega$the circular frequency.

The Hamiltonian of the particle is:

:$H = frac\left\{p^2\right\}\left\{2m\right\} + frac\left\{1\right\}\left\{2\right\} m omega^2 x^2$

where "x" is the position operator, and "p" is the momentum operator $left\left(p = -i hbar \left\{d over dx\right\} ight\right)$. 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 left| psi ight angle = E left| psi ight angle$.

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

:$leftlangle x | psi_n ight angle = sqrt\left\{frac\left\{1\right\}\left\{2^n,n! cdot left\left(frac\left\{momega\right\}\left\{pi hbar\right\} ight\right)^\left\{1/4\right\} cdot expleft\left(- frac\left\{momega x^2\right\}\left\{2 hbar\right\} ight\right) cdot H_nleft\left(sqrt\left\{frac\left\{momega\right\}\left\{hbar x ight\right)$:$n = 0, 1, 2, ldots$

The first eight solutions ("n" = 0 to 7) are shown on the right. The functions $H_n$ are the Hermite polynomials:

:$H_n\left(x\right)=\left(-1\right)^n e^\left\{x^2\right\}frac\left\{d^n\right\}\left\{dx^n\right\}e^\left\{-x^2\right\}$

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

:$E_n = hbar omega left\left(n + \left\{1over 2\right\} ight\right)$.

This energy spectrum is noteworthy for three reasons. Firstly, the energies are "quantized", and may only take the discrete values of $hbaromega$ 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 $hbaromega/2$, which is called the "ground state energy" or zero-point energy. In the ground state, according to quantum mechanics, an oscillator performs null oscillations and its average kinetic energy is positive. 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. The final reason is that the energy levels are equally spaced, unlike the Bohr model or the particle in a box.

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.

The spectral method 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"

:

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

The operator "a" and "a" have properties as below:

:

We can also define a number operator N which has the following property:

:$N = a^\left\{dagger\right\}a$:$Nleft| psi _n ight angle =nleft| psi _n ight angle$

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

:hbar}momega over 2} left( a^{dagger}-a ight)end{matrix}

The "x" and "p" operators obey the following identity, known as the canonical commutation relation:

:$left \left[x , p ight\right] = ihbar$.

The square brackets in this equation are a commonly-used notational device, known as the commutator, defined as

:$left \left[A , B ight\right] stackrel\left\{mathrm\left\{def${=} AB - BA.

Using the above, we can prove the identities

:$H = hbar omega left\left(a^\left\{dagger\right\}a + 1/2 ight\right)$:$left \left[a , a^\left\{dagger\right\} ight\right] = 1$.

Now, let $left|psi_E ight angle$ denote an energy eigenstate with energy "E". The inner product of any ket with itself must be non-negative, so

:$left\left(a left|psi_E ight angle, a left|psi_E ight angle ight\right) = leftlanglepsi_E ight| a^dagger a left| psi_E ight angle ge 0$.

Expressing "a""a" in terms of the Hamiltonian:

:$leftlanglepsi_E left| \left\{H over hbar omega\right\} - \left\{1 over 2\right\} ight|psi_E ight angle = left\left(\left\{E over hbar omega\right\} - \left\{1 over 2\right\} ight\right) ge 0$,

so that $E ge hbar omega / 2$. Note that when ($a left| psi_E ight angle$) is the zero ket (i.e. a ket with length zero), the inequality is saturated, so that $E = hbar omega / 2$. It is straightforward to check that there exists a state satisfying this condition; it is the ground ("n" = 0) state given in the preceding section.

Using the above identities, we can now show that the commutation relations of "a" and "a" with "H" are:

:.

Thus, provided ($a left| psi_E ight angle$) is not the zero ket,

:.

Similarly, we can show that

:$H \left(a^dagger left| psi_E ight angle\right) = \left(E + hbaromega\right) \left(a^dagger left| psi_E ight angle\right)$.

In other words, "a" acts on an eigenstate of energy "E" to produce, up to a multiplicative constant, another eigenstate of energy $E - hbar omega$, and "a" acts on an eigenstate of energy "E" to produce an eigenstate of energy $E + hbar omega$. For this reason, "a" is called a "lowering operator", and "a" a "raising operator". The two operators together are called ladder operators. In quantum field theory, "a" and "a" are alternatively called "annihilation" and "creation" operators because they destroy and create particles, which correspond to our quanta of energy.

Given any energy eigenstate, we can act on it with the lowering operator, "a", to produce another eigenstate with $hbar omega$-less energy. By repeated application of the lowering operator, it seems that we can produce energy eigenstates down to "E" = −∞. However, this would contradict our earlier requirement that $E ge hbar omega / 2$. Therefore, there must be a ground-state energy eigenstate, which we label $left| 0 ight angle$ (not to be confused with the zero ket), such that

:$a left| 0 ight angle = 0 hbox\left\{\left(zero ket\right)\right\}$.

In this case, subsequent applications of the lowering operator will just produce zero kets, instead of additional energy eigenstates. Furthermore, we have shown above that

:$H left|0 ight angle = \left(hbaromega/2\right) left|0 ight angle$

Finally, by acting on $left| 0 ight angle$ with the raising operator and multiplying by suitable normalization factors, we can produce an infinite set of energy eigenstates $left\left\{left| 0 ight angle, left| 1 ight angle, left| 2 ight angle, ... , left| n ight angle, ... ight\right\}$, such that

:$H left|n ight angle = hbaromega \left(n + 1/2\right) left|n ight angle$

which matches the energy spectrum which we gave in the preceding section.

This method can also be used to quickly find the ground state wave function of the quantum harmonic oscillator.Indeed $a left| 0 ight angle = 0$ becomes

:$xpsi_0\left(x\right) + frac\left\{hslash\right\}\left\{m omega\right\} frac\left\{d psi_0\right\}\left\{dx\right\}\left(x\right) = 0$

so that

:$frac\left\{dpsi_0\left(x\right)\right\}\left\{psi_0\left(x\right)\right\}=-frac\left\{momega\right\}\left\{hslash\right\}x;dx Rightarrow lnpsi_0\left(x\right)= -frac\left\{momega\right\}\left\{2hslash\right\}x^2 + ext\left\{ Const.\right\}$

After normalization this leads to the following position space representation of the ground state wave function.

:$psi_0\left(x\right)= left\left(\left\{momega over pihbar\right\} ight\right)^\left\{1 over 4\right\}e^\left\{-\left\{momega over 2hbar\right\}x^2\right\}.$

Natural length and energy scales

The quantum harmonic oscillator possesses natural scales for length and energy, which can be used to simplify the problem. These can be found by nondimensionalization. The result is that if we measure energy in units of $hbar omega$ and distance in units of $left\left(hbar / left\left(m omega ight\right) ight\right)^\left\{1/2\right\}$, then the Schrödinger equation becomes:

:$H = - \left\{1over2\right\} \left\{d^2 over du^2 \right\} + \left\{1 over 2\right\} u^2$,

and the energy eigenfunctions and eigenvalues become

:$leftlangle x | psi_n ight angle = \left\{1 over sqrt\left\{2^n n! pi^\left\{-1/4\right\} hbox\left\{exp\right\} \left(-u^2 / 2\right) H_n\left(u\right)$

:$E_n = n + \left\{1over 2\right\}$

where $H_n\left(u\right)$ are the Hermite polynomials.

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

Example: diatomic molecules

In diatomic molecules, the natural frequency can be found by: [http://hyperphysics.phy-astr.gsu.edu/hbase/quantum/hosc.html] :$omega = sqrt\left\{frac\left\{k\right\}\left\{m_r$where:$omega = 2 pi f$ is the angular frequency,:"k" is the bond force constant, and:$m_r$ is the reduced mass.

"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"1, ..., "x""N". Corresponding to each position coordinate is a momentum; we label these "p"1, ..., "p""N". The canonical commutation relations between these operators are

:.

The Hamiltonian for this system is

:$H = sum_\left\{i=1\right\}^N left\left( \left\{p_i^2 over 2m\right\} + \left\{1over 2\right\} m omega^2 x_i^2 ight\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"1, ..., "x""N" would refer to the positions of each of the "N" particles. This is a happy property of the "r"2 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\left\{x\right\}|psi_\left\{\left\{n angle=prod_\left\{i=1\right\}^Nlangle x_i|psi_\left\{n_i\right\} angle$

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

:.

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 Hamiltonian is:$H = hbar omega , sum_\left\{i=1\right\}^N left\left(a_i^dagger ,a_i + frac\left\{1\right\}\left\{2\right\} ight\right).$This Hamiltonian is invariant under the dynamic symmetry group "U(N)" (the unitary group in "N" dimensions), defined by:$U, a_i^dagger ,U^dagger = sum_\left\{j=1\right\}^N a_j^dagger,U_\left\{ji\right\}quadhbox\left\{for all\right\}quadU in U\left(N\right),$where $U_\left\{ji\right\}$ is an element in the defining matrix representation of "U(N)".

The energy levels of the system are

:$E = hbar omega left \left[\left(n_1 + cdots + n_N\right) + \left\{Nover 2\right\} ight\right]$.:$n_i = 0, 1, 2, dots quad \left(hbox\left\{the number of bosons in mode \right\} i\right).$

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"1 + "n"2 + "n"3. All states with the same "n" will have the same energy. For a given "n", we choose a particular "n"1. Then "n"2 + "n"3 = "n" − "n"1. There are "n" − "n"1 + 1 possible groups {"n"2, "n"3}. "n"2 can take on the values 0 to "n" − "n"1, and for each "n"2 the value of "n"3 is fixed. The degree of degeneracy therefore is:

:$g_n = sum_\left\{n_1=0\right\}^n n - n_1 + 1 = frac\left\{\left(n+1\right)\left(n+2\right)\right\}\left\{2\right\}$Formula for general "N" and "n" ["g"n being the dimension of the symmetric irreducible "n"th power representation of the unitary group "U(N)"] ::The special case "N = 3", given above, follows directly from this general equation.

Example: 3D isotropic harmonic oscillator

The Schrödinger equation of a spherically-symmetric three-dimensional harmonic oscillator can be solved explicitly by separation of variables, see this article for the present case. This procedure is analogous to the separation performed in the hydrogen-like atom problem, but with the spherically symmetric potential:$V\left(r\right) = \left\{1over 2\right\} mu omega^2 r^2,$where $mu$ is the mass of the problem. Because "m" will be used below for the magnetic quantum number, mass is indicated by $mu$, instead of "m", as earlier in this article.

The solution reads:$psi_\left\{klm\right\}\left(r, heta,phi\right) = N_\left\{kl\right\} r^\left\{l\right\}e^\left\{- u r^2\right\}\left\{L_k\right\}^\left\{\left(l+\left\{1over 2\right\}\right)\right\}\left(2 u r^2\right) Y_\left\{lm\right\}\left( heta,phi\right)$where:$N_\left\{kl\right\}=sqrt\left\{sqrt\left\{frac\left\{2 u ^\left\{3\left\{pi frac\left\{2^\left\{k+2l+3\right\};k!; u ^\left\{l\left\{\left(2k+2l+1\right)!!,$ is a normalization constant.:$u equiv \left\{mu omega over 2 hbar\right\}$:$\left\{L_k\right\}^\left\{\left(l+\left\{1over 2\right\}\right)\right\}\left(2 u r^2\right)$ are generalized Laguerre polynomials. The order "k" of the polynomial is a non-negative integer.:$Y_\left\{lm\right\}\left( heta,phi\right),$ is a spherical harmonic function.:$hbar$ is the reduced Planck constant: $hbarequivfrac\left\{h\right\}\left\{2pi\right\}$.

The energy eigenvalue is:$E=hbar omega \left(2k+l+\left\{3over 2\right\}\right)$The energy is usually described by the single quantum number:$nequiv 2k+l,$

Because "k" is a non-negative integer, for every even "n" we have $l=0,2...,n-2,n$ and for every odd "n" we have $l=1,3...,n-2,n$. The magnetic quantum number "m" is an integer satisfying $-l le m le l$, so for every "n" and "l" there are "2l+1" different quantum states, labeled by "m". Thus, the degeneracy at level "n" is:$sum_\left\{l=ldots,n-2,n\right\} \left(2l+1\right) = \left\{\left(n+1\right)\left(n+2\right)over 2\right\}$where the sum starts from 0 or 1, according to whether "n" is even or odd.This result is in accordance with the dimension formula above.

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"1, "x"2, ..., 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 "x"s are vector quantities. The Hamiltonian of the total system is

:$H = sum_\left\{i=1\right\}^N \left\{p_i^2 over 2m\right\} + \left\{1over 2\right\} m omega^2sum_\left\{\left\{ij\right\} \left(nn\right)\right\} \left(x_i - x_j\right)^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.

ee also

*Gas in a harmonic trap
*Creation and annihilation operators
*Coherent state
*Morse potential
*Hooke's atom

References

*cite book | author=Griffiths, David J. | title=Introduction to Quantum Mechanics (2nd ed.) | publisher=Prentice Hall |year=2004 |id=ISBN 0-13-805326-X
*cite book | author=Liboff, Richard L. | title=Introductory Quantum Mechanics | publisher=Addison-Wesley | year=2002 | id=ISBN 0-8053-8714-5

* [http://hyperphysics.phy-astr.gsu.edu/hbase/quantum/hosc.html Quantum Harmonic Oscillator]
*Calculation using a noncommutative free monoid: [http://www.schmarsow.net/oscillate.pdf (mathematical version)] / [http://www.schmarsow.net/oscillateSmall.pdf (abbreviated version)]

Wikimedia Foundation. 2010.

### Look at other dictionaries:

• Quantum mechanics — For a generally accessible and less technical introduction to the topic, see Introduction to quantum mechanics. Quantum mechanics …   Wikipedia

• Quantum field theory — In quantum field theory (QFT) the forces between particles are mediated by other particles. For instance, the electromagnetic force between two electrons is caused by an exchange of photons. But quantum field theory applies to all fundamental… …   Wikipedia

• Quantum state — In quantum physics, a quantum state is a mathematical object that fully describes a quantum system. One typically imagines some experimental apparatus and procedure which prepares this quantum state; the mathematical object then reflects the… …   Wikipedia

• Quantum wire — In condensed matter physics, a quantum wire is an electrically conducting wire, in which quantum effects are affecting transport properties. Due to the confinement of conduction electrons in the transverse direction of the wire, their transverse… …   Wikipedia

• Quantum vibration — A quantum vibration is a vibration of a chemical bond in a molecule that must be treated quantum mechanically. The low lying vibration energy states can be described as states of the quantum harmonic oscillator, and at higher vibrational states,… …   Wikipedia

• Mathematical formulation of quantum mechanics — Quantum mechanics Uncertainty principle …   Wikipedia

• P-adic quantum mechanics — One may compute the energy levels for a potential well like this one.[note 1] P adic quantum mechanics is a relatively recent approach to understanding the nature of fundamental physics. It is the application of p adic analysis to quantum… …   Wikipedia

• Gas in a harmonic trap — The results of the quantum harmonic oscillator can be used to look at the equilibrium situation for a quantum ideal gas in a harmonic trap, which is a harmonic potential containing a large number particles that do not interact with each other… …   Wikipedia

• Perturbation theory (quantum mechanics) — In quantum mechanics, perturbation theory is a set of approximation schemes directly related to mathematical perturbation for describing a complicated quantum system in terms of a simpler one. The idea is to start with a simple system for which a …   Wikipedia