Unit 4: Hydrogen Atom, Central Force Dynamics & Hydrogenic Orbitals

The hydrogen atom consists of a nucleus of charge \(+Z e\) (mass \(M\)) and an electron of charge \(-e\) (mass \(m_e\)) interacting via the attractive Coulomb electrostatic potential.

Separation of Center-of-Mass and Relative Coordinates

The total classical and quantum Hamiltonian for the two-body system is:

By introducing the center-of-mass vector \(\mathbf{R} = \frac{M \mathbf{r}_N + m_e \mathbf{r}_e}{M + m_e}\) with total mass \(M_{\text{tot}} = M + m_e\), and the relative position vector \(\mathbf{r} = \mathbf{r}_e - \mathbf{r}_N\) with reduced mass:

the kinetic energy operator factors exactly into:

where the Coulomb potential \(V(r) = -\frac{Z e^2}{4\pi \varepsilon_0 r}\) depends strictly on the scalar radial distance \(r = |\mathbf{r}|\), establishing central force symmetry. The center-of-mass motion describes a free particle of mass \(M_{\text{tot}}\), while the internal electronic dynamics are governed by the relative Hamiltonian:

Spherical Polar Coordinate Representation

In spherical polar coordinates \((r, \theta, \phi)\), where \(x = r \sin\theta \cos\phi\), \(y = r \sin\theta \sin\phi\), \(z = r \cos\theta\), the Laplacian operator transforms to:

Notice that the angular derivative operators are precisely related to the square of the orbital angular momentum operator \(\hat{L}^2\):

Consequently, the relative Hamiltonian is written compactly as:

Because \(\hat{H}\), \(\hat{L}^2\), and \(\hat{L}_z\) mutually commute:

they possess a complete simultaneous orthonormal eigenbasis. The full stationary wavefunctions \(\psi(r, \theta, \phi)\) separate into a radial function and angular Spherical Harmonics:

Substituting this separable product and utilizing \(\hat{L}^2 Y_l^m(\theta, \phi) = \hbar^2 l(l + 1) Y_l^m(\theta, \phi)\), the angular dependence cancels identically, yielding the Radial Schrödinger Equation:

The term \(\frac{\hbar^2 l(l+1)}{2\mu r^2}\) represents a repulsive centrifugal barrier arising from orbital angular momentum.

To solve the radial equation analytically for bound states (\(E < 0\)), we define the reduced radial function \(u(r) = r R(r)\).

Dimensionless Radial Equation

Substituting \(u(r)\) converts the equation into:

We introduce the parameter \(\kappa = \sqrt{-2\mu E}/\hbar\) and the dimensionless variable:

In terms of \(\rho\), the differential equation becomes:

where \(\lambda = \frac{Z e^2}{4\pi \varepsilon_0 \hbar} \sqrt{\frac{\mu}{-2 E}}\).

Asymptotic Analysis

1. As \(\rho \rightarrow \infty\): The equation reduces to \(\frac{d^2 u}{d\rho^2} - \frac{1}{4} u = 0\), whose physically well-behaved normalizable solution is \(u(\rho) \sim e^{-\rho / 2}\).

2. As \(\rho \rightarrow 0\): The dominant term is \(\frac{d^2 u}{d\rho^2} - \frac{l(l+1)}{\rho^2} u = 0\). The ansatz \(u(\rho) \sim \rho^s\) requires \(s(s - 1) = l(l + 1)\), which has roots \(s = l + 1\) and \(s = -l\). Since \(u(0)\) must vanish for \(R(0)\) to remain finite, we retain \(s = l + 1\). Thus, \(u(\rho) \sim \rho^{l+1}\), which implies \(R(\rho) \sim \rho^l\).

Series Expansion & Associated Laguerre Polynomials

Factoring out the asymptotic behavior, we set:

Substituting into the radial equation yields Kummer's confluent hypergeometric equation for \(v(\rho)\):

For the power series expansion \(v(\rho) = \sum_{k=0}^\infty a_k \rho^k\) to terminate into a polynomial of degree \(n_r\) (preventing \(v(\rho)\) from diverging as \(e^\rho\) at large distances), the numerator in the recurrence relation must vanish:

Defining the principal quantum number \(n = n_r + l + 1\), we have \(\lambda = n\), where \(n \in \{1, 2, 3, \dots\}\) and \(l \in \{0, 1, \dots, n-1\}\). The polynomial solutions are proportional to the Associated Laguerre Polynomials \(L_{n-l-1}^{2l+1}(\rho)\):

Normalized Hydrogenic Radial Wavefunctions

The fully normalized radial eigenfunction is:

where \(a_\mu = \frac{4\pi \varepsilon_0 \hbar^2}{\mu e^2} \approx a_0 = 0.529177 \text{ \AA}\) is the modified Bohr radius and \(\rho = \frac{2 Z r}{n a_\mu}\). Explicit formulas for the lowest states:

• 1s (\(n=1, l=0\)): \(R_{10}(r) = 2 \left(\frac{Z}{a_0}\right)^{3/2} e^{-Z r / a_0}\)
• 2s (\(n=2, l=0\)): \(R_{20}(r) = \frac{1}{\sqrt{2}} \left(\frac{Z}{a_0}\right)^{3/2} \left(1 - \frac{Z r}{2 a_0}\right) e^{-Z r / 2 a_0}\)
• 2p (\(n=2, l=1\)): \(R_{21}(r) = \frac{1}{2\sqrt{6}} \left(\frac{Z}{a_0}\right)^{3/2} \left(\frac{Z r}{a_0}\right) e^{-Z r / 2 a_0}\)
• 3s (\(n=3, l=0\)): \(R_{30}(r) = \frac{2}{81\sqrt{3}} \left(\frac{Z}{a_0}\right)^{3/2} \left(27 - 18\frac{Zr}{a_0} + 2\frac{Z^2 r^2}{a_0^2}\right) e^{-Zr/3a_0}\)

From the polynomial termination condition \(\lambda = n\), the bound-state energy eigenvalues are obtained directly:

where \(R_y = \frac{\mu e^4}{32 \pi^2 \varepsilon_0^2 \hbar^2} \approx 13.60569\text{ eV} = 1\text{ Ry} = \frac{1}{2} E_h\) (where \(E_h = 27.2114\text{ eV}\) is 1 Hartree).

Energy Levels and Quantum Numbers

The electronic energy depends strictly on the principal quantum number \(n\) and is independent of \(l\) and \(m\). This degeneracy is remarkable:

1. Principal quantum number \(n\): \(n \in \{1, 2, 3, \dots\}\) governs overall orbital scale and energy.

2. Azimuthal (orbital) quantum number \(l\): \(l \in \{0, 1, 2, \dots, n - 1\}\) designates orbital angular momentum magnitude \(|\mathbf{L}| = \hbar \sqrt{l(l+1)}\).

3. Magnetic quantum number \(m\): \(m \in \{-l, -l+1, \dots, +l\}\) governs the spatial projection \(L_z = m\hbar\).

Degeneracy Analysis

For a given \(n\), the orbital angular momentum can take \(n\) values (\(l = 0, 1, \dots, n-1\)). For each \(l\), there are \(2l + 1\) distinct \(m\) projections. The total spatial orbital degeneracy \(g_n\) is:

Including electron spin degeneracy (\(m_s = \pm 1/2\)), the total state degeneracy is \(2 n^2\):

• \(n = 1\): \(1^2 = 1\) orbital (\(1s\)) \(\rightarrow\) 2 quantum states
• \(n = 2\): \(2^2 = 4\) orbitals (\(2s, 2p_x, 2p_y, 2p_z\)) \(\rightarrow\) 8 quantum states
• \(n = 3\): \(3^2 = 9\) orbitals (\(3s, 3p (3), 3d (5)\)) \(\rightarrow\) 18 quantum states

The accidental degeneracy across different \(l\) values for a given \(n\) is a consequence of the higher dynamical symmetry group \(SO(4)\) generated by the conservation of the Laplace-Runge-Lenz vector in a pure \(1/r\) Coulomb potential.

The Rydberg Formula for Optical Transitions

A radiative transition between an upper state \(n_2\) and a lower state \(n_1\) involves photon emission or absorption with photon energy:

The transition wavenumber \(\tilde{\nu} = 1/\lambda\) is given by the general Rydberg formula:

where \(R_\infty = \frac{m_e e^4}{8 \varepsilon_0^2 h^3 c} \approx 109737.316\text{ cm}^{-1}\), and \(R_H = R_\infty \frac{M_p}{M_p + m_e} \approx 109677.583\text{ cm}^{-1}\) for atomic hydrogen. The classic spectral series of hydrogen:

• Lyman Series (\(n_1 = 1, n_2 \ge 2\)): Ultraviolet range (\(121.6\text{ nm} \rightarrow 91.2\text{ nm}\))
• Balmer Series (\(n_1 = 2, n_2 \ge 3\)): Visible range (\(656.3\text{ nm} \text{ [H}\alpha\text{]} \rightarrow 364.6\text{ nm}\))
• Paschen Series (\(n_1 = 3, n_2 \ge 4\)): Near-infrared range (\(1875\text{ nm} \rightarrow 820.4\text{ nm}\))
• Brackett Series (\(n_1 = 4\)) and Pfund Series (\(n_1 = 5\)): Mid-infrared range.

The physical probability of locating the electron inside an infinitesimal volume element \(d\tau = r^2 \sin\theta \, dr \, d\theta \, d\phi\) is:

The Radial Distribution Function \(P(r)\)

To find the probability of finding the electron at a radial distance between \(r\) and \(r + dr\) regardless of angle, we integrate over all angles \(\theta \in [0, \pi]\) and \(\phi \in [0, 2\pi]\):

Because the Spherical Harmonics are normalized (\(\int |Y_l^m|^2 d\Omega = 1\)), we obtain the Radial Probability Density:

Notice the crucial difference:

• The probability density per unit volume at the nucleus (\(r = 0\)) is \([R_{n0}(0)]^2 > 0\) for all \(s\)-orbitals (\(l = 0\)).
• The radial probability density \(P(r)\) vanishes at the nucleus: \(P_{n l}(0) = 0^2 [R_{n l}(0)]^2 = 0\) for all states because the spherical shell volume \(4\pi r^2 dr\) shrinks to zero as \(r \rightarrow 0\).

Most Probable Radius vs Mean Radius \(\langle r \rangle\)

For the ground state of hydrogen (\(1s, Z = 1\)):

To find the most probable radius \(r_{\text{mp}}\), we maximize \(P_{10}(r)\):

The most probable distance of the electron from the proton in the ground state matches the historical first Bohr radius! In contrast, the expectation value (average distance) \(\langle r \rangle\) is:

Because the exponential tail extends outward, \(\langle r \rangle > r_{\text{mp}}\). In general:

Node Topology

The total number of nodes in an atomic orbital wavefunction is \(n - 1\):

1. Radial nodes: Points where \(R_{n l}(r) = 0\) (excluding \(r = 0\) and \(r = \infty\)). The number of radial nodes is:

2. Angular nodal surfaces: Conical or planar surfaces where \(Y_l^m(\theta, \phi) = 0\). The number of angular nodes is:

3. Total nodes:

Examples:

• \(1s\): \(n=1, l=0 \implies 0\) radial nodes, 0 angular nodes.
• \(2s\): \(n=2, l=0 \implies 1\) radial node (\(r = 2 a_0 / Z\)), 0 angular nodes.
• \(2p\): \(n=2, l=1 \implies 0\) radial nodes, 1 angular nodal plane (e.g., \(xy\)-plane for \(p_z\)).
• \(3s\): \(n=3, l=0 \implies 2\) radial nodes, 0 angular nodes.
• \(3d\): \(n=3, l=2 \implies 0\) radial nodes, 2 angular nodal surfaces.

The angular parts of hydrogenic wavefunctions are the Spherical Harmonics \(Y_l^m(\theta, \phi)\), normalized over the unit sphere:

with \(Y_l^{-m}(\theta, \phi) = (-1)^m [Y_l^m(\theta, \phi)]^*\).

Real Orbitals in Chemistry

