Skip to content

Hydrogen Atom

The hydrogen atom is the canonical exactly solvable three-dimensional bound-state problem. In the nonrelativistic approximation, the electron and proton interact through an attractive Coulomb potential. After separating the center-of-mass motion, the relative coordinate obeys a one-body central-potential problem.

Hydrogen is important because it is both concrete and structurally rich: it introduces principal, orbital, and magnetic quantum numbers; it connects spherical harmonics to radial special functions; and it exposes degeneracies that later corrections split.

This page treats the spinless nonrelativistic Coulomb Hamiltonian for a point electron and point proton, with reduced mass. It does not include fine structure, hyperfine structure, the Lamb shift, finite nuclear size, external fields, radiative transitions, or many-electron effects.

Those corrections are physically essential for precision spectroscopy, but they should not be folded into the basic Coulomb solution. The purpose here is to solve the standard Schrödinger eigenvalue problem cleanly.

Let the electron and proton masses be mem_e and mpm_p. The reduced mass is

μ=mempme+mp.\mu = \frac{m_em_p}{m_e+m_p}.

After separating the free center-of-mass motion, the relative-coordinate Hamiltonian is

H^=−ℏ22μ∇2−e24πϵ0r.\hat H = -\frac{\hbar^2}{2\mu}\nabla^2 - \frac{e^2}{4\pi\epsilon_0r}.

The stationary Schrödinger equation is

H^ψ(r)=Eψ(r).\hat H\psi(\mathbf r)=E\psi(\mathbf r).

It is useful to introduce the reduced-mass Bohr radius

a0=4πϵ0ℏ2μe2.a_0 = \frac{4\pi\epsilon_0\hbar^2}{\mu e^2}.

If one uses the electron mass instead of μ\mu, this is the infinite-proton-mass approximation. The numerical difference is small for hydrogen but conceptually important.

Because the potential depends only on rr, write

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

The angular quantum numbers satisfy

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

The principal quantum number is

n=1,2,3,….n=1,2,3,\ldots.

The angular meaning of the n,ℓ,mn,\ell,m labels, parity, and shell multiplets is developed in Hydrogen Atom Angular Structure.

Within the spinless bound-state sector, the commuting observables HH, L2L^2, and LzL_z distinguish the states by nn, ℓ\ell, and mm. Their role as a sector-dependent complete set is explained in Complete Sets of Commuting Observables.

The radial equation is the Coulomb specialization of the general radial equation:

−ℏ22μd2udr2+[−e24πϵ0r+ℏ2ℓ(ℓ+1)2μr2]u=Eu,u=rR.-\frac{\hbar^2}{2\mu}\frac{d^2u}{dr^2} + \left[ -\frac{e^2}{4\pi\epsilon_0r} + \frac{\hbar^2\ell(\ell+1)}{2\mu r^2} \right]u =Eu, \qquad u=rR.

Normalizable bound-state solutions exist only for the allowed combinations of nn and ℓ\ell above. The angular part is supplied by spherical harmonics, and the radial polynomial part is supplied by Laguerre polynomials.

The bound-state energies are

En=−μe42(4πϵ0)2ℏ21n2=−ℏ22μa021n2.E_n = -\frac{\mu e^4}{2(4\pi\epsilon_0)^2\hbar^2} \frac{1}{n^2} = -\frac{\hbar^2}{2\mu a_0^2} \frac{1}{n^2}.

For ordinary hydrogen this is approximately

En≈−13.6 eVn2,E_n\approx-\frac{13.6\,\mathrm{eV}}{n^2},

with the reduced-mass correction included in high-precision values.

The energy depends only on nn, not on ℓ\ell or mm, in the ideal nonrelativistic Coulomb problem. The mm degeneracy follows from rotational symmetry for any central potential. The additional ℓ\ell degeneracy is special to the Coulomb potential and is explained in Degeneracy of the Hydrogen Atom.

For a fixed nn, the number of spatial bound states is

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

If spin is included without spin-dependent interactions, this doubles to 2n22n^2. Fine structure, hyperfine structure, Lamb-shift physics, and external fields split this ideal degeneracy.

The detailed reference for the functions in this section is Radial Wavefunctions.

With

ρ=2rna0,\rho=\frac{2r}{na_0},

one common normalized convention is

Rnℓ(r)=(2na0)3(n−ℓ−1)!2n[(n+ℓ)!]e−ρ/2ρℓLn−ℓ−1(2ℓ+1)(ρ).R_{n\ell}(r) = \sqrt{ \left( \frac{2}{na_0} \right)^3 \frac{(n-\ell-1)!}{2n[(n+\ell)!]} } e^{-\rho/2} \rho^\ell L_{n-\ell-1}^{(2\ell+1)}(\rho).

The full wavefunction is RnℓYℓmR_{n\ell}Y_\ell^m. The normalization convention is

∫0∞∣Rnℓ(r)∣2r2 dr=1,∫S2∣Yℓm∣2 dΩ=1.\int_0^\infty \lvert R_{n\ell}(r)\rvert^2r^2\,dr =1, \qquad \int_{S^2} \lvert Y_\ell^m\rvert^2\,d\Omega =1.

Equivalently, if unℓ=rRnℓu_{n\ell}=rR_{n\ell}, then

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

