White Dwarfs

A white dwarf is what is left when a star like the Sun runs out of ways to burn its nuclear fuel. Stars below roughly 8 M⊙8\,M_\odot never ignite carbon; after the asymptotic-giant-branch phase blows off the envelope as a planetary nebula, the naked carbon–oxygen core is exposed with no fuel left to burn. Ordinary stars hold themselves up by being hot — remove the heat and they shrink. A white dwarf does not care. It is supported by something that survives to absolute zero: electron degeneracy pressure, the flat refusal of the Pauli exclusion principle to let two electrons share the same quantum state. Squeeze the gas and the electrons are forced into ever-higher momentum states, and the momentum flux of that filled Fermi sea pushes back — no heat required. That single fact reorganizes everything. Because the pressure is set by density rather than temperature, the equation of state is a polytrope, P=KργP=K\rho^{\gamma}, and the structure obeys the Lane–Emden equation. The exponent is γ=5/3\gamma=5/3 while the electrons are slow, but softens to γ=4/3\gamma=4/3 once the Fermi momentum turns relativistic. From this one softening flow the two celebrated results of the subject: the inverted mass–radius relation R∝M−1/3R\propto M^{-1/3}, in which heavier white dwarfs are smaller, and a hard ceiling on the mass — the Chandrasekhar limit MCh≈1.4 M⊙M_{\rm Ch}\approx1.4\,M_\odot, built entirely from ℏ\hbar, cc, GG and the baryon mass, with no reference to how the star was made. Above it no equilibrium exists, and gravity wins. A white dwarf is thus the closest thing astronomy has to a cold, quantum-stiffened solid with a size of the Earth and the mass of the Sun, radiating its stored heat over aeons toward a final darkness.

Key takeaways
  • White dwarfs are the degenerate remnants of stars below ∼8 M⊙\sim8\,M_\odot: the exposed carbon–oxygen core, held up not by thermal pressure but by electron degeneracy pressure that persists at T=0T=0.
  • Their structure is a polytrope; the equation of state softens from P∝ρ5/3P\propto\rho^{5/3} (non-relativistic) to P∝ρ4/3P\propto\rho^{4/3} (ultra-relativistic), and that softening is the hinge of the whole subject.
  • The mass–radius relation is inverted, R∝M−1/3R\propto M^{-1/3}: adding mass makes a white dwarf smaller and denser, opposite to an ordinary star.
  • There is a maximum mass — the Chandrasekhar mass MCh≈1.46 (2/μe)2 M⊙≈1.4 M⊙M_{\rm Ch}\approx1.46\,(2/\mu_e)^2\,M_\odot\approx1.4\,M_\odot, where μe\mu_e is the mean molecular weight per electron — above which no equilibrium exists and collapse is inevitable.
  • A typical 0.6 M⊙0.6\,M_\odot white dwarf is Earth-sized (R≈R⊕R\approx R_\oplus), with mean density ∼106 g cm−3\sim10^6\,\mathrm{g\,cm^{-3}} and surface gravity ∼108 cm s−2\sim10^8\,\mathrm{cm\,s^{-2}}.
  • With no fuel they simply cool: the Mestel law tcool∝L−5/7t_{\rm cool}\propto L^{-5/7} turns the faint end of the luminosity function into a clock for stellar populations.
  • Key idea: a white dwarf is a star where gravity is balanced by electron degeneracy pressure — a pressure independent of temperature — and because that pressure softens from γ=5/3\gamma=5/3 to γ=4/3\gamma=4/3 as electrons turn relativistic, white dwarfs shrink as they gain mass and cannot exist above the Chandrasekhar limit.
The core insight

A white dwarf is a star in mechanical equilibrium in which gravity is balanced by electron degeneracy pressure — a pressure that persists at zero temperature. Because this pressure depends on density (not temperature), and because it softens from γ=5/3\gamma=5/3 to γ=4/3\gamma=4/3 as electrons turn relativistic, white dwarfs shrink as they gain mass and cannot exist above a threshold: the Chandrasekhar limit.

Where this lives today

White dwarfs have gone from curiosities to precision laboratories. The Gaia mission has measured parallaxes for hundreds of thousands of them, and the resulting mass–radius data confirm the degenerate equation of state to a few percent — the theory above is now checked star by star. The faint end of the white-dwarf luminosity function serves as a cosmochronometer: the coolest, faintest dwarfs are the oldest, and the sharp cutoff dates the Galactic disk at 88–1010 billion years and gives independent ages for globular clusters. At the high-mass extreme, a carbon–oxygen white dwarf driven toward MChM_{\rm Ch} by accretion or merger ignites explosive carbon burning and detonates as a Type Ia supernova; because the ignition mass is nearly universal, these are the calibratable standard candles that revealed dark energy in 1998. Gaia photometry has even caught the "cooling pile-up" predicted when the core crystallizes and releases latent heat. Push a degenerate core past the limit a different way — by electron capture rather than detonation — and it collapses to a neutron star, tying white-dwarf physics to the densest matter known.

