Creation and annihilation operators

From Wikipedia, the free encyclopedia
Jump to: navigation, search

Creation and annihilation operators are mathematical operators that have widespread applications in quantum mechanics, notably in the study of quantum harmonic oscillators and many-particle systems.[1] An annihilation operator lowers the number of particles in a given state by one. A creation operator increases the number of particles in a given state by one, and it is the adjoint of the annihilation operator. In many subfields of physics and chemistry, the use of these operators instead of wavefunctions is known as second quantization.

Creation and annihilation operators can act on states of various types of particles. For example, in quantum chemistry and many-body theory the creation and annihilation operators often act on electron states. They can also refer specifically to the ladder operators for the quantum harmonic oscillator. In the latter case, the raising operator is interpreted as a creation operator, adding a quantum of energy to the oscillator system (similarly for the lowering operator). They can be used to represent phonons.

The mathematics for the creation and annihilation operators for bosons is the same as for the ladder operators of the quantum harmonic oscillator.[2] For example, the commutator of the creation and annihilation operators that are associated with the same boson state equals one, while all other commutators vanish. However, for fermions the mathematics is different, involving anticommutators instead of commutators.[3]

Ladder operators for the quantum harmonic oscillator[edit]

In the context of the quantum harmonic oscillator, we reinterpret the ladder operators as creation and annihilation operators, adding or subtracting fixed quanta of energy to the oscillator system. Creation/annihilation operators are different for bosons (integer spin) and fermions (half-integer spin). This is because their wavefunctions have different symmetry properties.

First consider the simpler bosonic case of the phonons of the quantum harmonic oscillator.

Start with the Schrödinger equation for the one-dimensional time independent quantum harmonic oscillator

\left(-\frac{\hbar^2}{2m} \frac{d^2}{d x^2} + \frac{1}{2}m \omega^2 x^2\right) \psi(x) = E \psi(x)

Make a coordinate substitution to nondimensionalize the differential equation

x \ = \  \sqrt{ \frac{\hbar}{m \omega}} q.

and the Schrödinger equation for the oscillator becomes

 \frac{\hbar \omega}{2} \left(-\frac{d^2}{d q^2} + q^2 \right) \psi(q) = E \psi(q).

Note that the quantity  \hbar \omega = h \nu is the same energy as that found for light quanta and that the parenthesis in the Hamiltonian can be written as

 -\frac{d^2}{dq^2} + q^2 = \left(-\frac{d}{dq}+q \right) \left(\frac{d}{dq}+ q \right) + \frac {d}{dq}q - q \frac {d}{dq}

The last two terms can be simplified by considering their effect on an arbitrary differentiable function f(q),

\left(\frac{d}{dq} q- q \frac{d}{dq} \right)f(q) = \frac{d}{dq}(q  f(q)) - q  \frac{df(q)}{dq} = f(q)

which implies,

\frac{d}{dq} q- q \frac{d}{dq}  = 1


  -\frac{d^2}{dq^2} + q^2 = \left(-\frac{d}{dq}+q \right) \left(\frac{d}{dq}+ q \right) + 1

and the Schrödinger equation for the oscillator becomes, with substitution of the above and rearrangement of the factor of 1/2,

 \hbar \omega \left[\frac{1}{\sqrt{2}} \left(-\frac{d}{dq}+q \right)\frac{1}{\sqrt{2}} \left(\frac{d}{dq}+ q \right) + \frac{1}{2} \right] \psi(q) = E \psi(q).

If we define

a^\dagger \ = \  \frac{1}{\sqrt{2}} \left(-\frac{d}{dq} + q\right) as the "creation operator" or the "raising operator" and
 a \ \ = \  \frac{1}{\sqrt{2}} \left(\ \ \ \!\frac{d}{dq} + q\right) as the "annihilation operator" or the "lowering operator"

then the Schrödinger equation for the oscillator becomes

 \hbar \omega \left( a^\dagger a + \frac{1}{2} \right) \psi(q) = E \psi(q)

This is significantly simpler than the original form. Further simplifications of this equation enables one to derive all the properties listed above thus far.

Letting p = - i \frac{d}{dq}, where "p" is the nondimensionalized momentum operator then we have

 [q, p] = i \,


