Skip to content

Hyperfine structure is the splitting and mixing of atomic levels caused by electromagnetic interactions between the nucleus and the electrons. The leading terms couple the nuclear magnetic dipole moment and electric quadrupole moment to electronic fields at the nucleus. They resolve a fine-structure level of electronic angular momentum JJ into levels labeled by the total atomic angular momentum

F=I+J,\mathbf F=\mathbf I+\mathbf J,

where II is the nuclear spin.

The word hyperfine denotes a layer of structure finer than electronic fine structure. It is not another name for spin–orbit coupling: fine structure couples electronic degrees of freedom and is present even for a spinless nucleus, while hyperfine structure depends on the isotope’s nuclear spin and electromagnetic moments.

Required background. Hydrogenic states and alkali atoms supply the electronic levels. Total angular momentum and magnetic moments supply the coupling and sign conventions.

The formulas below first fix an electronic level ∣γJ⟩|\gamma J\rangle and resolve its nuclear-spin coupling. This is justified when hyperfine matrix elements are small compared with neighboring electronic separations. Otherwise off-diagonal hyperfine interactions mix different JJ levels and must also be included.

The nucleus has angular momentum operators

I2∣ImI⟩=ℏ2I(I+1)∣ImI⟩,Iz∣ImI⟩=ℏmI∣ImI⟩.\begin{aligned} I^2|I m_I\rangle &= \hbar^2I(I+1)|I m_I\rangle, \\ I_z|I m_I\rangle &= \hbar m_I|I m_I\rangle. \end{aligned}

Its magnetic dipole moment is conventionally written

μI=gIμNIℏ,μN=eℏ2mp,\boldsymbol\mu_I = g_I\mu_N\frac{\mathbf I}{\hbar}, \qquad \mu_N = \frac{e\hbar}{2m_p},

where gIg_I is the nuclear gg factor and μN\mu_N is the nuclear magneton. Nuclear magnetic moments are typically smaller than electronic moments because μN/μB=me/mp\mu_N/\mu_B=m_e/m_p.

The magnetic moment need not be parallel to I\mathbf I with a positive coefficient; the sign is carried by gIg_I. Isotopes of the same element can have different II, gIg_I, and quadrupole moment, so their hyperfine spectra can differ even when their electronic configurations are nearly identical.

A nucleus with I≥1I\ge1 may possess a nonspherical electric charge distribution characterized by a spectroscopic quadrupole moment QQ. A spin-00 or spin-1/21/2 nucleus cannot carry a permanent rank-2 moment in a stationary angular-momentum eigenstate.

The measured electric-quadrupole hyperfine constant is not QQ alone. It is proportional to the product of the nuclear quadrupole moment and the electronic electric-field gradient at the nucleus. Extracting one from the other requires either independent nuclear information or an atomic-structure calculation.

Higher nuclear multipoles, such as a magnetic octupole moment, can produce still smaller terms. They matter in selected precision studies but are outside the standard AA–BB model developed here.

Because I\mathbf I and J\mathbf J act on nuclear and electronic factors, the uncoupled basis is

∣ImI⟩⊗∣γJmJ⟩.|I m_I\rangle\otimes|\gamma J m_J\rangle.

For a rotationally invariant zero-field hyperfine Hamiltonian, the coupled basis is more useful:

∣γ;IJ;FmF⟩,|\gamma;I J;F m_F\rangle,

with

F=∣I−J∣,∣I−J∣+1,…,I+JF=|I-J|,|I-J|+1,\ldots,I+J

and

mF=−F,−F+1,…,F.m_F=-F,-F+1,\ldots,F.

The state transformation is a Clebsch–Gordan expansion,

∣IJ;FmF⟩=∑mI,mJCImI,JmJFmF×∣ImI⟩∣JmJ⟩.|I J;F m_F\rangle = \sum_{m_I,m_J} C^{F m_F}_{I m_I,J m_J} \times |I m_I\rangle |J m_J\rangle.

The dimensions agree:

(2I+1)(2J+1)=∑F=∣I−J∣I+J(2F+1).(2I+1)(2J+1) = \sum_{F=|I-J|}^{I+J}(2F+1).

At zero external field, rotational symmetry makes all 2F+12F+1 values of mFm_F degenerate. Hyperfine coupling correlates nuclear and electronic projections; it does not in general preserve mIm_I and mJm_J separately.

