Spin-Statistics Theorem

Every article in this chapter has leaned on one unexplained rule: half-integer-spin particles take the antisymmetric sector, integer-spin particles the symmetric one. Electrons are fermions, photons are bosons — and nonrelativistic quantum mechanics simply accepts this as input; the Schrödinger equation would be perfectly happy with symmetric electrons. Where does the rule come from? The spin-statistics theorem belongs to relativistic quantum field theory: in any Lorentz-invariant local theory with an energy floor and healthy probabilities, integer-spin fields must be quantized as bosons and half-integer-spin fields as fermions. We will not pretend to prove that theorem here — its full proof needs machinery this track builds much later. What we will do is exhibit both catastrophes exactly, in the two smallest laboratories that show them: a 2×2 mode-pair energy audit in which a fermionic spin-0 field or a bosonic spin-½ field loses its energy floor, and a causality audit in which the spin-0 equal-time commutator function vanishes identically outside the light cone — we prove the zero symbolically — while the anticommutator function, computed exactly in closed form, does not: quantize a scalar with anticommutators and its correlations leak outside the light cone. Both computations are verified numerically to eight digits. This article also settles a debt: two earlier articles, Rotation Operators and SU(2) and Spin-1/2 and the Pauli Matrices, explicitly defer the spin-statistics connection here — and this is where it is paid.

Key takeaways

In relativistic quantum field theory, integer-spin fields must be quantized with commutators (bosons) and half-integer-spin fields with anticommutators (fermions) — forced by four hypotheses (Lorentz invariance, microcausality, an energy floor, positive norms). The wrong choices fail visibly: a fermionic scalar's correlations leak outside the light cone (Δ1(r,0)=mK1(mr)/2π2r≠0\Delta_1(r,0)=mK_1(mr)/2\pi^2r\neq0), and a bosonic Dirac field's energy has no floor. Nonrelativistic quantum mechanics cannot prove the theorem — there the boson/fermion split is an empirical postulate, and the popular 2π2\pi-rotation argument is a mnemonic, not a proof. Every downstream consequence — exclusion, shells, degeneracy pressure, condensation — inherits its authority from this one relativistic result (Fierz 1939, Pauli 1940).

The core insight

Lorentz invariance + microcausality + energy bounded below force integer-spin fields to be bosons and half-integer-spin fields to be fermions — the boson/fermion postulate of nonrelativistic QM is a theorem of relativistic QFT.

Where this lives today

The theorem is under permanent experimental audit, because any violation would falsify quantum field theory itself. The cleanest tests are statistics tests: spectroscopy on ⁴He and ³He — one a boson, one a fermion, their wildly different low-temperature fates dissected in Bosons and Fermions — confirms the composite rule; searches for Pauli-forbidden atomic transitions (VIP-2 at Gran Sasso listens for X-rays from copper electrons dropping into filled shells) keep returning null. On the theory side, the dimensional fine print became a field: in two-dimensional electron gases the braid-group loophole is real, and the fractional-quantum-Hall anyons that exploit it were directly braided in 2020-era interferometry, with topological quantum computing betting on their non-Abelian cousins. The (−1)2j(-1)^{2j} ingredient itself is routine engineering — every qubit lab that drives a 2π2\pi rotation and watches the sign flip re-measures the theorem's key input (the fringe experiment lives in Rotation Operators and SU(2)).

The physical picture

Picture a quantum field as a machine distributed over space that can create or destroy one particle at each point. Relativity imposes a non-negotiable discipline: an intervention here must not disturb any measurement there if "there" lies outside the light cone. In operator language, observables built from the field at spacelike-separated points must commute — and whether they do depends on a single sign choice at quantization (commutators or anticommutators), which the field's spin decides.

Figure
Two spacelike points and the discipline of the light cone
t x light cone of x event x event y — spacelike right statistics: field functions at x and y cancel exactly — Δ(x−y) = 0 outside the cone wrong statistics: Δ₁(x−y) ≠ 0 leaks across — spacelike interventions interfere

