White Dwarfs and Electron Degeneracy Pressure

A white dwarf is the exposed, cooling core a Sun-like star leaves behind — the mass of a star crammed into the volume of a planet, roughly a solar mass inside an Earth-sized sphere. Once nuclear burning stops, there is no more heat to generate the thermal pressure that held the star up for billions of years, and nothing should stop gravity from crushing the core to a point. Yet it stops, held rigid by a pressure that owes nothing to temperature and everything to quantum mechanics: electron degeneracy pressure. The mechanism is the Pauli exclusion principle. Electrons are fermions, so no two may occupy the same quantum state. Squeeze the gas and you cannot pile electrons into the low-momentum states that are already full — they are forced up into ever-higher momentum levels, and that irreducible zero-point motion presses outward. The pressure depends only on how densely the electrons are packed, not on how hot they are: P∝ρ5/3P \propto \rho^{5/3} when they move slowly, softening to P∝ρ4/3P \propto \rho^{4/3} once the fastest approach the speed of light. Two consequences follow. First, balancing this pressure against gravity gives a mass–radius relation that runs backwards — R∝M−1/3R \propto M^{-1/3}, so a heavier white dwarf is a smaller one. Second, because the relativistic exponent 4/34/3 exactly matches what gravity demands, the pressure eventually loses the race: there is a maximum mass, the Chandrasekhar mass MCh≈1.44 M⊙M_{\mathrm{Ch}} \approx 1.44\,M_\odot, above which no white dwarf can exist. That single number, built from ℏ\hbar, cc, and GG alone, sets the fate of every dying low-mass star and lights the standard candles that measure the universe.

Key takeaways
  • A white dwarf is the degenerate remnant of a star born below ∼8 M⊙\sim 8\,M_\odot, held up not by heat but by the quantum degeneracy pressure of its electrons.
  • The Pauli exclusion principle forces compressed electrons into higher-momentum states, giving a pressure that depends on density alone: P∝ρ5/3P \propto \rho^{5/3} non-relativistic, P∝ρ4/3P \propto \rho^{4/3} ultra-relativistic.
  • Balancing degeneracy pressure against gravity gives the backwards relation R∝M−1/3R \propto M^{-1/3}: pile on mass and the star shrinks.
  • Because the ultra-relativistic exponent is exactly 4/34/3, pressure cannot outrun gravity above the Chandrasekhar mass MCh≈1.44 (2/μe)2 M⊙M_{\mathrm{Ch}} \approx 1.44\,(2/\mu_e)^2\,M_\odot.
  • A typical white dwarf has M∼0.6 M⊙M \sim 0.6\,M_\odot, R∼R⊕R \sim R_\oplus, central density ∼109 kg m−3\sim 10^9\,\mathrm{kg\,m^{-3}}, and Fermi energy ∼\sim MeV, cooling for billions of years along a track in the Hertzsprung–Russell diagram.
  • Key idea: the Pauli exclusion principle turns compression into pressure that depends only on density, giving white dwarfs the inverted relation R∝M−1/3R \propto M^{-1/3} and — once electrons turn relativistic — an absolute maximum mass fixed by fundamental constants.
The core insight

Electron degeneracy pressure arises because the Pauli exclusion principle forbids two electrons from sharing a quantum state; compression forces them into higher-momentum states, generating a pressure that depends on density alone — which gives white dwarfs the paradoxical relation R∝M−1/3R \propto M^{-1/3} and, once the electrons turn relativistic, an absolute maximum mass.

Where this lives today

Electron degeneracy is a working tool across modern astrophysics. Gaia's survey of billions of stars has resolved the white dwarf mass–radius relation to percent precision, and the sample now shows a "crystallization sequence": at low temperatures, the carbon–oxygen ions freeze into a Coulomb lattice, releasing latent heat that stalls the cooling exactly as degenerate-matter theory predicts. Because a white dwarf's cooling rate is well modelled, the faintest, oldest white dwarfs act as a chronometer — white-dwarf cosmochronology dates the Galactic disk and the oldest globular clusters independently of the cosmic expansion. The Chandrasekhar limit is the linchpin of precision cosmology: a carbon–oxygen white dwarf accreting toward MChM_{\mathrm{Ch}} ignites carbon and detonates as a Type Ia supernova, and because the trigger mass is nearly universal, these are the standardizable candles that revealed cosmic acceleration in 1998. Open questions remain lively: whether Type Ia progenitors are single- or double-degenerate, why some explode below MChM_{\mathrm{Ch}}, and how ultramassive white dwarfs near 1.35 M⊙1.35\,M_\odot probe the stiffest equation of state.

The physical picture

