Hydrogen atom quantum chemistry

Intuition [Beginner]

The hydrogen atom -- one proton, one electron -- is the only atom with an exact quantum-mechanical solution. Every orbital in every atom is named after a hydrogen orbital. Understanding hydrogen is the key to understanding all of chemistry.

The electron in hydrogen does not orbit like a planet. Instead, it occupies an orbital: a three-dimensional region where the electron is likely to be found. Each orbital is labelled by three quantum numbers. The principal quantum number $n=1,2,3,\dots$ sets the energy and the overall size. The angular momentum quantum number $l=0,1,\dots ,n-1$ sets the shape: $l=0$ is an $s$ orbital (spherical), $l=1$ is a $p$ orbital (dumbbell), $l=2$ is a $d$ orbital (cloverleaf). The magnetic quantum number $m_l=-l,\dots ,+l$ sets the orientation.

The fourth quantum number is $m_s=\pm 1/2$ (spin), which does not come from the orbital itself but from an intrinsic property of the electron. Together, these four numbers fully specify an electron's state.

Visual [Beginner]

The shapes of orbitals are among the most recognisable images in science. An $s$ orbital is a sphere centred on the nucleus. A $p$ orbital has two lobes with a node (a region of zero probability) at the nucleus. There are three $p$ orbitals ($p_x$, $p_y$, $p_z$) oriented along the three axes. The five $d$ orbitals have four lobes (most of them) or a unique torus-plus-lobe shape ($d_{z^2}$).

The radial distribution function $4\pi r^2|\psi {|}^2$ tells you how much probability lies in a thin spherical shell at distance $r$ from the nucleus. For the $1s$ orbital this peaks at $r=a_0=52.9$ pm (the Bohr radius). For the $2s$ orbital it peaks much further out, reflecting the larger size of the $n=2$ shell.

Worked example [Beginner]

Chlorine has $Z=17$. Describe its $3p$ valence orbitals.

Step 1. Chlorine's electron configuration is $1s^{2}2s^{2}2p^{6}3s^{2}3p^5$. The valence electrons are in $n=3$: two in $3s$ and five in $3p$.

Step 2. The $3p$ subshell has $l=1$, so $m_l=-1,0,+1$. That gives three $3p$ orbitals: $3p_x$, $3p_y$, $3p_z$. Each holds up to 2 electrons. With 5 electrons in three orbitals, Hund's rule gives one electron in each orbital and a paired electron in one of them.

Step 3. Nodes. The $3p$ orbital has $n-l-1=3-1-1=1$ radial node (a spherical shell where the wave function passes through zero) and $l=1$ angular node (a plane through the nucleus -- the $xy$-plane for $p_z$, etc.). Total nodes: $n-1=2$.

Step 4. Chlorine has 7 valence electrons ($3s^{2}3p^5$) because the $3p$ subshell is one electron short of being full ($3p^6$ would complete it). This near-full shell makes chlorine highly electronegative -- it strongly attracts one more electron to complete the octet.

Check your understanding [Beginner]

Formal definition [Intermediate+]

The hydrogenic atom (nuclear charge $Ze$, one electron) is governed by the Hamiltonian

$$\hat{H}=-\frac{{\hslash }^2}{2m_e}{\nabla }^2-\frac{Z{e}^2}{4\pi {\epsilon }_{0}r}.$$

The bound-state solutions are labelled by three quantum numbers $(n,l,m_l)$ with $n=1,2,3,\dots$, $l=0,1,\dots ,n-1$, $m_l=-l,\dots ,+l$. The energy eigenvalues depend only on $n$:

$$E_n=-\frac{Z^2{R}_H}{n^2},R_H=\frac{m_e{e}^4}{8{\epsilon }_0^2{h}^2}\approx 13.6\text{eV}.$$

The $n$-fold degeneracy (for each $n$, all allowed $l$ and $m_l$ give the same energy) is specific to the $1/r$ potential and is a consequence of the $\mathrm{SO}(4)$ symmetry of the Coulomb problem. This degeneracy is broken in multi-electron atoms by electron-electron repulsion.

The wave functions separate into radial and angular parts:

$${\psi }_{nl{m}_l}(r,\theta ,\phi )=R_{nl}(r)Y_l^{m_l}(\theta ,\phi ),$$

where $R_{nl}(r)$ is the radial wave function (expressed in terms of associated Laguerre polynomials) and $Y_l^{m_l}$ are the spherical harmonics.

The radial distribution function is

$$P(r)=r^2|R_{nl}(r){|}^2,$$

giving the probability per unit $r$ of finding the electron at distance $r$ from the nucleus (integrated over all angles). The factor $r^2$ accounts for the increasing volume of spherical shells at larger $r$. The most probable distance for the $1s$ electron is $r_{\max }=a_0/Z$, where $a_0=52.9$ pm is the Bohr radius.

Nodes. A radial node is a value of $r>0$ where $R_{nl}(r)=0$. The number of radial nodes is $n-l-1$. An angular node is a surface (plane or cone) where $Y_l^{m_l}=0$; the count is $l$. Total nodes: $n-1$.

Orbital notation

$l$ value	Letter	Name	Orbitals in subshell	Max electrons
0	$s$	sharp	1	2
1	$p$	principal	3	6
2	$d$	diffuse	5	10
3	$f$	fundamental	7	14

