The Bohr Radius and Characteristic Atomic Scales

Every solid object you have ever touched — this screen, your hand, a mountain range — is built from atoms roughly one ångström across. Not a millimetre, not a picometre: always about $10^{-10}$ metres. Chemistry tabulates thousands of bond lengths and they all crowd between 1 and 2 Å. Nature owns a ruler, and it stamps the same tick spacing on every material in the universe. Where is that ruler kept? The answer is a single combination of constants that falls out of the hydrogen atom: the Bohr radius $a_0 = 4\pi\varepsilon_0\hbar^2/(m_e e^2) \approx 0.529\,\text{Å}$ — the radius at which an electron's confinement-driven kinetic energy, the price the uncertainty principle charges for squeezing a quantum cloud, exactly balances the Coulomb attraction pulling it inward. Squeeze tighter and the kinetic bill explodes; relax and the binding goes slack. Because the ingredients are universal constants, every atom everywhere strikes roughly the same bargain. This article is the chapter's currency exchange. From $a_0$ flows a complete system of natural units: the hartree $E_h = 27.2\,\text{eV}$ , the attosecond clock tick, the $5 \times 10^{11}\,\text{V/m}$ field an electron feels one Bohr radius from a proton — and the exchange rate between them all, the fine-structure constant $\alpha \approx 1/137$ , which doubles as the atom's speedometer. Master these scales and you can estimate almost any atomic quantity before solving anything.

Key takeaways

The Bohr radius $a_0 = 4\pi\varepsilon_0\hbar^2/(m_e e^2) = 0.529177\,\text{Å}$ is the unique length constructible from $\hbar$ , $m_e$ , and the Coulomb coupling — the natural size of any atom.
It marks the balance point where confinement kinetic energy ( $\sim 1/r^2$ ) matches Coulomb attraction ( $\sim 1/r$ ); minimizing the two-term budget gives $a_0$ and $-13.6\,\text{eV}$ exactly.
The companion energy unit is the hartree $E_h = e^2/(4\pi\varepsilon_0 a_0) = \alpha^2 m_e c^2 = 27.2114\,\text{eV}$ — twice the Rydberg; the atomic time unit is $\hbar/E_h = 24.19$ attoseconds.
The fine-structure constant $\alpha = 1/137.036$ is the atom's speedometer: the ground-state electron moves at $v \approx \alpha c$ , so relativity enters only at the $\alpha^2 \sim 5\times10^{-5}$ level.
The scales flex predictably: charge $Z$ shrinks lengths by $1/Z$ , a heavier orbiting mass shrinks them by $m_e/\mu$ , excitation inflates them by $n^2$ .
Key idea: quantum kinetic pressure against Coulomb attraction sets one length, $a_0$ — and from it, every atomic scale follows by pure bookkeeping.

The core insight

The Bohr radius is the length at which confinement-driven kinetic energy balances Coulomb attraction; minimizing $E(r) \sim \hbar^2/(2m_e r^2) - e^2/(4\pi\varepsilon_0 r)$ yields $r = a_0$ and fixes the size and energy scale of every atom.

Where this lives today

The Bohr radius is quietly present in every quantitative statement about matter. Density-functional-theory codes — the workhorses that design battery cathodes and catalysts — operate natively in bohrs and hartrees, exactly as Hartree set them up; every simulated molecule's geometry file is a list of coordinates in $a_0$ . In semiconductors the effective Bohr radius (band mass, screened $\varepsilon$ ) balloons to nanometres — it sets exciton sizes and with them the colours of QLED quantum dots: same formula, different currency. At the precision frontier the scales themselves are instruments: attosecond streaking clocks electron release against the $24.19\,\text{as}$ tick; strong-field experiments grade laser power against the atomic intensity. And muonic hydrogen — the atom shrunk 186-fold by a heavy orbiter — turned the flexing of $a_0$ into a probe of the proton: its spectroscopy delivered a charge radius precise enough to spark (and settle) the decade-long "proton radius puzzle."

The physical picture

Figure

The energy budget — quantum pressure against Coulomb attraction

