Electron Degeneracy Pressure

Electron degeneracy pressure is the force that holds up a dead star. When a Sun-like star exhausts its nuclear fuel, gravity is left unopposed by any burning, and the core contracts until a quantum-mechanical wall stops it — not a force between particles, but a consequence of counting. Electrons are spin-½ fermions, and the Pauli exclusion principle forbids any two of them from sharing the same quantum state. Pack them densely enough and the low-momentum states fill up completely, out to a sharp cutoff called the Fermi momentum pFp_F. Electrons denied a slow state are forced to carry large momenta, and the flux of that momentum is a pressure that survives even at absolute zero. The astrophysical payoff is a single, strange fact: this pressure depends only on the electron number density nen_e, not on the temperature. An ordinary gas cools and its pressure fails, so a cooling star must shrink; a degenerate gas holds its pressure as the temperature drops to nothing. That decoupling is what lets a white dwarf radiate for billions of years without collapsing — it does not need to be hot to resist gravity. The same physics with neutrons supports neutron stars. But the wall has a flaw. As density climbs, the fastest electrons approach the speed of light, and a relativistic gas cannot stiffen as steeply against compression: the equation of state softens from P∝ne5/3P\propto n_e^{5/3} to P∝ne4/3P\propto n_e^{4/3}. A self-gravitating body needs P∝ρ4/3P\propto\rho^{4/3} just to break even, so the softening removes the star's margin of safety. Above a critical mass — the Chandrasekhar limit, about 1.4 M⊙1.4\,M_\odot — no degenerate equilibrium exists at all, and the remnant must collapse.

Key takeaways
  • Electron degeneracy pressure is the non-thermal pressure of a dense electron gas whose momentum states are filled to the Fermi momentum pF=ℏ(3π2ne)1/3p_F=\hbar(3\pi^2 n_e)^{1/3} because the Pauli principle forbids double occupancy.
  • It survives to T=0T=0 and depends only on the electron number density nen_e, not the temperature — unlike ordinary gas pressure P=nekBTP=n_e k_B T.
  • Non-relativistic: P=(3π2)2/35ℏ2mene5/3P=\frac{(3\pi^2)^{2/3}}{5}\frac{\hbar^2}{m_e}n_e^{5/3} (polytrope γ=5/3\gamma=5/3); ultra-relativistic: P=(3π2)1/34ℏc ne4/3P=\frac{(3\pi^2)^{1/3}}{4}\hbar c\,n_e^{4/3} (γ=4/3\gamma=4/3).
  • The relativistic softening to γ=4/3\gamma=4/3 is exactly why a maximum mass exists — the Chandrasekhar limit MCh≈1.44 (2/μe)2 M⊙M_\mathrm{Ch}\approx1.44\,(2/\mu_e)^2\,M_\odot.
  • A white dwarf of central density 106 g cm−310^6\ \mathrm{g\,cm^{-3}} has a Fermi temperature TF∼2×109 KT_F\sim2\times10^9\ \mathrm{K}, far above its interior temperature, so its electrons are strongly degenerate though the star glows white-hot.
  • Key idea: filling momentum space to pFp_F forces electrons into high-momentum states independent of temperature; the resulting momentum flux is a pressure P∝ne5/3P\propto n_e^{5/3} (non-relativistic) softening to ne4/3n_e^{4/3} (ultra-relativistic) that supports a cold, dense star up to a maximum mass.
The core insight

Filling momentum space to the Fermi momentum forces electrons into high-momentum states independent of temperature; the resulting momentum flux is a pressure P∝ne5/3P\propto n_e^{5/3} (non-relativistic) or P∝ne4/3P\propto n_e^{4/3} (ultra-relativistic) that supports a cold, dense star — and the relativistic softening caps its mass.

Where this lives today

The degenerate equation of state is now a precision tool, not a curiosity. Type Ia supernovae — the standardizable candles that revealed the accelerating Universe — are believed to be white dwarfs driven to or near the Chandrasekhar mass by accretion or merger, so the value 1.44 M⊙1.44\,M_\odot is quite literally built into the cosmic distance ladder. ESA's Gaia mission has measured parallaxes and fluxes for hundreds of thousands of white dwarfs, mapping the γ=5/3\gamma=5/3 mass–radius relation to per-cent accuracy and even resolving the pile-up of stars whose cooling has stalled as they crystallize into a quantum Coulomb lattice. At the extreme end, the same Fermi-gas logic with me→mnm_e\to m_n underlies neutron stars, where NICER's X-ray timing and LIGO–Virgo's tidal-deformability measurement in GW170817 now constrain the dense-matter equation of state directly. And in the laboratory it is not exotic at all: degeneracy pressure is what makes ordinary metals nearly incompressible and what stabilizes ultracold trapped Fermi gases.

