Penetration and Shielding

Sodium's outermost electron feels a pull from a nucleus with eleven protons. Yet it ionizes almost as easily as a hydrogen atom held by a single proton. Its first ionization energy is only 496 kJ/mol. Hydrogen's lone electron, in a much smaller shell, has an ionization energy of 1312 kJ/mol. The resolution is not a change in nuclear charge. Instead, it is a change in what the outer electron actually experiences. Ten inner electrons stand between it and the nucleus. They cancel nearly all of that extra charge. This is shielding. It is one of two ideas — the other is penetration — that explain almost every irregularity. These irregularities appear in orbital energies, ionization trends, and the entire structure of the periodic table. Penetration and shielding matter for a key reason. The simple picture from the hydrogen atom, where energy depends only on the principal quantum number nn, collapses instantly when a second electron is added. In hydrogen, 2s2s and 2p2p are exactly degenerate. In every other atom, 2s2s lies measurably below 2p2p. Also, 3d3d can lie above 4s4s depending on occupation. This ll-splitting within a shell generates the block structure of the periodic table. It explains why the transition metals appear where they do. It explains why lanthanide contraction happens. It also explains why the aufbau order is 4s4s before 3d3d, but 3d3d before 4s4s in ionization. None of this follows from Coulomb's law alone. Instead, it follows from the shapes of radial distribution functions. These shapes are built on Radial Distribution Functions and Nodes and Orbital Shapes and Boundary Surfaces. The big picture is this: every electron in a multielectron atom moves in an averaged, spherically symmetric potential. This potential is produced by the nucleus plus all the other electrons. An electron in an orbital with substantial probability density close to the nucleus is called a penetrating orbital. It spends part of its time inside the charge cloud of the other electrons. There, it feels more of the bare nuclear charge and less screening. Orbitals of the same nn but different ll penetrate to different extents. This happens because their radial distribution functions have different numbers of nodes. They also have different near-nucleus behavior. The result is an effective nuclear charge ZeffZ_{\text{eff}} that depends on both nn and ll, not just nn. This article builds the physical picture from radial probability distributions. It introduces Slater's rules as a quantitative, if approximate, recipe for ZeffZ_{\text{eff}}. It works through the energetics of ll-splitting with real numbers for sodium, chlorine, and a transition-metal ion. It closes by connecting penetration and shielding to the topics they set up. These are Effective Nuclear Charge as a systematic quantity, and the full Orbital Energy Ordering in Multielectron Atoms that governs the aufbau principle.

Key takeaways
  • Inner electrons shield outer electrons from the full nuclear charge, defining an effective nuclear charge Zeff<ZZ_{\text{eff}} < Z.
  • Penetration is the tendency of an orbital's radial density to have a peak or "spike" close to the nucleus, inside the shielding cloud of other electrons.
  • For a given nn, penetration decreases in the order s>p>d>fs > p > d > f, because ss orbitals have the most nodes and the largest near-origin density.
  • Greater penetration means less effective shielding felt by that electron, so it is bound more tightly — this splits the ll-degeneracy that exists in hydrogen.
  • Slater's rules give a working, hand-calculable estimate of ZeffZ_{\text{eff}} using empirical shielding constants σ\sigma grouped by orbital type.
  • Penetration/shielding together explain ionization energy trends, the 4s4s/3d3d energy crossover, and lanthanide contraction.
  • Key idea: electrons in more penetrating orbitals sample more of the bare nucleus and less of the other electrons' charge, so orbital energy within a shell splits as s<p<d<fs < p < d < f, purely because of orbital shape, not principal quantum number.

Where this lives today

Penetration and shielding are not a historical curiosity; they underpin working explanations across chemistry. In spectroscopy, X-ray photoelectron spectroscopy (XPS) resolves core-level binding energies that shift systematically with oxidation state precisely because changing the number of valence electrons changes the shielding felt by core electrons — the "chemical shift" in XPS is a direct, quantitative signature of penetration/shielding physics, used routinely to determine oxidation states of metals in catalysts and battery electrode materials. In coordination chemistry, the fact that 3d3d orbitals barely penetrate the argon-like core (as shown numerically above) is why they are strongly perturbed by ligand fields yet remain compact and localized — the entire crystal-field and ligand-field treatment of transition-metal color and magnetism assumes exactly this contracted, non-penetrating character of dd orbitals. Lanthanide contraction — the anomalously small ionic radii of the third-row transition metals following the lanthanides — is a direct penetration/shielding effect: poorly-shielding 4f4f electrons let the nucleus pull inward on the 5d5d/6s6s shells more than expected, with consequences for gold's chemistry, the density of tungsten, and the near-identical radii of zirconium and hafnium that famously complicates their separation in nuclear reactor materials. In drug design and bioinorganic chemistry, the effective charge felt by metal-bound electrons (understood via penetration arguments) governs Lewis acidity trends used to predict how zinc, iron, and copper centers in metalloenzymes activate substrates. Finally, every periodic-table trend taught at the introductory level — ionization energy, atomic radius, electron affinity, electronegativity — is, underneath the memorized arrows, a story about how penetration-dependent shielding reshapes ZeffZ_{\text{eff}} across rows and down columns, the systematic version of which is developed under Effective Nuclear Charge.

