Identical Particles and Many-Body Quantum Mechanics

Intuition [Beginner]

Pick up two electrons. Can you tell which is which? No -- and this is not a limitation of your measuring apparatus. It is a fundamental fact about nature. All electrons are truly identical: there is no label, no serial number, no hidden variable that distinguishes one from another. The same is true for photons, for protons, for every particle of a given type.

This indistinguishability has profound consequences. In classical physics, you could in principle paint one billiard ball red and the other blue and track them individually. In quantum mechanics, the wavefunction does not label particles. If you swap two identical particles, the physical state must be unchanged -- you cannot detect that a swap occurred. But "unchanged" in quantum mechanics means the probability density is unchanged, which allows the wavefunction itself to pick up a phase factor.

It turns out nature uses exactly two possible phase factors: $+1$ and $-1$.

• Bosons (integer spin: photons, pions, helium-4 atoms): the wavefunction is symmetric under exchange. $\psi (1,2)=+\psi (2,1)$. Bosons like to clump together. This is what makes lasers work (many photons in the same state) and what drives Bose-Einstein condensation.

• Fermions (half-integer spin: electrons, protons, neutrons, quarks): the wavefunction is antisymmetric under exchange. $\psi (1,2)=-\psi (2,1)$. Fermions refuse to share the same quantum state. This is the Pauli exclusion principle, and it explains the entire structure of the periodic table, the stability of matter, and the existence of white dwarf stars.

The connection between spin and exchange symmetry (integer spin $\to$ boson, half-integer spin $\to$ fermion) is the spin-statistics theorem, which can only be derived in relativistic quantum field theory.

For fermions, antisymmetry means if two particles occupy the same state, $\psi (1,2)=-\psi (2,1)=-\psi (1,2)$, which forces $\psi =0$. The state simply does not exist. This is the Pauli exclusion principle in its most fundamental form.

For many-particle systems, managing the symmetry constraints becomes combinatorially complex. Second quantization provides an elegant solution: instead of labeling particles and then symmetrizing/antisymmetrizing, you work with occupation numbers -- how many particles are in each single-particle state -- and creation/annihilation operators that add or remove particles while automatically enforcing the correct statistics.

Visual [Beginner]

EXCHANGE SYMMETRY
==================

  Bosons (symmetric):           Fermions (antisymmetric):

  psi(1,2) = +psi(2,1)         psi(1,2) = -psi(2,1)

    1   2      2   1              1   2      2   1
   (a) (b)    (a) (b)           (a) (b)    (a) (b)
    |   |      |   |              |   |      X   |
    +---+      +---+              +---+      +---+
  same sign                   sign flips on swap


PAULI EXCLUSION PRINCIPLE
==========================

  Two fermions in same state:

    psi(a,a) = -psi(a,a)
                |
    psi(a,a) = 0    (forbidden!)

  Allowed:                          Forbidden:
    +-  state 2                       +-  state 2
    | electron up                     | electron up + down
    +-  state 1                       +-  state 1
       electron up                       (would need both spins
                                         same state = antisymm = 0)


OCCUPATION NUMBER REPRESENTATION
=================================

  Single-particle states: |n1, n2, n3, ...>

  Bosons:  n_k = 0, 1, 2, 3, ...   (any number per state)
  Fermions: n_k = 0 or 1            (at most one per state)

  Example (3 particles, 4 states):

  Bosons:  |3,0,0,0> or |1,1,1,0> or |2,1,0,0>  ...
  Fermions: |1,1,1,0> or |1,1,0,1>  ...  (no state has n > 1)

Property	Bosons	Fermions
Spin	Integer ($0,1,2,\dots$)	Half-integer ($1/2,3/2,\dots$)
Exchange symmetry	$\psi (1,2)=+\psi (2,1)$	$\psi (1,2)=-\psi (2,1)$
Multiple occupancy	Allowed (any number)	Forbidden (Pauli exclusion)
Example	Photons, phonons, He-4	Electrons, protons, neutrons
Statistical distribution	Bose-Einstein	Fermi-Dirac

Worked example [Beginner]

Problem: Two non-interacting electrons are confined in a 1D infinite square well of width $a$. One is in the ground state $n=1$ and the other is in the first excited state $n=2$. Construct the total wavefunction including spin, and show that the spin-singlet and spin-triplet states lead to different spatial symmetries.

Solution:

Step 1: The single-particle spatial wavefunctions are ${\varphi }_n(x)=\sqrt{2/a}\sin (n\pi x/a)$.

Step 2: For two identical particles, we must form symmetric or antisymmetric spatial wavefunctions:

$${\psi }_S(x_1,x_2)=\frac{1}{\sqrt{2}}[{\varphi }_1(x_1){\varphi }_2(x_2)+{\varphi }_2(x_1){\varphi }_1(x_2)]$$

$${\psi }_A(x_1,x_2)=\frac{1}{\sqrt{2}}[{\varphi }_1(x_1){\varphi }_2(x_2)-{\varphi }_2(x_1){\varphi }_1(x_2)]$$

Step 3: For fermions, the total wavefunction (spatial $\times$ spin) must be antisymmetric. There are two ways to achieve this:

• Antisymmetric spatial + symmetric spin (triplet, $S=1$): $${\Psi }_{\text{total}}={\psi }_A(x_1,x_2)\otimes {\chi }_{S=1,M_S}$$ where ${\chi }_{1,1}=|\uparrow \uparrow ⟩$, ${\chi }_{1,0}=\frac{1}{\sqrt{2}}(|\uparrow \downarrow ⟩+|\downarrow \uparrow ⟩)$, ${\chi }_{1,-1}=|\downarrow \downarrow ⟩$.

• Symmetric spatial + antisymmetric spin (singlet, $S=0$): $${\Psi }_{\text{total}}={\psi }_S(x_1,x_2)\otimes {\chi }_{S=0}$$ where ${\chi }_{0,0}=\frac{1}{\sqrt{2}}(|\uparrow \downarrow ⟩-|\downarrow \uparrow ⟩)$.