The physical picture

Start with the balance itself. Gravity pulls every gram toward the centre; against it stands the degeneracy pressure of a cold, dense electron gas. Compress that gas and the exclusion principle forces electrons to occupy higher momentum states — the Fermi momentum climbs as pF∝ne1/3p_F\propto n_e^{1/3} — and the outward momentum flux stiffens, and hence the degeneracy pressure. This stiffness has nothing to do with the thermal heat: a cold electron gas resists compression exactly as fiercely as a hot one. The first figure fixes the scale of what that balance produces — the Sun's worth of material crushed into the volume of the Earth.

Figure
A 0.6 M⊙ white dwarf packs the Sun's material into an Earth-sized ball — R ≈ 6400 km, mean density ρ̄ ≈ 1×10⁶ g cm⁻³, surface gravity g ≈ 1×10⁸ cm s⁻². Inward gravity is balanced by outward electron degeneracy pressure P ∝ ρ^(5/3), not by thermal pressure
A star the mass of the Sun, the size of the Earth — held up by quantum pressure degenerate C/O interior nearly isothermal ρ_c ≈ 3×10⁶ g cm⁻³ thin non-degenerate envelope (insulating skin) gravity in ▸ degeneracy pressure out ◂ mass M ≈ 0.6 M⊙ radius R ≈ 6400 km ≈ R⊕ mean ρ̄ ≈ 1×10⁶ g cm⁻³ gravity g ≈ 1×10⁸ cm s⁻² support law P ∝ ρ^5/3 (T-independent) a million times denser than water, yet cold degeneracy holds it rigid

Now turn the mass up. In an ordinary gas star, more mass means more heat and a bigger star; here the opposite is true. Adding mass strengthens gravity, which compresses the star to higher density until the now-denser electron gas, occupying higher momentum states, pushes back hard enough — and because P∝ρ5/3P\propto\rho^{5/3} rises steeply, balance is restored at a smaller radius. Keep going and the electrons near the top of the Fermi sea approach cc (become relativistic); the pressure law softens toward ρ4/3\rho^{4/3}, exactly the scaling gravity itself demands, and the self-correcting balance fails. The radius plunges to zero at one special mass. The second figure is that story as a single curve — the inverted mass–radius relation, ending at the Chandrasekhar wall.

Figure
The white-dwarf mass–radius relation for μₑ = 2 — radius falls as R ∝ M^−1/3 in the non-relativistic regime, then collapses to zero at the Chandrasekhar mass M_Ch ≈ 1.4 M⊙. A 0.6 M⊙ dwarf sits at R ≈ R⊕; heavier means smaller radius, and no equilibrium exists beyond the wall
The punchline: heavier white dwarfs are smaller — R ∝ M^−1/3 up to the wall at 1.4 M⊙ 0 0.5 1.0 1.5 2.0 0.2 0.6 1.0 1.4 mass M (M⊙) radius R (R⊕) 0.6 M⊙ → R ≈ R⊕ M_Ch ≈ 1.4 M⊙ R → 0, no equilibrium R ∝ M^−1/3 (non-relativistic branch)

Two features complete the physical picture. First, the highly degenerate electrons conduct heat so efficiently that the stellar interior remains nearly isothermal. At the same time, a thin, non-degenerate envelope acts as an insulating layer that regulates the escape of energy from the surface. Second, having exhausted all nuclear energy sources, the white dwarf is powered solely by the gradual release of the thermal energy stored in its ions; the degenerate electrons contribute little to the heat capacity because only those near the Fermi surface can be thermally excited. Consequently, the star evolves primarily by cooling rather than by structural change, maintaining an almost constant radius as it descends the white dwarf cooling sequence in the Hertzsprung–Russell diagram. As cooling proceeds, the interior eventually crystallizes, and over timescales far exceeding the present age of the Universe, the remnant is expected to evolve toward a cold, inert black dwarf.

The core equations

The engine of everything is the pressure of a cold, filled Fermi sea. Fill momentum space with two spin states (up & down) per phase-space cell up to pFp_F; the pressure is one-third of the momentum-times-velocity flux over that sphere. In the non-relativistic limit v=p/mev=p/m_e, the integral gives a clean power law in the electron number density.

Pdeg=(3π2)2/35ℏ2me ne5/3(non-rel.)(1)P_{\rm deg}=\frac{(3\pi^2)^{2/3}}{5}\frac{\hbar^2}{m_e}\,n_e^{5/3}\quad(\text{non-rel.}) \qquad (1)

