Kramers' Opacity Law

Kramers' opacity law is a foundational quantum-mechanical scaling law in stellar astrophysics that describes how the Rosseland mean radiative opacity κ of an ionized stellar gas depends on mass density ρ, temperature T, and chemical composition when photon absorption is dominated by bound-free transitions (photoionization of bound atomic electrons) and free-free transitions (inverse bremsstrahlung absorption by free electrons passing near atomic nuclei). Derived in 1923 by Dutch theoretical physicist Hendrik Anthony Kramers using Niels Bohr's correspondence principle and semi-classical quantum electrodynamics, and introduced into stellar structure theory by Sir Arthur Stanley Eddington in 1926, the law establishes that radiative opacity scales directly with density and inversely with the 3.5 power of temperature:

κ = κ₀ ρ T^−7/2

Across astronomy, Kramers' opacity law governs radiative energy transport inside intermediate- and solar-mass main-sequence stars, triggers the radiative core nucleation that carries contracting pre-main-sequence stars from the Hayashi track onto the Henyey track, stabilizes radiative stellar envelopes against thermal runaways, and explains why partial ionization zones of hydrogen and helium invert the temperature exponent to drive the κ-mechanism pulsations of Cepheid and RR Lyrae variable stars along the Period-Luminosity relation.

Free-Free (Inverse Bremsstrahlung) and Bound-Free (Photoionization) Cross-Sections

In a stellar interior at temperatures between 10⁴ K and 3 × 10⁷ K, a photon of frequency ν can be absorbed by an electron through two primary coulombic interactions with an ion of nuclear charge Ze:

1. Free-Free Absorption (κ_ff, Inverse Bremsstrahlung): A free thermal electron of velocity v passing within impact parameter b of an ion of charge Ze absorbs a photon of energy hν to jump to a higher hyperbolic continuum orbit. By Kramers' semi-classical electrodynamic calculation, the absorption cross-section per ion–electron encounter scales with ion charge squared (Z²), inversely with electron velocity (v⁻¹), and inversely with photon frequency cubed (ν⁻³). Multiplying by the number density of target ions (n_i ∝ ρ) and the phase-space density of free thermal electrons (n_e ∝ ρ) having Maxwell–Boltzmann thermal speed ⟨v⟩ ∝ (k_BT / m_e)^1/2 yields the monochromatic mass absorption coefficient:
   • κ_ff,ν ∝ (n_i n_e / ρ) Z² ⟨v⟩⁻¹ ν⁻³ g_ff(ν, T) ∝ ρ T^−1/2 ν⁻³ g_ff(ν, T)
2. Bound-Free Absorption (κ_bf, Photoionization): A photon with energy hν exceeding the ionization potential χ_n = Z² R_∞ / n² ejects an electron bound in principal quantum shell n into the continuum. The Kramers photoionization cross-section for hydrogenic shell n is σ_bf,n(ν) = (64π⁴ m_e e¹⁰ Z⁴) / (3√3 c h⁶ n⁵ ν³) g_bf(ν, n), which likewise scales as ν⁻³ above the ionization edge. Because the Saha ionization equation dictates that the fraction of transiently recombined bound states in a mostly ionized gas scales proportionally to electron density and inversely with the thermal partition function (N_bound / N_ion ∝ n_e T^−3/2 exp(χ_n / k_BT) ∝ ρ T^−3/2), multiplying the bound-state population fraction by the ν⁻³ cross-section and correcting for stimulated emission [1 − exp(−hν / k_BT)] ≈ hν / (k_BT) yields the same monochromatic scaling:
   • κ_bf,ν ∝ ρ T^−1/2 ν⁻³ g_bf(ν, T)

Here g_ff and g_bf are the quantum-mechanical Gaunt factors (introduced by John Arthur Gaunt in 1930), which are dimensionless correction factors of order unity (g ≈ 0.8–1.2) that account for exact wave-mechanical Coulomb wavefunctions relative to Kramers' semi-classical parabolic trajectory approximation.

Integration to the Rosseland Mean Opacity

Deep inside a stellar radiative zone (far below the photosphere governed by the Eddington-Barbier relation), the photon mean free path ℓ = (κ ρ)⁻¹ is millimetres to centimetres—roughly 10⁻¹¹ of the stellar radius. Under this diffusion regime, radiative energy transport depends on the harmonic Rosseland mean opacity κ_R (formulated by Svein Rosseland in 1924), which weights monochromatic opacities by the temperature derivative of the Planck blackbody function ∂B_ν(T) / ∂T:

1 / κ_R = [ ∫₀^∞ κ_ν⁻¹ [1 − exp(−hν / k_BT)]⁻¹ (∂B_ν / ∂T) dν ] / [ ∫₀^∞ (∂B_ν / ∂T) dν ]