In physical chemistry, complex orbitals with \(m = \pm 1, \pm 2\) are generally replaced by real linear combinations that point along Cartesian coordinate axes. Since the Hamiltonian is degenerate in \(m\), any linear combination of degenerate eigenfunctions is also an exact eigenfunction.

\(p\)-Orbitals (\(l = 1\)):

• \(p_z = Y_1^0 = \sqrt{\frac{3}{4\pi}} \cos\theta = \sqrt{\frac{3}{4\pi}} \frac{z}{r}\)
• \(p_x = \frac{1}{\sqrt{2}} (-Y_1^1 + Y_1^{-1}) = \sqrt{\frac{3}{4\pi}} \sin\theta \cos\phi = \sqrt{\frac{3}{4\pi}} \frac{x}{r}\)
• \(p_y = \frac{i}{\sqrt{2}} (Y_1^1 + Y_1^{-1}) = \sqrt{\frac{3}{4\pi}} \sin\theta \sin\phi = \sqrt{\frac{3}{4\pi}} \frac{y}{r}\)

Each \(p\)-orbital consists of two lobes with opposite signs (phases) separated by an angular nodal plane through the nucleus:

• \(p_z\) has nodal plane \(z = 0\) (\(xy\)-plane).
• \(p_x\) has nodal plane \(x = 0\) (\(yz\)-plane).
• \(p_y\) has nodal plane \(y = 0\) (\(xz\)-plane).

\(d\)-Orbitals (\(l = 2\)):

There are five real \(d\)-orbitals formed from \(l = 2\) Spherical Harmonics:

• \(d_{z^2} = Y_2^0 = \sqrt{\frac{5}{16\pi}} (3\cos^2\theta - 1) = \sqrt{\frac{5}{16\pi}} \frac{3z^2 - r^2}{r^2}\) (two conical nodes at \(\cos\theta = \pm 1/\sqrt{3} \approx 54.74^\circ\))
• \(d_{xz} = \frac{1}{\sqrt{2}} (-Y_2^1 + Y_2^{-1}) = \sqrt{\frac{15}{4\pi}} \sin\theta \cos\theta \cos\phi = \sqrt{\frac{15}{4\pi}} \frac{xz}{r^2}\)
• \(d_{yz} = \frac{i}{\sqrt{2}} (Y_2^1 + Y_2^{-1}) = \sqrt{\frac{15}{4\pi}} \sin\theta \cos\theta \sin\phi = \sqrt{\frac{15}{4\pi}} \frac{yz}{r^2}\)
• \(d_{xy} = \frac{i}{\sqrt{2}} (-Y_2^2 + Y_2^{-2}) = \sqrt{\frac{15}{4\pi}} \sin^2\theta \sin\phi \cos\phi = \sqrt{\frac{15}{4\pi}} \frac{xy}{r^2}\)
• \(d_{x^2-y^2} = \frac{1}{\sqrt{2}} (Y_2^2 + Y_2^{-2}) = \sqrt{\frac{15}{16\pi}} \sin^2\theta \cos(2\phi) = \sqrt{\frac{15}{16\pi}} \frac{x^2 - y^2}{r^2}\)

Orbital Boundary Surfaces

In molecular chemistry, orbital shapes are conventionally visualized as 3D isosurfaces encompassing a specified probability (typically 90% or 95% of total electron density):

The spatial directional orientation of \(p\) and \(d\) orbitals governs chemical valence, hybridization (\(sp, sp^2, sp^3, d^2sp^3\)), crystal field splitting in transition metal complexes, and stereochemistry.

In high-resolution spectroscopy, hydrogenic spectral lines reveal small splittings called fine structure, with energy differences of order \(\alpha^2 E_n \sim 10^{-4}\text{ eV}\), where \(\alpha = \frac{e^2}{4\pi\varepsilon_0 \hbar c} \approx \frac{1}{137.036}\) is the fine structure constant.

The Spin-Orbit Interaction Hamiltonian

From the rest frame of the orbiting electron, the positively charged nucleus circulates with velocity \(-\mathbf{v}\), creating an effective internal magnetic field \(\mathbf{B}_{\text{int}}\):

