Skip to content

Bosonic Commutation Relations

Bosonic commutation relations are the operator algebra that makes bosonic occupation-number notation work. They encode two facts at once: a bosonic mode may be occupied by any nonnegative integer number of particles or excitations, and creation operators for different modes commute because exchanging identical bosons does not introduce a minus sign.

For modes labeled by i,ji,j, the canonical bosonic algebra is

[ai,aj†]=δijI,[ai,aj]=0,[ai†,aj†]=0.[a_i,a_j^\dagger] = \delta_{ij}I, \qquad [a_i,a_j] =0, \qquad [a_i^\dagger,a_j^\dagger] =0.

Here

[A,B]=AB−BA[A,B] = AB-BA

is the commutator, and II is the identity operator on the relevant Fock space. Many texts suppress II and write [ai,aj†]=δij[a_i,a_j^\dagger]=\delta_{ij}, but the right-hand side is still an operator identity.

For one bosonic mode, the algebra reduces to

[a,a†]=I.[a,a^\dagger]=I.

Equivalently,

aa†=a†a+I.aa^\dagger = a^\dagger a+I.

This identity is the algebraic source of the familiar one-step offset between creating first and annihilating first. If one creates a quantum and then removes it, there is one extra way to return to the original state compared with removing first and then creating.

For many modes,

[ai,aj†]=0when i≠j.[a_i,a_j^\dagger] =0 \qquad \text{when }i\ne j.

Thus operations on distinct bosonic modes commute. Creating a boson in mode ii and then one in mode jj gives the same state as doing those operations in the opposite order:

ai†aj†∣0⟩=aj†ai†∣0⟩.a_i^\dagger a_j^\dagger\lvert0\rangle = a_j^\dagger a_i^\dagger\lvert0\rangle.

This is the operator form of bosonic exchange symmetry.

Let ∣n⟩\lvert n\rangle be the normalized occupation-number state of one bosonic mode. The creation and annihilation operators act by

a†∣n⟩=n+1 ∣n+1⟩,a∣n⟩=n ∣n−1⟩,a^\dagger\lvert n\rangle = \sqrt{n+1}\,\lvert n+1\rangle, \qquad a\lvert n\rangle = \sqrt n\,\lvert n-1\rangle,

with

a∣0⟩=0.a\lvert0\rangle =0.

These formulas imply the commutation relation on every number state:

aa†∣n⟩=(n+1)∣n⟩,a†a∣n⟩=n∣n⟩.aa^\dagger\lvert n\rangle = (n+1)\lvert n\rangle, \qquad a^\dagger a\lvert n\rangle = n\lvert n\rangle.

Subtracting gives

[a,a†]∣n⟩=∣n⟩.[a,a^\dagger]\lvert n\rangle = \lvert n\rangle.

Since the number states form the standard basis of the single-mode bosonic Fock space, this is the identity operator on that space.

The square-root factors are forced by this algebra together with normalization. If one tried to use a†∣n⟩=∣n+1⟩a^\dagger\lvert n\rangle=\lvert n+1\rangle for every nn, the commutator would not equal the identity on the normalized number basis.

The number operator for mode ii is

Ni=ai†ai.N_i = a_i^\dagger a_i.

On the occupation-number basis,

Ni∣n1,n2,…⟩B=ni∣n1,n2,…⟩B.N_i \lvert n_1,n_2,\ldots\rangle_B = n_i \lvert n_1,n_2,\ldots\rangle_B.

The total number operator is

N=∑iNi.N = \sum_i N_i.

Using the commutation relations, one finds

[Ni,aj†]=δijaj†,[Ni,aj]=−δijaj.[N_i,a_j^\dagger] = \delta_{ij}a_j^\dagger, \qquad [N_i,a_j] = -\delta_{ij}a_j.

These identities say exactly what the words say: aj†a_j^\dagger raises the occupation of mode jj by one, while aja_j lowers it by one. For i≠ji\ne j, the occupation of mode ii is unchanged.

A bosonic occupation state is written

∣n1,n2,…,ni,…⟩B.\lvert n_1,n_2,\ldots,n_i,\ldots\rangle_B.

The mode-ii operators act only on the iith occupation number:

ai†∣…,ni,…⟩B=ni+1 ∣…,ni+1,…⟩B,a_i^\dagger \lvert\ldots,n_i,\ldots\rangle_B = \sqrt{n_i+1}\, \lvert\ldots,n_i+1,\ldots\rangle_B,