The letters originate from the appearance of spectral lines: sharp, principal, diffuse, and fine/fundamental series observed in alkali metal spectra.

Key theorem with proof [Intermediate+]

Theorem (Hydrogenic energy levels). The bound-state energies of the hydrogenic Hamiltonian $\hat{H}=-\frac{{\hslash }^2}{2m_e}{\nabla }^2-\frac{Z{e}^2}{4\pi {\epsilon }_{0}r}$ are

$$E_n=-\frac{m_e{Z}^2{e}^4}{2(4\pi {\epsilon }_0{)}^2{\hslash }^2{n}^2}=-\frac{Z^2{R}_H}{n^2},n=1,2,3,\dots$$

with degeneracy $n^2$ (ignoring spin) or $2n^2$ (including spin).

Proof sketch. The proof proceeds by separation of variables in spherical coordinates. The angular equation yields the spherical harmonics $Y_l^{m_l}$ with quantised angular momentum $L^2={\hslash }^{2}l(l+1)$ and $L_z=\hslash m_l$. The radial equation becomes

$$-\frac{{\hslash }^2}{2m_e}\frac{1}{r^2}\frac{d}{dr}\left(r^2\frac{dR}{dr}\right)+\left[\frac{{\hslash }^{2}l(l+1)}{2m_e{r}^2}-\frac{Z{e}^2}{4\pi {\epsilon }_{0}r}\right]R=ER.$$

The substitution $u(r)=rR(r)$ and the introduction of the effective potential $V_{\text{eff}}(r)=-Z{e}^2/(4\pi {\epsilon }_{0}r)+{\hslash }^{2}l(l+1)/(2m_e{r}^2)$ reduce this to a one-dimensional Schrodinger equation. The boundary conditions ($u(0)=0$, $u\to 0$ as $r\to \infty$) quantise the energy. The bound-state solutions exist only for $E<0$, and the requirement of normalisability (the wave function must vanish at infinity) restricts $n$ to positive integers. The detailed calculation involves the associated Laguerre polynomials and parallels the treatment in unit 12.06.01 pending. $\square$

Corollary. For multi-electron atoms, the energy depends on both $n$ and $l$. The $n$-fold degeneracy of hydrogen is lifted by electron-electron repulsion and shielding effects. The energy ordering follows the empirical $(n+l)$ rule (Madelung rule) rather than $n$ alone.

Worked example at intermediate level

Compute the number of angular and radial nodes for the $3d$ orbital of hydrogen, and verify that the most probable radius is $9a_0/Z$.

For $3d$: $n=3$, $l=2$. Angular nodes: $l=2$. Radial nodes: $n-l-1=0$. The $3d$ orbital has no radial nodes -- the radial function is nonzero for all $r>0$ except at $r=0$ and $r\to \infty$.

The radial distribution function for $3d$ peaks at $r_{\max }=n^2{a}_0/Z=9a_0/Z$. For hydrogen ($Z=1$), $r_{\max }=9\times 52.9=476$ pm. The general formula for the most probable radius of an orbital with quantum numbers $(n,l)$ in a hydrogenic atom is

$$r_{\max }=\frac{n^2{a}_0}{Z}\left[1+\frac{1}{2}\left(1-\frac{l(l+1)}{n^2}\right)\right],$$

which simplifies to $n^2{a}_0/Z$ when $l=n-1$ (the maximum-$l$ case, where there are no radial nodes and the distribution has a single peak).

Exercises [Intermediate+]

Relativistic corrections and fine structure [Master]

The Schrödinger treatment of hydrogen produces energies that depend only on $n$. High-resolution spectroscopy contradicts this: every spectral line resolves into a narrow multiplet at finer wavelength resolution. The Balmer-alpha line at 656 nm is not a single line but a doublet, and at still higher resolution a multiplet of seven components. The splittings are tiny -- of order $10^{-4}$ eV -- but they are real, and they encode physics that lies outside non-relativistic quantum mechanics.

The dimensionless small parameter is the fine-structure constant $$ \alpha = \frac{e^2}{4\pi\varepsilon_0\hbar c} \approx \frac{1}{137.036}. $$ For a hydrogenic atom of charge $Z$, the ratio of the electron's orbital speed to the speed of light is $v/c\sim Z\alpha$. For hydrogen this is one part in 137; for heavier hydrogenic ions like ${\text{Au}}^{78+}$ ($Z=79$) the ratio approaches one half, and relativistic effects dominate. The leading correction to the Schrödinger energies is of order ${\alpha }^2{Z}^4$ relative to the Rydberg energy, giving relative splittings $\sim 10^{-4}$ for hydrogen and growing as $Z^4$ for heavier nuclei. Relativistic chemistry is not academic: the colour of gold, the liquid state of mercury at room temperature, and the contracted $6s$ orbital that suppresses the basicity of the heavy alkali atoms are all consequences of $Z^4$ scaling of relativistic corrections.