Think of the two failure modes as two cliffs with one ridge between them. Cliff one is energy: quantize the wrong way and filling certain modes lowers the total energy without limit — no ground state to stand on. Cliff two is causality: the field's correlations refuse to vanish outside the light cone. Integer spin falls off cliff one under anticommutators; half-integer spin falls off cliff two (and one) under commutators. The single consistent path is the theorem, and the 2×2 audit is Derivation B below. The next figure is the article's punchline, and it is fully quantitative. For a free scalar field of mass mm we can compute both candidate equal-time correlation structures in closed form. The commutator function Δ(r,0)\Delta(r,0) — the one Bose quantization actually uses — vanishes identically at any spatial separation r>0r>0: we prove the zero symbolically below. The anticommutator function Δ1(r,0)\Delta_1(r,0) — the structure a hypothetical Fermi scalar would be built on — comes out as mK1(mr)/2π2rmK_1(mr)/2\pi^2 r: exponentially small beyond the Compton wavelength, but never zero. The wrong choice does not fail loudly; it fails by a whisper that never dies, and a whisper across the light cone is enough to wreck relativistic causality.

Figure
The leak outside the light cone — Δ vanishes exactly and Δ₁ never does
spatial separation r in Compton units r/λ̄_C log₁₀ of Δ₁(r)/Δ₁(λ̄_C) 0 −2 −4 1 2 3 4 5 6 Δ₁(r,0) = mK₁(mr)/2π²r — the Fermi-scalar leak, e^(−r/λ̄) small but never zero Δ(r,0) = 0 exactly, every r — the Bose choice (off the log scale) Compton scale λ̄_C = ħ/mc ×0.116 per λ̄ at r=2 · ×0.00134 at r=5

Crucially, the theorem is relativistic. Nonrelativistic quantum mechanics is perfectly consistent with either statistics for either spin — Permutation Symmetry and the Exchange Operator built both sectors without ever asking about spin. Experiment picks the pairing, and only relativistic field theory explains the pick. The rule you have been using all chapter is the shadow of a light-cone argument.

The theorem, stated with its hypotheses

Before any equation, fix what is being claimed. The theorem is a statement about fields — the operator machines of relativistic quantum field theory — and every hypothesis is load-bearing. In a quantum field theory with (i) Lorentz-invariant dynamics, (ii) microcausality — observables built from fields at spacelike-separated points are compatible, (iii) a Hamiltonian bounded below with a stable vacuum, and (iv) a positive-norm state space, the following pairing is forced:

s∈{0,1,2,… }  ⟹  Bose (commutators),s∈{12,32,… }  ⟹  Fermi (anticommutators).(1)s \in \{0,1,2,\dots\} \;\Longrightarrow\; \text{Bose (commutators)}, \qquad s \in \{\tfrac12,\tfrac32,\dots\} \;\Longrightarrow\; \text{Fermi (anticommutators)}. \qquad (1)

where: ss — the spin of the field, i.e. which representation of the rotation (Lorentz) group it transforms under; Bose/Fermi — whether mode operators obey [a,a†]=1[a,a^\dagger]=1 or {a,a†}=1\{a,a^\dagger\}=1, exactly the symmetric/antisymmetric sector split of Symmetric and Antisymmetric States. Each hypothesis has a job: drop Lorentz invariance and nonrelativistic mechanics happily hosts symmetric electrons; drop microcausality and the Fermi scalar below becomes legal; drop the energy floor and the Bose–Dirac disaster is admissible; drop positive norms and ghosts fake any statistics. In two spatial dimensions the exchange topology itself changes and the conclusion genuinely relaxes — that loophole is Anyons and Fractional Statistics' story.

The two correlation functions of a free scalar

The cleanest laboratory is the free real scalar field of mass mm — spin 0, no indices. Expand it in modes (the standard construction, quoted as the starting point, ℏ=c=1\hbar=c=1):

ϕ(x)=∫ ⁣d3p(2π)3 2Ep(ap e−ip⋅x+ap† e+ip⋅x),Ep=∣p∣2+m2.(2)\phi(x)=\int\!\frac{\mathrm{d}^3p}{(2\pi)^3\,2E_p}\Big(a_{\mathbf p}\,e^{-ip\cdot x}+a_{\mathbf p}^\dagger\,e^{+ip\cdot x}\Big), \qquad E_p=\sqrt{|\mathbf p|^2+m^2}. \qquad (2)