and

ai∣…,ni,…⟩B=ni ∣…,ni−1,…⟩B.a_i \lvert\ldots,n_i,\ldots\rangle_B = \sqrt{n_i}\, \lvert\ldots,n_i-1,\ldots\rangle_B.

For two distinct modes 11 and 22,

a1a2†∣n1,n2⟩B=n1n2+1 ∣n1−1,n2+1⟩B.a_1a_2^\dagger \lvert n_1,n_2\rangle_B = \sqrt{n_1}\sqrt{n_2+1}\, \lvert n_1-1,n_2+1\rangle_B.

Because [a1,a2†]=0[a_1,a_2^\dagger]=0, the same result is obtained if one applies a2†a1a_2^\dagger a_1 instead.

Repeated creation in the same mode is allowed. Starting from the vacuum,

(ai†)n∣0⟩≠0for every n=1,2,….(a_i^\dagger)^n\lvert0\rangle \ne 0 \qquad \text{for every }n=1,2,\ldots .

This is the sharp operator contrast with fermions, where the corresponding creation operator squares to zero.

The bosonic Fock space over a one-particle Hilbert space h\mathcal h is

FB(h)=⨁N=0∞Sym⁡Nh.\mathcal F_B(\mathcal h) = \bigoplus_{N=0}^{\infty} \operatorname{Sym}^N\mathcal h.

The commutation relations are not an extra rule pasted on top of this space. They are the efficient operator representation of the symmetric tensor powers.

The same algebra also determines the commutator of the position-space field operators once the field is expanded in a one-particle basis; see Mode Expansions.

Given an orthonormal mode basis, a normalized occupation state is

∣n1,n2,…⟩B=∏i(ai†)nini!∣0⟩,\lvert n_1,n_2,\ldots\rangle_B = \prod_i \frac{(a_i^\dagger)^{n_i}}{\sqrt{n_i!}} \lvert0\rangle,

where only finitely many nin_i are nonzero in the finite-particle sector. The product ordering is irrelevant for bosonic creation operators because

[ai†,aj†]=0.[a_i^\dagger,a_j^\dagger]=0.

For two distinct one-particle modes ∣u⟩\lvert u\rangle and ∣v⟩\lvert v\rangle,

au†av†∣0⟩⟷12(∣u⟩1∣v⟩2+∣v⟩1∣u⟩2).a_u^\dagger a_v^\dagger\lvert0\rangle \longleftrightarrow \frac{1}{\sqrt2} \bigl( \lvert u\rangle_1\lvert v\rangle_2 + \lvert v\rangle_1\lvert u\rangle_2 \bigr).

For two bosons in the same mode,

(au†)22!∣0⟩⟷∣u⟩1∣u⟩2.\frac{(a_u^\dagger)^2}{\sqrt{2!}} \lvert0\rangle \longleftrightarrow \lvert u\rangle_1\lvert u\rangle_2.

The factorial in the occupation-state formula exactly compensates for the repeated symmetric copies of the same mode.

For one oscillator-like bosonic mode,

N=a†a,N∣n⟩=n∣n⟩.N=a^\dagger a, \qquad N\lvert n\rangle = n\lvert n\rangle.

The Hamiltonian of a harmonic oscillator can be written

H=ℏω(N+12),H = \hbar\omega \left( N+\frac12 \right),

so the number basis is also the energy basis for that simple system.

For two modes, the normalized state with three bosons in mode 11 and one boson in mode 22 is

∣3,1⟩B=(a1†)3a2†3!∣0⟩.\lvert3,1\rangle_B = \frac{(a_1^\dagger)^3a_2^\dagger}{\sqrt{3!}} \lvert0\rangle.

Applying a1a_1 gives

a1∣3,1⟩B=3 ∣2,1⟩B.a_1\lvert3,1\rangle_B = \sqrt3\,\lvert2,1\rangle_B.

Applying a2†a_2^\dagger gives

a2†∣3,1⟩B=2 ∣3,2⟩B.a_2^\dagger\lvert3,1\rangle_B = \sqrt2\,\lvert3,2\rangle_B.