Hyperfine doublet with I equals three halves and J equals one half, showing the degeneracy-weighted centroid

For I=3/2I=3/2, J=1/2J=1/2, and A>0A>0, the fivefold F=2F=2 level lies at 3A/43A/4 and the threefold F=1F=1 level at −5A/4-5A/4. Their weighted centroid is zero, although their midpoint is not. Generally E+−E−=A(I+1/2)E_+-E_-=A(I+1/2); negative AA reverses the order, and the transition energy is the absolute difference.

Inside an isolated fine-structure level, rotational symmetry makes the leading magnetic-dipole interaction a scalar. With AA defined as an energy,

HM1=Aℏ2I⋅J.H_{\mathrm{M1}} = \frac{A}{\hbar^2} \mathbf I\cdot\mathbf J.

The angular identity

2I⋅J=F2−I2−J22\mathbf I\cdot\mathbf J = F^2-I^2-J^2

immediately gives

EFM1=A2K,E_F^{\mathrm{M1}} = \frac{A}{2}K,

where

K≡F(F+1)−I(I+1)−J(J+1).K \equiv F(F+1)-I(I+1)-J(J+1).

The constant AA includes the nuclear magnetic moment and the electronic magnetic field at the nucleus for the chosen state. It may be positive or negative. The formula predicts angular dependence after AA is known; it does not calculate AA from first principles.

For adjacent FF manifolds in the pure magnetic-dipole model,

EF−EF−1=A2[KF−KF−1],=AF.\begin{aligned} E_F-E_{F-1} &= \frac{A}{2} \left[K_F-K_{F-1}\right], \\ &=AF. \end{aligned}

Thus adjacent intervals are proportional to the larger FF. Electric-quadrupole terms, higher multipoles, second-order mixing, and external fields violate this simple interval rule.

The degeneracy-weighted magnetic-dipole centroid is unchanged:

∑F(2F+1)EFM1=0.\sum_F(2F+1)E_F^{\mathrm{M1}}=0.

This trace statement is useful when reconstructing a fine-structure centroid from resolved hyperfine levels. An unweighted average is generally wrong.

The effective constant AA can contain several electronic interactions:

Contribution Electronic structure Typical importance
orbital magnetic field electron orbital current coupled to μI\boldsymbol\mu_I states with nonzero orbital angular momentum
spin-dipolar field anisotropic dipole–dipole coupling of electron and nucleus non-SS density and open shells
Fermi contact term electron spin density at the nucleus dominant for many SS states
relativistic and correlation terms small components, core polarization, configuration mixing essential in heavy and many-electron atoms

For a one-electron SS state, the leading contact interaction can be written

HF=2μ03geμBgIμNδ3(r)I⋅Sℏ2,H_{\mathrm F} = \frac{2\mu_0}{3} g_e\mu_Bg_I\mu_N \delta^3(\mathbf r) \frac{\mathbf I\cdot\mathbf S}{\hbar^2},

where geg_e denotes the positive magnitude convention for the electron spin gg factor. Its first-order constant is proportional to

AnS(0)∝gIμNgeμB∣ψnS(0)∣2.A_{nS}^{(0)} \propto g_I\mu_N g_e\mu_B |\psi_{nS}(0)|^2.

This explains the strong contact sensitivity of hydrogenic and alkali SS states. It does not justify inserting a hydrogenic wavefunction into a heavy neutral atom. Core polarization, exchange, relativistic penetration, finite magnetization, and electronic correlation alter the spin density at the nucleus.

The electric-quadrupole interaction couples rank-2 nuclear and electronic tensors. Schematically,

HE2=∑q=−22(−1)qTq(2)(nuc)T−q(2)(el).H_{\mathrm{E2}} = \sum_{q=-2}^{2} (-1)^q T_q^{(2)}(\mathrm{nuc}) T_{-q}^{(2)}(\mathrm{el}).

Within an isolated I,JI,J manifold, the conventional first-order energy is

EFE2=BNIJFDIJ,NIJF=34K(K+1)−I(I+1)J(J+1),DIJ=2I(2I−1)J(2J−1).\begin{aligned} E_F^{\mathrm{E2}} &= B\frac{N_{IJF}}{D_{IJ}}, \\ N_{IJF} &= \frac34K(K+1) -I(I+1)J(J+1), \\ D_{IJ} &= 2I(2I-1)J(2J-1). \end{aligned}