where: ap,ap†a_{\mathbf p},a_{\mathbf p}^\dagger — annihilation/creation operators for momentum p\mathbf p; p⋅x=Ept−p⋅xp\cdot x = E_p t-\mathbf p\cdot\mathbf x — the Minkowski phase. Physically: ϕ(x)\phi(x) superposes "destroy a particle here" and "create a particle here" — the combination relativity itself forces, since a purely positive-frequency field cannot be local. Everything about statistics now sits in one question: does the field-pair bracket hand you the difference or the sum of the two frequency exponentials? Bose commutators give the difference — the commutator function; Fermi anticommutators give the sum. At equal times (tx=tyt_x=t_y, spatial separation rr) they are:

Δ(r,0)= ⁣∫ ⁣d3p(2π)32Ep(eip⋅r−e−ip⋅r),Δ1(r,0)= ⁣∫ ⁣d3p(2π)32Ep(eip⋅r+e−ip⋅r).(3)\Delta(r,0)=\!\int\!\frac{\mathrm{d}^3p}{(2\pi)^3 2E_p}\Big(e^{i\mathbf p\cdot\mathbf r}-e^{-i\mathbf p\cdot\mathbf r}\Big), \qquad \Delta_1(r,0)=\!\int\!\frac{\mathrm{d}^3p}{(2\pi)^3 2E_p}\Big(e^{i\mathbf p\cdot\mathbf r}+e^{-i\mathbf p\cdot\mathbf r}\Big). \qquad (3)

where: Δ\Delta — the structure [ϕ(x),ϕ(y)][\phi(x),\phi(y)] delivers under Bose quantization; Δ1\Delta_1 — the structure {ϕ(x),ϕ(y)}\{\phi(x),\phi(y)\} delivers under the (hypothetical) Fermi quantization; equal times is the most spacelike slice, so it is the decisive test. Read the sign: the difference can cancel, the sum cannot. That one sign is the entire causality half of the theorem for spin 0, and Derivation A does both integrals exactly. The sum evaluates in closed form to Δ1(r,0)=mK1(mr)/2π2r\Delta_1(r,0)=mK_1(mr)/2\pi^2 r (Derivation A, verified symbolically and numerically), with K1K_1 the modified Bessel function and mr=r/λ ⁣ ⁣ˉCmr=r/\lambda\!\!\bar{}_C — the punchline curve of the figure above. It is positive for every rr: exponentially suppressed beyond the Compton wavelength (∼e−mr\sim e^{-mr}), yet never zero. A Fermi-quantized scalar would carry irreducible operator structure connecting spacelike-separated points — measurements outside each other's light cones could not be independently prescribed.

The energy audit — one mode pair is enough

The second cliff needs even less machinery: a single (particle, antiparticle) mode pair of energy EE. Write the raw, unordered Hamiltonian each field theory delivers before any statistics choice is imposed. For the scalar, the energy density is a sum of squares, so both sectors enter positively; for the Dirac field, the antiparticle spin-sum flips the sign of its term — we quote that sign with its source (the vv-spinor completeness relation; Peskin & Schroeder §3.5) rather than re-deriving spinor algebra here:

Hscalarraw=E(a†a+b b†),HDiracraw=E(a†a−b b†).(5)H_{\text{scalar}}^{\text{raw}}=E\big(a^\dagger a+b\,b^\dagger\big), \qquad H_{\text{Dirac}}^{\text{raw}}=E\big(a^\dagger a-b\,b^\dagger\big). \qquad (5)

where: aa — particle mode, bb — antiparticle mode; the ordering b b†b\,b^\dagger is what the mode expansion literally produces. Everything now hinges on the reordering rule the statistics choice supplies: Bose gives b b†=b†b+1b\,b^\dagger=b^\dagger b+1, Fermi gives b b†=1−b†bb\,b^\dagger=1-b^\dagger b. Substituting each rule into each raw Hamiltonian generates the four-entry audit of the second figure — two constants shift harmlessly, one spectrum keeps a floor, and two lose it. Derivation B runs the table exactly. The mode-level table is exact and honest about its scope: it demonstrates the mechanism by which the wrong pairing destroys the energy floor, mode by mode. The general theorem — that no rescue is possible for any interacting, any-spin variant — is the axiomatic result of Fierz, Pauli, and the 1950s formalization; we cite, not derive.

