Chapter 16: Multi-Electron Atoms and the Building of the Periodic Table

"The underlying physical laws necessary for the mathematical theory of a large part of physics and the whole of chemistry are thus completely known, and the difficulty is only that the exact application of these laws leads to equations much too complicated to be soluble." --- Paul Dirac, 1929

"The periodic table is quantum mechanics' gift to chemistry --- and chemistry's greatest validation of quantum mechanics."

In Chapter 5, we solved the hydrogen atom exactly. One electron, one proton, one Coulomb potential --- and we obtained the energy levels $E_n = -13.6\;\text{eV}/n^2$, the wavefunctions $\psi_{nlm}(r,\theta,\phi) = R_{nl}(r)Y_l^m(\theta,\phi)$, and the quantum numbers $n$, $l$, $m$ that organize the entire structure. That solution was clean, exact, and beautiful.

Now we confront reality: every atom except hydrogen has more than one electron. Helium has two. Carbon has six. Iron has twenty-six. Uranium has ninety-two. And the moment we add a second electron, the problem becomes exactly unsolvable. The electron--electron repulsion term $e^2/4\pi\epsilon_0 r_{12}$ couples the two electrons together, destroying the separability that made hydrogen tractable.

Yet here is the remarkable fact: despite this mathematical intractability, quantum mechanics explains the entire periodic table. Every element's chemical behavior, every trend in ionization energy, atomic radius, electron affinity, and electronegativity --- all of it follows from the Schrodinger equation, the Pauli exclusion principle, and a handful of clever approximation schemes. This chapter develops those schemes and shows how the periodic table emerges as quantum mechanics' single greatest triumph in chemistry.

We will build on three prerequisite chapters: - Chapter 5 (hydrogen atom): the foundation --- quantum numbers, radial wavefunctions, the $1/r$ potential - Chapter 14 (addition of angular momentum): coupling $\hat{\mathbf{L}}$ and $\hat{\mathbf{S}}$, Clebsch--Gordan coefficients, term symbols - Chapter 15 (identical particles): the symmetrization postulate, the Pauli exclusion principle, Slater determinants

🏃 Fast Track: If your primary interest is in how quantum mechanics explains the periodic table, you may skim Sections 16.2 and 16.6 (central field approximation and Hartree--Fock) and focus on Sections 16.3--16.5 and 16.7--16.8, where the physical consequences are richest. However, the approximation sections contain essential conceptual machinery that reappears in Chapters 17--19.

16.1 Beyond Hydrogen: The Multi-Electron Challenge

The Hamiltonian We Cannot Solve

Consider an atom with nuclear charge $Z$ and $N$ electrons. In Dirac notation, the full Hamiltonian is:

$$\hat{H} = \sum_{i=1}^{N}\left[-\frac{\hbar^2}{2m_e}\nabla_i^2 - \frac{Ze^2}{4\pi\epsilon_0 r_i}\right] + \sum_{i

The first sum contains the kinetic energy and the nuclear attraction for each electron --- these are the one-body terms. If this were all there was, the problem would separate into $N$ independent hydrogen-like equations, each with nuclear charge $Z$. We could solve them immediately.

The second sum is the problem. The electron--electron repulsion $e^2/(4\pi\epsilon_0 r_{ij})$, where $r_{ij} = |\mathbf{r}_i - \mathbf{r}_j|$, couples every pair of electrons. For helium ($N = 2$), there is one such term. For carbon ($N = 6$), there are $\binom{6}{2} = 15$ terms. For uranium ($N = 92$), there are $\binom{92}{2} = 4186$ terms. These terms make exact analytical solutions impossible for any $N > 1$.

💡 Key Insight: The difficulty is not merely computational --- it is fundamental. The electron--electron repulsion term depends on the distance between electrons, $r_{ij} = |\mathbf{r}_i - \mathbf{r}_j|$. This means that the probability of finding electron 1 at position $\mathbf{r}_1$ depends on where electron 2 happens to be, which depends on where electron 1 is, which depends on where electron 2 is... The electrons are correlated. Separating the Schrodinger equation into independent single-particle equations requires approximation.

How Big Is the Electron--Electron Interaction?

To appreciate the challenge, consider helium ($Z = 2$, $N = 2$). If we ignore electron--electron repulsion entirely, each electron sees a hydrogen-like potential with $Z = 2$:

$$E_{\text{no repulsion}} = 2 \times \left(-\frac{Z^2 \times 13.6\;\text{eV}}{n^2}\right)\bigg|_{n=1,\,Z=2} = 2 \times (-54.4\;\text{eV}) = -108.8\;\text{eV}$$

The experimental ground state energy of helium is $-79.0\;\text{eV}$. The electron--electron repulsion contributes $+29.8\;\text{eV}$ --- roughly 27% of the total energy. This is not a small correction. It is comparable to the nuclear attraction itself. Any viable theory must account for it.

📊 By the Numbers: The relative importance of electron--electron repulsion: | Atom | $Z$ | Nuclear attraction (eV) | e-e repulsion (eV) | Ratio | |------|------|------------------------|---------------------|-------| | He | 2 | $-108.8$ | $+29.8$ | 27% | | Li | 3 | $-244.8$ | $+66.4$ | 27% | | C | 6 | $-979$ | $+250$ | 26% | | Fe | 26 | $-18{,}400$ | $+4{,}100$ | 22% |

The ratio decreases for heavy atoms (the nuclear attraction grows as $Z^2$ while the repulsion grows more slowly), but it never becomes negligible.

🔄 Check Your Understanding (Spaced Review --- Ch 5): The hydrogen atom energy formula is $E_n = -13.6\;\text{eV}/n^2$. For a hydrogen-like ion with nuclear charge $Z$ and one electron, what is the energy formula? What is the ground state energy of He$^+$ ($Z = 2$)? Of Li$^{2+}$ ($Z = 3$)? Why can we solve these exactly but not neutral helium?

The Helium Atom: Simplest Multi-Electron System

Before attacking the general case, let us examine helium ($Z = 2$, $N = 2$) more carefully. The Hamiltonian is:

$$\hat{H}_{\text{He}} = \underbrace{-\frac{\hbar^2}{2m_e}\nabla_1^2 - \frac{2e^2}{4\pi\epsilon_0 r_1}}_{\hat{h}_1\;\text{(hydrogen-like, }Z=2\text{)}} + \underbrace{-\frac{\hbar^2}{2m_e}\nabla_2^2 - \frac{2e^2}{4\pi\epsilon_0 r_2}}_{\hat{h}_2} + \underbrace{\frac{e^2}{4\pi\epsilon_0 |\mathbf{r}_1 - \mathbf{r}_2|}}_{\text{electron-electron repulsion}}$$

The first two terms are independent hydrogen-like Hamiltonians --- we know their eigenstates and eigenvalues exactly. The third term, the electron--electron repulsion, couples the two electrons. This coupling term has no preferred direction in space (it depends only on the distance $|\mathbf{r}_1 - \mathbf{r}_2|$), but it depends on the positions of both electrons simultaneously, preventing separation.

Three strategies exist for dealing with this coupling, and we will encounter all three in this textbook:

1. Perturbation theory (Chapter 17): Treat the electron--electron repulsion as a perturbation on the two independent hydrogen-like Hamiltonians. This gives a first-order energy correction of $+34.0\;\text{eV}$, yielding $E \approx -74.8\;\text{eV}$ (compare experimental $-79.0\;\text{eV}$). The error is substantial because the perturbation is not small.

2. Variational method (Chapter 19): Use a trial wavefunction with an adjustable effective nuclear charge $Z'$. Minimizing the energy gives $Z' = Z - 5/16 = 1.6875$, meaning each electron partially screens the nuclear charge from the other. This yields $E \approx -77.5\;\text{eV}$ --- much better than perturbation theory.

3. Hartree--Fock method (Section 16.6): Solve self-consistently for the optimal single-particle orbitals. This gives $E \approx -77.9\;\text{eV}$, capturing all but the correlation energy of $\sim 1.1\;\text{eV}$.

🧪 Experiment: The most precise measurement of helium's ionization energy gives $E_I = 24.587\,387\,7(2)\;\text{eV}$ (Kandula et al., 2010). Combined with the exact He$^+$ energy ($-54.4\;\text{eV}$), this gives a total ground state energy of $-79.005\;\text{eV}$, known to better than one part in a billion. Matching this precision requires relativistic corrections, quantum electrodynamic (QED) effects, and nuclear recoil corrections well beyond Hartree--Fock. The helium atom is one of the most precisely calculated and measured systems in all of physics.

The Pauli Constraint

From Chapter 15, we know that electrons are fermions: the total wavefunction (spatial $\times$ spin) must be antisymmetric under the exchange of any two electrons. For $N$ electrons, the wavefunction can be written as a Slater determinant:

$$\Psi(\mathbf{x}_1, \mathbf{x}_2, \ldots, \mathbf{x}_N) = \frac{1}{\sqrt{N!}}\begin{vmatrix} \chi_a(\mathbf{x}_1) & \chi_b(\mathbf{x}_1) & \cdots & \chi_N(\mathbf{x}_1) \\ \chi_a(\mathbf{x}_2) & \chi_b(\mathbf{x}_2) & \cdots & \chi_N(\mathbf{x}_2) \\ \vdots & & \ddots & \vdots \\ \chi_a(\mathbf{x}_N) & \chi_b(\mathbf{x}_N) & \cdots & \chi_N(\mathbf{x}_N) \end{vmatrix}$$

where each $\chi_i(\mathbf{x}) = \psi_i(\mathbf{r})\sigma_i(s)$ is a spin-orbital --- the product of a spatial orbital and a spin function ($|\!\uparrow\rangle$ or $|\!\downarrow\rangle$), and $\mathbf{x} = (\mathbf{r}, s)$ denotes both spatial and spin coordinates. The determinantal form automatically guarantees antisymmetry: if two spin-orbitals are identical ($\chi_a = \chi_b$), two columns are equal and the determinant vanishes. This is the Pauli exclusion principle in action.

🔗 Connection (Ch 15): In Chapter 15 we derived the Slater determinant for two particles. Here we need it for $N$ particles. The mathematical structure is identical --- only the bookkeeping is more elaborate. Review Section 15.4 if the determinantal form is unfamiliar.

16.2 The Central Field Approximation

The Key Physical Idea

The central field approximation is the single most important conceptual move in multi-electron atomic physics. The idea is this: each electron moves not in the bare nuclear potential $-Ze^2/(4\pi\epsilon_0 r)$ alone, but in an effective central potential $V_{\text{eff}}(r)$ that approximately accounts for the average repulsion from all the other electrons.

If this effective potential is spherically symmetric --- depending only on $r$, not on $\theta$ or $\phi$ --- then the angular part of the single-particle Schrodinger equation separates exactly, just as it did for hydrogen. The angular solutions are still spherical harmonics $Y_l^m(\theta,\phi)$. The quantum numbers $l$ and $m_l$ are still good quantum numbers. Only the radial equation changes, because the radial potential is no longer the simple $-Ze^2/(4\pi\epsilon_0 r)$ Coulomb potential.

Mathematically, we replace the full $N$-electron Hamiltonian with $N$ independent single-particle Hamiltonians:

$$\hat{H} \approx \sum_{i=1}^{N}\hat{h}_i, \qquad \hat{h}_i = -\frac{\hbar^2}{2m_e}\nabla_i^2 + V_{\text{eff}}(r_i)$$

Each single-particle equation is then:

$$\hat{h}_i\,\psi_{n_i l_i m_i}(\mathbf{r}_i) = \varepsilon_{n_i l_i}\,\psi_{n_i l_i m_i}(\mathbf{r}_i)$$

where $\varepsilon_{n_i l_i}$ is the single-particle orbital energy.

The Crucial Consequence: Lifting the $l$-Degeneracy

Here is the key point --- the one that makes the periodic table structurally richer than hydrogen would suggest.

In hydrogen, the energy depends only on $n$: $E_n = -13.6\;\text{eV}/n^2$. The states $2s$ and $2p$ have the same energy. The states $3s$, $3p$, and $3d$ have the same energy. This "$l$-degeneracy" is a special property of the $1/r$ Coulomb potential (it reflects the hidden $SO(4)$ symmetry we noted in Chapter 5).

In a multi-electron atom, $V_{\text{eff}}(r)$ is not a pure $1/r$ potential. It looks approximately Coulombic at very small $r$ (where the electron is close to the nucleus and sees nearly the full charge $Z$) and at very large $r$ (where the electron sees the nucleus plus all inner electrons as a net charge of $+1$). But in between, the potential is modified by the screening effect of other electrons. This modified potential breaks the $l$-degeneracy.

💡 Key Insight: In a multi-electron atom, the energy of an orbital depends on both $n$ and $l$:

$$\varepsilon_{nl} \quad\text{not just}\quad \varepsilon_n$$

Specifically, for a given $n$, states with lower $l$ penetrate closer to the nucleus (their radial wavefunctions have more probability density near $r = 0$), experience less screening, and therefore have lower (more negative) energy. The ordering is:

$$\varepsilon_{ns} < \varepsilon_{np} < \varepsilon_{nd} < \varepsilon_{nf}$$

This is why the $2s$ orbital fills before $2p$, and why $3d$ fills after $4s$. Without this $l$-dependence, chemistry would be radically different.

Physical Picture: Screening and Penetration

Consider an electron at distance $r$ from the nucleus. The electrons closer to the nucleus (those with $r' < r$) partially cancel the nuclear charge, while those farther away (with $r' > r$) create no net force (by the shell theorem). The effective nuclear charge seen by an electron at $r$ is:

$$Z_{\text{eff}}(r) = Z - \sigma(r)$$

where $\sigma(r)$ is the screening function --- the average number of electrons inside radius $r$. The effective potential is then:

$$V_{\text{eff}}(r) \approx -\frac{Z_{\text{eff}}(r)\,e^2}{4\pi\epsilon_0\,r}$$

An $s$-orbital ($l = 0$) has nonzero probability density at $r = 0$: the radial wavefunction $R_{ns}(0) \neq 0$. A $p$-orbital ($l = 1$) vanishes at the origin as $r^1$, a $d$-orbital as $r^2$, and an $f$-orbital as $r^3$. Because $s$-orbitals penetrate more deeply toward the nucleus, they experience a larger effective $Z_{\text{eff}}$ on average and are therefore more tightly bound.

⚠️ Common Misconception: "Screening means outer electrons don't feel the nucleus." Wrong. Every electron feels the nucleus --- just with reduced effective charge. Even the outermost valence electron in cesium ($Z = 55$) still feels a net attraction; it merely sees an effective charge of roughly $Z_{\text{eff}} \approx 1$ rather than 55. The attraction is weaker but never vanishes (for bound states).

The Quantum Defect: A Measurable Signature of Penetration

The breaking of $l$-degeneracy in multi-electron atoms has a beautiful experimental signature, particularly visible in alkali metals (Li, Na, K, ...), where a single valence electron orbits outside a closed-shell core.

For an alkali metal, the energy levels of the valence electron can be written as:

$$E_{nl} = -\frac{13.6\;\text{eV}}{(n - \delta_l)^2}$$

where $\delta_l$ is the quantum defect --- a measure of how much the electron "penetrates" the core and sees an enhanced nuclear charge. For hydrogen, $\delta_l = 0$ for all $l$ (no core to penetrate). For sodium:

The $d$ and $f$ electrons in sodium have nearly hydrogen-like energies because they never venture close to the nucleus. The $s$ electrons, by contrast, have wavefunctions with substantial probability density inside the core, feel a much larger effective charge, and are correspondingly more tightly bound. The quantum defect $\delta_l$ decreases monotonically with $l$, directly confirming our penetration argument.

📊 By the Numbers: The quantum defect for sodium's $3s$ valence electron is $\delta_0 = 1.37$, so $E_{3s} = -13.6/(3 - 1.37)^2 = -5.12\;\text{eV}$. The experimental value is $-5.14\;\text{eV}$ --- the agreement is remarkable for such a simple formula.

🔄 Check Your Understanding (Spaced Review --- Ch 5): In hydrogen, the $2s$ and $2p$ orbitals have the same energy. Look at the radial wavefunctions $R_{20}(r)$ and $R_{21}(r)$ from Chapter 5. Which has more probability density near the nucleus? Which would be more affected by the presence of an inner $1s$ electron screen? Use this to explain qualitatively why $\varepsilon_{2s} < \varepsilon_{2p}$ in a multi-electron atom.

A Note on the Born--Oppenheimer Approximation

Throughout this discussion, we have treated the nucleus as fixed at the origin. This is the Born--Oppenheimer approximation, which is justified by the large mass ratio $m_p/m_e \approx 1836$: the nucleus moves so slowly compared to the electrons that, to a good approximation, the electrons instantaneously adjust to the nuclear position. Corrections to this approximation are small (of order $m_e/m_p$) and are usually absorbed into the "reduced mass" correction that we already encountered for hydrogen in Chapter 5. For multi-electron atoms, the Born--Oppenheimer approximation is always used as the starting point.

16.3 Electron Configurations and the Aufbau Principle

The Rules of Orbital Filling

With the central field approximation in place, each electron occupies a single-particle orbital labeled by $(n, l, m_l, m_s)$. The Pauli exclusion principle demands that no two electrons share all four quantum numbers. The question becomes: in what order do the orbitals fill?

The Aufbau principle (from the German Aufbauprinzip, "building-up principle") states: electrons fill orbitals in order of increasing energy $\varepsilon_{nl}$, subject to the Pauli exclusion principle and Hund's rules.

The approximate energy ordering for multi-electron atoms is:

$$1s < 2s < 2p < 3s < 3p < 4s < 3d < 4p < 5s < 4d < 5p < 6s < 4f < 5d < 6p < 7s < 5f < 6d < 7p$$

The Madelung rule (also called the $(n + l)$ rule or Klechkowski rule) provides a useful mnemonic: orbitals fill in order of increasing $n + l$; for equal $n + l$, the orbital with lower $n$ fills first. This gives the ordering above.

📊 By the Numbers --- Subshell Capacity:

Building Up the Periodic Table: The First 18 Elements

Let us apply the Aufbau principle systematically:

The noble gas core notation $[\text{He}]$, $[\text{Ne}]$, $[\text{Ar}]$, etc. abbreviates the filled inner shells.

The First Complications: Potassium Through Krypton

After argon ($Z = 18$), the $3d$ subshell should be next if we follow the hydrogen-like ordering ($3d$ before $4s$). But the Madelung rule tells us that $4s$ ($n + l = 4$) fills before $3d$ ($n + l = 5$). This is a direct consequence of the screening-induced $l$-dependence of orbital energies.

The exceptions for chromium and copper arise from a subtle energetic preference for half-filled ($d^5$) and fully filled ($d^{10}$) $d$-subshells. The exchange interaction (a consequence of antisymmetry) provides extra stabilization when all spins in a subshell are aligned.

⚠️ Common Misconception: "The $4s$ orbital always has lower energy than $3d$." This is true for neutral potassium and calcium but becomes false partway across the transition metal series. In transition metal ions, $3d$ is always lower than $4s$, which is why transition metals lose $4s$ electrons first when ionized (e.g., Fe$^{2+}$ is $[\text{Ar}]\,3d^6$, not $[\text{Ar}]\,3d^4\,4s^2$). The Madelung rule is a guide for neutral atom ground states, not a universal energy ordering.

✅ Checkpoint: Write the ground-state electron configuration for (a) oxygen ($Z = 8$), (b) titanium ($Z = 22$), (c) bromine ($Z = 35$). Identify which subshell is being filled in each case. Check your answers against a periodic table.

16.4 Spectroscopic Notation: Understanding Term Symbols

From Configurations to Terms

Knowing the electron configuration (e.g., $1s^2\,2s^2\,2p^2$ for carbon) is not enough to fully specify the state. The two $2p$ electrons can have their orbital and spin angular momenta arranged in several different ways, each with a different total energy. To classify these arrangements, we use Russell--Saunders term symbols (also called LS-coupling term symbols).

The term symbol encodes three quantum numbers:

$${}^{2S+1}L_J$$

where: - $S$ is the total spin quantum number: $\hat{\mathbf{S}} = \sum_i \hat{\mathbf{s}}_i$, so $S = 0, 1, \ldots$ depending on how individual spins couple. - $L$ is the total orbital angular momentum quantum number: $\hat{\mathbf{L}} = \sum_i \hat{\boldsymbol{\ell}}_i$. It is written as a capital letter: $L = 0 \to S$, $L = 1 \to P$, $L = 2 \to D$, $L = 3 \to F$, $L = 4 \to G$, etc. - $J$ is the total angular momentum quantum number: $\hat{\mathbf{J}} = \hat{\mathbf{L}} + \hat{\mathbf{S}}$, so $|L - S| \leq J \leq L + S$. - $2S + 1$ is the spin multiplicity: the number of possible $M_S$ values (i.e., the number of ways total spin can project onto the $z$-axis). For $S = 0$, the multiplicity is 1 (singlet); for $S = 1$, it is 3 (triplet).

🔗 Connection (Ch 14): The quantum numbers $S$, $L$, $J$ arise from the addition of angular momentum, which we developed in detail in Chapter 14. The Clebsch--Gordan coefficients determine the allowed values of $J$ given $L$ and $S$. If you are uncertain about this coupling, review Section 14.3.

Letter Designation for $L$

(The first four letters come from spectroscopic terminology: Sharp, Principal, Diffuse, Fundamental. Beyond $F$, the letters follow alphabetical order, skipping $J$.)

Worked Example 1: Helium ($1s^2$)

Both electrons are in $1s$: $l_1 = l_2 = 0$ and the subshell is filled.

• $L = 0$ (only possibility when both $l = 0$)
• $S = 0$ (Pauli forces both $m_s$ values to be different: $+\tfrac{1}{2}$ and $-\tfrac{1}{2}$, giving total $S = 0$)
• $J = |L - S| = 0$
• Term symbol: ${}^1S_0$

A filled subshell always gives ${}^1S_0$. This is a general result worth remembering: filled subshells contribute $L = 0$, $S = 0$ to the total. Only the electrons in partially filled subshells matter for the term symbol.

Worked Example 2: Carbon ($1s^2\,2s^2\,2p^2$)

The filled $1s^2$ and $2s^2$ subshells contribute $L = S = 0$. We need only consider the two $2p$ electrons ($l_1 = l_2 = 1$).

Step 1: Find all allowed $L$ values.

Two electrons with $l = 1$ can couple to: $$L = l_1 + l_2, \; l_1 + l_2 - 1, \; \ldots, \; |l_1 - l_2| = 2, 1, 0$$

So $L = 0$ (S), $1$ (P), or $2$ (D).

Step 2: Find all allowed $S$ values.

Two spin-$\tfrac{1}{2}$ particles couple to $S = 0$ (singlet) or $S = 1$ (triplet).

Step 3: Apply Pauli exclusion.

Not all combinations of $L$ and $S$ are allowed for equivalent electrons (electrons in the same subshell). The wavefunction must be antisymmetric overall: since the spatial part has symmetry $(-1)^L$ under electron exchange and the spin part has symmetry $(-1)^{S+1}$, we need $(-1)^{L+S+1} = -1$, which requires $L + S$ to be even.

The allowed terms are ${}^1S$, ${}^3P$, and ${}^1D$.

Step 4: Find $J$ values for each term.

• ${}^1S$: $J = 0$ $\Rightarrow$ ${}^1S_0$
• ${}^3P$: $J = 2, 1, 0$ $\Rightarrow$ ${}^3P_2, {}^3P_1, {}^3P_0$
• ${}^1D$: $J = 2$ $\Rightarrow$ ${}^1D_2$

Verification: Count microstates. Two equivalent $p$-electrons give $\binom{6}{2} = 15$ microstates. Our terms account for $1 + (5 + 3 + 1) + 5 = 15$ states. $\checkmark$

Worked Example 3: Nitrogen ($1s^2\,2s^2\,2p^3$)

Three equivalent $p$-electrons. The number of microstates is $\binom{6}{3} = 20$. The systematic enumeration (using the $M_L, M_S$ table method from Chapter 14) yields the allowed terms:

$${}^4S, \quad {}^2D, \quad {}^2P$$

Verification: $(1) \times 4 + (5) \times 2 + (3) \times 2 = 4 + 10 + 6 = 20$ microstates. $\checkmark$

The key result is that the quartet state ($S = 3/2$, all spins aligned) only occurs with $L = 0$. This is because three electrons with $l = 1$ can have $m_l = +1, 0, -1$ (using all three values), which gives $M_L = 0$, hence $L = 0$. The requirement that all three spins are aligned ($S = 3/2$) forces all spatial quantum numbers to be different (Pauli principle), which forces $M_L = 0$.

🔵 Historical Note: The Russell--Saunders coupling scheme (also called LS-coupling) was developed by Henry Norris Russell and Frederick Saunders in 1925. It works well for light atoms ($Z \lesssim 30$) where the spin--orbit interaction is much weaker than the electrostatic electron--electron interaction. For heavy atoms, the alternative $jj$-coupling scheme (where each electron's $\hat{\boldsymbol{\ell}}_i$ and $\hat{\mathbf{s}}_i$ couple first to form $\hat{\mathbf{j}}_i$, then the individual $\hat{\mathbf{j}}_i$ couple to give $\hat{\mathbf{J}}$) becomes more appropriate. Most atoms fall somewhere between the two extremes.

16.5 Hund's Rules: Finding the Ground State

The Three Rules

Given a set of allowed term symbols for a configuration, which one is the ground state? The answer is provided by Hund's rules (Friedrich Hund, 1927), ordered by priority:

Rule 1 (Maximum $S$): The term with the largest $S$ (highest spin multiplicity $2S + 1$) has the lowest energy.

Physical basis: Electrons with parallel spins must have an antisymmetric spatial wavefunction, which tends to keep them farther apart, reducing the electron--electron repulsion energy. This is the exchange interaction --- a purely quantum-mechanical effect with no classical analogue.

Rule 2 (Maximum $L$, given Rule 1): Among terms with the same $S$, the one with the largest $L$ has the lowest energy.

Physical basis: Larger $L$ means the electrons are orbiting in the same direction, which tends to keep them on opposite sides of the nucleus, again reducing repulsion. The physical argument is more subtle than for Rule 1, but the empirical pattern is clear.

Rule 3 (Minimum $J$ for less than half-filled; Maximum $J$ for more than half-filled): - If the subshell is less than half-filled: $J = |L - S|$ (minimum $J$). - If the subshell is more than half-filled: $J = L + S$ (maximum $J$). - If exactly half-filled: $L = 0$, so $J = S$.

Physical basis: Rule 3 arises from the spin--orbit interaction $\hat{H}_{SO} \propto \hat{\mathbf{L}} \cdot \hat{\mathbf{S}}$. For less than half-filled subshells, the spin--orbit coupling constant is positive (energy increases with $J$), so the lowest $J$ wins. For more than half-filled subshells, the coupling constant is negative (energy decreases with $J$), so the highest $J$ wins.

🔴 Warning: Hund's rules are strictly valid only for the ground state of a free atom. They do not reliably predict the ordering of excited states, nor do they necessarily apply to atoms in molecules or solids where the crystal field or molecular bonding can override these atomic-scale tendencies.

The Fermi Hole: Physical Mechanism Behind Hund's First Rule

To understand why parallel spins lower the energy, consider two electrons in the same subshell. When both spins are parallel (say, both $|\!\uparrow\rangle$), the overall wavefunction must be antisymmetric, so the spatial part is antisymmetric:

$$\psi_{\text{spatial}}(\mathbf{r}_1, \mathbf{r}_2) = \frac{1}{\sqrt{2}}\left[\phi_a(\mathbf{r}_1)\phi_b(\mathbf{r}_2) - \phi_b(\mathbf{r}_1)\phi_a(\mathbf{r}_2)\right]$$

Setting $\mathbf{r}_1 = \mathbf{r}_2 = \mathbf{r}$, we get $\psi_{\text{spatial}}(\mathbf{r}, \mathbf{r}) = 0$. The probability of finding both electrons at the same point in space vanishes. This "exclusion zone" around each electron is called the Fermi hole (or exchange hole). It means that electrons with parallel spins tend to stay apart, reducing their Coulomb repulsion.

For electrons with antiparallel spins (singlet state), the spatial wavefunction is symmetric, and there is no such prohibition --- the electrons can be found at the same point. They are, on average, closer together and experience more repulsion.

This is the essence of the exchange interaction: it is not a new force but rather a correlation in the spatial distribution of electrons imposed by the Pauli principle. Electrons with parallel spins automatically avoid each other, lowering the repulsion energy. This is why Hund's first rule works: the maximum-spin state minimizes the electron--electron repulsion.

Applying Hund's Rules: Worked Examples

Carbon ($2p^2$): Allowed terms are ${}^1S_0$, ${}^3P_{2,1,0}$, ${}^1D_2$.

1. Maximum $S$: ${}^3P$ has $S = 1$ (triplet), highest among all terms. $\checkmark$
2. Only one term with $S = 1$, so Rule 2 is automatically satisfied.
3. Less than half-filled ($2$ of $6$ $p$-electrons): $J = |L - S| = |1 - 1| = 0$.

Ground state of carbon: ${}^3P_0$

Nitrogen ($2p^3$): Allowed terms are ${}^4S_{3/2}$, ${}^2D_{5/2,3/2}$, ${}^2P_{3/2,1/2}$.

1. Maximum $S$: ${}^4S$ has $S = 3/2$ (quartet). $\checkmark$
2. Only one term with $S = 3/2$. Rule 2 is trivially satisfied.
3. Exactly half-filled: $L = 0$, so $J = S = 3/2$.

Ground state of nitrogen: ${}^4S_{3/2}$

Oxygen ($2p^4$): We can think of $2p^4$ as $2p^6$ minus $2p^2$ --- the holes have the same term structure as $2p^2$. The allowed terms are again ${}^1S$, ${}^3P$, ${}^1D$.

1. Maximum $S$: ${}^3P$ with $S = 1$. $\checkmark$
2. Only one triplet term. $\checkmark$
3. More than half-filled ($4$ of $6$): $J = L + S = 1 + 1 = 2$.

Ground state of oxygen: ${}^3P_2$

Note the difference: carbon has ${}^3P_0$ (minimum $J$) and oxygen has ${}^3P_2$ (maximum $J$). The reversal occurs at the half-filled mark. This is a common source of errors on examinations --- always check whether the subshell is less than or more than half-filled before applying Rule 3.

Sulfur ($[\text{Ne}]\,3s^2\,3p^4$): Sulfur has exactly the same valence configuration as oxygen ($p^4$), so the terms and ground state are identical: ${}^3P_2$. This illustrates a general principle: elements in the same group of the periodic table have the same valence electron term symbol structure.

Iron ($[\text{Ar}]\,3d^6\,4s^2$): The filled $4s^2$ contributes $S = L = 0$. The six $3d$ electrons (more than half of 10) give:

1. Maximum $S$: With 6 electrons in 5 orbitals, at most 5 can have parallel spin. So maximum $S = 5/2 \times 2 - 6/2 = 2$. (More carefully: 5 spin-up, 1 spin-down, $M_S = 5/2 - 1/2 = 2$, so $S = 2$.)
2. Maximum $L$ consistent with $S = 2$: Distribute the 5 spin-up electrons among $m_l = 2, 1, 0, -1, -2$ (one each). The sixth electron (spin-down) goes in $m_l = 2$ for maximum $M_L$. So $M_L = (2+1+0-1-2) + 2 = 2$, hence $L = 2$ (D-term).
3. More than half-filled: $J = L + S = 2 + 2 = 4$.

Ground state of iron: ${}^5D_4$

A Summary Table of Ground State Terms

For reference, here are the ground state term symbols for the first-row transition metals, all determined by Hund's rules applied to the $3d^n$ configuration (the filled $4s^2$ contributes nothing):

Note the beautiful pattern: $S$ increases from Sc to Mn (half-filled), then decreases. The $L$ values follow a symmetric pattern around the half-filled point, and $J$ switches from minimum to maximum at the half-filled mark. The term symbol for $d^n$ is related to that for $d^{10-n}$ by the hole equivalence, with the only difference being the $J$ value (Rule 3 reversal).

✅ Checkpoint: Determine the ground state term symbol for (a) sulfur ($3p^4$), (b) vanadium ($3d^3\,4s^2$), (c) cobalt ($3d^7\,4s^2$). Verify your $J$ values using the half-filled rule.

16.6 The Hartree--Fock Method: Self-Consistent Fields

The Central Idea

How do we actually compute the effective potential $V_{\text{eff}}(r)$ used in the central field approximation? The answer, developed by Douglas Hartree (1928) and extended by Vladimir Fock (1930), is the Hartree--Fock method --- a self-consistent procedure that remains the foundation of modern atomic and molecular calculations.

The physical reasoning is beautifully circular:

1. Guess an initial set of single-particle orbitals $\{\psi_i(\mathbf{r})\}$.
2. Compute the average potential that each electron feels due to all the others. This includes: - The direct (Coulomb) interaction: each electron moves in the average electrostatic field of all other electrons. - The exchange interaction: a purely quantum-mechanical correction that arises from the antisymmetry (Slater determinant) requirement.
3. Solve the single-particle Schrodinger equation in this new potential to get updated orbitals $\{\psi_i'(\mathbf{r})\}$.
4. Iterate steps 2--3 until the orbitals stop changing (to within a chosen tolerance). The result is a self-consistent field (SCF).

The Hartree--Fock Equations (Conceptual Form)

For the $i$-th electron, the Hartree--Fock equation takes the form:

$$\left[-\frac{\hbar^2}{2m_e}\nabla^2 - \frac{Ze^2}{4\pi\epsilon_0 r} + \hat{V}_{\text{direct}} + \hat{V}_{\text{exchange}}\right]\psi_i(\mathbf{r}) = \varepsilon_i\,\psi_i(\mathbf{r})$$

The direct (Hartree) potential is a local potential:

$$V_{\text{direct}}(\mathbf{r}) = \frac{e^2}{4\pi\epsilon_0}\sum_{j \neq i}\int\frac{|\psi_j(\mathbf{r}')|^2}{|\mathbf{r} - \mathbf{r}'|}\,d^3\mathbf{r}'$$

This has a transparent interpretation: the charge density of electron $j$ is $-e|\psi_j(\mathbf{r}')|^2$, and the integral computes the electrostatic potential at $\mathbf{r}$ due to this charge distribution. The sum runs over all occupied orbitals except the $i$-th.

The exchange potential is nonlocal:

$$\hat{V}_{\text{exchange}}\,\psi_i(\mathbf{r}) = -\frac{e^2}{4\pi\epsilon_0}\sum_{j \neq i}^{\text{same spin}}\int\frac{\psi_j^*(\mathbf{r}')\psi_i(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|}\,d^3\mathbf{r}'\;\psi_j(\mathbf{r})$$

This term has no classical analogue --- it arises purely from the antisymmetry of the wavefunction (the Slater determinant). Notice that it is nonlocal: the operator acting on $\psi_i(\mathbf{r})$ depends on the value of $\psi_i$ at all points $\mathbf{r}'$, not just at $\mathbf{r}$.

💡 Key Insight: The exchange term is responsible for the extra stability of half-filled and filled subshells that we observed in Section 16.3 (the exceptions for Cr and Cu). Electrons with parallel spin experience an exchange interaction that lowers their combined energy relative to anti-parallel spin. This is the quantum-mechanical origin of Hund's first rule.

What Hartree--Fock Gets Right --- and What It Misses

Hartree--Fock captures approximately 99% of the total energy of an atom. For most purposes in atomic physics and chemistry, this is remarkably good. It correctly predicts:

• Electron configurations and the ordering of orbital energies
• Ionization energies (to within ~1--5%)
• Qualitative trends across the periodic table
• The shell structure of atoms

What it misses is the correlation energy --- the difference between the exact non-relativistic energy and the Hartree--Fock energy:

$$E_{\text{corr}} = E_{\text{exact}} - E_{\text{HF}}$$

For helium, $E_{\text{HF}} = -77.9\;\text{eV}$ and $E_{\text{exact}} = -79.0\;\text{eV}$, so $E_{\text{corr}} = -1.1\;\text{eV}$. This 1.4% error matters enormously for chemical bonding, which typically involves energy differences of a few eV. Capturing correlation energy requires methods beyond Hartree--Fock: configuration interaction (CI), coupled cluster theory (CCSD(T)), or density functional theory (DFT).

🧪 Experiment: The first successful Hartree--Fock calculations were performed by Hartree himself, for rubidium and other alkali metals, in the late 1920s --- using a hand-cranked mechanical calculator. Each iteration took weeks. Charlotte Froese Fischer, starting in the 1960s, developed the first general-purpose computer programs for multi-configuration Hartree--Fock calculations, and her work remains foundational.

The Self-Consistent Field in Pictures

The SCF procedure can be visualized as a convergent iteration:

Iteration 0:  Guess orbitals → Compute potential → Solve equations → New orbitals
Iteration 1:  New orbitals   → Compute potential → Solve equations → Newer orbitals
Iteration 2:  Newer orbitals → Compute potential → Solve equations → Even newer orbitals
      ...         ...              ...                  ...              ...
Convergence:  Final orbitals  = Self-consistent solution

Typically, convergence is achieved in 10--30 iterations for atoms. The initial guess is often a set of hydrogen-like orbitals with a screened nuclear charge (Slater-type orbitals).

Beyond Hartree--Fock: A Brief Roadmap

For completeness, we mention the major approaches that go beyond Hartree--Fock:

Configuration Interaction (CI): Instead of a single Slater determinant, the wavefunction is expanded as a linear combination of many Slater determinants built from different orbital occupations ("configurations"). Full CI is exact in a given basis set but scales exponentially with the number of electrons. Truncated CI (e.g., CISD, including only single and double excitations) is practical for small systems.

Coupled Cluster Theory (CCSD(T)): Widely regarded as the "gold standard" of quantum chemistry. It uses an exponential ansatz $|\Psi\rangle = e^{\hat{T}}|\Phi_0\rangle$ where $\hat{T}$ generates cluster excitations. CCSD(T) achieves chemical accuracy (errors below $\sim 1\;\text{kcal/mol} \approx 0.04\;\text{eV}$) for most molecular properties.

Density Functional Theory (DFT): Rather than computing the many-electron wavefunction, DFT works directly with the electron density $\rho(\mathbf{r})$. The Hohenberg--Kohn theorem (1964) guarantees that the ground state energy is a unique functional of the density. In practice, the Kohn--Sham formalism replaces the many-body problem with an auxiliary set of non-interacting electrons in an effective potential. DFT is by far the most widely used electronic structure method due to its favorable cost-to-accuracy ratio, though the exact density functional is unknown and must be approximated.

These methods are the subject of graduate courses in theoretical chemistry and condensed matter physics. For our purposes, the conceptual framework of Hartree--Fock is sufficient: a self-consistent mean-field theory that captures shell structure, orbital energies, and the qualitative physics of multi-electron atoms.

16.7 Building the Periodic Table from Quantum Mechanics

The Quantum Mechanical Origin of Periodic Structure

We now have all the ingredients to explain the periodic table from first principles:

1. Quantized orbitals: The central field approximation produces orbitals labeled by $(n, l, m_l)$ with energies $\varepsilon_{nl}$.
2. Pauli exclusion: Each orbital (including spin) can hold at most one electron, giving subshell capacities of $2(2l + 1)$.
3. Aufbau principle: Orbitals fill in order of increasing $\varepsilon_{nl}$.
4. Hund's rules: Within a partially filled subshell, electrons maximize $S$, then $L$, then follow the half-filling rule for $J$.

The periods of the periodic table correspond to the filling of successive shells:

The period lengths follow the pattern $2, 8, 8, 18, 18, 32, 32$ --- that is, $2n^2$ for each principal quantum number $n$, with each value appearing twice (because $s$-orbitals from two successive $n$ values bracket each group of higher-$l$ orbitals).

💡 Key Insight: The fact that the period lengths are $2, 8, 8, 18, 18, 32, 32$ rather than $2, 8, 18, 32, \ldots$ is a direct consequence of the Madelung rule ($n + l$ ordering). The "$l$-degeneracy lifting" interleaves orbitals from different principal quantum numbers, creating the characteristic "double period" structure of the periodic table. If atoms were truly hydrogen-like (with only $n$-dependent energies), the periodic table would have a completely different structure with periods of length $2, 8, 18, 32, \ldots$.

Explaining Chemical Periodicity

The groups (vertical columns) of the periodic table contain elements with the same valence electron configuration. This is why they share chemical properties:

• Group 1 (Alkali metals): $ns^1$ valence configuration. One loosely bound electron outside a noble gas core. Easily ionized; form $+1$ cations. Ionization energies decrease down the group as the valence electron is farther from the nucleus.

• Group 2 (Alkaline earth metals): $ns^2$. Two $s$-electrons. Less reactive than Group 1 (paired $s$ electrons are harder to remove).

• Groups 3--12 (Transition metals): Filling of $(n-1)d$ subshell while $ns$ is already occupied. The $d$-electrons have similar energies, leading to variable oxidation states, colored ions, and catalytic properties.

• Group 17 (Halogens): $ns^2\,np^5$. One electron short of a filled $p$ subshell. High electron affinities; form $-1$ anions.

• Group 18 (Noble gases): $ns^2\,np^6$ (except He: $1s^2$). Filled $s$ and $p$ subshells. Very stable; high ionization energies; low reactivity.

Ionization Energy Trends

The first ionization energy $E_I$ --- the energy required to remove the outermost electron --- provides a stringent test of the quantum mechanical picture. The general trend is:

• $E_I$ increases across a period (left to right): the nuclear charge increases while the screening is imperfect, so the effective nuclear charge on the valence electron grows.
• $E_I$ decreases down a group (top to bottom): the valence electron is in a higher $n$ orbital, farther from the nucleus.

But there are deviations from the smooth trend that quantum mechanics explains precisely:

• B ($Z = 5$) has lower $E_I$ than Be ($Z = 4$): Beryllium's $2s^2$ configuration has both electrons in a well-shielded $s$-orbital. Boron's extra electron goes into $2p$, which is higher in energy (less penetrating) and therefore easier to remove.

• O ($Z = 8$) has lower $E_I$ than N ($Z = 7$): Nitrogen has a half-filled $2p^3$ configuration (maximum exchange stabilization). Oxygen's fourth $2p$ electron must pair with an existing one, experiencing additional electron--electron repulsion that makes it easier to remove.

📊 By the Numbers --- First Ionization Energies (eV):

Note the dips at B and O, exactly as predicted by the quantum mechanical analysis.

Atomic Radius Trends

Atomic radii follow a complementary pattern:

• Decreasing across a period: More protons, imperfect screening → stronger attraction → smaller atom.
• Increasing down a group: Higher $n$ → larger orbitals → bigger atom.
• $d$-block contraction: The $3d$ electrons are poor screeners (low penetration), so $Z_{\text{eff}}$ increases significantly across the transition metals, causing atomic radii to decrease. This makes the $4d$ elements (Y through Cd) have similar sizes to the $3d$ elements.
• Lanthanide contraction: The $4f$ electrons are even poorer screeners than $3d$, producing a cumulative contraction across the lanthanide series that makes the $5d$ elements (La through Hg) have nearly the same radii as the $4d$ elements.

The lanthanide contraction has far-reaching chemical consequences: it explains why Zr and Hf have almost identical chemistry (despite being separated by 32 elements), why gold is dense and unreactive, and why mercury is liquid at room temperature.

Electron Affinity and Electronegativity

The electron affinity $E_A$ is the energy released when an atom gains an electron: $X + e^- \to X^-$. Quantum mechanics predicts:

• Halogens (group 17) have the highest electron affinities: gaining one electron completes the $p$-subshell, achieving a noble gas configuration. Fluorine: $E_A = 3.40\;\text{eV}$; chlorine: $E_A = 3.61\;\text{eV}$.
• Noble gases (group 18) have approximately zero electron affinities: an extra electron would go into the next shell ($n+1$), which is weakly bound.
• Group 2 (alkaline earths) and group 12 ($ns^2$ configurations) have very low electron affinities: the added electron goes into a $p$-orbital, which is screened and weakly bound.

A surprising prediction: nitrogen ($2p^3$, half-filled) has a lower electron affinity ($-0.07\;\text{eV}$, actually slightly negative, meaning N$^-$ is unstable) than carbon ($2p^2$, $E_A = 1.26\;\text{eV}$). Adding an electron to nitrogen forces spin pairing in the $2p$ subshell, destroying the favorable half-filled exchange stabilization. This subtlety is correctly predicted by the quantum mechanical model and is inexplicable by any classical theory.

Electronegativity --- the tendency of an atom in a molecule to attract electrons --- is not a ground-state atomic property but rather depends on the bonding environment. Linus Pauling's famous electronegativity scale, Mulliken's scale (based on the average of $E_I$ and $E_A$), and modern Allen electronegativity (based on average ionization energies of valence electrons) all reflect the underlying quantum mechanical orbital structure. The general trends --- increasing across a period, decreasing down a group --- mirror the trends in $Z_{\text{eff}}$ that we have been analyzing.

A Thought Experiment: The Periodic Table Without Spin

It is instructive to consider what the periodic table would look like if electrons had no spin (or equivalently, if we ignored the Pauli exclusion principle). Without the exclusion principle, every electron would fall into the $1s$ orbital. All atoms would be spherically symmetric, chemically inert, and essentially identical except for their mass and size. There would be no shell structure, no periodicity, no chemistry, and no life.

Even if we retained the exclusion principle but removed spin, each orbital could hold only one electron (since $m_s$ would not exist). The subshell capacities would halve: $s$ holds 1, $p$ holds 3, $d$ holds 5, $f$ holds 7. The periodic table would have periods of length 1, 4, 4, 9, 9, 16, 16. The first noble gas would be hydrogen ($1s^1$). Helium would be an alkali metal. Carbon would be a noble gas ($1s^1\,2s^1\,2p^3$). The entire landscape of chemistry would be unrecognizable.

⚖️ Interpretation: The periodic table --- and with it, all of chemistry, biology, and the material complexity of the universe --- depends critically on the existence of electron spin and the Pauli exclusion principle. These are not incidental features of quantum mechanics. They are the structural foundations upon which the diversity of matter is built.

16.8 Screening, Shielding, and the Failure of the Hydrogen-Like Model

Slater's Rules: A Quantitative Screening Model

In 1930, John C. Slater proposed a simple empirical scheme for estimating the effective nuclear charge experienced by each electron. Slater's rules assign screening constants $\sigma$ based on the orbital group of the electron:

1. Write the electron configuration in the following groups: $(1s)(2s,2p)(3s,3p)(3d)(4s,4p)(4d)(4f)\ldots$
2. For an electron in group $(ns, np)$: - Each other electron in the same group contributes $0.35$ (except $1s$, where it is $0.30$). - Each electron in the $(n-1)$ group contributes $0.85$. - Each electron in groups with $n' \leq n-2$ contributes $1.00$.
3. For an electron in group $(nd)$ or $(nf)$: - Each other electron in the same group contributes $0.35$. - Each electron in any inner group contributes $1.00$.

The effective nuclear charge is then $Z_{\text{eff}} = Z - \sigma$.

Example: Carbon ($Z = 6$, configuration $(1s^2)(2s^2\,2p^2)$):

For a $2p$ electron: - Same group: 3 other electrons $\times\,0.35 = 1.05$ - $(1s)$ group: 2 electrons $\times\,0.85 = 1.70$ - $\sigma = 2.75$; $Z_{\text{eff}} = 6 - 2.75 = 3.25$

For a $1s$ electron: - Same group: 1 other electron $\times\,0.30 = 0.30$ - $\sigma = 0.30$; $Z_{\text{eff}} = 6 - 0.30 = 5.70$

The inner $1s$ electrons see nearly the full nuclear charge, while the outer $2p$ electrons see an effective charge of about $3.25$. This quantifies the screening picture we have been using qualitatively.

Effective Quantum Numbers and Orbital Energies

Slater also proposed effective quantum numbers $n^*$:

The approximate orbital energy is then:

$$\varepsilon_{nl} \approx -13.6\;\text{eV}\times\frac{Z_{\text{eff}}^2}{(n^*)^2}$$

This gives a rough but physically motivated estimate of orbital energies across the periodic table. For example, the $2p$ energy in carbon is approximately:

$$\varepsilon_{2p} \approx -13.6\times\frac{3.25^2}{4} = -35.9\;\text{eV}$$

The experimental value (from photoelectron spectroscopy) is about $-11.3\;\text{eV}$ for the valence $2p$ electrons. The Slater estimate is too negative because the simple model overestimates the nuclear attraction felt by $2p$ electrons (Slater's rules ignore the difference between $2s$ and $2p$ screening). This illustrates both the utility and the limitations of empirical screening models.

Why the Hydrogen-Like Model Fails Quantitatively

The fundamental failure of the hydrogen-like model for multi-electron atoms can be summarized:

1. No $l$-dependence of energy: Hydrogen has $E_n$ depending only on $n$. Multi-electron atoms have $\varepsilon_{nl}$ depending on both $n$ and $l$. This changes the orbital filling order and the structure of the periodic table.

2. No exchange interaction: The hydrogen-like model treats each electron independently. It cannot capture the extra stability of parallel-spin configurations (Hund's first rule) or the subtle preference for half-filled subshells.

3. Correlation neglected entirely: Even the Hartree--Fock method, which includes exchange, misses correlation --- the instantaneous, dynamic avoidance between electrons that goes beyond the mean-field approximation.

4. Relativistic effects absent: For heavy atoms ($Z \gtrsim 50$), relativistic contraction of $s$ and $p$ orbitals and relativistic expansion of $d$ and $f$ orbitals become significant. Gold's color, mercury's liquidity, and lead's inertness (the "inert pair effect") all require relativistic corrections to explain.

⚠️ Common Misconception: "Quantum mechanics cannot explain chemistry --- it is just too complicated." This is backwards. Quantum mechanics explains chemistry completely in principle. The Schrodinger equation (or the Dirac equation for relativistic effects) contains all the information. The challenge is computational, not conceptual. Modern computational chemistry routinely achieves "chemical accuracy" (errors below 1 kcal/mol $\approx 0.04$ eV) for systems with hundreds of atoms.

Relativistic Effects: When $v/c$ Is Not Negligible

For heavy atoms ($Z \gtrsim 50$), the innermost electrons move at speeds comparable to the speed of light. A simple estimate: for a $1s$ electron in a hydrogen-like ion, the average speed is $v \sim Zc\alpha$, where $\alpha \approx 1/137$ is the fine-structure constant. For gold ($Z = 79$), $v_{1s}/c \approx 79/137 \approx 0.58$, and the relativistic Lorentz factor is $\gamma \approx 1.22$. The $1s$ orbital contracts by roughly 22% due to relativistic mass increase.

This contraction propagates outward: contracted inner $s$ and $p$ orbitals provide more effective screening of the nucleus from outer $d$ and $f$ orbitals, causing the latter to expand. The interplay between relativistic contraction (of $s, p$) and expansion (of $d, f$) has remarkable chemical consequences:

• Gold is yellow: The $5d \to 6s$ absorption edge shifts from the ultraviolet (in silver) into the blue due to the relativistic lowering of the $6s$ orbital. Gold absorbs blue light and reflects yellow.
• Mercury is liquid: Relativistic contraction of the $6s^2$ pair creates an exceptionally tight, inert orbital that resists metallic bonding.
• Lead is inert: The "inert pair effect" --- the reluctance of $6s^2$ electrons in Tl, Pb, Bi to participate in bonding --- is a direct consequence of relativistic stabilization.

These are not exotic fringe effects. They determine the color of gold, the state of mercury, and the chemistry of every $5d$ and $6p$ element. A fully non-relativistic periodic table would incorrectly predict gold to be silver-colored, mercury to be solid, and lead to be as reactive as tin.

🔗 Connection (Forward): We will develop the Dirac equation (the relativistic quantum theory of the electron) in Chapter 29. There, we will see how relativistic effects are not merely corrections to Schrodinger theory but arise naturally from requiring the quantum theory to be consistent with special relativity.

🔄 Check Your Understanding (Spaced Review --- Ch 14): In Chapter 14, we learned to add angular momenta using Clebsch--Gordan coefficients. Write the state $|L=1, M_L=0\rangle$ for two $p$-electrons as a linear combination of $|m_{l_1}, m_{l_2}\rangle$ states. What constraint does the Pauli principle impose on which combinations are allowed?

16.9 Summary and Project Checkpoint

What We Have Learned

This chapter has taken us from the exactly solvable hydrogen atom to the rich, approximate, but extraordinarily powerful theory of multi-electron atoms. The key ideas are:

1. The multi-electron Hamiltonian includes electron--electron repulsion terms that prevent exact solutions.
2. The central field approximation replaces the full many-body problem with effective single-particle equations in a spherically symmetric potential.
3. The $l$-degeneracy is broken: unlike hydrogen, multi-electron atoms have energies that depend on both $n$ and $l$, with lower $l$ orbitals having lower energy (more nuclear penetration).
4. Electron configurations follow the Aufbau principle and the Madelung $(n+l)$ rule, with exceptions for transition metals and rare earths.
5. Term symbols ${}^{2S+1}L_J$ classify the states of partially filled subshells by their total spin, orbital, and total angular momentum.
6. Hund's rules identify the ground state: maximize $S$, then $L$, then choose $J$ by the half-filling rule.
7. The Hartree--Fock method provides a systematic, self-consistent procedure for computing orbital energies and wavefunctions.
8. The periodic table --- its periods, groups, and chemical trends --- emerges directly from quantum mechanics.

The Bigger Picture

The periodic table is not an empirical curiosity that quantum mechanics happens to explain. It is the inevitable consequence of four quantum mechanical ingredients: (i) the Schrodinger equation with the Coulomb potential, (ii) quantization of angular momentum, (iii) electron spin, and (iv) the Pauli exclusion principle. Strip away any one of these four, and the periodic table would not exist.

Mendeleev arranged the elements by chemical properties in 1869, noting periodic patterns and predicting missing elements. His table was brilliant empiricism. Quantum mechanics, arriving half a century later, provided the explanation: the chemical periodicity of the elements is a manifestation of the shell structure of electrons in atoms, which in turn follows from the quantum numbers $n$, $l$, $m_l$, $m_s$ and the exclusion principle.

🔗 Connection (Forward): In Chapter 17, we will develop time-independent perturbation theory, which provides a systematic framework for computing corrections to the central field approximation. In Chapter 18, we will apply these tools to the fine structure of hydrogen --- the spin--orbit, relativistic, and Darwin corrections that split the $l$-degenerate levels. And in Chapter 19, the variational method will give us a powerful alternative approach, with the helium atom as our proving ground.

Project Checkpoint: Quantum Toolkit v1.6 --- Atomic Physics Module

What you are building: A multi-electron atom module that computes electron configurations, generates term symbols, applies Hund's rules to identify ground states, and implements a simplified Hartree solver for helium-like atoms.

New capabilities added in this checkpoint:

1. electron_config(Z) --- returns the ground state electron configuration for element with atomic number $Z$, following the Aufbau/Madelung rule.

2. term_symbols(electrons, l) --- given the number of equivalent electrons and their $l$ value, returns all allowed Russell--Saunders term symbols ${}^{2S+1}L_J$.

3. hund_ground_state(terms, n_electrons, max_electrons) --- selects the ground state term from the list using Hund's three rules.

4. slater_zeff(Z, orbital) --- computes the Slater effective nuclear charge for a specified orbital.

5. hartree_helium(Z, tol) --- simplified Hartree solver for helium-like atoms, iterating to self-consistency.

Acceptance criteria: - electron_config(26) returns "[Ar] 3d^6 4s^2" (iron) - term_symbols(2, 1) returns ['1S_0', '3P_2,1,0', '1D_2'] (carbon $p^2$) - hund_ground_state(...) selects ${}^3P_0$ for carbon - hartree_helium(2) converges to $E \approx -77.9\;\text{eV}$ (Hartree--Fock limit for He)

Run code/project-checkpoint.py to verify all tests pass.

Key Equations Summary

Multi-Electron Hamiltonian

$$\hat{H} = \sum_{i=1}^{N}\left[-\frac{\hbar^2}{2m_e}\nabla_i^2 - \frac{Ze^2}{4\pi\epsilon_0 r_i}\right] + \sum_{i

Central Field Single-Particle Equation

$$\left[-\frac{\hbar^2}{2m_e}\nabla^2 + V_{\text{eff}}(r)\right]\psi_{nlm}(\mathbf{r}) = \varepsilon_{nl}\,\psi_{nlm}(\mathbf{r})$$

Effective Nuclear Charge (Slater)

$$Z_{\text{eff}} = Z - \sigma$$

Term Symbol

$${}^{2S+1}L_J, \quad |L - S| \leq J \leq L + S$$

Hartree Direct Potential

$$V_{\text{direct}}(\mathbf{r}) = \frac{e^2}{4\pi\epsilon_0}\sum_{j \neq i}\int\frac{|\psi_j(\mathbf{r}')|^2}{|\mathbf{r} - \mathbf{r}'|}\,d^3\mathbf{r}'$$

Subshell Capacity

$$\text{Max electrons} = 2(2l + 1)$$

Orbital Energy Estimate (Slater)

$$\varepsilon_{nl} \approx -13.6\;\text{eV}\times\frac{Z_{\text{eff}}^2}{(n^*)^2}$$

In the next chapter, we develop the systematic approximation tools --- perturbation theory --- that let us move beyond the central field approximation and compute corrections to increasing precision.