Skip to content

Hydrogen Atom

The hydrogen atom is the nonrelativistic electron–proton Coulomb problem after separation of center-of-mass and relative motion. Its internal Hamiltonian is exactly soluble and supplies the canonical quantum numbers, orbitals, radial functions, degeneracies, and energy scale of one-electron atoms.

The exactness belongs to a sharply defined model: point charges, instantaneous Coulomb interaction, no relativistic terms, no coupling to the quantized electromagnetic field, no nuclear structure, and no external fields. Precision hydrogen requires a hierarchy of corrections beyond this baseline.

FieldStandard internal hydrogen model
Degrees of freedomElectron–proton relative coordinate r\mathbf r
Hilbert spaceL2(R3,d3r)L^2(\mathbb R^3,d^3r), with spin omitted
Hamiltonianp^ 2/(2μ)−e2/(4πϵ0r)\hat{\mathbf p}^{\,2}/(2\mu)-e^2/(4\pi\epsilon_0r)
Reduced massμ=memp/(me+mp)\mu=m_em_p/(m_e+m_p)
Bound-state labelsn,ℓ,mn,\ell,m
Bound energiesEn∝−1/n2E_n\propto-1/n^2
Bound spatial degeneracyn2n^2 in the ideal Coulomb model
Continuum thresholdE=0E=0
SolvabilityExact separation into spherical harmonics and associated Laguerre functions
Key special symmetryConserved Runge–Lenz structure in addition to rotations
Canonical homeHydrogen Atom

For a hydrogenic ion with nuclear charge +Ze+Ze, replace the Coulomb coupling by Ze2Z e^2. The Hydrogen Atom Hamiltonian is the compact operator card; this page identifies the model and its physical scope.

Before separating coordinates, the spinless nonrelativistic Hamiltonian is

H^tot=p^e 22me+p^N 22mN−Ze24πϵ0∣re−rN∣.\hat H_{\mathrm{tot}} =\frac{\hat{\mathbf p}_e^{\,2}}{2m_e} +\frac{\hat{\mathbf p}_N^{\,2}}{2m_N} -\frac{Ze^2} {4\pi\epsilon_0 \lvert\mathbf r_e-\mathbf r_N\rvert}.

Define total and reduced masses

M=me+mN,μ=memNme+mN,M=m_e+m_N, \qquad \mu=\frac{m_em_N}{m_e+m_N},

and use center-of-mass and relative coordinates. The Hamiltonian separates:

H^tot=P^ 22M+H^rel,\hat H_{\mathrm{tot}} =\frac{\hat{\mathbf P}^{\,2}}{2M} +\hat H_{\mathrm{rel}}, H^rel=p^ 22μ−Ze24πϵ0r.\hat H_{\mathrm{rel}} =\frac{\hat{\mathbf p}^{\,2}}{2\mu} -\frac{Ze^2}{4\pi\epsilon_0r}.

The free center-of-mass plane wave does not affect the internal spectrum. The model card concerns H^rel\hat H_{\mathrm{rel}}. For ordinary hydrogen, Z=1Z=1 and mN=mpm_N=m_p.

Using mem_e instead of μ\mu is the infinite-nuclear-mass approximation. It is often adequate for qualitative work but should not be silently mixed with reduced-mass spectroscopic formulas.

The internal spinless Hilbert space is

Hrel=L2(R3,d3r).\mathcal H_{\mathrm{rel}} =L^2(\mathbb R^3,d^3r).

The Coulomb singularity is part of a well-defined self-adjoint Schrödinger operator on its standard domain. Physical bound states are square integrable, regular at the origin in the full three-dimensional sense, and decay at infinity.

After angular separation,

ψnℓm(r,θ,ϕ)=Rnℓ(r)Yℓm(θ,ϕ).\psi_{n\ell m}(r,\theta,\phi) =R_{n\ell}(r)Y_\ell^m(\theta,\phi).

The radial normalization is