The electron possesses an intrinsic magnetic dipole moment \(\boldsymbol{\mu}_s = -g_s \frac{e}{2 m_e} \mathbf{S} \approx -\frac{e}{m_e} \mathbf{S}\). Incorporating the Thomas precession factor of \(1/2\) (due to the accelerating non-inertial reference frame of the electron), the spin-orbit Hamiltonian is:

Total Angular Momentum Coupling

Define total angular momentum \(\hat{\mathbf{J}} = \hat{\mathbf{L}} + \hat{\mathbf{S}}\). Squaring both sides:

The eigenvalues of \(\hat{\mathbf{L}} \cdot \hat{\mathbf{S}}\) in the coupled representation \(|j, m_j, l, s\rangle\) are:

where for a single electron \(s = 1/2\), so \(j = l + 1/2\) or \(j = l - 1/2\) (for \(l > 0\)).

Relativistic Corrections and Total Fine Structure

The full relativistic correction to order \(\alpha^2\) comprises three terms:

1. Relativistic kinetic energy correction: \(\hat{H}_{\text{rel}} = -\frac{\hat{p}^4}{8 m_e^3 c^2}\)

2. Spin-orbit coupling: \(\hat{H}_{\text{SO}}\)

3. Darwin term (for \(s\)-states, \(l=0\)): \(\hat{H}_D = \frac{\pi \hbar^2 Z e^2}{2 m_e^2 c^2 (4\pi \varepsilon_0)} \delta^3(\mathbf{r})\)

Remarkably, Paul Dirac's relativistic wave equation yields the exact combined energy formula:

States with the same principal quantum number \(n\) and the same total angular momentum \(j\) have identical energies in Dirac theory:

• The \(2s_{1/2}\) (\(n=2, l=0, j=1/2\)) and \(2p_{1/2}\) (\(n=2, l=1, j=1/2\)) levels are strictly degenerate in the Dirac equation.
• In 1947, Willis Lamb and Robert Retherford discovered that \(2s_{1/2}\) lies approximately \(1057.8\text{ MHz}\) above \(2p_{1/2}\). This Lamb shift arises from quantum electrodynamic (QED) vacuum fluctuations of the electromagnetic field and electron self-energy.

When an atom is placed in an external uniform magnetic field \(\mathbf{B} = B \hat{\mathbf{z}}\), its spectral lines split into closely spaced components.

Magnetic Interaction Hamiltonian

The total magnetic dipole moment of an atomic electron is:

where \(\mu_B = \frac{e \hbar}{2 m_e} \approx 9.274 \times 10^{-24}\text{ J}\cdot\text{T}^{-1}\) is the Bohr magneton, and \(g_e \approx 2.002319\) is the electron spin gyromagnetic \(g\)-factor (taken as \(2\) to first order). The Zeeman Hamiltonian is:

Normal Zeeman Effect (Zero Spin, \(S = 0\))

In singlet states where total spin \(S = 0\), \(\hat{\mathbf{J}} = \hat{\mathbf{L}}\), and the interaction Hamiltonian reduces to:

Under optical dipole selection rules \(\Delta m_l = 0, \pm 1\), any spectral transition of frequency \(\nu_0\) splits into exactly three equally spaced lines:

The central line (\(\Delta m = 0\), \(\pi\)-component) is linearly polarized parallel to \(\mathbf{B}\), while the side lines (\(\Delta m = \pm 1\), \(\sigma\)-components) are circularly polarized perpendicular to \(\mathbf{B}\).

Anomalous Zeeman Effect (Non-Zero Spin, Weak Field)

When spin is non-zero and the external magnetic field is weak compared to internal spin-orbit coupling (\(B \ll B_{\text{int}} \sim 10\text{ T}\)), \(j\) and \(m_j\) remain good quantum numbers. In first-order perturbation theory, the Zeeman energy shift is:

where \(g_J\) is the Landé \(g\)-factor:

Values of \(g_J\):

