← SciSim / Chemistry
Hydrogen Atom Orbitals (s, p, d) Topic #15

§1 Interactive Simulation

Three.js r128
Orbital
1s
E (eV)
-13.60
n, l
1, 0
Radial nodes
0
Angular nodes
0
⟨r⟩ (Å)
0.79
⚙ Controls
n (principal)
114
l (angular, 0=s 1=p 2=d)
003
Rotation angle
0°0°360°
Rotation Speed
0×1×4×
Phase coloring (±)
Show nodal planes
Show axes

§2 The Idea, Step by Step

From a fuzzy cloud to a precise shape

You can never pin down exactly where an atom's electron is — only where it is likely to be. Picture a swarm of bees around a hive: no single bee sits in a fixed spot, yet the swarm as a whole has a definite shape. An orbital is that shape: the fuzzy cloud showing where the electron probably is. The cloud above is exactly that probability map for a single electron around a hydrogen nucleus.

Just three whole numbers — the quantum numbers — fix the cloud. The principal number $n$ (1, 2, 3, …) sets the cloud's size and energy: bigger $n$ means a larger, higher-energy cloud. The angular number $l$ (0, 1, 2, called s, p, d) sets the shape: $l=0$ is a ball (s), $l=1$ a dumbbell (p), $l=2$ a cloverleaf (d). For hydrogen the energy depends only on $n$: $E_n = -\dfrac{13.6}{n^2}\text{ eV}$. Worked number: the ground state $n=1$ gives $E_1 = -13.6$ eV — exactly the energy it takes to tear hydrogen's electron away (its ionization energy). Step up to $n=2$ and $E_2 = -13.6/4 = -3.40$ eV, much closer to zero and far easier to remove.

At the precise level the full wavefunction splits into a radial part and an angular part, $\psi_{n,l,m}(r,\theta,\phi) = R_{n,l}(r)\,Y_l^m(\theta,\phi)$. The radial part $R_{n,l}$ controls how the density falls off with distance and how many spherical shells it has — the count of radial nodes is $n-l-1$. The angular part $Y_l^m$ sets the shape, contributing $l$ angular nodes (flat planes or cones of zero probability), so the total node count is always $n-1$. What you actually see drawn as the cloud is $|\psi|^2$, the probability density. In the sim, the $n$ slider grows the cloud and lifts the energy bars; the $l$ slider morphs ball → dumbbell → cloverleaf and adds nodal planes; the rotation slider spins the 3D isosurface so you can confirm it is a real solid shape, not a flat drawing.

Try this in the sim above: (1) Set $n=2$, $l=0$ and switch to the 4πr²R² graph — find the gap at $r=2a_0$ where the 2s probability drops to zero, a radial node you can see. (2) Keep $n=2$ and slide $l$ from 0 to 1: watch the ball become a two-lobed dumbbell with a nodal plane slicing through the nucleus. (3) Open the E levels graph and raise $n$ from 1 to 4 — the rungs crowd toward zero, showing why high-$n$ levels lie almost on top of one another.

§3 Equation Derivation

Hydrogen Atom Wavefunction — Separation of Variables
\[\psi_{n,l,m}(r,\theta,\phi) = R_{n,l}(r)\cdot Y_l^m(\theta,\phi)\] \[E_n = -\frac{13.6\,Z^2}{n^2}\text{ eV}\]
SymbolMeaningValues
\(n\)Principal quantum number1, 2, 3, …
\(l\)Angular momentum quantum number0 to n−1
\(m\)Magnetic quantum number−l to +l
\(R_{n,l}(r)\)Radial wavefunctionLaguerre polynomials × e^{−r/na₀}
\(Y_l^m(\theta,\phi)\)Spherical harmonic (shape)gives orbital shape
\(a_0\)Bohr radius = 52.9 pmm
Step-by-Step Derivation
Step 1 — Schrödinger equation in spherical coordinates \(-\dfrac{\hbar^2}{2m_e}\nabla^2\psi - \dfrac{e^2}{4\pi\varepsilon_0 r}\psi = E\psi\) In spherical coords, \(\nabla^2 = \frac{1}{r^2}\frac{\partial}{\partial r}\!\left(r^2\frac{\partial}{\partial r}\right) + \frac{1}{r^2\sin\theta}\frac{\partial}{\partial\theta}\!\left(\sin\theta\frac{\partial}{\partial\theta}\right) + \frac{1}{r^2\sin^2\theta}\frac{\partial^2}{\partial\phi^2}\)
Step 2 — Separation: ψ = R(r)Y(θ,φ) The equation separates into a radial equation (for R) and an angular equation (for Y). The angular equation gives spherical harmonics \(Y_l^m\) with eigenvalue \(l(l+1)\).
Step 3 — Radial wavefunctions (key examples) \(R_{10} = 2a_0^{-3/2}e^{-r/a_0}\) (1s, 0 radial nodes)
\(R_{20} = \frac{1}{2\sqrt{2}}a_0^{-3/2}(2-r/a_0)e^{-r/2a_0}\) (2s, node at r = 2a₀)
\(R_{21} = \frac{1}{2\sqrt{6}}a_0^{-3/2}(r/a_0)e^{-r/2a_0}\) (2p, no radial nodes)
Step 4 — Angular wavefunctions (shapes) \(Y_0^0 = \frac{1}{2\sqrt{\pi}}\) (s: spherical)
\(Y_1^0 = \sqrt{\frac{3}{4\pi}}\cos\theta\) (p_z: dumbbell along z)
\(Y_2^0 = \sqrt{\frac{5}{16\pi}}(3\cos^2\theta-1)\) (d_{z²}: lobe + torus)
Step 5 — Node counting Radial nodes = n − l − 1 · Angular nodes (nodal planes/cones) = l · Total nodes = n − 1
Example: 3d (n=3, l=2): radial nodes = 0, angular nodes = 2, total = 2
Step 6 — Radial distribution function \(P(r) = 4\pi r^2|R_{n,l}|^2\) — probability per unit radial distance
Peak of 1s: \(r_{max} = a_0\). Mean radius: \(\langle r\rangle_{ns} = \frac{a_0}{2}(3n^2-l(l+1))\)
Worked Example — Radial Node of 2s orbital

