White Dwarf

A white dwarf is a compact, electron-degenerate stellar remnant that represents the final evolutionary endpoint for all low- and intermediate-mass main-sequence stars with initial birth masses between ~0.08 M_☉ and ~8 M_☉—accounting for approximately 97% of all stars in the Milky Way. Lacking any internal thermonuclear fusion source, a typical white dwarf packs a mass comparable to that of the Sun (M ≈ 0.5–0.7 M_☉, with a narrow observational peak near 0.60 M_☉) into a volume comparable to that of the Earth (R ≈ 0.008–0.02 R_☉, or 5,500–14,000 km). This extreme compression yields mean bulk densities of ρ ≈ 10⁵ to 10⁷ g cm⁻³ (roughly one metric ton per cubic centimeter) and surface gravities (log g ≈ 7.5–8.5 in cgs units) roughly 100,000 to 300,000 times stronger than Earth's surface gravity.

Because a white dwarf no longer generates nuclear energy to sustain thermal gas pressure, its hydrostatic stability against gravitational collapse is maintained by electron degeneracy pressure, a quantum-mechanical manifestation of the Pauli exclusion principle first applied to stellar interiors by Ralph H. Fowler in 1926. In 1930, Indian astrophysicist Subrahmanyan Chandrasekhar incorporated special relativity into the degenerate equation of state and demonstrated that an electron-degenerate star possesses a maximum stable mass—the Chandrasekhar limit of M_Ch ≈ 1.44 M_☉—beyond which degeneracy pressure fails, triggering either a thermonuclear Type Ia supernova or accretion-induced collapse into a neutron star.

---

Asymptotic Giant Branch Mass Loss and Core Stratification

When an intermediate-mass star (0.8 M_☉ ≲ M_init ≲ 8 M_☉) exhausts hydrogen in its core and leaves the main sequence, it ascends the red giant branch (RGB), ignites core helium fusion via the triple-alpha process (3 ⁴He → ¹²C + γ and ¹²C + ⁴He → ¹⁶O + γ), and subsequently ascends the asymptotic giant branch (AGB) with an inert carbon–oxygen (C/O) core surrounded by concentrically burning helium and hydrogen shells. During the thermally pulsing AGB phase, radiation-driven dust winds (Ṁ ≈ 10⁻⁴ M_☉ yr⁻¹) eject the star's loosely bound outer envelope as an expanding planetary nebula, exposing the hot (T_eff > 100,000 K), bare degenerate core at the center, which then enters the white dwarf cooling track in the lower-left quadrant of the Hertzsprung–Russell diagram.

Depending on progenitor mass and binary interactions, white dwarfs divide into three internal compositional classes:

1. Helium White Dwarfs (M ≲ 0.45 M_☉): Because single stars below ~0.8 M_☉ have main-sequence lifetimes longer than the current 13.8-Gyr age of the universe, observed low-mass helium (⁴He) white dwarfs form exclusively in close binary systems where Roche-lobe overflow or common-envelope stripping removes the red giant's hydrogen envelope before the core reaches the ~10⁸ K required for helium ignition.
2. Carbon–Oxygen (C/O) White Dwarfs (0.45 M_☉ ≲ M ≲ 1.05 M_☉): The vast majority of white dwarfs (descending from 0.8–8 M_☉ progenitors, including the future fate of the Sun) consist of a degenerate core of roughly equal parts carbon (¹²C) and oxygen (¹⁶O), encased in a thin residual helium layer (M_He ≈ 10⁻² M_*) and an even thinner outermost hydrogen layer (M_H ≈ 10⁻⁴ M_*).
3. Oxygen–Neon–Magnesium (ONeMg) White Dwarfs (1.05 M_☉ ≲ M ≲ 1.38 M_☉): Stars in the narrow super-AGB transition window (M_init ≈ 8–10.5 M_☉) reach core temperatures high enough (T_c ≈ 6 × 10⁸ K) to ignite off-center carbon fusion (¹²C + ¹²C), leaving a high-mass degenerate core composed of ¹⁶O, ²⁰Ne, and ²⁴Mg.

---

Fermi–Dirac Statistics and Non-Relativistic Polytropes

At densities of ρ ≈ 10⁶ g cm⁻³, the mean inter-particle spacing between ionized atomic nuclei and free electrons is on the order of Δx ≈ 10⁻¹⁰ cm (1 picometer), hundredfolds smaller than a Bohr atomic radius (a₀ = 5.29 × 10⁻⁹ cm). Under the Pauli exclusion principle, no two identical spin-1/2 fermions (electrons) can occupy the same quantum state inside a phase-space volume cell of h³. Consequently, even as thermal temperature T → 0, electrons are forced to fill momentum states up to a high Fermi momentum (p_F):

