Radiation Transport

By the end of this reading, you will be able to:

• Define

  opacity $\kappa$ and optical depth $\tau$ and explain what each measures

• Derive

  the mean free path $\ell = 1/(\kappa\rho)$ from the microscopic relation $\ell = 1/(n\sigma)$

• Use

  a one-zone solar model to estimate the mean free path, optical depth, and diffusion time

• Explain

  why a random walk gives $d \sim \sqrt{N}\,\ell$ rather than $d \sim N\ell$

• Write

  the radiative-diffusion flux and luminosity equations and what controls the energy flow

• Explain

  why convection can take over when radiative transport becomes inefficient

• Compare

  radiation pressure with gas pressure in the Sun

• Derive

  the Eddington luminosity from force balance between radiation and gravity

Guiding question: if photons move at the speed of light, why does energy released in the solar core take roughly $10^5$ years to leak out?

Reading 3 answered the energy-generation question: fusion releases energy in the core. This reading answers a different one: how does that energy get out? The logic chain is the spine of the topic: matter blocks radiation, so stars are opaque; opacity and density set the photon mean free path; a tiny mean free path turns transport into a random walk; random walks are inefficient, so energy diffuses outward slowly; if diffusion becomes too inefficient, convection carries energy by bulk motion; and the same radiation field carries momentum, so radiation pressure can compete with gravity in luminous stars.

Part 1: Why Stars Are Opaque

The travel-time paradox

Start with the wrong model on purpose. If the Sun were transparent, energy from the core would stream straight to the surface in $t_\text{stream} = R_\odot/c = (7.0 \times 10^{10})/(3.0 \times 10^{10}) \approx 2.3~\text{s}$. That is obviously not how a real star works. The missing physics is opacity.

Opacity: how effectively matter blocks radiation

The opacity $\kappa$ measures how strongly matter absorbs or scatters radiation per unit mass, with units $\text{cm}^2\,\text{g}^{-1}$. A larger $\kappa$ means matter intercepts radiation more effectively; in stellar structure $\kappa$ is usually a mean opacity packaging many microscopic processes into one quantity.

Opacity

A measure of how strongly matter absorbs or scatters radiation per unit mass, $\kappa$ (units $\text{cm}^2\,\text{g}^{-1}$). Larger opacity means a shorter photon mean free path and more resistance to radiative energy flow.

For a fully ionized hydrogen-rich plasma, the cleanest baseline is electron scattering: $\kappa_\text{es} \approx 0.2(1 + X) \approx 0.34~\text{cm}^2\,\text{g}^{-1}$ with $X \approx 0.70$ the hydrogen mass fraction — the standard one-zone baseline for the solar interior.

Optical depth: how many mean free paths thick is the star?

Combine opacity with density. The optical-depth increment is $d\tau = \kappa\rho\,dr$; integrating across a region of size $R$ gives the optical depth $\tau \sim \kappa\rho R$. Using the mean free path $\ell = 1/(\kappa\rho)$, this is $\tau \sim R/\ell$ — literally how many mean free paths fit across the star.

Optical depth

The dimensionless thickness of a medium in mean free paths, $\tau \sim \kappa\rho R \sim R/\ell$. $\tau \ll 1$ is optically thin (photons stream freely); $\tau \gg 1$ is optically thick (photons interact many times).

Quick check

Suppose a star becomes denser while its opacity stays the same. (1) Does the mean free path get longer or shorter? (2) Does the optical depth get larger or smaller? (3) Should energy escape more easily or less easily?

Part 2: The Mean Free Path — How Far a Photon Travels

Start microscopically. If a photon moves through targets of number density $n$, each with cross section $\sigma$, the mean free path is $\ell = 1/(n\sigma)$. Since opacity is cross-section per unit mass, $\kappa = n\sigma/\rho$, so $n\sigma = \kappa\rho$ and:

Higher density means more targets per volume; higher opacity means more effective targets per gram. Both shorten the distance a photon travels before interacting.

Mean free path

The average distance a photon travels between interactions, $\ell = 1/(\kappa\rho)$ — about $0.02~\text{cm}$ in the solar core. When $\ell \ll R$, radiation cannot stream out; it diffuses.

Worked Example 1Mean Free Path in the Solar Core

Problem

