Energy Transport by Conduction

Energy transport by conduction is heat carried by the motion and collisions of massive particles — in a star, the free electrons — rather than by photons. In ordinary stellar plasma this channel is almost irrelevant: photons have far longer mean free paths and travel at $c$ , so radiative diffusion carries essentially all the flux and conduction is a rounding error. The situation inverts once matter is compressed to electron degeneracy. When the Pauli principle forces the electrons into a filled Fermi sea, only those within $\sim k_B T$ of the Fermi surface can scatter; the rest glide collisionlessly, their mean free paths grow enormous, and the gas conducts heat almost like a block of metal. Conduction then becomes the most efficient way to move energy — with only a whisper of a temperature gradient. The physical archetype is the white dwarf core, together with the degenerate helium and carbon–oxygen cores of evolved stars. There the conductive flux is so large that the interior is very nearly isothermal — a thermal short circuit — while all the temperature structure is squeezed into a thin, radiatively opaque envelope that throttles the cooling. White dwarf models therefore treat the interior as a uniform-temperature reservoir in an insulating skin, the setting for white dwarf cooling and the runaway ignition of the helium flash. Remarkably, none of this needs a new structural equation: conduction is folded in by writing its flux as a diffusion law with an effective conductive opacity, and because the two channels act in parallel their opacities combine harmonically.

Key takeaways

Conduction transports heat by particle collisions; in stars the carrier is the degenerate electron gas, whose long mean free paths make it a superb conductor in dense cores.
Conduction and radiation act in parallel under the same gradient, so their opacities add inversely: $1/\kappa = 1/\kappa_{\rm rad} + 1/\kappa_{\rm cond}$ , and the combined $\kappa$ slots into the radiative gradient.
In non-degenerate gas conduction is negligible because photons carry heat far more efficiently; it dominates only where degeneracy suppresses electron scattering.
In white dwarf and degenerate red-giant cores conduction is so efficient the interior is nearly isothermal, decoupling core temperature from the steep envelope.
Degenerate conductive opacity scales as $\kappa_{\rm cond}\propto T^2/\rho^2$ — falling as density rises, opposite to radiative opacity — so conduction inevitably wins at high $\rho$ .
Key idea: conduction and radiation are parallel diffusive channels combined as $1/\kappa = 1/\kappa_{\rm rad} + 1/\kappa_{\rm cond}$ ; because degenerate electron conductive opacity falls as $\rho^{-2}$ , conduction dominates dense cores and drives them to near-isothermality.

The core insight

Conduction and radiation are parallel diffusive channels; combining them as $1/\kappa = 1/\kappa_{\rm rad} + 1/\kappa_{\rm cond}$ lets a single effective opacity carry the flux, and because degenerate electron conductive opacity falls as $\rho^{-2}$ , conduction dominates and nearly isothermalizes dense degenerate cores.

Where this lives today

Conductive opacities are working inputs in every modern stellar-evolution code. In white dwarf cosmochronology, the age of the Galactic disk and of globular clusters is read off the faint end of the white dwarf luminosity function, and the inferred ages are only as accurate as the conductive opacities that set the core's cooling rate — the 2007 Cassisi tables shifted derived ages by hundreds of Myr. The same physics governs the helium flash in low-mass red giants, where efficient conduction keeps the degenerate core isothermal so ignition runs away nearly simultaneously. In cooling neutron stars, the relation between the observed surface temperature and the hot interior — what Chandra and XMM-Newton constrain — is fixed by electron conduction through the degenerate envelope, the same physics at nuclear density. And warm-dense-matter experiments at the National Ignition Facility and X-ray free-electron lasers now measure electron transport in compressed degenerate plasmas directly, testing the models stellar physics has relied on for decades.

The physical picture