Confine a cold gas of electrons and quantum mechanics changes the rules. Each electron occupies a phase-space cell of volume h3h^3, and only two electrons — one spin up and one spin down — may share one cell. At white dwarf densities, the electrons run out of low-momentum cells, so even as the temperature falls toward zero, they are forced to fill states up to a large Fermi momentum pFp_F. The topmost electrons move at a sizeable fraction of the speed of light, not because the gas is hot but because there is nowhere slower left to go. That relentless, temperature-independent motion is the pressure, and the first figure shows its origin: filling momentum space to the brim.

Figure
The degenerate Fermi sea in momentum space. At T = 0, electrons fill every state out to the Fermi momentum p_F = ħ(3π²n_e)^⅓ (solid blue sphere). Compressing the gas 8-fold doubles p_F, since p_F ∝ n_e^⅓; once p_F crosses m_e c (dashed amber) the fastest electrons turn relativistic and the equation of state softens from P ∝ n_e^(5/3) to P ∝ n_e^(4/3)
Compression fills higher momenta — the origin of degeneracy pressure p_F m_e c density n_e p_F below m_e c : non-relativistic P ∝ n_e^(5/3) p_F density 8 n_e p_F above m_e c : relativistic P ∝ n_e^(4/3) → ×8 denser

The mass–radius relation is the part that feels wrong. Add mass to a planet and it grows; add mass to a white dwarf and it shrinks. The reason is a mismatch of exponents: extra mass deepens the gravitational well, scaling the required central pressure as M2/R4M^2/R^4, while non-relativistic degeneracy pressure stiffens only as ρ5/3∼(M/R3)5/3\rho^{5/3} \sim (M/R^3)^{5/3} — equating the two forces the radius down as mass rises. Keep going and the interior densifies until the electrons turn relativistic; the pressure then softens to ρ4/3\rho^{4/3}, which scales with radius exactly as gravity does, and the star loses its ability to push back at all. The second figure is the punchline: the mass–radius curve diving to zero at the Chandrasekhar mass.

Figure
The white dwarf mass–radius relation: heavier means smaller. At low mass R ∝ M^(−1/3) (the non-relativistic polytrope), but as electrons turn relativistic, the radius plunges to zero at the Chandrasekhar mass M_Ch = 1.44 M_⊙. Sirius B, at M = 1.02 M_⊙ and R ≈ 0.9 Earth radii, sits on the curve
The punchline: R(M) for white dwarfs — R ∝ M^(−1/3), vanishing at M_Ch = 1.44 M_⊙ 0 0.5 1.0 1.5 2.0 0.4 0.8 1.2 1.6 white dwarf mass M (M_⊙) radius R (Earth radii) M_Ch = 1.44 M_⊙ R → 0 Sirius B 1.02 M_⊙, 0.9 R_⊕ R ∝ M^(−1/3)

The core equations

Pressure is the flux of momentum across a surface. For a gas filling every state up to pFp_F, we integrate that flux over the occupied momenta, isotropically averaged — hence the factor 1/31/3.

Pdeg=13∫0pFp v(p) n(p) dp(1)P_{\mathrm{deg}} = \dfrac{1}{3}\int_0^{p_F} p\, v(p)\, n(p)\, \mathrm{d}p \qquad (1)

where: PdegP_{\mathrm{deg}} — electron degeneracy pressure; pFp_F — Fermi momentum, the highest occupied momentum; v(p)v(p) — speed of an electron of momentum pp; n(p) dp=(8πp2/h3) dpn(p)\,\mathrm{d}p = (8\pi p^2/h^3)\,\mathrm{d}p — number density of electrons with momentum in [p, p+dp][p,\,p+\mathrm{d}p], two spins per phase-space cell h3h^3. Everything downstream is this integral evaluated with the right v(p)v(p). The Fermi momentum is fixed entirely by how densely the electrons are packed. Integrating the density of states up to pFp_F and inverting gives the single most important microscopic relation of the subject.

pF=ℏ (3π2ne)1/3(2)p_F = \hbar\,(3\pi^2 n_e)^{1/3} \qquad (2)

where: ne=ρ/(μemH)n_e = \rho/(\mu_e m_H) — electron number density; ρ\rho — mass density; μe\mu_e — mean molecular weight per electron (≈2\approx 2 for a carbon–oxygen white dwarf); ℏ\hbar — reduced Planck constant. Compress the gas and pFp_F rises as the cube root of density — slowly, but without bound. When the electrons move slowly, v=p/mev = p/m_e, and the integral yields a stiff polytrope of index 5/35/3 — the pressure supporting ordinary, lower-mass white dwarfs.

PNR=(3π2)2/35ℏ2me ne5/3(3)P_{\mathrm{NR}} = \dfrac{(3\pi^2)^{2/3}}{5}\dfrac{\hbar^2}{m_e}\,n_e^{5/3} \qquad (3)

where: mem_e — electron mass; nen_e — electron number density. No TT appears anywhere — the defining signature of degeneracy. When the fastest electrons approach the speed of light, v→cv \to c, the extra power of speed is lost and the integral softens to a polytrope of index 4/34/3 — the regime that breeds the mass limit.