Three distinct corrections combine to produce the fine structure. The relativistic kinetic-energy correction $-{\hat{p}}^4/(8m_e^3{c}^2)$ arises from expanding the relativistic kinetic energy $\sqrt{p^2{c}^2+m^2{c}^4}-m{c}^2$ to one order beyond $p^2/2m$. The spin-orbit coupling $\xi (r)\hat{\mathbf{L}}\cdot \hat{\mathbf{S}}$ arises from the magnetic field that an electron experiences in its own rest frame as it orbits the proton's electrostatic field. The Darwin term $({\hslash }^2/(8m_e^2{c}^2)){\nabla }^{2}V$ is a contact correction localised at the nucleus and affects only $s$-states. All three are equally large for hydrogen, and they conspire so that the total fine-structure correction depends only on the total angular momentum $j=l\pm 1/2$, not on $l$ separately. The remarkable result is that the energy $$ E_{n,j} = -\frac{Z^2 R_H}{n^2}\left[1 + \frac{(Z\alpha)^2}{n}\left(\frac{1}{j + 1/2} - \frac{3}{4n}\right)\right] $$ depends only on $n$ and $j$. States with the same $j$ but different $l$ remain degenerate at this order -- a cancellation that is mysterious from the Schrödinger viewpoint and obvious from the Dirac viewpoint.

The Dirac equation $(i\hslash {\gamma }^{\mu }{\partial }_{\mu }-mc)\psi =0$, when reduced to the non-relativistic limit by the Foldy-Wouthuysen transformation, reproduces all three corrections in a single stroke. Dirac (1928) solved the hydrogen atom exactly within his equation and obtained $$ E_{n,j}^{\text{Dirac}} = mc^2\left[1 + \left(\frac{Z\alpha}{n - (j + 1/2) + \sqrt{(j + 1/2)^2 - (Z\alpha)^2}}\right)^2\right]^{-1/2}, $$ whose expansion in $Z\alpha$ recovers the Schrödinger formula at zeroth order, the fine-structure formula above at $O({\alpha }^2)$, and additional corrections at higher orders. The Dirac formulation makes spin an intrinsic property of the four-component wave function rather than an ad hoc addition. The chemist usually meets these corrections only through their spectroscopic consequences, but the underlying mathematical structure is Dirac's, not Schrödinger's.

In spectroscopic notation, atomic states are labelled by term symbols ${}^{2S+1}{L}_J$, where $L$ is the total orbital angular momentum, $S$ is the total spin, and $J$ is the total angular momentum. The $2P$ level of hydrogen ($n=2$, $l=1$) splits into ${}^2{P}_{1/2}$ ($j=1/2$) and ${}^2{P}_{3/2}$ ($j=3/2$). The splitting is $\Delta E=(Z\alpha {)}^2{R}_H/(16)\approx 4.5\times 10^{-5}$ eV, corresponding to a wavelength shift of about 0.014 nm between the two components of Balmer-alpha. The $2S_{1/2}$ level lies at the same fine-structure energy as $2P_{1/2}$ by the $j$-only cancellation -- but the Lamb shift, a quantum-electrodynamic effect at order ${\alpha }^3$, breaks even this degeneracy, raising $2S_{1/2}$ above $2P_{1/2}$ by about $4\times 10^{-6}$ eV.

The chemist's view of this physics differs from the physicist's. A spectroscopist reading a sodium-lamp emission spectrum sees the famous yellow doublet -- the D${}_1$ line at 589.6 nm and the D${}_2$ line at 589.0 nm. The chemist labels these transitions as $3^2{P}_{1/2}\to 3^2{S}_{1/2}$ and $3^2{P}_{3/2}\to 3^2{S}_{1/2}$, using term symbols to encode the fact that the upper state is fine-structure-split. The same term-symbol notation labels every transition in the atomic absorption spectroscopy used in environmental analysis, every flame test in qualitative analysis, and every line in the periodic-table-organised spectra used for elemental identification. Behind that nomenclature lies the physics described above: spin-orbit coupling, the relativistic correction, and the Darwin term, fused by the Dirac equation. Chemistry textbooks treat $j$ as a useful label; physics textbooks derive its origin. Both communities use the same multiplet structure to read spectra.

For transition-metal chemistry, fine structure is small compared with the crystal-field splittings (of order 1 eV) that dominate $d$-orbital energies in coordination complexes. The chemist working on ${\text{Fe}}^{3+}$ or ${\text{Mn}}^{2+}$ rarely needs the explicit spin-orbit operator. For lanthanide and actinide chemistry, however, fine structure becomes essential: the $f$-electrons in Eu${}^{3+}$ or Tb${}^{3+}$ experience strong spin-orbit coupling, and the sharp emission lines of lanthanide phosphors (used in fluorescent lamps and security inks) are transitions between fine-structure-split $f$-manifolds. The colour of europium phosphors is set by ${\alpha }^2{Z}^4$ scaling rather than by the chemical environment. This is one of the few corners of chemistry where the physicist's view -- relativistic atomic structure as the dominant effect -- is the chemist's view.

Hyperfine structure and the role of nuclear spin [Master]

Even after fine structure is accounted for, the spectral lines of hydrogen retain residual splittings of order $10^{-6}$ eV. These come from the interaction of the electron's magnetic moment with the nuclear magnetic moment. The proton, despite its small size, is a spinning charged composite particle and carries an intrinsic magnetic moment $$ \boldsymbol{\mu}_N = g_p \mu_N \mathbf{I}/\hbar, \qquad \mu_N = \frac{e\hbar}{2m_p}, $$ where ${\mu }_N\approx {\mu }_B/1836$ is the nuclear magneton, $g_p\approx 5.586$ is the proton g-factor, and $\mathbf{I}$ is the nuclear spin angular momentum (here $I=1/2$). The electron magnetic moment is proportional to ${\mu }_B=e\hslash /(2m_e)$ -- about two thousand times larger because the proton mass appears in the denominator of ${\mu }_N$. The hyperfine interaction is the magnetic coupling between these two moments, encoded in the operator $A\mathbf{I}\cdot \mathbf{J}$ where $\mathbf{J}$ is the total electronic angular momentum and $A$ is the hyperfine coupling constant.