Derivations

Derivation AThe equal-time commutator vanishes exactly and the anticommutator does not

The goal is the causality half of the audit for spin 0: compute both candidate correlation structures of equation (3) at equal times, in closed form, and establish the sharp asymmetry — the Bose structure Δ(r,0)\Delta(r,0) is identically zero for every spatial separation, while the Fermi structure Δ1(r,0)\Delta_1(r,0) is strictly positive everywhere. The strategy is elementary: reduce each three-dimensional momentum integral to a radial one by doing the angular integral first, watch the two frequency exponentials produce identical radial kernels, and then either subtract them (exact zero) or add them (a known Bessel integral). The only care needed is convergence: the radial integrand decays like sin⁡(pr)\sin(pr) times a slowly varying function, so the sum case is a conditionally convergent oscillatory integral — which is why we verify it two independent ways.

Step 1. Angular integral. For either exponential, put u=cos⁡θu=\cos\theta between p\mathbf p and r\mathbf r; then ∫−11e±ipru du=2sin⁡(pr)/(pr)\int_{-1}^{1}e^{\pm i p r u}\,\mathrm{d}u = 2\sin(pr)/(pr) — the same kernel for both signs, because flipping p→−p\mathbf p\to-\mathbf p maps one exponential onto the other while d3p/Ep\mathrm{d}^3p/E_p is invariant. Why this is the whole story for the commutator: the two terms of Δ(r,0)\Delta(r,0) reduce to byte-identical radial integrals, so their difference is zero before any radial integration is attempted. The vanishing is symmetry, not cancellation-by-miracle.

Step 2. Assemble the commutator. Both terms equal 14π2∫0∞p dpEp sin⁡(pr)r\frac{1}{4\pi^2}\int_0^\infty \frac{p\,\mathrm{d}p}{E_p}\,\frac{\sin(pr)}{r}, hence

Δ(r,0)=0for every r>0,\Delta(r,0)=0 \qquad\text{for every } r>0,

exactly. Equal-time separations are spacelike, and a boost carries this zero across the whole spacelike region (each term is a Lorentz-invariant function of (x−y)2(x-y)^2) — the Bose-quantized scalar is causal.

Step 3. Assemble the anticommutator. The sum doubles the same kernel:

Δ1(r,0)=12π2r∫0∞p sin⁡(pr)p2+m2 dp.\Delta_1(r,0)=\frac{1}{2\pi^2 r}\int_0^\infty \frac{p\,\sin(pr)}{\sqrt{p^2+m^2}}\,\mathrm{d}p.

Step 4. Evaluate by the standard Bessel representation. Since ∫0∞cos⁡(pr)/p2+m2 dp=K0(mr)\int_0^\infty \cos(pr)/\sqrt{p^2+m^2}\,\mathrm{d}p=K_0(mr), differentiating with respect to rr gives ∫0∞psin⁡(pr)/p2+m2 dp=mK1(mr)\int_0^\infty p\sin(pr)/\sqrt{p^2+m^2}\,\mathrm{d}p=mK_1(mr), so

Δ1(r,0)=m K1(mr)2π2 r>0for every r>0.\Delta_1(r,0)=\frac{m\,K_1(mr)}{2\pi^2\,r}>0 \qquad\text{for every } r>0.

Why positivity kills the Fermi scalar: if the field's basic bracket is the anticommutator, this non-vanishing spacelike structure is irreducible — no operator ordering removes it, and local measurements at spacelike separation cease to be independently specifiable. The suppression scale is the Compton wavelength: the leak dies like e−mre^{-mr} but never reaches zero. *Verified: sympy reduces both angular integrals to the identical 2sin⁡(pr)/pr2\sin(pr)/pr kernel (difference exactly 0); the radial integral is computed numerically (oscillatory-weighted quadrature) and agrees with mK1(mr)/2π2rmK_1(mr)/2\pi^2r to 3×10⁻⁹ relative at r/λ̄ = 0.5–5; all quoted values recomputed.*