Use one-zone solar-core values $\rho_c \approx 150~\text{g cm}^{-3}$ and $\kappa_\text{es} \approx 0.34~\text{cm}^2\,\text{g}^{-1}$ to find the mean free path, and compare it to the solar radius.

Step·Compute the mean free path

$\kappa\rho_c = (0.34)(150) = 51~\text{cm}^{-1} \Rightarrow \ell = \frac{1}{\kappa\rho_c} \approx 1.96 \times 10^{-2}~\text{cm} \approx 0.20~\text{mm}$

Step·Count mean free paths across the Sun

$\frac{R_\odot}{\ell} \approx \frac{7.0 \times 10^{10}}{2.0 \times 10^{-2}} \approx 3.5 \times 10^{12} = \tau$

Dimensional check

$[\kappa\rho] = (\text{cm}^2\text{g}^{-1})(\text{g cm}^{-3}) = \text{cm}^{-1}$, so $[\ell] = 1/[\kappa\rho] = \text{cm}$ ✓.

Result

A photon travels only ~0.2 mm between interactions, and the Sun is $\sim 3.5 \times 10^{12}$ mean free paths thick — not just opaque, but enormously optically thick. The crucial result is $\ell \ll R_\odot$: the diffusion regime.

Quick check

Suppose $\kappa$ doubles while $\rho$ stays fixed. (1) What happens to $\ell$? (2) To $\tau$ across the same star? (3) Does diffusion become faster or slower?

Part 3: The Random Walk — Why It Takes Hundreds of Thousands of Years

A photon does not travel radially outward in one path. After each interaction, the next step goes in a new direction — a random walk. For step size $\ell$ and $N$ steps, the net displacement is:

The total path length is $d_\text{total} \sim N\ell$, but the net outward displacement is only $d \sim \sqrt{N}\,\ell$. So to cross a distance $R$ requires $N \sim (R/\ell)^2$ steps.

Random walk

A path of many steps each taken in an independent random direction. The net displacement grows only as $\sqrt{N}\,\ell$ (not $N\ell$), which is why crossing a distance $R$ needs $N \sim (R/\ell)^2$ steps and makes radiative transport slow.

Worked Example 2Photon Diffusion Time in the Sun

Problem

Use $R_\odot = 7.0 \times 10^{10}~\text{cm}$ and $\ell \approx 2.0 \times 10^{-2}~\text{cm}$ to estimate the number of scatterings and the diffusion time.

Step·Number of steps to cross the Sun

$N \sim \left(\frac{R_\odot}{\ell}\right)^2 \sim (3.5 \times 10^{12})^2 \approx 1.2 \times 10^{25}$

Step·Diffusion time (each step takes ell/c)

$t_\text{diff} \sim N\frac{\ell}{c} \sim \frac{R_\odot^2}{\ell c} = \frac{(7.0 \times 10^{10})^2}{(2.0 \times 10^{-2})(3.0 \times 10^{10})}~\text{s} \approx 2.6 \times 10^5~\text{yr}$

Dimensional check

$\left[R^2/(\ell c)\right] = \text{cm}^2 / (\text{cm} \cdot \text{cm s}^{-1}) = \text{s}$ ✓.

Result

About $10^{25}$ scatterings and $t_\text{diff} \sim 2.6 \times 10^5~\text{yr}$ — the right order of magnitude. More realistic stratified models stay in the same $10^5$-year range.

A useful alternate form: since $\tau \sim R/\ell$, the diffusion time is $t_\text{diff} \sim \tau\,(R/c)$. Here $R/c$ is the straight-line crossing time and $\tau$ is how many mean-free-path layers thick the star is — an optically thick star multiplies the crossing time by an enormous factor.

Problem

Use $t_\text{diff} \sim R^2/(\ell c)$. (1) If $\ell$ decreases by a factor of 10, what happens to $t_\text{diff}$? (2) If $R$ doubles while $\ell$ stays fixed?

Part 4: Radiative Diffusion — Energy Flow as a Leaking Fog

The random walk tells us why transport is slow; stellar structure needs a law for how much luminosity a temperature gradient can carry. A thermal radiation field has energy density $u_\text{rad} = aT^4$ (with $a = 7.56 \times 10^{-15}~\text{erg cm}^{-3}\text{K}^{-4}$). If the interior is hotter than the exterior, $u_\text{rad}$ decreases outward, and that gradient drives a net outward flux:

A steeper temperature gradient drives a larger flux; higher opacity or density suppress it; the minus sign means energy flows from hot to cool. Opacity acts like a resistance to energy flow — this is radiative diffusion.

Radiative diffusion

Energy transport by photons random-walking down a temperature gradient through an optically thick medium. The diffusive flux is $F_\text{rad} = -(c/3\kappa\rho)\,d(aT^4)/dr$ — throttled by opacity, driven by the gradient.

Luminosity is flux times area, $L(r) = 4\pi r^2 F_\text{rad}$:

This treatment assumes the medium is optically thick, in local thermodynamic equilibrium, radiation-dominated, and close to isotropic.

Reading the Math: the diffusion luminosity

The flux law carries a derivative, $d(aT^4)/dr$, so we run the signature move again — now on the temperature profile — and read off how a star’s luminosity scales with its mass and radius.

① Approximate the derivative. Temperature falls from $T_c$ at the center to essentially nothing at the surface over a distance $\sim R$:

$\frac{d(aT^4)}{dr} \;\approx\; \frac{aT_\text{surf}^4 - aT_c^4}{R - 0} \;\approx\; -\frac{aT_c^4}{R}.$

The coefficient out front, $c/(3\kappa\rho)$, is the radiative diffusion coefficient $D = c\ell/3$ — Part 2’s tiny mean free path $\ell$ in one symbol. A small $\ell$ makes $D$ small, which is the random walk of Part 3.

② Extract the scaling. Substitute into $F = -(c/3\kappa\rho)\,d(aT^4)/dr$, form the luminosity $L = 4\pi R^2 F$, and replace $\rho \sim M/R^3$:

$L \;\sim\; 4\pi R^2 \cdot \frac{c}{3\kappa\rho}\cdot\frac{aT_c^4}{R} \;\sim\; \frac{4\pi a c}{3\kappa}\,\frac{T_c^4 R^4}{M}.$

Read the exponents — but note $T_c$ is not free. Reading 5 substitutes the hydrostatic core temperature $T_c \sim \mu G M m_p/(k_B R)$; the four factors of $R$ cancel and one power of $M$ survives, leaving the headline mass–luminosity law $L \propto M^3$. Hydrostatic equilibrium (Reading 2) built the first half of that result; this diffusion law is the second.

③ Name the assumption. The audit gains a transport row:

Same deal as Reading 2: right exponents, approximate coefficient. The scaling $L \propto M^3$ is robust; the prefactor needs the true opacity and structure.

At fixed luminosity, high opacity makes radiative transport less efficient, so the star needs a steeper temperature gradient to carry the same energy. If the required gradient becomes too steep, radiation is no longer the preferred mechanism and convection can take over.

Multiple choice

At some radius, $L(r)$ is fixed but $\kappa$ increases by a factor of 10. To carry the same luminosity, what must the temperature gradient $|dT/dr|$ do?

Steepen — radiative transport is harderShallow out — radiative transport is easierStay the same

Part 5: When Radiation Isn’t Enough — Convection

Radiative diffusion has a built-in limit: photons move energy by many tiny random steps, high opacity makes those steps smaller, and energy flow becomes inefficient. To carry a fixed luminosity through a very opaque region, the star must steepen the temperature gradient. What if it becomes too steep?

Unlike photons, the stellar plasma is a fluid. A blob slightly hotter than its surroundings is less dense, so buoyancy pushes it upward, carrying energy; cooler, denser gas sinks. This circulating flow is convection.

Convection

Energy transport by the bulk motion of fluid, driven by buoyancy: hot material rises, cool material sinks, and energy is carried by moving matter rather than diffusing radiation. It takes over where the radiative temperature gradient would be too steep to be stable.

Compared mechanism to mechanism — radiation makes many tiny random-walk steps, convection moves energy by large-scale motion — convection transports energy much more efficiently where diffusion struggles. Conceptually: if radiation can carry the luminosity the star stays radiative; if it struggles, the gradient steepens; if the gradient becomes too steep, the fluid becomes unstable and convection turns on.

The Schwarzschild criterion

Make the buoyancy argument quantitative. Displace a blob of gas upward by $dr$ and let it stay in pressure balance with its new surroundings while expanding adiabatically — it rises too fast to exchange heat. The blob keeps going — the layer is unstable to convection — only if it remains less dense than the gas around it after the nudge. At fixed pressure, lower density means higher temperature, so the instability condition is that the surroundings cool faster with height than the adiabatic blob does:

$\left|\frac{dT}{dr}\right|_\text{rad} > \left|\frac{dT}{dr}\right|_\text{ad} \qquad\Longleftrightarrow\qquad \nabla_\text{rad} > \nabla_\text{ad},$

written compactly with the logarithmic gradient $\nabla \equiv d\ln T/d\ln P$. For an ideal monatomic gas $\nabla_\text{ad} = 1 - 1/\gamma \approx 0.4$ is essentially fixed, while the radiative gradient $\nabla_\text{rad}$ is whatever slope diffusion must adopt to carry the luminosity. From the diffusion law, $\nabla_\text{rad}$ grows with $\kappa$ and with the local luminosity-to-mass ratio $L/M$. So convection switches on exactly where the required radiative gradient is forced past the adiabatic value — in high-opacity zones (cool stellar envelopes) and where energy generation is fiercely concentrated (the CNO-burning cores of massive stars, Reading 3).

Name the assumption. This criterion assumes the displaced blob exchanges no heat (perfectly adiabatic), stays in instantaneous pressure equilibrium, and meets no composition gradient — refinements that real stellar-structure codes restore one at a time.

Quick check

If opacity increases in some region of a star, does radiative transport become easier or harder? What must happen to the temperature gradient? Would this make convection more or less likely? Explain in words.

Part 6: Radiation Pressure

Photons carry not only energy but momentum: $p = E/c$. So when radiation is absorbed or scattered, it exerts a force on matter. In an isotropic thermal radiation field, the radiation pressure is:

This is why radiation transport is also a force problem.

Radiation pressure

The pressure exerted by a photon field through its momentum, $P_\text{rad} = aT^4/3$ for an isotropic thermal field. It grows as $T^4$ — negligible in the Sun, but dominant in very massive, hot stars.

Worked Example 3Radiation Pressure vs. Gas Pressure in the Sun

Problem

Compare $P_\text{rad} = aT^4/3$ with the ideal-gas pressure $P_\text{gas} \approx \rho k_B T/(\mu m_H)$ at the solar core ($T_c \approx 1.5 \times 10^7~\text{K}$, $\rho_c \approx 150~\text{g cm}^{-3}$, mean molecular weight $\mu \approx 0.61$).

Step·Radiation pressure

$P_\text{rad} = \frac{aT_c^4}{3} \approx 1.3 \times 10^{14}~\text{dyn cm}^{-2}$

Step·Gas pressure

$P_\text{gas} \approx \frac{(150)(1.38 \times 10^{-16})(1.5 \times 10^7)}{(0.61)(1.67 \times 10^{-24})} \approx 3.0 \times 10^{17}~\text{dyn cm}^{-2}$

Dimensional check

Both reduce to $\text{erg cm}^{-3} = \text{dyn cm}^{-2}$ (a pressure) ✓.

Result

$P_\text{rad}/P_\text{gas} \approx 4 \times 10^{-4}$ — radiation pressure is tiny in the Sun. But $P_\text{rad} \propto T^4$ while $P_\text{gas} \propto T$ (at fixed $\rho$), so radiation pressure grows far faster with temperature — which is why it dominates in very massive stars.

Deriving the Eddington luminosity

How luminous can a star become before radiation pushes outward as strongly as gravity pulls inward? Assume spherical symmetry, steady state, and electron-scattering opacity. The outward radiative acceleration on a shell is $g_\text{rad} = \kappa F/c = \kappa L/(4\pi r^2 c)$; the inward gravitational acceleration is $g = GM/r^2$. Setting them equal (the $r^2$ cancels) gives the Eddington luminosity:

A force-balance ceiling: below $L_\text{Edd}$ gravity wins; near it, radiation reshapes the envelope and drives strong winds; above it, hydrostatic balance becomes very hard to maintain.

Eddington luminosity

The luminosity at which radiation’s outward force balances gravity, $L_\text{Edd} = 4\pi GMc/\kappa$ (about $3.8 \times 10^4\,L_\odot$ for the Sun). A structural ceiling, not a hard cutoff — stars near it drive winds and instabilities.

Worked Example 4The Sun's Eddington Luminosity

Problem