a = \frac{1}{\sqrt{2}}(q + i p) = \frac{1}{\sqrt{2}}\left( q + \frac{d}{dq}\right)
a^\dagger = \frac{1}{\sqrt{2}}(q - i p) = \frac{1}{\sqrt{2}}\left( q - \frac{d}{dq}\right).

Note that these imply that

 [a, a^\dagger ] = \frac{1}{2} [ q + ip , q-i p] = \frac{1}{2} ([q,-ip] + [ip, q]) = \frac{-i}{2} ([q, p] + [q, p]) = 1 .

The operators a and a^\dagger may be contrasted with normal operators, which commute with their adjoints. A normal operator has a representation A= B + i C where B,C are self-adjoint and commute, i.e. BC=CB. By contrast, a has the representation a=q+ip where p,q are self-adjoint but [p,q]=1. Then B and  C have a common set of eigenfunctions (and are simultaneously diagonalizable), whereas p and q famously don't and aren't.

Despite this, we go on. Using the commutation relations given above, the Hamiltonian operator can be expressed as

\hat H = \hbar \omega \left(  a \, a^\dagger - \frac{1}{2}\right) = \hbar \omega \left(  a^\dagger \, a + \frac{1}{2}\right).\qquad\qquad(*)

One can compute the commutation relations between the a and a^\dagger operators and the Hamiltonian:[4]

[\hat H, a ]  = [\hbar \omega( a a^\dagger - \frac{1}{2}),a] = \hbar \omega [ a a^\dagger, a] = \hbar \omega ( a [a^\dagger,a] + [a,a] a^\dagger) =  -\hbar \omega a.
[\hat H, a^\dagger ]  = \hbar \omega \, a^\dagger .

These relations can be used to easily find all the energy eigenstates of the quantum harmonic oscillator. Assuming that \psi_n is an eigenstate of the Hamiltonian \hat H \psi_n = E_n\, \psi_n. Using these commutation relations, it follows that[4]

\hat H\, a\psi_n = (E_n - \hbar \omega)\, a\psi_n .
\hat H\, a^\dagger\psi_n = (E_n + \hbar \omega)\, a^\dagger\psi_n .

This shows that a\psi_n and a^\dagger\psi_n are also eigenstates of the Hamiltonian, with eigenvalues E_n - \hbar \omega and E_n + \hbar \omega respectively. This identifies the operators a and a^\dagger as "lowering" and "raising" operators between the eigenstates. The energy difference between adjacent eigenstates is \Delta E = \hbar \omega.

The ground state can be found by assuming that the lowering operator possesses a nontrivial kernel, a\, \psi_0 = 0 with \psi_0\ne0. Using the formula above for the Hamiltonian, one obtains

0=\hbar \omega\, a^\dagger a \psi_0 = \left(\hat H - \frac{\hbar \omega}{2} \right) \,\psi_0.

so \psi_0 is an eigenfunction of the Hamiltonian. This gives the ground state energy E_0 = \hbar \omega /2. This allows one to identify the energy eigenvalue of any eigenstate \psi_n as[4]

E_n = \left(n + \frac{1}{2}\right)\hbar \omega.

Furthermore, it turns out that the first-mentioned operator in (*), the number operator N=a^\dagger a\,, plays a most important role in applications, while the second one, a \,a^\dagger\,, can simply be replaced by N +1\,. So one simply gets

\hbar\omega \,\left(N+\frac{1}{2}\right)\,\psi (q) =E\,\psi (q).

Explicit eigenfunctions[edit]

The ground state \ \psi_0(q) of the quantum harmonic oscillator can be found by imposing the condition that

 a \ \psi_0(q) = 0.

Written out as a differential equation, the wavefunction satisfies

q \psi_0 + \frac{d\psi_0}{dq} = 0

which has the solution

\psi_0(q) = C \exp\left(-{q^2 \over 2}\right).

The normalization constant C is found to be  1\over \sqrt[4]{\pi}  from \int_{-\infty}^\infty \psi_0^* \psi_0 \,dq = 1,  using the Gaussian integral. Explicit formulas for all the eigenfunctions can now be found by repeated application of  a^\dagger to  \psi_0. This, and further operator formalism, can be found in Glimm and Jaffe, Quantum Physics, pp. 12-20.