Picture two independent pipes carrying heat between the same pair of temperatures: photons through one, electrons through the other. If either pipe is wide — low resistance — heat flows easily. Radiative opacity $\kappa_{\rm rad}$ is the resistance of the photon pipe, conductive opacity $\kappa_{\rm cond}$ that of the electron pipe. Because both feel the same gradient and their fluxes add, the resistances combine like parallel resistors: the inverse opacities sum, and the wider pipe dominates. In an ordinary star the photon pipe is far wider; compress the gas into degeneracy and the electron pipe opens until it swamps the photon pipe. The microscopic reason is the Fermi sea. In a degenerate gas the electrons fill every momentum state up to the Fermi momentum, and the exclusion principle blocks almost all scattering: an electron can only be deflected into an empty state, and empty states exist only within about $k_B T$ of the Fermi surface. Most electrons therefore travel without scattering, their mean free paths grow very long, and the conductivity soars — the same reason a cold metal conducts so well. A white dwarf core is, in this respect, a ball of degenerate metal the size of the Earth. The two mechanisms respond oppositely to compression. Squeeze the gas and there are more electrons at a higher Fermi energy with longer relative mean free paths, so conduction gets stronger; radiation is throttled by the extra matter it must diffuse through. Deep in a compact star the channels swap dominance at a crossover density, below which radiation wins and above which conduction wins.

Figure

White dwarf temperature profile: conduction makes the degenerate core nearly isothermal with κ ≈ κ_cond and ΔT/T ≈ 2 percent, pushing the whole temperature drop into the thin envelope where radiation sets the cooling rate — 1/κ = 1/κ_rad + 1/κ_cond

Because conduction needs only a tiny gradient to carry a large luminosity, the degenerate interior spreads its heat until it sits at one temperature — the flat core and steep envelope above. The conductive opacity is not constant: it plunges as density climbs, and where it undercuts the radiative opacity marks the edge of the isothermal region, quantified next.

Figure

Opacity crossover at T = 10⁷ K for a carbon–oxygen core: radiative opacity stays near κ_rad ≈ 0.2 cm²/g while conductive opacity falls as κ_cond ∝ T²/ρ², so above ρ ≈ 1.6×10⁴ g/cm³ conduction short-circuits radiation and the effective 1/κ = 1/κ_rad + 1/κ_cond collapses toward κ_cond

The core equations

Conduction begins with the Fourier law: a heat flux proportional to the temperature gradient, set by the conductivity of the carrier — here the degenerate electron gas.

F_{\rm cond} = -K_{\rm cond}\,\frac{\mathrm{d}T}{\mathrm{d}r} \qquad (1)

where: $F_{\rm cond}$ — conductive energy flux (per area per time); $K_{\rm cond}$ — thermal conductivity of the electron gas; $\mathrm{d}T/\mathrm{d}r$ — radial temperature gradient. The minus sign sends heat from hot to cold; a larger conductivity carries the same flux with a smaller gradient. The radiative flux in an optically thick interior takes the diffusion form of the radiative-transfer article, photons random-walking out at a rate set by the Rosseland mean opacity. It is quoted here to match forms against conduction.

F_{\rm rad} = -\frac{4ac\,T^3}{3\kappa_{\rm rad}\rho}\,\frac{\mathrm{d}T}{\mathrm{d}r} \qquad (2)

where: $a$ — radiation constant, $a=7.566\times10^{-16}\ \mathrm{J\,m^{-3}\,K^{-4}}$ ; $c$ — speed of light; $\kappa_{\rm rad}$ — radiative (Rosseland mean) opacity; $\rho$ — density; $T$ — temperature. Both fluxes are linear in $\mathrm{d}T/\mathrm{d}r$ — the whole trick — so they share a shape. The masterstroke is to define an equivalent conductive opacity making Fourier's law identical to radiative diffusion; matching coefficients (1) and (2) fixes it.

\kappa_{\rm cond}\equiv\frac{4ac\,T^3}{3\rho\,K_{\rm cond}} \qquad (3)

where: $\kappa_{\rm cond}$ — equivalent conductive opacity; $K_{\rm cond}$ — electron thermal conductivity. A large conductivity means a small opacity — an efficient channel — packaging the electron-gas physics into one opacity-like number. With both fluxes in the same form under the same gradient they superpose; the common coefficient factors out, leaving the inverse opacities added.

F_{\rm tot}=F_{\rm rad}+F_{\rm cond}=-\frac{4ac\,T^3}{3\rho}\left(\frac{1}{\kappa_{\rm rad}}+\frac{1}{\kappa_{\rm cond}}\right)\frac{\mathrm{d}T}{\mathrm{d}r} \qquad (4)