The physical picture

Picture the electrons not in ordinary space but in momentum space, where each quantum state is a tiny cell of volume h3h^3 holding exactly two electrons — spin up and spin down. At absolute zero the electrons fill the cells from the center outward, taking the lowest-momentum states first. But cheap states run out. With enough electrons in a small box you must keep filling until the outermost occupied cells sit at a large radius pFp_F, the Fermi momentum. Those outermost electrons move fast even though the gas is stone cold, and their ceaseless drumming on the walls is the degeneracy pressure. Squeeze the box and the same electrons must fit into a smaller volume of real space, so by the uncertainty principle they spread over a larger volume of momentum space — pFp_F rises, and the pressure climbs steeply. That stiffness is what resists gravity.

Figure
Electrons fill momentum space to the Fermi radius p_F. Ordinary matter (left) fills only a small sphere; white-dwarf matter (right) at n_e ≈ 3×10³⁵ m⁻³ fills out to p_F ≈ 0.80 m_e c, close to the dashed relativistic boundary p = m_e c where the gas begins to soften
Momentum space fills to p_F — compress the gas and p_F climbs toward the relativistic wall p_F p_F p = m_e c compress n_e ↑ ⇒ p_F ↑ ordinary matter n_e ≈ 10²⁹ m⁻³, p_F ≪ m_e c white-dwarf core n_e ≈ 3×10³⁵ m⁻³, p_F ≈ 0.80 m_e c

The contrast with an ordinary gas is the whole point. A thermal gas obeys P=nekBTP=n_e k_B T, so as it cools its pressure fails and it contracts; a degenerate gas sets its pressure by density alone, so a white dwarf can radiate away its heat over the age of the Universe while holding its shape and its own weight. But there is a catch built into relativity. As density rises, pFp_F eventually exceeds mecm_e c, and the fastest electrons can no longer speed up in proportion to their momentum — their velocity saturates at cc. The pressure then grows more slowly with density, and the equation of state softens from a stiff γ=5/3\gamma=5/3 power law to a fragile γ=4/3\gamma=4/3. The second figure traces this on the pressure–density plane: two straight lines on a log–log plot, hinged where the electrons turn relativistic.

Figure
The degenerate equation of state on a log–log plane: stiff P ∝ ρ^(5/3) below the transition, softening to P ∝ ρ^(4/3) above it. The hinge sits where p_F = m_e c, at ρ ≈ 2×10⁹ kg m⁻³ (μ_e = 2). The softer 4/3 slope is exactly the slope gravity demands, which is why a maximum mass appears
The punchline: the equation of state softens from γ = 5/3 to γ = 4/3 at p_F = m_e c 10¹⁶ 10²⁰ 10²⁴ 10²⁸ 10⁶ 10⁹ 10¹² 10¹⁵ mass density ρ (kg m⁻³) pressure P (Pa) P ∝ ρ^(5/3) non-relativistic, stiff P ∝ ρ^(4/3) ultra-relativistic, soft p_F = m_e c ρ ≈ 2×10⁹

The core equations

Everything begins with counting states. At T=0T=0 the electrons occupy every momentum cell out to pFp_F; equating that count to the number density gives the Fermi momentum, the single quantity that controls the whole equation of state.

ne=8π3h3 pF3  ⟺  pF=ℏ (3π2ne)1/3(1)n_e=\frac{8\pi}{3h^3}\,p_F^3 \;\Longleftrightarrow\; p_F=\hbar\,(3\pi^2 n_e)^{1/3} \qquad (1)

where: nen_e — electron number density; pFp_F — Fermi momentum, the sharp cutoff of filled states; hh — Planck constant; ℏ=h/2π\hbar=h/2\pi. Reading: the Fermi surface swells as the cube root of density — there is no avoiding ever-faster states. Pressure is a flux of momentum. For an isotropic gas the pressure is one-third of the momentum carried across a surface, integrated over every occupied state — a formula that holds whether the electrons are slow or relativistic, because it never assumes a velocity law.

P=13∫0pFp v(p) n(p) dp,n(p) dp=8πh3p2 dp(2)P=\frac{1}{3}\int_0^{p_F} p\,v(p)\,n(p)\,\mathrm{d}p,\qquad n(p)\,\mathrm{d}p=\frac{8\pi}{h^3}p^2\,\mathrm{d}p \qquad (2)

