The Chandrasekhar Mass

The Chandrasekhar mass is not a historical curiosity; it is a live number at the centre of several frontier programmes. Its most consequential modern role is as the trigger of Type Ia supernovae. A carbon–oxygen white dwarf that accretes from a companion, or merges with another white dwarf, is driven toward $M_{\mathrm{Ch}}$; as it nears the limit the central density and temperature climb until carbon ignites, and because the electrons are degenerate the burning cannot regulate itself by expanding — it runs away, unbinding the entire star in seconds. Because the trigger mass is very nearly universal, the peak brightness is standardizable, and it was Type Ia supernovae, calibrated by exactly this physics, that revealed in 1998 that the expansion of the Universe is accelerating — the discovery behind the 2011 Nobel Prize and the whole industry of dark-energy cosmology. Wide surveys, from the Zwicky Transient Facility to the Vera Rubin Observatory's LSST, are now collecting these explosions by the thousands, while the question of exactly how a white dwarf reaches the limit — single-degenerate accretion versus double-degenerate merger — remains open. The same number governs the other end of stellar death: massive stars build degenerate iron cores that grow toward the temperature-corrected limit until electron capture removes the support and the core implodes into a neutron star or black hole, powering a core-collapse supernova. And the handful of reported "super-Chandrasekhar" Type Ia events are an active puzzle, most plausibly explained by rapid rotation, magnetic support, or mergers rather than any failure of the limit itself.

Imagine feeding mass onto a white dwarf, one spoonful at a time. Gravity squeezes harder, so the star shrinks, and the shrinking forces its electrons into ever higher-momentum quantum states — that is what degeneracy pressure is, the cost of stacking fermions into occupied phase space. As long as the electrons move slowly compared with light, they push back ferociously: the pressure climbs as $\rho^{5/3}$, faster than gravity's demand, and the star always finds a new, smaller equilibrium. This is the regime of the ordinary mass–radius relation, in which heavier white dwarfs are, paradoxically, smaller. But every extra spoonful drives the central density higher and the electrons faster, and there is a ceiling on speed. Once the most energetic electrons are moving at nearly $c$, extra momentum buys almost no extra speed, and the pressure gain slackens to $\rho^{4/3}$. That single change of exponent is the whole story. The first figure shows it directly: the equation of state bends from a steep $5/3$ slope to a shallower $4/3$ slope as the gas crosses into the relativistic regime.

Why does that shallower slope spell disaster? Because $\rho^{4/3}$ is precisely the scaling of gravity itself. For a self-gravitating ball the pressure gravity demands at the centre grows as $\rho^{4/3}$ too, so once the gas softens to that exponent, support and demand rise in lockstep and the radius drops out of the balance entirely. There is no longer a size the star can settle into to relieve the deficit; only the mass decides. Below a critical mass the gas wins by a constant factor at every radius and equilibrium exists; above it, gravity wins at every radius and the star has nowhere to stop. The second figure shows the consequence for real white dwarfs: as the mass climbs toward $M_{\mathrm{Ch}}=1.44\,M_\odot$, the equilibrium radius does not level off — it plunges toward zero.

A dimensional shortcut captures the whole result in one line. Balance the Fermi energy of a relativistic electron gas against the gravitational energy each electron carries, the radius cancels, and what survives is a pure condition on the number of electrons the star can hold: $N_{\max}\sim(\hbar c/Gm_H^2)^{3/2}\sim\alpha_G^{-3/2}$, where $\alpha_G=Gm_H^2/\hbar c\approx5.9\times10^{-39}$ is the gravitational fine-structure constant. The Chandrasekhar mass is just this enormous number of nucleons times the mass of one. The deepest lesson is the contrast with an ordinary star: a main-sequence star can always heat up to meet any gravitational demand, because thermal pressure rises with temperature. Degenerate matter cannot cheat that way — at zero temperature its pressure is fixed by density alone, and that rigidity is exactly what makes a maximum mass possible.

The headline result is a formula for the limiting mass that contains no property of any particular star — only constants of nature and the composition of the gas through $\mu_e$. It follows from treating the ultra-relativistic degenerate gas as an $n=3$ polytrope and solving the resulting stellar-structure equation once and for all.

where: $\hbar$ — reduced Planck constant; $c$ — speed of light; $G$ — gravitational constant; $\mu_e$ — mean molecular weight per electron (nucleons per electron), $\mu_e=2$ for fully ionized carbon–oxygen; $m_H$ — hydrogen/nucleon mass; $\omega_3^{0}=(-\xi^2\,\mathrm{d}\theta/\mathrm{d}\xi)_{\xi_1}\approx2.018$, the dimensionless surface constant of the $n=3$ Lane–Emden solution. Read it as: the maximum supportable mass is a fixed number of "gravitational quanta" $(\hbar c/G)^{3/2}$ divided by the square of the mass carried per electron. Plugging carbon–oxygen composition into the prefactor gives the canonical number. The small drop from $1.457$ to the commonly quoted $1.44$ comes from Coulomb-lattice and general-relativistic corrections applied on top of the ideal result.