R₂₀(r) = 0 when \(2 - r/a_0 = 0\) → r = 2a₀ = 2 × 52.9 pm = 105.8 pm

At this radius, the 2s electron has zero probability of being found. Inside this shell: electron density peaks near r ≈ a₀/2. Outside: second maximum near r ≈ 5a₀.

Reference: Atkins & de Paula — Physical Chemistry, 11th Ed., §9A "Hydrogenic Atoms" | McQuarrie & Simon — Physical Chemistry: A Molecular Approach, Chapter 6

§4 Frequently Asked Questions

🔬 SimulationWhat does each simulation mode show?▼
The s Orbitals tab shows 1s/2s/3s as 2D cross-sections: spherically symmetric electron density clouds, with the 2s and 3s showing concentric nodal shells. p Orbitals shows the dumbbell pz shape — two lobes with opposite phase (purple = positive, red = negative). d Orbitals displays all five d shapes: dz², dxz, dyz, dxy, dx²−y². All Subshells places all orbitals of a shell side by side. Dot Density renders random dots where density ∝ |ψ|², mimicking how electron density looks in real electron microscopy. Key: Shape comes from angular quantum number l; number of shells/layers comes from n.
🌍 Real LifeWhere do orbital shapes directly control chemistry?▼
Orbital shapes control every aspect of bonding geometry. sp³ hybridization gives methane's 109.5° tetrahedral shape; sp² gives benzene's flat hexagon with π electrons above and below; sp gives acetylene's linear structure. Transition metal catalysis — used in pharmaceutical synthesis, polymer production, and petroleum refining — depends critically on d orbital geometry determining how substrate molecules approach the metal. The colors of gems (ruby = Cr d-d transitions, sapphire = Fe/Ti d-d transitions) arise directly from crystal field splitting of d orbitals. Key: Every molecular shape, chemical bond, and spectroscopic color ultimately traces back to atomic orbital geometry.
🧪 ConceptualWhat physically is the difference between ψ and |ψ|²?▼
ψ (the wavefunction) is a complex mathematical amplitude — it can be positive, negative, or complex. It is not directly observable. |ψ|² is the probability density: the probability of finding the electron in volume dV is |ψ|² dV. This is what experiments measure. The sign of ψ matters for bonding: when two atomic ψ functions overlap with the same sign, they interfere constructively (bonding MO, increased electron density between nuclei). Opposite signs give destructive interference (antibonding MO, node between nuclei). Key: ψ encodes phase (±) needed for bonding; |ψ|² is the observable probability — both are essential in chemistry.
🧮 MathematicalHow many orbitals are in each shell and subshell?▼
For quantum number l, m runs from −l to +l giving 2l+1 orbitals: s (l=0) → 1 orbital; p (l=1) → 3 orbitals; d (l=2) → 5 orbitals; f (l=3) → 7 orbitals. Total orbitals in shell n = n² (since Σ(2l+1) for l=0 to n−1 = n²). With spin (ms = ±½), capacity per shell = 2n². This gives: n=1: 2e, n=2: 8e, n=3: 18e — the electron capacity pattern of the periodic table. Example: 4th period has s(2) + p(6) + d(10) = 18 elements. Key: Orbital count per subshell = 2l+1; total per shell = n² — this directly explains the periodic table's row lengths.
💡 Non-ObviousWhy does 2s penetrate closer to nucleus than 2p, and why does this matter?▼
The 2s radial wavefunction has a small inner peak near r ≈ 0.6a₀ (before its radial node at 2a₀), while 2p has R₂₁ ∝ r, going to zero at the nucleus. This "penetration" of 2s into the inner core means 2s electrons spend some time very close to the nucleus, experiencing nearly the full nuclear charge Z. This lowers 2s energy below 2p in multi-electron atoms (for Li: 2s is at −5.39 eV, 2p at −3.54 eV). This penetration and shielding effect is why 4s fills before 3d in the Aufbau principle. Key: Orbital penetration (s > p > d > f) creates energy differences between subshells — this is what drives the Aufbau order.
🌌 DeepWhy are there exactly 5 d orbitals and not 4 or 6?▼
For l=2, the magnetic quantum number m takes values −2, −1, 0, +1, +2 — exactly 5 values (2l+1=5). This is a mathematical requirement: these are the only values for which the spherical harmonic Y₂ᵐ(θ,φ) is single-valued and normalizable. The 5 real d orbitals (dxy, dxz, dyz, dx²−y², dz²) are linear combinations of the complex m = ±1, ±2, 0 solutions. Note that dz² looks different from the other four — it's actually a linear combination of m=0 and a "hidden" combination that keeps the count at 5, not 6. Key: The 2l+1 rule is a mathematical inevitability from solving spherical harmonics — it cannot be any other number.