There are two distinct mechanisms. For $s$-states, the dominant contribution is the Fermi contact term, which depends on the probability density of the electron at the nucleus, $|\psi (0){|}^2$. Only $s$-orbitals have nonzero wave function at the nucleus (as derived in Exercise 9 of unit 12.06.01 pending), so only $s$-state hyperfine couplings carry the contact contribution. For $p$, $d$, and $f$ states, the coupling is the magnetic dipole-dipole interaction between the electron orbital plus spin moments and the nuclear moment, modulated by the angular distribution of the electron. The Fermi contact term for the ground state of hydrogen gives a coupling constant $A=1420.405751$ MHz -- one of the most precisely measured numbers in all of physics.

The two coupled angular momenta $\mathbf{I}=1/2$ and $\mathbf{J}=1/2$ combine into a total $\mathbf{F}=\mathbf{I}+\mathbf{J}$ with $F=0$ (singlet) or $F=1$ (triplet). The hyperfine splitting raises the triplet above the singlet by $\Delta E=hA\approx 5.87\times 10^{-6}$ eV. The corresponding photon wavelength is $\lambda =c/A=21.106$ cm -- the 21-centimetre line of neutral hydrogen, the most important spectral line in radio astronomy. Atomic hydrogen pervades the interstellar medium of the galaxy, and the 21-cm line maps the distribution of cold neutral hydrogen, the structure of spiral arms, and the rotation curves of galaxies that gave the first direct evidence for dark matter. The transition is forbidden (a magnetic dipole transition rather than an electric dipole), so its lifetime is about 10 million years per atom -- but the sheer density of hydrogen in the galaxy gives a measurable signal. Every radio telescope built since the 1950s has the 21-cm line as a primary observing target.

For chemists, the 21-cm line opens a door into astrochemistry: by mapping the spatial distribution and Doppler shifts of neutral hydrogen, radio astronomers chart the gas clouds where stars and planets form. The chemistry of interstellar molecules -- water ice, formaldehyde, methanol, complex organics including amino acids -- happens in environments traced by hyperfine spectroscopy. The hyperfine line is the chemical mapping tool of the galaxy.

A more direct chemical application of hyperfine interactions is nuclear magnetic resonance (NMR). The NMR experiment places nuclei in an external magnetic field $B_0$ and measures the Larmor precession frequency of their magnetic moments, $\omega =\gamma B_0$ where $\gamma =g_N{\mu }_N/\hslash$ is the nuclear gyromagnetic ratio. The physics is the same as the hyperfine interaction in hydrogen, with the electron-nucleus coupling replaced by the external field. The crucial chemical observable is the chemical shift: nuclei in different molecular environments experience slightly different effective magnetic fields because the surrounding electron density partially shields the nucleus from $B_0$. The shielding constant $\sigma$ depends on the electronic structure -- specifically on the diamagnetic and paramagnetic responses of the electron cloud around the nucleus. Both responses are computed from hydrogenic-type orbitals: chemical shifts are derived from the same wave-function machinery as the hyperfine coupling, just with external instead of internal magnetic fields.

The $J$-coupling between adjacent magnetic nuclei in a molecule -- the splitting patterns that make NMR diagnostic -- is also a hyperfine-type interaction. Two nuclei interact through the electron cloud that connects them, with the through-bond coupling determined by the Fermi contact density at each nucleus from the bonding orbitals. A proton-proton coupling of 7 Hz in ethanol is a measurement of the contact density of the C-H bonding orbital at each carbon nucleus, transmitted through the C-C bond. NMR spectroscopy thus rests directly on the hydrogenic Fermi contact physics. A chemistry student running a ${}^1$H NMR experiment in an undergraduate lab is, without knowing it, exploiting the same physics that produces the 21-cm line in radio astronomy.

Beyond hydrogen, hyperfine splittings carry information about nuclear structure. The hyperfine anomaly between ${}^{205}$Tl and ${}^{203}$Tl gave one of the first measurements of nuclear finite-size effects. Hyperfine spectroscopy of singly-ionised atoms in ion traps now provides the most precise frequency standards: the caesium hyperfine transition ${\nu }_{Cs}=9192631770$ Hz defines the SI second, and the hydrogen maser uses the hyperfine transition we just derived as a stable frequency reference. The chemist working with atomic-clock-stabilised spectroscopy in metrology, or with NMR shimming in protein crystallography, is one short conceptual step from the hydrogen hyperfine structure derived from $|\psi (0){|}^2$ for the $1s$ orbital.

A practical chemical consequence appears in electron spin resonance (ESR), the electronic analogue of NMR. In radicals or transition-metal complexes with unpaired electrons, ESR detects the splitting of the electron spin in an external field, with hyperfine structure imprinted by surrounding nuclei. The hydrogen-atom radical ${\text{H}}^{\bullet }$ has a particularly clean ESR doublet with splitting $A=1420$ MHz -- the same hyperfine constant that gives the 21-cm astronomical line. The same number, used in two completely different contexts, by two different scientific communities, is one of the cleanest examples of how the hydrogenic Hamiltonian unifies physics and chemistry.

