Creation and Annihilation Operators

There is a bookkeeping catastrophe hiding inside identical-particle quantum mechanics. Write two electrons as $\Psi(\mathbf r_1,\mathbf r_2)$ and you must antisymmetrize by hand; write ten and you are juggling $10!=3628800$ signed permutations in a Slater determinant; write a metal's worth and the notation collapses under its own weight. Worse, the number of particles is often not even fixed — photons are created and destroyed every time an atom emits, and a superconductor's ground state is a superposition of different electron counts. Labelled coordinates are the wrong language for a world of anonymous, appearing-and-disappearing particles. Creation and annihilation operators are the right language. Instead of asking which particle sits where, you ask only how many particles occupy each single-particle mode, and you build every state by acting with $a^\dagger$ on the empty vacuum $|0\rangle$ . The miracle is that a handful of algebraic relations — commutators for bosons, anticommutators for fermions — carry the entire content of Bose and Fermi statistics. Antisymmetry stops being a constraint you impose and becomes a sign the algebra tracks for you. This article develops that algebra from the harmonic-oscillator ladder you already know, proves the ladder normalizations $\sqrt{n+1}$ and $\sqrt{n}$ , and shows that the single relation $\{c,c^\dagger\}=1$ forces fermion occupation into $\{0,1\}$ — the Pauli principle written as the operator identity $(c^\dagger)^2=0$ . The neighbouring articles Second Quantization and Occupation-Number Representation and Fock Space set the stage and hand the algebra to us as its rightful owner; here we make it work.

Key takeaways

Creation operators $a^\dagger$ add a particle to a mode; annihilation operators $a$ remove one; the vacuum $|0\rangle$ is emptied by every $a$ .
Bosons obey commutators $[a_i,a_j^\dagger]=\delta_{ij}$ ; fermions obey anticommutators $\{c_i,c_j^\dagger\}=\delta_{ij}$ — that single choice is the boson/fermion distinction at the operator level.
The commutator forces the boson ladder factors $a^\dagger|n\rangle=\sqrt{n+1}\,|n+1\rangle$ , $a|n\rangle=\sqrt{n}\,|n-1\rangle$ , with no ceiling on $n$ .
The anticommutator forces $(c^\dagger)^2=0$ and $\hat n^2=\hat n$ , so fermion occupation is exactly $\{0,1\}$ — the Pauli exclusion principle as one line of algebra.
Cross-mode fermion anticommutation makes $c_i^\dagger c_j^\dagger|0\rangle=-c_j^\dagger c_i^\dagger|0\rangle$ , reproducing Slater-determinant antisymmetry automatically.
Key idea: two (anti)commutation relations among $a$ and $a^\dagger$ replace the entire apparatus of explicit (anti)symmetrization — the algebra remembers the statistics so the states don't have to be symmetrized by hand.

The core insight

A handful of (anti)commutation relations among $a$ and $a^\dagger$ replace the entire machinery of explicit (anti)symmetrization: choosing the commutator or the anticommutator fixes the statistics, and the algebra then carries the signs — and the Pauli principle — for you.

Where this lives today

The algebra of $a$ and $a^\dagger$ is arguably the most reused notation in modern physics. In quantum optics every mode of light is a bosonic oscillator, and photon-number states, coherent states $|\alpha\rangle=e^{-|\alpha|^2/2}\sum_n\alpha^n/\sqrt{n!}\,|n\rangle$ , and squeezed states are all defined through these operators; the fact that $a^\dagger|0\rangle\ne0$ is precisely why spontaneous emission happens even in the dark. In condensed matter the Hubbard, BCS, and tight-binding models are written entirely as $c_i^\dagger c_j$ hopping and $U\hat n_{i\uparrow}\hat n_{i\downarrow}$ interaction terms, with $(c^\dagger)^2=0$ enforcing single occupancy automatically. The most consequential recent use is in quantum simulation and quantum chemistry on quantum computers: the Jordan–Wigner transform that we used above to build two-mode fermion operators is exactly how molecular fermionic Hamiltonians are mapped onto qubit gates, so that a quantum processor can compute ground-state energies of molecules. The 1928 sign bookkeeping of Jordan and Wigner turns out to be the bridge between electrons and qubits.