These square roots are often where arithmetic mistakes enter many-body calculations.

  • Writing [ai,aj†]=1[a_i,a_j^\dagger]=1 for all modes instead of δijI\delta_{ij}I.
  • Omitting the identity operator conceptually on the right-hand side of [a,a†]=I[a,a^\dagger]=I.
  • Replacing commutators by anticommutators for bosonic modes.
  • Dropping the factors n\sqrt n and n+1\sqrt{n+1} in actions on normalized number states.
  • Assuming that bosonic occupation numbers stop at 11.
  • Confusing a harmonic oscillator ladder operator for one particle in a potential with a Fock-space operator that creates or removes a mode excitation.
  • P. A. M. Dirac, The Principles of Quantum Mechanics, 4th ed., Oxford University Press, 1958.
  • A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, McGraw-Hill, 1971.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • L. E. Ballentine, Quantum Mechanics: A Modern Development, 2nd ed., World Scientific, 2014.
  • A. Altland and B. Simons, Condensed Matter Field Theory, 2nd ed., Cambridge University Press, 2010.
  1. Verify [a,a†]∣n⟩=∣n⟩[a,a^\dagger]\lvert n\rangle=\lvert n\rangle from the action of aa and a†a^\dagger on number states.
Solution

First,

aa†∣n⟩=an+1 ∣n+1⟩=(n+1)∣n⟩.aa^\dagger\lvert n\rangle = a\sqrt{n+1}\,\lvert n+1\rangle = (n+1)\lvert n\rangle.

Second,

a†a∣n⟩=a†n ∣n−1⟩=n∣n⟩.a^\dagger a\lvert n\rangle = a^\dagger\sqrt n\,\lvert n-1\rangle = n\lvert n\rangle.

Therefore

[a,a†]∣n⟩=((n+1)−n)∣n⟩=∣n⟩.[a,a^\dagger]\lvert n\rangle = \bigl((n+1)-n\bigr)\lvert n\rangle = \lvert n\rangle.
  1. Show that a1†a2†∣0⟩=a2†a1†∣0⟩a_1^\dagger a_2^\dagger\lvert0\rangle=a_2^\dagger a_1^\dagger\lvert0\rangle for two bosonic modes.
Solution

For bosons,

[a1†,a2†]=0.[a_1^\dagger,a_2^\dagger]=0.

Thus

a1†a2†=a2†a1†,a_1^\dagger a_2^\dagger = a_2^\dagger a_1^\dagger,

as operators, so they give the same result on the vacuum:

a1†a2†∣0⟩=a2†a1†∣0⟩.a_1^\dagger a_2^\dagger\lvert0\rangle = a_2^\dagger a_1^\dagger\lvert0\rangle.
  1. Compute [Ni,ai†]∣ni⟩[N_i,a_i^\dagger]\lvert n_i\rangle for one mode ii.
Solution

Using Ni∣ni⟩=ni∣ni⟩N_i\lvert n_i\rangle=n_i\lvert n_i\rangle,

Niai†∣ni⟩=Nini+1 ∣ni+1⟩=(ni+1)ni+1 ∣ni+1⟩.N_i a_i^\dagger\lvert n_i\rangle = N_i\sqrt{n_i+1}\,\lvert n_i+1\rangle = (n_i+1)\sqrt{n_i+1}\,\lvert n_i+1\rangle.

Also,

ai†Ni∣ni⟩=nini+1 ∣ni+1⟩.a_i^\dagger N_i\lvert n_i\rangle = n_i\sqrt{n_i+1}\,\lvert n_i+1\rangle.

Subtracting gives

[Ni,ai†]∣ni⟩=ni+1 ∣ni+1⟩=ai†∣ni⟩.[N_i,a_i^\dagger]\lvert n_i\rangle = \sqrt{n_i+1}\,\lvert n_i+1\rangle = a_i^\dagger\lvert n_i\rangle.
  1. Normalize (a†)3∣0⟩(a^\dagger)^3\lvert0\rangle.
Solution

By repeated use of the creation formula,

(a†)3∣0⟩=123 ∣3⟩=3! ∣3⟩.(a^\dagger)^3\lvert0\rangle = \sqrt{1}\sqrt{2}\sqrt{3}\,\lvert3\rangle = \sqrt{3!}\,\lvert3\rangle.

Therefore the normalized state is

∣3⟩=(a†)33!∣0⟩.\lvert3\rangle = \frac{(a^\dagger)^3}{\sqrt{3!}}\lvert0\rangle.