Breakdown for multi-electron atoms: why hydrogen is special [Master]

In hydrogen, the energy depends only on $n$: $E_{2s}=E_{2p}$, $E_{3s}=E_{3p}=E_{3d}$. In multi-electron atoms this degeneracy is broken. The mechanism is the interplay of penetration and shielding.

Penetration refers to the probability of an electron being found close to the nucleus. For orbitals with the same $n$, those with lower $l$ have greater penetration. The $s$ orbital ($l=0$) has nonzero probability density at the nucleus ($|{\psi }_{ns}(0){|}^2\ne 0$), while $p$, $d$, and $f$ orbitals all vanish at $r=0$. The physical reason is that $l=0$ means zero angular momentum, so there is no centrifugal barrier keeping the electron away from the nucleus.

Shielding refers to the reduction in effective nuclear charge felt by an electron due to other electrons between it and the nucleus. The quantitative measure is $Z_{\text{eff}}=Z-\sigma$, where $\sigma$ is the shielding constant.

The interplay produces the energy ordering. Consider the $n=3$ shell: $3s$ penetrates most and feels the largest $Z_{\text{eff}}$, giving the lowest energy. $3p$ penetrates less, feels a smaller $Z_{\text{eff}}$, and is higher in energy. $3d$ penetrates least and is highest. The same logic explains why $4s<3d$: the $4s$ orbital, despite having $n=4$, has substantial inner loops that penetrate the core, feeling a higher $Z_{\text{eff}}$ than the $3d$ orbital whose electron density is concentrated further from the nucleus.

The quantitative framework is the self-consistent field method (Hartree, 1928; Hartree-Fock, 1930). The orbital energies ${\epsilon }_{nl}$ are computed iteratively: guess initial orbitals, compute the effective potential from all other electrons, solve the one-electron Schrodinger equation in that potential, and iterate to self-consistency. The resulting orbital energies follow the Aufbau ordering for most neutral atoms.

Koptyug's 1982 tabulation of $Z_{\text{eff}}$ values for all neutral atoms (computed from Hartree-Fock orbital energies via ${\epsilon }_{nl}=-Z_{\text{eff}}^2{R}_H/n^{\ast 2}$ where $n^{\ast }$ is an effective quantum number) shows that $Z_{\text{eff}}$ for $4s$ in potassium ($Z=19$) is approximately 2.20, while $Z_{\text{eff}}$ for $3d$ is approximately 3.17 -- yet $4s$ is lower in energy because the effective principal quantum number $n^{\ast }$ for $4s$ ($\approx 2.26$) is much smaller than for $3d$ ($\approx 2.90$), reflecting the greater penetration of $4s$.

Hartree-Fock is not the whole story. The HF method approximates the many-electron wave function as a single Slater determinant of one-electron orbitals; each electron moves in an averaged field of the others. Real electrons are correlated: they avoid one another instantaneously, not just on average. The electron correlation energy -- the difference between the exact ground-state energy and the HF energy -- is small in absolute terms (typically a few percent of the total electronic energy) but it dominates the chemistry of bond formation. The He atom HF energy is $-77.87$ eV; the exact non-relativistic energy is $-79.01$ eV; the difference of $1.14$ eV is the correlation energy. This is more than chemical accuracy ($\sim 0.04$ eV per bond), so chemists cannot use bare HF for quantitative predictions.

Post-Hartree-Fock methods recover correlation through systematic additions. Configuration interaction (CI) writes the wave function as a sum of Slater determinants -- the HF ground state plus single, double, triple, and higher excitations. Full CI is exact within the chosen one-electron basis but scales factorially with the number of electrons. Møller-Plesset perturbation theory (MP2, MP3, MP4) treats correlation as a perturbation to HF; MP2 is the cheapest correlated method and recovers most of the correlation energy for small systems. Coupled-cluster theory (CCSD, CCSD(T)) is the modern gold standard, summing infinite classes of excitations through an exponential ansatz. The "(T)" in CCSD(T) marks the perturbative triple excitations and brings the method to chemical accuracy for many molecules. Density functional theory (DFT) is a different approach altogether: the Hohenberg-Kohn theorems show that the electron density $\rho (\mathbf{r})$, a function of three variables rather than $3N$, suffices to determine the ground state, and Kohn-Sham DFT computes $\rho$ self-consistently in an effective single-particle potential whose exchange-correlation contribution is approximated. DFT is the workhorse of computational chemistry today; the choice of exchange-correlation functional (B3LYP, PBE, M06-2X, etc.) is where chemists encode their experience.

Even with these tools, the labels $s$, $p$, $d$, $f$ from hydrogen do not transfer cleanly to multi-electron atoms. In transition metals, the $3d$ and $4s$ orbital energies are close enough that small chemical perturbations rearrange their order. The Madelung rule ($n+l$ filling) is empirical and fails for chromium ($[Ar]3d^{5}4s^1$ rather than $3d^{4}4s^2$), copper ($[Ar]3d^{10}4s^1$), and several other anomalies. Once a transition metal is incorporated into a complex, the molecular orbitals diagonalise mixed $s$/$p$/$d$ character and the original labels lose meaning. The chemist's claim that an iron atom has $3d^{6}4s^2$ valence electrons is a useful shorthand, not an exact statement -- the real wave function is a CI mixture of many configurations.