where: $M_\odot$ — solar mass, $1.989\times10^{30}\,\mathrm{kg}$. This is the number to remember: roughly one and a half Suns, the mass scale of every white dwarf and neutron star. Behind the polytrope lies the exact zero-temperature equation of state of a degenerate electron gas, valid at any degree of relativity. It is what interpolates smoothly between the two limiting exponents.

where: $m_e$ — electron mass; $h$ — Planck constant; $x_F=p_F/m_ec$ — Fermi momentum in units of $m_ec$, the single parameter measuring how relativistic the gas is; $p_F=\hbar(3\pi^2 n_e)^{1/3}$ — Fermi momentum. Reading: the bracket reduces to $\tfrac{8}{5}x_F^5$ for $x_F\ll1$ (giving $P\propto\rho^{5/3}$) and to $2x_F^4$ for $x_F\gg1$ (giving $P\propto\rho^{4/3}$). The gas density and the electron density are linked by the composition. Pressure comes from the electrons, but the mass — and therefore gravity — is carried almost entirely by the nucleons bound to them, which is why $\mu_e$ rather than the full molecular weight sets the scale.

where: $n_e$ — electron number density; $\rho$ — mass density; $X_i$ — mass fraction of species $i$; $Z_i,A_i$ — its charge and mass number. For any element heavier than hydrogen with $Z/A\approx1/2$ (helium, carbon, oxygen), $\mu_e=2$; pure hydrogen gives $\mu_e=1$. In the non-relativistic regime the same physics yields the famous inverted mass–radius relation, the low-mass branch of the second figure.

where: $R$ — white-dwarf radius; $M$ — its mass. Reading: pile on mass and the star shrinks, driving the density up and pushing the electrons toward the relativistic regime where equation (1) takes over — the mechanism that carries a star toward its own limit. The physical origin of the limit is cleanest in energy language. For a fixed mass, both the gravitational binding energy and the ultra-relativistic degeneracy (kinetic) energy scale as $1/R$.

where: $a,b$ — order-unity coefficients fixed by the density profile. Reading: because both terms share the same $1/R$, the sign of the total energy is set by $M$ alone. For $M<M_{\mathrm{Ch}}$ the positive term wins and a bound equilibrium exists; for $M>M_{\mathrm{Ch}}$ the total energy is negative and falls without bound as $R\to0$ — runaway collapse. Finally, the whole scale can be read as a statement about the weakness of gravity, packaged in the gravitational fine-structure constant.

where: $\alpha_G$ — gravitational coupling between two nucleons, the analogue of the electromagnetic $\alpha\approx1/137$. Reading: a star can hold roughly $\alpha_G^{-3/2}\approx2\times10^{57}$ nucleons against gravity before relativity defeats degeneracy — and that pure number, times a proton mass, is a solar mass.

How the star behaves depends entirely on how relativistic its electrons are, which is measured by the single parameter $x_F=p_F/m_ec$.

1 · Evaluating the limit and its composition scaling (confidence builder)

The prefactor in equation (1) is $M_{\mathrm{Ch}}=5.83\,\mu_e^{-2}\,M_\odot$. Find the limit for a carbon–oxygen white dwarf, for a pure-helium white dwarf, and for a hypothetical pure-hydrogen one, and say what the pattern means.

Step 1. Fix $\mu_e$ for each composition using $\mu_e=(\sum_i X_i Z_i/A_i)^{-1}$. Carbon ($Z/A=6/12$), oxygen ($8/16$), and helium ($2/4$) all have $Z/A=1/2$, so $\mu_e=2$. Hydrogen has $Z/A=1$, so $\mu_e=1$. C/O and He both give $\mu_e=2$; H gives $\mu_e=1$.

Step 2. Insert into $M_{\mathrm{Ch}}=5.83/\mu_e^2\,M_\odot$. For $\mu_e=2$: $M_{\mathrm{Ch}}=5.83/4=\mathbf{1.46\,M_\odot}$ (the standard $1.44\,M_\odot$ after corrections). For $\mu_e=1$: $M_{\mathrm{Ch}}=5.83/1=\mathbf{5.83\,M_\odot}$.