PUR=(3π2)1/34 ℏc ne4/3(4)P_{\mathrm{UR}} = \dfrac{(3\pi^2)^{1/3}}{4}\,\hbar c\,n_e^{4/3} \qquad (4)

where: cc — speed of light; nen_e — electron number density. The exponent has dropped from 5/35/3 to 4/34/3: the pressure grows more sluggishly with density, and that small change is decisive. The star holds itself up by balancing this pressure gradient against its own weight — the equation of hydrostatic equilibrium, unchanged from ordinary stellar structure.

dPdr=−G m(r) ρ(r)r2(5)\dfrac{\mathrm{d}P}{\mathrm{d}r} = -\dfrac{G\,m(r)\,\rho(r)}{r^2} \qquad (5)

where: m(r)m(r) — mass enclosed within radius rr; ρ(r)\rho(r) — local density; GG — gravitational constant. Combined with a polytropic equation of state, this becomes the Lane–Emden equation. Balancing the non-relativistic equation of state against gravity gives the mass–radius relation — inverted, so that heavier white dwarfs are smaller.

R∝M−1/3(non-relativistic)(6)R \propto M^{-1/3} \qquad (\text{non-relativistic}) \qquad (6)

where: RR — stellar radius; MM — stellar mass. A 1 M⊙1\,M_\odot white dwarf is smaller than a 0.5 M⊙0.5\,M_\odot one — the opposite of every planet and main-sequence star. Pushed into the relativistic regime, the mass becomes independent of radius and locks onto a single value built from fundamental constants — the Chandrasekhar mass.

MCh≈1.457(2μe)2M⊙≈1.44 M⊙(7)M_{\mathrm{Ch}} \approx 1.457\left(\dfrac{2}{\mu_e}\right)^{2} M_\odot \approx 1.44\,M_\odot \qquad (7)

where: μe\mu_e — mean molecular weight per electron (=2=2 for a C–O composition); M⊙=1.989×1030 kgM_\odot = 1.989\times10^{30}\,\mathrm{kg}. The coefficient 1.457 M⊙1.457\,M_\odot is the ideal n=3n=3 polytrope value; general-relativistic and Coulomb corrections trim it to the familiar 1.44 M⊙1.44\,M_\odot. Finally, whether a white dwarf is non-relativistic or relativistic is decided by the kinetic energy at the Fermi surface, which interpolates smoothly between the two limits.

EF=(pFc)2+(mec2)2−mec2(8)E_F = \sqrt{(p_F c)^2 + (m_e c^2)^2} - m_e c^2 \qquad (8)

where: EFE_F — kinetic energy at the Fermi surface; mec2=0.511 MeVm_e c^2 = 0.511\,\mathrm{MeV} — electron rest energy. When pFc≪mec2p_F c \ll m_e c^2 the star is non-relativistic; when pFc≳mec2p_F c \gtrsim m_e c^2 it is on the road to the Chandrasekhar limit.

Derivation ADegeneracy pressure and the two equations of state

The goal is to turn the single momentum-flux integral of equation (1) into the two polytropic equations of state that govern white dwarfs — the stiff P∝ne5/3P \propto n_e^{5/3} that supports light ones and the soft P∝ne4/3P \propto n_e^{4/3} that dooms heavy ones — by substituting the correct electron speed v(p)v(p) and eliminating the Fermi momentum in favour of density. Both limits come from the same integral; only the kinematics of the fastest electrons differ.

Step 1. Count the states. With two electrons per phase-space cell of volume h3h^3, the density of electrons with momentum up to pFp_F is ne=∫0pF(8πp2/h3) dp=(8π/3h3) pF3n_e = \int_0^{p_F} (8\pi p^2/h^3)\,\mathrm{d}p = (8\pi/3h^3)\,p_F^3; inverting gives pF=ℏ(3π2ne)1/3p_F = \hbar(3\pi^2 n_e)^{1/3}, equation (2). Why: this is the only place the Pauli principle enters — it forces a nonzero pFp_F at T=0T = 0 and hence a nonzero pressure.

Step 2. Do the non-relativistic case. Put v=p/mev = p/m_e into equation (1): PNR=13∫0pFp (p/me) (8πp2/h3) dp=(8π/15meh3) pF5P_{\mathrm{NR}} = \tfrac{1}{3}\int_0^{p_F} p\,(p/m_e)\,(8\pi p^2/h^3)\,\mathrm{d}p = (8\pi/15 m_e h^3)\,p_F^5. Why: for slow electrons, momentum and speed are proportional, so the integrand is ∝p4\propto p^4 and the integral scales as pF5p_F^5 — the extra powers of pFp_F are what make this equation of state stiff.

