Skip to content

Second Quantization

Second quantization is the occupation-number and operator language for identical-particle quantum mechanics and quantum fields. It replaces explicit symmetrization or antisymmetrization of coordinate wavefunctions with bosonic commutators or fermionic anticommutators on Fock space.

The name is historical. In nonrelativistic many-body theory, one is usually rewriting an already quantum NN-particle theory in a representation that allows variable occupation. In canonical field theory, one instead begins with a classical field phase space and quantizes its modes. Those routes meet in similar creation, annihilation, and Fock-space structures, but they are not the same logical operation.

Before using this bridge, be comfortable with:

  • one-particle Hilbert spaces and orthonormal bases;
  • tensor products and permutation symmetry;
  • identical bosons and fermions;
  • the harmonic oscillator’s ladder algebra;
  • direct sums and occupation-number notation;
  • commutators, anticommutators, and operator ordering.

The canonical many-particle construction is developed in Second Quantization: Bridge to QFT. This page emphasizes the dictionary and its field-theory boundary.

Let h\mathcal h be the one-particle Hilbert space. The unsymmetrized NN-particle space is h⊗N\mathcal h^{\otimes N}. Identical particles occupy its symmetric or antisymmetric subspace:

h+(N)=Sym⁡Nh,h−(N)=⋀Nh.\mathcal h_+^{(N)} =\operatorname{Sym}^N\mathcal h, \qquad \mathcal h_-^{(N)} =\bigwedge^N\mathcal h.

Bosonic and fermionic Fock spaces collect all particle-number sectors:

F±(h)=⨁N=0∞h±(N).\mathcal F_\pm(\mathcal h) = \bigoplus_{N=0}^{\infty} \mathcal h_\pm^{(N)}.

The N=0N=0 sector is one dimensional:

h±(0)≅C,\mathcal h_\pm^{(0)} \cong\mathbb C,

and its normalized basis vector is the vacuum ∣0⟩\lvert0\rangle.

Fock space does not require the Hamiltonian to change particle number. It provides a common arena in which number-conserving and number-changing operators can both be represented.

Choose an ordered orthonormal basis {∣i⟩}\lbrace\lvert i\rangle\rbrace of h\mathcal h.

Bosonic operators satisfy

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

Their occupation basis is

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

The action on mode jj is

aj∣…,nj,…⟩=nj ∣…,nj−1,…⟩,aj†∣…,nj,…⟩=nj+1 ∣…,nj+1,…⟩.\begin{aligned} a_j\lvert\ldots,n_j,\ldots\rangle &= \sqrt{n_j}\, \lvert\ldots,n_j-1,\ldots\rangle, \\ a_j^\dagger\lvert\ldots,n_j,\ldots\rangle &= \sqrt{n_j+1}\, \lvert\ldots,n_j+1,\ldots\rangle. \end{aligned}

Fermionic operators satisfy

{ci,cj†}=δijI,{ci,cj}={ci†,cj†}=0.\lbrace c_i,c_j^\dagger\rbrace =\delta_{ij}I, \qquad \lbrace c_i,c_j\rbrace =\lbrace c_i^\dagger,c_j^\dagger\rbrace=0.

Because (ci†)2=0(c_i^\dagger)^2=0, each occupation is ni=0n_i=0 or 11. Fix the basis ordering through

∣n1,n2,…⟩=(c1†)n1(c2†)n2⋯∣0⟩.\lvert n_1,n_2,\ldots\rangle = (c_1^\dagger)^{n_1} (c_2^\dagger)^{n_2}\cdots \lvert0\rangle.

Then

cj∣n⟩=(−1)∑i<jninj∣n1,…,0j,…⟩,cj†∣n⟩=(−1)∑i<jni(1−nj)∣n1,…,1j,…⟩.\begin{aligned} c_j\lvert\boldsymbol n\rangle &= (-1)^{\sum_{i<j}n_i} n_j \lvert n_1,\ldots,0_j,\ldots\rangle, \\ c_j^\dagger\lvert\boldsymbol n\rangle &= (-1)^{\sum_{i<j}n_i} (1-n_j) \lvert n_1,\ldots,1_j,\ldots\rangle. \end{aligned}

The sign is not an optional convention once the mode order has been fixed. Changing the mode order changes basis-state phases and must be propagated consistently through every operator.

The total number operator is

N=∑iai†aiN =\sum_i a_i^\dagger a_i

for bosons, or the analogous expression N=∑ici†ciN=\sum_i c_i^\dagger c_i for fermions. It satisfies

[N,ai†]=ai†,[N,ai]=−ai,[N,a_i^\dagger]=a_i^\dagger, \qquad [N,a_i]=-a_i,

