Slater Determinants

A Slater determinant is an antisymmetric wavefunction for a system of identical fermions, built as the determinant of a matrix of single-particle spin-orbitals. Key result: it automatically satisfies the Pauli exclusion principle — interchanging any two particles flips the sign, and two identical occupied states give a zero wavefunction. Where it applies: quantum chemistry, atomic and molecular physics, Hartree–Fock theory, and condensed-matter many-body theory. Who/when: introduced by John C. Slater in 1929.

A Slater determinant is the simplest many-fermion wavefunction that respects the antisymmetry required by the spin-statistics theorem. Given $N$ orthonormal single-particle states (spin-orbitals) $χ_{1}, \dots, χ_{N}$ and $N$ identical fermions, one arranges the values $χ_{i} (x_{j})$ into an $N \times N$ matrix and takes its determinant, normalized by $1/ N!$ . Because a determinant changes sign when two columns (or rows) are swapped, the resulting wavefunction is automatically antisymmetric under particle exchange — and vanishes identically if two spin-orbitals coincide, encoding the Pauli exclusion principle. This matters because it is the foundational building block of nearly all of electronic-structure theory; the Hartree–Fock approximation is precisely the variational best single Slater determinant. The surprising fact: even a single Slater determinant, despite being a product-like object, builds in genuine quantum correlation between same-spin electrons (the "exchange hole"), purely from antisymmetry and with no interaction assumed.

The Core Insight

Identical fermions are indistinguishable, and their joint wavefunction must be antisymmetric under exchange of any two of them. A determinant is, by definition, the unique (up to scale) totally antisymmetric function of its columns. So if you treat each particle as a column and each occupied state as a row, the determinant is the antisymmetrization machine — automatic, exact, and economical.

Formal Statement

Spin-orbital

A spin-orbital $χ_{i} (x)$ is a single-particle wavefunction depending on the combined coordinate $x = (r, ω)$ , where $r$ is the spatial position and $ω$ is the spin variable. Typically $χ_{i} (x) = ψ_{i} (r) σ (ω)$ , a spatial orbital times a spin function $σ \in {α, β}$ (spin up/down). We assume the set is orthonormal:

⟨ χ_{i} ∣ χ_{j} ⟩ = \int χ_{i}^{*} (x) χ_{j} (x) d x = δ_{ij}

⟳

\int d x

denotes integration over space and summation over the discrete spin variable; orthonormality is a defining convention, choosable by Gram–Schmidt.

Slater Determinant

For $N$ identical fermions occupying orthonormal spin-orbitals $χ_{1}, \dots, χ_{N}$ , the Slater determinant is

Ψ (x_{1}, \dots, x_{N}) = \frac{1}{N !} χ_{1} (x_{1}) χ_{1} (x_{2}) ⋮ χ_{1} (x_{N}) χ_{2} (x_{1}) χ_{2} (x_{2}) ⋮ χ_{2} (x_{N}) \dots \dots ⋱ \dots χ_{N} (x_{1}) χ_{N} (x_{2}) ⋮ χ_{N} (x_{N})

⟳row

j

holds particle

x_{j}

; column

i

holds spin-orbital

χ_{i}

. The prefactor

1/ N!

normalizes

Ψ

when the

χ_{i}

are orthonormal.

Expanding the determinant via the permutation (Leibniz) formula gives the equivalent and often more useful expression:

Ψ (x_{1}, \dots, x_{N}) = \frac{1}{N !} P \in S_{N} \sum sgn (P) i = 1 \prod N χ_{i} (x_{P (i)})

⟳the Leibniz formula for a determinant sums over all permutations

P

of the symmetric group

S_{N}

, with each term weighted by the sign (parity) of

P

A compact operator form uses the antisymmetrizer $\hat{A}$ :

Ψ = N! \hat{A} [χ_{1} (x_{1}) χ_{2} (x_{2}) \dots χ_{N} (x_{N})], \hat{A} = \frac{1}{N !} P \in S_{N} \sum sgn (P) \hat{P}

⟳

\hat{A}

is a projector (

\hat{A}^{2} = \hat{A}

\hat{A}^{†} = \hat{A}

); applying it to the Hartree product and rescaling by

