KNOWPIA
WELCOME TO KNOWPIA

The **quantum harmonic oscillator** is the quantum-mechanical analog of the classical harmonic oscillator. Because an arbitrary smooth potential can usually be approximated as a harmonic potential at the vicinity of a stable equilibrium point, it is one of the most important model systems in quantum mechanics. Furthermore, it is one of the few quantum-mechanical systems for which an exact, analytical solution is known.^{[1]}^{[2]}^{[3]}

The Hamiltonian of the particle is:

The time-independent Schrödinger equation is,

Then solve the differential equation representing this eigenvalue problem in the coordinate basis, for the wave function , using a spectral method. It turns out that there is a family of solutions. In this basis, they amount to Hermite functions,

The functions *H _{n}* are the physicists' Hermite polynomials,

The corresponding energy levels are

The variance in both position and momentum are observed to increase for higher energy levels. The lowest energy level has value of which is its minimum value due to uncertainty relation and also corresponds to a gaussian wavefunction.

This energy spectrum is noteworthy for three reasons. First, the energies are quantized, meaning that only discrete energy values (integer-plus-half multiples of *ħω*) are possible; this is a general feature of quantum-mechanical systems when a particle is confined. Second, these discrete energy levels are equally spaced, unlike in the Bohr model of the atom, or the particle in a box. Third, the lowest achievable energy (the energy of the *n* = 0 state, called the ground state) is not equal to the minimum of the potential well, but *ħω*/2 above it; this is called zero-point energy. Because of the zero-point energy, the position and momentum of the oscillator in the ground state are not fixed (as they would be in a classical oscillator), but have a small range of variance, in accordance with the Heisenberg uncertainty principle.

The ground state probability density is concentrated at the origin, which means the particle spends most of its time at the bottom of the potential well, as one would expect for a state with little energy. As the energy increases, the probability density peaks at the classical "turning points", where the state's energy coincides with the potential energy. (See the discussion below of the highly excited states.) This is consistent with the classical harmonic oscillator, in which the particle spends more of its time (and is therefore more likely to be found) near the turning points, where it is moving the slowest. The correspondence principle is thus satisfied. Moreover, special nondispersive wave packets, with minimum uncertainty, called coherent states oscillate very much like classical objects, as illustrated in the figure; they are *not* eigenstates of the Hamiltonian.

The "ladder operator" method, developed by Paul Dirac, allows extraction of the energy eigenvalues without directly solving the differential equation. It is generalizable to more complicated problems, notably in quantum field theory. Following this approach, we define the operators a and its adjoint *a*^{†},

These operators lead to the following representation of and ,

The operator a is not Hermitian, since itself and its adjoint *a*^{†} are not equal. The energy eigenstates |*n*⟩, when operated on by these ladder operators, give

From the relations above, we can also define a number operator N, which has the following property:

The following commutators can be easily obtained by substituting the canonical commutation relation,

and the Hamilton operator can be expressed as

so the eigenstate of N is also the eigenstate of energy.

The commutation property yields

and similarly,

This means that a acts on |*n*⟩ to produce, up to a multiplicative constant, |*n*–1⟩, and *a*^{†} acts on |*n*⟩ to produce |*n*+1⟩. For this reason, a is called an **annihilation operator** ("lowering operator"), and *a*^{†} a **creation operator** ("raising operator"). The two operators together are called ladder operators.

Given any energy eigenstate, we can act on it with the lowering operator, a, to produce another eigenstate with *ħω* less energy. By repeated application of the lowering operator, it seems that we can produce energy eigenstates down to *E* = −∞. However, since

the smallest eigenvalue of the number operator is 0, and

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

Finally, by acting on |0⟩ with the raising operator and multiplying by suitable normalization factors, we can produce an infinite set of energy eigenstates

such that

Arbitrary eigenstates can be expressed in terms of |0⟩,

The preceding analysis is algebraic, using only the commutation relations between the raising and lowering operators. Once the algebraic analysis is complete, one should turn to analytical questions. First, one should find the ground state, that is, the solution of the equation . In the position representation, this is the first-order differential equation

Explicitly connecting with the previous section, the ground state |0⟩ in the position representation is determined by ,

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 *energy* is measured in units of *ħω* and *distance* in units of √*ħ*/(*mω*), then the Hamiltonian simplifies to

To avoid confusion, these "natural units" will mostly not be adopted in this article. However, they frequently come in handy when performing calculations, by bypassing clutter.