where: PdegP_{\rm deg} — electron degeneracy pressure; ℏ\hbar — reduced Planck constant; mem_e — electron mass; nen_e — electron number density. Reading: pressure grows as the 5/35/3 power of electron number density and is independent of temperature — a cold, stiff support. Push the density up and the fastest electrons reach relativistic speeds. Their velocity saturates at cc, so a factor of p/mep/m_e in the pressure integral becomes a factor of cc, and the exponent drops by one-third — the seed of the mass limit.

Pdeg=(3π2)1/34 ℏc ne4/3(ultra-rel.)(2)P_{\rm deg}=\frac{(3\pi^2)^{1/3}}{4}\,\hbar c\,n_e^{4/3}\quad(\text{ultra-rel.}) \qquad (2)

where: cc — speed of light; nen_e — electron number density. Reading: the support softens to P∝ρ4/3P\propto\rho^{4/3} — the stiffness weakens just when gravity needs it most. To connect the electron gas to the mass density, count baryons per electron. Full ionization of carbon, oxygen or helium gives one electron per two nucleons, so μe≈2\mu_e\approx2; hydrogen would give μe≈1\mu_e\approx1. This one number sets the scale of the mass limit.

ne=ρμemu,μe=ρnemu≈2(3)n_e=\frac{\rho}{\mu_e m_u},\qquad \mu_e=\frac{\rho}{n_e m_u}\approx2 \qquad (3)

where: ρ\rho — mass density; mum_u — atomic mass unit; μe\mu_e — baryons per electron (≈2\approx2 for fully ionized C, O, He). Reading: a white dwarf carries two baryon masses per supporting electron. The star is a fluid in balance: the pressure gradient holds up the weight of every shell, and mass accumulates outward. Closed by a polytropic P(ρ)P(\rho), these two structure equations become the Lane–Emden equation.

dPdr=−G m(r) ρr2,dmdr=4πr2ρ(4)\frac{\mathrm{d}P}{\mathrm{d}r}=-\frac{G\,m(r)\,\rho}{r^2},\qquad \frac{\mathrm{d}m}{\mathrm{d}r}=4\pi r^2\rho \qquad (4)

where: m(r)m(r) — mass enclosed within radius rr; GG — gravitational constant. Reading: hydrostatic equilibrium plus mass continuity; with P=Kρ1+1/nP=K\rho^{1+1/n} the pair reduces to one dimensionless second-order ODE of index nn. Solving the n=3/2n=3/2 polytrope (the non-relativistic case) and restoring dimensions gives the hallmark inverse relation; the prefactor collects the Lane–Emden boundary constants and K1K_1.

R=(2.0) ℏ2G me (μemu)5/3 1M1/3  ∝  M−1/3(5)R=(2.0)\,\frac{\hbar^2}{G\,m_e\,(\mu_e m_u)^{5/3}}\,\frac{1}{M^{1/3}}\;\propto\;M^{-1/3} \qquad (5)

where: RR — stellar radius; MM — total mass; the numerical factor collects the n=3/2n=3/2 Lane–Emden constants. Reading: doubling the mass shrinks the radius by 2−1/3≈0.792^{-1/3}\approx0.79 — heavier white dwarfs are smaller, unlike main-sequence stars. At the relativistic n=3n=3 polytrope the mass formula loses all dependence on central density, leaving one critical mass built from constants alone — the Chandrasekhar mass.

MCh=ω3 03π2(ℏcG)3/21(μemu)2≈1.46(2μe)2M⊙(6)M_{\rm Ch}=\frac{\omega_3^{\,0}\sqrt{3\pi}}{2}\left(\frac{\hbar c}{G}\right)^{3/2}\frac{1}{(\mu_e m_u)^{2}}\approx1.46\left(\frac{2}{\mu_e}\right)^{2}M_\odot \qquad (6)

where: ω3 0≈2.018\omega_3^{\,0}\approx2.018 — structure constant of the n=3n=3 polytrope; mum_u — atomic mass unit; μe\mu_e — baryons per electron. Reading: only ℏ\hbar, cc, GG, mum_u and μe\mu_e appear — no temperature, no history; evaluating gives ≈1.4 M⊙\approx1.4\,M_\odot for μe=2\mu_e=2. Finally the observable life of the star: it shines as a blackbody of fixed radius, and its luminosity is fed by the leaking heat of its ions through the insulating envelope. The result is the Mestel cooling law.

L=4πR2σTeff4,L∝Tc7/2,tcool∝L−5/7(7)L=4\pi R^2\sigma T_{\rm eff}^4,\qquad L\propto T_c^{7/2},\qquad t_{\rm cool}\propto L^{-5/7} \qquad (7)

where: LL — luminosity; σ\sigma — Stefan–Boltzmann constant; TeffT_{\rm eff} — surface temperature; TcT_c — nearly uniform interior temperature. Reading: with RR frozen by degeneracy, luminosity tracks core temperature, and a white dwarf dims ever more slowly — the faintest are the oldest. The density at which the electrons turn relativistic and this softening begins, ρcrit=me3c3μemu/(3π2ℏ3)∼2×106 g cm−3\rho_{\rm crit}=m_e^3c^3\mu_e m_u/(3\pi^2\hbar^3)\sim2\times10^{6}\,\mathrm{g\,cm^{-3}}, sits right at typical central densities — which is why a real white dwarf is already mildly relativistic.