p_F = ℏ (3π² n_e)^(1/3) where n_e = (Z / A) (ρ / m_u) = ρ / (μ_e m_u)

Here n_e is the free electron number density and μ_e = A / Z ≈ 2.0 is the mean molecular weight per electron for helium, ¹²C, or ¹⁶O. In the non-relativistic regime (p_F ≪ m_e c, valid for ρ ≪ 10⁶ g cm⁻³ and M ≲ 0.7 M_☉), the electron translational kinetic energy is E = p² / (2 m_e), and the resulting integral over the momentum distribution yields the non-relativistic electron degeneracy pressure (P_e):

P_e = (1 / 5 m_e) (3π²)^(2/3) ℏ² n_e^(5/3) = K_NR (ρ / μ_e)^(5/3)

Remarkably, this equation of state (P ∝ ρ^(5/3), corresponding to a Lane–Emden polytrope of index n = 1.5) is independent of temperature T. Balancing the central degeneracy pressure P_c ∝ (M / R³)^(5/3) = M^(5/3) / R⁵ against the hydrostatic gravitational pressure P_grav ∝ G M² / R⁴ yields the counterintuitive white dwarf mass–radius relation:

R ∝ μ_e^(−5/3) M^(−1/3) and Volume V ∝ R³ ∝ M⁻¹

Thus, adding mass to a white dwarf causes it to shrink in radius and increase sharply in density: a 0.4 M_☉ white dwarf has a radius of ~11,000 km, a 0.6 M_☉ white dwarf (such as Sirius B as measured at ~0.98 M_☉, R ≈ 5,800 km, or 40 Eridani B at 0.57 M_☉, R ≈ 9,500 km) is roughly Earth-sized, and a 1.3 M_☉ white dwarf compresses to a radius of just ~3,000 km.

Relativistic Transition and the Chandrasekhar Limit

As a white dwarf's mass increases above ~0.8 M_☉ and its central density exceeds ρ_c ≈ 2 × 10⁶ g cm⁻³, the Fermi momentum p_F approaches and exceeds m_e c, meaning the degenerate electrons move at speeds approaching the speed of light c. In the ultra-relativistic limit (E → p c), the electron velocity can no longer increase proportionally with momentum (v ≤ c), softening the equation of state from P_e ∝ n_e^(5/3) to:

P_e,rel = (¼) (3π²)^(1/3) ℏ c n_e^(4/3) = K_UR (ρ / μ_e)^(4/3)

For an n = 3 polytrope (γ = 4/3), both the central degeneracy pressure (P_e ∝ (M / R³)^(4/3) = M^(4/3) / R⁴) and the inward gravitational pressure (P_grav ∝ G M² / R⁴) scale identically with radius as R⁻⁴. Radius cancels out of the equilibrium condition, yielding a unique maximum mass—the Chandrasekhar limit:

M_Ch = (ω₃⁰ √(3π) / 2) (ℏ c / G)^(3/2) (1 / μ_e m_u)² ≈ 5.836 / μ_e² M_☉ ≈ 1.44 M_☉ (for μ_e = 2)

When general relativistic instabilities and inverse beta decay (p + e⁻ → n + ν_e, neutronization) are included, the actual maximum stable mass for a non-rotating C/O white dwarf is ~1.38 M_☉, at which point escalating core density triggers runaway carbon fusion and a Type Ia supernova.

---

Gravitational Settling and Exoplanetary Debris Pollution

Because surface gravity on a white dwarf exceeds g ≈ 10⁸ cm s⁻², gravitational settling (sedimentation) acts within weeks to thousands of years to stratify elements by atomic weight: trace heavy elements sink into the interior while the lightest remaining element floats to the photosphere. Consequently, ~80% of white dwarfs display pure hydrogen atmospheres (spectral type DA), characterized by stark-broadened Balmer lines (Hα, Hβ, Hγ), while ~20% lost their residual hydrogen during a late AGB thermal pulse and exhibit helium-dominated atmospheres (DB, DO, DC).

Remarkably, high-resolution spectroscopy reveals that 25% to 50% of white dwarfs (DZ, DAZ, DBZ) display photospheric absorption lines of heavy rock-forming elements (Ca, Mg, Fe, Si, O). Because the gravitational sinking timescale of heavy metals in a DA atmosphere is only a few days to millennia—orders of magnitude shorter than the white dwarf's cooling age—these metals must be actively accreting today from the tidal disruption of rocky asteroids, comets, and minor planets scattered inward across the white dwarf's Roche limit (exemplified by the transiting disintegrating planetesimals and dusty gas disk orbiting WD 1145+017). White dwarf spectroscopy thus provides direct bulk geochemical assays of extrasolar planetary interiors.