The ground state has n=1n=1, ℓ=0\ell=0, and m=0m=0. Its radial and full wavefunctions are

R10(r)=2a0−3/2e−r/a0,ψ100(r,θ,ϕ)=1πa03e−r/a0.R_{10}(r) = 2a_0^{-3/2}e^{-r/a_0}, \qquad \psi_{100}(r,\theta,\phi) = \frac{1}{\sqrt{\pi a_0^3}}e^{-r/a_0}.

The wavefunction is largest at the origin, but the radial probability density includes the volume factor r2r^2:

P10(r)=∣R10(r)∣2r2=4r2a03e−2r/a0.P_{10}(r) = \lvert R_{10}(r)\rvert^2r^2 = \frac{4r^2}{a_0^3}e^{-2r/a_0}.

This radial probability density is maximal at r=a0r=a_0. This is one precise sense in which the Bohr radius sets the size of the ground-state atom.

Hydrogenic orbitals are stationary one-electron wavefunctions, not classical electron trajectories. The labels s,p,d,…s,p,d,\ldots correspond to ℓ=0,1,2,…\ell=0,1,2,\ldots. Complex orbitals are natural angular-momentum eigenstates; real orbitals are particular linear combinations within a degenerate mm subspace. The notation, real-versus-complex basis choice, nodal surfaces, and visualization cautions are developed in Atomic Orbitals.

For fixed nn and ℓ\ell, the number of radial nodes is

n−ℓ−1.n-\ell-1.

Angular nodes come from the spherical harmonics. The total nodal structure is therefore a product of radial and angular information. A separate page can treat orbital visualization in more detail; this page uses orbitals only as labels for the exact eigenfunctions.

The Coulomb potential supplies an infinite ladder of bound states accumulating at E=0E=0 from below. States with larger nn are less tightly bound and extend to larger radii. The continuum above E=0E=0 describes ionized scattering states; see Continuum States of the Coulomb Problem for the positive-energy sector.

The quantum numbers have distinct meanings:

  • nn controls the energy and overall radial scale in the ideal Coulomb problem.
  • ℓ\ell controls orbital angular momentum and the centrifugal barrier.
  • mm controls the projection of orbital angular momentum along the chosen zz axis.

The zz axis is a convention unless an external field or measurement apparatus selects a direction.

  • Using mem_e without noticing the reduced-mass approximation.
  • Treating orbitals as paths followed by an electron.
  • Forgetting the r2r^2 measure when interpreting radial probabilities.
  • Assuming the ideal n2n^2 degeneracy survives all physical corrections.
  • Confusing the radial wavefunction RnℓR_{n\ell} with the radial probability density ∣Rnℓ∣2r2\lvert R_{n\ell}\rvert^2r^2.
  • Thinking that mm changes the radial equation in a rotationally invariant Coulomb problem.
  1. Normalize the ground-state radial wavefunction R10(r)=2a0−3/2e−r/a0R_{10}(r)=2a_0^{-3/2}e^{-r/a_0} using the radial measure.
Solution

Compute

∫0∞∣R10(r)∣2r2 dr=4a03∫0∞r2e−2r/a0 dr.\int_0^\infty \lvert R_{10}(r)\rvert^2r^2\,dr = \frac{4}{a_0^3} \int_0^\infty r^2e^{-2r/a_0}\,dr.

Using

∫0∞r2e−αr dr=2α3,\int_0^\infty r^2e^{-\alpha r}\,dr = \frac{2}{\alpha^3},

with α=2/a0\alpha=2/a_0, the integral becomes

4a032(2/a0)3=1.\frac{4}{a_0^3} \frac{2}{(2/a_0)^3} =1.
  1. Show that the spatial degeneracy of a fixed principal shell is n2n^2.
Solution

For fixed nn, the allowed orbital quantum numbers are ℓ=0,1,…,n−1\ell=0,1,\ldots,n-1. For each ℓ\ell, there are 2ℓ+12\ell+1 allowed values of mm. Therefore

gn=∑ℓ=0n−1(2ℓ+1)=2∑ℓ=0n−1ℓ+n=n(n−1)+n=n2.g_n = \sum_{\ell=0}^{n-1}(2\ell+1) = 2\sum_{\ell=0}^{n-1}\ell +n = n(n-1)+n = n^2.
  1. Find the most probable radius in the ground state.
Solution

The radial probability density is

P10(r)=4r2a03e−2r/a0.P_{10}(r) = \frac{4r^2}{a_0^3}e^{-2r/a_0}.

Its derivative is proportional to

ddr(r2e−2r/a0)=2re−2r/a0(1−ra0).\frac{d}{dr} \left( r^2e^{-2r/a_0} \right) = 2re^{-2r/a_0} \left( 1-\frac{r}{a_0} \right).

The critical points are r=0r=0 and r=a0r=a_0. The nonzero maximum occurs at

r=a0.r=a_0.
  • D. J. Griffiths and D. F. Schroeter, Introduction to Quantum Mechanics, 3rd ed., Cambridge University Press, 2018.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • C. Cohen-Tannoudji, B. Diu, and F. Laloë, Quantum Mechanics, Wiley, 1977.
  • L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Non-Relativistic Theory, 3rd ed., Pergamon, 1977.