The BB term exists only when both angular momenta can support the required rank-2 coupling:

I≥1,J≥1.I\ge1, \qquad J\ge1.

For J=1/2J=1/2, the electronic quadrupole expectation vanishes, even if the nucleus has I≥1I\ge1. This is why alkali nS1/2nS_{1/2} and nP1/2nP_{1/2} levels have no first-order electric-quadrupole BB term, whereas nP3/2nP_{3/2} levels can.

The standard dipole-plus-quadrupole model is therefore

EFhfs=A2K+EFE2.E_F^{\mathrm{hfs}} = \frac{A}{2}K +E_F^{\mathrm{E2}}.

Different databases may tabulate AA and BB in frequency units, writing A/hA/h and B/hB/h implicitly. Always inspect the convention before substituting a number into an energy formula.

The hydrogen 1S1/21S_{1/2} ground state has

I=12,J=12,I=\frac12, \qquad J=\frac12,

so the allowed total angular momenta are F=0F=0 and F=1F=1. For the pure dipole Hamiltonian,

EF=1=A4,EF=0=−3A4.\begin{aligned} E_{F=1} &=\frac{A}{4}, \\ E_{F=0} &=-\frac{3A}{4}. \end{aligned}

Their difference is

ΔE1Shfs=A.\Delta E_{1S}^{\mathrm{hfs}}=A.

For ordinary hydrogen, A>0A>0, so the F=1F=1 triplet lies above the F=0F=0 singlet. The evaluated transition frequency quoted by NIST is

ν21=1420.405 751 766 7(10) MHz.\nu_{21} = 1420.405\,751\,766\,7(10) \ \mathrm{MHz}.

Using λ=c/ν\lambda=c/\nu gives approximately

λ21≃21.106 cm,\lambda_{21}\simeq21.106\ \mathrm{cm},

the celebrated 21 cm line. It is a weak magnetic-dipole transition, but the vast amount of neutral hydrogen in astrophysical environments makes it an exceptionally important probe.

The contact model gives the leading scale, not the final prediction. Precision hydrogen theory includes:

  • the electron anomalous magnetic moment and other QED corrections;
  • relativistic and recoil terms;
  • finite proton magnetization and charge distributions;
  • the Zemach correction, which couples charge and magnetic form factors;
  • proton polarizability and other two-photon-exchange effects;
  • weak and hadronic vacuum-polarization contributions at much smaller scales.

The measured interval is extraordinarily precise, while theoretical interpretation is limited more strongly by proton structure. Hyperfine spectroscopy therefore probes a different nuclear combination from the Lamb shift: the Lamb shift is primarily sensitive to charge-radius effects, whereas hyperfine structure is especially sensitive to magnetic structure and spin-dependent two-photon exchange.

The words singlet and triplet describe the coupling of electron and proton spins here. The particles are distinguishable, so this terminology does not impose the exchange-symmetry constraints that apply to two identical electrons.

An alkali ground state is approximately nS1/2nS_{1/2}, so J=1/2J=1/2. If I≠0I\ne0, it splits into

F−=I−12,F+=I+12.F_-=I-\frac12, \qquad F_+=I+\frac12.

The dipole energies are

EF+=AI2,EF−=−A(I+1)2,\begin{aligned} E_{F_+} &=\frac{AI}{2}, \\ E_{F_-} &=-\frac{A(I+1)}{2}, \end{aligned}

so the signed difference EF+−EF−E_{F_+}-E_{F_-} is

ΔEhfs=A(I+12).\Delta E_{\mathrm{hfs}} = A\left(I+\frac12\right).

Representative microwave scales are:

Isotope and state Nuclear spin Hyperfine levels Ground interval
1^1H 1S1/21S_{1/2} 1/21/2 F=0,1F=0,1 1.4204057517667 GHz1.4204057517667\ \mathrm{GHz}
85^{85}Rb 5S1/25S_{1/2} 5/25/2 F=2,3F=2,3 about 3.036 GHz3.036\ \mathrm{GHz}
87^{87}Rb 5S1/25S_{1/2} 3/23/2 F=1,2F=1,2 about 6.835 GHz6.835\ \mathrm{GHz}
133^{133}Cs 6S1/26S_{1/2} 7/27/2 F=3,4F=3,4 exactly 9.192631770 GHz9.192631770\ \mathrm{GHz} as an SI defining frequency

