Identical Particles, Permutation Symmetry & Many-Body Systems

1. Quantum Indistinguishability

In classical mechanics, identical particles can always be distinguished by tracking their continuous trajectories through phase space. In quantum mechanics, the Heisenberg uncertainty principle ($\Delta x \Delta p \ge \hbar/2$) and the spatial overlap of wavepackets render identical particles fundamentally indistinguishable. No physical measurement can determine which particle is which.

2. The Permutation (Exchange) Operator

For a system of $N$ identical particles, the permutation operator $\hat{P}_{ij}$ interchanges all coordinates (spatial and spin) of particles $i$ and $j$:

$$\hat{P}_{ij}\psi(\xi_1, \dots, \xi_i, \dots, \xi_j, \dots, \xi_N) = \psi(\xi_1, \dots, \xi_j, \dots, \xi_i, \dots, \xi_N)$$

where $\xi_k = (\vec{r}_k, \sigma_k)$ represents the combined spatial and spin coordinates of particle $k$. Applying $\hat{P}_{ij}$ twice returns the original configuration:

$$\hat{P}_{ij}^2 = \hat{I} \implies \text{eigenvalues of } \hat{P}_{ij} \text{ are } \pm 1$$

Because the particles are physically identical, the Hamiltonian is completely invariant under particle interchange:

$$[\hat{P}_{ij}, \hat{H}] = 0 \quad \forall i, j$$

3. The Symmetrization Postulate & Spin-Statistics Theorem

The Symmetrization Postulate states that all physical quantum states must be either completely symmetric or completely antisymmetric under the interchange of any pair of identical particles:

1. Bosons (Integer Spin: $s = 0, 1, 2, \dots$): The state vector is strictly symmetric under interchange: $$\hat{P}_{ij}|\psi\rangle = +|\psi\rangle$$ Bosons obey Bose-Einstein statistics and can condense into the identical single-particle state (Bose-Einstein Condensation). Examples: photons ($s=1$), gluons ($s=1$), $^{4}\text{He}$ atoms ($s=0$).
2. Fermions (Half-Integer Spin: $s = 1/2, 3/2, \dots$): The state vector is strictly antisymmetric under interchange: $$\hat{P}_{ij}|\psi\rangle = -|\psi\rangle$$ Fermions obey Fermi-Dirac statistics and the Pauli Exclusion Principle. Examples: electrons ($s=1/2$), protons ($s=1/2$), neutrons ($s=1/2$).

The deep connection between spin and permutation symmetry is proven rigorously in relativistic quantum field theory via the Spin-Statistics Theorem (Pauli, 1940), stemming from causality, Lorentz invariance, and positive energy stability.

1. The Pauli Exclusion Principle

Consider a system of $N$ non-interacting fermions. If two fermions were to occupy the identical single-particle quantum state $\phi_\alpha$, interchanging them would yield:

$$\psi(\xi_1, \xi_2) = -\psi(\xi_2, \xi_1)$$

If $\xi_1 = \xi_2$, then $\psi(\xi_1, \xi_1) = -\psi(\xi_1, \xi_1) \implies \psi = 0$. Two identical fermions cannot occupy the exact same quantum state simultaneously—this is the Pauli Exclusion Principle.

2. The Slater Determinant

To construct a totally antisymmetric $N$-particle state from a set of single-particle orthonormal orbitals $\{\phi_{\alpha_1}, \phi_{\alpha_2}, \dots, \phi_{\alpha_N}\}$, John Slater introduced the Slater determinant:

$$\Psi(\xi_1, \xi_2, \dots, \xi_N) = \frac{1}{\sqrt{N!}} \begin{vmatrix} \phi_{\alpha_1}(\xi_1) & \phi_{\alpha_2}(\xi_1) & \cdots & \phi_{\alpha_N}(\xi_1) \\ \phi_{\alpha_1}(\xi_2) & \phi_{\alpha_2}(\xi_2) & \cdots & \phi_{\alpha_N}(\xi_2) \\ \vdots & \vdots & \ddots & \vdots \\ \phi_{\alpha_1}(\xi_N) & \phi_{\alpha_2}(\xi_N) & \cdots & \phi_{\alpha_N}(\xi_N) \end{vmatrix}$$

