PHYSICAL CHEMISTRY 2 • ATOMIC AND MOLECULAR STRUCTURE

Atomic Orbitals & Electron Configurations — Atomic orbitals and electron configurations (QM framing)

How quantum mechanics replaces classical orbits with probability distributions to predict atomic structure and chemical behavior.

Historical Context & Motivation

The quest to understand how electrons inhabit atoms spans more than a century. Classical physics predicted that accelerating charges should radiate energy continuously, implying that atoms ought to collapse—a catastrophe that demanded a radically new theoretical framework. The Bohr model of 1913 introduced the idea of quantized energy levels and successfully reproduced the hydrogen emission spectrum, but it could not explain multi-electron atoms, fine spectral structure, or the chemical bond. These shortcomings compelled physicists to develop wave mechanics, a framework in which the electron is described not by a definite trajectory but by a wavefunction whose square modulus yields a probability density. This paradigm shift gave rise to the concept of atomic orbitals and the systematic rules for electron configurations that underpin all of modern chemistry.

1913
Bohr's Quantized Orbits
Niels Bohr postulated stationary orbits for hydrogen, deriving the Rydberg formula from angular-momentum quantization. The model succeeded for one-electron systems but failed for helium and beyond.
1925
Heisenberg & Pauli
Werner Heisenberg published matrix mechanics, and Wolfgang Pauli introduced the exclusion principle, stipulating that no two electrons in an atom can share the same set of four quantum numbers.
1926
Schrödinger's Wave Equation
Erwin Schrödinger formulated the time-independent wave equation for the hydrogen atom, yielding exact solutions—spherical harmonics multiplied by radial functions—that define atomic orbitals.
1928
Hartree's Self-Consistent Field
Douglas Hartree proposed an iterative method to approximate multi-electron wavefunctions by treating each electron as moving in the average potential of all others, laying the groundwork for modern computational chemistry.
1930s
Slater & the Aufbau Framework
John C. Slater introduced determinantal wavefunctions and orbital-energy ordering rules, enabling the systematic construction of ground-state electron configurations via the Aufbau principle.

The central question that this lesson addresses is: How do the quantum-mechanical solutions for hydrogen-like atoms generalize to many-electron systems, and what principles govern the filling of orbitals to determine ground-state electron configurations? Answering this question requires understanding the quantum numbers that label each orbital, the spatial form of the associated wavefunctions, and the rules—Aufbau, Pauli, and Hund—that dictate how electrons populate these orbitals.

Core Principles & Definitions

An atomic orbital is a one-electron wavefunction ψn,l,m(r, θ, φ) obtained by solving the Schrödinger equation for an electron in the Coulomb potential of a nucleus. Each orbital is labeled by three quantum numbers—n, l, and ml—and a fourth quantum number ms specifies the electron's intrinsic spin projection. The square modulus |ψ|² gives the probability density for finding the electron at a given point in space, replacing the Bohr model's deterministic trajectory with a statistical cloud. Five foundational ideas govern how these orbitals are populated in multi-electron atoms.

1

Quantum Numbers

Each orbital is uniquely specified by n (principal, 1, 2, 3 …), l (azimuthal, 0 to n−1), and ml (magnetic, −l to +l). The spin quantum number ms = ±½ completes the description.
2

Pauli Exclusion Principle

No two electrons in an atom may possess an identical set of all four quantum numbers. Consequently, each spatial orbital accommodates at most two electrons, which must have opposite spin projections.
3

Aufbau Principle

In the ground state, electrons fill orbitals in order of increasing effective energy. For multi-electron atoms, orbital energies depend on both n and l due to electron–electron repulsion and shielding effects.
4

Hund's Rule

Among degenerate orbitals (same n and l), electrons occupy them singly with parallel spins before pairing. This maximizes total spin and minimizes electron–electron repulsion.
5

Shielding & Penetration