where: $F_{\rm tot}$ — total diffusive energy flux. The parenthesis sums two conductances; the effective opacity is its reciprocal — the harmonic combination at the heart of the topic.

\frac{1}{\kappa}=\frac{1}{\kappa_{\rm rad}}+\frac{1}{\kappa_{\rm cond}} \qquad (5)

where: $\kappa$ — effective opacity inserted into the standard gradient. The smaller opacity dominates $\kappa$ , as the smaller resistor dominates a parallel pair — the more efficient channel wins. Because the total flux keeps the pure-radiative form with $\kappa_{\rm rad}\to\kappa$ , no new equation is needed: the standard gradient equation absorbs conduction automatically.

\frac{\mathrm{d}T}{\mathrm{d}r}=-\frac{3\kappa\rho}{4ac\,T^3}\,\frac{L_r}{4\pi r^2} \qquad (6)

where: $L_r$ — luminosity crossing the sphere of radius $r$ ; $r$ — radial coordinate. The gradient carrying a fixed $L_r$ is proportional to $\kappa$ : shrink the opacity and the gradient collapses. In the degenerate, non-relativistic regime the conductive opacity has a telling scaling — it falls as density rises, opposite to a Kramers radiative opacity.

\kappa_{\rm cond}\propto \frac{T^2}{\rho^2}\qquad\text{(degenerate, non-relativistic)} \qquad (7)

where: composition enters through the ionic charge $Z$ and $\mu_e$ (electron–ion scattering). This power law, derived below, is why conduction always overtakes radiation as matter is compressed. Whether the gas is degenerate is set by comparing the thermal energy $k_B T$ with the Fermi energy $E_F$ ; conduction dominates only well inside the degenerate regime.

k_B T \lesssim E_F,\qquad E_F=\frac{\hbar^2}{2m_e}\left(3\pi^2 n_e\right)^{2/3} \qquad (8)

where: $E_F$ — electron Fermi energy; $n_e$ — electron number density; $m_e$ — electron mass; $\hbar$ — reduced Planck constant; $k_B$ — Boltzmann constant. The parameter $E_F/k_B T$ measures how deep into the Fermi sea the gas sits; in a white dwarf core it reaches the hundreds.

Derivation AThe refractive-index trick — casting conduction as a diffusion equation

The goal is to show why conduction requires no new structural equation: because radiative and conductive fluxes are both linear in the same temperature gradient, an equivalent conductive opacity lets them add as parallel channels, and the sum of their inverse opacities defines one effective opacity that carries the flux through the ordinary gradient equation with $\kappa_{\rm rad}\to\kappa$ . We match forms, superpose, and read off the combination.

Step 1. Write the radiative flux in its diffusion form, $F_{\rm rad}=-\dfrac{4ac\,T^3}{3\kappa_{\rm rad}\rho}\dfrac{\mathrm{d}T}{\mathrm{d}r}$ , valid in the optically thick interior. Why: this is the reference shape into which we force conduction, so both channels compare coefficient by coefficient.

Step 2. Give Fourier's law $F_{\rm cond}=-K_{\rm cond}\,\mathrm{d}T/\mathrm{d}r$ the same form by defining $\kappa_{\rm cond}\equiv 4ac\,T^3/(3\rho K_{\rm cond})$ , so $F_{\rm cond}=-\dfrac{4ac\,T^3}{3\kappa_{\rm cond}\rho}\dfrac{\mathrm{d}T}{\mathrm{d}r}$ . Why: a pure re-labelling — no physics added — converting a conductivity into an opacity, the currency the structure equations speak.

Step 3. Add the two fluxes; both respond to the identical gradient, so they superpose: $F_{\rm tot}=-\dfrac{4ac\,T^3}{3\rho}\left(\kappa_{\rm rad}^{-1}+\kappa_{\rm cond}^{-1}\right)\dfrac{\mathrm{d}T}{\mathrm{d}r}$ . Why: heat leaves through both pipes at once, so it is the fluxes that add, not the opacities.

