Skip to content

Creation, Annihilation, and Second Quantization

Second quantization is the operator language of Fock space. Creation and annihilation operators add or remove one quantum from a mode, their algebra encodes bosonic or fermionic statistics, and one- and two-body operators become sums of mode-operator products.

The historical name can mislead. The procedure does not quantize an already quantized particle a second time. It rewrites many-particle quantum mechanics in a representation where exchange symmetry and variable particle number are built into the state space and operator algebra.

Use the volume orientation to choose the physical subsystem or mode decomposition, and Fock Space and Occupation Number to construct the state space. This chapter owns the operator construction that acts on that space.

The construction sequence is

choose one-particle modes↓define creation and annihilation maps↓impose the statistics algebra↓lift one- and two-body operators.\begin{gathered} \text{choose one-particle modes} \\ \downarrow \\ \text{define creation and annihilation maps} \\ \downarrow \\ \text{impose the statistics algebra} \\ \downarrow \\ \text{lift one- and two-body operators}. \end{gathered}
TaskCanonical pageMain output
add and remove mode quantaCreation and Annihilation Operatorsnormalized action on bosonic and fermionic number states
encode bosonic statisticsBosonic Commutation Relationscanonical commutators and symmetric Fock-space action
encode fermionic statisticsFermionic Anticommutation Relationscanonical anticommutators, exclusion, and ordering signs
count particles and excitationsNumber Operatorsmode number, total number, and conservation commutators
change between mode and position languageMode Expansionscompleteness, basis changes, momentum modes, and spin labels
use position-space creation and annihilation fieldsField Operatorsoperator-valued distributions, density, and field algebra
lift single-particle observablesOne-Body Operatorsbilinear mode and field expressions
lift pair interactionsTwo-Body Operatorsquartic expressions, matrix elements, and double-counting factor
assemble physical modelsMany-Particle Hamiltoniansgases, lattices, interactions, and pairing preview
reorder operator productsNormal Orderingcreation-left convention, fermionic signs, and reference vacuum
preview contraction bookkeepingWick’s Theorem PreviewGaussian-state contractions and scope limits
identify what changes in relativistic theoryBridge to QFTfree modes, fields, locality, antiparticles, and QFT boundaries

Let ∣n1,n2,…⟩\lvert n_1,n_2,\ldots\rangle be an occupation-number state. A bosonic creation operator acts as

ar†∣…,nr,…⟩=nr+1×∣…,nr+1,…⟩,\begin{aligned} a_r^\dagger \lvert\ldots,n_r,\ldots\rangle &= \sqrt{n_r+1} \\ &\quad\times \lvert\ldots,n_r+1,\ldots\rangle, \end{aligned}

while annihilation acts as

ar∣…,nr,…⟩=nr ∣…,nr−1,…⟩.a_r \lvert\ldots,n_r,\ldots\rangle = \sqrt{n_r}\, \lvert\ldots,n_r-1,\ldots\rangle.

The square-root factors preserve normalization and produce the canonical algebra. In particular, ar∣0⟩=0a_r\lvert0\rangle=0.

Fermionic operators change occupations between zero and one. With mode order 1,2,…1,2,\ldots fixed, define

ηr=(−1)∑s<rns.\eta_r = (-1)^{\sum_{s<r}n_s}.

Then

cr†∣…,nr,…⟩=ηr(1−nr)∣…,1,…⟩,cr∣…,nr,…⟩=ηrnr∣…,0,…⟩.\begin{aligned} c_r^\dagger \lvert\ldots,n_r,\ldots\rangle &= \eta_r(1-n_r) \lvert\ldots,1,\ldots\rangle, \\ c_r \lvert\ldots,n_r,\ldots\rangle &= \eta_r n_r \lvert\ldots,0,\ldots\rangle. \end{aligned}

The sign records how many occupied modes the operator passes in the chosen ordering. Different consistent conventions represent the same physics, but mixing conventions produces errors.

Bosonic mode operators satisfy