N!

reproduces equation [3].

Antisymmetry and Pauli Exclusion

The Slater determinant $Ψ$ satisfies:

Antisymmetry: $Ψ (\dots, x_{a}, \dots, x_{b}, \dots) = - Ψ (\dots, x_{b}, \dots, x_{a}, \dots)$ for any pair $a, b$ .
Pauli exclusion: if $χ_{i} = χ_{j}$ for $i \neq = j$ , then $Ψ \equiv 0$ .
Normalization: $⟨ Ψ∣Ψ ⟩ = 1$ when the $χ_{i}$ are orthonormal.

Why It Holds / Derivation

Antisymmetry, Pauli, and Normalization

Step 1 — Antisymmetry. Swapping particles

x_{a}

and

x_{b}

in equation [2] interchanges rows

a

and

b

of the determinant. A determinant changes sign under interchange of two rows. Hence

Ψ \to - Ψ

, proving property (1). Step 2 — Pauli exclusion (from equal orbitals). If

χ_{i} = χ_{j}

, two columns of the determinant are identical. A determinant with two equal columns vanishes, so

Ψ \equiv 0

. No two fermions can occupy the same spin-orbital. Step 3 — Pauli exclusion (from equal coordinates). Set

x_{a} = x_{b}

. Then rows

a

and

b

are equal, so

Ψ = 0

. Two same-spin electrons cannot be at the same point — the "Fermi/exchange hole." Step 4 — Normalization. Using the expanded form [3]:

⟨ Ψ∣Ψ ⟩ = \frac{1}{N !} P, Q \sum sgn (P) sgn (Q) i \prod ⟨ χ_{i} (x_{P (i)}) ∣ χ_{i} (x_{Q (i)})⟩

⟳take the bra–ket of two sums over permutations and integrate over all

x_{j}

Each integral factorizes; orthonormality [1] forces each one-particle overlap to a Kronecker delta, so only terms with $P = Q$ survive (others contain at least one $⟨ χ_{i} ∣ χ_{k} ⟩ = 0$ ). There are $N!$ surviving terms, each equal to $1$ , with $sgn (P)^{2} = 1$ :

⟨ Ψ∣Ψ ⟩ = \frac{1}{N !} P \sum 1 = \frac{N !}{N !} = 1

⟳count the diagonal

P = Q

permutations; there are exactly

∣ S_{N} ∣ = N!

of them.

This proves all three properties. $■$

John Clarke Slater | 1900–1976

An American physicist and a towering figure in early quantum theory of atoms and solids. In a 1929 Physical Review paper, "The Theory of Complex Spectra," Slater introduced the determinantal wavefunction to handle multi-electron atoms, replacing Heisenberg's and Dirac's cumbersome explicit antisymmetrization. The construction unified spin and the Pauli principle in one elegant object. Slater went on to develop the X $α$ method, the augmented-plane-wave method in solid-state physics, and trained generations at MIT. The "Slater–Condon rules" for matrix elements between determinants, and Slater-type orbitals (STOs), also bear his name.

Building Intuition

Start with the naive guess: $N$ independent particles, each in its own state, give a product $χ_{1} (x_{1}) χ_{2} (x_{2}) \dots$ — the Hartree product. But this labels the particles ("particle 1 is in state 1"), which is forbidden: identical fermions are indistinguishable. To fix it, we must add all relabelings. Bosons add them all with $+$ signs (symmetric). Fermions add them with alternating signs determined by the parity of the swap (antisymmetric). The alternating-sign sum over all permutations is exactly the algebraic definition of a determinant. So the determinant is not a clever trick — it is the only object that is built from single-particle states and is totally antisymmetric. The $1/ N!$ is bookkeeping: there are $N!$ terms, each normalized, so dividing by $N!$ keeps total probability at one.

The exchange hole: probability density of finding a second same-spin electron near a reference electron at the origin. The dip to zero at contact comes purely from antisymmetry.

Worked Examples

Two-electron determinant: helium ground state (model)

Given: Two electrons in the same spatial orbital

ψ_{1 s} (r)

but opposite spins,

χ_{1} = ψ_{1 s} α

χ_{2} = ψ_{1 s} β