Inner electrons screen the nuclear charge felt by outer electrons. Orbitals with lower l penetrate closer to the nucleus and experience less shielding, which explains why 4s fills before 3d in potassium.
KEY TAKEAWAY
Think of atomic orbitals as the resonance modes of a three-dimensional drum membrane clamped around a nucleus. Just as a drumhead can vibrate only in specific standing-wave patterns characterized by the number of nodal lines, an electron in an atom adopts only certain standing-wave patterns—each labeled by quantum numbers. The Aufbau principle is the rule that tells you which drum modes get 'excited' first when you add energy (electrons) to the system, while Hund's rule reflects the physical preference for distributing energy evenly across degenerate modes rather than concentrating it in one.

Visual Explanation — Orbital Shapes & Nodes

The angular part of each orbital wavefunction is a spherical harmonic Ylm(θ, φ), whose shape depends on l and ml. The following diagram presents cross-sectional sketches of the 1s, 2s, 2p, and 3d orbitals, highlighting radial and angular nodes. Recall that an orbital possesses n − l − 1 radial nodes (spherical surfaces where ψ = 0) and l angular nodes (planes or cones), yielding a total of n − 1 nodes.

Cross-sectional sketches of the 1s, 2s, 2p, and 3d orbitals. Dashed pink circles indicate radial nodes; dashed yellow or green lines indicate angular nodal planes. The probability density is highest near the nucleus for s orbitals and concentrated along lobes for p and d orbitals.

In the diagram, the 1s orbital is a simple sphere with no nodes; its probability density peaks at the nucleus and decays exponentially. The 2s orbital retains spherical symmetry but introduces one radial node—a spherical shell where the wavefunction changes sign and ψ = 0. The 2p orbital breaks spherical symmetry with a nodal plane through the nucleus, producing its characteristic dumbbell shape oriented along one Cartesian axis. Finally, the 3d orbital (shown here as dxy) features two nodal planes and a four-lobed cloverleaf pattern. Recognizing these patterns is essential for predicting orbital overlap and bond directionality in molecular orbital theory.

Mathematical Framework

The hydrogen atom Hamiltonian in atomic units (ℏ = me = e = 4πε₀ = 1) is Ĥ = −½∇² − Z/r, and the time-independent Schrödinger equation Ĥψ = Eψ is exactly solvable in spherical coordinates by writing ψ(r, θ, φ) = Rnl(r) × Ylm(θ, φ). The separation of variables yields a radial equation whose solutions involve associated Laguerre polynomials and an angular equation solved by spherical harmonics. The resulting energy eigenvalues depend only on n for hydrogen, but in multi-electron atoms the degeneracy in l is lifted by electron–electron interactions.

HYDROGEN-LIKE ENERGY EIGENVALUES
Eₙ = −Z² / (2n²) [atomic units] or Eₙ = −Z² × 13.6 eV / n²
Where Z is the nuclear charge, n is the principal quantum number (1, 2, 3, …), and 13.6 eV is the Rydberg energy (half a Hartree). For hydrogen, Z = 1 and E₁ = −13.6 eV.
RADIAL WAVEFUNCTION (GENERAL FORM)
Rₙₗ(r) = Nₙₗ × (2Zr / na₀)ˡ × e^(−Zr / na₀) × L^(2l+1)_(n−l−1)(2Zr / na₀)
Here Nnl is a normalization constant, a0 = 0.529 Å is the Bohr radius, and L denotes an associated Laguerre polynomial of degree n − l − 1. The factor (2Zr/na₀)ˡ ensures the wavefunction vanishes at the origin for l > 0.
RADIAL PROBABILITY DISTRIBUTION
P(r) dr = r² |Rₙₗ(r)|² dr
The factor r² arises from the volume element in spherical coordinates (dV = r² sin θ dr dθ dφ). The most probable radius for the 1s orbital of hydrogen is rmp = a₀, recovered by setting dP/dr = 0.
ANGULAR WAVEFUNCTIONS — SPHERICAL HARMONICS
Yₗᵐ(θ, φ) = Θₗₘ(θ) × Φₘ(φ) with Φₘ(φ) = (1/√2π) e^(imφ)
The polar part Θ involves associated Legendre polynomials. Real combinations (e.g., px ∝ sin θ cos φ, py ∝ sin θ sin φ) are used in chemistry because they align along Cartesian axes and facilitate visualization of bonding directionality.
⚛️ Multi-Electron Degeneracy Lifting
In hydrogen, orbitals with the same n are degenerate regardless of l. In multi-electron atoms, electron–electron repulsion and the resulting shielding break this degeneracy: within a given shell, lower-l orbitals (which penetrate closer to the nucleus) are stabilized relative to higher-l orbitals. This is the physical origin of the (n + l) rule used in the Aufbau ordering.

