Madelung Rules: Orbital Filling Order

• Madelung Rules are a set of principles that define the electron orbital filling order based on the sum (n+ℓ) and the principal quantum number n.
• Analytic derivations reveal their foundation in central potentials with hidden O(4) symmetry, which explains observed energy degeneracies and periodic trends.
• Recent studies incorporate topological perspectives and account for electron correlations and relativistic effects, clarifying exceptions in d- and f-block elements.

The Madelung rules, also known as the (n + ℓ, n) principle or Aufbau principle, formalize the empirically observed ordering in which electronic subshells (atomic orbitals) are filled in atoms. This filling sequence underpins the structure of the periodic table, the energetic ordering of orbitals, and a broad range of atomic properties. The Madelung rule asserts: orbitals are filled in order of increasing sum n + ℓ, where n is the principal quantum number and ℓ is the azimuthal quantum number; when two orbitals have the same value of n + ℓ, the orbital with the lower n is filled first. Recent research has provided analytic proofs of the rule as an emergent property of certain mean-field atomic Hamiltonians and explicated the group-theoretic and topological structures underpinning its validity and exceptions.

1. Formal Statement and Empirical Manifestation

The Madelung rule specifies that atomic orbitals are occupied in the order determined by the two integers (n + ℓ, n). Specifically,

• Orbitals are filled with increasing values of N = n + ℓ.
• For fixed N, lower values of n take precedence. This ordering directly yields the observed periodicities in the lengths of periods of the periodic table (2, 8, 8, 18, 18, 32, ...), and the familiar sequence of subshell filling (1s, 2s, 2p, 3s, 3p, 4s, 3d, ...).

Mathematically, for two orbitals (n₁, ℓ₁) and (n₂, ℓ₂):

• If (n₁ + ℓ₁) < (n₂ + ℓ₂), fill (n₁, ℓ₁) before (n₂, ℓ₂).
• If (n₁ + ℓ₁) = (n₂ + ℓ₂), fill the one with smaller n (Loza et al., 2012).

This sequence corresponds to the following sequence of shells:

$1s < 2s < 2p < 3s < 3p < 4s < 3d < 4p \ldots$

2. Analytic Derivation and the Role of Central Potentials

The Madelung rule can be derived as a consequence of the energy spectrum of electron states in certain effective central potentials that admit a hidden dynamical O(4) symmetry (Belokolos, 2017, Kholodenko et al., 2019, Kholodenko, 2020). Specifically, the mean-field atomic Hamiltonian,

$$H = \sum_{k=1}^Z\left(-\frac12 \Delta_k + V(r_k)\right)$$

(with central potential V) admits exact O(4) symmetry only if V takes a specific form. This form, identified as the Tietz potential,

$V(r) = -\frac{Z}{r (1 + r/R)^2}$

where $R \sim Z^{-1/3}$, yields a spectrum for single-particle energies depending solely on the combination $n + \ell$ (Belokolos, 2017).

Semiclassical analysis (using the Bohr–Sommerfeld quantization near $E = 0$) and quantum algebraic methods (exploiting the commutation relations of angular momentum $\mathbf L$ and a Runge–Lenz-type vector $\mathbf A$), reveal that the spectrum is governed by the Casimir invariant

$$C_1 = L^2 + M^2 \;\longrightarrow\; p\,n_r + q\,\ell$$

with $q/p = 2$, thus giving $$H = \sum_{k=1}^Z\left(-\frac12 \Delta_k + V(r_k)\right)$$0, exactly reproducing the Madelung rule and period lengths in the Mendeleev table (Belokolos, 2017).

A full quantum mechanical treatment shows that for the Tietz-type or "deformed fish-eye" potentials—derivable as slight modifications to the Coulomb potential and justifiable by geometric (Bertrand spacetime) and conformal invariance techniques—the eigenvalues take the explicit form (Kholodenko et al., 2019):

$$H = \sum_{k=1}^Z\left(-\frac12 \Delta_k + V(r_k)\right)$$1

confirming that all states with the same $$H = \sum_{k=1}^Z\left(-\frac12 \Delta_k + V(r_k)\right)$$2 are (accidentally) degenerate.

3. Beyond Spherical Symmetry: Topology, Conformal Invariance, and Hopf Mapping

Advances in mathematical physics have clarified the geometric and topological underpinnings of the Madelung rule. The effective many-electron potential, when expressed as a "type I" Bertrand potential, can be interpreted as encoding the metric structure of a curved, spherically symmetric, conformally invariant 3-manifold (Kholodenko et al., 2019). Lifting the problem to the three-sphere $$H = \sum_{k=1}^Z\left(-\frac12 \Delta_k + V(r_k)\right)$$3 via stereographic projection and implementing the Hopf fibration $$H = \sum_{k=1}^Z\left(-\frac12 \Delta_k + V(r_k)\right)$$4 allows for a treatment of the Schrödinger equation that recasts the spectrum as that of a Laplacian on $$H = \sum_{k=1}^Z\left(-\frac12 \Delta_k + V(r_k)\right)$$5.

Crucially, this geometric framework:

• Enforces conformal invariance via the Yamabe Laplacian.
• Explains the $$H = \sum_{k=1}^Z\left(-\frac12 \Delta_k + V(r_k)\right)$$6 degeneracy rigorously via the unique properties of the three-sphere's symmetry generators.
• Clarifies the correspondence between the spectral ordering and the Madelung rule.