For example, the fundamental solution (propagator) of *H* − *i∂ _{t}*, the time-dependent Schrödinger operator for this oscillator, simply boils down to the Mehler kernel,

The coherent states (also known as Glauber states) of the harmonic oscillator are special nondispersive wave packets, with minimum uncertainty *σ _{x}*

The coherent states are indexed by *α* ∈ **C** and expressed in the |*n*⟩ basis as

Since coherent states are not energy eigenstates, their time evolution is not a simple shift in wavefunction phase. The time-evolved states are, however, also coherent states but with phase-shifting parameter α instead: .

Because and via the Kermack-McCrae identity, the last form is equivalent to a unitary displacement operator acting on the ground state: . Calculating the expectation values:

where is the phase contributed by complex α. These equations confirm the oscillating behavior of the particle.

The uncertainty calculated using the numeric method are:

which gives . Since the only wavefunction that can have lowest position-momentum uncertainty, , is a gaussian wavefunction, and since the coherent state wavefunction has minimum position-momentum uncertainty, we note that the general gaussian wavefunction in quantum mechanics has the form:

The probability of each energy eigenstates can be calculated to find the energy distribution of the wavefunction:

which corresponds to Poisson distribution.

When n is large, the eigenstates are localized into the classical allowed region, that is, the region in which a classical particle with energy *E*_{n} can move. The eigenstates are peaked near the turning points: the points at the ends of the classically allowed region where the classical particle changes direction. This phenomenon can be verified through asymptotics of the Hermite polynomials, and also through the WKB approximation.

The frequency of oscillation at x is proportional to the momentum *p*(*x*) of a classical particle of energy *E*_{n} and position x. Furthermore, the square of the amplitude (determining the probability density) is *inversely* proportional to *p*(*x*), reflecting the length of time the classical particle spends near x. The system behavior in a small neighborhood of the turning point does not have a simple classical explanation, but can be modeled using an Airy function. Using properties of the Airy function, one may estimate the probability of finding the particle outside the classically allowed region, to be approximately

In the phase space formulation of quantum mechanics, eigenstates of the quantum harmonic oscillator in several different representations of the quasiprobability distribution can be written in closed form. The most widely used of these is for the Wigner quasiprobability distribution.

The Wigner quasiprobability distribution for the energy eigenstate |*n*⟩ is, in the natural units described above,^{[citation needed]}

Meanwhile, the Husimi Q function of the harmonic oscillator eigenstates have an even simpler form. If we work in the natural units described above, we have

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

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 convenient 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 the energy eigenfunctions for the *N*-dimensional oscillator are expressed in terms of the 1-dimensional eigenfunctions as:

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

By an analogous procedure to the one-dimensional case, we can then show that each of the *a _{i}* and

The energy levels of the system are

As in the one-dimensional case, the energy is quantized. The ground state energy is *N* times the one-dimensional ground 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 pairs {*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:

This arises due to the constraint of putting *N* quanta into a state ket where and , which are the same constraints as in integer partition.

The Schrödinger equation for a particle in 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 a different spherically symmetric potential

The solution to the equation is:^{[7]}

- is a normalization constant; ;

are generalized Laguerre polynomials; The order k of the polynomial is a non-negative integer;

- is a spherical harmonic function;
- ħ is the reduced Planck constant:

The energy eigenvalue is

Because k is a non-negative integer, for every even n we have *ℓ* = 0, 2, …, *n* − 2, *n* and for every odd n we have *ℓ* = 1, 3, …, *n* − 2, *n* . The magnetic quantum number m is an integer satisfying −*ℓ* ≤ *m* ≤ *ℓ*, so for every n and *ℓ* there are 2*ℓ* + 1 different quantum states, labeled by m . Thus, the degeneracy at level n is

The notation of a harmonic oscillator can be extended to a one-dimensional lattice of many particles. Consider a one-dimensional quantum mechanical *harmonic chain* of *N* identical atoms. This is the simplest quantum mechanical model of a lattice, and we will see how phonons arise from it. The formalism that we will develop for this model is readily generalizable to two and 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 _{i}* = 0 if the particle i is at its equilibrium position). In two or more dimensions, the

We introduce, then, a set of N "normal coordinates" *Q _{k}*, defined as the discrete Fourier transforms of the xs, and N "conjugate momenta" Π defined as the Fourier transforms of the ps,