For lanthanides, the situation is worse: $4f$, $5d$, $6s$, and $6p$ orbitals overlap energetically, and lanthanide chemistry requires explicit relativistic and correlation treatment. For actinides, the $5f$ orbitals are sufficiently extended that they participate in covalent bonding (uranyl ${\text{UO}}_2^{2+}$ has $5f$-character in the U-O bonds), breaking the textbook story that $f$-electrons are core-like. These are not exotic anomalies -- they are the chemistry of fluorescent phosphors, of nuclear fuel cycles, of magnetic permanent magnets (neodymium-iron-boron) and of medical imaging contrast agents (gadolinium).

The deepest breakdown of the hydrogenic picture occurs in strongly correlated systems: high-temperature superconductors, Mott insulators, heavy-fermion compounds, and the metal-insulator transitions in transition-metal oxides. In these systems, single-particle pictures fail entirely; the relevant low-energy degrees of freedom are collective excitations (Cooper pairs, magnons, Kondo singlets) that have no description in terms of independent hydrogenic orbitals. Yet even here the language inherited from hydrogen -- $d$-band, $f$-electron, crystal-field-split -- remains the chemist's starting vocabulary. The hydrogenic orbitals are the basis we expand around, and the chemistry of correlated systems is the description of how that basis fails.

The takeaway: hydrogen is special because it has exactly one electron, so its solution is exact, its degeneracies are protected by SO(4) symmetry, and its orbital labels are unambiguous. Every other atom is a perturbation away from this ideal, and the size of the perturbation determines how much of hydrogen's clean structure survives. The art of modern computational chemistry is knowing which hydrogenic features to keep, which to modify, and which to discard. The next unit on multi-electron atomic structure quantifies these modifications systematically.

Cross-domain bridge: connection to the physics-side hydrogen unit [Master]

This unit and unit 12.06.01 pending cover the same physical system from different starting points. The same Coulomb Hamiltonian governs both treatments; the same wave functions ${\psi }_{nlm}$ describe the same electron-proton bound states; the same energy levels $E_n=-R_H/n^2$ emerge from the same separation of variables. Yet the two units exist in parallel, neither delegating to the other, because the framing differs and the framing determines what is worth proving, what is worth labelling, and what counts as the natural language of the subject.

The physics unit 12.06.01 frames hydrogen as a bound-state problem in a central force. The starting question is: given a $1/r$ potential, what is the spectrum of the corresponding Hamiltonian as a self-adjoint operator on $L^2({\mathbb{R}}^3)$? The answer involves separation of variables, the radial equation, the associated Laguerre equation, and the SO(4) symmetry that explains the accidental $l$-degeneracy. The natural objects in this framing are operators, commutators, and spectral decompositions. The Pauli algebraic solution using the Runge-Lenz vector is the conceptual centrepiece because it makes the symmetry explicit. The next steps from the physics unit are perturbation theory (Stark effect, fine structure, hyperfine splitting as perturbations of the unperturbed Hamiltonian) and the algebraic structure of more general central-force problems.

The chemistry unit 14.04.01 -- the one you are reading -- frames hydrogen as the template for the periodic table. The starting question is: given that hydrogen has the only exactly solvable atom, how do its orbital shapes and quantum-number labels organise the chemistry of every other element? The answer involves orbital shapes ($s$, $p$, $d$, $f$ visualised as boundary surfaces), radial distribution functions (where electrons live in space), node-counting rules (what makes $4s$ penetrate more than $3d$), and the Aufbau principle (how to fill orbitals). The natural objects in this framing are orbitals as spatial distributions, electron configurations as filling rules, and trends across the periodic table. The next steps from the chemistry unit are the construction of molecular orbitals as linear combinations of atomic orbitals, the explanation of ionisation energies and electronegativities from orbital structure, and the empirical periodicities that gave chemistry its organising principles before quantum mechanics existed.

Concretely, where the units overlap and where they diverge:

• Schrödinger equation and its solution. Both units derive the same energy levels and wave functions. 12.06.01 carries the proof through to completion (including the SO(4) symmetry analysis); 14.04.01 sketches the proof and refers to 12.06.01 for the details. The chemistry version emphasises the form of the radial functions $R_{nl}(r)$ and the radial distribution function $r^2{R}_{nl}^2$, because these are the quantities chemists need for visualising orbital shapes and computing electron probability densities. The physics version emphasises the operator structure and the spectral theorem.

• Quantum numbers. Both units use $(n,l,m_l)$. The physics unit emphasises that $l$ labels representations of SO(3) and $n$ labels representations of SO(4); the chemistry unit emphasises that $l$ encodes orbital shape ($s$, $p$, $d$, $f$) and $n$ encodes shell. The Pauli exclusion principle and electron spin $m_s$ appear in both but play different roles: in chemistry they generate the Aufbau filling rule, while in physics they connect to the spin-statistics theorem.

• Fine and hyperfine structure. This unit derives the qualitative chemical consequences (term symbols, the sodium D-line, the 21-cm line, the connection to NMR). Unit 12.06.01 derives the quantitative perturbation-theoretic structure (Dirac equation, spin-orbit operator, hyperfine Hamiltonian). The Lamb shift is mentioned in both but proved fully in neither; quantum electrodynamics lies past both units.