Matrix representation[edit]

The matrix expression of the creation and annihilation operators of the quantum harmonic oscillator with respect to the above orthonormal basis is

0 & 0 & 0 & \dots & 0 &\dots \\
\sqrt{1} & 0 & 0 & \dots & 0 & \dots\\
0 & \sqrt{2} & 0 & \dots & 0 & \dots\\
0 & 0 & \sqrt{3} & \dots & 0 & \dots\\
\vdots & \vdots & \vdots & \ddots  & \vdots  & \dots\\
0 & 0 & 0 & 0 & \sqrt{n} &\dots &  \\
\vdots & \vdots & \vdots & \vdots & \vdots  &\ddots \end{pmatrix}
0 & \sqrt{1} & 0 & 0 & \dots & 0 & \dots \\
0 & 0 & \sqrt{2} & 0 & \dots & 0 & \dots \\
0 & 0 & 0 & \sqrt{3} & \dots & 0 & \dots \\
0 & 0 & 0 & 0 & \ddots & \vdots & \dots \\
\vdots & \vdots & \vdots & \vdots & \ddots & \sqrt{n} & \dots \\
0 & 0 & 0 & 0 & \dots & 0 & \ddots \\
\vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix}

These can be obtained via the relationships a^\dagger_{ij} = \langle\psi_i | {a}^\dagger | \psi_j\rangle and a_{ij} = \langle\psi_i | {a} | \psi_j\rangle. The eigenvectors \psi_i are those of the quantum harmonic oscillator, and are sometimes called the "number basis".

Generalized creation and annihilation operators[edit]

Main article: CCR_and_CAR_algebras

The operators derived above are actually a specific instance of a more generalized notion of creation and annihilation operators. The more abstract form of the operators are constructed as follows. Let H be a one-particle Hilbert space (that is, any Hilbert space, viewed as representing the state of a single particle).

The (bosonic) CCR algebra over H is the algebra-with-conjugation-operator (called *) abstractly generated by elements a(f), where f ranges freely over H, subject to the relations

[a(f),a^\dagger(g)]=\langle f|g \rangle,

where we are using bra–ket notation. The map a : fa(f) from H to the bosonic CCR algebra is required to be complex antilinear (this adds more relations). Its adjoint is a(f), and the map fa(f) is complex linear in H. Thus H embeds as a complex vector subspace of its own CCR algebra. In a representation of this algebra, the element a(f) will be realized as an annihilation operator, and a(f) as a creation operator.

In general, the CCR algebra is infinite dimensional. If we take a Banach space completion, it becomes a C* algebra. The CCR algebra over H is closely related to, but not identical to, a Weyl algebra.

For fermions, the (fermionic) CAR algebra over H is constructed similarly, but using anticommutator relations instead, namely

\{a(f),a^\dagger(g)\}=\langle f|g \rangle.

The CAR algebra is finite dimensional only if H is finite dimensional. If we take a Banach space completion (only necessary in the infinite dimensional case), it becomes a C* algebra. The CAR algebra is closely related to, but not identical to, a Clifford algebra.

Physically speaking, a(f) removes (i.e. annihilates) a particle in the state | f\scriptstyle \rangle whereas a(f) creates a particle in the state | f\scriptstyle \rangle .

The free field vacuum state is the state | 0 \scriptstyle \rangle with no particles, characterized by


If | f\scriptstyle \rangle is normalized so that \scriptstyle \langle f | f\scriptstyle \rangle = 1, then N = a(f) a(f) gives the number of particles in the state | f\scriptstyle \rangle .

Creation and annihilation operators for reaction-diffusion equations[edit]

The annihilation and creation operator description has also been useful to analyze classical reaction diffusion equations, such as the situation when a gas of molecules A diffuse and interact on contact, forming an inert product: A + A → ∅ . To see how this kind of reaction can be described by the annihilation and creation operator formalism, consider n_{i} particles at a site i on a 1-d lattice. Each particle moves to the right or left with a certain probability, and each pair of particles at the same site annihilates each other with a certain other probability.