Step 3. Eliminate pFp_F. Substitute pF=ℏ(3π2ne)1/3p_F = \hbar(3\pi^2 n_e)^{1/3} and collect constants to reach equation (3), PNR=(3π2)2/35ℏ2mene5/3P_{\mathrm{NR}} = \tfrac{(3\pi^2)^{2/3}}{5}\tfrac{\hbar^2}{m_e} n_e^{5/3}. Why: writing pressure in terms of density exposes the polytropic index γ=5/3\gamma = 5/3 (Lane–Emden n=3/2n = 3/2) — the input to the mass–radius relation.

Step 4. Do the ultra-relativistic case. Now the fastest electrons satisfy p≫mecp \gg m_e c, so v→cv \to c: PUR=13∫0pFp c (8πp2/h3) dp=(2πc/3h3) pF4P_{\mathrm{UR}} = \tfrac{1}{3}\int_0^{p_F} p\,c\,(8\pi p^2/h^3)\,\mathrm{d}p = (2\pi c/3h^3)\,p_F^4, and pF∝ne1/3p_F \propto n_e^{1/3} gives equation (4). Why: replacing p/mep/m_e by the constant cc removes one power of pp, so pressure scales as pF4p_F^4 not pF5p_F^5 — the softening to γ=4/3\gamma = 4/3 that eventually lets gravity win.

Go deeper · Derivation B · 1 minThe inverse mass–radius relation from the Lane–Emden polytrope
Derivation BThe inverse mass–radius relation from the Lane–Emden polytrope

The goal is to show why heavier white dwarfs are smaller, deriving R∝M−1/3R \propto M^{-1/3} from the non-relativistic equation of state. A polytrope P=Kρ1+1/nP = K\rho^{1+1/n} in hydrostatic equilibrium reduces to the dimensionless Lane–Emden equation, whose solution ties both radius and mass to the central density; eliminating that shared density produces the inverse scaling — a result a quick scaling argument, run first, already reveals.

Step 1. Scaling shortcut. The central pressure needed to support the star is, from hydrostatic balance, Pc∼GM2/R4P_c \sim GM^2/R^4; the equation of state supplies Pc∼Kρc5/3∼K(M/R3)5/3P_c \sim K\rho_c^{5/3} \sim K(M/R^3)^{5/3}. Equating, KM5/3/R5∼GM2/R4K M^{5/3}/R^5 \sim GM^2/R^4, which rearranges to R∝KM−1/3/GR \propto K M^{-1/3}/G. Why: with no calculus this already reveals the sign — more mass, smaller radius — because the pressure demanded by gravity (M2/R4M^2/R^4) outruns what a ρ5/3\rho^{5/3} gas can supply.

Step 2. Set up the polytrope properly. Write ρ=ρc θn\rho = \rho_c\,\theta^n with n=3/2n = 3/2 (since γ=1+1/n=5/3\gamma = 1 + 1/n = 5/3), and r=αξr = \alpha\xi with α2=(n+1)Kρc(1/n)−1/(4πG)\alpha^2 = (n+1)K\rho_c^{(1/n)-1}/(4\pi G). Hydrostatic equilibrium plus Poisson's equation collapse to the Lane–Emden equation 1ξ2ddξ ⁣(ξ2dθdξ)=−θn\tfrac{1}{\xi^2}\tfrac{\mathrm{d}}{\mathrm{d}\xi}\!\left(\xi^2\tfrac{\mathrm{d}\theta}{\mathrm{d}\xi}\right) = -\theta^n. Why: the substitution absorbs every physical constant into α\alpha, leaving one universal dimensionless equation whose solution θ(ξ)\theta(\xi) is computed once and reused for all polytropes of a given nn.