• For pure orbital state (\(s=0\)): \(g_J = 1\)
• For pure spin state (\(l=0, j=1/2, s=1/2\)): \(g_J = 2\)
• For \(^2P_{1/2}\) (\(l=1, s=1/2, j=1/2\)): \(g_J = 1 + \frac{3/4 + 3/4 - 2}{2 \times 3/4} = 1 + \frac{-1/2}{3/2} = 1 - \frac{1}{3} = \frac{2}{3}\)
• For \(^2P_{3/2}\) (\(l=1, s=1/2, j=3/2\)): \(g_J = 1 + \frac{15/4 + 3/4 - 2}{2 \times 15/4} = 1 + \frac{5/2}{15/2} = 1 + \frac{1}{3} = \frac{4}{3}\)

Because different states possess different \(g_J\) factors, transitions split into complex multiplet patterns with more than 3 lines (the anomalous Zeeman effect).

Strong-Field Paschen-Back Limit

When the magnetic field is very strong (\(B \gg B_{\text{int}}\)), the external field decouples \(\mathbf{L}\) and \(\mathbf{S}\). Both \(\hat{L}_z\) and \(\hat{S}_z\) become independently conserved:

With electric dipole selection rules \(\Delta m_s = 0\) and \(\Delta m_l = 0, \pm 1\), the spectrum collapses back into a triplet resembling the normal Zeeman effect.

When an atom is exposed to an external uniform electrostatic field \(\mathbf{E}_{\text{ext}} = \mathcal{E} \hat{\mathbf{z}}\), its spectral lines shift and split—the Stark Effect.

The Stark Interaction Hamiltonian

The electric dipole moment of an atomic electron is \(\mathbf{d} = -e \mathbf{r}\). The electrostatic perturbation Hamiltonian is:

Quadratic Stark Effect (Non-Degenerate Ground State)

For the ground state of hydrogen (\(1s\), \(n = 1\)), the first-order energy correction is:

This vanishes identically because \(|\psi_{1s}|^2\) is symmetric under spatial inversion (\(z \rightarrow -z\)) while \(z\) is antisymmetric (Laporte parity selection rule). The lowest non-vanishing shift occurs at second order:

where \(\alpha_{\text{pol}}\) is the static atomic dipole polarizability. For ground-state hydrogen:

The energy shift is strictly proportional to \(\mathcal{E}^2\) (Quadratic Stark Effect).

Linear Stark Effect (Degenerate Excited States)

For the first excited state (\(n = 2\)), there are 4 degenerate orbitals: \(2s, 2p_0, 2p_{+1}, 2p_{-1}\). Because the perturbation \(\hat{H}' = e \mathcal{E} z\) preserves cylindrical symmetry around \(z\), \([\hat{H}', \hat{L}_z] = 0\). Matrix elements between states with different \(m\) vanish:

• \(2p_{+1}\) (\(m = +1\)) and \(2p_{-1}\) (\(m = -1\)) have no partners with the same \(m\) and different parity, so their first-order shifts are zero.
• The \(2s\) (\(m = 0\)) and \(2p_z\) (\(m = 0\)) states possess different parities and couple strongly:

The secular determinant within the \(\{2s, 2p_z\}\) subspace is:

The degenerate \(n = 2\) level splits linearly with field strength (Linear Stark Effect):

• \(E_+^{(1)} = +3 e a_0 \mathcal{E}\) (eigenstate \(\frac{1}{\sqrt{2}}(2s - 2p_z)\), dipole pointing against field)
• \(E_0^{(1)} = 0\) (two-fold degenerate: \(2p_x, 2p_y\))
• \(E_-^{(1)} = -3 e a_0 \mathcal{E}\) (eigenstate \(\frac{1}{\sqrt{2}}(2s + 2p_z)\), dipole pointing along field)

The linear Stark effect is unique to hydrogenic systems due to their accidental \(l\)-degeneracy.

## Advanced Mathematical Supplement: SO(4) Symmetry & Runge-Lenz Invariant

Dynamical SO(4) Symmetry of the Hydrogen Atom & The Laplace-Runge-Lenz Vector

The accidental degeneracy of hydrogenic energy levels with respect to orbital angular momentum \(l\) (where \(E_n\) depends strictly on \(n\)) is not an accident. It is the direct consequence of an underlying four-dimensional rotational symmetry \(\mathbf{SO}(4)\).

The Quantum Laplace-Runge-Lenz Operator