. Build the

2 \times 2

Slater determinant. Step 1. Write the determinant from [2] with

N = 2

Ψ (x_{1}, x_{2}) = \frac{1}{2} χ_{1} (x_{1}) χ_{1} (x_{2}) χ_{2} (x_{1}) χ_{2} (x_{2})

⟳expand a

2 \times 2

determinant as

a d - b c

Step 2. Expand:

Ψ = \frac{1}{2} [ψ_{1 s} (r_{1}) ψ_{1 s} (r_{2})] [α (1) β (2) - β (1) α (2)]

⟳factor out the common spatial part

ψ_{1 s} (r_{1}) ψ_{1 s} (r_{2})

since both orbitals share it.

Step 3. The spatial part is symmetric; the spin part is the antisymmetric singlet. Their product is overall antisymmetric, as required. Boxed answer:

Ψ = ψ_{1 s} (r_{1}) ψ_{1 s} (r_{2}) \cdot \frac{1}{2} [α (1) β (2) - β (1) α (2)]

Two electrons can share a spatial orbital precisely because their spins differ — that is the closed-shell structure of helium. (Dimensionless wavefunction; spatial part carries units of length

^{- 3}

each.)

Energy expectation: Slater–Condon for two orbitals

Given: A two-electron determinant of orthonormal spin-orbitals

χ_{a}, χ_{b}

, with Hamiltonian

\hat{H} = \hat{h} (1) + \hat{h} (2) + 1/ r_{12}

(atomic units). Find

⟨ Ψ∣ \hat{H} ∣Ψ ⟩

. Step 1. One-electron part. The Slater–Condon rules give the sum of orbital energies:

⟨ Ψ∣ \hat{h} (1) + \hat{h} (2) ∣Ψ ⟩ = h_{aa} + h_{bb}, h_{ii} = ⟨ χ_{i} ∣ \hat{h} ∣ χ_{i} ⟩

⟳cross terms vanish by orthonormality; each electron contributes its own one-body integral.

Step 2. Two-electron part gives a Coulomb

J

minus exchange

K

⟨ Ψ \frac{1}{r _{12}} Ψ ⟩ = J_{ab} - K_{ab}

where

J_{ab} = ⟨ χ_{a} χ_{b} ∣ \frac{1}{r _{12}} ∣ χ_{a} χ_{b} ⟩

and

K_{ab} = ⟨ χ_{a} χ_{b} ∣ \frac{1}{r _{12}} ∣ χ_{b} χ_{a} ⟩

⟳the antisymmetric pairs in [3] generate direct (

P = Q

) and swapped terms; the swapped one supplies the minus exchange integral.

Step 3. Total energy:

E = h_{aa} + h_{bb} + J_{ab} - K_{ab}

⟳add [9] and [10]; this is the Hartree–Fock energy for a two-electron determinant.

Boxed answer:

E = h_{aa} + h_{bb} + J_{ab} - K_{ab} (hartree)

The exchange term

K_{ab}

— which is purely a consequence of the determinant's antisymmetry — lowers the energy whenever the spins are parallel (Hund's first rule).

Common Misconceptions

"A Slater determinant has no electron correlation." False in general. It contains exchange (Fermi) correlation between same-spin electrons exactly. What it lacks is Coulomb correlation between opposite-spin electrons; that omission is "the correlation energy" in the Hartree–Fock context. Saying "uncorrelated" is sloppy shorthand.
"The $1/ N!$ is arbitrary." It is fixed by requiring $⟨ Ψ∣Ψ ⟩ = 1$ given orthonormal orbitals; with non-orthogonal orbitals the normalization changes.
"Each electron is in a definite orbital." No — the determinant is a fully entangled object; you cannot say "particle 1 is in $χ_{1}$ ." Only the set of occupied orbitals is meaningful.
"Any antisymmetric wavefunction is a single Slater determinant." No. A general antisymmetric state is a linear combination of many determinants (configuration interaction).

An antisymmetric state that is NOT a single determinant

Consider two electrons and four spin-orbitals. The state

Φ = \frac{1}{2} (∣ χ_{1} χ_{2} ⟩ + ∣ χ_{3} χ_{4} ⟩)