Step 3. Interpret. A pure-helium and a carbon–oxygen white dwarf share the same limit because both carry two nucleons per electron — the electrons that provide pressure are equally "diluted" by weight. A hydrogen white dwarf would formally support four times as much mass, but no such object exists: hydrogen ignites long before it can degenerate at that density, so nature never realizes $\mu_e=1$. The limit is set by the electron-to-nucleon ratio, not by which specific ions are present — helium and carbon–oxygen white dwarfs cap at the same $1.44\,M_\odot$.

2 · Central density and Fermi energy at the limit

A white dwarf approaching $M_{\mathrm{Ch}}$ has a central density around $\rho_c\approx2\times10^{9}\,\mathrm{kg/m^3}$ ($2\times10^{6}\,\mathrm{g/cm^3}$) with $\mu_e=2$. Find the electron number density, the Fermi momentum in units of $m_ec$, and confirm the electrons are relativistic.

Step 1. Electron number density from equation (4): $n_e=\rho/(\mu_e m_H)=(2\times10^{9})/(2\times1.67\times10^{-27})=\mathbf{6.0\times10^{35}\,m^{-3}}$. That is about $10^{9}$ times denser than the conduction electrons in copper.

Step 2. Fermi momentum $p_F=\hbar(3\pi^2 n_e)^{1/3}$. Compute $3\pi^2 n_e=29.6\times6.0\times10^{35}=1.78\times10^{37}$, cube root $=2.6\times10^{12}\,\mathrm{m^{-1}}$, so $p_F=1.055\times10^{-34}\times2.6\times10^{12}=\mathbf{2.7\times10^{-22}\,kg\,m/s}$.

Step 3. Compare with $m_ec=9.11\times10^{-31}\times3.0\times10^{8}=2.7\times10^{-22}\,\mathrm{kg\,m/s}$. Therefore $x_F=p_F/m_ec\approx\mathbf{1.0}$. The Fermi momentum equals $m_ec$: the most energetic electrons are exactly transrelativistic, right where the equation of state is bending from slope $5/3$ to $4/3$ in the first figure.

Step 4. Interpret. This is not a coincidence — it is the definition of the limit. A star reaches $M_{\mathrm{Ch}}$ precisely when its central electrons become relativistic, because that is when the supporting pressure softens to gravity's scaling. Equivalently, in the language of equation (7), it can hold about $\alpha_G^{-3/2}\approx2\times10^{57}$ nucleons — roughly $1.85\,M_\odot$ before the $\mu_e^{-2}$ and polytrope factors sharpen it — after which relativity defeats degeneracy. At the limit the central $x_F\sim1$: the star sits exactly at the crossover density where degeneracy pressure loses its edge over gravity, and the mass of a dying star is $\alpha_G^{-3/2}$ protons.

The Chandrasekhar mass is inseparable from the young man who found it. It is one of the rare fundamental results in astrophysics attached to a specific voyage, a famous act of resistance, and a Nobel Prize awarded half a century later.

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., §16.4 (White Dwarfs) — the standard undergraduate treatment with the mass–radius relation · Shapiro & Teukolsky, Black Holes, White Dwarfs, and Neutron Stars, Ch. 3–4 — the definitive graduate derivation of the degenerate EOS and the limit · Chandrasekhar, An Introduction to the Study of Stellar Structure (1939) — the original polytrope theory from its author · Dodelson, Modern Cosmology — Type Ia supernovae as standard candles · Mo, van den Bosch & White, Galaxy Formation and Evolution — stellar remnants in cosmic context · Binney & Tremaine, Galactic Dynamics — degenerate stars within stellar-dynamical systems. The key numerical values used throughout this article are: prefactor $M_{\mathrm{Ch}}=5.83\,\mu_e^{-2}\,M_\odot$; $n=3$ Lane–Emden surface constants $\xi_1=6.897$ and $(-\xi^2\theta')_{\xi_1}=2.018$; canonical value $1.44\,M_\odot$ for $\mu_e=2$ (equation (1) in SI units gives $2.9\times10^{30}$ kg $=1.46\,M_\odot$ before corrections); gravitational coupling $\alpha_G\approx5.9\times10^{-39}$; and $\alpha_G^{-3/2}m_H\approx1.85\,M_\odot$. The Fermi-gas pressure integral of equation (3) reduces to $\rho^{5/3}$ for $x_F\ll1$ and to $\rho^{4/3}$ for $x_F\gg1$.

Prerequisites. White Dwarfs and Electron Degeneracy Pressure · Asymptotic Giant Branch and Thermal Pulses · Advanced Burning Stages in Massive Stars. Leads to. Type Ia Supernovae · Neutron Stars and the Tolman–Oppenheimer–Volkoff Limit.

Ask the tutor: Ask me to walk through why the $\rho^{4/3}$ equation of state makes the balance mass independent of radius, or to work a Type Ia ignition estimate as a white dwarf approaches the limit.