The rubidium entries are rounded orientation values. For numerical hyperfine constants, use isotope- and state-specific measurements compiled by Allegrini, Arimondo, and Orozco or the convention-consistent data sheets by Steck; a general electronic-level database need not resolve hyperfine components.

For nP1/2nP_{1/2}, J=1/2J=1/2 again eliminates the first-order BB term. For nP3/2nP_{3/2}, J=3/2J=3/2 permits electric-quadrupole structure when I≥1I\ge1. Optical DD lines therefore split into multiple F→F′F\to F' components with isotope-dependent spacings and line strengths.

The contact interaction is no longer the entire story in excited or non-SS states. Orbital, spin-dipolar, relativistic, and many-body contributions can dominate. Hyperfine constants are stringent tests of electronic wavefunctions near the nucleus and of core-polarization calculations.

For identical electronic structure and pointlike nuclei, AA is proportional to gI=μI/(IμN)g_I=\mu_I/(I\mu_N), not simply to the tabulated magnetic moment μI\mu_I when isotopes have different spins. Finite nuclear magnetization and charge distributions produce departures from the gIg_I ratio known as hyperfine anomalies; ordinary mass-dependent electronic changes must also be controlled. They are physical isotope-dependent corrections, not necessarily inconsistent measurements.

Within one electronic state, transitions between hyperfine manifolds are typically driven by microwave magnetic fields. In the weak-field coupled basis, the leading magnetic-dipole rules are

ΔF=0,±1,F=0↮F′=0,\Delta F=0,\pm1, \qquad F=0\not\leftrightarrow F'=0,

with

ΔmF=0,±1\Delta m_F=0,\pm1

selected by the microwave polarization relative to the quantization axis. These are angular selection rules, not guarantees of equal line strength.

For optical electric-dipole transitions between hyperfine-resolved electronic levels, the same FF and mFm_F rules appear after the electronic reduced matrix element is recoupled to nuclear spin. Atomic Selection Rules gives the operational 6j6j recoupling and polarization rules; their general derivation belongs to the Wigner–Eckart Theorem and Selection Rules pages.

Hyperfine transitions can be exceptionally narrow because the two states share the same electronic configuration and the transition is not an allowed optical electric-dipole decay. Narrowness alone does not make a transition an ideal clock: it must also be preparable, interrogable, reproducible, and insensitive to environmental perturbations.

For 133^{133}Cs, I=7/2I=7/2 and the 6S1/26S_{1/2} ground state has F=3F=3 and F=4F=4. The microwave clock transition uses the states

∣F=3,mF=0⟩⟷∣F=4,mF=0⟩.|F=3,m_F=0\rangle \longleftrightarrow |F=4,m_F=0\rangle.

At weak magnetic field, a hyperfine Zeeman sublevel has the leading shift

ΔEZ(1)≃gFμBmFB.\Delta E_Z^{(1)} \simeq g_F\mu_Bm_FB.

Each clock state has mF=0m_F=0, so its diagonal first-order shift vanishes. This makes the transition first-order insensitive near zero field. It is not magnetically insensitive to all orders: the field mixes the two mF=0m_F=0 hyperfine states and produces a quadratic Zeeman shift.

The SI second is defined by fixing the numerical value of the unperturbed 133^{133}Cs ground-state hyperfine transition frequency:

ΔνCs=9 192 631 770 Hz\Delta\nu_{\mathrm{Cs}} = 9\,192\,631\,770\ \mathrm{Hz}

exactly. The word unperturbed is crucial. A laboratory clock measures atoms in finite fields, at finite temperature, with collisions, motion, microwave interrogation, and gravitational potential. A primary frequency standard evaluates and corrects its realization toward the defining ideal.

For a cesium standard, the main corrections include quadratic Zeeman, blackbody Stark, collisional, and microwave interrogation shifts. Comparing geographically separated clocks also requires gravitational-redshift corrections. Atomic Clocks explains how interrogation and the apparatus determine the measured frequency.

In a static axially symmetric magnetic field, mFm_F remains exact while FF is an approximate label when the Zeeman interaction is weak compared with the hyperfine interval. A practical comparison is

gJμBB≪ΔEhfs.g_J\mu_BB \ll \Delta E_{\mathrm{hfs}}.

In an intermediate field, states with the same conserved mF=mI+mJm_F=m_I+m_J mix. One should diagonalize