• Multi-electron atoms. The chemistry unit treats Hartree-Fock and electron correlation in the Master tier; the physics unit treats multi-electron atoms as a brief application of perturbation theory and refers forward to many-body physics. The chemistry framing here is essential -- the configuration $[Ar]3d^{6}4s^2$ is a chemical label that has no analogue in the operator-theoretic physics formulation.

• Spectroscopy. Both units mention the Rydberg formula and the Balmer/Lyman/Paschen series. The chemistry unit organises spectroscopy by chemical application (atomic absorption, emission, ESR, NMR); the physics unit organises spectroscopy by theoretical apparatus (transition rates, selection rules, gauge invariance of the dipole approximation).

Why both units? Because a chemistry student approaching the periodic table needs the orbital-shape framing to make sense of why argon is inert and chlorine is reactive, while a physics student approaching quantum mechanics needs the operator framing to make sense of why the spectrum is discrete and the bound-state subspace is complete. Neither framing is reducible to the other in any useful pedagogical sense. The chemistry student who knows only the orbital shapes will not be able to derive the SO(4) symmetry or generalise to non-Coulomb central potentials. The physics student who knows only the operator structure will not be able to predict ionic radii or explain the colour of transition-metal complexes. The two framings are complementary, and a graduate student in either field needs the other.

The hooks_out chain from this unit reflects the chemistry framing: 14.02.01 (Lewis structures, using valence orbitals), 14.05.01 (molecular orbital theory, building MOs from atomic orbitals), 16.01.01 (periodic trends, organised by orbital energetics). The hooks_out chain from 12.06.01 reflects the physics framing: 12.07.01 (perturbation theory, treating fine/hyperfine structure), 14.12.01 (UV-Vis spectroscopy, applying the dipole transition rate to atomic transitions). The two hook chains intersect at spectroscopy and atomic structure, providing natural cross-links for students and researchers moving between the two communities.

The chemistry framing has one feature absent from the physics treatment: it explicitly mentions the shapes of orbitals as visualisation aids. The cartoon $p_x$ dumbbell, the $d_{xy}$ cloverleaf, the $d_{z^2}$ donut-and-lobe -- these are mnemonic devices that help chemists predict bonding geometry, but they are not in any sense the only or even the most natural representation of the wave function. A physicist would more commonly use the complex spherical harmonics $Y_l^m$, with $p_x$ and $p_y$ appearing only as the real linear combinations $(Y_1^{-1}-Y_1^1)/\sqrt{2}$ and $(Y_1^{-1}+Y_1^1)/(i\sqrt{2})$. The chemist's choice of real orbitals is dictated by chemistry: $p_x$ points along the $x$-axis, and that is the orientation in which a $\sigma$-bond forms. The physicist's choice of complex orbitals is dictated by symmetry: $Y_1^m$ are eigenfunctions of ${\hat{L}}_z$, and that is the operator whose conservation matters. Reading both units side by side lets a student see the same wave function from two complementary angles.

Connections [Master]

• Physics hydrogen atom 12.06.01 pending provides the full mathematical solution to the Schrodinger equation for the hydrogen atom: the radial functions, spherical harmonics, and energy eigenvalues. That unit develops the SO(4) symmetry, the Pauli algebraic solution, and the spectral structure of the Coulomb Hamiltonian as a self-adjoint operator. The present unit adapts that treatment to the chemical context, emphasising orbital shapes, quantum numbers, term symbols, and the connection to the periodic table. The two units are complementary rather than redundant -- see the cross-domain bridge sub-section for an explicit comparison. A student doing computational chemistry uses both: the physics formalism to understand which approximations are valid, and the chemistry framing to predict bonding and reactivity.

• Atomic structure and electron configurations 14.01.01 uses the orbitals defined here to build up electron configurations via the Aufbau principle. The quantum numbers $(n,l,m_l,m_s)$ that label orbitals are the same ones that fill according to Aufbau ordering.

• Lewis structures 14.02.01 simplify orbital descriptions to dot diagrams of valence electrons. The valence orbitals identified here (the outermost $s$ and $p$ orbitals) are the ones that participate in Lewis bonding pictures.

• Molecular orbital theory 14.05.01 pending builds molecular orbitals as linear combinations of the atomic orbitals whose shapes and energies are established here. The symmetry of atomic orbitals ($s$ is spherical, $p$ has directional lobes) determines which atomic orbitals can overlap to form molecular orbitals.

• Periodic trends 16.01.01 pending -- ionisation energy, atomic radius, electronegativity -- are direct consequences of the orbital energies and radial distributions quantified here. The increase in ionisation energy across a period reflects increasing $Z_{\text{eff}}$ in orbitals of the same $n$.

Historical & philosophical context [Master]

The Bohr model (1913) introduced quantised orbits for the hydrogen atom and reproduced the Balmer series formula, but could not account for the shapes of orbitals or the fine structure of spectral lines. Schrodinger's wave equation (1926) provided the full solution, introducing the concept of the wave function as a probability amplitude. The interpretation of $|\psi {|}^2$ as a probability density was proposed by Max Born (1926), fundamentally changing the ontological status of the electron from a particle with a definite position to a delocalised probability distribution.

The term "orbital" was coined by Robert Mulliken in 1932 to distinguish the wave-mechanical description from the Bohr "orbit." Mulliken and Friedrich Hund independently developed molecular orbital theory in the late 1920s, building directly on the hydrogenic orbital shapes.