The quantity *k _{n}* will turn out to be the wave number of the phonon, i.e. 2

This preserves the desired commutation relations in either real space or wave vector space

From the general result

The Hamiltonian may be written in wave vector space as

Note that the couplings between the position variables have been transformed away; if the Qs and Πs were hermitian (which they are not), the transformed Hamiltonian would describe N *uncoupled* harmonic oscillators.

The form of the quantization depends on the choice of boundary conditions; for simplicity, we impose *periodic* boundary conditions, defining the (*N* + 1)-th atom as equivalent to the first atom. Physically, this corresponds to joining the chain at its ends. The resulting quantization is

The upper bound to n comes from the minimum wavelength, which is twice the lattice spacing a, as discussed above.

The harmonic oscillator eigenvalues or energy levels for the mode *ω _{k}* are

If we ignore the zero-point energy then the levels are evenly spaced at

So an **exact** amount of energy *ħω*, must be supplied to the harmonic oscillator lattice to push it to the next energy level. In analogy to the photon case when the electromagnetic field is quantised, the quantum of vibrational energy is called a phonon.

All quantum systems show wave-like and particle-like properties. The particle-like properties of the phonon are best understood using the methods of second quantization and operator techniques described elsewhere.^{[9]}

In the continuum limit, *a* → 0, *N* → ∞, while *Na* is held fixed. The canonical coordinates *Q _{k}* devolve to the decoupled momentum modes of a scalar field, , whilst the location index i (

- The vibrations of a diatomic molecule are an example of a two-body version of the quantum harmonic oscillator. In this case, the angular frequency is given by
^{[10]} - The Hooke's atom is a simple model of the helium atom using the quantum harmonic oscillator.
- Modelling phonons, as discussed above.
- A charge with mass in a uniform magnetic field is an example of a one-dimensional quantum harmonic oscillator: Landau quantization.

- Quantum pendulum
- Quantum machine – human-made device whose collective motion follows the laws of quantum mechanics
- Gas in a harmonic trap – Quantum mechanical model
- Creation and annihilation operators – Operators useful in quantum mechanics
- Coherent state – Specific quantum state of a quantum harmonic oscillator
- Morse potential – Model for the potential energy of a diatomic molecule
- Bertrand's theorem – Physics theorem
- Mehler kernel
- Molecular vibration – Periodic motion of the atoms of a molecule

**^**The normalization constant is , and satisfies the normalization condition .

**^**Griffiths, David J. (2004).*Introduction to Quantum Mechanics*(2nd ed.). Prentice Hall. ISBN 978-0-13-805326-0.**^**Liboff, Richard L. (2002).*Introductory Quantum Mechanics*. Addison–Wesley. ISBN 978-0-8053-8714-8.**^**Rashid, Muneer A. (2006). "Transition amplitude for time-dependent linear harmonic oscillator with Linear time-dependent terms added to the Hamiltonian" (PDF-Microsoft PowerPoint).*M.A. Rashid – Center for Advanced Mathematics and Physics*. National Center for Physics. Retrieved 19 October 2010.**^**Hall, Brian C. (2013),*Quantum Theory for Mathematicians*, Graduate Texts in Mathematics, vol. 267, Springer, Theorem 11.4, ISBN 978-1461471158**^**Pauli, W. (2000),*Wave Mechanics: Volume 5 of Pauli Lectures on Physics*(Dover Books on Physics). ISBN 978-0486414621 ; Section 44.**^**Condon, E. U. (1937). "Immersion of the Fourier transform in a continuous group of functional transformations",*Proc. Natl. Acad. Sci. USA***23**, 158–164. online**^**Albert Messiah,*Quantum Mechanics*, 1967, North-Holland, Ch XII, § 15, p 456.online**^**Fradkin, D. M. "Three-dimensional isotropic harmonic oscillator and SU3."*American Journal of Physics***33**(3) (1965) 207–211.**^**Mahan, GD (1981).*Many particle physics*. New York: Springer. ISBN 978-0306463389.**^**"Quantum Harmonic Oscillator".*Hyperphysics*. Retrieved 24 September 2009.

- Quantum Harmonic Oscillator
- Rationale for choosing the ladder operators
- Live 3D intensity plots of quantum harmonic oscillator
- Driven and damped quantum harmonic oscillator (lecture notes of course "quantum optics in electric circuits")