Step 4. Define the effective opacity $\kappa^{-1}\equiv\kappa_{\rm rad}^{-1}+\kappa_{\rm cond}^{-1}$ , the harmonic combination. Why: one coefficient reproduces the whole bracket, so a single opacity stands in for both mechanisms — parallel conductances, not series resistors.

Step 5. Set the flux to the local value $F_{\rm tot}=L_r/(4\pi r^2)$ and rearrange to $\dfrac{\mathrm{d}T}{\mathrm{d}r}=-\dfrac{3\kappa\rho}{4ac\,T^3}\dfrac{L_r}{4\pi r^2}$ , identical to the pure radiative equation under $\kappa_{\rm rad}\to\kappa$ . Why: the transport equation of stellar structure is untouched in form, so a single combined opacity table is all a model needs to include conduction. *Two limiting checks: when $\kappa_{\rm cond}\to\infty$ the harmonic combination gives $\kappa\to\kappa_{\rm rad}$ , recovering the pure radiative gradient; when $\kappa_{\rm cond}\to 0$ it gives $\mathrm{d}T/\mathrm{d}r\to 0$ , the isothermal limit. In both cases the combined opacity locks to the smaller of the two values.*

Go deeper · Derivation B · 1 minWhere the T-squared over rho-squared scaling comes from

Derivation BWhere the T-squared over rho-squared scaling comes from

The goal is to derive the density and temperature dependence of the degenerate conductive opacity, $\kappa_{\rm cond}\propto T^2/\rho^2$ , from the kinetic theory of an electron gas whose mean free path is set by electron–ion Coulomb scattering — and so explain the central fact of the topic, that conduction gets better with compression while radiation gets worse. We build the conductivity from a Wiedemann–Franz estimate, insert the degenerate kinematics, and feed it into $\kappa_{\rm cond}$ .

Step 1. Estimate the conductivity kinetically as $K_{\rm cond}\sim \dfrac{\pi^2}{3}\dfrac{n_e k_B^2 T\,\tau}{m_e}$ , with $\tau=\ell/v_F$ , $\ell=1/(n_i\sigma)$ , and $v_F=p_F/m_e$ . Why: only electrons within $k_B T$ of the Fermi surface carry heat, so the conductivity takes the metallic Wiedemann–Franz form, not the classical $\sqrt{T}$ law.

Step 2. Insert the degenerate kinematics: $p_F\propto n_e^{1/3}$ , so $E_F\propto n_e^{2/3}$ ; electron–ion Coulomb scattering gives $\sigma\propto Z^2/E_F^2\propto Z^2 n_e^{-4/3}$ (up to a Coulomb log), and $n_i=n_e/Z$ . Why: the exclusion principle fixes scattering at the Fermi energy, so $v_F$ and $\sigma$ are frozen functions of density, not temperature.

Step 3. Combine the denominator: $m_e v_F\,n_i\,\sigma\propto n_e^{1/3}\cdot(n_e/Z)\cdot(Z^2 n_e^{-4/3})=Z\,n_e^{0}$ , so $K_{\rm cond}\propto n_e k_B^2 T/Z\propto \rho\,T/(\mu_e Z)$ . Why: the density powers cancel to leave $K_{\rm cond}$ rising linearly with $\rho$ — squeezing the gas improves conduction.

Step 4. Feed this into $\kappa_{\rm cond}=4ac\,T^3/(3\rho K_{\rm cond})\propto T^3/(\rho\cdot \rho T)=T^2/\rho^2$ . Why: the $T^3$ of radiative diffusion over a conductivity growing as $\rho T$ leaves the signature $T^2/\rho^2$ — opposite in sign to Kramers ( $\kappa_{\rm rad}\propto\rho\,T^{-3.5}$ ), which guarantees a crossover. *The $\rho^{-2}$ exponent is consistent with the crossover density shown in the figure: at $T=10^7$ K and $\kappa_{\rm rad}\approx0.2\ \mathrm{cm^2\,g^{-1}}$ the two opacities meet at $\rho\approx1.6\times10^4\ \mathrm{g\,cm^{-3}}$ . Because $\kappa_{\rm cond}$ falls with density while $\kappa_{\rm rad}$ rises, conduction inevitably wins above that crossover. Detailed tables add relativistic and Coulomb-log corrections but preserve the leading powers.*