where: PP — pressure; v(p)v(p) — speed of an electron of momentum pp; n(p) dpn(p)\,\mathrm{d}p — number of electrons per unit volume with momentum in [p,p+dp][p,p+\mathrm{d}p], from the same phase-space counting as (1). Reading: to get a number you need only the velocity–momentum relation, and that is where relativity enters. The velocity law is the hinge of the whole subject. At low momentum it is Newtonian, v=p/mev=p/m_e; near and beyond pF∼mecp_F\sim m_e c it saturates at the speed of light.

v(p)=pc2p2c2+me2c4  → p≪mec   pme,→ p≫mec   c(3)v(p)=\frac{p c^2}{\sqrt{p^2c^2+m_e^2c^4}}\;\xrightarrow{\,p\ll m_e c\,}\;\frac{p}{m_e},\qquad \xrightarrow{\,p\gg m_e c\,}\;c \qquad (3)

where: cc — speed of light; mem_e — electron mass. Reading: the crossover at pF≈mecp_F\approx m_e c turns the stiff non-relativistic gas into the soft relativistic one — the physical origin of the mass limit. Carrying out the integral with the two limiting velocity laws gives the two workhorse equations of state. Below, the non-relativistic law that supports ordinary white dwarfs.

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

where: ℏ\hbar — reduced Planck constant; mem_e — electron mass. Reading: the exponent 5/35/3 makes the gas stiff — double the density and the pressure rises by 25/3≈3.172^{5/3}\approx3.17. At very high density the velocity saturates and the same integral softens by one power.

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

where: cc — speed of light. Reading: the electron mass has dropped out entirely — the ultra-relativistic pressure knows only ℏ\hbar, cc, and density, and this scale-free form is what fixes a unique mass. The stars are made of nuclei but supported by electrons, so we need the bookkeeping that links mass density to electron density, and the polytropic packaging that feeds directly into stellar-structure theory.

ρ=μemune,P=Kργ,γ=53 or 43(6)\rho=\mu_e m_u n_e,\qquad P=K\rho^{\gamma},\quad \gamma=\tfrac{5}{3}\ \text{or}\ \tfrac{4}{3} \qquad (6)

where: ρ\rho — total mass density; μe\mu_e — mean molecular weight per electron (μe≈2\mu_e\approx2 for helium, carbon, and oxygen); mum_u — atomic mass unit; KK — polytropic constant; γ\gamma — polytropic index. Reading: nucleons carry the mass, electrons carry the pressure, and μe\mu_e is the exchange rate between them. Finally, the criterion that decides whether a gas is degenerate at all: compare the Fermi energy to the thermal energy.

EF=pF2c2+me2c4−mec2  →NR  ℏ2(3π2ne)2/32me,TF=EFkB(7)E_F=\sqrt{p_F^2c^2+m_e^2c^4}-m_e c^2\;\xrightarrow{\text{NR}}\;\frac{\hbar^2(3\pi^2 n_e)^{2/3}}{2m_e},\qquad T_F=\frac{E_F}{k_B} \qquad (7)

where: EFE_F — Fermi (kinetic) energy of the topmost filled state; TFT_F — degeneracy temperature; kBk_B — Boltzmann constant. Reading: when T≪TFT\ll T_F the thermal energy cannot lift electrons across the Fermi surface and the gas is degenerate.

Derivation ANon-relativistic degeneracy pressure and the five-thirds law

The goal is to turn the abstract pressure integral of equation (2) into the concrete power law P∝ne5/3P\propto n_e^{5/3} that supports ordinary white dwarfs, using nothing but phase-space counting, the Newtonian velocity v=p/mev=p/m_e, and the Fermi momentum from equation (1). We fill the Fermi sphere at zero temperature, insert the low-momentum velocity law, do a one-line integral, and eliminate pFp_F in favor of the observable density — and the temperature never appears, which is the whole miracle.

Step 1. Substitute the momentum distribution n(p)=(8π/h3)p2n(p)=(8\pi/h^3)p^2 into the pressure integral, giving P=13∫0pFp v(p) 8πh3p2 dpP=\frac{1}{3}\int_0^{p_F} p\,v(p)\,\frac{8\pi}{h^3}p^2\,\mathrm{d}p. Why: this collapses the problem to a single integral over the filled sphere; all that remains is to choose the velocity law.

Step 2. In the non-relativistic regime p≤pF≪mecp\le p_F\ll m_e c, so v(p)=p/mev(p)=p/m_e and the integrand becomes p4/mep^4/m_e: P=8π3h3me∫0pFp4 dpP=\frac{8\pi}{3h^3 m_e}\int_0^{p_F}p^4\,\mathrm{d}p. Why: the fastest electron is still slow compared to light, so the simple Newtonian relation is accurate to better than a percent.