H=Hhfs−μJ⋅B−μI⋅BH = H_{\mathrm{hfs}} -\boldsymbol\mu_J\cdot\mathbf B -\boldsymbol\mu_I\cdot\mathbf B

in fixed-mFm_F blocks. In the strong-field hyperfine Paschen–Back regime, the uncoupled labels mI,mJm_I,m_J become more useful. Zeeman Effect in Atoms derives the Breit–Rabi crossover and polarization-resolved spectra; Zeeman Effect as a Perturbation Example develops the general perturbative field logic.

Avoid assigning a single field-independent gFg_F through this crossover. The eigenvectors change with BB, and avoided crossings occur between states with the same exact symmetries.

A J=1/2J=1/2 doublet determines the magnitude of AA from one interval; assigning the FF levels also determines its sign. When both AA and BB contribute, fit several levels with the complete formula. Reconstruct the electronic centroid with degeneracy weights, rather than averaging observed optical peak positions: peaks can share different upper and lower endpoints and have population-dependent intensities.

The first-order AA–BB model presumes an isolated JJ level. Nearby electronic levels produce second-order shifts through off-diagonal hyperfine matrix elements, so precision fits must distinguish those shifts from effective fitted constants. Fields, optical pumping, collisions, and unresolved Zeeman components can also alter the measured resonance.

1. Separate dipole and quadrupole constants

Section titled “1. Separate dipole and quadrupole constants”

For I=J=1I=J=1, measured signed intervals are (E1−E0)/h=5  MHz(E_1-E_0)/h=5\;\mathrm{MHz} and (E2−E1)/h=46  MHz(E_2-E_1)/h=46\;\mathrm{MHz}. Neglect higher-order terms. Find A/hA/h and B/hB/h, and verify the level centroid.

Solution

For F=0,1,2F=0,1,2, the values of KK are −4,−2,2-4,-2,2. The level formula gives

E0=−2A+52B,E1=−A−54B,E2=A+14B.E_0=-2A+\frac52B,\quad E_1=-A-\frac54B,\quad E_2=A+\frac14B.

Writing d1=E1−E0d_1=E_1-E_0 and d2=E2−E1d_2=E_2-E_1 gives

B=d2−2d19,A=d1+154B.B=\frac{d_2-2d_1}{9},\qquad A=d_1+\frac{15}{4}B.

Thus B/h=4  MHzB/h=4\;\mathrm{MHz}, A/h=20  MHzA/h=20\;\mathrm{MHz}, and the three levels are EF/h=(−30,−25,21)  MHzE_F/h=(-30,-25,21)\;\mathrm{MHz}. Their weighted sum is 1(−30)+3(−25)+5(21)=01(-30)+3(-25)+5(21)=0. Neither measured gap alone equals A/hA/h.

For general nuclear spin I≥1/2I\ge1/2 and J=1/2J=1/2, use EF=AK/2E_F=AK/2 to derive the two energies, their separation, and their degeneracy-weighted centroid.

Solution

For F+=I+1/2F_+=I+1/2,

K+=I,E+=AI2.K_+=I, \qquad E_+=\frac{AI}{2}.

For F−=I−1/2F_-=I-1/2,

K−=−(I+1),E−=−A(I+1)2.K_-=-(I+1), \qquad E_-=-\frac{A(I+1)}{2}.

Therefore

E+−E−=A(I+12).E_+-E_- = A\left(I+\frac12\right).

The degeneracies are 2I+22I+2 and 2I2I. Their weighted sum is

(2I+2)E++(2I)E−=AI(I+1)−AI(I+1),=0.\begin{aligned} (2I+2)E_+ +(2I)E_- &=AI(I+1)-AI(I+1), \\ &=0. \end{aligned}

Thus the pure dipole interaction splits the level about its degeneracy-weighted electronic centroid.

Decide whether a first-order electric-quadrupole hyperfine constant BB can occur for each pair: (a) I=1/2,J=3/2I=1/2,J=3/2; (b) I=3/2,J=1/2I=3/2,J=1/2; (c) I=3/2,J=3/2I=3/2,J=3/2; (d) I=0,J=2I=0,J=2. Explain each answer without evaluating the energy formula.

Solution

A diagonal rank-2 coupling requires both I≥1I\ge1 and J≥1J\ge1.