The probability that one particle leaves the site during the short time period dt is proportional to n_{i}dt, let us say a probability \alpha n_{i}dt to hop left and \alpha n_{i}dt to hop right. All n_i particles will stay put with a probability 1-2\alpha n_{i}dt. (Since dt is so short, the probability that two or more will leave during dt is very small and will be ignored.)

We can now describe the occupation of particles on the lattice as a `ket' of the form | ..., n-1, n0, n1, ...\scriptstyle \rangle . It represents the juxtaposition (or conjunction, or tensor product) of the number states ..., | n-1\scriptstyle \rangle , | n0\scriptstyle \rangle , | n1\scriptstyle \rangle , ... located at the individual sites of the lattice. A slight modification[clarification needed] of the annihilation and creation operators is needed so that

a|n\rangle= \sqrt{n} \ |n-1\rangle


a^{\dagger}|n\rangle= \sqrt{n+1}\ |n+1\rangle,

for all n≥0. This modification preserves the commutation relation


Now let ai=a πi, where πi selects the ith component of ψ. That is, ai makes a copy of the state | ni\scriptstyle \rangle in an abstract place and then applies a to it. Then aii a, where ιi inserts an abstract state at the ith site. Thus, for example, the net effect of ai-1ai is to move an eigenstate from the ith to the (i-1)th site while multiplying with the appropriate factor.

This allows us to write the pure diffusive behaviour of the particles as


where the sum is over i.

The reaction term can be deduced by noting that n particles can interact in n(n-1) different ways, so that the probability that a pair annihilates is \lambda n(n-1)dt, yielding a term


where number state n is replaced by number state n-2 at site i at a certain rate. Thus the state evolves by

\partial_{t}|\psi\rangle=-\alpha\sum(a_{i}^{\dagger}-a_{i-1}^{\dagger})(a_{i}-a_{i-1})|\psi\rangle+\lambda\sum(a_{i}^{2}-a_{i}^{\dagger 2}a_{i}^{2})|\psi\rangle

Other kinds of interactions can be included in a similar manner.

This kind of notation allows the use of quantum field theoretic techniques to be used in the analysis of reaction diffusion systems.

Creation and annihilation operators in quantum field theories[edit]

In quantum field theories and many-body problems one works with creation and annihilation operators of quantum states, a^\dagger_i and a^{\,}_i. These operators change the eigenvalues of the number operator,

N = \sum_i n_i = \sum_i a^\dagger_i a^{\,}_i,

by one, in analogy to the harmonic oscillator. The indices (such as i) represent quantum numbers that label the single-particle states of the system; hence, they are not necessarily single numbers. For example, a tuple of quantum numbers (n, l, m, s) is used to label states in the hydrogen atom.

The commutation relations of creation and annihilation operators in a multiple-boson system are,

[a^{\,}_i, a^\dagger_j] \equiv a^{\,}_i a^\dagger_j - a^\dagger_ja^{\,}_i = \delta_{i j},
[a^\dagger_i, a^\dagger_j] = [a^{\,}_i, a^{\,}_j] = 0,

where [\ \ , \ \ ] is the commutator and \delta_{i j} is the Kronecker delta.

For fermions, the commutator is replaced by the anticommutator \{\ \ , \ \ \},

\{a^{\,}_i, a^\dagger_j\} \equiv a^{\,}_i a^\dagger_j +a^\dagger_j a^{\,}_i = \delta_{i j},
\{a^\dagger_i, a^\dagger_j\} = \{a^{\,}_i, a^{\,}_j\} = 0.

Therefore, exchanging disjoint (i.e. i \ne j) operators in a product of creation of annihilation operators will reverse the sign in fermion systems, but not in boson systems.

If the states labelled by i are an orthonormal basis of a Hilbert space H, then the result of this construction coincides with the CCR algebra and CAR algebra construction in the previous section but one. If they represent "eigenvectors" corresponding to the continuous spectrum of some operator, as for unbound particles in QFT, then the interpretation is more subtle.

See also[edit]



  1. ^ (Feynman 1998, p. 151)
  2. ^ (Feynman 1998, p. 167)
  3. ^ (Feynman 1998, pp. 174–5)
  4. ^ a b c Branson, Jim. "Quantum Physics at UCSD". Retrieved 16 May 2012.