Step 3. The integral gives pF5/5p_F^5/5, so P=8π15 h3mepF5P=\frac{8\pi}{15\,h^3 m_e}p_F^5. Why: a pure power of pFp_F — the pressure is set entirely by where the Fermi surface sits, and nothing else.

Step 4. Eliminate pFp_F using pF=ℏ(3π2ne)1/3p_F=\hbar(3\pi^2 n_e)^{1/3} from equation (1), and write h3=(2πℏ)3=8π3ℏ3h^3=(2\pi\hbar)^3=8\pi^3\hbar^3. Collecting the constants, P=(3π2)2/35ℏ2mene5/3P=\frac{(3\pi^2)^{2/3}}{5}\frac{\hbar^2}{m_e}n_e^{5/3}. Why: this trades the invisible pFp_F for the measurable density and exposes the 5/35/3 exponent that makes the gas stiff.

Step 5. Notice what is absent: there is no TT. The pressure is a pure function of density, so a star built on this law can cool to absolute zero without losing an ounce of support. Why: this temperature-independence is the physical foundation of white-dwarf stability and of white-dwarf cooling. *Note: dimensionally [ℏ2/me][ne5/3][\hbar^2/m_e][n_e^{5/3}] is an energy density, i.e. a pressure. Numerically at ne=3.01×1035 m−3n_e=3.01\times10^{35}\ \mathrm{m^{-3}} this gives PNR=3.2×1021 PaP_\mathrm{NR}=3.2\times10^{21}\ \mathrm{Pa}.*

Go deeper · Derivation B · 1 minUltra-relativistic softening and the origin of a maximum mass
Derivation BUltra-relativistic softening and the origin of a maximum mass

The aim is to show why a maximum mass must exist — not to compute its precise value, but to expose the mechanism. We repeat the pressure integral with the relativistic velocity saturated at cc, obtain the softer γ=4/3\gamma=4/3 law in which the electron mass vanishes, and then set that pressure against the gravitational demand of a self-gravitating body. The radius cancels, and a single critical mass is left standing, independent of how hard you squeeze.

Step 1. In the ultra-relativistic limit v→cv\to c, the integrand in equation (2) becomes c p3c\,p^3: P=8πc3h3∫0pFp3 dp=2πc3h3pF4P=\frac{8\pi c}{3h^3}\int_0^{p_F}p^3\,\mathrm{d}p=\frac{2\pi c}{3h^3}p_F^4. Why: replacing v=p/mev=p/m_e by v=cv=c removes one power of pp and one factor of 1/me1/m_e — the softening and the disappearance of the electron mass happen in the same stroke.

Step 2. Insert pF=ℏ(3π2ne)1/3p_F=\hbar(3\pi^2 n_e)^{1/3} to get PUR=(3π2)1/34ℏc ne4/3∝ρ4/3P_\mathrm{UR}=\frac{(3\pi^2)^{1/3}}{4}\hbar c\,n_e^{4/3}\propto\rho^{4/3}. Why: the exponent has dropped from 5/35/3 to 4/34/3, and the constant KK now contains only ℏ\hbar and cc — a warning sign, because 4/34/3 is precisely the gravitational threshold.

Step 3. Estimate the gravitational demand. Hydrostatic equilibrium gives a central pressure Pc∼GM2/R4P_c\sim GM^2/R^4, while the mean density is ρ∼M/R3\rho\sim M/R^3, so the pressure needed to hold up the star scales as GM2/R4GM^2/R^4. Why: this is the load the electron gas must carry; comparing its supply to this demand decides whether equilibrium is possible.

Step 4. Set supply equal to demand: K(M/R3)4/3∼GM2/R4K(M/R^3)^{4/3}\sim GM^2/R^4. The left side scales as M4/3/R4M^{4/3}/R^4 and the right as M2/R4M^2/R^4 — the radius R4R^4 cancels on both sides. Why: because both pressure and gravity scale as 1/R41/R^4 for a γ=4/3\gamma=4/3 gas, the star has no lever to adjust; changing RR changes nothing.

Step 5. What remains is a condition on mass alone: M4/3∼M2M^{4/3}\sim M^2 can hold for only one value, Mmax⁡∼(ℏcG)3/21(μemu)2M_{\max}\sim\left(\frac{\hbar c}{G}\right)^{3/2}\frac{1}{(\mu_e m_u)^2}. Why: with the radius powerless, exactly one critical mass balances the two effects — the Chandrasekhar mass. Below it a stiffening correction saves the star; above it, gravity always wins. *Note: the scaling (ℏc/G)3/2/(μemu)2(\hbar c/G)^{3/2}/(\mu_e m_u)^2 evaluates to 1.9 M⊙1.9\,M_\odot for μe=2\mu_e=2; the exact Lane–Emden n=3n=3 coefficient 2.018/(4π)2.018/(4\pi) sharpens this to 1.44 M⊙1.44\,M_\odot.*