Derivation AThe Fermi pressure of a cold electron gas

The goal is to show, from nothing but the Pauli exclusion principle, that a zero-temperature electron gas has a real, calculable pressure depending only on its density — and to obtain the exponents 5/35/3 and 4/34/3 that control every result in this article. We fill momentum space up to a Fermi momentum, apply elementary kinetic theory for the pressure, and evaluate the integral in both the slow and fast limits.

Step 1. Count the states. In a volume VV the number of quantum cells up to momentum pp is 2⋅43πp3⋅V/h32\cdot\frac{4}{3}\pi p^3\cdot V/h^3, the factor 2 for spin. At T=0T=0 electrons fill every cell up to pFp_F, so ne=8π3h3pF3n_e=\frac{8\pi}{3h^3}p_F^3, which inverts to pF=ℏ(3π2ne)1/3p_F=\hbar(3\pi^2 n_e)^{1/3}. Why: degeneracy is just the ground state of a fermion gas, and pFp_F is set purely by how tightly the electrons are packed.

Step 2. Write the pressure. Kinetic theory gives P=13∫0pFp v 8πp2h3 dpP=\frac{1}{3}\int_0^{p_F} p\,v\,\frac{8\pi p^2}{h^3}\,\mathrm{d}p, one-third of the flux of momentum pp carried at velocity vv across a surface, summed over the filled sphere. Why: pressure is momentum transport, and the only inputs are the occupation (full up to pFp_F) and the relation between pp and vv.

Step 3. Do the non-relativistic case, v=p/mev=p/m_e. The integrand is p4/mep^4/m_e, and ∫0pFp4dp=pF5/5\int_0^{p_F}p^4\mathrm{d}p=p_F^5/5, giving P=8π15meh3pF5P=\frac{8\pi}{15 m_e h^3}p_F^5. Substituting pFp_F yields P=(3π2)2/35ℏ2mene5/3P=\frac{(3\pi^2)^{2/3}}{5}\frac{\hbar^2}{m_e}n_e^{5/3} — equation (1). Why: the ne5/3n_e^{5/3} law follows from pF∝ne1/3p_F\propto n_e^{1/3}; the exponent is pure phase-space counting.

Step 4. Do the ultra-relativistic case, v→cv\to c. Now the integrand is c p3c\,p^3, ∫0pFp3dp=pF4/4\int_0^{p_F}p^3\mathrm{d}p=p_F^4/4, and P=2πc3h3pF4=(3π2)1/34ℏc ne4/3P=\frac{2\pi c}{3h^3}p_F^4=\frac{(3\pi^2)^{1/3}}{4}\hbar c\,n_e^{4/3} — equation (2). Why: replacing v=p/mev=p/m_e by v=cv=c removes one power of pp, softening the exponent from 5/35/3 to 4/34/3 — the microscopic origin of the mass limit.

Go deeper · Derivation B · 1 minThe inverted mass–radius relation from pressure balance
Derivation BThe inverted mass–radius relation from pressure balance

The goal is to derive the counter-intuitive result that a white dwarf shrinks as it gains mass, R∝M−1/3R\propto M^{-1/3}, using only dimensional pressure balance — no full Lane–Emden solution required — and to recover the Earth-sized radius for a solar-mass-scale object. We estimate the central pressure gravity demands, set it equal to the non-relativistic degeneracy pressure the electrons supply, and solve for the equilibrium radius.

Step 1. Gravity's demand. From hydrostatic equilibrium, dP/dr∼−GMρ/R2\mathrm{d}P/\mathrm{d}r\sim -GM\rho/R^2 with ρ∼M/R3\rho\sim M/R^3, so the central pressure scales as Pc∼GM2/R4P_c\sim GM^2/R^4. Why: this is the weight of the column per unit area, rising steeply as the star is squeezed.

Step 2. The electrons' supply. Insert ne=ρ/(μemu)n_e=\rho/(\mu_e m_u) and ρ∼M/R3\rho\sim M/R^3 into equation (1): Pdeg∼ℏ2me(μemu)5/3M5/3R5P_{\rm deg}\sim\frac{\hbar^2}{m_e(\mu_e m_u)^{5/3}}\frac{M^{5/3}}{R^5}. Why: compression raises the degeneracy pressure too, but at a different power of RR — and that mismatch fixes equilibrium.