Case BB term? Reason
(a) no the spin-1/21/2 nucleus has no permanent quadrupole moment
(b) no a J=1/2J=1/2 electronic manifold cannot support the required rank-2 expectation
(c) yes both angular momenta are at least one
(d) no an I=0I=0 nucleus has neither magnetic-dipole nor electric-quadrupole orientation

Hyperfine mixing with other JJ levels can produce higher-order effects even when a first-order diagonal BB constant vanishes, but that does not change this selection rule.

Using c=2.99792458×108 m s−1c=2.99792458\times10^8\ \mathrm{m\,s^{-1}} and ν=1.4204057518×109 s−1\nu=1.4204057518\times10^9\ \mathrm{s^{-1}}, compute the vacuum wavelength of the hydrogen ground-state hyperfine transition.

Solution

The wavelength is

λ=cν,≃2.99792458×1081.4204057518×109 m,≃0.21106 m.\begin{aligned} \lambda &=\frac{c}{\nu}, \\ &\simeq \frac{2.99792458\times10^8} {1.4204057518\times10^9} \ \mathrm m, \\ &\simeq0.21106\ \mathrm m. \end{aligned}

Thus

λ≃21.106 cm.\lambda\simeq21.106\ \mathrm{cm}.

The nickname “21 cm line” is a rounded wavelength label; precision work uses the transition frequency and its stated convention.

5. Why clock states still need a bias field

Section titled “5. Why clock states still need a bias field”

The caesium clock transition uses mF=0m_F=0 states, for which the first-order Zeeman formula vanishes. Why is a small bias magnetic field still useful, and why must its magnitude be measured?

Solution

A bias field defines a quantization axis and separates the mF≠0m_F\ne0 Zeeman components. That makes state preparation, polarization control, and identification of the clock resonance more robust.

The field cannot be ignored because it mixes the two mF=0m_F=0 hyperfine states and produces a quadratic Zeeman shift. Spatial inhomogeneity also makes different atomic trajectories sample different shifts. A primary standard therefore measures the field, often using field-sensitive Zeeman transitions, and corrects the clock frequency to the unperturbed limit with an uncertainty.

  • M. Allegrini, E. Arimondo, and L. A. Orozco, “Survey of Hyperfine Structure Measurements in Alkali Atoms,” Journal of Physical and Chemical Reference Data 51, 043102 (2022), doi:10.1063/5.0098061.
  • E. Arimondo, M. Inguscio, and P. Violino, “Experimental Determinations of the Hyperfine Structure in the Alkali Atoms,” Reviews of Modern Physics 49, 31–75 (1977), doi:10.1103/RevModPhys.49.31.
  • G. Breit and I. I. Rabi, “Measurement of Nuclear Spin,” Physical Review 38, 2082–2083 (1931), doi:10.1103/PhysRev.38.2082.2.
  • Bureau International des Poids et Mesures, SI Base Unit: Second, defining ΔνCs=9 192 631 770 Hz\Delta\nu_{\mathrm{Cs}}=9\,192\,631\,770\ \mathrm{Hz} exactly.
  • C. E. Carlson, V. Nazaryan, and K. Griffioen, “Proton Structure Corrections to Electronic and Muonic Hydrogen Hyperfine Splitting,” Physical Review A 78, 022517 (2008), doi:10.1103/PhysRevA.78.022517.
  • E. Fermi, “Über die magnetischen Momente der Atomkerne,” Zeitschrift für Physik 60, 320–333 (1930), doi:10.1007/BF01339933.
  • C. J. Foot, Atomic Physics (Oxford University Press, 2005), Chapters 2, 6, and 7.
  • A. Kramida, Y. Ralchenko, J. Reader, and the NIST ASD Team, NIST Atomic Spectra Database, National Institute of Standards and Technology, doi:10.18434/T4W30F.
  • National Institute of Standards and Technology, Atomic Spectroscopy: Atomic States, Shells, and Configurations, including the evaluated 1^1H ground-state hyperfine interval.
  • I. I. Sobelman, Atomic Spectra and Radiative Transitions, 2nd ed. (Springer, 1992), doi:10.1007/978-3-642-76907-8.
  • D. A. Steck, Alkali D Line Data, isotope-resolved reference sheets with constants, conventions, and source histories.
  • G. K. Woodgate, Elementary Atomic Structure, 2nd ed. (Oxford University Press, 1980), Chapters 9–12.