Ask the tutor: Why does the 3d orbital drop below 4s in energy for transition-metal cations even though Slater's rules treat the core screening of 3d as essentially unchanged from the neutral atom?

The physical picture

Start from the one-electron atom, where the Schrödinger equation gives energies depending only on nn:

En=−Z2n2 RyE_n = -\frac{Z^2}{n^2}\,\text{Ry}

For hydrogen, 2s2s and 2p2p have identical energy. Add a second electron — helium, or better, think of sodium's valence electron moving in the field of a Na+\mathrm{Na}^{+} core — and the potential the outer electron feels is no longer a clean −Ze2/r-Ze^2/r. Close to the nucleus, where inner-shell electron density is sparse, the outer electron sees nearly the full nuclear charge ZZ. Far from the nucleus, the inner electrons' charge cloud sits between the outer electron and the nucleus, and by Gauss's law that inner charge cancels part of the nuclear attraction, so the outer electron sees an effective charge close to Z−(number of inner electrons)Z - (\text{number of inner electrons}). The crucial question is: how much time does an orbital spend in the "sees full ZZ" region versus the "screened" region? That is exactly what the radial distribution function 4πr2Rnl2(r)4\pi r^2 R_{n l}^2(r) answers. Recall that nsns orbitals have n−1n-1 radial nodes, npnp have n−2n-2, ndnd have n−3n-3, and so on — for fixed nn, lower ll means more radial nodes and (crucially) a larger-amplitude inner lobe close to r=0r=0. An electron in 3s3s has real, non-negligible probability density very near the nucleus (inside the 1s1s and 2s,2p2s,2p shells), while a 3p3p electron has a smaller inner lobe, and a 3d3d electron — with no radial nodes at all for n=3n=3, just the pure rlr^{l} falloff at small rr — has essentially none. This is penetration made quantitative: it is literally area under the radial distribution curve at small rr, inside the core.

Figure
Radial distribution functions for n = 3: penetration order 3s > 3p > 3d
r (distance from nucleus) 4πr²R²(r) core region 3s (2 inner peaks + outer) 3p (1 inner peak + outer) 3d (no inner peak) outer maxima nearly coincide near same r

The figure makes the essential point visually: all three orbitals (3s3s, 3p3p, 3d3d) have their main, outermost lobe at roughly similar radius (they belong to the same shell), but only 3s3s and 3p3p have extra inner lobes that dip deep into the core region shaded on the left. The 3s3s curve has two small inner humps (from its two inner radial nodes) plus the outer lobe; 3p3p has one inner hump; 3d3d has none. That inner probability is what "penetration" means physically, and it directly controls how much nuclear charge is screened away. Now connect this to energy. Define, for an electron in orbital ψnl\psi_{nl}, the effective nuclear charge

Zeff(n,l)=Z−σnlZ_{\text{eff}}(n,l) = Z - \sigma_{nl}

where σnl\sigma_{nl} is the shielding constant, the amount of nuclear charge cancelled by the other electrons as experienced by this particular orbital. Because σnl\sigma_{nl} depends on how much the electron penetrates the core, σ3s<σ3p<σ3d\sigma_{3s} < \sigma_{3p} < \sigma_{3d}, so ZeffZ_{\text{eff}} decreases in that same order, and by the hydrogen-like energy formula (now with Z→ZeffZ\to Z_{\text{eff}}),

Enl≈−(Zeff(n,l)n)2RyE_{nl} \approx -\left(\frac{Z_{\text{eff}}(n,l)}{n}\right)^{2}\text{Ry}

a larger ZeffZ_{\text{eff}} means a more negative (more stable) energy. Hence within the same nn: Ens<Enp<End<EnfE_{ns} < E_{np} < E_{nd} < E_{nf}. This is the ll-splitting the hydrogen atom does not show, and it is entirely a penetration/shielding effect.

Figure
Energy-level splitting within n = 3 and n = 4 due to penetration
Hydrogen (l-degenerate) vs. multielectron atom (l-split) n=3, all l (H atom) 3d 3p 3s Turning on electron–electron repulsion + penetration splits l

Note this is purely a shielding/penetration effect within one shell; it says nothing yet about cross-shell competition such as 4s4s vs. 3d3d, which involves comparing different nn and is treated with real numbers in the worked examples below and formally in Orbital Energy Ordering in Multielectron Atoms.

The core ideas

The governing relation is the modified hydrogenic energy formula:

Enl=−(Zeff(n,l)n∗)2×13.6 eVE_{nl} = -\left(\frac{Z_{\text{eff}}(n,l)}{n^{*}}\right)^{2} \times 13.6\ \text{eV}

where:

  • Zeff(n,l)=Z−σnlZ_{\text{eff}}(n,l) = Z - \sigma_{nl} is the effective nuclear charge for that specific orbital,
  • n∗n^{*} is the effective principal quantum number (equal to nn for n≤3n \le 3; reduced for n=4,5,6n = 4,5,6 in Slater's scheme),
  • σnl\sigma_{nl} is built from other electrons' contributions, weighted by how much they shield this orbital.

Reading: the formula says an electron's binding energy scales with the square of the charge it actually experiences, not the bare nuclear charge — this is why even a modest change in shielding (a change in σ\sigma of order 1) produces a large, chemically significant change in energy, since the effect is squared. The shielding constant itself, in Slater's rules, is built additively:

σnl=∑iki\sigma_{nl} = \sum_{i} k_i

where:

  • the sum runs over all other electrons in the atom,
  • kik_i is the individual screening contribution of electron ii, which depends on whether ii is in the same group, the shell just inside, or further inside, and on whether the orbital being screened is s/ps/p-type or d/fd/f-type.

Reading: this equation is simply "shielding is additive over the other electrons," but the values of kik_i encode all the penetration physics — an inner dd electron screens an outer ss electron much more effectively (nearly full k=1k=1) than an outer dd electron screens another dd electron in the same subshell (partial k≈0.35k \approx 0.35), because ss orbitals penetrate deep inside dd shells while same-subshell dd–dd interactions are mutually poorly screening. Finally, effective nuclear charge connects to observable ionization energy through

IE1≈(Zeffn∗)2×13.6 eVIE_1 \approx \left(\frac{Z_{\text{eff}}}{n^{*}}\right)^{2}\times 13.6\ \text{eV}

where:

  • IE1IE_1 is the first ionization energy (removing the outermost, highest-n∗n^* electron),
  • the same Zeff,n∗Z_{\text{eff}}, n^{*} as defined above apply to the electron being removed.

Reading: this is the empirical bridge between the abstract shielding constant and something measured in a mass spectrometer or photoelectron spectrometer — every trend in ionization energy across a period or down a group is, at bottom, a story about how ZeffZ_{\text{eff}} for the outermost electron changes.

Derivation AWhy s beats p beats d in penetration: near-nucleus behavior of radial wavefunctions

Goal: show, from the mathematical form of hydrogenic radial functions, that at small rr the penetration ranking s>p>d>fs > p > d > f is guaranteed by the power of rr with which each orbital vanishes at the origin.

Step 1 — General small-rr form. Every hydrogenic radial function behaves near the origin as

Rnl(r)→r→0Cnl r lR_{nl}(r) \xrightarrow{r\to 0} C_{nl}\, r^{\,l}

This comes directly from solving the radial equation: the centrifugal term l(l+1)/r2l(l+1)/r^2 in the effective potential forces Rnl∼rlR_{nl}\sim r^l as the only solution that stays finite (and normalizable) at r=0r=0. An ss orbital (l=0l=0) approaches a finite, nonzero value at the nucleus; a pp orbital (l=1l=1) vanishes linearly; a dd orbital (l=2l=2) vanishes quadratically; an ff orbital (l=3l=3) vanishes cubically.

Step 2 — Radial probability density near the origin. The probability of finding the electron in a thin shell at radius rr is

P(r) dr=4πr2Rnl2(r) dr∼4πCnl2 r2l+2 drP(r)\,dr = 4\pi r^2 R_{nl}^2(r)\,dr \sim 4\pi C_{nl}^2\, r^{2l+2}\,dr

For ss: P(r)∼r2P(r)\sim r^2. For pp: P(r)∼r4P(r)\sim r^4. For dd: P(r)∼r6P(r)\sim r^6. For ff: P(r)∼r8P(r)\sim r^8. All of these vanish at r=0r=0 (as they must, since the volume element 4πr24\pi r^2 itself vanishes), but they vanish at very different rates. For small but nonzero rr inside the core, r2≫r4≫r6≫r8r^2 \gg r^4 \gg r^6 \gg r^8 whenever r<1r<1 (in atomic units, core radii are indeed sub-unity), so the ss electron has dramatically more density in the near-nucleus zone than pp, which has more than dd, which has more than ff.

Step 3 — Integrate over the core region. Define the penetration fraction as the probability found inside some core radius rcr_c (roughly the radius of the inner shells):

fpen(nl)=∫0rc4πr2Rnl2(r) drf_{\text{pen}}(nl) = \int_0^{r_c} 4\pi r^2 R_{nl}^2(r)\,dr

Given the Step 2 scaling, and the fact that the outer lobes of nsns, npnp, ndnd all sit at roughly comparable radius (same principal shell, so similar total extent), the extra inner lobes of low-ll orbitals contribute disproportionately to fpenf_{\text{pen}}. Explicit evaluation for hydrogenic n=3n=3 functions gives fpen(3s)>fpen(3p)>fpen(3d)f_{\text{pen}}(3s) > f_{\text{pen}}(3p) > f_{\text{pen}}(3d) by roughly an order of magnitude between successive ll, consistent with the radial-node count argument: 3s3s has 2 nodes (3 lobes total, 2 of them inner), 3p3p has 1 node (2 lobes, 1 inner), 3d3d has 0 nodes (1 lobe, no inner probability at all).

Step 4 — Translate penetration into shielding felt. Inside radius rcr_c, the inner electrons contribute little cancelling charge (Gauss's law: charge outside your radius does not shield you), so an electron sampling that region via fpenf_{\text{pen}} feels close to the bare ZZ. Averaging over the whole orbital, the effective shielding constant is reduced in proportion to how much probability sits inside rcr_c:

σnl≈σmax⁡[1−fpen(nl)]\sigma_{nl} \approx \sigma_{\max}\left[1 - f_{\text{pen}}(nl)\right]

where σmax⁡\sigma_{\max} would be the shielding felt by a completely non-penetrating point charge at the outer radius. Since fpen(ns)>fpen(np)>fpen(nd)f_{\text{pen}}(ns) > f_{\text{pen}}(np) > f_{\text{pen}}(nd), it follows immediately that

σns<σnp<σnd\sigma_{ns} < \sigma_{np} < \sigma_{nd}

which is the shielding-order result used throughout this article, derived here purely from the power-law behavior of hydrogenic radial functions at small rr — no empirical input required.

Limiting case check: as l→n−1l\to n-1 (the maximum possible angular momentum for given nn, e.g. 3d3d, 4f4f), the radial function has zero nodes and the smallest possible near-origin exponent is maximal (rlr^l with the largest ll), so fpen→f_{\text{pen}}\to its minimum value for that shell — consistent with dd and ff orbitals being the least penetrating, most hydrogen-like (highest symmetry, most "outer-shell-like") members of their shell.

Go deeper · Derivation B · 2 minSlater's rules: building Z_eff by hand
Derivation BSlater's rules: building Z_eff by hand

Goal: derive a working numerical estimate of ZeffZ_{\text{eff}} for a chosen electron using Slater's empirical grouping scheme, and show explicitly how it encodes the penetration ranking derived above.

Step 1 — Group the electrons. Write the electron configuration in Slater groups:

(1s)(2s,2p)(3s,3p)(3d)(4s,4p)(4d)(4f)(5s,5p)…(1s)(2s,2p)(3s,3p)(3d)(4s,4p)(4d)(4f)(5s,5p)\ldots

Each parenthesis is a group; nsns and npnp are grouped together (they screen each other with the same weight) while ndnd and nfnf are always separate groups from everything at the same nn. This asymmetry is itself an encoding of penetration: ns/npns/np penetrate similarly enough (both have some inner density) that Slater treats them as equivalent screeners of each other, whereas ndnd is different enough (no inner lobes) that it is kept separate.

Step 2 — Assign screening contributions kik_i. For an electron being screened in an nsns or npnp orbital:

  • each other electron in the same group contributes k=0.35k = 0.35 (0.30 if the group is 1s1s),
  • each electron in a group with principal quantum number one less (n−1n-1 shell) contributes k=0.85k = 0.85,
  • each electron in any group with principal quantum number two or more less contributes k=1.00k = 1.00 (complete screening).

For an electron being screened in an ndnd or nfnf orbital:

  • each other electron in the same group contributes k=0.35k = 0.35,
  • every electron in any lower group (lower nn, or same nn but lower ll) contributes the full k=1.00k = 1.00 — there is no partial-screening (0.85) tier for d/fd/f electrons, because dd and ff orbitals do not penetrate the core enough to justify anything less than full screening from everything inside them.

This asymmetry — d/fd/f electrons get fully screened by everything inside, with no partial credit — is the direct numerical fossil of the penetration argument in the first derivation: since d,fd,f orbitals barely penetrate, essentially all inner charge shields them completely.

Step 3 — Sum and subtract.

σ=∑iki,Zeff=Z−σ\sigma = \sum_i k_i, \qquad Z_{\text{eff}} = Z - \sigma

Step 4 — Worked instance (kept general here; full numeric example below). Consider comparing a 3p3p and a 3d3d electron in the same hypothetical atom with core 1s2 2s22p6 3s21s^2\,2s^2 2p^6\,3s^2 and one additional electron either in 3p3p or (hypothetically) 3d3d. For the 3p3p electron: same-group (3s,3p3s,3p) partners contribute 0.350.35 each; the n=2n=2 shell (2s,2p2s,2p, 8 electrons) contributes 0.850.85 each; the 1s1s shell contributes 1.001.00. For a 3d3d electron in the same atom: the 3s,3p3s,3p electrons (lower ll, same nn) now contribute the full 1.001.00 each (not 0.850.85), and so do the n=2n=2 and n=1n=1 shells. The 3d3d electron is screened more (larger σ\sigma, smaller ZeffZ_{\text{eff}}) purely because of the ll-dependent screening tier — exactly reproducing σ3d>σ3p\sigma_{3d} > \sigma_{3p} from the physical derivation.

Limiting case check: for hydrogen itself (Z=1Z=1, one electron, no other electrons to sum over), σ=0\sigma = 0 trivially for any orbital, so Zeff=1Z_{\text{eff}} = 1 regardless of ll — correctly recovering the exact hydrogenic ll-degeneracy Slater's rules are built to reduce to when there is nothing to screen.

Worked examples

Example 1 — Foundational: Z_eff for the sodium valence electron. Sodium, Z=11Z=11, configuration 1s2 2s22p6 3s11s^2\,2s^2 2p^6\,3s^1. We want ZeffZ_{\text{eff}} felt by the lone 3s3s electron. Group structure: (1s2)(2s22p6)(3s1)(1s^2)(2s^2 2p^6)(3s^1). Screeners of the 3s3s electron:

  • Same group (3s,3p3s,3p): no other electrons here (the group has only this one electron), contributes 00.
  • n=2n=2 shell (2s22p62s^2 2p^6, 8 electrons), one shell in: each contributes 0.850.85:
8×0.85=6.808 \times 0.85 = 6.80
  • n=1n=1 shell (1s21s^2, 2 electrons), two shells in: each contributes 1.001.00:
2×1.00=2.002 \times 1.00 = 2.00

Total shielding:

σ=6.80+2.00=8.80\sigma = 6.80 + 2.00 = 8.80
Zeff=11−8.80=2.20Z_{\text{eff}} = 11 - 8.80 = 2.20

Compare with the naive (wrong) picture that ten core electrons cancel ten units of charge, giving Zeff=1Z_{\text{eff}}=1: Slater's rules show the real number is over twice that, because the 1s1s electrons screen almost perfectly but the n=2n=2 shell screens only 85% as effectively — some 3s3s density penetrates inside the n=2n=2 shell. Using n∗=3n^*=3 (Slater's rule: n∗=nn^*=n for n≤3n\le 3):

IE1≈(2.203)2×13.6 eV=(0.733)2×13.6=0.5379×13.6=7.32 eVIE_1 \approx \left(\frac{2.20}{3}\right)^2 \times 13.6\ \text{eV} = (0.733)^2\times 13.6 = 0.5379\times13.6 = 7.32\ \text{eV}

Converting, 7.32 eV×96.485 kJ mol−1eV−1=706 kJ/mol7.32\ \text{eV} \times 96.485\ \text{kJ mol}^{-1}\text{eV}^{-1} = 706\ \text{kJ/mol}. The experimental value is 496 kJ/mol. Slater's rules overestimate here (a known limitation for alkali metals, where the simple k=0.85k=0.85 tier undershoots how much the n=2n=2 shell actually screens an outer nsns electron), but the qualitative conclusion — ZeffZ_{\text{eff}} is small (around 2, not 11 or 1) — is exactly the physically correct statement, and the order-of-magnitude drop in ionization energy compared to hydrogen's 1312 kJ/mol is captured. This is the essence of the article's opening claim: sodium ionizes easily not because its nucleus is weak, but because ten inner electrons collectively (and imperfectly) cancel nine-elevenths of the nuclear charge as seen by the valence electron.

Example 2 — Applied: comparing 3s, 3p ionization within chlorine, and rationalizing electron affinity trends. Chlorine, Z=17Z=17, configuration 1s2 2s22p6 3s23p51s^2\,2s^2 2p^6\,3s^2 3p^5. Compute ZeffZ_{\text{eff}} for a 3p3p electron. Group structure: (1s2)(2s22p6)(3s23p5)(1s^2)(2s^2 2p^6)(3s^2 3p^5). The electron of interest is one of the five 3p3p electrons. Screeners:

  • Same group (3s,3p3s,3p): total electrons in group =2+5=7=2+5=7; screening electron itself excluded, so 66 others, each contributing 0.350.35:
6×0.35=2.106 \times 0.35 = 2.10
  • n=2n=2 shell (8 electrons), each contributing 0.850.85:
8×0.85=6.808\times 0.85 = 6.80
  • n=1n=1 shell (2 electrons), each contributing 1.001.00:
2×1.00=2.002\times 1.00 = 2.00
σ=2.10+6.80+2.00=10.90\sigma = 2.10+6.80+2.00 = 10.90
Zeff=17−10.90=6.10Z_{\text{eff}} = 17-10.90 = 6.10

Now do the same for a 3s3s electron in chlorine (same configuration, screening a 3s3s instead of 3p3p — Slater's rules treat 3s,3p3s,3p identically as screeners and as the screened electron, so the same-group count is identical: 66 other electrons in the combined 3s,3p3s,3p group):

σ3s=6×0.35+8×0.85+2×1.00=2.10+6.80+2.00=10.90\sigma_{3s} = 6\times0.35 + 8\times0.85+2\times1.00 = 2.10+6.80+2.00=10.90

Slater's rules predict Zeff(3s)=Zeff(3p)=6.10Z_{\text{eff}}(3s) = Z_{\text{eff}}(3p) = 6.10 for chlorine — this is the well-known blind spot of Slater's original scheme: by lumping nsns and npnp into one group, it cannot resolve the ss–pp splitting that penetration actually produces. Physically, the 3s3s electron does penetrate more than 3p3p (as shown in the first derivation) and so should have a somewhat larger true ZeffZ_{\text{eff}} — more refined parameter sets (Clementi–Raimondi, from self-consistent-field calculations) give Zeff(3s)≈6.12Z_{\text{eff}}(3s)\approx 6.12 vs. Zeff(3p)≈6.76Z_{\text{eff}}(3p)\approx 6.76 for chlorine, restoring the physically correct order Zeff(3p)<Zeff(3s)Z_{\text{eff}}(3p) < Z_{\text{eff}}(3s). Chemically, this matters for chlorine's very negative electron affinity (−349-349 kJ/mol, the most exothermic of any element): the incoming electron enters a 3p3p orbital that is reasonably well shielded from other 3p3p electrons individually (each same-subshell partner only screens at k=0.35k=0.35) yet the whole atom still exerts a respectably large Zeff≈6.1Z_{\text{eff}} \approx 6.1 on it, giving strong net attraction with only modest electron–electron repulsion penalty — exactly the balance that makes halogens such effective electron acceptors compared with, say, oxygen or noble gases (full shell, extra electron forced into next nn, tiny ZeffZ_{\text{eff}}).

Example 3 — Exam/competition level: why 4s empties before 3d fills — Ca vs. Ca²⁺ vs. Sc, using explicit numbers. This is the classic "which orbital is lower, 4s4s or 3d3d?" puzzle, and penetration/shielding resolves the apparent contradiction between neutral-atom filling order and ion configurations.

Neutral calcium, Z=20Z=20: 1s2 2s22p6 3s23p6 4s21s^2\,2s^2 2p^6\,3s^2 3p^6\,4s^2 (not 3d23d^2). Compute ZeffZ_{\text{eff}} for the two competing orbitals as if populated, in the same underlying core 1s22s22p63s23p61s^2 2s^2 2p^6 3s^2 3p^6 (18 electrons): For a hypothetical 3d3d electron added to this core (Slater group rules: n∗=3n^*=3 for 3d3d, and using the d/fd/f tier where every electron in a lower group screens fully):

σ3d=18×1.00=18.00(all 18 core electrons are "lower" than 3d)\sigma_{3d} = 18\times 1.00 = 18.00 \quad(\text{all 18 core electrons are "lower" than } 3d)
Zeff(3d)=20−18.00=2.00,n∗=3Z_{\text{eff}}(3d) = 20-18.00 = 2.00,\qquad n^*=3
E3d∝(2.003)2=0.444E_{3d}\propto \left(\frac{2.00}{3}\right)^2 = 0.444

For a 4s4s electron added to this same core:

  • Same group (4s,4p4s,4p): only this one electron, contributes 00.
  • n=3n=3 shell (3s23p63s^23p^6, 8 electrons), one shell in, k=0.85k=0.85:
8×0.85=6.808\times0.85=6.80
  • n≤2n\le2 shells (1s22s22p61s^22s^22p^6, 10 electrons), two-or-more shells in, k=1.00k=1.00:
10×1.00=10.0010\times1.00=10.00
σ4s=6.80+10.00=16.80\sigma_{4s} = 6.80+10.00=16.80
Zeff(4s)=20−16.80=3.20Z_{\text{eff}}(4s) = 20-16.80=3.20

Slater's rule uses n∗=3.7n^{*}=3.7 for n=4n=4 (the empirical Slater correction for n≥4n\ge4: n∗=3.7n^*=3.7 at n=4n=4, 4.04.0 at n=5n=5, 4.24.2 at n=6n=6), so:

E4s∝(3.203.7)2=(0.865)2=0.748E_{4s}\propto\left(\frac{3.20}{3.7}\right)^2 = (0.865)^2=0.748

Comparing the two proportional energies (more negative = lower energy = more stable, so we want the larger magnitude of (Zeff/n∗)2(Z_{\text{eff}}/n^*)^2 since E=−(Zeff/n∗)2×E=-(Z_{\text{eff}}/n^*)^2\timesRy):

0.748 (for 4s)  >  0.444 (for 3d)0.748\ (\text{for }4s) \;>\; 0.444\ (\text{for }3d)

The 4s4s orbital is predicted (correctly) to be lower in energy than 3d3d for building up neutral calcium — even though 4s4s has a smaller ZeffZ_{\text{eff}} (3.20 vs. 2.00), its much smaller effective quantum number penalty (n∗=3.7n^*=3.7, close to 3d3d's n∗=3n^*=3) is not enough to overcome 4s4s's superior penetration advantage relative to how poorly-shielded it manages to stay. The n2n^2 denominator working against 4s4s is outweighed by 3d3d's far smaller ZeffZ_{\text{eff}} (since every one of the 18 core electrons fully shields 3d3d, versus only the outer eight doing so at 85% for 4s4s) — this is penetration in action: 4s4s has small inner lobes that let it "peek" past part of the core, while 3d3d has no such lobes and is shielded almost completely.

Now ionize: for Ca2+\mathrm{Ca}^{2+} (remove the two 4s4s electrons, giving argon-core configuration 1s22s22p63s23p61s^22s^22p^63s^23p^6), and ask which orbital an added electron would prefer — this models Sc2+→Sc+\mathrm{Sc}^{2+}\to\mathrm{Sc}^+ or explains why, once 3d3d starts filling in the transition series, 3d3d becomes lower than 4s4s (transition-metal cations have configurations like [Ar]3dn[\mathrm{Ar}]3d^n, never [Ar]3dn−24s2[\mathrm{Ar}]3d^{n-2}4s^2-for-the-ion). With the 4s4s electrons already gone, the 3d3d orbital's screening environment hasn't changed (still 18 electrons all fully screening it in this Slater approximation, since 3d3d's formula doesn't distinguish "core" further), but critically, once even one electron occupies 3d3d in the neutral transition metals, additional 3d3d electrons screen each other only at k=0.35k=0.35 (poor, same-subshell screening) while continuing to be almost perfectly shielded from outside by the core — so as ZZ increases across the first transition series, Zeff(3d)Z_{\text{eff}}(3d) rises faster than Zeff(4s)Z_{\text{eff}}(4s), and by scandium the 3d3d orbital has already dropped below 4s4s in orbital energy for the ion, even though the neutral atom is still built by filling 4s4s first due to interelectronic repulsion subtleties beyond simple Slater accounting. This is precisely the resolved paradox behind the famous "4s4s fills first, empties first" rule, and it is treated in full in Orbital Energy Ordering in Multielectron Atoms.

Limiting cases, scope & common traps

RegimeWhat happensWhy
Hydrogen / one-electron ionsNo ll-splitting at all; σ=0\sigma=0 alwaysNo other electrons exist to screen; Zeff=ZZ_{\text{eff}}=Z for every ll
Same nn, increasing ll (e.g. 3s,3p,3d3s,3p,3d)Energy rises (Ens<Enp<EndE_{ns}<E_{np}<E_{nd})Penetration falls with ll (fewer inner radial lobes), so σ\sigma rises and ZeffZ_{\text{eff}} falls
Across a period, same orbital typeZeffZ_{\text{eff}} rises smoothly (roughly +0.65+0.65 per element for outer s/ps/p in Slater's scheme, since only k=0.35k=0.35 same-group screening partly offsets each new proton)Each added proton is only ~35–65% cancelled by the new electron in the same shell
Down a group, same orbital typeZeffZ_{\text{eff}} rises only slightly; ionization energy usually fallsNew shells add both protons and near-complete inner screening from the previous shell; the dominant change is a larger nn (bigger orbital, more distance)
dd and ff orbitals specificallyAlmost no penetration; σ\sigma close to "one electron per inner charge"Zero radial nodes at minimum nn for that ll; explains why 4f4f orbitals in lanthanides are so contracted and non-bonding
Very high ZZ (heavy elements)Simple non-relativistic Slater/hydrogenic treatment breaks downRelativistic contraction of s,p1/2s,p_{1/2} orbitals becomes comparable in size to penetration effects (e.g. gold's color, mercury's liquidity)

The chemist behind it

Origins · 1 minJohn Clarke Slater and how this was discovered
JC
John Clarke Slater
1900–1976 · American theoretical physicist and chemist, MIT

Slater trained as a physicist under the early quantum mechanics revolution, working alongside Bohr and contributing to the Slater determinant formalism that enforces the Pauli exclusion principle in multielectron wavefunctions. His 1930 rules for estimating shielding constants were meant as a practical, hand-computable substitute for the full variational self-consistent-field calculations that were, at the time, prohibitively laborious without computers. Slater's broader career spanned band theory of solids, the theory of insulators and semiconductors, and some of the earliest systematic quantum-chemical calculations on molecules — his shielding rules survive in every introductory chemistry course as the simplest quantitative bridge between orbital shape and orbital energy.

History
- 1913: Bohr model gives orbit energies depending only on nn, with no mechanism for ll-splitting.
- 1926: Schrödinger equation solved exactly for hydrogen, confirming exact ll-degeneracy for one-electron systems.
- 1928–1930: Hartree develops the self-consistent field method, numerically computing true multielectron orbital energies and revealing the s<p<d<fs<p<d<f ordering within a shell.
- 1930: J. C. Slater publishes his empirical rules for estimating shielding constants and effective nuclear charge by hand.
- 1963: Clementi and Raimondi publish refined, orbital-resolved ZeffZ_{\text{eff}} values from accurate Hartree–Fock calculations, correcting known Slater-rule failures such as the 3s3s/3p3p degeneracy problem.
What exactly is the difference between shielding and screening — are they different words for the same thing?

In standard usage they are synonyms; "shielding constant" and "screening constant" both refer to σ\sigma, the amount of nuclear charge cancelled by other electrons as felt by a given electron. Some texts use "screening" for the general phenomenon and "shielding constant" specifically for the numerical parameter in Slater-type formulas, but there is no accepted distinction beyond stylistic preference. What must not be confused with either term is "penetration," which describes the orbital's own radial density near the nucleus rather than the effect other electrons have on it — penetration is a cause, shielding/screening is (partly) its consequence.

Does penetration mean an electron in a 3s orbital literally passes "through" the 1s and 2s electrons?

Loosely, yes, in the sense that the 3s3s radial distribution function has genuine, calculable probability density at radii smaller than the typical radius of the 1s1s and 2s2s shells — the electron is not confined to a shell "outside" the inner electrons the way a planetary orbit would be. There is no literal collision or trajectory crossing since these are probability densities, not classical paths, but the physical content is real: quantum mechanically, a 3s3s electron spends a measurable fraction of its time closer to the nucleus than the average radius of the inner shells, and during that time it experiences less screening and more of the bare nuclear attraction.

Why do d orbitals penetrate so little compared to s and p orbitals of the same shell?

Because the minimum principal quantum number for a dd orbital is n=3n=3, and at that minimum nn, dd orbitals have zero radial nodes — their radial wavefunction is a single lobe that behaves as r2r^2 near the origin (since l=2l=2), the steepest near-origin falloff of any orbital in that shell. Contrast with 3s3s (two radial nodes, finite density at r=0r=0) and 3p3p (one radial node, linear falloff). The absence of radial nodes for maximum-ll orbitals in a given shell is a general feature of hydrogenic solutions, and it is the direct mathematical reason dd (and even more so ff) orbitals are the least penetrating, most strongly shielded orbital type in any shell.

How does penetration explain why noble gas configurations are so stable, beyond just "full shell"?

A full ns2np6ns^2np^6 configuration means the next electron must go into a new shell with a much larger nn, and critically, that new electron's orbital — however penetrating in principle — must now get past a complete, tightly-packed inner shell rather than a partial one. The jump in shielding when moving from filling an existing shell to starting a new one is unusually large, because the new outer electron's small residual penetration into the just-completed shell contributes only a small correction to what is otherwise near-total screening by that full shell. This produces the sharp drop in ionization energy from (say) neon to sodium, and it is why the noble-gas boundary marks the clearest discontinuity in atomic properties across the periodic table.

Is the s < p < d < f energy ordering within a shell absolute, or can it be violated?

Within a single shell of a single, fixed atom or ion in its ground-state electron configuration, the ordering Ens<Enp<End<EnfE_{ns}<E_{np}<E_{nd}<E_{nf} from penetration/shielding is essentially universal and not violated — it is a robust consequence of the radial node-counting argument. What does get more subtle, and can appear to "violate" simple expectations, is comparing orbitals from different shells, such as 4s4s versus 3d3d: there, penetration (favoring 4s4s) competes with the smaller principal quantum number of 3d3d (favoring 3d3d through the 1/n21/n^2 scaling), and the outcome can flip depending on how many electrons already occupy each orbital and the overall charge state, as shown numerically for calcium versus its ions in the worked examples.

Why does Slater's original scheme fail to distinguish 3s from 3p, and does that matter chemically?

Slater's rules lump nsns and npnp into a single group for screening purposes, assigning them identical shielding contributions both as screeners and as the screened orbital — a simplification made for hand-calculation convenience in 1930, before computers made more refined self-consistent-field calculations routine. It does matter for precision: real ZeffZ_{\text{eff}} values (Clementi–Raimondi) show 3s3s experiences measurably less shielding than 3p3p, consistent with 3s3s's greater penetration. For most qualitative chemical arguments — predicting which element has a higher ionization energy, rationalizing electron affinity trends — the ss/pp distinction Slater misses is a second-order correction, but for precise spectroscopic term energies or subtle periodic anomalies, the refined values are necessary.

Before this

Where this leads