Electron Configuration — Aufbau Ordering & Notation

The ground-state electron configuration of an atom is constructed by placing electrons into orbitals in order of increasing energy, following the Aufbau principle, the Pauli exclusion principle, and Hund's rule. In the commonly used (n + l) rule (also called the Madelung rule), orbitals are filled in order of increasing n + l; for two subshells with the same n + l, the one with the smaller n fills first. This ordering is approximate—notable exceptions occur at Cr (Z = 24) and Cu (Z = 29) in the first transition series, where half-filled or fully filled d subshells confer extra stability through exchange energy. The diagram below illustrates the Aufbau filling order as an energy-level diagram.

Aufbau energy-level diagram for multi-electron atoms. Each horizontal line represents one orbital (capable of holding two electrons). Subshells are labeled on the right; the (n + l) value governs the filling order. Note that 4s fills before 3d because both have n + l = 4, but 4s has the lower n.
Aufbau filling order with quantum number details
Subshellnln + lm_l valuesMax electrons
1s10102
2s20202
2p213−1, 0, +16
3s30302
3p314−1, 0, +16
4s40402
3d325−2, −1, 0, +1, +210
4p415−1, 0, +16
⚠️ Exceptions to the Madelung Rule
Chromium (Z = 24) adopts the configuration [Ar] 3d⁵ 4s¹ instead of the expected [Ar] 3d⁴ 4s², and copper (Z = 29) is [Ar] 3d¹⁰ 4s¹ rather than [Ar] 3d⁹ 4s². These anomalies reflect the favorable exchange energy gained when a subshell is exactly half-filled or fully filled—a quantum-mechanical effect not captured by simple orbital-energy ordering. Similar exceptions occur in heavier transition metals and the lanthanides/actinides.

Worked Example — Electron Configuration of Iron (Z = 26)

Let us construct the ground-state electron configuration of iron (Z = 26) step by step, applying the Aufbau principle, Pauli exclusion, and Hund's rule. We will also identify the term symbol quantum numbers for the ground state.

