Chandrasekhar Limit

The Chandrasekhar limit (denoted M_Ch ≈ 1.44 M_☉, or approximately 2.86 × 10³⁰ kg for an ideal carbon–oxygen or helium composition with mean molecular weight per electron μ_e = 2.0) is the maximum mass that a non-rotating, unmagnetized white dwarf star—or the electron-degenerate core of an evolved star—can sustain in hydrostatic equilibrium against gravitational self-collapse via quantum-mechanical electron degeneracy pressure. Expressed in terms of the fundamental constants of quantum mechanics (ℏ), special relativity (c), gravitation (G), and the proton/hydrogen mass (m_H), the limit scales with the Planck mass cubed divided by the nucleon mass squared: M_Ch ∝ (ℏc/G)^3/2 m_H⁻² μ_e⁻².

The limit was calculated in July 1930 by the 19-year-old Indian astrophysicist Subrahmanyan Chandrasekhar during his voyage from Madras to England aboard the steamship SS Pilsna to begin graduate studies at Trinity College, Cambridge (work for which he was awarded the 1983 Nobel Prize in Physics, following related independent polytropic calculations by Wilhelm Anderson [1929] and Edmund C. Stoner [1930]). Earlier, in 1926, Ralph H. Fowler had resolved the paradox of white dwarf stars such as Sirius B—which compresses a solar mass into a planetary volume comparable to Earth (ρ̄ ~ 10⁹ kg m⁻³), left behind after a solar-type star sheds a planetary nebula such as the Ring Nebula—by applying the newly formulated Fermi–Dirac statistics and Wolfgang Pauli's exclusion principle (P_deg ∝ ρ^5/3). Fowler's non-relativistic treatment predicted that a white dwarf of any mass could reach a stable equilibrium radius by contracting to a higher density.

Chandrasekhar realized that inside white dwarfs approaching 1 M_☉, the electron Fermi momentum p_F approaches and exceeds m_ec, meaning that the degenerate electrons move at velocities approaching the speed of light c. Incorporating Albert Einstein's special relativity into Fermi–Dirac statistical mechanics softened the equation of state at high densities from P ∝ ρ^5/3 to P ∝ ρ^4/3, proving the existence of a finite upper mass ceiling above which no stable white dwarf configuration exists. Despite public opposition from Sir Arthur Stanley Eddington at the January 1935 meeting of the Royal Astronomical Society—who dismissed relativistic degeneracy as "stellar buffoonery" because it implied the gravitational collapse of massive stars into what are now known as neutron stars and black holes—astronomical observations and quantum electrodynamics vindicated Chandrasekhar's derivation.

---

Fermi Momentum and the Pauli Exclusion Principle

Consider an ionized, cold (T → 0 K) stellar plasma of mass density ρ, composed of atomic nuclei with mass number A and atomic number Z alongside free electrons. Electrical neutrality requires an electron number density n_e related to mass density by:

n_e = ( Z / A ) ( ρ / m_H ) = ρ / ( μ_e m_H )

where μ_e ≡ A/Z is the mean molecular weight per electron (μ_e = 2.00 for ⁴He, ¹²C, and ¹⁶O; μ_e = 56/26 = 2.154 for ⁵⁶Fe). By the Pauli exclusion principle and the four quantum numbers of an electron, at most two electrons with opposite spin projections (m_s = +½, −½) can occupy each elementary phase-space cell of volume h³ = (2*πℏ*)³. At zero temperature, electrons fill momentum states uniformly inside a sphere in momentum space up to the Fermi momentum p_F:

n_e = ( 2 / h³ ) ∫₀^p_F 4*πp*² dp = ( 8*π* p_F³ ) / ( 3 h³ ) ⇒ p_F = ℏ ( 3*π*² n_e )^1/3

From kinetic theory, the isotropic pressure exerted by an ideal Fermi gas of particles with momentum p and velocity v(p) = ∂E_k/∂p is:

P_deg = ( 1 / 3 ) ∫₀^p_F p v(p) ( 8*πp*² / h³ ) dp

Regime 1: Non-Relativistic Degeneracy (p_F ≪ m_ec, Low-Mass White Dwarfs)

In low-mass white dwarfs (M ≲ 0.5 M_☉, ρ ≲ 10⁹ kg m⁻³), the Fermi momentum is small compared to m_ec, so the electron kinetic energy is Newtonian, E_k = p²/(2*m*_e), and the velocity is v(p) = p/m_e. Integrating from 0 to p_F yields the non-relativistic degenerate equation of state (a polytrope P = K₁ρ^γ with adiabatic index γ = 5/3 and polytropic index n = 1/(γ − 1) = 3/2):

P_deg,NR = [ ℏ² / ( 5 m_e ) ] ( 3*π*² )^2/3 n_e^5/3 = K₁ ( ρ / μ_e )^5/3