Go deeper · Derivation C · 1 minThe degeneracy temperature and why hot stars can be cold gases
Derivation CThe degeneracy temperature and why hot stars can be cold gases

The purpose here is to make quantitative the loose phrase "cold, dense gas" and to resolve an apparent paradox: white dwarfs glow at tens of thousands of kelvin, yet we call their electrons degenerate. The resolution is that "cold" means cold relative to the Fermi temperature, and for stellar densities that temperature is astronomically high. We compute EFE_F, define TFT_F, and read off the criterion.

Step 1. Take the non-relativistic Fermi energy, EF=pF2/2me=ℏ2(3π2ne)2/3/2meE_F=p_F^2/2m_e=\hbar^2(3\pi^2 n_e)^{2/3}/2m_e, the kinetic energy of the most energetic filled state. Why: this is the energy an electron would need to be promoted just above the Fermi surface — the natural yardstick for thermal disturbance.

Step 2. Define the degeneracy temperature TF=EF/kBT_F=E_F/k_B. The gas is degenerate when kBT≪EFk_B T\ll E_F, i.e. T≪TFT\ll T_F; then thermal agitation can shuffle only the thin sliver of electrons within kBTk_B T of the Fermi surface, leaving the deep sea untouched. Why: exclusion, not temperature, then sets the pressure, and finite-TT corrections enter only at order (kBT/EF)2(k_B T/E_F)^2 through the Sommerfeld expansion.

Step 3. For a white-dwarf core at ne≈3×1035 m−3n_e\approx3\times10^{35}\ \mathrm{m^{-3}} (ρ≈106 g cm−3\rho\approx10^6\ \mathrm{g\,cm^{-3}}), EF≈160 keVE_F\approx160\ \mathrm{keV} and TF≈1.9×109 KT_F\approx1.9\times10^9\ \mathrm{K}, while the interior temperature is ∼107 K\sim10^7\ \mathrm{K}. Why: with T/TF∼10−2T/T_F\sim10^{-2} the electrons are strongly degenerate even though the star is hot by any human standard — the paradox dissolves. *Note: TF=1.90×109 KT_F=1.90\times10^9\ \mathrm{K} at ρ=106 g cm−3\rho=10^6\ \mathrm{g\,cm^{-3}}, and the Sommerfeld corrections are of order (kBT/EF)2∼3×10−5(k_B T/E_F)^2\sim3\times10^{-5}, wholly negligible.*

Limiting cases

The degenerate equation of state has two limiting forms bracketing the transition at pF∼mecp_F\sim m_e c, plus the classical and finite-temperature limits that mark its edges.

CaseResultMeaning
Non-relativistic degeneracyP=(3π2)2/35ℏ2mene5/3P=\frac{(3\pi^2)^{2/3}}{5}\frac{\hbar^2}{m_e}n_e^{5/3}Valid for pF≪mecp_F\ll m_e c (moderate density). Polytrope γ=5/3\gamma=5/3; stiff — pressure rises fast under compression, giving stable white dwarfs.
Ultra-relativistic degeneracyP=(3π2)1/34ℏc ne4/3P=\frac{(3\pi^2)^{1/3}}{4}\hbar c\,n_e^{4/3}Valid for pF≫mecp_F\gg m_e c (very high density). Polytrope γ=4/3\gamma=4/3; soft — the electron mass drops out and a maximum mass appears.
Classical (non-degenerate) limitP=nekBTP=n_e k_B TWhen kBT≫EFk_B T\gg E_F the Pauli principle is inoperative and ordinary ideal-gas pressure returns; degeneracy is negligible.
Finite-temperature correctionP=P0 ⁣[1+5π212 ⁣(kBTEF)2+⋯ ]P=P_0\!\left[1+\frac{5\pi^2}{12}\!\left(\frac{k_B T}{E_F}\right)^2+\cdots\right]Sommerfeld expansion: the leading thermal correction is tiny when T≪TFT\ll T_F, so real white dwarfs are very nearly fully degenerate.

Worked examples

1 · The central conditions of a white dwarf (fully worked)

A carbon–oxygen white dwarf has central mass density ρ=106 g cm−3=109 kg m−3\rho=10^6\ \mathrm{g\,cm^{-3}}=10^9\ \mathrm{kg\,m^{-3}} and μe=2\mu_e=2. Find the electron number density, the Fermi momentum, whether the gas is relativistic, the degeneracy pressure, and the Fermi temperature. Its interior temperature is ∼107 K\sim10^7\ \mathrm{K} — is it degenerate?