and likewise for ci†,cic_i^\dagger,c_i. Thus creation and annihilation operators map

h±(N)⟶h±(N±1).\mathcal h_\pm^{(N)} \longrightarrow \mathcal h_\pm^{(N\pm1)}.

A Hamiltonian conserves total particle number exactly when

[H,N]=0.[H,N]=0.

Pairing terms such as ci†cj†+h.c.c_i^\dagger c_j^\dagger+\text{h.c.} do not conserve NN, though they may conserve fermion-number parity.

Suppose a new one-particle basis is

∣α⟩=∑i∣i⟩Uiα,U†U=I.\lvert\alpha\rangle =\sum_i\lvert i\rangle U_{i\alpha}, \qquad U^\dagger U=I.

The corresponding operators transform as

bα†=∑iUiαai†,bα=∑iUiα∗ai.b_\alpha^\dagger =\sum_i U_{i\alpha}a_i^\dagger, \qquad b_\alpha =\sum_i U_{i\alpha}^\ast a_i.

Unitarity preserves the canonical algebra:

[bα,bβ†]=δαβI[b_\alpha,b_\beta^\dagger] =\delta_{\alpha\beta}I

for bosons, with the analogous fermionic anticommutator. Creation and annihilation operators therefore depend on a chosen mode basis; the underlying many-particle operator does not.

Let hh be a one-particle operator with matrix elements

hij=⟨i∣h∣j⟩.h_{ij}=\langle i\vert h\vert j\rangle.

On the NN-particle sector, the corresponding additive operator is

h(1)+⋯+h(N).h^{(1)}+\cdots+h^{(N)}.

Its Fock-space lift is

dΓ(h)=∑ijhijai†ajd\Gamma(h) =\sum_{ij} h_{ij}a_i^\dagger a_j

for bosons, and

dΓ(h)=∑ijhijci†cjd\Gamma(h) =\sum_{ij} h_{ij}c_i^\dagger c_j

for fermions. The notation dΓ(h)d\Gamma(h) emphasizes that this is the differential second quantization of a one-particle operator.

Important examples are

One-particle operatorFock-space operator
Identity IhI_{\mathcal h}Number operator N=dΓ(Ih)N=d\Gamma(I_{\mathcal h})
One-particle Hamiltonian hhFree many-body Hamiltonian ∑ijhijai†aj\sum_{ij}h_{ij}a_i^\dagger a_j
Momentum p\boldsymbol pTotal momentum ∑ij⟨i∣p∣j⟩ai†aj\sum_{ij}\langle i\vert\boldsymbol p\vert j\rangle a_i^\dagger a_j
One-particle symmetry generator ggAdditive many-body generator dΓ(g)d\Gamma(g)

Every number-conserving one-body lift commutes with NN.

For a symmetric pair interaction vv acting on two one-particle slots, define ordered-basis matrix elements

Vij;kl=⟨i,j∣v∣k,l⟩.V_{ij;kl} = \langle i,j\vert v\vert k,l\rangle.

With unsymmetrized matrix elements and the displayed operator order, the second-quantized interaction is

V^=12∑ijklVij;klai†aj†alak,\widehat V =\frac12 \sum_{ijkl} V_{ij;kl} a_i^\dagger a_j^\dagger a_l a_k,

or the same expression with cc operators for fermions.

The factor 1/21/2 removes double counting of particle pairs. Some sources use antisymmetrized fermionic matrix elements, change index order, or absorb combinatorial factors into Vij;klV_{ij;kl}. Copy the matrix-element definition together with the operator formula.

The operator contains two creators and two annihilators, so it conserves total particle number.

Field Operators in Nonrelativistic Many-Body Theory

Section titled “Field Operators in Nonrelativistic Many-Body Theory”

Let φi(x)=⟨x∣i⟩\varphi_i(\boldsymbol x)=\langle\boldsymbol x\vert i\rangle. Define

ψ(x)=∑iφi(x)ai,ψ†(x)=∑iφi∗(x)ai†.\psi(\boldsymbol x) =\sum_i \varphi_i(\boldsymbol x)a_i, \qquad \psi^\dagger(\boldsymbol x) =\sum_i \varphi_i^\ast(\boldsymbol x)a_i^\dagger.

For fermions, replace aia_i by cic_i. Completeness of the one-particle basis gives

[ψ(x),ψ†(y)]=δ(d)(x−y)I[\psi(\boldsymbol x),\psi^\dagger(\boldsymbol y)] = \delta^{(d)} (\boldsymbol x-\boldsymbol y)I

for bosons, and

{ψ(x),ψ†(y)}=δ(d)(x−y)I\lbrace \psi(\boldsymbol x),\psi^\dagger(\boldsymbol y) \rbrace = \delta^{(d)} (\boldsymbol x-\boldsymbol y)I