Step 3. Read off radius and mass. The surface is the first zero, θ(ξ1)=0\theta(\xi_1) = 0, so R=αξ1R = \alpha\xi_1 and M=4πα3ρc (−ξ2 dθ/dξ)ξ1M = 4\pi\alpha^3\rho_c\,(-\xi^2\,\mathrm{d}\theta/\mathrm{d}\xi)_{\xi_1}, both depending on central density only through α∝ρc(1/n−1)/2\alpha \propto \rho_c^{(1/n - 1)/2}. Why: ξ1\xi_1 and (−ξ2θ′)ξ1(-\xi^2\theta')_{\xi_1} are pure numbers for n=3/2n = 3/2, so the physics of a given star is carried by ρc\rho_c alone.

Step 4. Eliminate the central density. For n=3/2n = 3/2, tracking powers of ρc\rho_c gives R∝ρc−1/6R \propto \rho_c^{-1/6} and M∝ρc1/2M \propto \rho_c^{1/2}; solving the second for ρc∝M2\rho_c \propto M^2 and inserting into the first yields R∝M−1/3R \propto M^{-1/3}. Why: the shared ρc\rho_c is the hidden variable linking MM and RR; removing it leaves the observable relation, confirmed numerically at R≈0.013 R⊙R \approx 0.013\,R_\odot for M=0.6 M⊙M = 0.6\,M_\odot.

Go deeper · Derivation C · 1 minThe Chandrasekhar mass from the relativistic limit
Derivation CThe Chandrasekhar mass from the relativistic limit

The goal is to explain why a white dwarf cannot exceed about 1.44 M⊙1.44\,M_\odot — not as an empirical ceiling but as an inevitability of the relativistic equation of state. The ultra-relativistic polytrope has index n=3n = 3, for which the mass depends on the central density not at all: the Lane–Emden machinery returns a single mass, built from fundamental constants, at which pressure and gravity scale identically and equilibrium becomes a knife-edge.

Step 1. Use the soft equation of state. As density climbs and the electrons become relativistic across most of the star, the pressure follows equation (4): P=K2ρ4/3P = K_2\rho^{4/3} with K2=(3π2)1/3ℏc/4 (μemH)−4/3K_2 = (3\pi^2)^{1/3}\hbar c/4\,(\mu_e m_H)^{-4/3}, i.e. γ=4/3\gamma = 4/3 and Lane–Emden index n=3n = 3. Why: this is the equation of state a massive white dwarf is forced into, because its high density has pushed pFp_F well past mecm_e c throughout the interior.

Step 2. Track the central-density dependence of the mass. For a general polytrope M∝α3ρc∝ρc(3/2)(1/n−1)+1M \propto \alpha^3\rho_c \propto \rho_c^{(3/2)(1/n - 1) + 1}. Setting n=3n = 3 makes the exponent (3/2)(1/3−1)+1=0(3/2)(1/3 - 1) + 1 = 0, so M∝ρc0M \propto \rho_c^{0} — the mass is independent of central density and hence of radius. Why: this is the mathematical heart of the limit; for n=3n = 3 every equilibrium configuration, whatever its size, has exactly the same mass, so that mass is uniquely determined.

Step 3. Evaluate the constant. With n=3n = 3 the mass is M=4π (K2/πG)3/2(−ξ2θ′)ξ1M = 4\pi\,(K_2/\pi G)^{3/2}(-\xi^2\theta')_{\xi_1}; inserting K2K_2, the n=3n=3 Lane–Emden number (−ξ2θ′)ξ1=2.018(-\xi^2\theta')_{\xi_1} = 2.018, and μe=2\mu_e = 2 gives MCh=1.457 (2/μe)2 M⊙≈1.44 M⊙M_{\mathrm{Ch}} = 1.457\,(2/\mu_e)^2\,M_\odot \approx 1.44\,M_\odot. Why: the result scales as (ℏc/G)3/2/(μemH)2(\hbar c/G)^{3/2}/(\mu_e m_H)^2 — quantum, relativistic, and gravitational constants with no free parameters, which is why the limit is universal.

Step 4. Interpret the stability. At γ=4/3\gamma = 4/3 the star is marginally stable: compression changes pressure and gravitational demand in lockstep, so nothing restores it. For M>MChM > M_{\mathrm{Ch}} gravity wins outright and the core collapses to a neutron star or, if carbon ignites first, detonates as a Type Ia supernova. Why: the exponent 4/34/3 is the boundary between stable and unstable stars — what makes the Chandrasekhar mass a hard wall, not a soft trend.

Limiting cases

The degenerate equation of state has two clean regimes; the classical limits confirm degeneracy is a genuinely distinct source of pressure.

CaseResultMeaning
Non-relativistic degenerateP=K1(ρμemH)5/3, K1=(3π2)2/3ℏ25meP = K_1\left(\dfrac{\rho}{\mu_e m_H}\right)^{5/3},\ K_1 = \dfrac{(3\pi^2)^{2/3}\hbar^2}{5m_e}Low-mass white dwarfs, pF≪mecp_F \ll m_e c; polytrope n=3/2n = 3/2 gives R∝M−1/3R \propto M^{-1/3}.
Ultra-relativistic degenerateP=K2(ρμemH)4/3, K2=(3π2)1/3ℏc4P = K_2\left(\dfrac{\rho}{\mu_e m_H}\right)^{4/3},\ K_2 = \dfrac{(3\pi^2)^{1/3}\hbar c}{4}High-mass white dwarfs, pF≫mecp_F \gg m_e c; polytrope n=3n = 3, mass independent of radius — the Chandrasekhar limit.
Strong-degeneracy validitykBT≪EFk_B T \ll E_FDegeneracy dominates when the Fermi energy swamps thermal energy; in white dwarfs EF∼E_F \sim MeV, kBT∼k_B T \sim keV, holding by ∼102\sim 10^{2}–10310^{3}.
Classical ideal-gas limitP=ρkBTμmHP = \dfrac{\rho k_B T}{\mu m_H}Recovered at low density or high temperature; this supports main-sequence stars, not white dwarfs.

Worked examples

1 · Is Sirius B relativistic? (fully worked)

Sirius B has M=1.02 M⊙M = 1.02\,M_\odot and radius R=0.0084 R⊙≈R = 0.0084\,R_\odot \approx one Earth radius. Estimate its mean density, electron number density, Fermi momentum, the relativity parameter x=pF/mecx = p_F/m_e c, the Fermi kinetic energy, and confirm that the gas is strongly degenerate at T≈2×107 KT \approx 2\times10^7\,\mathrm{K}. Take μe=2\mu_e = 2.

Step 1. Mean density. With M=2.03×1030 kgM = 2.03\times10^{30}\,\mathrm{kg} and R=0.0084×6.957×108=5.84×106 mR = 0.0084\times6.957\times10^{8} = 5.84\times10^{6}\,\mathrm{m}, the volume is 43πR3=8.34×1020 m3\tfrac{4}{3}\pi R^3 = 8.34\times10^{20}\,\mathrm{m^3}, so ρˉ=M/V=2.4×109 kg m−3\bar\rho = M/V = \mathbf{2.4\times10^{9}\ kg\,m^{-3}} — a million times denser than any laboratory solid (the central density is ∼6×\sim 6\times higher).

Step 2. Electron density. ne=ρˉ/(μemH)=2.4×109/(2×1.673×10−27)=7.2×1035 m−3n_e = \bar\rho/(\mu_e m_H) = 2.4\times10^{9}/(2\times1.673\times10^{-27}) = \mathbf{7.2\times10^{35}\ m^{-3}}.

Step 3. Fermi momentum. pF=ℏ(3π2ne)1/3p_F = \hbar(3\pi^2 n_e)^{1/3}; with 3π2ne=2.1×10373\pi^2 n_e = 2.1\times10^{37} and cube root 2.77×1012 m−12.77\times10^{12}\,\mathrm{m^{-1}}, pF=1.055×10−34×2.77×1012=2.9×10−22 kg m s−1p_F = 1.055\times10^{-34}\times2.77\times10^{12} = \mathbf{2.9\times10^{-22}\ kg\,m\,s^{-1}}.

Step 4. Relativity parameter. x=pF/(mec)=2.9×10−22/2.73×10−22=1.07x = p_F/(m_e c) = 2.9\times10^{-22}/2.73\times10^{-22} = \mathbf{1.07}. Because x≈1x \approx 1, Sirius B is mildly relativistic — its electrons straddle the ρ5/3\rho^{5/3} and ρ4/3\rho^{4/3} regimes, which is why a ∼1 M⊙\sim 1\,M_\odot white dwarf sits where the mass–radius curve begins to bend toward the Chandrasekhar limit.

Step 5. Fermi energy. pFc=8.76×10−14 J=0.547 MeVp_F c = 8.76\times10^{-14}\,\mathrm{J} = 0.547\,\mathrm{MeV}, so EF=0.5472+0.5112−0.511=0.24 MeVE_F = \sqrt{0.547^2 + 0.511^2} - 0.511 = \mathbf{0.24\ MeV}.

Step 6. Degeneracy check. At T=2×107 KT = 2\times10^7\,\mathrm{K}, kBT=1.7×10−3 MeVk_B T = 1.7\times10^{-3}\,\mathrm{MeV}, so kBT/EF≈0.007k_B T/E_F \approx \mathbf{0.007}. The gas is strongly degenerate and the pressure temperature-independent: a "cold" quantum object even though its surface glows at 25,000 K25{,}000\,\mathrm{K}.

2 · The backwards mass–radius relation

A 0.6 M⊙0.6\,M_\odot carbon–oxygen white dwarf has R≈0.013 R⊙R \approx 0.013\,R_\odot. Using the non-relativistic scaling, predict the radius of a 1.0 M⊙1.0\,M_\odot white dwarf, and explain why the real value is even smaller.

Step 1. Apply R∝M−1/3R \propto M^{-1/3}: R1.0=0.013 R⊙×(1.667)−1/3=0.013×0.843=0.011 R⊙R_{1.0} = 0.013\,R_\odot \times (1.667)^{-1/3} = 0.013\times0.843 = \mathbf{0.011\ R_\odot}, about 1.21.2 Earth radii — a 67%67\% increase in mass shrinks the star by roughly 16%16\%.

Step 2. Central density scaling. Since ρc∝M/R3∝M2\rho_c \propto M/R^3 \propto M^2, going from 0.60.6 to 1.0 M⊙1.0\,M_\odot raises the central density by (1.0/0.6)2=2.8×(1.0/0.6)^2 = \mathbf{2.8\times} — denser cores mean faster Fermi electrons.

Step 3. Why reality beats the prediction. The measured radius of Sirius B (≈0.0084 R⊙\approx 0.0084\,R_\odot) is smaller than the 0.011 R⊙0.011\,R_\odot the non-relativistic law predicts. From Example 1, at ∼1 M⊙\sim 1\,M_\odot the electrons are already relativistic (x≈1x \approx 1), so the equation of state has begun softening toward ρ4/3\rho^{4/3}; the softer pressure lets gravity compress the star further, pulling the radius below the M−1/3M^{-1/3} curve toward MChM_{\mathrm{Ch}}, where R→0R \to 0.

3 · How composition shifts the mass limit

A carbon–oxygen white dwarf has μe=2\mu_e = 2. An iron core would have μe≈2.15\mu_e \approx 2.15. By what factor does the Chandrasekhar mass change?

Step 1. From equation (7), MCh∝μe−2M_{\mathrm{Ch}} \propto \mu_e^{-2}, so the ratio is (2.15/2.00)−2=(1.075)−2=0.865(2.15/2.00)^{-2} = (1.075)^{-2} = \mathbf{0.865}.

Step 2. An iron white dwarf's limit is 1.44×0.865≈1.25 M⊙1.44\times0.865 \approx \mathbf{1.25\ M_\odot} — about 13%13\% lower. Fewer electrons per nucleon means less degeneracy pressure per unit mass, so gravity wins sooner; μe\mu_e is the only material property entering the mass limit.

Physics of the portrait

The white dwarf is where quantum statistics first reshaped astronomy — a story running from an unexplained wobble in the brightest star in the sky to a young man's shipboard calculation that fixed, from three constants of nature, the heaviest a dead star can be.

Origins · 1 minSubrahmanyan Chandrasekhar and how this was discovered
SC
Subrahmanyan Chandrasekhar
1910–1995 · relativistic degeneracy and the maximum mass of white dwarfs

White dwarfs had been a puzzle since Ralph Fowler showed in 1926 that their staggering densities could be explained by the new Fermi–Dirac statistics: the electrons form a degenerate gas whose quantum pressure, not heat, holds the star up. Fowler's picture seemed to permit a white dwarf of any mass. In 1930, on the long sea voyage from Madras to Cambridge to begin his doctorate, the nineteen-year-old Subrahmanyan Chandrasekhar carried it one crucial step further. He asked what happens when the density is so high that the Fermi electrons move near the speed of light — and found the equation of state softens from P∝ρ5/3P \propto \rho^{5/3} to P∝ρ4/3P \propto \rho^{4/3}. With that softer exponent, the mass of a degenerate star becomes independent of its radius and locks onto a single value, MCh≈1.4 M⊙M_{\mathrm{Ch}} \approx 1.4\,M_\odot, above which a star has no stable endpoint as a white dwarf. Arthur Eddington, the most influential astrophysicist of the age, ridiculed the result at the 1935 Royal Astronomical Society meeting, insisting nature would "find a way," and his prestige delayed acceptance for years. But the mathematics was right. The limit now bears Chandrasekhar's name, underpins the theory of neutron stars and black holes, and sets the trigger mass for Type Ia supernovae. He received the Nobel Prize in Physics in 1983.

History
1844
— Bessel infers that Sirius has an unseen massive companion from its wobble
1926
— Fowler explains white dwarf matter as a degenerate Fermi–Dirac electron gas
1930
— Chandrasekhar derives the relativistic maximum mass on his voyage to England
1935
— Eddington publicly disputes the mass limit; the community sides with him for years
1983
— Chandrasekhar receives the Nobel Prize for his work on stellar structure

Going further

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., Ch. 16 (Degenerate Remnants of Stars) · Chandrasekhar, An Introduction to the Study of Stellar Structure (1939) · Shapiro & Teukolsky, Black Holes, White Dwarfs, and Neutron Stars, Ch. 3–4 · Kippenhahn, Weigert & Weiss, Stellar Structure and Evolution, 2nd ed., Ch. 15 · Padmanabhan, Theoretical Astrophysics, Vol. II. The key results used throughout: pF=ℏ(3π2ne)1/3p_F = \hbar(3\pi^2 n_e)^{1/3}; equations of state P∝ρ5/3P \propto \rho^{5/3} (Lane–Emden n=3/2n = 3/2) and P∝ρ4/3P \propto \rho^{4/3} (Lane–Emden n=3n = 3); mass–radius scaling R∝M−1/3R \propto M^{-1/3}; Chandrasekhar's 1931 polytrope coefficient MCh=1.457 (2/μe)2 M⊙M_{\mathrm{Ch}} = 1.457\,(2/\mu_e)^2\,M_\odot (equivalently 5.83/μe2 M⊙5.83/\mu_e^2\,M_\odot), trimmed to ≈1.44 M⊙\approx 1.44\,M_\odot by general-relativistic and Coulomb corrections. Physical constants used: mec2=0.511 MeVm_e c^2 = 0.511\,\mathrm{MeV}, M⊙=1.989×1030 kgM_\odot = 1.989\times10^{30}\,\mathrm{kg}, R⊙=6.957×108 mR_\odot = 6.957\times10^8\,\mathrm{m}.

Builds on. Asymptotic Giant Branch and Thermal Pulses · Post-Main-Sequence Evolution and the Red Giant Branch · The Hertzsprung–Russell Diagram. Leads to. The Chandrasekhar Mass · Advanced Burning Stages in Massive Stars.

Ask the tutor: Ask me why the relativistic index 4/34/3 turns a stable star unstable, or to compute the Fermi energy and radius of any white dwarf from its mass.

Why does degeneracy pressure not depend on temperature?

Because the electrons fill every momentum state up to the Fermi momentum pFp_F regardless of how hot the gas is, as long as kBT≪EFk_B T \ll E_F. The pressure is set by that density-determined Fermi sea — the population of fast electrons the exclusion principle forces into existence — not by thermal agitation. Heating the gas only nudges a few electrons near the Fermi surface into higher states, a negligible change when EF∼E_F \sim MeV and kBT∼k_B T \sim keV. So the pressure persists unchanged as the star cools toward T=0T = 0, which is why a white dwarf does not shrink as it radiates away its heat over billions of years.

Why are more massive white dwarfs smaller?

It comes from a mismatch of exponents. Adding mass deepens the gravitational well, so the central pressure gravity demands scales as M2/R4M^2/R^4, while non-relativistic degeneracy pressure supplies only ρ5/3∼(M/R3)5/3\rho^{5/3} \sim (M/R^3)^{5/3}. Setting the two equal gives R∝M−1/3R \propto M^{-1/3}: to make the extra pressure a heavier star needs, the gas must compress to higher density, meaning a smaller radius — the opposite of a planet, where added mass just piles on volume. Degeneracy pressure is "soft," stiffening more slowly with density than gravity strengthens, so equilibrium is reached only by shrinking. Push far enough and the electrons turn relativistic, softening the pressure further and driving the radius to zero at the Chandrasekhar mass.

What is the mean molecular weight per electron μ_e, and why is it 2?

μe\mu_e is the gas mass, in atomic mass units, per free electron: ne=ρ/(μemH)n_e = \rho/(\mu_e m_H). For any fully ionized element heavier than hydrogen, each nucleus contributes ZZ electrons and a mass of roughly A≈2ZA \approx 2Z, so μe=A/Z≈2\mu_e = A/Z \approx 2. Carbon (A=12A=12, Z=6Z=6), oxygen (A=16A=16, Z=8Z=8), and helium (A=4A=4, Z=2Z=2) all give exactly 22; only hydrogen gives μe=1\mu_e = 1 — but white dwarfs have burned their hydrogen. Because the Chandrasekhar mass scales as μe−2\mu_e^{-2}, this single number ties the mass limit to composition.

Is a white dwarf held up like a solid, or is it still a gas?

It is a gas — but a quantum-degenerate one, unlike any classical gas. The electrons are not bound to atoms; they roam freely through a lattice of ionized nuclei, forming a Fermi sea that provides essentially all the pressure. In that sense a white dwarf is closer to a giant metal than to a rock: the electrons are delocalized and conduct heat superbly, so the interior is nearly isothermal. The ions carry the star's thermal energy and at low temperatures can crystallize into a Coulomb lattice, but the support against gravity is electron degeneracy pressure — a property of the quantum gas, not the mechanical strength of any solid.

What happens if a white dwarf gains mass beyond the Chandrasekhar limit?

It cannot remain a white dwarf. As accretion pushes the mass toward MChM_{\mathrm{Ch}}, the central density climbs until, in a carbon–oxygen star, carbon ignites — and because degeneracy pressure is temperature-independent, the burning does not trigger the expansion that would normally throttle it. The result is a thermonuclear runaway that unbinds the entire star as a Type Ia supernova. In other cases (an oxygen–neon core, or rapid accretion) electron captures soften the pressure and trigger collapse to a neutron star instead. Either way, MChM_{\mathrm{Ch}} is a genuine wall: above it, degeneracy pressure scales exactly as gravity does, so no stable equilibrium exists.

Why does the relativistic exponent 4/3 cause instability while 5/3 does not?

Stability against collapse requires the pressure to stiffen faster than gravity when the star is squeezed — quantitatively, an adiabatic index γ>4/3\gamma > 4/3. The non-relativistic gas has γ=5/3>4/3\gamma = 5/3 > 4/3, so a small compression raises the pressure more than it raises gravity's demand, and the star springs back: stable. But when electrons turn relativistic, γ→4/3\gamma \to 4/3 exactly, and now compression raises pressure and gravitational demand in perfect step — no restoring force, marginal equilibrium. Any additional mass tips the balance and gravity wins. The value 4/34/3 is the precise dividing line between a star that can support itself and one that must collapse, and it is what makes the Chandrasekhar mass a hard limit.

Before this

Where this leads