Go deeper · Derivation B · 2 minThe 2 by 2 energy audit — which pairings keep a floor
Derivation BThe 2 by 2 energy audit — which pairings keep a floor

The goal is the energy half of the audit, run exactly on one particle–antiparticle mode pair: substitute each statistics reordering rule into each raw Hamiltonian of equation (5) and read off the four spectra. The strategy is bookkeeping made honest — every spectrum is assembled on the exact occupation-number basis, with the fermionic rules first proven as operator identities on explicit 2×2 matrices, and the sign that distinguishes scalar from Dirac quoted from the spinor completeness sums with its source rather than smuggled in. Watch for the two distinct ways a floor can die: per-mode collapse (each filled antiparticle mode subtracts energy, so infinitely many modes dig an infinite hole) and per-occupation collapse (one mode alone runs away as its occupation grows).

Step 1. Certify the reordering rules. On the two-state fermion mode, {a,a†}=1\{a,a^\dagger\}=\mathbb 1 and a2=0a^2=0 hold as matrix identities, giving Fermi's rule b b†=1−b†bb\,b^\dagger=1-b^\dagger b; the bosonic ladder gives Bose's rule b b†=b†b+1b\,b^\dagger=b^\dagger b+1.

Step 2. Scalar with Bose (the right pairing): H=E(na+nb+1)H=E(n_a+n_b+1). Floor +E+E at na=nb=0n_a=n_b=0; every excitation costs +E+E. Consistent.

Step 3. Scalar with Fermi (wrong): H=E(na+(1−nb))H=E\big(n_a+(1-n_b)\big). The spectrum on {na,nb}∈{0,1}2\{n_a,n_b\}\in\{0,1\}^2 is {0,E,E,2E}\{0,E,E,2E\} — filling the antiparticle mode lowers the energy by EE. Why this is fatal even though four states look bounded: the field has a mode pair at every momentum, and each filled bb-mode subtracts its own EpE_p. Two modes with E1=1,E2=3E_1=1,E_2=3 already drop the energy by 4 below the b-vacuum; the sum over all momenta is −∑pEp-\sum_p E_p, bottomless. The floor dies by a thousand cuts.

Step 4. Dirac with Fermi (the right pairing): H=E(na−(1−nb))=E(na+nb)−EH=E\big(n_a-(1-n_b)\big)=E(n_a+n_b)-E. Spectrum {−E,0,0,E}\{-E,0,0,E\}: a harmless constant shift, absorbed into the vacuum energy; excitations cost +E+E. Consistent — this sign repair is precisely why the Dirac sea never materializes in the modern formalism.

Step 5. Dirac with Bose (wrong): H=E(na−(nb+1))H=E\big(n_a-(n_b+1)\big). Now nbn_b is bosonic — unbounded — and the spectrum contains −(nb+1)E-(n_b+1)E for every nbn_b: floors −6E,−11E,−21E-6E,-11E,-21E at nb=5,10,20n_b=5,10,20 and no bottom. One single mode destroys the theory. Why one mode suffices here but not in Step 3: bosonic occupation is unlimited, so the runaway needs no help from other modes — the catastrophe is visible in a two-level table, which is why Pauli could run this argument in 1940 without a computer. *Verified: fermionic identities {a,a†}=1, a²=0, n²=n checked on explicit matrices; all four spectra assembled on exact number bases — floors +E, per-mode −E loss, −E constant, and −(n_b+1)E → −∞ reproduced; boson commutator certified on a 30-level ladder.*

The bridge back to the antisymmetric sector