Balancing outward pressure gradient P/R against inward Newtonian gravity GMρ/R² with ρ ∝ M/R³ gives:

P_deg,NR ∝ M^5/3 / R⁵ ≈ G M² / R⁴ ⇒ R ∝ M^−1/3 μ_e^−5/3

This is the inverted white dwarf mass–radius relation: adding mass to a white dwarf causes its equilibrium radius R to shrink (V ∝ R³ ∝ 1/M) and its central density to climb as ρ_c ∝ M².

Regime 2: Ultra-Relativistic Degeneracy (p_F ≫ m_ec) and the Chandrasekhar Limit

As mass increases and the white dwarf contracts to densities ρ ≫ 2 × 10⁹ kg m⁻³, p_F surpasses m_ec. In the ultra-relativistic limit (p_F ≫ m_ec), the relativistic energy–momentum relation E = √(p²c² + m_e²c⁴) approaches E ≈ pc, and the electron velocity saturates at the speed of light, v(p) = ∂E/∂p ≈ c. Substituting v(p) = c into the pressure integral softens the equation of state from n_e^5/3 to n_e^4/3 (adiabatic index γ = 4/3, polytropic index n = 3):

P_deg,UR = ( ℏc / 4 ) ( 3*π*² )^1/3 n_e^4/3 = K₂ ( ρ / μ_e )^4/3

Now consider virial hydrostatic equilibrium between the central ultra-relativistic degeneracy pressure (P_deg,UR ∝ ρ^4/3 ∝ [M/R³]^4/3 = M^4/3/R⁴) and the self-gravitational pressure (P_grav ∝ GM²/R⁴):

P_deg,UR / P_grav ∝ ( ℏc M^4/3 / R⁴ ) / ( G M² / R⁴ ) ∝ ( ℏc / G ) M^−2/3

Notably, the stellar radius R cancels out on both sides of the equation. Because both outward ultra-relativistic degeneracy pressure and inward gravity scale as 1/R⁴, a star whose mass exceeds the threshold where P_grav > P_deg,UR cannot halt its contraction by shrinking to a smaller radius; it will contract toward R → 0.

Solving the Lane–Emden differential equation (1/ξ²) d/dξ (ξ² dθ/dξ) = −*θ*ⁿ for polytropic index n = 3 yields the dimensionless homology constant ω₃⁰ = −*ξ*²(dθ/dξ)|_ξ1 = 2.018236 at the surface zero ξ₁ = 6.89685. Substituting K₂ into the Lane–Emden mass integral gives the Chandrasekhar mass:

M_Ch = [ ω₃⁰ √(3*π*) / 2 ] ( ℏc / G )^3/2 ( 1 / [ μ_e m_H ] )² ≈ 5.836 / μ_e² M_☉

For a pure helium, carbon, or oxygen white dwarf (μ_e = 2.00), the ideal relativistic Chandrasekhar limit is:

M_Ch (μ_e = 2.00) = 5.836 / 4.00 = 1.459 M_☉ (~1.44 M_☉ using ¹²C atomic mass unit convention)

Using the exact relativistic Fermi-gas parametric function f(x) (where x = p_F/(m_ec)), Chandrasekhar (1935) derived the complete white-dwarf mass–radius curve interpolating smoothly from R ∝ M^−1/3 at low mass to R → 0 as M → M_Ch:

R(M) ≈ R₀ ( M / M_Ch )^−1/3 [ 1 − ( M / M_Ch )^4/3 ]^1/2

---

Coulomb Lattice, Electron Capture, and General Relativistic Corrections

In real astrophysical white dwarfs and evolved stellar cores, three physical effects modify the ideal 1.459 M_☉ value:

1. Electrostatic Coulomb Lattice Attraction (Salpeter 1961): Because positive C/O nuclei are localized in a body-centered cubic crystal lattice rather than spread uniformly, electrostatic attraction slightly lowers the net pressure by ~1.5–2%, reducing M_Ch from 1.459 M_☉ to ~1.43 M_☉.
2. Inverse Beta Decay / Pycnonuclear Electron Capture: At central densities ρ_c > 10¹³ kg m⁻³, the electron Fermi energy exceeds the threshold to trigger endothermic electron capture onto nuclei (p + e⁻ → n + ν_e, such as ¹⁶O → ¹⁶N or ²⁴Mg → ²⁴Na). Removing degenerate electrons raises μ_e and drops the effective adiabatic index below γ < 4/3, triggering instability before R = 0 is reached.
3. General Relativistic Hydrostatic Instability (Chandrasekhar 1964): Because general relativity steepens gravity relative to Newtonian mechanics, dynamical instability sets in slightly before the adiabatic index drops all the way to 4/3 (at γ_crit = 4/3 + 1.12 r_s/R), capping the maximum stable mass of a cold non-rotating C–O white dwarf at M_Ch,eff ≈ 1.38–1.40 M_☉ (at a finite minimum radius of R_min ≈ 1,800–2,000 km and central density ρ_c ≈ 2.5 × 10¹³ kg m⁻³).
4. Iron Core Composition (μ_e = 2.154): Inside the evolved iron core (⁵⁶Fe, Z = 26, A = 56) of a massive supergiant star, the higher neutron-to-proton ratio yields μ_e = 56/26 = 2.154, lowering the bare Chandrasekhar limit to 5.836 / (2.154)² = 1.258 M_☉ (though finite thermal entropy in hot supergiant cores can temporarily raise the effective core mass to 1.35–1.60 M_☉ prior to collapse).