The physical picture

Picture each single-particle mode — a momentum state, an atomic orbital, a photon polarization — as a numbered box. Two operators act on the boxes. The annihilation operator $a$ removes one particle from its box; the creation operator $a^\dagger$ drops one in. The state of the whole system is nothing more than the list of occupation numbers, $|n_1,n_2,\dots\rangle$ : how many particles sit in box $1$ , box $2$ , and so on. You never say which particle is in which box, because identical particles wear no labels — and that omission is exactly right, not a loss of information. The two statistics are two rules for stacking. For bosons a box holds any number of particles, and climbing it is literally the harmonic-oscillator ladder you met one mode at a time in Ladder Operators and the Algebraic Oscillator — one oscillator per mode, rungs at $n=0,1,2,\dots$ forever. For fermions a box holds at most one particle; asking $c^\dagger$ to add a second returns exactly zero. That hard stop is the Pauli principle, and it is not a separate postulate glued on — it drops straight out of the anticommutator, as we will prove. The signs are the memory of antisymmetry. Reordering two fermion creation operators flips the sign, $c_i^\dagger c_j^\dagger=-c_j^\dagger c_i^\dagger$ , exactly as swapping two columns flips a Slater determinant. Nothing is lost in translation; the algebra does the antisymmetrizing that Slater Determinants does by hand.

Figure

Boson ladder (no ceiling) versus fermion ladder (two rungs)

Assembling a many-body state is then just a sequence of $a^\dagger$ 's on the vacuum. Fock space is the direct sum of all particle-number sectors, and $a^\dagger$ walks you up the tower from one sector to the next — no relabelling of particles, ever.

Figure

Building Fock space by raising the vacuum

The core equations

The whole formalism grows from one relation per statistics. For a single bosonic mode the commutator $[a,a^\dagger]=1$ is identical to the harmonic-oscillator algebra, and it alone fixes how the operators move between Fock states. The number operator $\hat n=a^\dagger a$ has nonnegative-integer eigenvalues, and the raising and lowering factors are not free constants — they are forced by normalization, as the first derivation shows.

a^\dagger|n\rangle=\sqrt{n+1}\,|n+1\rangle,\qquad a|n\rangle=\sqrt{n}\,|n-1\rangle \qquad (1)

where: $|n\rangle$ is the normalized occupation-number (Fock) state of a single mode holding $n$ quanta; $a^\dagger,a$ raise and lower that count; the square-root factors are the unique normalization keeping every $|n\rangle$ of unit length. Read physically: the ground rung is the vacuum $|0\rangle$ with $a|0\rangle=0$ (you cannot remove what is not there), and there is no top rung — a bosonic mode can be filled without limit. The commutation relations that generate this ladder, and their fermionic partners, are the algebra's entire axiom set. Everything else is a theorem.

[a_i,a_j^\dagger]=\delta_{ij},\qquad [a_i,a_j]=[a_i^\dagger,a_j^\dagger]=0 \qquad (2)

where: $i,j$ label single-particle modes (momentum, orbital, spin-polarization); $\delta_{ij}$ picks out the same mode. Distinct bosonic modes commute freely, which is why bosonic Fock states are fully symmetric under relabelling. For fermions the same structural relations hold with anticommutators replacing commutators. This one substitution is the whole of Fermi–Dirac statistics.

\{c_i,c_j^\dagger\}=\delta_{ij},\qquad \{c_i,c_j\}=\{c_i^\dagger,c_j^\dagger\}=0 \qquad (3)

where: $\{A,B\}=AB+BA$ is the anticommutator; $c_i,c_i^\dagger$ are the fermionic mode operators. Setting $i=j$ in the second relation gives $\{c_i^\dagger,c_i^\dagger\}=2(c_i^\dagger)^2=0$ , and the consequences are immediate and physical.

(c_i^\dagger)^2=0,\qquad \hat n_i^2=\hat n_i\ \Rightarrow\ n_i\in\{0,1\} \qquad (4)