Because the Planck weighting function ∂B_ν / ∂T peaks at photon energies hν ≈ 3.83 k_BT (Wien's law: ν_peak ∝ T), the frequency factor ν⁻³ in Kramers' monochromatic cross-sections evaluates over the thermal spectrum as ⟨ν⁻³⟩ ∝ T⁻³. Combining the electron velocity/Saha factor (ρ T^−1/2) with the thermal photon frequency factor (T⁻³) produces Kramers' opacity law:

κ_R = κ₀ ρ T^−1/2 T⁻³ = κ₀ ρ T^−7/2

Dependence on Hydrogen (X) and Metal (Z) Mass Fractions

In stellar structure models parameterized by hydrogen mass fraction X, helium mass fraction Y, and heavy-element ("metal", Z > 2) mass fraction Z (such that X + Y + Z = 1), free electrons are supplied predominantly by ionized hydrogen and helium: n_e = (ρ / m_p)(X + Y/2) = (ρ / 2m_p)(1 + X). Standard calibration of the bound-free and free-free Kramers prefactors in SI units (m² kg⁻¹, with ρ in kg m⁻³ and T in K; divide prefactors by 100 for CGS cm² g⁻¹) yields:

• Bound-Free Kramers Opacity (Metals Photoionization):

κ_bf ≈ 4.34 × 10²¹ (g_bf / t) Z (1 + X) ρ T^−7/2 m² kg⁻¹

• Free-Free Kramers Opacity (Hydrogen/Helium Bremsstrahlung):

κ_ff ≈ 3.68 × 10¹⁸ g_ff (1 − Z) (1 + X) ρ T^−7/2 m² kg⁻¹

where t ≈ 1–10 is the guillotine cut-off factor accounting for the progressive stripping of inner K- and L-shell electrons from C, N, O, Ne, and Fe nuclei at high temperatures. Even in solar-metallicity stars (Z ≈ 0.014–0.02), bound-free photoionization of heavy elements is comparable to or exceeds free-free absorption because the Z⁴ nuclear-charge enhancement of bound-free cross-sections compensates for the ~10⁻² mass abundance of metals.

Across the density–temperature plane of stellar interiors, Kramers' opacity operates between low-temperature H⁻ absorption and high-temperature Thomson electron scattering (which sets the Eddington luminosity):

Radiative Core Formation and the Hayashi-to-Henyey Transition

During pre-main-sequence Kelvin-Helmholtz contraction, a young T Tauri star starts as a cool, distended, fully convective polytrope on the Hayashi track because low interior temperatures (T < 10⁶ K) yield immense Kramers opacity, trapping photons and forcing the radiative temperature gradient ∇_rad = (3 κ L P) / (16π a c G m T⁴) above the adiabatic threshold ∇_ad ≈ 0.4. As gravitational contraction shrinks the radius R, the virial theorem dictates that central density scales as ρ_c ∝ M R⁻³ and central temperature scales as T_c ∝ M R⁻¹. Substituting these homologous scalings into Kramers' opacity law reveals how the central opacity evolves during contraction:

κ_c ∝ ρ_c T_c^−7/2 ∝ (R⁻³) (R⁻¹)^−7/2 = R^+1/2

Thus, as the protostar contracts (R decreases), its central Kramers opacity decreases as R^1/2! Once ∇_rad drops below ∇_ad at the center, convection shuts off in the core and a stable radiative core expands outward, deflecting the star onto the nearly horizontal Henyey track on the Hertzsprung–Russell diagram. In stars more massive than ~1.3 M_⊙, subsequent ignition of the highly temperature-sensitive CNO cycle at the Gamow peak (ε_CNO ∝ T¹⁶) concentrates nuclear luminosity at the very center, driving ∇_rad back above ∇_ad and creating a convective core subject to convective overshooting surrounded by a Kramers/electron-scattering radiative envelope.

Thermal Stability and the κ-Mechanism Pulsation Valve

In normal radiative zones governed by Kramers' law, if a spherical shell is adiabatically compressed (δρ / ρ > 0 and δT / T = (γ − 1) δρ / ρ = (2/3) δρ / ρ for an ideal monatomic gas with γ = 5/3), the fractional change in Kramers opacity is:

δκ / κ = (∂lnκ / ∂lnρ) (δρ / ρ) + (∂lnκ / ∂lnT) (δT / T) = [ 1 + (−7/2)(2/3) ] (δρ / ρ) = − (4/3) (δρ / ρ) < 0

Because compression heats the gas and reduces Kramers opacity (δκ / κ = −1.33 δρ / ρ), compressed radiative layers become more transparent and leak excess heat away, strongly damping mechanical oscillations.

However, within the narrow He II → He III partial ionization zone (T ≈ 40,000 K) and the combined H I → H II / He I → He II ionization zone (T ≈ 10,000–15,000 K) of stars crossing the Cepheid instability strip, work done on the gas during compression goes into ionizing atoms rather than raising kinetic temperature. The adiabatic exponent γ − 1 plunges toward 0.1–0.15 while the partial ionization edge causes opacity to rise with temperature (∂lnκ / ∂lnT > 0). Consequently, compression increases the opacity (δκ / κ > 0), damming up outgoing radiation like a thermodynamic piston during maximum compression and releasing the trapped heat during expansion. First identified by Arthur Eddington and quantified by Sergei Zhevakin and John P. Cox, this κ-mechanism (supplemented by the γ-mechanism) drives the coherent radial pulsations of classical Cepheids, RR Lyrae stars, and δ Scuti variables.

1. Kramers, H. A. (1923). "On the theory of X-ray absorption and of the continuous X-ray spectrum." The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, Series 6, 46(275), 836–871. https://doi.org/10.1080/14786442308565244
2. Rosseland, S. (1924). "Note on the absorption of radiation within a star." Monthly Notices of the Royal Astronomical Society, 84(7), 525–528. https://doi.org/10.1093/mnras/84.7.525
3. Iglesias, C. A., & Rogers, F. J. (1996). "Updated Opal Opacities." The Astrophysical Journal, 464, 943–953. https://doi.org/10.1086/177381
4. Kippenhahn, R., Weigert, A., & Weiss, A. (2012). Stellar Structure and Evolution (2nd ed.). Springer. https://doi.org/10.1007/978-3-642-30304-3