∫0∞∣Rnℓ(r)∣2r2 dr=1.\int_0^\infty \lvert R_{n\ell}(r)\rvert^2r^2\,dr=1.

If unℓ(r)=rRnℓ(r)u_{n\ell}(r)=rR_{n\ell}(r), then

∫0∞∣unℓ(r)∣2 dr=1,\int_0^\infty \lvert u_{n\ell}(r)\rvert^2\,dr=1,

and the regular radial boundary condition is u(0)=0u(0)=0. Confusing RR, uu, and the radial probability density is one of the most common errors in this model.

SymbolMeaningRole
mem_eElectron massSets the infinite-nuclear-mass scale
mNm_NNuclear massEnters the reduced mass
μ\muElectron–nucleus reduced massControls orbital length and binding energy
ZZNuclear charge numberZ=1Z=1 for hydrogen
eeElementary charge magnitudeCoulomb coupling
ϵ0\epsilon_0Vacuum permittivitySI-unit convention
αfs\alpha_{\mathrm{fs}}Fine-structure constantNonrelativistic expansion parameter through ZαfsZ\alpha_{\mathrm{fs}}

Define the reduced-mass Bohr radius

aμ=4πϵ0ℏ2μe2a_\mu =\frac{4\pi\epsilon_0\hbar^2}{\mu e^2}

and the hydrogenic length scale

aZ=aμZ.a_Z=\frac{a_\mu}{Z}.

The corresponding binding-energy scale is

EC=ℏ22μaμ2=μc2αfs22.E_{\mathrm C} =\frac{\hbar^2}{2\mu a_\mu^2} =\frac{\mu c^2\alpha_{\mathrm{fs}}^2}{2}.

For a hydrogenic ion, the leading binding scale grows as μZ2\mu Z^2, while the orbital scale shrinks as 1/(μZ)1/(\mu Z).

In Hartree atomic units with an infinitely heavy nucleus, the internal Hamiltonian is

H^=−12∇2−Zr,\hat H =-\frac12\nabla^2-\frac{Z}{r},

and energies are

En=−Z22n2E_n=-\frac{Z^2}{2n^2}

in hartrees. Do not combine this compact form with SI parameters without translating units.

Rotational symmetry permits simultaneous eigenstates of

H^,L^2,L^z.\hat H, \qquad \hat L^2, \qquad \hat L_z.

The bound-state quantum numbers obey

n=1,2,3,…,n=1,2,3,\ldots, ℓ=0,1,…,n−1,\ell=0,1,\ldots,n-1, m=−ℓ,−ℓ+1,…,ℓ.m=-\ell,-\ell+1,\ldots,\ell.

The angular functions are Spherical Harmonics, and the radial polynomial factors use Laguerre Polynomials. The complete separation and normalization are owned by Radial Schrödinger Equation and Radial Wavefunctions.

Exact solvability includes:

  • the complete discrete bound spectrum;
  • normalized bound eigenfunctions;
  • analytic continuum Coulomb functions and scattering phases;
  • separation into radial and angular sectors;
  • algebraic symmetry methods for the bound spectrum.

Closed form does not mean every observable is a one-line expression. Continuum normalization, matrix elements, sums over intermediate states, and precision corrections can remain technically demanding.

The nonrelativistic hydrogenic bound energies are

En=−μZ2e42(4πϵ0)2ℏ21n2,n=1,2,3,….E_n =-\frac{\mu Z^2e^4} {2(4\pi\epsilon_0)^2\hbar^2} \frac{1}{n^2}, \qquad n=1,2,3,\ldots.

Equivalently,

En=−μc2(Zαfs)22n2.E_n =-\frac{\mu c^2(Z\alpha_{\mathrm{fs}})^2}{2n^2}.

The zero of energy is the separated electron–nucleus threshold. Bound states have E<0E<0, while the continuum has E>0E>0. The bound levels accumulate at zero from below.

For fixed nn, the number of spatial states is