where: $\hat n_i=c_i^\dagger c_i$ is the fermionic occupation operator. The left identity says a mode cannot be doubly created; the right says its occupation solves $n(n-1)=0$ , i.e. $n\in\{0,1\}$ . This is the Pauli exclusion principle — the same content that Pauli Exclusion Principle states in three languages, here as one operator identity. With the moves in hand, any many-body state is built from the vacuum. For bosons repeated raising on a mode needs a $\sqrt{n_i!}$ to stay normalized:

|n_1,n_2,\dots\rangle=\prod_i\frac{(a_i^\dagger)^{n_i}}{\sqrt{n_i!}}\,|0\rangle \qquad (5)

where: $|0\rangle$ is the vacuum with $a_i|0\rangle=0$ for every $i$ ; each $n_i$ is any nonnegative integer. For fermions the product of distinct $c_k^\dagger$ 's in a fixed order builds an antisymmetric state — the operator form of a Slater determinant:

|k_1,\dots,k_N\rangle=c_{k_1}^\dagger c_{k_2}^\dagger\cdots c_{k_N}^\dagger|0\rangle,\qquad k_1<k_2<\cdots<k_N \qquad (6)

where: the increasing order fixes the overall sign, and $N$ is the particle number with each mode occupied at most once. The continuum limit of these operators — the field operators $\psi(\mathbf r)=\sum_k\phi_k(\mathbf r)\,a_k$ that create and destroy particles at a point — is the subject of Field Operators and Many-Body Wavefunctions, which owns their $\delta$ -function (anti)commutators; here we supply only the per-mode algebra those field relations are built from.

Derivation AThe ladder factors from the commutator

We want to prove that the numbers in equation (1) — the $\sqrt{n+1}$ and $\sqrt{n}$ — are not chosen but forced, given only $[a,a^\dagger]=1$ and the definition $\hat n=a^\dagger a$ . The strategy is to show first that $a^\dagger$ and $a$ shift the eigenvalue of $\hat n$ by $\pm1$ , so they map Fock states to neighbouring Fock states up to a constant, and then pin that constant by demanding unit normalization. Watch that the argument also proves the spectrum of $\hat n$ is exactly the nonnegative integers, with a bottom rung annihilated by $a$ — there is no separate assumption of "integer quanta".

Step 2. Compute the commutators of $\hat n$ with the ladder operators. From $[a,a^\dagger]=1$ , $[\hat n,a^\dagger]=a^\dagger[a,a^\dagger]+[a^\dagger,a^\dagger]a=a^\dagger$ and likewise $[\hat n,a]=-a$ . Why this matters: it means $\hat n(a^\dagger|n\rangle)=(a^\dagger\hat n+a^\dagger)|n\rangle=(n+1)a^\dagger|n\rangle$ , so $a^\dagger|n\rangle$ is an eigenstate of $\hat n$ with eigenvalue $n+1$ ; hence $a^\dagger|n\rangle=\lambda_n|n+1\rangle$ for some constant $\lambda_n$ , and symmetrically $a|n\rangle=\mu_n|n-1\rangle$ .

Step 3. Fix the constants by normalization. $\|a^\dagger|n\rangle\|^2=\langle n|a a^\dagger|n\rangle=\langle n|(a^\dagger a+1)|n\rangle=n+1$ , so $|\lambda_n|=\sqrt{n+1}$ . Similarly $\|a|n\rangle\|^2=\langle n|a^\dagger a|n\rangle=n$ , so $|\mu_n|=\sqrt{n}$ . Choosing the standard phases gives exactly equation (1).

Step 4. The bottom rung. Lowering repeatedly must terminate or it would produce negative-norm states, contradicting Step 1; termination happens when $a|0\rangle=0$ , i.e. $\|a|0\rangle\|^2=0=n_{\min}$ . So the spectrum is $\{0,1,2,\dots\}$ with no ceiling. Verified: numerically, a truncated $12\times12$ ladder reproduces $a^\dagger|n\rangle=\sqrt{n+1}\,|n+1\rangle$ and $a|n\rangle=\sqrt{n}\,|n-1\rangle$ for all interior $n$ , with $\hat n$ eigenvalues $0,\dots,11$ and $[\hat n,a^\dagger]=a^\dagger$ ; see the verifier.