Step 3. Balance them. Setting GM2/R4∼ℏ2me(μemu)5/3M5/3/R5GM^2/R^4\sim \frac{\hbar^2}{m_e(\mu_e m_u)^{5/3}}M^{5/3}/R^5, one power of RR survives: R∼ℏ2Gme(μemu)5/3M−1/3R\sim\frac{\hbar^2}{G m_e(\mu_e m_u)^{5/3}}M^{-1/3}. Why: gravity scales as R−4R^{-4} but non-relativistic pressure as R−5R^{-5}, so pressure wins at small RR — there is always a stable equilibrium, and it moves inward as MM grows.

Step 4. Put numbers in. For M∼0.6 M⊙M\sim0.6\,M_\odot and μe=2\mu_e=2 the formula returns R∼10−2 R⊙∼7000 kmR\sim10^{-2}\,R_\odot\sim7000\,\mathrm{km}, an Earth-sized star. Why: the combination ℏ2/(Gmemu5/3)\hbar^2/(G m_e m_u^{5/3}) sets the natural length of degenerate matter, and it lands squarely at planetary radii.

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

The goal is to expose why a maximum mass must exist: when the electrons turn ultra-relativistic the pressure law softens to P∝ρ4/3P\propto\rho^{4/3}, the equilibrium radius drops out of the balance entirely, and only one special mass can satisfy it. We repeat the pressure balance with the relativistic equation of state, watch the radius cancel, and identify the surviving combination of constants as MChM_{\rm Ch}.

Step 1. Relativistic supply. Using equation (2) with ρ∼M/R3\rho\sim M/R^3, Pdeg∼ℏc ne4/3∼ℏc(μemu)4/3M4/3R4P_{\rm deg}\sim\hbar c\,n_e^{4/3}\sim\frac{\hbar c}{(\mu_e m_u)^{4/3}}\frac{M^{4/3}}{R^4}. Why: the softened exponent makes the pressure fall as R−4R^{-4} — the same power as gravity's demand.

Step 2. Balance against gravity. Setting GM2/R4∼ℏc(μemu)4/3M4/3/R4GM^2/R^4\sim\frac{\hbar c}{(\mu_e m_u)^{4/3}}M^{4/3}/R^4, both sides carry R−4R^{-4} and the radius cancels identically. Why: with pressure and gravity scaling the same way in RR, the star can no longer adjust its size to find balance.

Step 3. Solve for the mass. What remains is GM2∼ℏc(μemu)4/3M4/3GM^2\sim\frac{\hbar c}{(\mu_e m_u)^{4/3}}M^{4/3}, i.e. MCh∼(ℏcG)3/21(μemu)2M_{\rm Ch}\sim\left(\frac{\hbar c}{G}\right)^{3/2}\frac{1}{(\mu_e m_u)^2}. Why: equilibrium exists only at one mass; below it stiffness rescues the star, above it gravity always wins.

Step 4. The exact treatment. The full n=3n=3 Lane–Emden solution supplies the structure constant ω3 0=2.018\omega_3^{\,0}=2.018 and gives MCh=ω3 03π2(ℏcG)3/2(μemu)−2≈1.46 (2/μe)2 M⊙M_{\rm Ch}=\frac{\omega_3^{\,0}\sqrt{3\pi}}{2}\left(\frac{\hbar c}{G}\right)^{3/2}(\mu_e m_u)^{-2}\approx1.46\,(2/\mu_e)^2\,M_\odot. Why: the dimensional estimate already had the right constants; the polytrope only fixes the pure number, landing at ≈1.4 M⊙\approx1.4\,M_\odot. *Note: the n=3n=3 polytrope mass is independent of central density, and numerical evaluation gives MCh=1.456 M⊙M_{\rm Ch}=1.456\,M_\odot for μe=2\mu_e=2.*

Limiting cases

The degenerate equation of state and the structure it supports change qualitatively between regimes:

CaseResultMeaning
Non-relativistic degeneracy (low mass)P=K1ρ5/3,  K1=(3π2)2/3ℏ25me(μemu)5/3P=K_1\rho^{5/3},\ \ K_1=\dfrac{(3\pi^2)^{2/3}\hbar^2}{5m_e(\mu_e m_u)^{5/3}}polytrope index n=3/2n=3/2; stable equilibrium with R∝M−1/3R\propto M^{-1/3}; valid for M≲1 M⊙M\lesssim1\,M_\odot
Ultra-relativistic degeneracy (near limit)P=K2ρ4/3,  K2=(3π2)1/3ℏc4(μemu)4/3P=K_2\rho^{4/3},\ \ K_2=\dfrac{(3\pi^2)^{1/3}\hbar c}{4(\mu_e m_u)^{4/3}}polytrope index n=3n=3; mass fixed independent of radius, giving MChM_{\rm Ch}; no equilibrium for M>MChM>M_{\rm Ch}
Finite temperature (young hot dwarf)kBT/EF≪1k_BT/E_F\ll1thermal corrections to PP are of order (kBT/EF)2(k_BT/E_F)^2 and small; the star is well approximated as fully degenerate
GR and inverse-beta decay (extreme density)ρ≳109 g cm−3\rho\gtrsim10^{9}\,\mathrm{g\,cm^{-3}}general-relativistic corrections and electron capture onto nuclei destabilize the star just below the Newtonian limit, triggering collapse

