Accidental Degeneracy and the SO(4) Symmetry

Line up the energy levels of any atom except hydrogen and you see a staircase: in each shell the ss level sits lowest, then pp, then dd — sodium's 3ss–3pp split is a very visible 2.1 eV, the yellow of a streetlamp. Now look at hydrogen: 2ss and 2pp, dead even; 3ss, 3pp, 3dd, dead even. Every level of the same nn coincides exactly, packing n2n^2 states into each rung. Rotational symmetry explains why the 2ℓ+12\ell+1 states of one ℓ\ell agree (Degeneracy and the Role of the Magnetic Quantum Number proves it) but is silent on why different ℓ\ell's should. Textbooks called this "accidental degeneracy." The name is a confession, not an explanation. There are no accidents of this size: a spectrum more degenerate than the visible symmetry demands means something else is conserved. Here it is the Laplace–Runge–Lenz vector — classically the arrow that pins a Kepler ellipse's orientation and forbids precession, quantum mechanically a conserved operator that shuffles ℓ\ell without touching the energy. Add its three components to L\mathbf{L}'s three and the algebra that closes is not SO(3)SO(3) but SO(4)SO(4): hydrogen secretly carries the symmetry of rotations in four dimensions. This article cashes the observation all the way: the enlarged algebra, its split into two su(2)s, and Pauli's 1926 tour de force — the full Bohr spectrum, degeneracy n2n^2 included, extracted from commutators alone, before wave mechanics touched the atom.