The two derivations settle the physics; one line settles the bookkeeping that makes the theorem matter for this chapter. Fermionic creation operators anticommute, {a^i†,a^j†}=0\{\hat a_i^\dagger,\hat a_j^\dagger\}=0, so two-particle states are automatically antisymmetric, a^i†a^j†∣0⟩=−a^j†a^i†∣0⟩\hat a_i^\dagger\hat a_j^\dagger|0\rangle=-\hat a_j^\dagger\hat a_i^\dagger|0\rangle, and coincident labels give (a^i†)2=0(\hat a_i^\dagger)^2=0 — no double occupancy, the seed Pauli Exclusion Principle grows into the periodic table and the white dwarf. Bosonic operators commute and populate the symmetric sector the same way. The relativistic sign choice thereby selects between the two sectors that Permutation Symmetry and the Exchange Operator proved were the only options — it never re-derives them, and the nonrelativistic postulate (fermions antisymmetric, bosons symmetric) is inherited, not assumed. (Verified: (a†)2=0(a^\dagger)^2=0 and n2=nn^2=n are exact matrix identities in the verify script.)

Limiting cases

caseresultmeaning
r≫λ ⁣ ⁣ˉCr \gg \lambda\!\!\bar{}_CΔ1∼e−mr\Delta_1 \sim e^{-mr}, still ≠0\neq 0the acausal leak of wrong statistics is exponentially small — and still fatal, because exactly zero is what causality demands
r→λ ⁣ ⁣ˉCr \to \lambda\!\!\bar{}_C and belowΔ1=O(m2)\Delta_1 = O(m^2)within a Compton wavelength the leak is order one: a Fermi scalar is not even approximately causal at its own natural scale
m→∞m\to\inftyΔ1→0\Delta_1\to0 pointwise for r>0r>0an infinitely heavy wrong-statistics field hides its sickness — the nonrelativistic limit, where the theorem loses its teeth
nonrelativistic QMeither statistics formally consistentno light cone, no energy-positivity constraint on antiparticles: the pairing must be postulated, as this chapter has done
two spatial dimensionsexchange phase eiθe^{i\theta} allowedthe configuration-space topology changes (braid group replaces SNS_N): Anyons and Fractional Statistics — the theorem's honest fine print
spin ladder upwardalternation persistsspin 1 (photon): Bose; spin 32\tfrac32: Fermi; spin 2 (graviton): Bose — the audit's pattern holds for every rung

Worked examples

Example 1 — how big is the causality leak of a Fermi-quantized electron-mass scalar? Task: a hypothetical spin-0 field with the electron's mass is quantized with anticommutators. Quantify Δ1(r,0)\Delta_1(r,0) at 1, 2 and 5 reduced Compton wavelengths (λ ⁣ ⁣ˉC=ℏ/mec=0.386 pm\lambda\!\!\bar{}_C=\hbar/m_ec=0.386\,\mathrm{pm}), in units of m2m^2.

Step 1. Enumerate. The closed form is Δ1(nλ ⁣ ⁣ˉC,0)=m2K1(n)/2π2n\Delta_1(n\lambda\!\!\bar{}_C,0)=m^2K_1(n)/2\pi^2n; everything reduces to Bessel values K1(1)=0.60191K_1(1)=0.60191, K1(2)=0.13987K_1(2)=0.13987, K1(5)=4.0446×10−3K_1(5)=4.0446\times10^{-3}.

Step 2. Apply. At n=1n=1: Δ1=0.60191/(2π2)=0.030493 m2\Delta_1=0.60191/(2\pi^2)=\mathbf{0.030493\,m^2}. At n=2n=2: 0.13987/(4π2)=0.0035430 m20.13987/(4\pi^2)=\mathbf{0.0035430\,m^2}. At n=5n=5: 4.0446×10−3/(10π2)=4.098×10−5 m24.0446\times10^{-3}/(10\pi^2)=\mathbf{4.098\times10^{-5}\,m^2}.

Step 3. Boundary number — the suppression ladder. Relative to one Compton wavelength, the leak falls to 11.6%11.6\% at 2λ ⁣ ⁣ˉC2\lambda\!\!\bar{}_C and 0.134%0.134\% at 5λ ⁣ ⁣ˉC5\lambda\!\!\bar{}_C. At atomic distances (r≈137 λ ⁣ ⁣ˉCr\approx137\,\lambda\!\!\bar{}_C) the leak is ∼e−137\sim e^{-137}: hopeless to measure, yet nonzero — and a single nonzero value anywhere outside the cone already contradicts the causality axiom. The wrong theory is not approximately right; it is exactly wrong.