Worked examples

1 · Mean density and the degeneracy check (fully worked)

A typical white dwarf has M≈0.6 M⊙M\approx0.6\,M_\odot packed into R≈R⊕R\approx R_\oplus. Find its mean density and confirm the electrons are strongly degenerate.

Step 1. Convert. M=0.6×1.989×1033 g=1.19×1033 gM=0.6\times1.989\times10^{33}\,\mathrm{g}=1.19\times10^{33}\,\mathrm{g}, and R=R⊕=6.37×108 cmR=R_\oplus=6.37\times10^{8}\,\mathrm{cm}.

Step 2. Mean density. ρˉ=3M4πR3=3(1.19×1033)4π(6.37×108)3\bar\rho=\dfrac{3M}{4\pi R^3}=\dfrac{3(1.19\times10^{33})}{4\pi(6.37\times10^{8})^3}. The denominator is 4π(2.58×1026)=3.24×1027 cm34\pi(2.58\times10^{26})=3.24\times10^{27}\,\mathrm{cm^3}, so ρˉ=1.1×106 g cm−3\bar\rho=\mathbf{1.1\times10^{6}\ g\,cm^{-3}} — about a million times denser than water.

Step 3. Electron density. With μe=2\mu_e=2, ne=ρˉ/(μemu)=1.1×106/(2×1.66×10−24)=3.3×1029 cm−3n_e=\bar\rho/(\mu_e m_u)=1.1\times10^6/(2\times1.66\times10^{-24})=\mathbf{3.3\times10^{29}\ cm^{-3}}.

Step 4. Degeneracy test. pF=ℏ(3π2ne)1/3p_F=\hbar(3\pi^2 n_e)^{1/3} gives EF≈pF2/2me∼0.2 MeVE_F\approx p_F^2/2m_e\sim0.2\ \mathrm{MeV}, whereas even at Tc∼107 KT_c\sim10^7\,\mathrm{K} the thermal energy is kBT∼1 keVk_BT\sim1\ \mathrm{keV}. Since kBT/EF∼10−2≪1k_BT/E_F\sim10^{-2}\ll1, the gas is strongly degenerate, and the zero-temperature equation of state is an excellent approximation.

2 · Is a real white dwarf relativistic? (fully worked)

Take ρc≈3×106 g cm−3\rho_c\approx3\times10^{6}\,\mathrm{g\,cm^{-3}}, a few times the mean. Compute pFp_F and compare pFcp_F c to mec2=0.511 MeVm_e c^2=0.511\,\mathrm{MeV}.

Step 1. Electron density at the centre: ne=ρc/(μemu)=3×106/(3.32×10−24)=9.0×1029 cm−3n_e=\rho_c/(\mu_e m_u)=3\times10^6/(3.32\times10^{-24})=9.0\times10^{29}\,\mathrm{cm^{-3}}.

Step 2. Fermi momentum: pF=ℏ(3π2ne)1/3p_F=\hbar(3\pi^2 n_e)^{1/3}. With 3π2ne=2.66×10313\pi^2 n_e=2.66\times10^{31}, its cube root is 3.0×1010 cm−13.0\times10^{10}\,\mathrm{cm^{-1}}, so pF=1.055×10−27×3.0×1010=3.15×10−17 g cm s−1p_F=1.055\times10^{-27}\times3.0\times10^{10}=3.15\times10^{-17}\,\mathrm{g\,cm\,s^{-1}}.

Step 3. Relativity parameter: pFc=3.15×10−17×3×1010=9.4×10−7 erg=0.59 MeVp_F c=3.15\times10^{-17}\times3\times10^{10}=9.4\times10^{-7}\,\mathrm{erg}=\mathbf{0.59\ MeV}. Against mec2=0.511 MeVm_e c^2=0.511\,\mathrm{MeV}, the ratio is x=pF/mec≈1.2x=p_F/m_e c\approx1.2.

Step 4. Interpretation. A typical white dwarf is mildly relativistic — its most energetic electrons already move near cc, so the equation of state sits partway between γ=5/3\gamma=5/3 and γ=4/3\gamma=4/3. This is why real white dwarfs feel the Chandrasekhar limit even at ordinary masses, and why the mass–radius curve steepens well before 1.4 M⊙1.4\,M_\odot.

3 · Evaluating the Chandrasekhar mass from constants (fully worked)

Compute MChM_{\rm Ch} for μe=2\mu_e=2 from equation (6) and confirm the famous ≈1.4 M⊙\approx1.4\,M_\odot.