Mestel's Cooling Law and Core Crystallization

In 1952, Leon Mestel established the analytical theory of white dwarf cooling. Because degenerate electrons have long mean free paths, thermal conduction keeps the interior degenerate core (>99% of the star's mass) nearly isothermal (T_core ≈ 10⁶–10⁸ K), while heat transport to space is throttled by a thin, non-degenerate outer radiative envelope (<1% of the mass) obeying Kramer's opacity (κ ∝ ρ T^(−3.5)). Connecting core thermal energy (U_th = (3/2) k_B T_core M / (A m_u)) to surface luminosity (L ∝ M T_core^(7/2)) yields Mestel's cooling law, where luminosity declines with cooling time t as L ∝ t^(−7/5) (T_eff ∝ t^(−7/20)).

As the core cools below T_core ≈ (2–10) × 10⁶ K, the electrostatic potential energy between positively charged ¹²C⁶⁺ and ¹⁶O⁸⁺ ions exceeds their thermal vibrational kinetic energy by a critical Coulomb coupling parameter Γ = (Z e)² / (a_i k_B T) ≈ 175. At this threshold, the liquid-like ionic plasma undergoes a first-order phase transition and crystallizes from the center outward into a solid body-centered cubic (bcc) lattice—releasing latent heat of fusion (~1 k_B T per ion) and gravitational sedimentation energy (as oxygen-rich crystals sink), which pauses cooling for 1 to 2 billion years and produces a distinct pile-up on the Gaia HR diagram. Because cooling to ambient cosmic background temperatures (~3 K) requires >10¹³–10¹⁵ years, no cooled "black dwarf" yet exists in our 13.8-billion-year-old universe.

---

Why does a more massive white dwarf have a smaller radius?

In a normal main-sequence star, adding mass increases internal thermal pressure and makes the star larger. In a white dwarf, gravity is balanced by non-relativistic electron degeneracy pressure (P ∝ ρ^(5/3)), which depends only on particle density and not on temperature. When mass is added, gravity increases (P_grav ∝ M² / R⁴), forcing the star to contract to a higher density (R ∝ M^(−1/3)) so that its tightly packed electrons acquire enough quantum Fermi momentum to balance the stronger gravitational pull.

How can a white dwarf contain heavy metals in its atmosphere if its gravity pulls them down in days?

The presence of calcium, iron, silicon, and magnesium in DZ and DAZ white dwarfs is direct evidence of surviving exoplanetary systems. When a star becomes a white dwarf, outer planets and asteroid belts survive at wider orbits and dynamically scatter rocky planetesimals inward. Passing inside the white dwarf's tidal disruption radius (~1 R_☉), these asteroids are shredded into a dusty debris disk that rains rock-forming elements onto the star's surface.

Will the Sun explode as a supernova when it becomes a white dwarf?

No. The Sun has an initial mass of 1.0 M_☉ and is an isolated single star. After exhausting core hydrogen and helium in roughly 5 to 7 billion years, it will shed ~45% of its mass in asymptotic giant branch winds and leave behind a stable ~0.54 M_☉ carbon–oxygen white dwarf—well below the 1.44 M_☉ Chandrasekhar limit required for a Type Ia supernova—which will slowly crystallize and cool over trillions of years.

---

1. Chandrasekhar, S. (1931). "The maximum mass of ideal white dwarfs." The Astrophysical Journal, 74, 81–82. https://doi.org/10.1086/143324
2. Mestel, L. (1952). "On the theory of white dwarf stars. I. The energy sources of white dwarfs." Monthly Notices of the Royal Astronomical Society, 112(6), 583–597. https://doi.org/10.1093/mnras/112.6.583
3. Koester, D., & Chanmugam, G. (1990). "Physics of white dwarf stars." Reports on Progress in Physics, 53(7), 837–915. https://doi.org/10.1088/0034-4885/53/7/001
4. Tremblay, P.-E., Fontaine, G., Fusillo, N. P. G., Dunlap, B. H., Gänsicke, B. T., Hollands, M. A., Hermes, J. J., Marsh, T. R., Cukanovaite, E., & Cunningham, T. (2019). "Core crystallization and pile-up in the cooling sequence of evolving white dwarfs." Nature, 565(7738), 202–205. https://doi.org/10.1038/s41586-018-0791-x