Step 1. Electron density from equation (6): ne=ρ/(μemu)=109/(2×1.6605×10−27)=3.01×1035 m−3n_e=\rho/(\mu_e m_u)=10^9/(2\times1.6605\times10^{-27})=\mathbf{3.01\times10^{35}\ m^{-3}}.

Step 2. Fermi momentum from equation (1): pF=ℏ(3π2ne)1/3=(1.0546×10−34)(3π2×3.01×1035)1/3=2.19×10−22 kg m s−1p_F=\hbar(3\pi^2 n_e)^{1/3}=(1.0546\times10^{-34})(3\pi^2\times3.01\times10^{35})^{1/3}=\mathbf{2.19\times10^{-22}\ kg\,m\,s^{-1}}.

Step 3. Relativistic? Compare to mec=9.109×10−31×2.998×108=2.73×10−22m_e c=9.109\times10^{-31}\times2.998\times10^8=2.73\times10^{-22}. The ratio is x=pF/mec=0.80x=p_F/m_e c=\mathbf{0.80} — the fastest electrons are mildly relativistic, sitting right at the hinge of the two figures, so the true pressure lies between the NR and UR laws.

Step 4. Non-relativistic estimate from equation (4): PNR=(3π2)2/35ℏ2mene5/3=3.2×1021 PaP_\mathrm{NR}=\frac{(3\pi^2)^{2/3}}{5}\frac{\hbar^2}{m_e}n_e^{5/3}=\mathbf{3.2\times10^{21}\ Pa}. The ultra-relativistic law (5) gives PUR=4.9×1021 PaP_\mathrm{UR}=4.9\times10^{21}\ \mathrm{Pa}; the exact pressure at x=0.80x=0.80 falls between, around 3×1021 Pa3\times10^{21}\ \mathrm{Pa} — some 10510^5 times the pressure at the center of the Sun.

Step 5. Fermi temperature from equation (7): EF≈160 keVE_F\approx160\ \mathrm{keV}, so TF=EF/kB=1.9×109 KT_F=E_F/k_B=\mathbf{1.9\times10^{9}\ K}. Since the interior temperature ∼107 K\sim10^7\ \mathrm{K} gives T/TF∼0.005≪1T/T_F\sim0.005\ll1, the electrons are strongly degenerate — the star can cool by two more orders of magnitude in temperature without its support weakening at all.

2 · The Chandrasekhar mass from dimensional balance

Estimate the maximum mass of a white dwarf (μe=2\mu_e=2) from the ultra-relativistic equation of state, and compare to the exact value.

Step 1. Set the degenerate supply against the gravitational demand, as in Derivation 2: Kρ4/3∼GM2/R4K\rho^{4/3}\sim GM^2/R^4 with ρ∼M/R3\rho\sim M/R^3. Because both sides scale as R−4R^{-4}, the radius cancels and only mass survives: Mmax⁡∼(ℏc/G)3/2/(μemu)2M_{\max}\sim(\hbar c/G)^{3/2}/(\mu_e m_u)^2.

Step 2. Put in numbers: ℏc/G=4.74×10−16 kg2\hbar c/G=4.74\times10^{-16}\ \mathrm{kg^2}, so (ℏc/G)3/2=1.03×10−23 kg3(\hbar c/G)^{3/2}=1.03\times10^{-23}\ \mathrm{kg^3}. Dividing by (μemu)2=(2×1.66×10−27)2=1.10×10−53 kg2(\mu_e m_u)^2=(2\times1.66\times10^{-27})^2=1.10\times10^{-53}\ \mathrm{kg^2} gives Mmax⁡≈9.4×1029 kg≈1.9 M⊙M_{\max}\approx\mathbf{9.4\times10^{29}\ kg}\approx1.9\,M_\odot.

Step 3. The order-of-magnitude balance overshoots because it ignores the run of density through the star. Solving the γ=4/3\gamma=4/3 Lane–Emden (n=3n=3) equation properly supplies the dimensionless coefficient 2.018/(4π)2.018/(4\pi), which sharpens the estimate to MCh=1.457 (2/μe)2 M⊙=1.44 M⊙M_\mathrm{Ch}=1.457\,(2/\mu_e)^2\,M_\odot=\mathbf{1.44\ M_\odot} for μe=2\mu_e=2 — the observed upper mass of white dwarfs and the fuse that lights Type Ia supernovae. The dimensional estimate lands within 30 % of the exact limit, and, crucially, predicts that the answer is a pure mass independent of radius.