Step 4: The spin-singlet state (with symmetric spatial part) has a non-zero probability of finding both electrons at the same location: $|{\psi }_S(x,x){|}^2>0$. The spin-triplet state (with antisymmetric spatial part) has $|{\psi }_A(x,x){|}^2=0$: the two electrons avoid each other. This is an effective "exchange interaction" -- not a real force, but a consequence of antisymmetry. It underlies ferromagnetism (Hund's rules).

Check your understanding [Beginner]

Formal definition [Intermediate+]

Indistinguishability postulate. For a system of $N$ identical particles, the Hilbert space is not the full tensor product ${\mathcal{H}}^{\otimes N}$ but rather the symmetric or antisymmetric subspace. The permutation operator $P_{ij}$ that exchanges particles $i$ and $j$ acts on physical states as:

$$P_{ij}|\psi ⟩=\pm |\psi ⟩$$

with $+1$ for bosons and $-1$ for fermions.

Symmetrized states. Given single-particle states $|\alpha ⟩,|\beta ⟩,\dots$, the properly symmetrized/antisymmetrized $N$-particle state is:

$$|\alpha ,\beta ,\dots {⟩}_{\pm }=\frac{1}{\sqrt{N!n_1!n_2!\cdots }}\sum _P(\pm 1{)}^{P}P|{\alpha }_{P(1)}⟩\otimes |{\alpha }_{P(2)}⟩\otimes \cdots$$

where the sum runs over all $N!$ permutations, $(\pm 1{)}^P$ is $+1$ for bosons and $\text{sgn}(P)$ for fermions, and $n_i$ are the occupation numbers (relevant for bosonic normalization).

Slater determinants. For fermions, the antisymmetrized state can be written as a Slater determinant:

$$\Psi (1,2,\dots ,N)=\frac{1}{\sqrt{N!}}\left|\begin{array}{ccc}{\psi }_a(1)& {\psi }_b(1)& \cdots \\ {\psi }_a(2)& {\psi }_b(2)& \cdots \\ ⋮& ⋮& \ddots \end{array}\right|$$

If any two single-particle states are identical, two columns are equal and the determinant vanishes -- this is the Pauli exclusion principle.

Helium atom. The Hamiltonian for the helium atom (two electrons, nuclear charge $Z=2$) is:

$$H=\left(-\frac{{\hslash }^2}{2m}{\nabla }_1^2-\frac{Z{e}^2}{4\pi {ϵ}_0{r}_1}\right)+\left(-\frac{{\hslash }^2}{2m}{\nabla }_2^2-\frac{Z{e}^2}{4\pi {ϵ}_0{r}_2}\right)+\frac{e^2}{4\pi {ϵ}_0|{\mathbf{r}}_1-{\mathbf{r}}_2|}$$

The last term is the electron-electron repulsion. In zeroth order (ignoring repulsion), the ground state has both electrons in the $1s$ orbital with opposite spins (singlet). First-order perturbation theory gives a ground-state energy correction. The exchange integral, arising from the antisymmetry requirement, leads to different energies for the singlet (para-helium) and triplet (ortho-helium) configurations.

Second quantization. The occupation number (Fock) representation labels states by the number of particles in each mode:

$$|n_1,n_2,n_3,\dots ⟩,n_k=0,1,2,\dots \text{ (bosons)};n_k=0,1\text{ (fermions)}$$

Creation ($a_k^{†}$) and annihilation ($a_k$) operators change occupation numbers. For bosons:

$$[a_k,a_{k^{\prime }}^{†}]={\delta }_{k{k}^{\prime }},[a_k,a_{k^{\prime }}]=[a_k^{†},a_{k^{\prime }}^{†}]=0$$

For fermions, the commutators are replaced by anticommutators:

$$\{a_k,a_{k^{\prime }}^{†}\}={\delta }_{k{k}^{\prime }},\{a_k,a_{k^{\prime }}\}=\{a_k^{†},a_{k^{\prime }}^{†}\}=0$$

The anticommutation relations automatically enforce the Pauli principle: $(a_k^{†}{)}^2=0$ for fermions (you cannot create two fermions in the same state).

The field operator creates/annihilates a particle at a spatial point:

$${\hat{\Psi }}^{†}(\mathbf{x})=\sum _k{\varphi }_k^{\ast }(\mathbf{x})a_k^{†}$$

where ${\varphi }_k(\mathbf{x})$ are the single-particle basis wavefunctions.

Many-body Hamiltonians are expressed in second-quantized form:

$$\hat{H}=\int {\hat{\Psi }}^{†}(\mathbf{x})h(\mathbf{x})\hat{\Psi }(\mathbf{x})d^{3}x+\frac{1}{2}\int \int {\hat{\Psi }}^{†}(\mathbf{x}){\hat{\Psi }}^{†}({\mathbf{x}}^{\prime })V(\mathbf{x}-{\mathbf{x}}^{\prime })\hat{\Psi }({\mathbf{x}}^{\prime })\hat{\Psi }(\mathbf{x})d^{3}x{d}^3{x}^{\prime }$$

where $h(\mathbf{x})$ is the single-particle Hamiltonian and $V$ is the two-body interaction.

Key results [Intermediate+]

1. Spin-statistics theorem. Particles with integer spin are bosons (symmetric wavefunctions); particles with half-integer spin are fermions (antisymmetric wavefunctions). This theorem requires relativistic quantum field theory for its proof; in nonrelativistic quantum mechanics it is assumed as a postulate.

2. Hartree-Fock approximation. The ground state of an $N$-fermion system is approximated as a single Slater determinant. Minimizing the energy leads to the Hartree-Fock equations:

$$h{\varphi }_i(\mathbf{x})+\sum _{j\ne i}\int \frac{|{\varphi }_j({\mathbf{x}}^{\prime }){|}^2{e}^2}{|\mathbf{x}-{\mathbf{x}}^{\prime }|}{\varphi }_i(\mathbf{x})d^3{x}^{\prime }-\sum _{j\ne i}{\delta }_{{\sigma }_i{\sigma }_j}\int \frac{{\varphi }_j^{\ast }({\mathbf{x}}^{\prime }){\varphi }_i({\mathbf{x}}^{\prime })e^2}{|\mathbf{x}-{\mathbf{x}}^{\prime }|}{\varphi }_j(\mathbf{x})d^3{x}^{\prime }={ϵ}_i{\varphi }_i(\mathbf{x})$$

The last term (exchange term) has no classical analog -- it arises purely from antisymmetry.

Advanced treatment [Master]

Correlation energy. The Hartree-Fock approximation captures the exchange energy exactly but misses the correlation energy -- the energy lowering due to correlated motion of electrons beyond what the mean field and exchange can account for. The correlation energy is defined as $E_c=E_{\text{exact}}-E_{\text{HF}}$ and is always negative (correlation always lowers the energy). Post-Hartree-Fock methods (configuration interaction, coupled cluster, Moller-Plesset perturbation theory) systematically recover this energy.

Green's functions and the self-energy. The single-particle Green's function:

$$G(\mathbf{x},t;{\mathbf{x}}^{\prime },t^{\prime })=-i⟨{\Psi }_0|T\{\hat{\Psi }(\mathbf{x},t){\hat{\Psi }}^{†}({\mathbf{x}}^{\prime },t^{\prime })\}|{\Psi }_0⟩$$

encodes the full spectrum and all single-particle properties. Dyson's equation relates the interacting Green's function to the noninteracting one through the self-energy $\Sigma$:

$$G=G_0+G_0\Sigma G$$

The self-energy contains all exchange and correlation effects. Diagrammatically, it is the sum of all one-particle irreducible diagrams.

Landau Fermi liquid theory. For a system of interacting fermions at low temperature, Landau postulated that the low-lying excitations are quasiparticles that behave like weakly interacting fermions, with a one-to-one correspondence to the noninteracting Fermi gas states. The quasiparticle effective mass $m^{\ast }$, the Landau parameters $F_{\ell }^s$, $F_{\ell }^a$, and the quasiparticle lifetime $\tau \propto (E-E_F{)}^{-2}$ characterize the system. This framework explains the specific heat, magnetic susceptibility, and transport properties of metals.

Bose-Einstein condensation. For noninteracting bosons in 3D, below a critical temperature $T_c$:

$$T_c=\frac{2\pi {\hslash }^2}{m{k}_B}{\left(\frac{n}{\zeta (3/2)}\right)}^{2/3}$$

a macroscopic fraction of particles occupies the ground state. The condensate fraction is $N_0/N=1-(T/T_c{)}^{3/2}$. For interacting bosons, Bogoliubov theory treats the condensate as a classical field (macroscopic occupation) plus small quantum fluctuations:

$$\hat{\Psi }(\mathbf{x})=\sqrt{n_0}+\delta \hat{\Psi }(\mathbf{x})$$

The fluctuations are diagonalized by Bogoliubov transformations, yielding quasiparticles with the spectrum $E_k=\sqrt{{ϵ}_k({ϵ}_k+2g{n}_0)}$ where $g$ is the interaction strength and ${ϵ}_k={\hslash }^2{k}^2/2m$. At low $k$, this is linear (phonons), not quadratic (free particles), demonstrating that interactions fundamentally change the excitation spectrum.

Quantum many-body entanglement. The ground state of a many-body system is typically highly entangled. For gapped local Hamiltonians in 1D, the area law states that the entanglement entropy of a subsystem scales with the boundary, not the volume. This underlies the success of matrix product states (MPS) and the density matrix renormalization group (DMRG) for 1D systems.

Jordan-Wigner transformation. In 1D, fermions can be mapped to spins:

$$a_j^{†}={\sigma }_j^{+}\prod _{k<j}{\sigma }_k^z,a_j={\sigma }_j^{-}\prod _{k<j}{\sigma }_k^z$$

This string operator accounts for the anticommutation relations. The resulting spin models (e.g., the transverse-field Ising model maps to free fermions) provide exactly solvable examples of many-body fermion systems.

Connections [Master]

• 12.05.01 -- Angular momentum and spin provide the single-particle quantum numbers that label the states in the occupation number representation.
• 12.08.01 -- Scattering amplitudes must be symmetrized/antisymmetrized for identical particles, modifying cross-sections.
• 12.10.01 -- Quantum statistical mechanics builds on the occupation number formalism to derive Bose-Einstein and Fermi-Dirac distributions.
• 12.12.01 -- Canonical QFT elevates the creation/annihilation operators to fundamental objects; particle number is no longer fixed but emerges from field quantization.
• 12.13.01 -- Quantum electrodynamics involves fermion (electron) and boson (photon) fields interacting via second-quantized operators.
• 09.07.01 -- Continuum mechanics fields are classical precursors to quantum fields; the transition from discrete to continuous degrees of freedom mirrors the many-body to QFT progression.

Bibliography [Master]

1. Fetter, A. L. and Walecka, J. D. Quantum Theory of Many-Particle Systems, Dover, 2003. The standard graduate text for many-body quantum mechanics with second quantization, Green's functions, and applications.

2. Abrikosov, A. A., Gorkov, L. P. and Dzyaloshinski, I. E. Methods of Quantum Field Theory in Statistical Physics, Dover, 1975. Classic text on diagrammatic methods for many-body systems.

3. Mattuck, R. D. A Guide to Feynman Diagrams in the Many-Body Problem, 3rd ed., Dover, 1992. An accessible introduction to diagrammatic techniques for many-body physics.

4. Negele, J. W. and Orland, H. Quantum Many-Particle Systems, Westview Press, 1998. Comprehensive treatment of path integrals, Green's functions, and diagrammatic methods.

5. Pitaevskii, L. and Stringari, S. Bose-Einstein Condensation and Superfluidity, Oxford University Press, 2016. Authoritative treatment of BEC theory and experiments.

6. Coleman, P. Introduction to Many-Body Physics, Cambridge University Press, 2015. Modern treatment emphasizing the renormalization group and emergent phenomena in many-body systems.