Step 1. Core combination: ℏcG=(1.055×10−27)(3.0×1010)6.674×10−8=4.74×10−10 g2\dfrac{\hbar c}{G}=\dfrac{(1.055\times10^{-27})(3.0\times10^{10})}{6.674\times10^{-8}}=4.74\times10^{-10}\,\mathrm{g^2}, so (ℏc/G)3/2=1.03×10−14 g3(\hbar c/G)^{3/2}=1.03\times10^{-14}\,\mathrm{g^3}.

Step 2. Divide by the baryon-mass square: (μemu)2=(2×1.66×10−24)2=1.10×10−47 g2(\mu_e m_u)^2=(2\times1.66\times10^{-24})^2=1.10\times10^{-47}\,\mathrm{g^2}, giving 1.03×10−14/1.10×10−47=9.4×1032 g1.03\times10^{-14}/1.10\times10^{-47}=9.4\times10^{32}\,\mathrm{g}.

Step 3. Apply the structure prefactor ω3 03π2=2.018×3.0702=3.10\dfrac{\omega_3^{\,0}\sqrt{3\pi}}{2}=\dfrac{2.018\times3.070}{2}=3.10: MCh=3.10×9.4×1032=2.9×1033 gM_{\rm Ch}=3.10\times9.4\times10^{32}=2.9\times10^{33}\,\mathrm{g}.

Step 4. In solar units: MCh=2.9×1033/1.989×1033=1.46 M⊙M_{\rm Ch}=2.9\times10^{33}/1.989\times10^{33}=\mathbf{1.46\ M_\odot}. The compact catalogue form MCh=5.83 μe−2 M⊙M_{\rm Ch}=5.83\,\mu_e^{-2}\,M_\odot returns the same — the number that caps every white dwarf.

Physics of the portrait

The white dwarf is where quantum statistics first collided with gravity — a collision inseparable from the young man who worked it out at sea.

Origins · 1 minSubrahmanyan Chandrasekhar and how this was discovered
SC
Subrahmanyan Chandrasekhar
1910–1995 · Indian-American astrophysicist who derived the maximum mass of a white dwarf

In the summer of 1930, a nineteen-year-old Subrahmanyan Chandrasekhar boarded a ship from Madras to England to begin doctoral work at Cambridge. On the voyage he combined Fowler's 1926 idea that white dwarfs are held up by electron degeneracy with special relativity — and found that when the electrons are forced to relativistic speeds, the pressure softens from P∝ρ5/3P\propto\rho^{5/3} to P∝ρ4/3P\propto\rho^{4/3}, and the supporting mass can no longer grow without bound. There was a maximum: about 1.41.4 times the mass of the Sun, built from nothing but ℏ\hbar, cc, GG and the mass of a nucleon. The implication was seismic — a heavy enough stellar core has no stable endpoint as a white dwarf and must collapse to something denser still. Arthur Eddington, the towering figure of the field, rejected the result publicly, most famously at a 1935 Royal Astronomical Society meeting where he insisted there ought to be "a law of Nature to prevent a star from behaving in this absurd way." Eddington was wrong, but his prestige delayed acceptance for years. Chandrasekhar moved to Chicago and poured the calculation into his 1939 monograph An Introduction to the Study of Stellar Structure. The limiting mass now bears his name, underpins the theory of neutron stars, black holes and Type Ia supernovae, and won him the 1983 Nobel Prize — a vindication half a century in the making.

History
1862
— Alvan Clark resolves Sirius B, the first white dwarf ever identified
1926
— Fowler explains white-dwarf support by Fermi–Dirac electron degeneracy
1930
— Chandrasekhar derives the relativistic limiting mass en route to England
1935
— Chandrasekhar and Eddington clash publicly over the mass limit
1983
— Chandrasekhar receives the Nobel Prize for the theory of stellar structure

Going further

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 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. 2–4 · Kippenhahn, Weigert & Weiss, Stellar Structure and Evolution, Ch. 15, 35 · Mo, van den Bosch & White, Galaxy Formation and Evolution (stellar-remnant and cosmochronology context) · Dodelson, Modern Cosmology (standard candles and dark energy). Key numerical benchmarks: the Fermi-momentum integral yields P∝ne5/3P\propto n_e^{5/3} and ne4/3n_e^{4/3}; the n=3/2n=3/2 polytrope gives R∝M−1/3R\propto M^{-1/3} and an Earth-sized radius for 0.6 M⊙0.6\,M_\odot; the n=3n=3 mass is independent of central density, giving MCh=1.456 (2/μe)2 M⊙M_{\rm Ch}=1.456\,(2/\mu_e)^2\,M_\odot; the anchor densities are ρˉ≈1.1×106 g cm−3\bar\rho\approx1.1\times10^{6}\,\mathrm{g\,cm^{-3}} and ρcrit≈2×106 g cm−3\rho_{\rm crit}\approx2\times10^{6}\,\mathrm{g\,cm^{-3}}, with pFc≈0.59 MeVp_F c\approx0.59\,\mathrm{MeV} at ρc=3×106 g cm−3\rho_c=3\times10^{6}\,\mathrm{g\,cm^{-3}}; and the Mestel cooling law gives L∝Tc7/2L\propto T_c^{7/2}, tcool∝L−5/7t_{\rm cool}\propto L^{-5/7}.