Physics of the portrait

Electron degeneracy pressure is where quantum mechanics reached out and touched the sky. Within four years of Fermi and Dirac writing down the statistics of identical fermions, a teenager on a steamship worked out that those statistics set an absolute ceiling on the mass of a dead star — a claim so unwelcome it triggered one of the most famous feuds in astrophysics.

Origins · 1 minChandrasekhar and how this was discovered
C
Chandrasekhar
1910–1995 · the relativistic degenerate star and its maximum mass

Subrahmanyan Chandrasekhar built the theory of the degenerate star before he was twenty. Fermi–Dirac statistics had appeared only in 1926, and that same year Ralph Fowler used them to explain why white dwarfs — puzzlingly dense, since Sirius B packs a solar mass into an Earth-sized volume — do not collapse: their electrons are degenerate, held apart by exclusion rather than heat. Fowler's gas was non-relativistic and seemed to allow a white dwarf of any mass. In 1930, on the sea voyage from Madras to Cambridge to begin his doctorate, Chandrasekhar noticed what Fowler had missed: at the densities of a massive white dwarf the electrons must be relativistic, the equation of state softens to γ=4/3\gamma=4/3, and the radius shrinks to zero at a finite mass of about 1.4 M⊙1.4\,M_\odot, above which no equilibrium exists. Lev Landau reached a similar conclusion in 1932. Arthur Eddington, the towering figure of British astrophysics, found the result absurd and attacked it publicly at the 1935 Royal Astronomical Society meeting, declaring there should 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. The collapse he could not stomach is now understood to produce neutron stars, black holes, and — when a white dwarf is pushed over the edge — the Type Ia supernovae used to measure the expansion of the Universe. Chandrasekhar shared the 1983 Nobel Prize for this work, half a century after a nineteen-year-old first sketched it in a ship's notebook.

History
1844
— Bessel infers an unseen companion of Sirius from its proper-motion wobble
1926
— Fowler applies the new Fermi–Dirac statistics to the dense matter of white dwarfs
1930
— Chandrasekhar derives a maximum mass for a relativistic degenerate star at age nineteen
1935
— Eddington disputes the limiting mass at the Royal Astronomical Society
1983
— Chandrasekhar receives the Nobel Prize in Physics for the theory of stellar structure

Going further

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., Ch. 16 (Degenerate Remnants of Stars) · Shapiro & Teukolsky, Black Holes, White Dwarfs, and Neutron Stars, Ch. 2–3 · Chandrasekhar, An Introduction to the Study of Stellar Structure (1939) · Kippenhahn, Weigert & Weiss, Stellar Structure and Evolution, Ch. 15 · Binney & Tremaine, Galactic Dynamics, 2nd ed. (Fermi-gas appendix). The key results follow from phase-space counting and the pressure integral: ne=8πpF3/3h3n_e=8\pi p_F^3/3h^3 and pF=ℏ(3π2ne)1/3p_F=\hbar(3\pi^2 n_e)^{1/3}; PNR=(3π2)2/3ℏ2ne5/3/5meP_\mathrm{NR}=(3\pi^2)^{2/3}\hbar^2 n_e^{5/3}/5m_e and PUR=(3π2)1/3ℏc ne4/3/4P_\mathrm{UR}=(3\pi^2)^{1/3}\hbar c\,n_e^{4/3}/4 from the pressure integral with the correct v(p)v(p); the RR-cancellation gives MCh∝(ℏc/G)3/2/(μemu)2M_\mathrm{Ch}\propto(\hbar c/G)^{3/2}/(\mu_e m_u)^2, with the Lane–Emden n=3n=3 coefficient fixing 1.44 M⊙1.44\,M_\odot for μe=2\mu_e=2. The worked-example anchors are ne=3.01×1035 m−3n_e=3.01\times10^{35}\ \mathrm{m^{-3}}, pF/mec=0.80p_F/m_e c=0.80, P∼3×1021 PaP\sim3\times10^{21}\ \mathrm{Pa}, TF=1.9×109 KT_F=1.9\times10^9\ \mathrm{K} at ρ=106 g cm−3\rho=10^6\ \mathrm{g\,cm^{-3}}.

Prerequisites. Fermi–Dirac Statistics · Pauli Exclusion Principle · Ideal Gas and Phase Space. Leads to. White Dwarfs · The Chandrasekhar Mass Limit · White Dwarf Cooling and Crystallization · Neutron Stars · The Equation of State of Dense Matter.