Go deeper · Derivation C · 1 minWhy a degenerate core goes nearly isothermal

Derivation CWhy a degenerate core goes nearly isothermal

The goal is to show quantitatively that a degenerate core cannot sustain more than a tiny temperature gradient, so it is modelled as an isothermal reservoir with the whole temperature drop pushed into the thin envelope. We solve the transport equation for the gradient using the small conductive opacity, evaluate it for realistic white dwarf numbers, and compare against what radiation alone would demand.

Step 1. In the degenerate interior $\kappa_{\rm cond}\ll\kappa_{\rm rad}$ , so the harmonic combination locks onto the smaller value, $\kappa\approx\kappa_{\rm cond}$ . Why: the effective opacity is dominated by the more efficient channel, so conduction alone sets transport.

Step 2. The gradient carrying $L_r$ is $\left|\dfrac{\mathrm{d}T}{\mathrm{d}r}\right|=\dfrac{3\kappa_{\rm cond}\rho}{4ac\,T^3}\dfrac{L_r}{4\pi r^2}$ , proportional to $\kappa_{\rm cond}$ . Why: a low-resistance channel moves a fixed flux with almost no push, so the gradient shrinks in lockstep with the opacity.

Step 3. Evaluate for a carbon–oxygen core with $\rho=10^9\ \mathrm{kg\,m^{-3}}$ , $T=10^7$ K, $\kappa_{\rm cond}\approx5\times10^{-6}\ \mathrm{m^2\,kg^{-1}}$ , $L_r\approx10^{-3}L_\odot$ at $r\approx3\times10^6$ m: $|\mathrm{d}T/\mathrm{d}r|\approx0.055\ \mathrm{K\,m^{-1}}$ (full arithmetic in Worked example 3). Why: over the whole core radius this integrates to $\Delta T\approx1.7\times10^5$ K, only about 2% of the $10^7$ K interior — the core really is nearly isothermal.

Step 4. Compression drives it flatter still, since $\kappa_{\rm cond}\propto T^2/\rho^2$ falls toward the centre where $\rho$ is largest. The temperature drop is banished to the thin envelope, whose radiative opacity insulates the reservoir and fixes the cooling rate. Why: this two-zone structure — isothermal core, insulating radiative skin — is the basis of Mestel's cooling theory and of the near-simultaneous ignition that becomes the helium flash. *For comparison, radiation alone — with $\kappa$ 4000 times larger — would require a gradient of $\approx220\ \mathrm{K\,m^{-1}}$ and an integrated temperature drop $\Delta T\sim7\times10^8$ K, larger than the core temperature itself and therefore impossible. Conduction is what makes a nearly isothermal core physically achievable.*

Limiting cases

How conduction competes with radiation across the regimes a star's matter passes through:

Case	Result	Meaning
Non-degenerate gas	$\kappa_{\rm cond}\gg\kappa_{\rm rad}\ \Rightarrow\ \kappa\approx\kappa_{\rm rad}$	Ordinary plasma (the Sun's interior): photons win, conduction is a rounding error, and the effective opacity is the radiative one.
Strongly degenerate core	$\kappa_{\rm cond}\ll\kappa_{\rm rad}\ \Rightarrow\ \kappa\approx\kappa_{\rm cond}\propto T^2/\rho^2$	White dwarf and degenerate red-giant cores: conduction dominates, flattening the gradient toward isothermality.
Transition (onset of degeneracy)	$\kappa_{\rm cond}\sim\kappa_{\rm rad}$	Both channels comparable; the full harmonic combination is needed. This crossover density marks the edge of the isothermal region.
Isothermal limit	$\kappa\to 0\ \Rightarrow\ \mathrm{d}T/\mathrm{d}r\to 0$	When conduction is overwhelmingly efficient the gradient carrying $L_r$ vanishes: a uniform-temperature reservoir.
Ultra-relativistic degeneracy	scaling softens toward $\kappa_{\rm cond}\propto T^2/\rho^{4/3}$	In the densest cores $v_F\to c$ ; conduction stays dominant but its density dependence weakens.

Common trapDo not add the opacities — the fluxes add, so the opacities add inversely

The single most common error is to combine the channels by summing their opacities, $\kappa=\kappa_{\rm rad}+\kappa_{\rm cond}$ . That is backwards. The two channels carry heat in parallel under the same gradient, so it is the fluxes that add, and each flux is proportional to $1/\kappa$ — hence the inverse opacities sum: $1/\kappa=1/\kappa_{\rm rad}+1/\kappa_{\rm cond}$ , like parallel conductances. Adding the opacities directly treats the channels as series resistances, making the more efficient channel raise the total resistance instead of lowering it. Two further traps: conduction is negligible in normal stars — invoke it only where the electron gas is degenerate — and the conductivity here is the degenerate electron conductivity, which rises with density, not the classical Spitzer conductivity of a dilute plasma, which has the opposite density dependence and gives a qualitatively wrong crossover.

Worked examples

1 · Is a white dwarf core actually degenerate? (confidence builder)

A carbon–oxygen white dwarf core has $\rho=10^6\ \mathrm{g\,cm^{-3}}=10^9\ \mathrm{kg\,m^{-3}}$ , $\mu_e=2$ , and $T=10^7$ K. Compute the electron Fermi energy and the degeneracy parameter to confirm the electrons are strongly degenerate.

Step 1. Electron number density: $n_e=\rho/(\mu_e m_u)=10^9/(2\times1.6605\times10^{-27})=\mathbf{3.01\times10^{35}\ m^{-3}}$ — a hundred million times denser than a metal.

Step 2. Fermi energy from equation (8), with $(3\pi^2 n_e)^{1/3}=2.07\times10^{12}\ \mathrm{m^{-1}}$ : $E_F=2.63\times10^{-14}\ \mathrm{J}=\mathbf{164\ keV}$ . (The Fermi momentum gives $p_F/m_ec=0.80$ — mildly relativistic, hinting the $\kappa_{\rm cond}$ scaling softens at higher density.)

Step 3. Degeneracy parameter: $E_F/k_BT=2.63\times10^{-14}/(1.381\times10^{-23}\times10^7)=\mathbf{190}$ . With $E_F\gg k_BT$ by two orders of magnitude the gas is strongly degenerate. The core is degenerate with $E_F/k_BT\approx190$ , so conduction, not radiation, governs its heat transport.

2 · The combined opacity in the degenerate interior

In the same core the radiative opacity is set by electron scattering, $\kappa_{\rm rad}\approx0.2\ \mathrm{cm^2\,g^{-1}}$ , and the tabulated conductive opacity is $\kappa_{\rm cond}\approx5\times10^{-5}\ \mathrm{cm^2\,g^{-1}}$ . Find the effective opacity.

Step 1. Harmonic combination from equation (5): $\dfrac{1}{\kappa}=\dfrac{1}{0.2}+\dfrac{1}{5\times10^{-5}}=5+20000=20005\ \mathrm{g\,cm^{-2}}$ .

Step 2. Invert: $\kappa=1/20005=\mathbf{5.0\times10^{-5}\ cm^2\,g^{-1}}$ — indistinguishable from $\kappa_{\rm cond}$ ; the radiative term adds only $5$ parts in $20005$ ( $0.025\%$ ).

Step 3. The ratio $\kappa_{\rm rad}/\kappa=0.2/5\times10^{-5}=\mathbf{4000}$ : the effective opacity is 4000 times smaller than radiation alone, so the interior needs 4000 times less gradient to carry the same luminosity — the meaning of "conduction short-circuits radiation". $\kappa\approx\kappa_{\rm cond}=5.0\times10^{-5}\ \mathrm{cm^2\,g^{-1}}$ , set by the conductive channel.

3 · How isothermal is the core? (fully worked)

Using the numbers above, estimate the gradient conduction must sustain to carry $L_r\approx10^{-3}L_\odot=3.8\times10^{23}$ W at $r\approx3\times10^6$ m, and the temperature drop across the core.

Step 1. Local flux: $F=L_r/(4\pi r^2)=3.8\times10^{23}/[4\pi(3\times10^6)^2]=\mathbf{3.4\times10^{9}\ W\,m^{-2}}$ .

Step 2. Gradient from equation (6) with $\kappa=\kappa_{\rm cond}=5\times10^{-6}\ \mathrm{m^2\,kg^{-1}}$ , $\rho=10^9\ \mathrm{kg\,m^{-3}}$ , $T=10^7$ K: with $3\kappa\rho=1.5\times10^{4}$ and $4ac\,T^3=9.1\times10^{14}$ , $\left|\dfrac{\mathrm{d}T}{\mathrm{d}r}\right|=\dfrac{1.5\times10^4}{9.1\times10^{14}}\times3.4\times10^9=\mathbf{0.055\ K\,m^{-1}}$ .

Step 3. Integrated over the core radius, $\Delta T\approx0.055\times3\times10^6=\mathbf{1.7\times10^5\ K}$ , only $\Delta T/T\approx\mathbf{2\%}$ of the $10^7$ K interior. Radiation alone (4000 times larger $\kappa$ ) would demand $\Delta T\sim7\times10^8$ K, larger than $T$ itself and impossible. The degenerate core is nearly isothermal ( $\Delta T/T\approx2\%$ ); conduction is what makes such a flat interior possible.

Physics of the portrait

Conduction in stars is the quantum theory of metals transplanted to the cosmos; its decisive consequence — the isothermal degenerate core and its cooling law — was crystallized by one theorist in the 1950s.

Origins · 1 minLeon Mestel and how this was discovered

Leon Mestel

1927–2017 · white dwarf cooling and the isothermal degenerate core

Leon Mestel was a British astrophysicist whose 1952 paper "On the theory of white dwarf stars" turned qualitative ideas about degenerate matter into a working cooling theory — and it turned on conduction. Ralph Fowler had shown in 1926 that white dwarf matter is a degenerate electron gas, and Arnold Sommerfeld's quantum free-electron theory of metals (1927–28) had given the heat conductivity of such a gas, but no one had assembled these into a predictive model of how a white dwarf loses its heat. Mestel's insight was that the degenerate interior conducts so efficiently it is essentially isothermal — a vast reservoir of thermal energy stored in the ions, whose heat leaks out only through a thin, radiatively opaque envelope. That envelope, not the interior, controls the luminosity, giving the Mestel cooling law in which $L$ declines as the reservoir cools, so the coolest, faintest white dwarfs become clocks for dating stellar populations. The same isothermal-core physics governs the onset of the helium flash. Mestel later did foundational work on stellar magnetism and star formation, but his cooling theory remains the frame in which conduction earns its place in stellar structure.

History

1926

— Fowler shows dense white dwarf matter is a degenerate electron gas

1928

— Sommerfeld's quantum theory of metals yields the conductivity of a degenerate gas

1952

— Mestel derives white dwarf cooling from an isothermal, conduction-dominated core

1969

— Hubbard and Lampe tabulate conductive opacities for degenerate stellar matter

2007

— Cassisi and collaborators publish modern relativistic conductive opacity tables

Going further

Texts. Carroll and Ostlie, An Introduction to Modern Astrophysics, 2nd ed., §10.4 (energy transport) and §16.2 (white dwarfs) · Kippenhahn, Weigert and Weiss, Stellar Structure and Evolution, 2nd ed., Ch. 5 · Hansen, Kawaler and Trimble, Stellar Interiors, 2nd ed., Ch. 4 · Prialnik, An Introduction to the Theory of Stellar Structure and Evolution, 2nd ed., Ch. 5. Key values used throughout: $a=7.566\times10^{-16}\ \mathrm{J\,m^{-3}\,K^{-4}}$ ; the parallel-flux argument yields $1/\kappa=1/\kappa_{\rm rad}+1/\kappa_{\rm cond}$ ; the degenerate scaling is $\kappa_{\rm cond}\propto T^2/\rho^2$ ; the Fermi energy is $E_F=\tfrac{\hbar^2}{2m_e}(3\pi^2 n_e)^{2/3}$ . For a $10^6\ \mathrm{g\,cm^{-3}}$ , $10^7$ K carbon–oxygen core these give $E_F/k_BT\approx190$ , combined $\kappa=5.0\times10^{-5}\ \mathrm{cm^2\,g^{-1}}$ , and $\Delta T/T\approx2\%$ . Conductive opacity tables are from Hubbard and Lampe (1969) and Cassisi et al. (2007).

Builds on: Opacity and Radiative Transfer in Stars · Radiative Temperature Gradient · Equation of State for Stellar Matter · Equations of Stellar Structure. Leads to: Convective Energy Transport · Equations of Stellar Structure · Energy Generation Equation.

Ask the tutor: Ask me to walk through why a white dwarf core is nearly isothermal, or to compare the conductive and radiative opacities at any density and temperature you choose.

Why do the fluxes add but the opacities combine inversely?

Because radiation and conduction are two parallel routes for the same heat, driven by the same gradient. Heat can leave a shell by either route, so the total flux is the sum of the two fluxes. The subtlety is that each flux is proportional to $1/\kappa$ , not $\kappa$ : opacity is a resistance, so a bigger opacity means less flux. Adding two fluxes each proportional to an inverse opacity means $1/\kappa$ is what adds — like resistors in parallel, where conductances add. Summing the opacities directly is the series rule, the wrong picture here.

Why is conduction negligible in the Sun but dominant in a white dwarf?

It comes down to mean free paths. In the Sun's non-degenerate interior, photons travel much farther between interactions than electrons and move at $c$ , so radiative diffusion carries essentially all the flux. In a white dwarf the gas is degenerate: the Pauli principle blocks almost all scattering, so electron mean free paths become very long and the conductivity soars while the radiative opacity stays high. The electron channel then opens far wider than the photon channel and conduction takes over. The switch is tied to the onset of degeneracy, $E_F\gtrsim k_BT$ .

What does it mean that the conductive opacity falls with density?

An opacity is a resistance to heat flow, so a falling $\kappa_{\rm cond}$ means the gas conducts better the denser it is. Compressing a degenerate gas adds electrons and raises the Fermi energy, so the carriers move faster and their relative mean free paths grow; the conductivity rises roughly linearly with density, and its inverse, the opacity, falls. This is the opposite of a Kramers radiative opacity, which rises with density — and the two opposite slopes guarantee that at some crossover density conduction overtakes radiation.

Why does efficient conduction make the core isothermal rather than just cooler?

Because the gradient a star needs is set by how hard the heat is to move: the transport equation says $\mathrm{d}T/\mathrm{d}r\propto\kappa$ , so pushing a fixed luminosity through a low-opacity medium takes only a shallow gradient. In a degenerate core $\kappa_{\rm cond}$ is thousands of times below the radiative value, so the gradient collapses by the same factor — a drop of a percent or two across the whole core. The core still cools over time, but as a unit: a well-mixed reservoir whose slow leak through the insulating envelope is what the cooling law tracks.

Is this the same conductivity as in an ordinary metal or a laboratory plasma?

It is the same physics as a metal — the degenerate free-electron gas — but not the classical Spitzer conductivity of a hot dilute plasma. In a metal or a white dwarf, degeneracy restricts scattering to electrons near the Fermi surface, giving Wiedemann–Franz behaviour with conductivity rising with density. In a dilute plasma, Spitzer conductivity scales as $T^{5/2}$ and decreases with density. Using the Spitzer formula in a degenerate core gets the wrong sign for how conduction responds to compression and misses the crossover. The correct starting point is the degenerate electron gas, corrected for electron–ion scattering and, at high density, relativity.

How is conduction put into a stellar-structure code in practice?

Very cleanly, which is the point. Because the conductive flux is written in the same diffusion form as radiation, a modeller never adds a new equation. At each point the code looks up $\kappa_{\rm rad}$ and $\kappa_{\rm cond}$ from precomputed tables (functions of $\rho$ , $T$ , composition), combines them harmonically into a single $\kappa$ , and inserts it into the ordinary gradient equation. The isothermal core then emerges wherever $\kappa_{\rm cond}$ is tiny. The only real work is generating accurate conductive-opacity tables — a hard dense-plasma many-body problem — which is why references like Hubbard–Lampe and Cassisi et al. matter so much.