Prerequisites. Electron Degeneracy Pressure · Equation of State of Dense Matter · Hydrostatic Equilibrium and the Lane–Emden Equation. Leads to. The Chandrasekhar Mass Limit · White-Dwarf Cooling and Crystallization · Neutron Stars · Type Ia Supernovae.

Ask the tutor: Ask me why the radius literally cancels out of the balance at γ=4/3\gamma=4/3, or to estimate the cooling age of a white dwarf from its luminosity.

Why is a white dwarf held up by electrons and not by the nuclei?

Both electrons and ions can be degenerate in principle, but degeneracy pressure scales inversely with particle mass: for a given momentum a lighter particle moves faster and exerts more pressure. Electrons are about 18001800 times lighter than nucleons, so at the same number density their degeneracy pressure utterly dominates. The heavy carbon and oxygen ions remain an almost classical gas that provides the star's mass and — through their thermal motion — its slowly draining heat reservoir. Support and mass thus come from different components: electrons hold the star up, ions weigh it down and store its warmth.

Does the Chandrasekhar mass depend on what the white dwarf is made of?

Only through μe\mu_e, the number of baryons per electron, and even that dependence is weak in practice. For any fully ionized mix of carbon, oxygen or helium there is one electron per two nucleons, so μe=2\mu_e=2 and MCh≈1.46 M⊙M_{\rm Ch}\approx1.46\,M_\odot. A pure-hydrogen composition would give μe=1\mu_e=1 and quadruple the limit — but hydrogen-dominated cores never survive to become white dwarfs, so that case is academic. What the limit does not depend on is temperature, rotation, or assembly history: it is built from ℏ\hbar, cc, GG and the nucleon mass alone, which is why the same number governs both Type Ia ignition and the white-dwarf–neutron-star boundary.

What is the difference between a white dwarf and a neutron star?

Both are degenerate compact remnants, but supported by different particles at different densities. A white dwarf is held up by electron degeneracy at ∼106 g cm−3\sim10^6\,\mathrm{g\,cm^{-3}} and is Earth-sized; a neutron star is held up by neutron degeneracy (stiffened by nuclear forces) at ∼1014 g cm−3\sim10^{14}\,\mathrm{g\,cm^{-3}} and is only about 10 km10\,\mathrm{km} across. The bridge is the Chandrasekhar limit: push a core past ≈1.4 M⊙\approx1.4\,M_\odot and electron pressure fails, electrons are captured onto protons by inverse beta decay, and the core collapses until neutron degeneracy halts it — governed then by the Tolman–Oppenheimer–Volkoff equation, not the Newtonian balance used here.

If white dwarfs have no fuel, why are some of them so hot?

Because they are born hot and cool slowly. A freshly exposed core inherits the temperature of the stellar interior that made it — surface temperatures can exceed 100,000 K100{,}000\,\mathrm{K} — with an enormous reservoir of thermal energy in the motion of its ions. With no nuclear source it can only radiate that heat away, but the degenerate electrons conduct so efficiently that the interior stays nearly isothermal while a thin envelope throttles the leak. The result is the Mestel law tcool∝L−5/7t_{\rm cool}\propto L^{-5/7}: the star fades ever more slowly, taking billions of years to grow faint. The hottest white dwarfs are simply the youngest.

What happens when a white dwarf is pushed over the Chandrasekhar mass?

One of two dramatic things, depending on composition. A carbon–oxygen white dwarf accreting toward the limit heats up until carbon ignites under degenerate conditions; because degeneracy pressure is temperature-independent, the burning cannot expand and cool the star, so it runs away and detonates the whole object as a Type Ia supernova, leaving nothing behind. An oxygen–neon–magnesium dwarf instead tends to lose its electron support to inverse beta decay before carbon ignites, collapsing to a neutron star ("accretion-induced collapse"). Either way, ≈1.4 M⊙\approx1.4\,M_\odot is a genuine end of the line — it cannot remain a white dwarf.

Will a white dwarf ever truly disappear?

Not on any timescale that has yet elapsed. As a white dwarf cools, its ions eventually arrange into a crystal lattice — the core solidifies, releasing latent heat that briefly slows the cooling, a delay now detected in Gaia data. It then fades toward a cold, dark "black dwarf" radiating essentially nothing. But reaching that state takes far longer than the current age of the Universe, so no black dwarf exists yet. Every white dwarf ever formed is still, in astronomical terms, glowing warm — the oldest are merely the dimmest, which is exactly what makes them such good clocks.

Before this

Where this leads