for fermions.

The field operator has dimensions

[ψ]=L−d/2[\psi]=L^{-d/2}

under this delta-function convention. It annihilates a particle at position x\boldsymbol x in a distributional sense; it is not a complex-valued single-particle wavefunction.

For a one-particle kernel h(x,y)h(\boldsymbol x,\boldsymbol y),

H^1=∫ddx ddy ψ†(x)h(x,y)ψ(y).\widehat H_1 = \int d^dx\,d^dy\, \psi^\dagger(\boldsymbol x) h(\boldsymbol x,\boldsymbol y) \psi(\boldsymbol y).

For a local differential one-particle Hamiltonian,

H^1=∫ddx ψ†(x)[−ℏ2∇22m+U(x)]ψ(x).\widehat H_1 = \int d^dx\, \psi^\dagger(\boldsymbol x) \left[ -\frac{\hbar^2\nabla^2}{2m} +U(\boldsymbol x) \right] \psi(\boldsymbol x).

A two-body potential has the standard form

V^=12∫ddx ddy ψ†(x)ψ†(y)V(x,y)×ψ(y)ψ(x).\begin{aligned} \widehat V =\frac12 \int d^dx\,d^dy\, &\psi^\dagger(\boldsymbol x) \psi^\dagger(\boldsymbol y) V(\boldsymbol x,\boldsymbol y) \\ &\times \psi(\boldsymbol y) \psi(\boldsymbol x). \end{aligned}

For fermions, the written order fixes the signs. Equal-point products can also require regularization in continuum models.

Let ΨN(x1,…,xN)\Psi_N(\boldsymbol x_1,\ldots,\boldsymbol x_N) be normalized and symmetric for bosons or antisymmetric for fermions. Its Fock-space vector is

∣ΨN⟩=1N!∫ddx1⋯ddxN ΨN(x1,…,xN)×ψ†(x1)⋯ψ†(xN)∣0⟩.\begin{aligned} \lvert\Psi_N\rangle =\frac{1}{\sqrt{N!}} \int d^dx_1\cdots d^dx_N\, &\Psi_N(\boldsymbol x_1,\ldots,\boldsymbol x_N) \\ &\times \psi^\dagger(\boldsymbol x_1)\cdots \psi^\dagger(\boldsymbol x_N) \lvert0\rangle. \end{aligned}

Conversely,

ΨN(x1,…,xN)=1N!⟨0∣ψ(xN)⋯ψ(x1)∣ΨN⟩.\begin{aligned} \Psi_N(\boldsymbol x_1,\ldots,\boldsymbol x_N) =\frac{1}{\sqrt{N!}} \langle0\vert &\psi(\boldsymbol x_N)\cdots \\ &\psi(\boldsymbol x_1) \vert\Psi_N\rangle. \end{aligned}

This equivalence shows what second quantization does: it packages exchange symmetry into the operator algebra. It does not introduce a different empirical theory for a fixed nonrelativistic particle number.

  • Exchange symmetry is automatic once the operator algebra is chosen.
  • Variable particle number is represented without changing Hilbert spaces.
  • Additive observables and few-body interactions have compact formulas.
  • Lattice, momentum, orbital, band, and position bases are related by unitary mode transformations.
  • Collective excitations and quasiparticles can be introduced by new operator bases.
  • Perturbations that create or destroy excitations can be written directly.
  • The same algebra prepares the mode language of quantum field theory.

Nonrelativistic field operators and relativistic quantum fields look similar, but several new ingredients enter.

For a relativistic field, canonical quantization begins from a classical field and its conjugate momentum. Normal modes become operator modes after imposing commutation or anticommutation relations. The Harmonic Oscillator to Fields card develops the bosonic free-field case.

A complex relativistic field generally contains distinct particle and antiparticle creation operators. A Dirac field combines electron and positron modes. This structure is not captured by merely allowing the particle number of a nonrelativistic Schrödinger field to vary.

Relativistic fields must satisfy locality conditions at spacelike separation. The connection between integer spin and bosonic commutation, and between half-integer spin and fermionic anticommutation, is a QFT theorem with assumptions including Lorentz invariance, locality, and positive energy. In nonrelativistic many-body mechanics, Bose or Fermi statistics is imposed as the particle type.

Quantum fields are operator-valued distributions. Products at the same point can be singular, so expressions such as ϕ4(x)\phi^4(x) or ψ‾(x)ψ(x)\overline\psi(x)\psi(x) require a regulator and renormalization prescription.

Free Fock space supplies perturbative and asymptotic particle language, but an interacting continuum QFT need not be globally equivalent to the free representation. “Create a particle by acting once with the field” is therefore an approximation or spectral statement, not a universal exact identity.