[ar,as†]=δrs,[ar,as]=0,[ar†,as†]=0.\begin{aligned} [a_r,a_s^\dagger]&=\delta_{rs}, \\ [a_r,a_s]&=0, \\ [a_r^\dagger,a_s^\dagger]&=0. \end{aligned}

Fermionic mode operators satisfy

{cr,cs†}=δrs,{cr,cs}=0,{cr†,cs†}=0.\begin{aligned} \{c_r,c_s^\dagger\}&=\delta_{rs}, \\ \{c_r,c_s\}&=0, \\ \{c_r^\dagger,c_s^\dagger\}&=0. \end{aligned}

Setting r=sr=s in the final relation gives

(cr†)2=0,(c_r^\dagger)^2=0,

which is Pauli exclusion in operator form. Bosonic creation operators commute, while fermionic creation operators anticommute and change sign under reordering.

For either statistics, use the appropriate mode operators to define

n^r=dr†dr,N^=∑rn^r,\hat n_r=d_r^\dagger d_r, \qquad \hat N=\sum_r\hat n_r,

where dr=ard_r=a_r for bosons and dr=crd_r=c_r for fermions.

The universal number commutators are

[N^,dr†]=dr†,[N^,dr]=−dr.\begin{aligned} [\hat N,d_r^\dagger]&=d_r^\dagger, \\ [\hat N,d_r]&=-d_r. \end{aligned}

Thus dr†d_r^\dagger raises total number by one and drd_r lowers it by one. Bilinears dr†dsd_r^\dagger d_s conserve total number because they remove one quantum and add one.

For fermions, n^r2=n^r\hat n_r^2=\hat n_r, so its eigenvalues are zero and one. Bosonic n^r\hat n_r has all nonnegative integer eigenvalues.

Let {ϕr(x)}\{\phi_r(\mathbf x)\} be a complete orthonormal basis of the one-particle Hilbert space. The nonrelativistic annihilation and creation fields are expanded as

ψ^(x)=∑rϕr(x)dr,ψ^†(x)=∑rϕr∗(x)dr†.\begin{aligned} \hat\psi(\mathbf x) &= \sum_r \phi_r(\mathbf x)d_r, \\ \hat\psi^\dagger(\mathbf x) &= \sum_r \phi_r^*(\mathbf x)d_r^\dagger. \end{aligned}

The inverse relation is

dr=∫d3x ϕr∗(x)ψ^(x).d_r = \int d^3x\, \phi_r^*(\mathbf x) \hat\psi(\mathbf x).

Completeness turns the discrete algebra into the equal-time field algebra. For bosons,

[ψ^(x),ψ^†(y)]=δ(3)(x−y),\left[ \hat\psi(\mathbf x), \hat\psi^\dagger(\mathbf y) \right] = \delta^{(3)}(\mathbf x-\mathbf y),

and for fermions the commutator is replaced by an anticommutator.

Field operators at sharp points are operator-valued distributions. Meaningful operators are obtained by smearing them with suitable wavefunctions or by placing the theory on a lattice or other regulated basis.

Basis Changes Preserve the Abstract Operator

Section titled “Basis Changes Preserve the Abstract Operator”

Suppose the one-particle bases are related by

∣α⟩=∑rUαr∣r⟩.\lvert\alpha\rangle = \sum_r U_{\alpha r}\lvert r\rangle.

The mode operators transform contragrediently so that the field and physical observables remain unchanged. With the convention above,

bα†=∑rUαrdr†,bα=∑rUαr∗dr.b_\alpha^\dagger = \sum_r U_{\alpha r}d_r^\dagger, \qquad b_\alpha = \sum_r U_{\alpha r}^*d_r.

Matrix elements and operator coefficients must be transformed together. A diagonal one-body Hamiltonian in one basis may contain off-diagonal hopping or mixing terms in another.

Let AA be a one-particle operator with matrix elements

Ars=⟨r∣A∣s⟩.A_{rs} = \langle r|A|s\rangle.

