Slater Determinants

⚡ Key Takeaways

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.

💡 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.

DefinitionSpin-orbital

A spin-orbital $\chi_i(\mathbf{x})$ is a single-particle wavefunction depending on the combined coordinate $\mathbf{x} = (\mathbf{r}, \omega)$ , where $\mathbf{r}$ is the spatial position and $\omega$ is the spin variable. Typically $\chi_i(\mathbf{x}) = \psi_i(\mathbf{r})\,\sigma(\omega)$ , a spatial orbital times a spin function $\sigma \in \{\alpha, \beta\}$ (spin up/down). We assume the set is orthonormal:

1

\quad \langle \chi_i | \chi_j \rangle = \int \chi_i^*(\mathbf{x})\,\chi_j(\mathbf{x})\,d\mathbf{x} = \delta_{ij}

Derivation hint: $\int d\mathbf{x}$ denotes integration over space and summation over the discrete spin variable; orthonormality is a defining convention, choosable by Gram–Schmidt.

DefinitionSlater Determinant

For $N$ identical fermions occupying orthonormal spin-orbitals $\chi_1,\dots,\chi_N$ , the Slater determinant is

2

\quad \Psi(\mathbf{x}_1,\dots,\mathbf{x}_N) = \frac{1}{\sqrt{N!}}\begin{vmatrix} \chi_1(\mathbf{x}_1) & \chi_2(\mathbf{x}_1) & \cdots & \chi_N(\mathbf{x}_1) \\ \chi_1(\mathbf{x}_2) & \chi_2(\mathbf{x}_2) & \cdots & \chi_N(\mathbf{x}_2) \\ \vdots & \vdots & \ddots & \vdots \\ \chi_1(\mathbf{x}_N) & \chi_2(\mathbf{x}_N) & \cdots & \chi_N(\mathbf{x}_N) \end{vmatrix}

Derivation hint: row $j$ holds particle $\mathbf{x}_j$ ; column $i$ holds spin-orbital $\chi_i$ . The prefactor $1/\sqrt{N!}$ normalizes $\Psi$ when the $\chi_i$ are orthonormal.

TheoremAntisymmetry and Pauli Exclusion

The Slater determinant $\Psi$ satisfies: 1. Antisymmetry: $\Psi(\dots,\mathbf{x}_a,\dots,\mathbf{x}_b,\dots) = -\,\Psi(\dots,\mathbf{x}_b,\dots,\mathbf{x}_a,\dots)$ for any pair $a,b$ . 2. Pauli exclusion: if $\chi_i = \chi_j$ for $i\ne j$ , then $\Psi \equiv 0$ . 3. Normalization: $\langle \Psi | \Psi \rangle = 1$ when the $\chi_i$ are orthonormal.

ProofAntisymmetry, Pauli, and Normalization

Step 1 — Antisymmetry. Swapping particles $\mathbf{x}_a$ and $\mathbf{x}_b$ in equation [2] interchanges rows $a$ and $b$ of the determinant. A determinant changes sign under interchange of two rows. Hence $\Psi \to -\Psi$ , proving property (1).

Step 2 — Pauli exclusion (from equal orbitals). If $\chi_i = \chi_j$ , two columns of the determinant are identical. A determinant with two equal columns vanishes, so $\Psi \equiv 0$ . No two fermions can occupy the same spin-orbital.

Step 3 — Pauli exclusion (from equal coordinates). Set $\mathbf{x}_a = \mathbf{x}_b$ . Then rows $a$ and $b$ are equal, so $\Psi = 0$ . Two same-spin electrons cannot be at the same point — the "Fermi/exchange hole."

Step 4 — Normalization. Using the expanded form [3]:

5

\quad \langle \Psi | \Psi \rangle = \frac{1}{N!}\sum_{P,Q}\operatorname{sgn}(P)\operatorname{sgn}(Q)\prod_i \langle \chi_i(\mathbf{x}_{P(i)}) | \chi_i(\mathbf{x}_{Q(i)})\rangle

Derivation hint: take the bra–ket of two sums over permutations and integrate over all $\mathbf{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 $\langle\chi_i|\chi_k\rangle = 0$ ). There are $N!$ surviving terms, each equal to $1$ , with $\operatorname{sgn}(P)^2 = 1$ :

6

\quad \langle \Psi | \Psi \rangle = \frac{1}{N!}\sum_{P}1 = \frac{N!}{N!} = 1