Example 2 — the energy books for one mode pair, all four ways. Task: a particle–antiparticle mode pair has E=1E=1 (units of mc2mc^2). Tabulate the exact spectra of the four spin-statistics pairings and identify each verdict.

Step 1. Enumerate the raw Hamiltonians and reordering rules: Hscalarraw=E(a†a+bb†)H^{\text{raw}}_{\text{scalar}}=E(a^\dagger a+bb^\dagger), HDiracraw=E(a†a−bb†)H^{\text{raw}}_{\text{Dirac}}=E(a^\dagger a-bb^\dagger); Bose bb†=nb+1bb^\dagger=n_b+1, Fermi bb†=1−nbbb^\dagger=1-n_b.

Step 2. Apply each rule to each raw form:

pairingspectrumfloorverdict
scalar + BoseE(na+nb+1)E(n_a+n_b+1)+E+Estable
scalar + FermiE(na+1−nb)E(n_a+1-n_b): {0,E,E,2E}\{0,E,E,2E\}−Ep-E_p per filled bb-modecollapses with mode count
Dirac + FermiE(na+nb)−EE(n_a+n_b)-E: {−E,0,0,E}\{-E,0,0,E\}−E-E constantstable (vacuum shift)
Dirac + BoseE(na−nb−1)E(n_a-n_b-1)−(nb+1)E→−∞-(n_b+1)E\to-\inftycollapses in one mode

Step 3. Boundary number. For the wrong-Dirac case the floor at bosonic occupation nb=20n_b=20 is already −21E\mathbf{-21E}; nothing stops nb=106n_b=10^6. For the wrong-scalar case, two modes with E1=1,E2=3E_1=1,E_2=3 sit 4E\mathbf{4E} below their own bb-vacuum — and the real field has a continuum of modes. Both wrong cells have no ground state; both right cells do.

Origins · 1 minSpin-Statistics Theorem and how this was discovered
ST
Spin-Statistics Theorem
1939–1958 · Zurich, Princeton

The theorem was assembled in the shadow of war. Markus Fierz — Pauli's assistant in Zurich — proved the free-field case in 1939 (Helv. Phys. Acta 12, 3): integer-spin fields demand commutators, half-integer anticommutators, on pain of losing causality or the energy floor. Pauli himself, marooned at Princeton's Institute for Advanced Study in 1940, distilled it into "The Connection Between Spin and Statistics" (Phys. Rev. 58, 716) — the paper that turned his own 1925 exclusion rule (told in Pauli Exclusion Principle) from decree into consequence. Lüders and Zumino, and independently Burgoyne, re-proved it in 1958 from field-theory axioms alone; Streater and Wightman canonized it in 1964. The coda is cautionary: Duck and Sudarshan's 1997 survey examined seventy years of "elementary" proofs — including Feynman's 1986 belt-trick argument — and found every nonrelativistic shortcut wanting.

History
1925
— Pauli's exclusion rule and, months later, spin itself enter atomic physics as separate postulates
1926
— Fermi and Dirac work out the statistics of antisymmetric ensembles; Dirac's paper splits the two sectors
1939
— Fierz proves the free-field spin-statistics connection in Helvetica Physica Acta
1940
— Pauli's Phys. Rev. 58, 716, written in Princeton exile, gives the theorem its definitive wartime form
1958
— Lüders–Zumino and Burgoyne re-derive it from axiomatic field theory; PCT and spin-statistics become siblings
1964
— Streater and Wightman canonize the axiomatic proof in PCT, Spin and Statistics, and All That
1975
— neutron interferometers watch a spinor change sign under 2π rotation — the theorem's key ingredient, measured
1997
— Duck and Sudarshan audit all "simple proofs" and conclude none succeeds without relativity

Going further

References. Pauli, "The Connection Between Spin and Statistics", Phys. Rev. 58, 716 (1940) · Fierz, Helv. Phys. Acta 12, 3 (1939) · Streater & Wightman, PCT, Spin and Statistics, and All That (1964) · Weinberg, The Quantum Theory of Fields Vol. I, Ch. 5 · Peskin & Schroeder, An Introduction to Quantum Field Theory, Sec. 3.5 · Duck & Sudarshan, Pauli and the Spin-Statistics Theorem (1997) · Griffiths, Introduction to Quantum Mechanics, Ch. 5 · Sakurai & Napolitano, Modern Quantum Mechanics, Ch. 7.