Go deeper · Derivation B · 1 minPauli exclusion from the anticommutator

Derivation BPauli exclusion from the anticommutator

Here we prove that fermion occupation is restricted to $\{0,1\}$ using nothing but the anticommutation relation $\{c,c^\dagger\}=1$ together with $\{c,c\}=0$ . The route is to first show the creation operator is nilpotent, $(c^\dagger)^2=0$ — you literally cannot write down a doubly-occupied fermion mode — and then to show the number operator is a projector, $\hat n^2=\hat n$ , whose only eigenvalues are $0$ and $1$ . The subtle move is the middle rearrangement, where the anticommutator is used to slide a $c$ past a $c^\dagger$ ; keep track of the sign there and everything falls out.

Step 1. Nilpotency. From $\{c^\dagger,c^\dagger\}=2(c^\dagger)^2=0$ we get $(c^\dagger)^2=0$ directly, and identically $c^2=0$ . Creating (or annihilating) twice in one mode gives zero.

Step 2. The number operator is idempotent. Write $\hat n^2=c^\dagger c\,c^\dagger c$ . Use $\{c,c^\dagger\}=1$ in the form $c c^\dagger=1-c^\dagger c$ :

\hat n^2=c^\dagger(c c^\dagger)c=c^\dagger(1-c^\dagger c)c=c^\dagger c-(c^\dagger)^2 c^2=\hat n-0=\hat n.

Why this is the whole story: an operator satisfying $\hat n^2=\hat n$ is a projector, and projectors have eigenvalues satisfying $\lambda^2=\lambda$ , i.e. $\lambda\in\{0,1\}$ .

Step 3. Read off the spectrum. $\hat n(\hat n-1)=0$ means every occupation number solves $n(n-1)=0$ , so $n\in\{0,1\}$ : at most one fermion per single-particle mode. This is the Pauli exclusion principle, obtained without ever antisymmetrizing a wavefunction. Verified: in the $2\times2$ fermion representation $c=\begin{pmatrix}0&1\\0&0\end{pmatrix}$ , sympy confirms $\{c,c^\dagger\}=\mathbb 1$ , $(c^\dagger)^2=0$ , $\hat n=\mathrm{diag}(0,1)$ , $\hat n^2=\hat n$ , and $\mathrm{solve}(n^2=n)=\{0,1\}$ .

Go deeper · Derivation C · 1 minCross-mode antisymmetry and the Slater determinant

Derivation CCross-mode antisymmetry and the Slater determinant

The single-mode algebra says nothing about exchange; that lives in the relations between different modes. We want to show that two fermion creation operators on different modes anticommute, $c_i^\dagger c_j^\dagger=-c_j^\dagger c_i^\dagger$ , and that projecting the two-fermion state onto position space reproduces a normalized $2\times2$ Slater determinant. The point of interest is that antisymmetry under particle exchange is identical to the sign flip under operator reordering — one fact, two costumes — so the determinant of Slater Determinants emerges here for free.

Step 1. The exchange sign. Setting $i\ne j$ in $\{c_i^\dagger,c_j^\dagger\}=0$ gives $c_i^\dagger c_j^\dagger=-c_j^\dagger c_i^\dagger$ . Acting on the vacuum, $c_i^\dagger c_j^\dagger|0\rangle=-c_j^\dagger c_i^\dagger|0\rangle$ : relabelling the two occupied modes flips the state's sign.

Step 2. Position projection. Overlapping the state $c_a^\dagger c_b^\dagger|0\rangle$ with the two-particle position bra $\langle \mathbf r_1\mathbf r_2|$ yields

\langle\mathbf r_1\mathbf r_2|c_a^\dagger c_b^\dagger|0\rangle=\tfrac1{\sqrt2}\big[\phi_a(\mathbf r_1)\phi_b(\mathbf r_2)-\phi_a(\mathbf r_2)\phi_b(\mathbf r_1)\big].

Why the minus sign is inevitable: it is the same anticommutation sign from Step 1, now expressed on the coordinate wavefunction — the row/column-swap antisymmetry of a determinant.