In classical Keplerian planetary orbits, the Laplace-Runge-Lenz vector points from the focus along the major axis toward perihelion, remaining strictly constant in time. In quantum mechanics, the Hermitian Laplace-Runge-Lenz operator is:

Evaluating the commutator with the Coulomb Hamiltonian \(\hat{H} = \frac{\hat{\mathbf{p}}^2}{2\mu} - \frac{Z e^2}{4\pi\varepsilon_0 r}\):

Because both \(\hat{\mathbf{L}}\) and \(\hat{\mathbf{M}}\) commute with \(\hat{H}\), the hydrogen atom possesses six independent conserved continuous generators.

The SO(4) Lie Algebra

For bound states (\(E < 0\)), define the scaled vector:

The commutation relations between \(\hat{\mathbf{L}}\) and \(\hat{\mathbf{K}}\) are:

This is precisely the Lie algebra of the four-dimensional orthogonal rotation group \(\mathbf{SO}(4)\). By defining two decoupled commuting angular momentum vectors:

they satisfy two independent \(\mathfrak{su}(2)\) algebras:

Because \(\hat{\mathbf{L}} \cdot \hat{\mathbf{K}} = 0\), the Casimirs are equal: \(\hat{\mathbf{I}}_1^2 = \hat{\mathbf{I}}_2^2 = j(j + 1)\hbar^2\). The total Casimir relates directly to the Hamiltonian:

Setting \(2j + 1 = n\), the bound-state energy eigenvalues are obtained algebraically without ever solving differential equations:

The state degeneracy is \((2j + 1)^2 = n^2\), matching the observed spatial orbital degeneracy.

## Research Monograph: Relativistic Dirac Equation & Quantum Electrodynamic Corrections

Dirac's Relativistic Wave Mechanics

Paul Dirac sought a first-order wave equation that is linear in both time and space derivatives to ensure positive-definite probability densities while satisfying relativistic invariance \(E^2 = p^2 c^2 + m_0^2 c^4\):

where \(\Psi\) is a 4-component spinor (bispinor), and \(\boldsymbol{\alpha} = (\alpha_x, \alpha_y, \alpha_z)\) and \(\beta\) are \(4 \times 4\) Dirac matrices:

where \(\sigma_i\) are the \(2 \times 2\) Pauli spin matrices.

Consequences of the Dirac Equation

1. Intrinsic Spin \(s = 1/2\) and \(g = 2\): Electron spin and the gyromagnetic ratio \(g_s = 2\) emerge automatically from relativistic invariance without empirical introduction!

2. Antimatter (The Positron): Negative energy solutions with \(E < -m_e c^2\) led to Dirac's hole theory, predicting the existence of the positron, experimentally discovered by Carl Anderson in 1932.

3. Zitterbewegung: Rapid trembling motion of the electron position operator at frequency \(\omega \sim 2 m_e c^2 / \hbar \approx 10^{21}\text{ Hz}\) over a Compton wavelength \(\hbar / m_e c \approx 3.86 \times 10^{-13}\text{ m}\).

Quantum Electrodynamic (QED) Radiative Corrections

While Dirac's equation predicts exact degeneracy for \(2s_{1/2}\) and \(2p_{1/2}\), Quantum Electrodynamics (QED) accounts for interaction with the quantized electromagnetic vacuum:

1. Electron Self-Energy & Vacuum Fluctuations: Zero-point fluctuations of the vacuum electric field jiggle the bound electron, smearing its interaction with the nuclear Coulomb potential. Because the \(s\)-electron spends time at the nucleus while \(p\)-electrons do not, this shifts the \(2s_{1/2}\) level upward by \(+1010\text{ MHz}\).

2. Vacuum Polarization: Virtual electron-positron pairs in the vacuum shield the bare nuclear charge at distances \(r \gg \hbar / m_e c\). At very short distances, the electron penetrates the screening cloud, feeling a slightly stronger effective attraction, which lowers the \(2s_{1/2}\) level by \(-27\text{ MHz}\).

The net QED shift (\(+1057.8\text{ MHz}\)) matches the experimental Lamb shift with ten-figure accuracy, representing the pinnacle of modern theoretical physics.