Derivation hint: count the diagonal $P=Q$ permutations; there are exactly $|S_N| = N!$ of them. This proves all three properties. $\blacksquare$

∎

JC

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 $\alpha$ 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.

🧠 Intuition

Start with the naive guess: $N$ independent particles, each in its own state, give a product $\chi_1(\mathbf{x}_1)\chi_2(\mathbf{x}_2)\cdots$ — 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/\sqrt{N!}$ is bookkeeping: there are $N!$ terms, each normalized, so dividing by $\sqrt{N!}$ keeps total probability at one.

Figure

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.

ExampleTwo-electron determinant: helium ground state (model)

Given: Two electrons in the same spatial orbital $\psi_{1s}(\mathbf{r})$ but opposite spins, $\chi_1 = \psi_{1s}\alpha$ , $\chi_2 = \psi_{1s}\beta$ . Build the $2\times2$ Slater determinant.

Step 1. Write the determinant from [2] with $N=2$ :

7

\quad \Psi(\mathbf{x}_1,\mathbf{x}_2) = \frac{1}{\sqrt{2}}\begin{vmatrix} \chi_1(\mathbf{x}_1) & \chi_2(\mathbf{x}_1) \\ \chi_1(\mathbf{x}_2) & \chi_2(\mathbf{x}_2)\end{vmatrix}

Derivation hint: expand a $2\times2$ determinant as $ad - bc$ .

Step 2. Expand:

8

\quad \Psi = \frac{1}{\sqrt{2}}\big[\psi_{1s}(\mathbf{r}_1)\psi_{1s}(\mathbf{r}_2)\big]\big[\alpha(1)\beta(2) - \beta(1)\alpha(2)\big]

Derivation hint: factor out the common spatial part $\psi_{1s}(\mathbf{r}_1)\psi_{1s}(\mathbf{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:

\boxed{\Psi = \psi_{1s}(\mathbf{r}_1)\psi_{1s}(\mathbf{r}_2)\cdot\tfrac{1}{\sqrt2}\big[\alpha(1)\beta(2)-\beta(1)\alpha(2)\big]}

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.)

ExampleEnergy expectation: Slater–Condon for two orbitals

Given: A two-electron determinant of orthonormal spin-orbitals $\chi_a,\chi_b$ , with Hamiltonian $\hat{H} = \hat{h}(1)+\hat{h}(2) + 1/r_{12}$ (atomic units). Find $\langle\Psi|\hat H|\Psi\rangle$ .

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

9

\quad \langle\Psi|\hat h(1)+\hat h(2)|\Psi\rangle = h_{aa} + h_{bb}, \qquad h_{ii}=\langle\chi_i|\hat h|\chi_i\rangle

Derivation hint: 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$ :

10

\quad \Big\langle\Psi\Big|\tfrac{1}{r_{12}}\Big|\Psi\Big\rangle = J_{ab} - K_{ab}

where $J_{ab}=\langle\chi_a\chi_b|\tfrac{1}{r_{12}}|\chi_a\chi_b\rangle$ and $K_{ab}=\langle\chi_a\chi_b|\tfrac{1}{r_{12}}|\chi_b\chi_a\rangle$ . Derivation hint: the antisymmetric pairs in [3] generate direct ( $P=Q$ ) and swapped terms; the swapped one supplies the minus exchange integral.

Step 3. Total energy:

11

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

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

Boxed answer:

\boxed{E = h_{aa}+h_{bb}+J_{ab}-K_{ab}\ \text{(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 mistake

"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/\sqrt{N!}$ is arbitrary." It is fixed by requiring $\langle\Psi|\Psi\rangle = 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 $\chi_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).

CounterexampleAn antisymmetric state that is NOT a single determinant

Consider two electrons and four spin-orbitals. The state

\Phi = \tfrac{1}{\sqrt2}\big(|\chi_1\chi_2\rangle + |\chi_3\chi_4\rangle\big)

(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.

connection

The Slater determinant is the natural object in second quantization: a determinant of $\chi_{k_1},\dots,\chi_{k_N}$ corresponds to the occupation-number state $\hat a^\dagger_{k_1}\cdots\hat a^\dagger_{k_N}|0\rangle$ , 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$ .

🔬 Where this appears

⏳ History

- 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.

open-problem

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.

Formal Statement

Why It Holds / Derivation

Building Intuition

Worked Examples

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also