Ground-State Configuration of Fe (Z = 26)
1
Step 1 — Count the electronsIron is a neutral atom with Z = 26, so there are 26 electrons to place. We fill orbitals in the Aufbau order: 1s, 2s, 2p, 3s, 3p, 4s, 3d, …
2
Step 2 — Fill the core shells (1s through 3p)1s² accounts for 2 electrons, 2s² for 2 more (total 4), 2p⁶ for 6 more (total 10), 3s² for 2 more (total 12), and 3p⁶ for 6 more (total 18). These 18 electrons constitute the argon core, [Ar].
[Ar] = 1s² 2s² 2p⁶ 3s² 3p⁶ — 18 electrons placed
3
Step 3 — Fill 4s before 3dThe 4s subshell (n + l = 4) fills before 3d (n + l = 5) in the Aufbau scheme. Two electrons enter 4s with opposite spins, bringing the total to 20.
[Ar] 4s² — 20 electrons placed
4
Step 4 — Place remaining electrons in 3d26 − 20 = 6 electrons remain. The 3d subshell has five orbitals (ml = −2, −1, 0, +1, +2). By Hund's rule, we place one electron in each of the five orbitals with parallel spins (↑ in each), consuming 5 electrons. The sixth electron pairs with one of the existing electrons (say ml = −2), giving ↑↓ in one orbital and ↑ in the remaining four.
[Ar] 3d⁶ 4s² — all 26 electrons placed
5
Step 5 — Determine the ground-state term symbolThe 3d⁶ configuration has four unpaired electrons, so S = 4 × ½ = 2, giving 2S + 1 = 5. The total orbital angular momentum L = (−1) + 0 + 1 + 2 = 2, corresponding to a D term. Since the subshell is more than half-filled, J = L + S = 4. The ground-state Russell–Saunders term symbol is therefore ⁵D₄.
Fe ground state: [Ar] 3d⁶ 4s², term symbol ⁵D₄
💡 Why write 3d before 4s?
Although 4s fills before 3d, spectroscopic convention lists subshells in order of n (so 3d before 4s) to emphasize the spatial proximity to the nucleus. When ionizing transition metals, 4s electrons are removed first because, in the cation, 3d is more tightly bound than 4s. Thus Fe²⁺ is [Ar] 3d⁶, not [Ar] 3d⁴ 4s².

Strengths & Limitations of the Orbital Approximation

The orbital model—assigning each electron to a single-particle wavefunction—is an approximation. It works remarkably well for predicting ground-state configurations, ionization energies, and periodic trends, but it has inherent limitations tied to the neglect of instantaneous electron–electron correlation. Understanding where the model succeeds and fails is crucial for knowing when to invoke more sophisticated methods.

Strengths and limitations of the orbital approximation
AspectStrengthsLimitations
Periodic trendsCorrectly predicts ionization energy, electron affinity, and atomic radius trends across periods and groups.Fails to quantitatively reproduce ionization energies without empirical shielding corrections (e.g., Slater rules).
SpectroscopyExplains gross features of atomic emission/absorption spectra and term symbol labeling.Cannot explain fine structure (spin–orbit coupling) or Lamb shift without relativistic and QED corrections.
Multi-electron atomsHartree–Fock provides reasonable total energies (≈ 99% of exact) and qualitative orbital pictures.Misses electron correlation energy (~1 eV per electron pair), critical for chemical accuracy.
Chemical bondingServes as the foundation for molecular orbital theory and LCAO approximations.Static correlation in bond-breaking processes requires multi-reference methods (CASSCF, MRCI).
Configuration exceptionsThe Madelung rule covers > 80% of elements correctly.Breaks down for Cr, Cu, Mo, Pd, and many lanthanides/actinides; requires explicit exchange-energy analysis.
KEY TAKEAWAY
The orbital picture is like a map projection: it faithfully represents the topology of electronic structure (which orbitals fill, how they order) even though it distorts the fine metric details (correlation energies, multiplet splittings). Just as a Mercator projection is indispensable for navigation despite its area distortion, the orbital approximation remains the essential starting point for understanding atomic and molecular electronic structure—one that post-Hartree–Fock methods refine rather than replace.

Connection to Advanced Theory

While the single-determinant orbital picture captures most of atomic structure, modern computational chemistry systematically improves upon it. The Hartree–Fock (HF) method replaces the intuitive Aufbau ordering with a self-consistent variational procedure that optimizes orbital shapes for a given electron configuration. Beyond HF, post-Hartree–Fock methods such as configuration interaction (CI), coupled cluster (CC), and multi-configurational self-consistent field (MCSCF) recover the missing correlation energy by mixing multiple Slater determinants. Density functional theory (DFT) offers an alternative route, replacing the many-electron wavefunction with the electron density as the fundamental variable while incorporating exchange-correlation effects through approximate functionals. The table below summarizes how the orbital picture connects to these more rigorous treatments.