(a sum of two determinants) is perfectly antisymmetric and normalized, yet it cannot be written as any single Slater determinant of two spin-orbitals. This is the heart of why exact wavefunctions require multi-determinant expansions (CI, CASSCF), and why the single-determinant Hartree–Fock approximation is only a starting point.

Connections

The Slater determinant is the natural object in second quantization: a determinant of $χ_{k_{1}}, \dots, χ_{k_{N}}$ corresponds to the occupation-number state $\overset{a}{^}_{k_{1}}^{†} \dots \overset{a}{^}_{k_{N}}^{†} ∣0 ⟩$ , where the fermionic creation operators automatically build in antisymmetry through their anticommutation relations Second Quantization. It is the trial wavefunction of Hartree–Fock Theory, whose variational minimization yields the canonical orbitals. It is the reference state for post-HF correlation methods like Configuration Interaction and coupled cluster. Its emergence is mandated by the Pauli Exclusion Principle and, more deeply, the Spin-Statistics Theorem. The mathematics it relies on is the Determinant and the antisymmetrizing properties of the symmetric group $S_{N}$ .

Real-World Applications

Quantum chemistry (Hartree–Fock / SCF): every molecular electronic-structure calculation starts from one or more Slater determinants; programs like Gaussian, ORCA, and PySCF build them from STOs or Gaussian basis sets.
Atomic spectroscopy: Slater's original use — computing term energies and multiplet splittings of complex atoms via Slater–Condon parameters $F^{k}$ , $G^{k}$ .
Density Functional Theory: the Kohn–Sham scheme uses a single Slater determinant of fictitious non-interacting orbitals to reproduce the true density.
Solid-state band theory: Bloch-orbital Slater determinants describe metals and semiconductors; the exchange hole underlies the cohesive energy of metals.
Quantum Monte Carlo: Slater–Jastrow trial wavefunctions multiply a determinant by a correlation factor; the determinant supplies the nodal structure (the "fixed-node" surface).
Nuclear physics: the shell model and Hartree–Fock–Bogoliubov theory build nucleon wavefunctions as Slater determinants of single-particle nuclear states.
Quantum dots and cold atoms: few-fermion systems in traps are modeled by small Slater determinants, directly probing exchange effects experimentally.

Historical Development

1925 — Pauli formulates the exclusion principle for electron quantum numbers.
1926 — Fermi & Dirac develop the statistics of indistinguishable antisymmetric particles (Fermi–Dirac statistics).
1926 — Heisenberg & Dirac show many-electron wavefunctions must be antisymmetric, using explicit but unwieldy combinations.
1927 — Hartree introduces the self-consistent field using simple product (Hartree) wavefunctions — antisymmetry not yet built in.
1929 — John C. Slater publishes "The Theory of Complex Spectra" (Phys. Rev. 34, 1293), introducing the determinantal wavefunction that elegantly enforces antisymmetry and incorporates spin.
1930 — Fock generalizes Hartree's method to antisymmetric determinantal wavefunctions, giving the Hartree–Fock equations.
1931–1932 — Slater & Condon derive the rules ("Slater–Condon rules") for evaluating matrix elements between determinants.
1950s — Roothaan & Hall recast Hartree–Fock in a finite basis (the Roothaan equations), making determinantal methods computationally practical.
1965 — Kohn & Sham repurpose a single determinant of non-interacting orbitals as the cornerstone of practical DFT.

Subtleties & Limits

A single Slater determinant captures exchange but not the full Coulomb correlation, and recovering "the correlation energy" — formally the difference between the exact and Hartree–Fock energies — remains the central computational challenge of electronic-structure theory. The exact many-fermion wavefunction generally requires exponentially many determinants, and no polynomial-time algorithm is known for the general interacting-fermion ground-state problem (it is QMA-hard in the worst case). Open practical issues include: the fixed-node problem in diffusion Monte Carlo, where the unknown exact nodal surface limits accuracy; strong (static) correlation in bond-breaking and transition metals, where a single determinant qualitatively fails and multireference methods are mandatory; and the design of density functionals that mimic correlation cheaply yet reliably. Whether a compact, systematically improvable representation of correlated fermionic states exists in general is still unresolved.