gn=∑ℓ=0n−1(2ℓ+1)=n2.g_n =\sum_{\ell=0}^{n-1}(2\ell+1) =n^2.

Rotational symmetry explains degeneracy in mm for any central potential. The additional independence from ℓ\ell is special to the 1/r1/r potential and is tied to the conserved quantum Runge–Lenz structure. Degeneracy of the Hydrogen Atom develops this hidden symmetry.

Appending electron spin while omitting all spin-dependent terms doubles the spatial count to 2n22n^2. Real fine, hyperfine, radiative, nuclear-size, and external-field terms split parts of this ideal degeneracy.

The normalized ground state is

ψ100(r)=1πaZ3exp⁡(−raZ).\psi_{100}(\mathbf r) =\frac{1}{\sqrt{\pi a_Z^3}} \exp\left(-\frac{r}{a_Z}\right).

Its probability density is largest at the origin, but the probability per unit radius is

P100(r)=4r2aZ3exp⁡(−2raZ),P_{100}(r) =4\frac{r^2}{a_Z^3} \exp\left(-\frac{2r}{a_Z}\right),

which is maximal at r=aZr=a_Z.

For a general bound state:

  • the number of radial nodes is n−ℓ−1n-\ell-1;
  • angular nodes are determined by YℓmY_\ell^m;
  • parity is (−1)ℓ(-1)^\ell;
  • s,p,d,…s,p,d,\ldots label ℓ=0,1,2,…\ell=0,1,2,\ldots;
  • complex mm eigenstates and real orbital combinations are different bases within degenerate subspaces.

An orbital is a stationary one-electron wavefunction, not a path followed by an electron. Atomic Orbitals owns the basis, node, and visualization conventions.

The most common observables and diagnostics are:

QuantityWhat it probes
H^\hat HInternal binding energy
L^2\hat L^2 and L^z\hat L_zOrbital angular momentum labels
ParitySpatial inversion symmetry and selection structure
rr and radial projectorsRadial size and probability
Electric dipole ∝r^\propto\hat{\mathbf r}Leading optical transition amplitudes
L^\hat{\mathbf L} and spin magnetic momentZeeman response after spin is included
Continuum fluxIonization and Coulomb scattering

For the ground state,

⟨r⟩100=32aZ,⟨1r⟩100=1aZ.\langle r\rangle_{100} =\frac32a_Z, \qquad \left\langle\frac{1}{r}\right\rangle_{100} =\frac{1}{a_Z}.

The Coulomb virial theorem gives

2⟨T⟩=−⟨V⟩,2\langle T\rangle=-\langle V\rangle,

so for any stationary bound state,

⟨T⟩=−En,⟨V⟩=2En.\langle T\rangle=-E_n, \qquad \langle V\rangle=2E_n.

Electric-dipole matrix elements obey Δℓ=±1\Delta\ell=\pm1 and polarization-dependent Δm=0,±1\Delta m=0,\pm1 in the ideal basis. The full tensor-operator statement belongs at Dipole Transitions.

Hydrogen is the canonical model for:

  • exact two-body reduction and reduced mass;
  • central-potential separation in three dimensions;
  • spherical harmonics and radial special functions;
  • discrete and continuum sectors in one Hamiltonian;
  • principal, orbital, and magnetic quantum numbers;
  • radial probability measures and nodal structure;
  • ordinary rotational degeneracy versus hidden Coulomb degeneracy;
  • selection rules and spectroscopic line organization;
  • symmetry breaking by controlled corrections;
  • radial-solver and basis-set benchmarks.

It is also the bridge from one-particle quantum mechanics to atomic structure. The labels survive in approximate form for more complex atoms even when the exact one-electron Coulomb solution does not.

The exact card Hamiltonian is not a precision theory of real hydrogen.