4. Exceptions: Correlation, Relativistic Effects, and Statistical Models

Empirically, the Madelung rule exhibits systematic violations, predominantly among d- and f-block elements. Of the ~104 elements with established ground states, 21 (~20%) violate the rule (Loza et al., 2012). These violations cluster among late transition metals, lanthanides, and actinides.

Several analytic and computational approaches delineate the boundaries of Madelung validity:

• Real atomic Hamiltonians deviate from idealized spherically symmetric one-body scenarios due to nonsphericity, spin-orbit coupling, and many-body electron correlation effects. Such deviations break the exact O(4) symmetry, lifting degeneracies and inducing orbital re-orderings (Kholodenko et al., 2019, Kholodenko, 2020).
• Relativistic corrections, most importantly those incorporated via the Dirac–Coulomb equation, generate fine-structure splittings and, for sufficiently high $$H = \sum_{k=1}^Z\left(-\frac12 \Delta_k + V(r_k)\right)$$7, quantitatively reverse the energy ordering of shells, producing Madelung-exceptional atoms (Kholodenko, 2020).
• Empirical violation statistics reveal that s- and p-block elements invariably obey the rule (SP-rule), as do closed shells (C-rule) and half-filled subshells (M-rule); any subshell with eight electrons always violates (8-rule) (Loza et al., 2012).
• The frequency and nature of violations can be modeled statistically via two independent factors: electron count $$H = \sum_{k=1}^Z\left(-\frac12 \Delta_k + V(r_k)\right)$$8 in the open subshell and hole count $$H = \sum_{k=1}^Z\left(-\frac12 \Delta_k + V(r_k)\right)$$9. The probability that an element with $V(r) = -\frac{Z}{r (1 + r/R)^2}$0 will obey Madelung is

$V(r) = -\frac{Z}{r (1 + r/R)^2}$1

where $V(r) = -\frac{Z}{r (1 + r/R)^2}$2 are empirical counts for obeying/violating elements, providing high correlation ($V(r) = -\frac{Z}{r (1 + r/R)^2}$3) with observed configurations (Loza et al., 2012).

Subshell	Simple Violation	Double Violation	Complex Violation	Always Obey
s, p blocks	0	0	0	Yes
d block	38.5%	3.8%	3.8%
f block	33.3%	4.1%	0
Closed shell	0	0	0	Yes

5. Extensions: Regularized Spectral Ordering and Madelung’s Mathematical Origin

Recent research employing canonical regularization in the de Broglie–Bohm representation has shown that even for single-electron (hydrogenic) atoms, Sturm-Liouville boundary conditions and Langer-type corrections lift Coulomb degeneracy and analytically reproduce a spectral sequence conforming to Madelung’s principle (Kumar, 15 Jun 2026). The analytic solutions involve generalized Laguerre, Legendre, and Bessel functions with non-integer parameters. The resulting orbital-dependent shifts yield the $V(r) = -\frac{Z}{r (1 + r/R)^2}$4 ordering observed in empirical filling sequences, with the canonical regularization acting as a minimal mathematical mechanism underlying the Aufbau principle, independent of conventional many-body screening arguments.

6. Physical and Chemical Implications

The Madelung rule encodes not only the periodic repetition of elemental properties, but also the observed energetic regularities governing chemical reactivity, ionization series, and bulk band-structure trends. The departures (Madelung-exceptional atoms) are especially salient in d- and f-block series where exchange, relativistic, and crystal field effects alter ground-state configurations. Notably, it has been shown that Madelung-exceptional atoms—those for which single-atom Dirac-Coulomb spectra cross at physical Z—play a disproportionately large role in the design of high-T_c superconductors, as the topological and group-theoretic nature of their ground states maps onto favourable conditions for BCS-like condensation when assembled into hydrides or complex lattices (Kholodenko, 2020).

The analytic frameworks recently developed also enable statistically robust predictions for superheavy elements (Z > 100), guiding both atomic spectroscopy and quantum-chemical modeling (Loza et al., 2012).

7. Mathematical and Topological Perspectives

The most recent theoretical advances interpret the appearance of Madelung regularity as intimately connected to the topology of the underlying atomic manifold. The transition from Madelung-regular to Madelung-exceptional atoms is cast as a topological phase change, detectable via the properties of the associated Seiberg–Witten monopole equations and spin^c structures (Kholodenko, 2020). The presence of a nontrivial moduli space of solutions for specific Z maps onto the existence of atomic superconducting order parameters at the single-atom level, with implications extending to solid-state high-T_c phenomena.

In summary, the Madelung rules arise from a confluence of dynamical symmetry (O(4)), analytic properties of central potentials (Tietz or deformed fish-eye forms), conformal geometry, and topological field theory. Their predictive power remains central to both atomic structure theory and the rational design of advanced materials. Exceptions and anomalies are analytically understood as consequences of symmetry breaking, relativistic reordering, and many-body effects, setting the stage for precise modeling of both known and as-yet-unknown elements (Belokolos, 2017, Kholodenko et al., 2019, Kumar, 15 Jun 2026, Loza et al., 2012, Kholodenko, 2020).