Ask the tutor: Want to see how the γ=4/3\gamma=4/3 softening turns into a hard mass limit, or how the same Fermi-gas logic rebuilds itself for neutron stars? Ask me to walk through either.

Is degeneracy pressure a new fundamental force?

No — it is ordinary kinetic pressure, computed with quantum bookkeeping. Pressure is momentum flux, and the only new ingredient is the Pauli exclusion principle, which forbids two electrons from sharing a quantum state and therefore forces them into high-momentum states once the low ones are full. Those large momenta produce a large flux, which we feel as pressure. There is no attraction or repulsion involved; in the ideal Fermi gas the electrons do not interact at all. So degeneracy pressure is not a fifth force but a statistical consequence of the exclusion principle acting on the electrons' kinetic energy — the same physics that makes metals hard also holds up white dwarfs.

Why does the pressure not depend on temperature?

Because at T=0T=0 the electrons already fill every momentum state up to pFp_F, and pFp_F is fixed by density alone through pF=ℏ(3π2ne)1/3p_F=\hbar(3\pi^2 n_e)^{1/3}. Raising the temperature can only nudge the thin layer of electrons within kBTk_B T of the Fermi surface; the vast degenerate sea beneath is frozen in place, since there are no empty lower states for those electrons to move into. When kBT≪EFk_B T\ll E_F the thermal correction to the pressure is of order (kBT/EF)2(k_B T/E_F)^2, utterly negligible for real white dwarfs. The pressure is therefore set by the density, and a degenerate star can cool indefinitely while keeping its structure — the reason white dwarfs are stable fossils rather than slowly deflating balloons.

Why does the equation of state soften at high density?

Because electrons cannot travel faster than light. In the pressure integral the momentum flux depends on the product p v(p)p\,v(p). At low density the fastest electrons obey v=p/mev=p/m_e, so flux grows as p2p^2 and the pressure rises steeply, P∝ne5/3P\propto n_e^{5/3}. But as density climbs, pFp_F eventually exceeds mecm_e c and the velocity saturates at cc; now the flux grows only as pp, one power slower, and the pressure softens to P∝ne4/3P\propto n_e^{4/3}. The electron mass drops out of the relativistic law entirely. This softening is not a minor correction — it is the difference between a gas that can always out-stiffen gravity and one that cannot, which is why it dictates the existence of a maximum mass.

Where does the Chandrasekhar mass come from physically?

From a coincidence of exponents. A self-gravitating body in hydrostatic balance needs its pressure to scale as ρ4/3\rho^{4/3} at the margin of stability, because both the central pressure and the gravitational load scale as M2/R4M^2/R^4. The non-relativistic degenerate gas is stiffer than this (γ=5/3\gamma=5/3), so it can always find a stable radius. But as the star grows more massive its center turns relativistic and the equation of state softens to exactly γ=4/3\gamma=4/3 — matching gravity's demand. At that point the radius cancels out of the force balance, and only one mass can satisfy it, MCh≈1.44 (2/μe)2 M⊙M_\mathrm{Ch}\approx1.44\,(2/\mu_e)^2\,M_\odot. Above this mass, gravity wins at every radius and collapse is inevitable.

How is neutron degeneracy pressure related?

It is the identical physics with the electron replaced by the neutron. When a white dwarf is pushed past the Chandrasekhar mass, electron captures (e−+p→n+νee^-+p\to n+\nu_e) convert the matter to neutrons, and neutron degeneracy pressure — the same Fermi-gas formulae with me→mnm_e\to m_n and ne→nnn_e\to n_n — takes over. Because the neutron is about 1839 times heavier, the non-relativistic pressure (carrying a factor 1/m1/m) is far weaker at a given density, so equilibrium is reached only at nuclear densities and radii of ∼10 km\sim10\ \mathrm{km}. A neutron star has its own maximum mass (the Tolman–Oppenheimer–Volkoff limit, ∼2\sim2–2.3 M⊙2.3\,M_\odot), but computing it requires general relativity and strong interactions, not just the ideal Fermi gas.

What is the helium flash, and how does degeneracy cause it?

In a low-mass red giant the inert helium core becomes electron-degenerate before it is hot enough to fuse helium, and degenerate pressure does not respond to temperature. In a normal gas a burst of nuclear heating raises the pressure, the region expands and cools, and burning is throttled — a thermostat. In a degenerate core there is no such relief valve: heating raises the temperature but not the pressure, so the reaction rate soars in a thermonuclear runaway called the helium flash, which briefly burns at a luminosity comparable to an entire galaxy though almost none escapes. Only when the temperature climbs enough to lift the degeneracy (T→TFT\to T_F) does the gas behave thermally again and the thermostat switch back on.

Where this leads