The shapes of orbitals -- spherical for $s$, dumbbell for $p$, cloverleaf for $d$ -- are boundary-surface representations of the probability density. These images, ubiquitous in chemistry textbooks, are somewhat misleading: the true probability density extends to infinity (decaying exponentially) and has no sharp boundary. The "surface" is typically drawn at the radius enclosing 90% of the probability, but this is a convention, not a physical boundary. The philosophical point: orbitals are mathematical constructs that make accurate predictions, but they are not directly observable. Quantum mechanics predicts only the results of measurements, and no measurement can directly image an orbital. Scanning tunnelling microscopy (STM) images of electron density on surfaces come closest, but they measure a related quantity (local density of states) rather than $|\psi {|}^2$ directly.

A deeper philosophical question concerns the status of orbitals in a many-electron atom. In hydrogen, the orbital is uniquely defined because there is exactly one electron and its wave function is a one-particle quantity. In any other atom, the true wave function is an antisymmetrised many-body object $\Psi ({\mathbf{r}}_1,{\mathbf{r}}_2,\dots ,{\mathbf{r}}_N)$ depending on $3N$ spatial coordinates plus spin, and there is no rigorous decomposition into individual orbitals. The Hartree-Fock orbitals are the best single-determinant approximation but are not unique -- any unitary mixing of occupied orbitals leaves the determinant unchanged, so the labels "$2s$", "$2p_x$", and so on become matters of convention rather than physical reality. The natural orbitals derived from the exact density matrix are unique but their occupation numbers are non-integer, reflecting electron correlation. The chemist's habit of speaking about "the $4s$ electron" in potassium is a useful shorthand, but the underlying quantum-mechanical object is the full antisymmetrised wave function, not a one-electron orbital.

This ambiguity surfaces in chemistry through the question of orbital "reality". When ARPES (angle-resolved photoemission spectroscopy) maps Fermi-surface states in transition-metal compounds, what is being measured is the spectral function $A(\mathbf{k},\omega )$ -- a many-body quantity that reduces to single-particle orbitals only in the non-interacting limit. When density functional theory gives Kohn-Sham orbitals, the orbitals themselves have no rigorous physical meaning beyond reproducing the correct ground-state density; their orbital energies are not exact ionisation potentials. The 2018 Nature paper claiming "observation of single hydrogen molecular orbitals" by STM, and the controversy it sparked, illustrates how blurred the line is between mathematical artefact and physical observable. The chemist must be aware that the standard orbital picture is a representation, not a microscope image.

For the philosophy of science, the hydrogen atom is the most thoroughly studied physical system in human history. The first 18 digits of the Lyman-alpha wavelength have been measured by frequency-comb laser spectroscopy, and they agree with QED predictions to a precision unmatched in any other domain. The proton-radius puzzle -- a 4-sigma discrepancy in proton charge radius measurements between electronic hydrogen and muonic hydrogen spectroscopy, resolved in 2019 by improved measurements -- shows that hydrogen continues to test fundamental physics at the boundaries. The same atom that organises the chemist's periodic table is also the most sensitive probe we have for new physics beyond the Standard Model. Quantum chemistry inherits this precision, and the chemistry built on hydrogenic orbitals is, at its foundation, in remarkable agreement with experiment.

Bibliography [Master]

• Schrodinger, E., "Quantisierung als Eigenwertproblem", Ann. Phys. 384 (1926), 361-376. The hydrogen atom solution.
• Born, M., "Zur Quantenmechanik der Stossvorgange", Z. Physik 37 (1926), 863-867. The probabilistic interpretation.
• Mulliken, R. S., "Electronic Structures of Polyatomic Molecules and Valence", Phys. Rev. 40 (1932), 55-62. The term "orbital."
• McQuarrie, D. A., Quantum Chemistry, 2e (University Science Books, 2008), Ch. 6. The hydrogen atom.
• Levine, I. N., Quantum Chemistry, 7e (Pearson, 2014), Ch. 6-7. Hydrogen and many-electron atoms.
• Tro, N. J., Chemistry: A Molecular Approach, 6e (Pearson, 2023), Ch. 7. Introductory treatment.
• Atkins, P. & Friedman, R., Molecular Quantum Mechanics, 5e (Oxford, 2011). Advanced reference for orbital theory.
• Pyykkö, P., "Relativistic Effects in Structural Chemistry", Chem. Rev. 88 (1988), 563-594. The standard chemistry-oriented review of how relativistic corrections shape periodic trends and reactivity in heavy elements.
• Schwerdtfeger, P. & Nagle, J. K., "2018 Table of Static Dipole Polarizabilities of the Neutral Elements", Mol. Phys. 117 (2019), 1200-1225. Reference table for atomic polarisabilities computed from accurate many-electron wave functions on the hydrogenic-orbital basis.
• Helgaker, T., Jørgensen, P. & Olsen, J., Molecular Electronic-Structure Theory (Wiley, 2000). The authoritative reference for Hartree-Fock, coupled cluster, and configuration interaction starting from one-electron orbitals.
• Drake, G. W. F., ed., Springer Handbook of Atomic, Molecular, and Optical Physics (Springer, 2006). Comprehensive reference covering hyperfine structure, fine structure, and QED corrections from hydrogen to multi-electron systems.

---

Wave 3 chemistry unit. All hooks_out targets are proposed; no receiving unit yet exists to confirm them. Status remains draft pending Tyler's review and external chemistry reviewer per CHEMISTRY_PLAN §6.