Step 3. Recognize the determinant. The bracket is $\det\begin{pmatrix}\phi_a(\mathbf r_1)&\phi_b(\mathbf r_1)\\\phi_a(\mathbf r_2)&\phi_b(\mathbf r_2)\end{pmatrix}/\sqrt2$ , the normalized two-electron Slater determinant. Swapping $\mathbf r_1\leftrightarrow\mathbf r_2$ changes its sign; setting $a=b$ makes it vanish, recovering exclusion in coordinate language. Verified: sympy confirms the bracketed combination is antisymmetric under $\mathbf r_1\leftrightarrow\mathbf r_2$ ; the two-mode Jordan–Wigner construction $c_1=c\otimes\mathbb 1$ , $c_2=Z\otimes c$ gives $\{c_1,c_2\}=0$ and $c_1^\dagger c_2^\dagger|0\rangle=-c_2^\dagger c_1^\dagger|0\rangle$ exactly.

The second-quantized Hamiltonian assembles these operators into a form valid in every particle-number sector at once. One-body terms take the form $\sum_{ij}\langle i|T|j\rangle\,a_i^\dagger a_j$ and two-body interactions $\tfrac12\sum_{ijkl}\langle ij|V|kl\rangle\,a_i^\dagger a_j^\dagger a_l a_k$ , with the reversed inner ordering $a_l a_k$ chosen so that direct and exchange matrix elements come out with the correct sign for either statistics. The equivalence of the one-body string to the first-quantized operator is proved and owned by Second Quantization; we supply the algebra it runs on.

Limiting cases

The two relations specialize cleanly to every regime you will meet.

Case	Result	Meaning
Single bosonic mode	$[a,a^\dagger]=1,\ \hat n=a^\dagger a$	Exactly the harmonic-oscillator ladder; eigenvalues $n=0,1,2,\dots$ unbounded above.
Single fermionic mode	$\{c,c^\dagger\}=1,\ c^2=0$	A two-state system: $\hat n^2=\hat n$ gives $n\in\{0,1\}$ only — Pauli exclusion.
Many bosonic modes	$[a_i,a_j^\dagger]=\delta_{ij}$	Independent oscillators; Fock states symmetric, normalized by $\sqrt{n_i!}$ .
Many fermionic modes	$\{c_i,c_j^\dagger\}=\delta_{ij}$	Cross-mode anticommutation ⇒ ordered products are Slater determinants.
Classical / large occupation	$a\to\sqrt{n}\,e^{i\theta}$	Coherent states with $\langle\hat n\rangle=\|\alpha\|^2\gg1$ recover a classical field amplitude.

Worked examples

Example 1 — Boson ladder normalization (a real number). Take a single bosonic mode in the Fock state $|3\rangle$ . Compute $a^\dagger|3\rangle$ and $a|4\rangle$ explicitly, and confirm the normalization by hand.

Step 1. Apply equation (1): $a^\dagger|3\rangle=\sqrt{3+1}\,|4\rangle=\sqrt4\,|4\rangle=\mathbf{2\,|4\rangle}$ .

Step 3. Now lower from $|4\rangle$ : $a|4\rangle=\sqrt4\,|3\rangle=\mathbf{2\,|3\rangle}$ . Applying $a^\dagger$ then $a$ returns $a\,a^\dagger|3\rangle=2\cdot2\,|3\rangle=4\,|3\rangle=(\hat n+1)|3\rangle$ , as required.

Result: $a^\dagger|3\rangle=2|4\rangle$ and $a|4\rangle=2|3\rangle$ ; the factor $2=\sqrt4$ appears in both, and $\|a^\dagger|3\rangle\|^2=\mathbf{4}$ . There is no ceiling — the next step gives $a^\dagger|4\rangle=\sqrt5\,|5\rangle$ , and so on forever.

Example 2 — Two-fermion occupation and total number (exam level). Build the two-mode fermion Fock space and find the eigenvalues of the total number operator $\hat N=c_1^\dagger c_1+c_2^\dagger c_2$ ; then verify the exchange sign numerically.