Why doesn't the electron simply fall onto the proton? Classically it should — an orbiting charge radiates, spirals in, and crashes in about sixteen picoseconds. Quantum mechanics forbids the crash with an accounting rule: confining the electron's cloud to size $r$ forces a momentum spread $\Delta p \sim \hbar/r$ , hence a kinetic energy $\sim \hbar^2/(2m_e r^2)$ . That cost climbs as $1/r^2$ while the Coulomb reward only deepens as $1/r$ : near the origin the kinetic wall always wins. The bottom of the combined curve sits at one radius, and evaluating the constants there gives $0.529\,\text{Å}$ — the atom's size is an output of quantum mechanics, not an input. The same balance explains why all atoms are ångström-sized: a uranium atom binds its outermost electrons with an effective charge not much larger than hydrogen's (inner electrons screen the nucleus — machinery in Hydrogenic Ions and Scaling with Nuclear Charge), so the outermost shell always lands within a factor of a few of $a_0$ . Nature's ruler is universal because its tick spacing is built from universal constants.

Figure

The ×137 ladder — three lengths, one constant apart

There is a deeper way to read $a_0$ , through the one dimensionless number in the problem. The fine-structure constant $\alpha = e^2/(4\pi\varepsilon_0\hbar c) \approx 1/137$ organizes three famous lengths into a geometric ladder: the classical electron radius $r_e = 2.82\,\text{fm}$ , the reduced Compton wavelength $\lambda\!\!\!^{-}_C = \hbar/m_e c = 386\,\text{fm}$ , and the Bohr radius — each a factor $1/\alpha \approx 137$ above the last. The atom lives two factors of 137 above the scale of relativistic pair physics. That margin is why the non-relativistic Schrödinger equation describes hydrogen to parts in $10^4$ — and why the corrections, when they come, arrive ordered in powers of $\alpha$ .

The one length the constants allow

Before any minimization, a bookkeeping question: what lengths can be built at all from the theory's ingredients? The non-relativistic hydrogen Hamiltonian contains exactly three constants — $\hbar$ , the electron mass $m_e$ , and the Coulomb coupling $k = e^2/4\pi\varepsilon_0$ . The unique length they form is

a_0 = \frac{4\pi\varepsilon_0\hbar^2}{m_e e^2} = \frac{\hbar^2}{m_e k} = 5.29177 \times 10^{-11}\,\text{m} = 0.529177\,\text{Å} \qquad (1)

where: $\varepsilon_0$ — vacuum permittivity; $\hbar$ — reduced Planck constant; $m_e$ — electron mass (strictly the reduced mass $\mu$ — machinery in The Bohr Energy Spectrum and Rydberg Constant); $e$ — elementary charge. Read the structure: $\hbar$ appears squared in the numerator — erase quantum mechanics ( $\hbar \to 0$ ) and the only available length shrinks to zero. Note what is absent: no $c$ (non-relativistic physics), no proton mass (the nucleus is just an anchor). Every length is fixed once $a_0$ is; the same three constants also single out one energy. Two equivalent expressions for it are worth memorizing:

E_h = \frac{e^2}{4\pi\varepsilon_0 a_0} = \alpha^2 m_e c^2 = 27.2114\,\text{eV} \qquad (2)

where: $E_h$ — the hartree, the atomic unit of energy. The first form says: $E_h$ is the Coulomb energy of two elementary charges one Bohr radius apart. The second: the electron's rest energy taxed twice by the coupling — $511\,\text{keV} \times (1/137)^2$ . Beware the factor of two: the hydrogen binding energy is $E_h/2 = 13.6057\,\text{eV}$ , the Rydberg — half because the virial theorem leaves the bound electron with kinetic $+E_h/2$ against potential $-E_h$ . The dimensionless key linking all scales deserves its own display:

\alpha = \frac{e^2}{4\pi\varepsilon_0\hbar c} = \frac{1}{137.036}, \qquad v_1 = \alpha c = 2.188 \times 10^{6}\,\text{m/s} \qquad (3)