Compute $L_\text{Edd}$ for the Sun with $\kappa \approx 0.34~\text{cm}^2\text{g}^{-1}$, $M_\odot = 2.0 \times 10^{33}~\text{g}$.

Step·Plug into 4πGMc/κ

$L_\text{Edd} = \frac{4\pi (6.674 \times 10^{-8})(2.0 \times 10^{33})(3.0 \times 10^{10})}{0.34} \approx 1.5 \times 10^{38}~\text{erg s}^{-1}$

Dimensional check

$\left[\dfrac{GMc}{\kappa}\right] = \dfrac{(\text{cm}^3\text{g}^{-1}\text{s}^{-2})(\text{g})(\text{cm s}^{-1})}{\text{cm}^2\text{g}^{-1}} = \text{g cm}^2\text{s}^{-3} = \text{erg s}^{-1}$ ✓.

Result

$L_\text{Edd} \approx 1.5 \times 10^{38}~\text{erg s}^{-1} \approx 3.8 \times 10^4\,L_\odot$. The Sun is far below this, so radiation pressure is dynamically negligible for solar structure — but massive stars come much closer because their luminosities rise faster than linearly with mass.

Problem

$L_\text{Edd} \propto M/\kappa$. (1) At fixed opacity, if stellar mass doubles, what happens to $L_\text{Edd}$? (2) Why do the most massive stars come closest to this limit even though $L_\text{Edd}$ rises with mass?

Reference Tables

Radiation Transport at a Glance

Symbol Legend

Summary: Energy’s Tortuous Journey

1. Stars are opaque — the photon mean free path in the solar core is $\sim 0.2~\text{mm}$, about a hundred billion times smaller than the solar radius.
2. The random walk explains the diffusion time: a one-zone estimate gives $t_\text{diff} \sim 2.6 \times 10^5~\text{yr}$, and realistic models stay in the same $10^5$-year range.
3. Radiative diffusion transports energy down the temperature gradient; the random walk explains why it is slow, the diffusion equation how much luminosity a gradient can carry.
4. Convection takes over when radiative diffusion becomes too inefficient — energy is then carried by rising and sinking fluid.
5. Radiation pressure ($P_\text{rad} = aT^4/3$) is negligible in the Sun but grows as $T^4$; in very luminous stars it competes with gravity and sets the Eddington luminosity.

Retrieval — equation: mean-free-path▸

Glossary

Convection

Energy transport by the bulk motion of fluid, driven by buoyancy: hot material rises, cool material sinks, and energy is carried by moving matter rather than diffusing radiation. It takes over where the radiative temperature gradient would be too steep to be stable.

Eddington luminosity

The luminosity at which radiation’s outward force balances gravity, $L_\text{Edd} = 4\pi GMc/\kappa$ (about $3.8 \times 10^4\,L_\odot$ for the Sun). A structural ceiling, not a hard cutoff — stars near it drive winds and instabilities.

Mean free path

The average distance a photon travels between interactions, $\ell = 1/(\kappa\rho)$ — about $0.02~\text{cm}$ in the solar core. When $\ell \ll R$, radiation cannot stream out; it diffuses.

Opacity

A measure of how strongly matter absorbs or scatters radiation per unit mass, $\kappa$ (units $\text{cm}^2\,\text{g}^{-1}$). Larger opacity means a shorter photon mean free path and more resistance to radiative energy flow.

Optical depth

The dimensionless thickness of a medium in mean free paths, $\tau \sim \kappa\rho R \sim R/\ell$. $\tau \ll 1$ is optically thin (photons stream freely); $\tau \gg 1$ is optically thick (photons interact many times).

Radiation pressure

The pressure exerted by a photon field through its momentum, $P_\text{rad} = aT^4/3$ for an isotropic thermal field. It grows as $T^4$ — negligible in the Sun, but dominant in very massive, hot stars.

Radiative diffusion

Energy transport by photons random-walking down a temperature gradient through an optically thick medium. The diffusive flux is $F_\text{rad} = -(c/3\kappa\rho)\,d(aT^4)/dr$ — throttled by opacity, driven by the gradient.

Random walk

A path of many steps each taken in an independent random direction. The net displacement grows only as $\sqrt{N}\,\ell$ (not $N\ell$), which is why crossing a distance $R$ needs $N \sim (R/\ell)^2$ steps and makes radiative transport slow.