Step 1. Enumerate the states. Each mode holds $0$ or $1$ fermion, so the Fock space is spanned by $|00\rangle,|10\rangle,|01\rangle,|11\rangle$ — dimension $2^2=\mathbf4$ . This is the fermion counting theorem in miniature: $M$ modes give a $2^M$ -dimensional Fock space.

Step 2. Apply $\hat N$ . The occupations are $N=0,1,1,2$ on the four states, so the eigenvalues of $\hat N$ are $\{0,1,1,2\}$ — the middle value is doubly degenerate (either mode singly occupied).

Step 3. The exchange sign. Using the Jordan–Wigner representation $c_1=c\otimes\mathbb 1$ , $c_2=Z\otimes c$ with $Z=\mathrm{diag}(1,-1)$ , one finds $c_1^\dagger c_2^\dagger|00\rangle=-\,c_2^\dagger c_1^\dagger|00\rangle$ , both of unit norm.

Result: the two-mode fermion Fock space has dimension $\mathbf4$ , total-number eigenvalues $\mathbf{\{0,1,1,2\}}$ , and the doubly-occupied state $c_1^\dagger c_2^\dagger|0\rangle$ changes sign under mode exchange — antisymmetry with no wavefunction ever written. Contrast the bosonic two-mode space, which is infinite-dimensional because each mode can hold any number of quanta.

Origins · 1 mincreation operator and how this was discovered

creation operator

formalism: 1927–1932 · The operators were born in field theory before they became the language of the many-body problem

The creation and annihilation operators arrived with quantum field theory rather than with atomic physics. In 1927 Dirac quantized the electromagnetic field by treating photons as excitations created and destroyed by $a^\dagger$ and $a$ — the paper usually taken to mark the birth of quantum electrodynamics. The same year Jordan and Klein wrote the boson field in second-quantized form; in 1928 Jordan and Wigner introduced the anticommutator $\{c,c^\dagger\}=1$ that makes the formalism describe electrons, encoding antisymmetry as an algebraic sign. Fock organized the resulting states in 1932 into the direct-sum "Fock space" that carries his name.

History

1925

— Pauli states the exclusion principle; later rendered as the operator identity (c†)² = 0

1927

— Dirac quantizes the electromagnetic field with a and a†, founding QED (Proc. R. Soc. A 114, 243)

1927

— Jordan and Klein write the second-quantized boson field (Z. Phys. 45, 751)

1928

— Jordan and Wigner introduce fermionic anticommutation {c, c†} = 1 (Z. Phys. 47, 631)

1932

— Fock assembles occupation-number states into Fock space (Z. Phys. 75, 622)

Going further

References. Sakurai & Napolitano, Modern Quantum Mechanics, 2nd ed., §7.5–7.6 · Cohen-Tannoudji, Diu & Laloë, Quantum Mechanics Vol. II, chapter on second quantization · Fetter & Walecka, Quantum Theory of Many-Particle Systems, Ch. 1 · Negele & Orland, Quantum Many-Particle Systems, Ch. 1 · Shankar, Principles of Quantum Mechanics, 2nd ed., identical-particles chapter.

Provenance. All algebra was recomputed independently: the boson ladder factors $\sqrt{n+1},\sqrt{n}$ from $[a,a^\dagger]=1$ on a truncated $12$ -dimensional ladder; the fermion identities $(c^\dagger)^2=0$ , $\hat n^2=\hat n$ , and $n\in\{0,1\}$ from $\{c,c^\dagger\}=1$ in the $2\times2$ representation; the two-mode Jordan–Wigner construction giving $\{c_1,c_2\}=0$ and the exchange sign; and the total-number spectrum $\{0,1,1,2\}$ . Historical dates and journal references were verified (Dirac 1927 Proc. R. Soc. A 114, 243; Jordan–Wigner 1928 Z. Phys. 47, 631; Fock 1932 Z. Phys. 75, 622). No base physics error was found; the export's algebra is standard and reproduced faithfully.

Prerequisites: Ladder Operators and the Algebraic Oscillator · Second Quantization · Occupation-Number Representation and Fock Space · Symmetric and Antisymmetric States · Bosons and Fermions. Leads to: Field Operators and Many-Body Wavefunctions · Slater Determinants · Pauli Exclusion Principle · Anyons and Fractional Statistics.