Mathematical properties:

1. Antisymmetry: Interchanging two particles corresponds to swapping two rows of the determinant, multiplying the determinant by $-1$.
2. Pauli Exclusion: If two orbitals are identical ($\alpha_i = \alpha_j$), two columns are identical, causing the determinant to vanish identically.
3. Normalization: The factor $1/\sqrt{N!}$ ensures $\langle \Psi | \Psi \rangle = 1$ when single-particle orbitals are orthonormal.

1. Total Wavefunction Factorization

For a system of two electrons (spin-1/2), the total wavefunction factors into a spatial part $\psi(\vec{r}_1, \vec{r}_2)$ and a spin part $\chi(\sigma_1, \sigma_2)$:

$$\Psi(1, 2) = \psi(\vec{r}_1, \vec{r}_2)\chi(\sigma_1, \sigma_2)$$

The total state must be antisymmetric under overall permutation $\hat{P}_{12} = \hat{P}_{\text{space}}\hat{P}_{\text{spin}} = -1$. As proven in Chapter 5:

• Spin Triplet ($S = 1$, $\chi_S$ symmetric): Requires an antisymmetric spatial wavefunction $\psi_A(\vec{r}_1, \vec{r}_2) = -\psi_A(\vec{r}_2, \vec{r}_1)$. This is termed Ortho-Helium.
• Spin Singlet ($S = 0$, $\chi_A$ antisymmetric): Requires a symmetric spatial wavefunction $\psi_S(\vec{r}_1, \vec{r}_2) = +\psi_S(\vec{r}_2, \vec{r}_1)$. This is termed Para-Helium.

2. The Direct and Exchange Integrals

Let the two electrons occupy single-particle spatial orbitals $\phi_a(\vec{r})$ and $\phi_b(\vec{r})$ ($a \ne b$). The normalized symmetric and antisymmetric spatial wavefunctions are:

$$\psi_{S/A}(\vec{r}_1, \vec{r}_2) = \frac{1}{\sqrt{2}}\left[ \phi_a(\vec{r}_1)\phi_b(\vec{r}_2) \pm \phi_b(\vec{r}_1)\phi_a(\vec{r}_2) \right]$$

Evaluating the expectation value of the Coulomb repulsion potential $V_{12} = \frac{e^2}{4\pi\varepsilon_0 |\vec{r}_1 - \vec{r}_2|}$:

$$\langle V_{12} \rangle_{S/A} = J \pm K$$

where:

1. $J \equiv \int |\phi_a(\vec{r}_1)|^2 \frac{e^2}{4\pi\varepsilon_0 r_{12}} |\phi_b(\vec{r}_2)|^2 \, d^3r_1 d^3r_2 > 0$ is the classical direct Coulomb integral (repulsion between charge clouds).
2. $K \equiv \int \phi_a^*(\vec{r}_1)\phi_b^*(\vec{r}_2) \frac{e^2}{4\pi\varepsilon_0 r_{12}} \phi_b(\vec{r}_1)\phi_a(\vec{r}_2) \, d^3r_1 d^3r_2 > 0$ is the purely quantum mechanical exchange integral.

3. Physical Consequences: Exchange Energy and Hund's Rule

The energy difference between singlet and triplet states is:

$$\Delta E = E_{\text{singlet}} - E_{\text{triplet}} = (E_0 + J + K) - (E_0 + J - K) = 2K > 0$$

The triplet state (ortho-helium) always lies lower in energy than the singlet state (para-helium). Because fermions in an antisymmetric spatial state have zero probability of occupying the same point in space ($\psi_A(\vec{r}, \vec{r}) = 0$), they are surrounded by a "Fermi exchange hole", keeping them further apart on average and dramatically reducing their electrostatic Coulomb repulsion energy. This is the microscopic physical origin of ferromagnetism and Hund's First Rule.

1. Statistical Mechanics of Non-Interacting Fermions in 3D

Consider a system of $N$ non-interacting spin-1/2 electrons confined in a volume $V = L^3$. Imposing periodic boundary conditions $\psi(x+L, y, z) = \psi(x, y, z)$, the allowed momentum states form a discrete grid in $k$-space:

$$\vec{k} = \frac{2\pi}{L}(n_x, n_y, n_z) \implies \Delta k^3 = \frac{(2\pi)^3}{V}$$

At absolute zero ($T = 0$), the electrons fill all available single-particle states from $k = 0$ up to a maximum momentum surface in $k$-space termed the Fermi sphere of radius $k_F$.

2. Fermi Momentum and Fermi Energy

Accounting for the spin degeneracy factor $g_s = 2s+1 = 2$:

$$N = 2 \frac{V}{(2\pi)^3} \left(\frac{4}{3}\pi k_F^3\right) = \frac{V k_F^3}{3\pi^2}$$

Solving for the Fermi wavevector $k_F$ in terms of the electron density $n = N/V$:

$$k_F = (3\pi^2 n)^{1/3}, \quad p_F = \hbar k_F = \hbar(3\pi^2 n)^{1/3}$$

The Fermi energy $E_F$ is the kinetic energy of the highest occupied state at $T = 0$:

$$E_F = \frac{\hbar^2 k_F^2}{2m} = \frac{\hbar^2}{2m}(3\pi^2 n)^{2/3}$$

The density of states $g(E) \equiv \frac{dN}{dE}$ evaluates to:

$$g(E) = \frac{V}{2\pi^2}\left(\frac{2m}{\hbar^2}\right)^{3/2}\sqrt{E}$$

3. Total Internal Energy and Quantum Degeneracy Pressure

The total ground-state kinetic energy $E_{\text{total}}$ is obtained by integrating over the Fermi sphere:

$$E_{\text{total}} = \int_0^{E_F} E g(E)\,dE = \frac{V}{2\pi^2}\left(\frac{2m}{\hbar^2}\right)^{3/2} \int_0^{E_F} E^{3/2}\,dE = \frac{V}{2\pi^2}\left(\frac{2m}{\hbar^2}\right)^{3/2} \frac{2}{5} E_F^{5/2} = \frac{3}{5} N E_F$$

Even at absolute zero temperature, quantum fermions possess enormous kinetic energy. The system exerts a macroscopic outward quantum degeneracy pressure:

$$P = -\left(\frac{\partial E_{\text{total}}}{\partial V}\right)_N = -\frac{3}{5}N \frac{\partial E_F}{\partial V} = -\frac{3}{5}N \left(-\frac{2}{3}\frac{E_F}{V}\right) = \frac{2}{5}n E_F = \frac{(3\pi^2)^{2/3}\hbar^2}{5m} n^{5/3}$$

Degeneracy pressure is purely quantum mechanical in origin (independent of temperature and electrostatic forces), providing the outward mechanical force that stabilizes white dwarf stars against gravitational collapse (up to the Chandrasekhar mass limit).

1. Minimal Coupling and the Landau Gauge

Consider an electron of mass $m$ and charge $q = -e$ confined to the $xy$-plane in a uniform perpendicular magnetic field $\vec{B} = B\hat{z}$. Under minimal electromagnetic coupling, the kinetic momentum operator is $\vec{\Pi} = \vec{p} - q\vec{A} = \vec{p} + e\vec{A}$. The Hamiltonian is:

$$\hat{H} = \frac{1}{2m}(\hat{\vec{p}} + e\vec{A})^2$$

Choosing the Landau gauge $\vec{A} = (0, Bx, 0)$, which yields $\nabla \times \vec{A} = B\hat{z}$:

$$\hat{H} = \frac{\hat{p}_x^2}{2m} + \frac{(\hat{p}_y + eBx)^2}{2m}$$

2. Harmonic Mapping and Landau Energy Spectrum

Since $\hat{H}$ does not contain $y$, the momentum operator $\hat{p}_y$ commutes with the Hamiltonian: $[\hat{H}, \hat{p}_y] = 0$. We can replace $\hat{p}_y$ by its continuous eigenvalue $\hbar k_y$:

$$\hat{H} = \frac{\hat{p}_x^2}{2m} + \frac{1}{2}m\omega_c^2 \left(x + \frac{\hbar k_y}{eB}\right)^2$$

where $\omega_c = \frac{eB}{m}$ is the classical cyclotron frequency. This Hamiltonian is mathematically identical to a 1D quantum harmonic oscillator shifted in equilibrium position to:

$$x_0 = -\frac{\hbar k_y}{eB} = -k_y l_B^2$$

where $l_B \equiv \sqrt{\frac{\hbar}{eB}}$ is the magnetic length.

The energy eigenvalues are the famous discrete Landau levels:

$$E_n = \hbar\omega_c\left(n + \frac{1}{2}\right) \quad (n = 0, 1, 2, \dots)$$

3. Macroscopic Degeneracy of Landau Levels

Notice that the energy $E_n$ is completely independent of $k_y$. For a sample of dimensions $L_x \times L_y$, the center of the harmonic oscillator must lie within the physical boundaries of the sample ($0 \le x_0 \le L_x$):

$$0 \le \frac{\hbar k_y}{eB} \le L_x \implies 0 \le k_y \le \frac{eBL_x}{\hbar}$$

With periodic boundary conditions along $y$, $k_y = \frac{2\pi n_y}{L_y}$, the number of allowed $k_y$ values in this interval is:

$$g = \frac{\Delta k_y}{2\pi / L_y} = \frac{eBL_x}{\hbar}\frac{L_y}{2\pi} = \frac{eB(L_x L_y)}{2\pi\hbar} = \frac{\Phi}{\Phi_0}$$

where $\Phi = B(L_x L_y)$ is the total magnetic flux through the sample and $\Phi_0 = \frac{h}{e}$ is the magnetic flux quantum. Each Landau level possesses a gigantic macroscopic degeneracy equal to the number of flux quanta piercing the sample, forming the physical basis of the Integer Quantum Hall Effect.

1. Nuclear Charge Independence and the Isospin Formalism

High-energy scattering experiments reveal that the strong nuclear interaction is charge-independent: the nuclear force between two protons ($p-p$), two neutrons ($n-n$), and a proton and neutron ($p-n$) in identical spatial and spin states is virtually identical.

To exploit this symmetry, Werner Heisenberg introduced the isospin formalism. The proton and neutron are viewed as two orthogonal isospin states of a single particle, the nucleon, with isospin $I = 1/2$:

$$|p\rangle = |1/2, +1/2\rangle = \begin{pmatrix} 1 \\ 0 \end{pmatrix}, \quad |n\rangle = |1/2, -1/2\rangle = \begin{pmatrix} 0 \\ 1 \end{pmatrix}$$

The electric charge of a nucleon is given by the Gell-Mann-Nishijima formula:

$$Q = e\left(I_3 + \frac{1}{2}\right)$$

2. The Two-Nucleon System and the Deuteron

Coupling the isospins of two nucleons ($I_1 = 1/2, I_2 = 1/2$):

$$\frac{1}{2} \otimes \frac{1}{2} = 1 \oplus 0$$

• Isospin Triplet ($I = 1$, symmetric in isospin): Includes $|1, 1\rangle = |pp\rangle$, $|1, -1\rangle = |nn\rangle$, and $|1, 0\rangle = \frac{1}{\sqrt{2}}(|pn\rangle + |np\rangle)$.
• Isospin Singlet ($I = 0$, antisymmetric in isospin): $|0, 0\rangle = \frac{1}{\sqrt{2}}(|pn\rangle - |np\rangle)$.

Under the Generalized Pauli Principle, the total state of two nucleons must be completely antisymmetric under overall interchange of spatial, spin, and isospin coordinates:

$$\hat{P}_{\text{space}} \hat{P}_{\text{spin}} \hat{P}_{\text{isospin}} = (-1)^L (-1)^{S+1} (-1)^{I+1} = -1 \implies (-1)^{L + S + I} = -1$$

For the bound ground state of the deuteron (a proton-neutron bound state): experimentally, $L = 0$ (mostly $s$-wave) and total spin $J = 1$ ($S = 1$). Substituting into the condition:

$$(-1)^{0 + 1 + I} = -1 \implies (-1)^I = +1 \implies I = 0$$

The deuteron bound state is an isospin singlet ($I = 0$). This explains why there are no bound diproton ($pp$) or dineutron ($nn$) states in nature: both require $I = 1$, which has higher energy due to spin-dependent nuclear tensor forces!