Before using a second-quantized formula, record:

  1. bosonic commutators or fermionic anticommutators;
  2. ordered mode basis, especially for fermions;
  3. normalization of continuum modes and delta functions;
  4. whether matrix elements are symmetrized or antisymmetrized;
  5. index order in two-body coefficients;
  6. whether total number, a species charge, or only parity is conserved;
  7. normal-ordering convention;
  8. whether the field is nonrelativistic, relativistic real, relativistic complex, spinor, or gauge;
  9. finite-volume versus continuum normalization.
  • Describing second quantization as literally quantizing the same theory twice.
  • Forgetting the vacuum sector in Fock space.
  • Treating creation operators as basis-independent objects without mode labels.
  • Reordering fermionic operators without the required sign.
  • Using bosonic occupation factors for fermions.
  • Combining antisymmetrized matrix elements with a prefactor meant for unsymmetrized elements.
  • Calling ψ(x)\psi(\boldsymbol x) a many-particle wavefunction.
  • Assuming Fock-space notation implies that the Hamiltonian changes particle number.
  • Inferring relativistic spin–statistics from nonrelativistic notation alone.
  • Treating products of continuum field operators at one point as automatically finite.

Show that a one-body lift

A^=∑ijAijai†aj\widehat A=\sum_{ij}A_{ij}a_i^\dagger a_j

commutes with N=∑kak†akN=\sum_k a_k^\dagger a_k.

Solution

Use

[N,ai†]=ai†,[N,aj]=−aj.[N,a_i^\dagger]=a_i^\dagger, \qquad [N,a_j]=-a_j.

The product rule gives

[N,ai†aj]=[N,ai†]aj+ai†[N,aj]=ai†aj−ai†aj=0.\begin{aligned} [N,a_i^\dagger a_j] &= [N,a_i^\dagger]a_j +a_i^\dagger[N,a_j] \\ &= a_i^\dagger a_j-a_i^\dagger a_j =0. \end{aligned}

Therefore every term commutes with NN, so [N,A^]=0[N,\widehat A]=0. The same argument works for fermions because NN and the bilinear are even operators and the ordinary commutator is used.

Using the convention

∣n1,n2,n3⟩=(c1†)n1(c2†)n2(c3†)n3∣0⟩,\lvert n_1,n_2,n_3\rangle =(c_1^\dagger)^{n_1} (c_2^\dagger)^{n_2} (c_3^\dagger)^{n_3} \lvert0\rangle,

calculate c2†∣1,0,1⟩c_2^\dagger\lvert1,0,1\rangle.

Solution

One occupied mode lies before mode 22, so

(−1)∑i<2ni=(−1)n1=−1.(-1)^{\sum_{i<2}n_i}=(-1)^{n_1}=-1.

Because n2=0n_2=0, creation is allowed:

c2†∣1,0,1⟩=−∣1,1,1⟩.c_2^\dagger\lvert1,0,1\rangle =-\lvert1,1,1\rangle.

The minus sign records the reordering needed to place c2†c_2^\dagger after c1†c_1^\dagger in the chosen canonical basis order.

Verify that a unitary mode transformation preserves the bosonic canonical commutator.

Solution

With

bα=∑iUiα∗ai,bβ†=∑jUjβaj†,b_\alpha =\sum_iU_{i\alpha}^\ast a_i, \qquad b_\beta^\dagger =\sum_jU_{j\beta}a_j^\dagger,

one obtains

[bα,bβ†]=∑ijUiα∗Ujβ[ai,aj†]=∑iUiα∗Uiβ=(U†U)αβ=δαβ.\begin{aligned} [b_\alpha,b_\beta^\dagger] &= \sum_{ij} U_{i\alpha}^\ast U_{j\beta} [a_i,a_j^\dagger] \\ &= \sum_i U_{i\alpha}^\ast U_{i\beta} \\ &= (U^\dagger U)_{\alpha\beta} =\delta_{\alpha\beta}. \end{aligned}

The fermionic calculation is identical with anticommutators.

Continue with canonical quantization of scalar and spinor fields, relativistic normalization, particles and antiparticles, Wick’s theorem, interacting fields, and the spin–statistics theorem at QFT.org.

  • A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, McGraw-Hill, 1971.
  • J. W. Negele and H. Orland, Quantum Many-Particle Systems, Addison-Wesley, 1988.
  • A. Altland and B. Simons, Condensed Matter Field Theory, 2nd ed., Cambridge University Press, 2010.
  • M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory, Addison-Wesley, 1995.
  • S. Weinberg, The Quantum Theory of Fields, Vol. I, Cambridge University Press, 1995.