Added physicsWhat changesCanonical route
Relativistic kinematics and spin-orbit termsFine-structure splitting and corrected energiesFine Structure
Quantized-field radiative correctionsLamb shift and radiative level shiftsLamb Shift Overview
Nuclear magnetic moment and electron–nucleus spin couplingHyperfine multipletsHyperfine Structure
Finite nuclear charge radiusShort-distance level shifts, strongest for states sampling the originHydrogen prototype and precision atomic pages
External electric fieldStark shifts and mixed parityStark Effect in Atoms
External magnetic fieldZeeman splitting and angular-momentum recouplingZeeman Effect in Atoms
Coupling to radiationTransition rates, linewidths, and finite excited-state lifetimesAtomic and light–matter interaction pages
Additional electronsScreening, exchange, correlation, and antisymmetryMany-electron atomic structure methods

The nonrelativistic expansion requires Zαfs≪1Z\alpha_{\mathrm{fs}}\ll1. The leading binding scale is of order μc2(Zαfs)2\mu c^2(Z\alpha_{\mathrm{fs}})^2, while relativistic fine-structure corrections enter at a smaller absolute scale of order μc2(Zαfs)4\mu c^2(Z\alpha_{\mathrm{fs}})^4. High-ZZ one-electron ions require increasingly relativistic and finite-nuclear-size treatment.

The static Coulomb Hamiltonian gives infinitely lived excited eigenstates. Spontaneous emission appears only after coupling the atom to quantized electromagnetic modes; it is not hidden inside the Schrödinger eigenvalue problem.

ModelChange from ordinary hydrogenMain caution
Hydrogenic IonNuclear charge Z>1Z>1Relativistic and nuclear-size effects grow with ZZ
Deuterium or tritiumDifferent nuclear mass and structureReduced-mass and isotope shifts are not the whole precision correction
PositroniumEqual electron and positron massesAnnihilation, spin, and QED effects are essential
Muonic atomElectron replaced by a muonMuch smaller orbital scale strongly probes nuclear size
Coulomb scatteringPositive-energy sectorRequires continuum and long-range scattering conventions
Alkali valence electronEffective central field with a coreQuantum defects replace exact hydrogenic degeneracy
Helium and many-electron atomsElectron–electron interaction addedNo separation into independent exact hydrogen orbitals

Hydrogen as Atomic Prototype organizes how the ideal labels and correction hierarchy are used in atomic physics.

Hydrogen is a demanding numerical benchmark because of the Coulomb singularity, the half-line radial domain, and the continuum threshold.

CheckExact target
Bound energiesEn∝−1/n2E_n\propto-1/n^2
Spatial shell degeneracyn2n^2
Ground radial maximumr=aZr=a_Z
Ground mean radius3aZ/23a_Z/2
Virial theorem2⟨T⟩=−⟨V⟩2\langle T\rangle=-\langle V\rangle
Angular eigenvaluesℏ2ℓ(ℓ+1)\hbar^2\ell(\ell+1) and ℏm\hbar m
Radial boundaryu(0)=0u(0)=0 and square-integrable decay for bound states

A radial solver should state the origin treatment, outer boundary, grid or basis, angular sector, energy zero, and convergence with resolution and box size. Radial Schrödinger Solvers gives the computational route.

  • Replacing the reduced mass by mem_e without declaring the fixed-nucleus approximation.
  • Forgetting that the card Hamiltonian describes relative motion after center-of-mass separation.
  • Mixing SI, Gaussian, natural, and atomic-unit Coulomb conventions.
  • Using the 13.6 eV13.6\,\mathrm{eV} hydrogen scale for a hydrogenic ion without the Z2Z^2 and reduced-mass factors.
  • Treating the n2n^2 spatial degeneracy as generic to all central potentials.
  • Assuming the ideal degeneracy survives fine, hyperfine, radiative, nuclear-size, or field corrections.
  • Confusing R(r)R(r), u(r)=rR(r)u(r)=rR(r), ∣R(r)∣2\lvert R(r)\rvert^2, and the radial probability ∣R(r)∣2r2dr\lvert R(r)\rvert^2r^2dr.
  • Drawing orbitals as electron trajectories or rigid material surfaces.
  • Forgetting the positive-energy Coulomb continuum.
  • Using nonrelativistic hydrogenic scaling at large ZZ without checking ZαfsZ\alpha_{\mathrm{fs}}.
  • Treating stationary excited states of the isolated Coulomb Hamiltonian as a theory of radiative lifetimes.
  • Applying the one-electron solution directly to many-electron atoms.