Provenance. Field expansions, (anti)commutators, and spin sums follow Peskin & Schroeder Sec. 3.5 and Weinberg QFT Vol. I Ch. 5. A verifier should confirm: the scalar commutator Δ(x−y)\Delta(x-y) vanishes for spacelike separation while the anticommutator does not; the Dirac spin sums ∑uuˉ=p+m\sum u\bar u=\cancel{p}+m, ∑vvˉ=p−m\sum v\bar v=\cancel{p}-m; and the energy-positivity argument requiring anticommutators for spin-1/2. Historical attribution to Fierz (1939) and Pauli (Phys. Rev. 58, 716, 1940) per Duck & Sudarshan (1997). Nonrelativistic consequences cross-checked against Griffiths Ch. 5 and Sakurai Ch. 7. (This lane's verify script executes the first and third confirmations exactly — the symbolic spacelike zero, the K₁ closed form to 3×10⁻⁹, and the four mode-level spectra; the spinor completeness sums are quoted with their source above.)

Builds on: Indistinguishability of Identical Particles · Permutation Symmetry and the Exchange Operator · Symmetric and Antisymmetric States · Bosons and Fermions. Leads to: Pauli Exclusion Principle · Slater Determinants · Aufbau Principle and Electron Configurations · The Helium Atom.

Is the spin-statistics connection an experimental fact or a theorem?

Both, at different layers. Within nonrelativistic quantum mechanics it is a postulate read off from experiment — helium spectroscopy, electron shells, photon bunching. Within relativistic quantum field theory it is a theorem: given Lorentz invariance, microcausality, an energy floor and positive norms, no other pairing is consistent. The experiments then double as tests of those axioms, which is why Pauli-violation searches are really tests of field theory's foundations.

Why does the anticommutator function Δ₁ being nonzero matter if it is exponentially small?

Because causality is an exact structural requirement, not an approximate one. The independence of spacelike measurements is what lets a relativistic theory assign consistent probabilities across reference frames; an acausal correlation of size e⁻¹³⁷ is as disqualifying as one of size 1 — it means the operator algebra itself connects regions no signal can connect. Small numbers can be negligible in predictions; they cannot be negligible in axioms.

Where exactly does the minus sign in the Dirac raw Hamiltonian come from?

From the completeness sums of the negative-frequency spinors: the antiparticle projector carries the opposite sign of the mass term relative to the particle projector. When the mode expansion is inserted into the energy density, that projector sign lands on the antiparticle term as −bb†. We quote this with its source (Peskin & Schroeder §3.5) rather than derive it, because spinor algebra belongs to the relativistic chapter of this track; the audit here needs only the sign.

Does the theorem say anything about composite particles like helium atoms?

Indirectly and decisively. The theorem fixes the statistics of the constituents (quarks, electrons: fermions); the exchange bookkeeping of Bosons and Fermions then makes a composite a boson or fermion according to the parity of its fermion count. That is why ⁴He condenses and ³He must pair first — a one-neutron difference flipping the sector, exactly as the composite-exchange rule demands.

What would it even look like experimentally if spin-statistics failed a little bit?

The standard parametrization imagines a small admixture of the wrong sector — a tiny probability that a new electron entering an atom binds symmetrically with an existing one. That would let it fall into a filled 1s shell and emit a shifted X-ray; VIP-2 at Gran Sasso looks for exactly those photons in current-carrying copper and sees none, bounding any violation below one part in 102910^{29}. Symmetric-sector photons would likewise alter two-photon selection rules; spectroscopy sees no such lines.

Do anyons violate the spin-statistics theorem?

No — they sidestep its hypotheses. The theorem's proof machinery assumes three or more spatial dimensions, where exchanging two particles twice is topologically trivial and the permutation group governs. In two dimensions the double exchange can wind, the braid group takes over, and a continuous exchange phase becomes consistent. Anyons in quantum Hall systems live exactly in that loophole; in 3D the theorem stands.

Before this

Where this leads