Orbital model versus advanced quantum-chemical methods
FeatureOrbital / Aufbau ModelAdvanced QM Methods
Wavefunction formSingle Slater determinant built from hydrogen-like orbitalsLinear combination of many determinants (CI) or cluster expansion (CC)
Electron interactionMean-field approximation; each electron sees averaged potentialExplicit electron correlation included perturbatively or variationally
Orbital energiesApproximate; follow (n + l) ordering with empirical correctionsSelf-consistently optimized; Koopmans' theorem links to ionization energies
AccuracyQualitative; errors of ~1–5 eV in total energiesChemical accuracy (~1 kcal/mol) achievable with CCSD(T) or high-level DFT
Computational costNegligible (pen-and-paper)Scales as N⁴ (HF) to N⁷ (CCSD(T)); DFT scales as N³–N⁴

Looking ahead, the concepts developed in this lesson—quantum numbers, orbital symmetries, and configuration rules—carry directly into molecular orbital theory, where atomic orbitals combine via the LCAO approximation to form bonding and antibonding molecular orbitals. Understanding atomic electron configurations is also prerequisite for term symbols and selection rules that govern atomic spectroscopy. The quantum-mechanical framing established here—wavefunctions, probability densities, quantum numbers, and variational optimization—is the conceptual bedrock upon which all of computational and theoretical chemistry is built.

Practice Problems

PROBLEM 1CONCEPTUAL
Explain why the 2s orbital in a multi-electron atom is lower in energy than the 2p orbital, even though both have the same principal quantum number n = 2. Your answer should reference the concepts of penetration and shielding.
PROBLEM 2BASIC CALCULATION
Calculate the energy of an electron in the n = 3 level of the Li²⁺ ion (Z = 3) using the hydrogen-like energy formula Eₙ = −Z² × 13.6 eV / n². How does this compare to the n = 1 level of hydrogen?
PROBLEM 3INTERMEDIATE
Write the full ground-state electron configuration and the Russell–Saunders term symbol for vanadium (V, Z = 23). How many unpaired electrons does it have, and what is the predicted spin-only magnetic moment μ = √(n(n+2)) μ_B?
PROBLEM 4APPLIED
The radial probability distribution for a hydrogen 2s orbital is P(r) = r² |R₂₀(r)|², where R₂₀(r) = (1/4√2)(1/a₀)^(3/2)(2 − r/a₀)e^(−r/2a₀). Find the radial distance at which the radial node occurs, and verify that P(r) = 0 at that point.
PROBLEM 5CRITICAL THINKING
Copper (Z = 29) has the ground-state configuration [Ar] 3d¹⁰ 4s¹ rather than the Aufbau-predicted [Ar] 3d⁹ 4s². Provide a quantum-mechanical argument based on exchange energy and electron–electron repulsion that rationalizes this exception. Why doesn't nickel (Z = 28) similarly adopt [Ar] 3d¹⁰ 4s⁰?

Lesson Summary

This lesson developed the quantum-mechanical description of atomic structure from first principles. We began with the historical progression from Bohr's quantized orbits to Schrödinger's wave equation, whose solutions define atomic orbitals—one-electron wavefunctions labeled by three quantum numbers (n, l, ml). The shapes of s, p, and d orbitals were visualized through their radial and angular node patterns, governed by the counting rule: n − l − 1 radial nodes and l angular nodes. The radial probability distribution P(r) = r²|R(r)|² determines the most probable electron–nucleus distance and underpins concepts of penetration and shielding in multi-electron atoms.

Ground-state electron configurations are constructed using three rules: the Aufbau principle (fill in order of increasing energy, guided by the n + l rule), the Pauli exclusion principle (at most two electrons per orbital, with opposite spins), and Hund's rule (maximize spin within degenerate subshells). Notable exceptions such as Cr and Cu arise from exchange-energy stabilization of half-filled and fully filled d subshells. The orbital approximation serves as the indispensable foundation upon which Hartree–Fock theory, post-HF methods, and molecular orbital theory are constructed.

Varsity Tutors • Physical Chemistry 2 • Atomic Orbitals & Electron Configurations