Skip to content

Identical Particles and Exchange Symmetry

9.1 Symmetrisation Postulate

For a system of NN identical particles, the wavefunction must satisfy:

ψ(,ri,,rj,)=±ψ(,rj,,ri,)\psi(\ldots, \mathbf{r}_i, \ldots, \mathbf{r}_j, \ldots) = \pm\psi(\ldots, \mathbf{r}_j, \ldots, \mathbf{r}_i, \ldots)

  • Bosons (integer spin): symmetric (++ sign). Any number can occupy the same state.
  • Fermions (half-integer spin): antisymmetric (- sign). Pauli exclusion: no two fermions can occupy the same state.

For two particles, the properly symmetrised states are:

ψS=12[ψa(1)ψb(2)+ψb(1)ψa(2)](bosons)\psi_S = \frac{1}{\sqrt{2}}\left[\psi_a(1)\psi_b(2) + \psi_b(1)\psi_a(2)\right] \quad \text{(bosons)}

ψA=12[ψa(1)ψb(2)ψb(1)ψa(2)](fermions)\psi_A = \frac{1}{\sqrt{2}}\left[\psi_a(1)\psi_b(2) - \psi_b(1)\psi_a(2)\right] \quad \text{(fermions)}

9.2 Exchange Interaction

Even without an explicit interaction potential, the requirement of (anti)symmetry leads to an effective exchange interaction. For two electrons in a box, the probability of finding them close together differs between the triplet (spatially antisymmetric, spin symmetric) and singlet (spatially symmetric, spin antisymmetric) states:

ψtriplet2=0whenr1=r2|\psi_{\text{triplet}|^2 = 0 \quad \text{when} \mathbf{r}_1 = \mathbf{r}_2}

ψsinglet2>0whenr1=r2|\psi_{\text{singlet}|^2 > 0 \quad \text{when} \mathbf{r}_1 = \mathbf{r}_2}

The triplet state keeps electrons apart (effective repulsion), while the singlet allows them to be close. This is the origin of the Hund”s first rule: parallel spins are energetically favourable for atoms because the exchange interaction lowers the Coulomb repulsion.

9.3 The Helium Atom

The helium Hamiltonian (ignoring nuclear motion):

H^=22me(12+22)2e24πε0r12e24πε0r2+e24πε0r1r2\hat{H} = -\frac{\hbar^2}{2m_e}\left(\nabla_1^2 + \nabla_2^2\right) - \frac{2e^2}{4\pi\varepsilon_0 r_1} - \frac{2e^2}{4\pi\varepsilon_0 r_2} + \frac{e^2}{4\pi\varepsilon_0|\mathbf{r}_1 - \mathbf{r}_2|}

Ground state (parahelium): Both electrons in the 1s1s orbital with opposite spins (singlet). The spatial part is symmetric: ψ100(r1)ψ100(r2)\psi_{100}(\mathbf{r}_1)\psi_{100}(\mathbf{r}_2).

First-order perturbation theory for the electron-electron repulsion:

E(1)=54e24πε0a0=52×13.6 eV=34.0 eVE^{(1)} = \frac{5}{4}\frac{e^2}{4\pi\varepsilon_0 a_0} = \frac{5}{2}\times 13.6\ \text{eV} = 34.0\ \text{eV}

The unperturbed ground state energy is E(0)=2×(54.4 eV)=108.8E^{(0)} = 2 \times (-54.4\ \text{eV}) = -108.8 eV (two electrons in Z=2Z = 2 Coulomb potential). Including perturbation: E108.8+34.0=74.8E \approx -108.8 + 34.0 = -74.8 eV. The experimental value is 79.0-79.0 eV.

Excited states: When one electron is excited to 1snl1s\,nlThe spin configuration matters:

  • Parahelium (singlet, S=0S = 0): symmetric spatial, antisymmetric spin. Lower energy for given configuration.
  • Orthohelium (triplet, S=1S = 1): antisymmetric spatial, symmetric spin. Higher energy.

The exchange integral KK and direct integral JJ:

J=ψa(1)2e24πε0r12ψb(2)2d3r1d3r2J = \iint |\psi_a(1)|^2\frac{e^2}{4\pi\varepsilon_0 r_{12}}|\psi_b(2)|^2\, d^3r_1 d^3r_2

K=ψa(1)ψb(2)e24πε0r12ψb(1)ψa(2)d3r1d3r2K = \iint \psi_a^*(1)\psi_b^*(2)\frac{e^2}{4\pi\varepsilon_0 r_{12}}\psi_b(1)\psi_a(2)\, d^3r_1 d^3r_2

The energy splitting between singlet and triplet is 2K2KWith the triplet lower by 2K2K.

Worked Example 9.1: Helium $1s2s$ States

For the 1s2s1s\,2s configuration of helium:

J1s,2s=e24πε0ψ1s(1)21r12ψ2s(2)2d3r1d3r2J_{1s,2s} = \frac{e^2}{4\pi\varepsilon_0}\int |\psi_{1s}(1)|^2\frac{1}{r_{12}}|\psi_{2s}(2)|^2\, d^3r_1 d^3r_2

K1s,2s=e24πε0ψ1s(1)ψ2s(2)1r12ψ2s(1)ψ1s(2)d3r1d3r2K_{1s,2s} = \frac{e^2}{4\pi\varepsilon_0}\int \psi_{1s}^*(1)\psi_{2s}^*(2)\frac{1}{r_{12}}\psi_{2s}(1)\psi_{1s}(2)\, d^3r_1 d^3r_2

Evaluating these (using the multipole expansion 1/r12=lr<l/r>l+1Pl(cosθ)1/r_{12} = \sum_l r_<^l/r_>^{l+1}\,P_l(\cos\theta)):

J1s,2s0.42 Ry=5.7 eVJ_{1s,2s} \approx 0.42\ \text{Ry} = 5.7\ \text{eV}

K1s,2s0.032 Ry=0.43 eVK_{1s,2s} \approx 0.032\ \text{Ry} = 0.43\ \text{eV}

The singlet (parahelium) has energy E=E0+J+KE = E_0 + J + KAnd the triplet (orthohelium) has E=E0+JKE = E_0 + J - K.

The splitting: EsingletEtriplet=2K0.86E_{\text{singlet} - E_{\text{triplet} = 2K \approx 0.86}} eV. This is the exchange splitting.

The orthohelium 23S2^3S state is metastable: it cannot decay to the ground state by electric dipole transition (because ΔS=0\Delta S = 0 for E1 transitions, and the ground state is a singlet). Its lifetime is 104\sim 10^4 s.

9.4 Slater Determinants

For NN fermions, the antisymmetric wavefunction is efficiently written as a Slater determinant:

Ψ(1,2,,N)=1N!ϕ1(1)ϕ2(1)ϕN(1)ϕ1(2)ϕ2(2)ϕN(2)ϕ1(N)ϕ2(N)ϕN(N)\Psi(1, 2, \ldots, N) = \frac{1}{\sqrt{N!}}\begin{vmatrix} \phi_1(1) & \phi_2(1) & \cdots & \phi_N(1) \\ \phi_1(2) & \phi_2(2) & \cdots & \phi_N(2) \\ \vdots & \vdots & \ddots & \vdots \\ \phi_1(N) & \phi_2(N) & \cdots & \phi_N(N) \end{vmatrix}

Properties:

  • Swapping any two rows (particles) changes the sign
  • If any two columns (orbitals) are identical, the determinant vanishes (Pauli exclusion)
  • The normalisation is correct if the spin-orbitals ϕi\phi_i are orthonormal