Its number-conserving many-body lift is

A^=∑r,sArsdr†ds.\widehat A = \sum_{r,s} A_{rs}d_r^\dagger d_s.

In position language,

A^=∫d3x ψ^†(x)Axψ^(x),\widehat A = \int d^3x\, \hat\psi^\dagger(\mathbf x) A_{\mathbf x} \hat\psi(\mathbf x),

where AxA_{\mathbf x} acts on the one-particle coordinate and any suppressed internal indices.

Examples include kinetic energy, external potentials, spin operators, mode energies, and hopping. The bilinear form is the canonical home of one-body physics in second quantization.

Let a pair interaction have matrix elements

Vrs;tu=⟨r,s∣V∣t,u⟩.V_{rs;tu} = \langle r,s|V|t,u\rangle.

For an interaction summed over unordered particle pairs, the second-quantized operator is

V^=12∑r,s,t,uVrs;tudr†ds†dudt.\widehat V = \frac12 \sum_{r,s,t,u} V_{rs;tu} d_r^\dagger d_s^\dagger d_u d_t.

The factor 1/21/2 removes double counting when the first-quantized interaction is ∑i<jV(i,j)\sum_{i<j}V(i,j). Index order matters, especially for fermions; changing the convention for Vrs;tuV_{rs;tu} requires changing the operator order consistently.

For a position-space pair potential,

V^=12∫d3x d3y ψ^†(x)ψ^†(y)×V(x,y)ψ^(y)ψ^(x).\begin{aligned} \widehat V =\frac12 \int d^3x\,d^3y\, &\hat\psi^\dagger(\mathbf x) \hat\psi^\dagger(\mathbf y) \\ &\times V(\mathbf x,\mathbf y) \hat\psi(\mathbf y) \hat\psi(\mathbf x). \end{aligned}

Contact interactions, Coulomb interactions, and Hubbard on-site terms are special cases after choosing modes and regularization.

A general number-conserving Hamiltonian with one- and two-body terms has the form

H=∑r,shrsdr†ds+12∑r,s,t,uVrs;tudr†ds†dudt.\begin{aligned} H &= \sum_{r,s}h_{rs}d_r^\dagger d_s \\ &\quad+ \frac12 \sum_{r,s,t,u} V_{rs;tu} d_r^\dagger d_s^\dagger d_u d_t. \end{aligned}

Every term contains equal numbers of creation and annihilation operators, so

[H,N^]=0.[H,\hat N]=0.

Pairing terms such as cr†cs†c_r^\dagger c_s^\dagger and cscrc_s c_r change particle number by two. They may appear in effective mean-field descriptions and preserve number parity even when they do not preserve N^\hat N. The physical interpretation and approximation must be stated.

Normal ordering moves creation operators to the left of annihilation operators. It is denoted by colons:

:drds†:={ds†dr,bosons,−ds†dr,fermions.:d_r d_s^\dagger: = \begin{cases} d_s^\dagger d_r,&\text{bosons}, \\ -d_s^\dagger d_r,&\text{fermions}. \end{cases}

Commutators or anticommutators generate the difference between an original product and its normal-ordered form. For one bosonic mode,

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

Normal ordering is defined relative to a chosen vacuum or reference state. In interacting many-body theory and QFT, changing the reference can change what “normal ordered” means.

Wick’s theorem reorganizes products or time-ordered products into normal-ordered terms plus contractions when the reference state has the required Gaussian structure. A contraction is bookkeeping for a two-point expectation value, not a new physical interaction.

The theorem is powerful because Gaussian states reduce higher correlators to two-point data. It does not apply unchanged to arbitrary interacting or non-Gaussian states. This volume provides only Wick’s Theorem Preview; full many-body and relativistic formulations belong in their canonical settings.

Quantum mechanics already supplies:

  • Fock space and number sectors;
  • creation and annihilation operators;
  • bosonic and fermionic algebras;
  • mode expansions;
  • operator-valued fields in position space;
  • normal ordering and contraction language.

Relativistic QFT adds essential structure: Lorentz symmetry, local relativistic fields, antiparticles, microcausality, interacting vacua, renormalization, and particle concepts that may depend on background or asymptotic regime.

A nonrelativistic field ψ^(x)\hat\psi(\mathbf x) annihilates a particle at equal time within the many-body model. It should not be identified without qualification with a relativistic scalar, spinor, or gauge field. Use Bridge to QFT for the handoff.

State the one-particle basis, internal labels, and fermionic ordering convention.

Use commutators for bosons and anticommutators for fermions. Do not mix the resulting reordering rules.

Compute hrsh_{rs} and Vrs;tuV_{rs;tu} from the one-particle basis before writing operator sums. Record index and complex-conjugation conventions.

One-body terms are bilinear. Two-body terms are quartic. Check the pair-counting factor and operator order.

5. Check Hermiticity and conservation laws

Section titled “5. Check Hermiticity and conservation laws”

Verify H†=HH^\dagger=H, particle-number or parity commutators, spin or lattice symmetries, and the action on simple occupation states.

Transform both mode operators and coefficient matrices. The abstract observable must not depend on notation.

Identify cutoffs, truncated modes, contact regularization, mean-field replacements, reference vacua, or nonrelativistic assumptions.

  • Omitting bosonic square-root factors. Number-state normalization then fails.
  • Ignoring fermionic mode ordering. Anticommutation signs depend on a consistent convention.
  • Using commutators for fermions or anticommutators for bosons. Statistics is encoded in the algebra.
  • Forgetting the pair factor 1/21/2. This double counts unordered interactions under the stated convention.
  • Reversing annihilation indices in a fermionic quartic term. Operator order and matrix-element convention must match.
  • Treating ψ^(x)\hat\psi(\mathbf x) as an ordinary function. Sharp-point fields are operator-valued distributions.
  • Calling every quartic term an interaction. Normal ordering and basis transformations can also produce quartic expressions in effective descriptions; inspect the model.
  • Assuming normal ordering is reference independent. The vacuum or reference state matters.
  • Calling second quantization a second physical quantization step. It is a many-body representation and operator formalism.
  • Equating the nonrelativistic field language with full QFT. Relativistic locality and particle interpretation add essential content.

Operator foundation: Creation and Annihilation Operators → Bosonic Commutation Relations or Fermionic Anticommutation Relations → Number Operators.

Mode and field language: Mode Expansions → Field Operators → Bridge to QFT.

Hamiltonian construction: One-Body Operators → Two-Body Operators → Many-Particle Hamiltonians.

Reordering and contractions: Normal Ordering → Wick’s Theorem Preview → Harmonic-Oscillator Ladder Operators.

Continuous-variable bridge: Mode Expansions → Continuous Variables and Modes → Gaussian States Preview.

Across-fields bridge: Many-Particle Hamiltonians → Entanglement Across Fields → Entanglement in Many-Body Physics.

Practice bridge: Reference, Problems, and Notebooks → Fock Space Exercises → Computational Notebooks.

  • A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, Dover, 2003.
  • J. W. Negele and H. Orland, Quantum Many-Particle Systems, Westview Press, 1998.
  • A. Altland and B. Simons, Condensed Matter Field Theory, 3rd ed., Cambridge University Press, 2023.
  • H. Bruus and K. Flensberg, Many-Body Quantum Theory in Condensed Matter Physics, Oxford University Press, 2004.
  • A. F. G. Wyatt, Excitations in Quantum Fluids, Cambridge University Press, 1999.
  • M. O. Scully and M. S. Zubairy, Quantum Optics, Cambridge University Press, 1997.
  • M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory, Westview Press, 1995.

Exercise 1: Recover the bosonic commutator

Section titled “Exercise 1: Recover the bosonic commutator”

Act with [a,a†][a,a^\dagger] on a one-mode number state ∣n⟩\lvert n\rangle and show that the result is ∣n⟩\lvert n\rangle.

Solution

First,

aa†∣n⟩=(n+1)∣n⟩,a†a∣n⟩=n∣n⟩.\begin{aligned} a a^\dagger\lvert n\rangle &= (n+1)\lvert n\rangle, \\ a^\dagger a\lvert n\rangle &= n\lvert n\rangle. \end{aligned}

Subtracting gives

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

Because the number states span the Fock space on the natural dense domain, this realizes [a,a†]=I[a,a^\dagger]=I there.

Exercise 2: Fermionic exclusion from the algebra

Section titled “Exercise 2: Fermionic exclusion from the algebra”

Use {cr†,cr†}=0\{c_r^\dagger,c_r^\dagger\}=0 to show that the same fermionic mode cannot be created twice.

Solution

The anticommutator gives

2(cr†)2=0,2(c_r^\dagger)^2=0,

so (cr†)2=0(c_r^\dagger)^2=0. Acting twice on any state therefore gives the zero vector. In particular,

(cr†)2∣0⟩=0.(c_r^\dagger)^2\lvert0\rangle=0.

This is the operator-algebra form of Pauli exclusion for a complete one-particle mode.

Exercise 3: A diagonal one-body Hamiltonian

Section titled “Exercise 3: A diagonal one-body Hamiltonian”

Let the one-particle Hamiltonian satisfy h∣r⟩=ϵr∣r⟩h\lvert r\rangle=\epsilon_r\lvert r\rangle. Write its second-quantized form and find its action on ∣n1,n2,…⟩\lvert n_1,n_2,\ldots\rangle.

Solution

The matrix elements are hrs=ϵrδrsh_{rs}=\epsilon_r\delta_{rs}, so

H0=∑rϵrdr†dr=∑rϵrn^r.H_0 = \sum_r\epsilon_r d_r^\dagger d_r = \sum_r\epsilon_r\hat n_r.

Therefore

H0∣n⟩=En∣n⟩,En=∑rϵrnr,\begin{aligned} H_0\lvert\mathbf n\rangle &= E_{\mathbf n}\lvert\mathbf n\rangle, \\ E_{\mathbf n} &= \sum_r\epsilon_r n_r, \end{aligned}

Here ∣n⟩=∣n1,n2,…⟩\lvert\mathbf n\rangle=\lvert n_1,n_2,\ldots\rangle. The result applies to bosons or fermions with the corresponding allowed occupations.

Suppose a first-quantized interaction is ∑i<jV(i,j)\sum_{i<j}V(i,j). Explain why a symmetric sum over ordered mode indices in second quantization carries a factor 1/21/2.

Solution

The unordered pair sum counts each physical pair once. A fully expanded operator sum naturally includes both orderings associated with the two particle slots or field coordinates. For a symmetric pair potential, the contributions related by exchanging the two particles are equal.

Multiplying the ordered sum by 1/21/2 removes this double counting. If the matrix elements are defined with antisymmetrization or another restricted convention already built in, the prefactor must be checked rather than copied mechanically.

Exercise 5: Number conservation of a quartic term

Section titled “Exercise 5: Number conservation of a quartic term”

Show that

Qrstu=dr†ds†dudtQ_{rstu} = d_r^\dagger d_s^\dagger d_u d_t

commutes with total number N^\hat N for bosons or fermions.

Solution

Use [N^,dj†]=dj†[\hat N,d_j^\dagger]=d_j^\dagger, [N^,dj]=−dj[\hat N,d_j]=-d_j, and the product rule for commutators. Each creation operator contributes +1+1 and each annihilation operator contributes −1-1:

[N^,Qrstu]=(1+1−1−1)Qrstu=0.\begin{aligned} [\hat N,Q_{rstu}] &= (1+1-1-1)Q_{rstu} \\ &=0. \end{aligned}

The argument uses ordinary commutators with the even operator N^\hat N and works for either underlying statistics algebra.