Key takeaways
  • Hydrogen's energy depends only on nn: each level holds ∑ℓ=0n−1(2ℓ+1)=n2\sum_{\ell=0}^{n-1}(2\ell+1) = n^2 states — more degeneracy than rotational symmetry can enforce.
  • The excess is policed by a hidden conserved operator, the Laplace–Runge–Lenz vector A^\hat{\mathbf{A}} — credentials (A^†=A^\hat{\mathbf{A}}^\dagger = \hat{\mathbf{A}}, [H^,A^]=0[\hat H, \hat{\mathbf{A}}] = 0, A^⋅L^=0\hat{\mathbf{A}}\cdot\hat{\mathbf{L}} = 0) proven in The Laplace–Runge–Lenz Vector in Quantum Mechanics.
  • On the bound spectrum, L^\hat{\mathbf{L}} and the rescaled A^′=−m/2H A^\hat{\mathbf{A}}' = \sqrt{-m/2H}\,\hat{\mathbf{A}} close the six-generator algebra so(4)\mathfrak{so}(4).
  • so(4)\mathfrak{so}(4) splits into two commuting su(2)s, J^(1,2)=12(L^±A^′)\hat{\mathbf{J}}^{(1,2)} = \tfrac12(\hat{\mathbf{L}} \pm \hat{\mathbf{A}}'); the constraint A^⋅L^=0\hat{\mathbf{A}}\cdot\hat{\mathbf{L}} = 0 forces j1=j2=jj_1 = j_2 = j.
  • Casimir bookkeeping then derives En=−mk2/2ℏ2n2E_n = -mk^2/2\hbar^2 n^2 with n=2j+1n = 2j+1 and degeneracy (2j+1)2=n2(2j+1)^2 = n^2 — Pauli's 1926 algebraic solution.
  • Key idea: "accidental" degeneracy is hidden symmetry — hydrogen's spectrum is an SO(4)SO(4) representation theory problem, and the 1/r1/r potential is exactly the potential for which that works.
The core insight

Energy depending only on nn (not ℓ\ell) signals a conserved quantity beyond angular momentum. That quantity is the Runge–Lenz vector A\mathbf{A}; with L\mathbf{L} it generates SO(4)SO(4), whose irreducible representations have dimension n2n^2 — exactly the hydrogenic degeneracy.

Where this lives today

Rydberg-atom laboratories work inside a single nn-shell the way this article does. At n∼50n \sim 50 the 2500 nearly degenerate states form a playground steered with small fields, in exactly this language: circular-state preparation (the ℓ=n−1\ell = n-1 rungs of cavity-QED) is adiabatic steering through the (j,j)(j, j) multiplet, the Runge–Lenz operator generating the ℓ\ell-transfers, and the linear Stark maps used to field-ionize Rydberg atoms read off the SO(4)SO(4) quantum numbers. The idea of a spectrum-generating symmetry — solve a system by closing the algebra of its conserved and ladder operators — became a working method across physics: the SO(4,2)SO(4,2) treatment of all hydrogen states at once, the interacting boson model's dynamical symmetries, supersymmetric factorization chains, the integrable-model toolbox — all descend from Pauli's 1926 template. And "accidental degeneracy" survives as the standard alarm phrase: when a numerical spectrum shows more coincidence than the visible symmetry explains, the modern reflex is to hunt the hidden conserved operator — in cold-atom traps, photonic lattices, and band structures alike.

The physical picture

Figure
Closed orbit, fixed arrow — the classical face of the hidden symmetry
A electron V = −k/r exactly ellipse closes; A frozen → extra conserved vector any other V(r) rosette precesses; A rotates → conservation lost the Laplace–Runge–Lenz vector points along the major axis — conserved only for 1/r

The classical tell has been in the sky all along. A planet in a pure 1/r1/r potential retraces the same ellipse forever, its long axis pointing at a fixed star. That frozen direction is a conserved vector — the Laplace–Runge–Lenz vector A\mathbf{A}, focus to perihelion, magnitude proportional to eccentricity. Perturb the potential even slightly and the ellipse becomes a turning rosette: A\mathbf{A} rotates, conservation gone (Mercury's perihelion precession is this diagnostic in reverse). Conservation of A\mathbf{A} fingerprints the exact Coulomb form. Quantum mechanics inherits the fingerprint with upgraded consequences. Rotations cannot change ℓ\ell — the L^\hat L's only reshuffle a multiplet internally. The components of A^\hat{\mathbf{A}} are different animals: they commute with H^\hat H yet carry vector (k=1k=1) character, so they connect ℓ\ell to ℓ±1\ell \pm 1. A conserved operator that changes ℓ\ell without changing energy forces degeneracy: apply A^z\hat A_z to 2ss and out comes (a multiple of) 2p0p_0 at the same energy, by conservation. The "accident" is this machine at work.

Figure
Hydrogen against a generic atom — the degeneracy the rotations cannot explain
generic atom (e.g. Na valence) 3s 3p 3d 2.1 eV ℓ levels split: core electrons spoil the pure 1/r potential hydrogen 2s 2p n=2: 4 states, one E 3s 3p 3d n=3: 9 states, one E 1s n=1: 1 state g(n) = n² — more than SO(3) can enforce: the extra generators are A's

Six conserved generators — three L^\hat L's, three A^\hat A's — are too many for the rotation group: SO(3)SO(3) has three; six is the count for SO(4)SO(4), rotations of four-dimensional Euclidean space. In 1935 Fock showed the four-dimensionality is literal: in momentum space, bound hydrogen is free motion on a 4-sphere (the portrait tells the story). Below is the operator version — close the algebra, split it into two su(2)s, let the Casimir arithmetic hand over the spectrum.

The cast: credentials quoted, then the algebra closed

Everything rests on the quantum Laplace–Runge–Lenz operator. Its construction is delicate — the naive p^×L^\hat{\mathbf p}\times\hat{\mathbf L} is not even hermitian — and the symmetrized operator, with proofs of all credentials, is the business of The Laplace–Runge–Lenz Vector in Quantum Mechanics. We quote what we need:

A^=12m(p^×L^−L^×p^)−k r^r,k=e24πε0(1)\hat{\mathbf{A}} = \frac{1}{2m}\left(\hat{\mathbf{p}}\times\hat{\mathbf{L}} - \hat{\mathbf{L}}\times\hat{\mathbf{p}}\right) - k\,\frac{\hat{\mathbf{r}}}{r}, \qquad k = \frac{e^2}{4\pi\varepsilon_0} \qquad (1)

where: mm — the (reduced) electron mass; the symmetrization makes A^\hat{\mathbf{A}} hermitian. Its credentials, all cited from that article: [H^,A^i]=0[\hat H, \hat A_i] = 0 (conserved for the Coulomb H^=p^2/2m−k/r\hat H = \hat p^2/2m - k/r, and only for it), A^⋅L^=L^⋅A^=0\hat{\mathbf{A}}\cdot\hat{\mathbf{L}} = \hat{\mathbf{L}}\cdot\hat{\mathbf{A}} = 0 (the arrow lies in the orbital plane), and the squared-length identity

A^2=2H^m(L^2+ℏ2)+k2(2)\hat{\mathbf{A}}^2 = \frac{2\hat H}{m}\left(\hat{\mathbf{L}}^2 + \hbar^2\right) + k^2 \qquad (2)

where: the lone +ℏ2+\hbar^2, absent from the classical A2=2EL2/m+k2A^2 = 2EL^2/m + k^2, is the residue of operator ordering — and precisely the term that will make nn an integer. A classical physicist would delete it as clutter; deleting it misquantizes hydrogen. The commutators among the six generators are the whole story, so they get a derivation card.

Derivation AClosing the algebra — six generators, one catch, one fix

The goal is the commutator table of {L^i,A^i}\{\hat L_i, \hat A_i\} and the discovery that it almost — then, after one rescaling, exactly — forms so(4)\mathfrak{so}(4). The strategy: two of the three bracket families cost nothing (pure transformation-character statements); the third, [A^i,A^j][\hat A_i, \hat A_j], is a genuine computation whose result we state, machine-verify, and domesticate by rescaling. Watch the catch: the [A^,A^][\hat A, \hat A] bracket closes onto L^\hat{\mathbf{L}} only through the Hamiltonian — which is what confines the clean group structure to one energy shell at a time.

Step 1. The rotation family: [L^i,L^j]=iℏ ϵijkL^k[\hat L_i, \hat L_j] = i\hbar\,\epsilon_{ijk}\hat L_k — the angular-momentum algebra itself, owned by Angular Momentum Commutation Relations, quoted verbatim.

Step 2. The vector-character family. A^\hat{\mathbf{A}} is built from vectors (p^\hat{\mathbf{p}}, L^\hat{\mathbf{L}}, r^/r\hat{\mathbf{r}}/r), so it transforms as a vector operator:

[L^i,A^j]=iℏ ϵijkA^k.[\hat L_i, \hat A_j] = i\hbar\,\epsilon_{ijk}\hat A_k.

Why free: the k=1k{=}1 grammar of "Angular Momentum Commutation Relations" and "Spherical Tensor Operators" — any vector built from vectors obeys it; nothing Coulomb-specific yet.

Step 3. The new family — the honest computation. Expanding [A^i,A^j][\hat A_i, \hat A_j] term by term (the p^×L^\hat p\times\hat L pieces against each other, the cross terms against r^/r\hat{\mathbf r}/r) yields, after substantial cancellation:

[A^i,A^j]=iℏ ϵijk(−2H^m)L^k.[\hat A_i, \hat A_j] = i\hbar\,\epsilon_{ijk}\left(-\frac{2\hat H}{m}\right)\hat L_k.

Why the Hamiltonian appears: dimensional bookkeeping demands it — the bracket has dimensions of A^2∼\hat A^2 \sim (velocity)², and energy/mass is the only conserved quantity that supplies them. The expansion is long but mechanical; we verify the identity symbolically as differential operators on test wavefunctions (footnote) rather than reproduce three pages of algebra.

Step 4. The catch and the fix. With H^\hat H inside the bracket, the six operators do not form a Lie algebra over the numbers — the structure "constants" are operators. But H^\hat H commutes with everything in sight, so on a fixed bound eigenspace (E<0E < 0) it becomes its eigenvalue. Define there

A^′=m−2E  A^⟹[A^i′,A^j′]=iℏ ϵijkL^k,\hat{\mathbf{A}}' = \sqrt{\frac{m}{-2E}}\;\hat{\mathbf{A}} \qquad\Longrightarrow\qquad [\hat A'_i, \hat A'_j] = i\hbar\,\epsilon_{ijk}\hat L_k,

Why legal: −2E/m>0-2E/m > 0 on bound states, so the root is real; on the E>0E > 0 continuum the same move needs m/2E\sqrt{m/2E} and flips the bracket's sign — the algebra becomes non-compact so(3,1)\mathfrak{so}(3,1), matching the infinite-dimensional scattering spectrum. The sign of the energy chooses the geometry.

Step 5. Read off the structure. The closed table — [L^,L^]∼L^[\hat L,\hat L] \sim \hat L, [L^,A^′]∼A^′[\hat L, \hat A'] \sim \hat A', [A^′,A^′]∼L^[\hat A', \hat A'] \sim \hat L — is the commutation pattern of the six rotation generators of 4D Euclidean space: L^\hat{\mathbf{L}} rotates the three familiar axes among themselves, A^′\hat{\mathbf{A}}' rotates each into the fourth. Hydrogen's bound shell carries so(4)\mathfrak{so}(4). Verified two independent ways: symbolically — [L^z,A^x]=iℏA^y[\hat L_z, \hat A_x] = i\hbar\hat A_y and [A^x,A^y]=−2iℏH^L^z/m[\hat A_x, \hat A_y] = -2i\hbar\hat H\hat L_z/m applied as differential operators to e−re^{-r} and ze−r/2z e^{-r/2} (sympy, exact zero residuals); and numerically in the n=2n{=}2 matrix representation, all nine brackets of each family to 10−1310^{-13}.

Two su(2)s and Pauli's spectrum

The algebra so(4)\mathfrak{so}(4) has a lovely secret: unlike its three-dimensional cousin it is not simple — it falls apart into two independent copies of the angular-momentum algebra. The right change of variables makes the split visible.

Go deeper · Derivation B · 1 minThe decoupling — so(4) is two angular momenta in disguise
Derivation BThe decoupling — so(4) is two angular momenta in disguise

The goal is to replace the six generators by two commuting triples, each a bona-fide angular momentum, so the entire su(2) representation theory imports wholesale. The strategy is the sum-and-difference combination — with a bonus: the constraint A^⋅L^=0\hat{\mathbf{A}}\cdot\hat{\mathbf{L}} = 0 becomes a statement about the two new spins that collapses the representation zoo to hydrogen's single ladder.

Step 1. Define the combinations

J^(1)=12(L^+A^′),J^(2)=12(L^−A^′).\hat{\mathbf{J}}^{(1)} = \tfrac12\left(\hat{\mathbf{L}} + \hat{\mathbf{A}}'\right), \qquad \hat{\mathbf{J}}^{(2)} = \tfrac12\left(\hat{\mathbf{L}} - \hat{\mathbf{A}}'\right).

Step 2. Check they commute:

[J^i(1),J^j(2)]=14([L^i,L^j]−[L^i,A^j′]+[A^i′,L^j]−[A^i′,A^j′])=14 iℏ ϵijk(L^k−A^k′+A^k′−L^k)=0.[\hat J^{(1)}_i, \hat J^{(2)}_j] = \tfrac14\left([\hat L_i,\hat L_j] - [\hat L_i,\hat A'_j] + [\hat A'_i,\hat L_j] - [\hat A'_i,\hat A'_j]\right) = \tfrac14\,i\hbar\,\epsilon_{ijk}\left(\hat L_k - \hat A'_k + \hat A'_k - \hat L_k\right) = 0.

Why exact: the [A^′,A^′]∼L^[\hat A',\hat A'] \sim \hat L closure — the previous card's hard-won Step 3 — enters with precisely the sign that kills the [L^,L^][\hat L, \hat L] term. The rescaling was chosen to make this happen.

Step 3. Each triple closes on itself the same way: [J^i(a),J^j(a)]=iℏ ϵijkJ^k(a)[\hat J^{(a)}_i, \hat J^{(a)}_j] = i\hbar\,\epsilon_{ijk}\hat J^{(a)}_k. Two independent angular momenta — so(4)≅su(2)⊕su(2)\mathfrak{so}(4) \cong \mathfrak{su}(2) \oplus \mathfrak{su}(2). Why this is a gift: the representation theory of "Ladder Operators and the Angular Momentum Spectrum" now applies twice over — eigenvalues ℏ2ja(ja+1)\hbar^2 j_a(j_a+1), ja∈{0,12,1,… }j_a \in \{0, \tfrac12, 1, \dots\}, half-integers included because nothing here is orbital.

Step 4. Import the constraint. Expanding the squares, [J^(1)]2−[J^(2)]2=L^⋅A^′[\hat{\mathbf J}^{(1)}]^2 - [\hat{\mathbf J}^{(2)}]^2 = \hat{\mathbf{L}}\cdot\hat{\mathbf{A}}' — and A^⋅L^=0\hat{\mathbf{A}}\cdot\hat{\mathbf{L}} = 0 is a proven credential of the operator:

[J^(1)]2=[J^(2)]2⟹j1=j2≡j.[\hat{\mathbf{J}}^{(1)}]^2 = [\hat{\mathbf{J}}^{(2)}]^2 \qquad\Longrightarrow\qquad j_1 = j_2 \equiv j.

Why this halves the world: a generic so(4)\mathfrak{so}(4) representation is an independent pair (j1,j2)(j_1, j_2). The orbit-geometry fact that A\mathbf{A} lies in the orbital plane strikes out everything off the diagonal: hydrogen only builds (j,j)(j,j). *Verified in the explicit (12,12)(\tfrac12,\tfrac12) matrix representation: all nine cross-commutators vanish, both triples close su(2), L^⋅A^′=0\hat{\mathbf L}\cdot\hat{\mathbf A}' = 0 identically (numpy, 10−1310^{-13}).*

With the machinery assembled, the spectrum falls out of bookkeeping. This is the calculation Wolfgang Pauli published in January 1926 — solving hydrogen by matrix mechanics while the differential-equation route did not yet exist.

Go deeper · Derivation C · 1 minPauli 1926 — the Bohr spectrum from Casimir arithmetic
Derivation CPauli 1926 — the Bohr spectrum from Casimir arithmetic

The goal is hydrogen's complete bound spectrum and degeneracies from the algebra alone: no radial equation, no special functions, no boundary conditions. The strategy: evaluate the remaining Casimir, [J^(1)]2+[J^(2)]2[\hat{\mathbf J}^{(1)}]^2 + [\hat{\mathbf J}^{(2)}]^2, two ways — by representation theory (2ℏ2j(j+1)2\hbar^2 j(j+1)) and by the operator identity (2) — and set them equal. The energy has nowhere to hide.

Step 1. Express the Casimir sum in the original generators, using L^⋅A^′+A^′⋅L^=0\hat{\mathbf L}\cdot\hat{\mathbf A}' + \hat{\mathbf A}'\cdot\hat{\mathbf L} = 0:

[J^(1)]2+[J^(2)]2=12(L^2+A^′2).[\hat{\mathbf{J}}^{(1)}]^2 + [\hat{\mathbf{J}}^{(2)}]^2 = \tfrac12\left(\hat{\mathbf{L}}^2 + \hat{\mathbf{A}}'^2\right).

Step 2. Evaluate A^′2\hat{\mathbf{A}}'^2 on the energy-EE shell via identity (2):

A^′2=m−2E A^2=m−2E[2Em(L^2+ℏ2)+k2]=−L^2−ℏ2−mk22E.\hat{\mathbf{A}}'^2 = \frac{m}{-2E}\,\hat{\mathbf{A}}^2 = \frac{m}{-2E}\left[\frac{2E}{m}(\hat{\mathbf{L}}^2 + \hbar^2) + k^2\right] = -\hat{\mathbf{L}}^2 - \hbar^2 - \frac{mk^2}{2E}.

Why the signs work: E<0E < 0 on bound states, so −mk2/2E>0-mk^2/2E > 0.

Step 3. Add L^2\hat{\mathbf{L}}^2 and watch it cancel:

12(L^2+A^′2)=12(−ℏ2−mk22E).\tfrac12\left(\hat{\mathbf{L}}^2 + \hat{\mathbf{A}}'^2\right) = \tfrac12\left(-\hbar^2 - \frac{mk^2}{2E}\right).

Why this matters: the Casimir — hence the energy — is blind to how the shell's angular momentum is distributed. The independence of EE from ℓ\ell stops being an output to explain; it is built into the bookkeeping at this line.

Step 4. Set the two evaluations equal, with j1=j2=jj_1 = j_2 = j from the decoupling card:

2ℏ2j(j+1)=12(−ℏ2−mk22E)⟹−mk22E=ℏ2[4j(j+1)+1]=ℏ2(2j+1)2.2\hbar^2 j(j+1) = \tfrac12\left(-\hbar^2 - \frac{mk^2}{2E}\right) \quad\Longrightarrow\quad -\frac{mk^2}{2E} = \hbar^2\left[4j(j+1) + 1\right] = \hbar^2(2j+1)^2.

Why the square is perfect: 4j(j+1)+1=(2j+1)24j(j+1) + 1 = (2j+1)^2 — exactly where the quantum +ℏ2+\hbar^2 of identity (2) earns its keep. Without it the bracket is 4j(j+1)4j(j+1), no perfect square, and the spectrum does not come out as 1/integer21/\text{integer}^2.

Step 5. Solve, defining n≡2j+1n \equiv 2j+1:

En=−mk22ℏ2n2=−13.6057 eVn2,n=1,2,3,…E_n = -\frac{mk^2}{2\hbar^2 n^2} = -\frac{13.6057\,\text{eV}}{n^2}, \qquad n = 1, 2, 3, \dots

Since jj runs over 0,12,1,…0, \tfrac12, 1, \dots, nn runs over all positive integers. Half-integer jj — illegal for any single orbital angular momentum — is not merely allowed but required: it builds the even-nn shells.

Step 6. Count the states. The (j,j)(j,j) representation has dimension (2j+1)2=n2(2j+1)^2 = n^2 — the degeneracy, derived rather than summed. Decomposing L^=J^(1)+J^(2)\hat{\mathbf{L}} = \hat{\mathbf{J}}^{(1)} + \hat{\mathbf{J}}^{(2)} by the Clebsch–Gordan series (Addition of Angular Momenta): j⊗j=0⊕1⊕⋯⊕2jj \otimes j = 0 \oplus 1 \oplus \cdots \oplus 2j — exactly one ℓ\ell-multiplet for each ℓ=0,…,n−1\ell = 0, \dots, n-1. The constraint ℓ≤n−1\ell \le n-1, a series termination in The Hydrogen Atom: Coulomb Radial Equation, is here a line of group theory. Verified: Casimir route at j=0,12,1j = 0, \tfrac12, 1 == the ladder −13.6057/n2-13.6057/n^2 eV exactly; perfect-square identity symbolic; CG counts ∑(2ℓ+1)=n2\sum(2\ell+1) = n^2 for n≤6n \le 6; matches Pauli (1926), Sakurai §4.1.

En=−mk22ℏ2n2,gn=n2,n=2j+1(3)E_n = -\frac{mk^2}{2\hbar^2 n^2}, \qquad g_n = n^2, \qquad n = 2j+1 \qquad (3)

where: gng_n — states per shell (spin doubles it to 2n22n^2 but adds no ℓ\ell-physics). Compare the chapter's two derivations of hydrogen: the analytic route needed a differential equation, an asymptotic peel, a series, and a termination demand; the algebraic route needed six commutators and one constraint. The price of elegance is fragility — the algebra knows only the exact 1/r1/r problem. One more reading deserves a paragraph. Fock's 1935 construction: stereographically project each bound momentum state, from the momentum sphere of radius p0=−2mEp_0 = \sqrt{-2mE}, onto the unit sphere S3S^3. There hydrogen's Schrödinger equation becomes a free particle's on S3S^3 — symmetry manifest — and the 4D spherical harmonics come in families of dimension exactly n2n^2: hydrogen's degeneracy is the statement that a hypersphere looks the same from every side. (Stated without derivation; the portrait band tells the story.)

Limiting cases

casesymmetry and spectrummeaning
exact Coulomb, E<0E<0SO(4)SO(4): En=−Ry/n2E_n = -\text{Ry}/n^2, gn=n2g_n = n^2full accidental degeneracy — this article
exact Coulomb, E>0E>0SO(3,1)SO(3,1) continuumrescaling flips sign; non-compact, infinite-dimensional representations
fine structure onEE depends on (n,j)(n, j)SO(4)SO(4) broken; 2s1/22s_{1/2}–2p1/22p_{1/2} survives until the Lamb shift — approximation chapters
screened potential (alkali)Enℓ=−Ry/(n−δℓ)2E_{n\ell} = -\text{Ry}/(n-\delta_\ell)^2core penetration spoils 1/r1/r: ℓ\ell-degeneracy lifts, mm survives
isotropic 3D oscillatorSU(3)SU(3): gN=12(N+1)(N+2)g_N = \tfrac12(N{+}1)(N{+}2)the other closed-orbit potential — Degeneracy and SU(3) Symmetry of the 3D Oscillator
uniform external fieldSO(4)→SO(4) \to lessa chosen axis; hydrogen's linear Stark effect is the degeneracy made visible

Worked examples

1 · The n = 2 shell, assembled from two spin-halves (confidence builder)

Work in the smallest nontrivial representation, n=2n = 2: j1=j2=12j_1 = j_2 = \tfrac12. Build the generators, verify the algebra, and watch the shell's orbital content — one 2s, three 2p — drop out of angular-momentum addition.

Step 1. Represent the two spins on C2⊗C2\mathbb{C}^2\otimes\mathbb{C}^2: J^(1)=ℏ2σ⊗1\hat{\mathbf{J}}^{(1)} = \tfrac{\hbar}{2}\boldsymbol\sigma\otimes\mathbb{1}, J^(2)=1⊗ℏ2σ\hat{\mathbf{J}}^{(2)} = \mathbb{1}\otimes\tfrac{\hbar}{2}\boldsymbol\sigma — 4×44\times4 matrices. (Symmetry spins, not electron spin: the four dimensions are the n=2n=2 shell's four states.)

Step 2. Machine-check the structure (this article's verify script): all nine [J^i(1),J^j(2)]=0[\hat J^{(1)}_i, \hat J^{(2)}_j] = 0; each triple closes su(2); L^=J^(1)+J^(2)\hat{\mathbf{L}} = \hat{\mathbf{J}}^{(1)} + \hat{\mathbf{J}}^{(2)} and A^′=J^(1)−J^(2)\hat{\mathbf{A}}' = \hat{\mathbf{J}}^{(1)} - \hat{\mathbf{J}}^{(2)} reproduce the full so(4)\mathfrak{so}(4) table; and L^⋅A^′=0\hat{\mathbf{L}}\cdot\hat{\mathbf{A}}' = 0 exactly — automatic in a (j,j)(j,j) representation.

Step 3. Decompose by total angular momentum: diagonalizing the 4×44\times4 L^2\hat{\mathbf{L}}^2 gives eigenvalues {0,2ℏ2,2ℏ2,2ℏ2}\{0, 2\hbar^2, 2\hbar^2, 2\hbar^2\} — one ℓ=0\ell = 0 state and three ℓ=1\ell = 1 states, the 2s and the 2p triplet: the 12⊗12=0⊕1\tfrac12\otimes\tfrac12 = 0 \oplus 1 Clebsch–Gordan series in the flesh.

Step 4. Confirm the energy label: n=2j+1=2n = 2j+1 = 2, and the Casimir arithmetic gives E2=−13.6057/4=−3.401 eVE_2 = -13.6057/4 = \mathbf{-3.401\,eV} for all four states. The 2s–2p degeneracy is not two integrals that happened to agree; it is the statement that both multiplets live in one irreducible (12,12)(\tfrac12, \tfrac12) box.

2 · Running the Casimir constraint up the ladder (exam level)

Evaluate the spectrum at j=0,12,1j = 0, \tfrac12, 1 against the radial-equation ladder; then find which states the Runge–Lenz operator annihilates.

Step 1. The constraint reads −mk2/2E=ℏ2(2j+1)2-mk^2/2E = \hbar^2(2j+1)^2, with mk2/2ℏ2=13.6057 eVmk^2/2\hbar^2 = 13.6057\,\text{eV} (the Rydberg — assembled from CODATA in The Bohr Radius and Characteristic Atomic Scales): j=0⇒−13.606 eVj = 0 \Rightarrow \mathbf{-13.606\,eV}; j=12⇒−3.401 eVj = \tfrac12 \Rightarrow \mathbf{-3.401\,eV}; j=1⇒−1.512 eVj = 1 \Rightarrow \mathbf{-1.512\,eV}. All three match the −Ry/n2-\text{Ry}/n^2 ladder of The Hydrogen Atom: Coulomb Radial Equation exactly — two derivations, zero daylight.

Step 2. Multiplet sizes: (2j+1)2=1,4,9(2j+1)^2 = 1, 4, 9 — with electron spin, 2n2=2,8,182n^2 = 2, 8, 18: the periodic table's first row lengths emerge as SO(4)SO(4) representation dimensions (the filling order beyond hydrogen involves screening, which breaks the symmetry being counted).

Step 3. Where does A^′\hat{\mathbf{A}}' run out of work? On an (n,ℓ)(n, \ell) level, Derivation C's Step 2 gives A^′2=ℏ2[n2−ℓ(ℓ+1)−1]\hat{\mathbf{A}}'^2 = \hbar^2[n^2 - \ell(\ell+1) - 1]: at the top rung ℓ=n−1\ell = n-1 this is ℏ2(n−1)\hbar^2(n-1), and it vanishes only at the ground state (1,0)(1, 0) — the one shell with no other ℓ\ell to shuffle to.

Step 4. Positivity sanity check: n2−ℓ(ℓ+1)−1≥0n^2 - \ell(\ell+1) - 1 \ge 0 demands ℓ≤n−1\ell \le n-1 — the operator identity knows the allowed quantum numbers. (Verified for n≤7n \le 7 in the script.)

Origins · 1 minVladimir Fock and how this was discovered
VF
Vladimir Fock
1898–1974 · Leningrad

Fock's 1935 paper "Zur Theorie des Wasserstoffatoms" answered the question Pauli's 1926 algebra left hanging: the commutators said hydrogen carries four-dimensional rotation symmetry — but what, physically, rotates? Fock's move was to leave position space entirely. Write the bound-state problem in momentum space at energy E=−p02/2mE = -p_0^2/2m, and stereographically project the momentum sphere of radius p0p_0 onto the unit sphere S3S^3. Under that map the Schrödinger integral equation becomes that of a free particle on the hypersphere: the potential dissolves into geometry, and the "hidden" SO(4)SO(4) is the statement that a sphere looks the same from every direction. The four-dimensional harmonics fall into families of dimension 1,4,9,16,…1, 4, 9, 16, \dots — the degeneracies, now as transparent as the 2ℓ+12\ell+1 of ordinary rotations. A year later Bargmann closed the loop, proving Fock's geometric group and Pauli's algebra identical. Working in Leningrad through the purge years (arrested twice, briefly; defender of relativity and quantum theory against ideological attack), Fock left his name on much of later many-body theory — other articles' stories. Here he is the man who found the fourth dimension hiding in the hydrogen atom.

History
1873
— Bertrand proves only 1/r and r² close every bounded orbit — the classical shadow of hidden symmetry.
1926
— Pauli derives the full hydrogen spectrum from the Runge–Lenz algebra, weeks before Schrödinger's wave-equation solution appears.
1935
— Fock maps bound hydrogen onto a free particle on the 4-sphere: the symmetry becomes geometry.
1936
— Bargmann proves Fock's geometric SO(4) and Pauli's operator algebra coincide.
1947
— Lamb and Retherford measure the 2s₁/₂–2p₁/₂ splitting: the last "accidental" degeneracy falls, to QED.

Going further

Texts. Sakurai & Napolitano, 3rd ed., §4.1 + Runge–Lenz/SO(4) supplement · Shankar, 2nd ed., Ch. 13 · Cohen-Tannoudji Vol. I, Runge–Lenz complement · Griffiths & Schroeter §4.2 · Pauli, Z. Phys. 36, 336 (1926) · Fock, Z. Phys. 98, 145 (1935) · Bargmann, Z. Phys. 99, 576 (1936).

Provenance. Spectrum and degeneracy: En=−13.6057 eV/n2E_n=-13.6057\,\text{eV}/n^2, gn=n2g_n=n^2 — verify against Bohr/Schrödinger solution (Griffiths §4.2, Sakurai §4). Runge–Lenz operator, its Hermitian symmetrization, conservation [H,A]=0[H,\mathbf{A}]=0, orthogonality A⋅L=0\mathbf{A}\cdot\mathbf{L}=0, and the identity A2=2Hμ(L2+ℏ2)+k2\mathbf{A}^2=\tfrac{2H}{\mu}(\mathbf{L}^2+\hbar^2)+k^2 — recompute by direct commutator algebra (Shankar Ch. 13; Cohen-Tannoudji complement). SO(4)≅SU(2)×SU(2)SO(4)\cong SU(2)\times SU(2) decomposition with j1=j2=(n−1)/2j_1=j_2=(n-1)/2 and the Casimir derivation of EnE_n — reproduce Pauli (1926). Clebsch–Gordan series j⊗j=⨁ℓ=02jℓj\otimes j=\bigoplus_{\ell=0}^{2j}\ell gives ℓ=0…n−1\ell=0\dots n-1. Bertrand's theorem (closed orbits only for 1/r1/r and r2r^2) — classical mechanics (Goldstein). No nonstandard numbers used; the only numeric value is the Rydberg energy 13.605713.6057 eV. Edition note (2026-07-08). The Runge–Lenz operator's construction, hermiticity, conservation, orthogonality, and squared-length identity are quoted with their owner named — The Laplace–Runge–Lenz Vector in Quantum Mechanics proves them; this article owns the closure, the su(2)⊕su(2) decoupling, and the Casimir spectrum. The base's second derivation (a matching-condition sketch that only V∝1/rV \propto 1/r conserves A\mathbf{A}) is not reproduced: its conclusion is stated via Bertrand's theorem, and operator-level uniqueness stays with the operator's article. All algebra machine-verified (symbolic differential operators + the explicit n=2n=2 representation).

Prerequisites. The hydrogen atom: Coulomb radial equation · the Bohr energy spectrum and Rydberg constant · degeneracy and the role of the magnetic quantum number · the Laplace–Runge–Lenz vector in quantum mechanics. Leads to. Hydrogenic wavefunctions and quantum numbers · hydrogenic ions and scaling with nuclear charge.

Is anything actually rotating in four dimensions, or is SO(4) just bookkeeping?

Something concrete rotates — in momentum space. Fock's map sends each bound state onto the sphere S3S^3, where hydrogen's dynamics is free motion on the hypersphere: three generators are ordinary rotations of p\mathbf{p}-space (L^\hat{\mathbf{L}}), the other three tilt the projection (A^′\hat{\mathbf{A}}'). No fourth spatial dimension exists in the lab — the geometry lives in the state space.

Why does the same trick not work for other atoms — helium, say?

Because the other electrons wreck the 1/r1/r form. Helium's electron feels the nucleus plus the other electron's cloud: central, but not Coulombic — so [H^,A^]≠0[\hat H, \hat{\mathbf{A}}] \ne 0 and the ℓ\ell-ladder collapses. That is why alkali spectra show the staircase hydrogen lacks, quantified by the quantum defect in Enℓ=−Ry/(n−δℓ)2E_{n\ell} = -\text{Ry}/(n-\delta_\ell)^2: δℓ\delta_\ell measures how deeply orbit ℓ\ell penetrates the core — how badly the hidden symmetry is broken for that ℓ\ell.

The two J spins can be half-integer — but orbital angular momentum forbids half-integers. Is there a contradiction?

No — single-valuedness constrains operators that differentiate real angles, which is what L^\hat{\mathbf{L}} is. The combinations J^(1,2)=12(L^±A^′)\hat{\mathbf{J}}^{(1,2)} = \tfrac12(\hat{\mathbf{L}} \pm \hat{\mathbf{A}}') are not orbital angular momenta of anything — no angle generates them — so the full algebraic spectrum applies, half-integers included. The observable L^=J^(1)+J^(2)\hat{\mathbf{L}} = \hat{\mathbf{J}}^{(1)} + \hat{\mathbf{J}}^{(2)} still comes out integer: two equal spins add to integer totals. The half-integers live in the scaffolding, never in the observable.

Does electron spin change the SO(4) story?

Only by a factor of two in the counting: spin is absent from the non-relativistic Coulomb Hamiltonian, so each of the n2n^2 orbital states carries both projections, gn=2n2g_n = 2n^2. Spin's real entrance is relativistic — spin–orbit coupling breaks SO(4)SO(4) and produces the fine-structure jj-dependence of the table above. The good bookkeeping becomes J^=L^+S^\hat{\mathbf{J}} = \hat{\mathbf{L}} + \hat{\mathbf{S}}, and the accidental degeneracy is over.

What exactly breaks first: the algebra, the conservation, or the degeneracy?

The conservation. Any δV\delta V not ∝1/r\propto 1/r gives [H^,A^]≠0[\hat H, \hat{\mathbf{A}}] \ne 0 immediately — the operator still exists, it just stops commuting. The shell-by-shell closure then fails, and the degeneracy lifts at whatever order δV\delta V connects the ℓ\ell-multiplets. The breaking is graded: fine structure preserves rotations and parity, so it splits only by jj; a uniform field preserves one rotation axis, so it splits in ∣m∣|m| patterns. Reading splittings against surviving subgroups is the spectroscopist's version of this article.

Is the n² pattern why the periodic table has rows of length 2, 8, 18, 32?

Half of why. The magic numbers 2n22n^2 are the SO(4)SO(4) dimensions doubled by spin — the shell capacities. The table's actual row lengths (2, 8, 8, 18, 18, 32…) reflect filling order, where screening breaks the ℓ\ell-degeneracy and interleaves shells (4s before 3d). SO(4)SO(4) sizes the boxes; symmetry breaking orders their use.

Before this

Where this leads