Are

a

and

a^\dagger

Hermitian, and if not, are they observable?

Neither is Hermitian — they are adjoints of each other, $a^\dagger=(a)^\dagger$ , not self-adjoint — so neither is directly an observable. What is observable is the number operator $\hat n=a^\dagger a$ , which is Hermitian with real eigenvalues $0,1,2,\dots$ (or $0,1$ for fermions). Physical measurables are always built from Hermitian combinations: the quadratures $x\propto a+a^\dagger$ and $p\propto i(a^\dagger-a)$ in quantum optics, or the density $\psi^\dagger\psi$ . The non-Hermitian operators are computational tools that move you between eigenstates of the Hermitian ones.

Why the

\sqrt{n!}

for bosons but no factorial for fermions?

Because a bosonic mode can be raised many times, and each raising on an already-occupied mode brings a $\sqrt{n+1}$ factor from equation (1); stacking $n$ of them gives $\sqrt{n!}$ , which you divide out to normalize. A fermionic mode is raised at most once — $(c^\dagger)^2=0$ — so there is never repeated raising on one mode and no factorial appears. What fermions have instead is the ordering sign from anticommuting distinct creation operators past one another, which bosons lack.

Is the "vacuum" the same as "nothing"?

Not quite. The vacuum $|0\rangle$ is the unique state annihilated by every $a_i$ (or $c_i$ ): $a_i|0\rangle=0$ . It has zero particles but is a perfectly good, normalized quantum state, not the absence of a state. In field theory it carries zero-point energy and supports fluctuations; the fact that $a^\dagger|0\rangle\ne0$ is why you can create particles out of it. "Empty of excitations" is the right reading, not "empty of physics".

How do these operators know about spin-statistics? Why not use commutators for everything?

The algebra by itself does not derive spin-statistics — it implements it. You are free to write commutators or anticommutators, but only one choice yields states with the correct symmetry: symmetric (Bose) states need commutators, antisymmetric (Fermi) states need anticommutators. The theorem that ties this to spin — integer spin must be bosonic, half-integer fermionic — is a result of relativistic quantum field theory (microcausality plus positive energy) and is owned by Spin-Statistics Theorem. Non-relativistically, the choice is an input matched to observed statistics.

What goes wrong if I forget the normal-ordering in the two-body term?

The two-body interaction is written $\tfrac12\sum\langle ij|V|kl\rangle\,a_i^\dagger a_j^\dagger a_l a_k$ with both creation operators to the left. If you instead sandwich $V$ without this ordering you generate spurious self-interaction terms (a particle interacting with itself) and double-count pairs. The reversed inner labels $a_l a_k$ (not $a_k a_l$ ) are what make the direct term $\langle ij|V|ij\rangle$ and the exchange term $\langle ij|V|ji\rangle$ come out with the physically correct relative sign — attractive-looking bunching for bosons, the exchange hole for fermions.

Does the Jordan–Wigner string mean fermions are secretly nonlocal?

The string $Z\otimes Z\otimes\cdots\otimes c$ that represents a fermion operator on mode $k$ looks nonlocal — it touches every earlier mode. But this is an artifact of forcing anticommuting objects into a tensor-product (qubit) space that is built from commuting factors; the physics is perfectly local. The string is bookkeeping for the exchange sign, and all physical observables (densities, currents, energies) are built from even numbers of fermion operators, on which the strings cancel. It becomes practically important in quantum computing, where the string length sets the gate cost of simulating fermions on qubits.

Can occupation numbers be non-integer, e.g. for anyons?

For ordinary bosons and fermions the number operator has strictly integer eigenvalues — that is forced by the (anti)commutator, as derived above. Anyons, which exist only in two dimensions, generalize the exchange phase (from $\pm1$ to a general $e^{i\theta}$ ) but their number operators still count whole quasiparticles; what becomes fractional is the statistics and, in the fractional quantum Hall effect, the electric charge of the quasiparticle. That story — braid groups replacing the permutation group — is the subject of Anyons and Fractional Statistics; here the occupation is always a nonnegative integer.