Reference: LibreTexts Chemistry — Atomic Orbitals https://chem.libretexts.org/Bookshelves/Physical_and_Theoretical_Chemistry_Textbook_Maps/Supplemental_Modules/Quantum_Mechanics | MIT OCW 5.61 Physical Chemistry, Lectures 10-12

§5 Common Misconceptions

❌ Misconception: "The two lobes of a p orbital mean the electron goes back and forth between them."
✅ Correction: The electron exists simultaneously throughout both lobes as a quantum mechanical standing wave. The sign difference (+ and −) between lobes represents the phase of ψ, not electron location. The probability density |ψ|² is identical in both lobes. The electron has zero probability of being in the nodal plane at the center — but this does not mean it "crosses" from one lobe to another; such a trajectory picture is meaningless in quantum mechanics.
📖 Reference: Atkins & de Paula — Physical Chemistry, 11th Ed., §9A.2: "The interpretation of the wavefunctions"
❌ Misconception: "The boundary surface drawn for an orbital is a hard wall — the electron cannot go outside it."
✅ Correction: The 90% contour surface conventionally shown for orbitals means there is 90% probability of finding the electron within that surface — 10% of the time the electron is outside it. The wavefunction decays exponentially but extends to infinity. This is not merely theoretical: quantum tunneling (essential in nuclear fusion, enzyme catalysis, and tunnel diodes) occurs precisely because electrons can exist in classically "forbidden" regions outside any classical boundary.
📖 Reference: Silberberg — Chemistry: The Molecular Nature, 9th Ed., Chapter 7.4: "The Quantum-Mechanical Model"
❌ Misconception: "The d orbital notation like dz² means the orbital is only along the z-axis — it's a line, not a 3D shape."
✅ Correction: The subscript z² refers to the angular function (3cos²θ−1), which depends on the angle from z — it creates a 3D shape with two large lobes along z plus a toroidal ring (doughnut) in the xy plane. All orbitals are three-dimensional probability distributions. The notation dxy, dxz, dyz, dx²−y², dz² describes the orientation of the lobes in space, not that the electron is confined to those axes.
📖 Reference: Housecroft & Sharpe — Inorganic Chemistry, 5th Ed., Chapter 1.7: "Atomic orbitals"
❌ Misconception: "3d and 4s orbitals have the same energy in all atoms because they fill together."
✅ Correction: In hydrogen, all orbitals with the same n have the same energy. In multi-electron atoms, 4s is lower than 3d for the neutral atoms K and Ca (which is why 4s fills first in Aufbau). However, once transition metals begin filling 3d, the 3d orbitals drop below 4s due to increased nuclear charge. This is why ions of transition metals (e.g., Fe²⁺, Fe³⁺) lose 4s electrons first — at that nuclear charge, 4s > 3d in energy.
📖 Reference: Housecroft & Sharpe — Inorganic Chemistry, 5th Ed., Chapter 1.9: "Many-electron atoms" — 3d/4s energy crossover explicitly addressed
❌ Misconception: "s orbitals are simple uniform spheres of constant electron density."
✅ Correction: s orbitals are spherically symmetric in angle but not uniform with distance r. For 1s: |ψ|² ∝ e^{−2r/a₀}, peaking at the nucleus and decaying outward. The radial distribution function 4πr²|R|² peaks at r = a₀ (not at the nucleus) because the volume element 4πr² grows with r. For 2s and 3s, |ψ|² shows radial nodes — shells of zero probability — making them multi-shell structures, not uniform spheres.
📖 Reference: Atkins & de Paula — Physical Chemistry, 11th Ed., Figure 9A.6: "Radial distribution functions" for 1s, 2s, 3s

Section 4 reference: Taber, K.S. — Chemical Misconceptions (RSC, 2002) | Tsaparlis, G. — "Atomic orbitals: a survey of misconceptions" J. Chem. Educ. 2001, 78, 1432 | Barke, H-D. et al. — Misconceptions in Chemistry (Springer, 2009) Ch. 4