where: $v_1$ — the characteristic ground-state electron speed, $\hbar/(m_e a_0)$ . This is the atom's speedometer reading: about $0.7\%$ of light speed. Kinetic energy compares to rest energy as $(v_1/c)^2 = \alpha^2 \approx 5 \times 10^{-5}$ — the promised smallness that justifies non-relativistic quantum mechanics and simultaneously predicts the size of the fine-structure corrections that eventually force relativity back in. How those corrections escalate for heavier one-electron ions ( $v \sim Z\alpha c$ — at uranium the speedometer reads $0.67c$ ) is the business of Hydrogenic Ions and Scaling with Nuclear Charge.

Derivation AThe stability argument — one minimization, the whole atom

The goal is the most physical derivation of $a_0$ in existence: show that an electron confined near a proton has a minimum-energy configuration at exactly $r = a_0$ with energy exactly $-13.6\,\text{eV}$ , using nothing but the uncertainty principle and one derivative. The strategy is to write the energy as a function of a single number — the cloud size $r$ — by pricing confinement with $\Delta p \sim \hbar/r$ , then minimize. Watch the powers of $r$ : the whole result is the fight between $1/r^2$ and $1/r$ . And keep one honesty flag flying: this is an estimate whose landing on the exact answer is a special grace of the Coulomb problem, not a general property of such arguments.