Starting from the two-particle kinetic energy, show why the internal kinetic term contains the reduced mass.

Solution

Define

R=mere+mNrNM,r=re−rN,\mathbf R =\frac{m_e\mathbf r_e+m_N\mathbf r_N}{M}, \qquad \mathbf r=\mathbf r_e-\mathbf r_N,

where M=me+mNM=m_e+m_N. The conjugate total and relative momenta can be chosen so that the canonical transformation gives

p^e 22me+p^N 22mN=P^ 22M+p^ 22μ,\frac{\hat{\mathbf p}_e^{\,2}}{2m_e} +\frac{\hat{\mathbf p}_N^{\,2}}{2m_N} =\frac{\hat{\mathbf P}^{\,2}}{2M} +\frac{\hat{\mathbf p}^{\,2}}{2\mu},

with

1μ=1me+1mN.\frac{1}{\mu} =\frac{1}{m_e}+\frac{1}{m_N}.

Therefore

μ=memNme+mN.\mu=\frac{m_em_N}{m_e+m_N}.

Because the Coulomb potential depends only on r\mathbf r, center-of-mass and internal motion separate.

Ignoring the small reduced-mass difference, compare the n=1n=1 and n=2n=2 energies of He+\mathrm{He}^+ with those of hydrogen. Also compare the hydrogenic length scales.

Solution

For a one-electron ion,

En(Z)=Z2En(Z=1).E_n(Z)=Z^2E_n(Z=1).

Thus the He+\mathrm{He}^+ ground state with Z=2Z=2 is four times as deeply bound as the hydrogen ground state:

E1(He+)≈−4(13.6 eV).E_1(\mathrm{He}^+) \approx-4(13.6\,\mathrm{eV}).

For n=2n=2,

E2(He+)=−4(13.6 eV)4≈−13.6 eV.E_2(\mathrm{He}^+) =-\frac{4(13.6\,\mathrm{eV})}{4} \approx-13.6\,\mathrm{eV}.

The basic Coulomb length is

aZ=aμZ,a_Z=\frac{a_\mu}{Z},

so the He+\mathrm{He}^+ orbital scale at fixed quantum numbers is half the corresponding hydrogenic scale, before reduced-mass and relativistic corrections.

Show that the ideal spatial degeneracy at fixed nn is n2n^2. Explain why this is not the observed degeneracy of precision hydrogen.

Solution

For each ℓ=0,…,n−1\ell=0,\ldots,n-1, there are 2ℓ+12\ell+1 magnetic quantum numbers. Therefore

gn=∑ℓ=0n−1(2ℓ+1)=n(n−1)+n=n2.\begin{aligned} g_n &=\sum_{\ell=0}^{n-1}(2\ell+1)\\ &=n(n-1)+n\\ &=n^2. \end{aligned}

This counts spatial states in the spinless nonrelativistic point-Coulomb model. Electron spin adds a trivial factor of two only if spin-dependent interactions are still omitted. Fine structure, the Lamb shift, hyperfine coupling, finite nuclear size, and external fields resolve or reorganize parts of the ideal degeneracy.

  • H. A. Bethe and E. E. Salpeter, Quantum Mechanics of One- and Two-Electron Atoms, Springer, 1957.
  • C. Cohen-Tannoudji, B. Diu, and F. Laloë, Quantum Mechanics, 2 vols., Wiley, 1977.
  • L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Non-Relativistic Theory, 3rd ed., Pergamon, 1977.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • D. J. Griffiths and D. F. Schroeter, Introduction to Quantum Mechanics, 3rd ed., Cambridge University Press, 2018.