Role in Type Ia Supernovae and Core-Collapse Supernovae

The Chandrasekhar limit serves as the primary bifurcation boundary of stellar evolution:

• Low- and Intermediate-Mass Stars (M_ZAMS ≲ 8 M_☉): Shed their outer envelopes on the Asymptotic Giant Branch to leave behind C–O white dwarfs with M ≈ 0.55–1.05 M_☉ < M_Ch, which cool eternally along the white dwarf sequence without further nuclear burning.
• Accreting or Merging White Dwarfs (Type Ia Supernovae): When a C–O white dwarf in a binary system accretes mass toward ~1.38 M_☉, core compression to ρ_c ≈ 2–3 × 10¹² kg m⁻³ ignites a thermonuclear runaway of carbon fusion before collapse can occur, exploding the star as a standardizable Type Ia supernova. (By contrast, an oxygen–neon–magnesium [O–Ne–Mg] white dwarf reaching M_Ch undergoes electron capture onto ²⁴Mg and ²⁰Ne, collapsing via Accretion-Induced Collapse [AIC] into a neutron star).
• Massive Stars (M_ZAMS ≳ 8 M_☉): Fuse elements all the way to an iron-peak core; once the iron core exceeds M_Ch,Fe ≈ 1.26–1.40 M_☉, endothermic photodisintegration and electron capture trigger gravitational core collapse past white-dwarf densities into a neutron star or, if the remnant exceeds the Tolman–Oppenheimer–Volkoff limit (M_TOV ≈ 2.15–2.30 M_☉), past its Schwarzschild radius into a black hole.

---

Related concepts include Hertzsprung–Russell Diagram and Dark Matter.

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

In a degenerate white dwarf, outward pressure is provided by the quantum-mechanical Pauli exclusion principle rather than thermal heat: electrons are squeezed into momentum states up to the Fermi momentum p_F ∝ n_e^1/3. When mass is added to a white dwarf, its self-gravity increases as GM²/R⁴. To generate enough additional electron degeneracy pressure (P_deg,NR ∝ ρ^5/3 ∝ M^5/3/R⁵) to balance the stronger gravity, the star must contract to a smaller radius according to R ∝ M^−1/3.

Why did incorporating special relativity introduce a maximum mass limit for white dwarfs?

In non-relativistic quantum mechanics, compressing a star by a factor R increases electron velocities without bound (v ∝ p_F ∝ 1/R), causing degeneracy pressure to grow as 1/R⁵—always outpacing gravity (which grows as 1/R⁴). In special relativity, however, electron velocities cannot exceed the speed of light c. As v → c, the equation of state softens from P ∝ ρ^5/3 (1/R⁵) to P ∝ ρ^4/3 (1/R⁴). Because pressure and gravity now scale with the exact same power of radius, any star heavier than ~1.44 M_☉ cannot halt collapse at any radius.

Can rapidly rotating or strongly magnetized white dwarfs exceed 1.44 solar masses?

Yes. The classical 1.44 M_☉ Chandrasekhar limit applies strictly to a non-rotating, spherically symmetric, unmagnetized star. Differential centrifugal force in rapidly spinning white dwarfs—commonly formed by the double-degenerate merger of two white dwarfs—can support super-Chandrasekhar white dwarfs with masses up to ~2.0–2.4 M_☉, explaining overluminous Type Ia supernovae such as SN 2003fg and SN 2009dc.

---

1. Fowler, R. H. (1926). "On dense matter." Monthly Notices of the Royal Astronomical Society, 87(2), 114–122. https://doi.org/10.1093/mnras/87.2.114
2. Chandrasekhar, S. (1931). "The maximum mass of ideal white dwarfs." The Astrophysical Journal, 74, 81–82. https://doi.org/10.1086/143324
3. Chandrasekhar, S. (1935). "The highly collapsed configurations of a stellar mass (Second paper)." Monthly Notices of the Royal Astronomical Society, 95(3), 207–225. https://doi.org/10.1093/mnras/95.3.207
4. Chandrasekhar, S. (1984). "On stars, their evolution and their stability" (Nobel Lecture). Reviews of Modern Physics, 56(2), 137–147. https://doi.org/10.1103/RevModPhys.56.137