Step 1. Price the confinement. A cloud of size $r$ has momentum spread $p \sim \hbar/r$ (the uncertainty floor — the rigorous inequality machinery is the Generalized Uncertainty Principle's), so its kinetic energy is $p^2/2m_e = \hbar^2/(2m_e r^2)$ . Add the Coulomb energy at distance $r$ :

E(r) = \frac{\hbar^2}{2m_e r^2} - \frac{e^2}{4\pi\varepsilon_0 r}.

Why this is legitimate: we are not solving anything — we are pricing the two competing effects at a common scale $r$ and trusting that the true ground state lives near the cheapest one.

Step 2. Minimize. Setting $dE/dr = 0$ :

-\frac{\hbar^2}{m_e r^3} + \frac{e^2}{4\pi\varepsilon_0 r^2} = 0 \quad\Longrightarrow\quad r_{\min} = \frac{4\pi\varepsilon_0\hbar^2}{m_e e^2} = a_0.

Why a minimum and not a maximum: as $r \to 0$ the positive $1/r^2$ term dominates and $E \to +\infty$ ; as $r \to \infty$ both terms die and $E \to 0^-$ . The stationary point between is the unique minimum ( $d^2E/dr^2 > 0$ there, verified below).

Step 3. Evaluate the energy at the minimum:

E(a_0) = \frac{\hbar^2}{2m_e a_0^2} - \frac{e^2}{4\pi\varepsilon_0 a_0} = \frac{E_h}{2} - E_h = -\frac{E_h}{2} = -13.6057\,\text{eV}.

Why the terms land in a 1:2 ratio: at the balance radius the kinetic term is exactly half the magnitude of the potential term — the virial theorem's $2\langle T\rangle = -\langle V\rangle$ appearing already in the estimate.

Step 4. State the honesty clause. The estimate reproduces the exact Schrödinger ground state — radius scale and energy — with no fudge factor. That precision is luck: writing $\Delta p = \hbar/r$ rather than $\hbar/2r$ was a choice, and for other potentials the same argument lands only within factors of order one. What is not luck is the scaling: any argument balancing $\hbar^2/m r^2$ against $k/r$ must produce $r \sim \hbar^2/mk$ and $E \sim -mk^2/\hbar^2$ . The exact solve that certifies the numbers is The Hydrogen Atom: Coulomb Radial Equation. Verified: minimization done symbolically (sympy) — $r_{\min} = \hbar^2/m_e k$ exactly, $E_{\min} = -m_e k^2/2\hbar^2$ exactly, $d^2E/dr^2 = k^4 m_e^3/\hbar^6 > 0$ ; CODATA values give $a_0 = 5.29177\times10^{-11}$ m and $E_{\min} = -13.6057$ eV.

Go deeper · Derivation B · 1 minWhy a₀ and nothing else — the dimensional bookkeeping

Derivation BWhy a₀ and nothing else — the dimensional bookkeeping

The goal is to prove the claim "the unique length constructible from $\hbar$ , $m_e$ , $k$ " rather than assert it — and, in passing, to see that the whole system of atomic units is forced, not chosen. The strategy is exponent bookkeeping: demand $\hbar^{x} m_e^{y} k^{z}$ have the dimensions of a length and show the exponents are uniquely determined. This is a three-equations-in-three-unknowns problem, and its non-degeneracy is the entire content.

Step 1. Tabulate dimensions in mass–length–time: $[\hbar] = ML^2T^{-1}$ , $[m_e] = M$ , $[k] = [\text{energy} \times \text{length}] = ML^3T^{-2}$ .

Step 2. Demand $[\hbar^x m_e^y k^z] = L$ . Matching exponents of $M$ , $L$ , $T$ in turn:

x + y + z = 0, \qquad 2x + 3z = 1, \qquad -x - 2z = 0.

Why three equations: one per base dimension — the system is square because we have exactly three constants.

Step 3. Solve: $x = 2$ , $y = -1$ , $z = -1$ , uniquely (the coefficient determinant is nonzero). Hence

a_0 = \frac{\hbar^2}{m_e k}

is the only length in non-relativistic hydrogen — no alternative combination exists, and no dimensionless freedom remains except an overall pure number, which the Schrödinger solution fixes to exactly 1. Why uniqueness matters physically: it means the theory has no size knob. Once $\hbar$ , $m_e$ , $e$ are measured anywhere in the universe, atomic sizes everywhere are fixed.

Step 4. Repeat for energy ( $ML^2T^{-2}$ ): uniquely $x = -2$ , $y = 1$ , $z = 2$ , giving $E_h = m_e k^2/\hbar^2$ — and for time, $t_0 = \hbar/E_h$ . The whole currency system is thereby derived: setting $\hbar = m_e = e = 4\pi\varepsilon_0 = 1$ is legitimate because those assignments exhaust the freedom. Verified: both linear systems solved numerically (numpy) — exponents $(2,-1,-1)$ and $(-2,1,2)$ , determinant nonzero (uniqueness); dimensional identities cross-checked against CODATA values to $10^{-9}$ relative.

The atomic currency system

Hartree's proposal — work in the atom's own units and the equations shed their constants — turns equation (1) into a full exchange table. In atomic units ( $\hbar = m_e = e = 4\pi\varepsilon_0 = 1$ ) the hydrogen Hamiltonian is just $H = -\tfrac12\nabla^2 - 1/r$ ; results convert back through this ladder:

quantity	atomic unit	SI / practical value	remember it as
length	$a_0$	$0.529177\,\text{Å}$	the atom's radius tick
energy	$E_h = m_e e^4/(4\pi\varepsilon_0)^2\hbar^2$	$27.2114\,\text{eV}$	twice the ionization energy of H
time	$\hbar/E_h$	$24.19\,\text{as}$	the electron's clock tick
velocity	$v_1 = \alpha c$	$2.188\times10^6\,\text{m/s}$	$c/137$ — the speedometer
electric field	$E_h/(e a_0)$	$5.142\times10^{11}\,\text{V/m}$	the field the electron actually feels
intensity	$\tfrac12 c\varepsilon_0 (E_h/ea_0)^2$	$3.51\times10^{16}\,\text{W/cm}^2$	where lasers rival the nucleus

The time entry rewards a pause. An electron at $v_1$ crosses one Bohr radius in $a_0/v_1 = 24.19$ attoseconds; a full "orbit" $2\pi a_0/v_1$ takes $152\,\text{as}$ . Chemistry's fastest events — charge migration after ionization, the first femtosecond of a breaking bond — play out in a few of these ticks; the 2023 Nobel Prize went to the laser techniques that built a camera with attosecond shutter speed, and its reference frame is this table's number.

\rho = \frac{r}{a_0}: \qquad \frac{d^2u}{d\rho^2} = \left[\frac{\ell(\ell+1)}{\rho^2} - \frac{2}{\rho} + \frac{|E|}{\text{Ry}}\right]u \qquad (4)

where: $u = rR$ — the reduced radial function (the reduction machinery is Effective Potential and the Centrifugal Barrier's); $\text{Ry} = E_h/2$ . This is the shape the radial equation takes when lengths are measured in $a_0$ : every physical constant has vanished. The next derivation shows that demanding this disappearance is itself a route to $a_0$ — the rigorous counterpart of the stability estimate.

Go deeper · Derivation C · 1 mina₀ as the scale that nondimensionalizes the radial equation

Derivation Ca₀ as the scale that nondimensionalizes the radial equation

The goal is to establish $a_0$ the honest way: not from a heuristic energy budget but from the exact dynamics. We show that the hydrogen radial equation forces the combination $4\pi\varepsilon_0\hbar^2/m_e e^2$ on us as the only length that strips it of all dimensional clutter — so $a_0$ is intrinsic to the Schrödinger problem, present before any solution is written down. The strategy: substitute $r = a\rho$ with $a$ unknown, and demand the kinetic and Coulomb coefficients merge into pure numbers.

Step 1. Start from the reduced radial equation (quoted from The Radial Equation and Effective Potential and the Centrifugal Barrier, which own its derivation):

-\frac{\hbar^2}{2m_e}\frac{d^2u}{dr^2} + \left[\frac{\hbar^2\ell(\ell+1)}{2m_e r^2} - \frac{e^2}{4\pi\varepsilon_0 r}\right]u = E\,u.

Step 2. Substitute $r = a\rho$ , so $d^2/dr^2 = a^{-2}\,d^2/d\rho^2$ , and divide through by the kinetic prefactor $\hbar^2/(2m_e a^2)$ :

-\frac{d^2u}{d\rho^2} + \left[\frac{\ell(\ell+1)}{\rho^2} - \frac{2m_e a}{\hbar^2}\cdot\frac{e^2}{4\pi\varepsilon_0}\cdot\frac{1}{\rho}\right]u = \frac{2m_e a^2 E}{\hbar^2}\,u.

Why divide by the kinetic coefficient: the centrifugal term then automatically becomes the pure number $\ell(\ell+1)/\rho^2$ — only the Coulomb coefficient still carries dimensions, and one choice of $a$ remains to fix it.

Step 3. Demand the Coulomb coefficient equal the pure number $2$ (the conventional choice that makes the bound-state condition come out in integers):

\frac{2m_e a}{\hbar^2}\cdot\frac{e^2}{4\pi\varepsilon_0} = 2 \quad\Longrightarrow\quad a = \frac{4\pi\varepsilon_0\hbar^2}{m_e e^2} = a_0.

Why this is unique: any other pure-number choice rescales $a$ by that number — but the combination of constants is forced. The dimensional-bookkeeping card proved no second length exists; here the dynamics selects the conventional normalization.

Step 4. The energy coefficient then reads $2m_e a_0^2 E/\hbar^2 = E/\text{Ry}$ with $\text{Ry} = \hbar^2/(2m_e a_0^2) = 13.6057\,\text{eV}$ , and the equation lands in the constant-free form (4). Everything dimensional about hydrogen now lives in two conversion factors, $a_0$ and Ry; the differential equation itself is arithmetic. Solving it — series, termination, the $-1/n^2$ ladder — is the business of The Hydrogen Atom: Coulomb Radial Equation. Verified: substitution algebra checked symbolically; $\text{Ry} = \hbar^2/(2m_e a_0^2) = 13.6057$ eV and $\text{Ry} = E_h/2$ confirmed with CODATA constants to $10^{-6}$ relative.

Limiting cases

case	scales become	meaning
$\hbar \to 0$	$a_0 \to 0$ , $E_h \to \infty$ (as $1/\hbar^2$ )	classical physics has no atomic size — matter collapses
nuclear charge $Z$	$a \to a_0/Z$ , $E \to Z^2 E_h$	stronger pull, tighter and hotter atom — machinery at Hydrogenic Ions and Scaling with Nuclear Charge
finite nuclear mass $M$	$m_e \to \mu = m_e M/(m_e{+}M)$ : $a_0$ grows by $\approx 0.05\%$ for H	the nucleus recoils; the correction that split hydrogen from deuterium — The Bohr Energy Spectrum and Rydberg Constant
heavier orbiter (muon)	all lengths $\times\, m_e/\mu_{\text{red}} \approx 1/186$	muonic hydrogen: a 285 fm atom that grazes the proton
high excitation $n$	size $\sim n^2 a_0$ , spacing $\sim E_h/n^3$	Rydberg atoms: micron-sized, nearly classical
$\alpha$ larger (thought experiment)	$a_0 \propto 1/\alpha$ shrinks, $v_1/c = \alpha$ grows	chemistry becomes relativistic; at $\alpha \sim 1$ atoms are unrecognizable

Common trap“a₀ is the radius of the electron's orbit”

There is no orbit. The ground-state electron is a standing cloud $\psi_{100} \propto e^{-r/a_0}$ with nonzero amplitude at every radius — including inside the nucleus. What $a_0$ actually marks: the most probable radius (peak of the shell probability $P(r) = 4\pi r^2|\psi|^2$ ), while the mean sits farther out at $\langle r\rangle = \tfrac32 a_0$ — different "sizes" for one state, disentangled in Radial Probability Distribution and Expectation Values. The Bohr-model planetary electron circling at $a_0$ gets the length right (hence the constant's name) but the kinematics wrong: nothing circulates in an $s$ -state — its angular momentum is zero, where Bohr's model demanded $\hbar$ . Keep the number, drop the orbit.

Worked examples

1 · Building the currency table from scratch (confidence builder)

Using only CODATA constants — $\hbar = 1.0546\times10^{-34}\,\text{J}\,\cdot\,\text{s}$ , $m_e = 9.1094\times10^{-31}\,\text{kg}$ , $e = 1.6022\times10^{-19}\,\text{C}$ , $\varepsilon_0 = 8.8542\times10^{-12}\,\text{F/m}$ , $c = 2.9979\times10^8\,\text{m/s}$ — reconstruct every atomic unit.

Step 1. The length. $a_0 = 4\pi\varepsilon_0\hbar^2/(m_e e^2)$ : numerator $4\pi(8.8542\times10^{-12})(1.0546\times10^{-34})^2 = 1.2374\times10^{-78}$ , denominator $(9.1094\times10^{-31})(1.6022\times10^{-19})^2 = 2.3383\times10^{-68}$ ; ratio $\mathbf{a_0 = 5.2918\times10^{-11}\,m}$ .

Step 2. The energy, first route: $E_h = e^2/(4\pi\varepsilon_0 a_0) = (1.6022\times10^{-19})^2/[4\pi(8.8542\times10^{-12})(5.2918\times10^{-11})] = 4.3597\times10^{-18}\,\text{J} = \mathbf{27.211\,eV}$ .

Step 3. The energy, second route as a cross-check: $\alpha = e^2/(4\pi\varepsilon_0\hbar c) = 1/137.036$ , so $E_h = \alpha^2 m_e c^2 = (511.00\,\text{keV})/137.036^2 = \mathbf{27.211\,eV}$ — the two routes agree to five digits, confirming the identity $a_0 = \hbar/(\alpha m_e c)$ that welds the table together.

Step 4. Time and speed. $t_0 = \hbar/E_h = \mathbf{24.19\,as}$ ; $v_1 = \alpha c = \mathbf{2.188\times10^6\,m/s}$ . Sanity check: $a_0/v_1 = 24.19\,\text{as} = t_0$ — self-consistent by construction.

Step 5. The field: $\mathcal{E}_0 = E_h/(e a_0) = \mathbf{5.142\times10^{11}\,V/m}$ . Laboratory insulators break down near $10^8$ V/m — the inside of an atom sustains five thousand times that, permanently.

2 · When does a laser out-shout the nucleus? (exam level)

An intense femtosecond laser is focused onto hydrogen atoms. At what intensity does the laser's peak electric field equal the atomic field $\mathcal{E}_0 = 5.142\times10^{11}$ V/m — the regime where the laser competes with the nucleus for custody of the electron?

Step 1. Relate intensity to peak field. A linearly polarized wave of peak field $E_0$ carries cycle-averaged intensity $I = \tfrac12 c\varepsilon_0 E_0^2$ (the Poynting-flux average — the factor $\tfrac12$ from $\langle\cos^2\rangle$ ).

Step 2. Insert the atomic field: $I = \tfrac12 (2.9979\times10^8)(8.8542\times10^{-12})(5.142\times10^{11})^2 = 3.51\times10^{20}\,\text{W/m}^2 = \mathbf{3.51\times10^{16}\,W/cm^2}$ .

Step 3. Interpret. At this benchmark the "perturbation" equals the binding field: perturbation theory is dead, the Coulomb barrier bends down, and the electron escapes by field ionization within an optical cycle. Amplified femtosecond systems reach $10^{18}$ – $10^{22}$ W/cm² — routinely above the atomic unit — which is why strong-field physics (tunnel ionization, high-harmonic generation, attosecond pulses) is mainstream experiment, not thought experiment.

Step 4. Corollary: since $I \propto E_0^2$ , a laser four orders below the benchmark still exerts $1\%$ of the atomic field — enough to Stark-shift levels measurably. The atomic intensity is the ruler against which all laser–atom physics is graded.

Origins · 1 minDouglas Rayner Hartree and how this was discovered

Douglas Rayner Hartree

1897–1958 · Cambridge and Manchester

Hartree's problem was drowning in exponents. Computing an atom's electron distribution by hand — his life's project from 1926 on — meant marching the radial Schrödinger equation outward step by numerical step, every line dragging factors of $10^{-34}$ , $10^{-31}$ , $10^{-19}$ . His fix was disarmingly practical: let the atom set the units. Measure length in $a_0$ and energy in $e^2/(4\pi\varepsilon_0 a_0)$ — now called the hartree — and every constant becomes exactly 1; the surviving numbers are order-one, and a human with a desk calculator can trust their own arithmetic for the months a single atom required. The 1928 papers introducing his self-consistent-field method open with precisely this currency reform, and the human computer cranking many of the integrations was his father, William Hartree, a retired engineer who took up his son's arithmetic as a second career. The units outlived the hand computation: every modern electronic-structure code still reports its energies in hartrees. The self-consistent-field method itself — quantum mechanics' first industrial-scale approximation scheme — has its own story in the many-electron chapters.

History

1913

— Bohr's model of hydrogen introduces the radius a₀ = 0.53 Å as the innermost allowed orbit, deriving atomic size from constants for the first time.

1926

— Schrödinger's wave mechanics re-derives a₀ with the orbit removed: it re-emerges as the decay scale of the ground-state cloud e^(−r/a₀).

1928

— Hartree, computing atoms numerically, proposes atomic units: lengths in a₀, energies in what is now the hartree — the constants vanish from the working equations.

1959

— Shull and Hall propose the name "hartree"; the currency becomes standard vocabulary.

2018

— CODATA fixes a₀ = 0.529177210903 Å to twelve digits, its precision riding on measurements of α — the atom's size is now among the best-known lengths in physics.

Going further

Texts. Griffiths & Schroeter, Introduction to Quantum Mechanics, 3rd ed., §4.2 · Sakurai & Napolitano, Modern Quantum Mechanics, 2nd ed., §3.7 · Shankar, Principles of Quantum Mechanics, 2nd ed., Ch. 13 · Cohen-Tannoudji, Diu & Laloë, Vol. 1, Ch. VII · CODATA 2018 recommended values (NIST).

Provenance. All constants are CODATA-recommended standard values: $a_0=5.29177\times10^{-11}$ m, $\mathrm{Ry}=13.6057$ eV, $\alpha^{-1}=137.036$ , verifiable in Griffiths §4.2 and NIST tables. A verifier should recompute $a_0=4\pi\varepsilon_0\hbar^2/(m_e e^2)$ from SI constants, confirm $\mathrm{Ry}=\hbar^2/(2m_e a_0^2)=m_e e^4/[2(4\pi\varepsilon_0)^2\hbar^2]=13.6$ eV, verify the minimization $r_{\min}=a_0$ , and check the integrals $r_{\text{peak}}=a_0$ , $\langle r\rangle=\tfrac{3}{2}a_0$ against Griffiths/Shankar tables of $\langle r^k\rangle$ . Edition note (2026-07-08). The base's radial-distribution derivation ( $r_{\text{peak}} = a_0$ , $\langle r\rangle = \tfrac32 a_0$ ) is stated here with its owner named: the probability craft belongs to Radial Probability Distribution and Expectation Values. All numbers in this article were machine-reverified against CODATA 2018.

Prerequisites. The radial equation · the hydrogen atom: Coulomb radial equation · effective potential and the centrifugal barrier. Leads to. The Bohr energy spectrum and Rydberg constant · hydrogenic wavefunctions and quantum numbers · radial probability distribution and expectation values · hydrogenic ions and scaling with nuclear charge.

Is a₀ something we measure, or something we compute?

Both, and they agree — that is the point. Computed, $a_0$ combines $\hbar$ , $m_e$ , $e$ , $\varepsilon_0$ , each measured independently; measured, atomic sizes appear in scattering cross-sections, X-ray bond lengths, and Rydberg spectroscopy, all consistent with it. Its twelve-digit CODATA precision actually rides on the fine-structure constant via $a_0 = \hbar/(\alpha m_e c)$ — so today the "size of the atom" is known better than almost any directly ruled length.

The stability argument used Δp ~ ħ/r. Wouldn't ħ/2r (the actual uncertainty bound) change the answer?

It would move the estimate — that choice is exactly the argument's slack. With $p = \hbar/2r$ the minimum lands at $4a_0$ with a quarter of the energy: right scaling, wrong prefactor. The choice $p = \hbar/r$ happens to reproduce the exact answer — a fluke of the Coulomb problem, not a derivation. Trust the argument for scales; prefactors belong to the honest solve in The Hydrogen Atom: Coulomb Radial Equation.

Why is the hartree, not the rydberg, the "atomic unit" of energy?

Because the hartree is what the unit system produces: with $\hbar = m_e = e = 4\pi\varepsilon_0 = 1$ , the energy unit is $m_e e^4/(4\pi\varepsilon_0)^2\hbar^2 = E_h$ , full stop. The rydberg $= E_h/2$ owes its factor to hydrogen's virial split and survives for spectroscopy's historical reasons. The safe habit: check whether a paper's "a.u." means hartrees (almost always, in quantum chemistry) before trusting any factor of 2.

If the electron has nonzero probability inside the nucleus, why doesn't it react with the proton?

Presence is not reaction. The integrated probability of being inside the proton is of order $10^{-14}$ (computed honestly in Radial Probability Distribution and Expectation Values), and more importantly there is no energetically allowed channel at eV scales: electron capture ( $p + e^- \to n + \nu$ ) costs 0.78 MeV. In hydrogen the contact amplitude shows up instead as hyperfine structure — the 21 cm line — not as decay.

Does a₀ change in a molecule or a solid?

The constant doesn't, but its role is played by a rescaled understudy. In a semiconductor an electron of effective mass $m^*$ moves through dielectric $\varepsilon_r$ , so bound states (excitons, dopant levels) have radius $a^* = a_0\,\varepsilon_r\,(m_e/m^*)$ — commonly 1–10 nm. That inflation is why hydrogen-like physics reappears at the nanoscale and why quantum-dot colours are size-tunable. The bare $a_0$ still sets the underlying chemistry — bond lengths, core levels — inside any material.

Why does α keep appearing in a non-relativistic article — isn't it a relativity constant?

$\alpha$ is a pure coupling strength. Here it appears only as a ratio of scales — $v_1/c$ , $a_0/\lambda\!\!\!^{-}_C$ — locating hydrogen relative to relativistic physics without invoking relativistic dynamics. The dynamics enters later: fine-structure corrections of relative size $\alpha^2$ , and for nuclear charge $Z$ the escalation $v \sim Z\alpha c$ that eventually breaks the non-relativistic treatment — continued in Hydrogenic Ions and Scaling with Nuclear Charge.

Is there a gravitational Bohr radius?

Formally yes: replace the Coulomb coupling by $Gm_e m_p$ and the same formula gives $a_0^{\text{grav}} \approx 1.2\times10^{29}$ m — larger than the observable universe. That absurdity is a real statement: gravity is $\sim10^{39}$ times weaker than electromagnetism at particle level, so gravitational "atoms" cannot exist. It also flags what $a_0$ silently assumes — electromagnetism alone shapes atoms.