Chandrasekhar Mass Limit

Every star that fuses hydrogen leaves a corpse, and for the great majority — the Sun included — that corpse is a white dwarf: the naked, cooling core exposed when the envelope drifts away. Nothing burns in it any longer, so the one thing between the star and gravitational ruin is a quantum effect. The Pauli exclusion principle forbids two electrons from sharing a state, and packing them together forces them into ever higher momenta; the resulting electron degeneracy pressure holds the star up with no need for heat at all. Non-relativistically that pressure stiffens as $P\propto\rho^{5/3}$ , steeply enough to guarantee an equilibrium at any mass, and it produces the strange signature relation $R\propto M^{-1/3}$ — heavier white dwarfs are smaller. That guarantee has a hidden expiry date. Squeeze the star hard enough and the electrons' Fermi momentum climbs past $m_e c$ ; they turn relativistic, their energy grows as $pc$ rather than $p^2/2m_e$ , and the equation of state softens to $P\propto\rho^{4/3}$ — exactly the power law at which self-gravity demands support. When both energies scale identically with size, balance can be met at one and only one mass, the Chandrasekhar mass $M_{\mathrm{Ch}}\approx1.4\,M_\odot$ . Above it, no cold equilibrium exists and gravity always wins. The number is built almost entirely from constants of Nature — $\hbar$ , $c$ , $G$ , and the nucleon mass — the only material input being $\mu_e$ , the nucleons per electron. That universality is why Type Ia supernovae, triggered when a white dwarf is pushed to $M_{\mathrm{Ch}}$ , all detonate at nearly the same mass and shine with nearly the same peak brightness: the standard candles that revealed the accelerating Universe.

Key takeaways

The Chandrasekhar mass is the maximum mass a white dwarf's electron degeneracy pressure can support against gravity: $M_{\mathrm{Ch}}\approx1.4\,M_\odot$ for a carbon–oxygen core.
It exists because at high density the electrons turn ultra-relativistic, softening the equation of state from $P\propto\rho^{5/3}$ to $P\propto\rho^{4/3}$ — the exact scaling of gravity's demand.
Its value depends only on fundamental constants and the nucleons-per-electron $\mu_e$ : $M_{\mathrm{Ch}}\propto(\hbar c/G)^{3/2}/(\mu_e m_H)^2=1.457\,(2/\mu_e)^2\,M_\odot$ .
The equilibrium radius formally shrinks to zero at $M_{\mathrm{Ch}}$ ; real stars are destabilised slightly earlier by inverse beta decay and general relativity.
White dwarfs driven toward $M_{\mathrm{Ch}}$ by accretion or merger ignite carbon and explode as Type Ia supernovae — the cosmological standard candles.
Key idea: in the ultra-relativistic limit both degeneracy pressure and gravity scale as $\rho^{4/3}$ , so balance can hold at only one mass, $M_{\mathrm{Ch}}\approx1.4\,M_\odot$ , above which pressure can never again catch gravity.

The core insight

In the ultra-relativistic limit both degeneracy pressure and gravity scale as $\rho^{4/3}$ ; equality can then hold at only one mass, $M_{\mathrm{Ch}}\approx1.4\,M_\odot$ , above which no compression lets electron pressure balance gravity and the star must collapse.

Where this lives today

The Chandrasekhar mass is a live parameter in today's astrophysics. Type Ia supernovae, whose progenitor white dwarfs are driven near $M_{\mathrm{Ch}}$ , remain the premier tool for measuring cosmic distances; the 2011 Nobel Prize for dark energy rested on their near-uniform luminosity, and surveys like the Zwicky Transient Facility and Vera Rubin Observatory's LSST now find thousands per year. How they ignite is still contested — a single-degenerate channel (a dwarf accreting to the limit) competes with a double-degenerate channel (two dwarfs merging), and "super-Chandrasekhar" events like SN 2009dc have forced models with rapid rotation or magnetic support that raise the effective limit. Gaia's astrometry has delivered masses and radii for hundreds of thousands of white dwarfs, testing the mass–radius relation directly and revealing the sharp cutoff near $1.4\,M_\odot$ . And the same instability, applied to an iron stellar core, triggers core-collapse supernovae and the birth of every neutron star and stellar black hole — continued by the Tolman–Oppenheimer–Volkoff limit.

The physical picture

Imagine slowly adding mass to a cold white dwarf and watching the energy bookkeeping. Two energies compete: the degenerate electrons carry a kinetic energy that resists compression, and gravity carries a binding energy that rewards it. Non-relativistically the kinetic energy scales as $R^{-2}$ while gravity scales only as $R^{-1}$ ; the steeper term wins at small $R$ , so a genuine energy minimum exists and the star settles into a stable radius. But as the star is loaded up and the electrons turn relativistic, the kinetic energy softens until it too scales as $R^{-1}$ — the same power as gravity. Now the whole energy is a single constant times $1/R$ , and the star's fate hangs on the sign of that constant.

Figure

Ultra-relativistic energy balance E_tot = C/R: because degeneracy and gravity both scale as 1/R, only the sign of the constant C survives. C > 0 below the limit and the star expands until non-relativistic stiffening halts it at a finite radius, C = 0 exactly at M_Ch (marginal, radius undetermined), and C < 0 above it, where the energy falls without bound as R → 0 and the star collapses

If the star is lighter than $M_{\mathrm{Ch}}$ the constant is positive, the energy falls as the star swells, and expansion continues until the lower-density electrons become non-relativistic again and the stiffer $R^{-2}$ term reasserts itself and pins the radius. If it is heavier, the constant is negative and the energy plunges without bound as $R\to0$ . Exactly at $M_{\mathrm{Ch}}$ the constant vanishes: the star is marginally stable at every radius. This is why $M_{\mathrm{Ch}}$ is a mass limit — the balance depends on mass alone. The second figure shows the observable fingerprint: white dwarfs marching toward zero radius as their mass approaches $1.44\,M_\odot$ .

Figure

Zero-temperature mass–radius relation for a carbon–oxygen white dwarf (μ_e = 2): heavier means smaller, following R ∝ M^(−1/3) at low mass, and the radius plunges to zero as the mass reaches M_Ch = 1.44 M⊙. Measured degenerate stars — Sirius B (1.02 M⊙, ≈ 5800 km), 40 Eridani B, Procyon B, Stein 2051 B — lie on the curve, and not one exceeds the limit

The core equations

The headline result is a mass assembled almost entirely from fundamental constants. Everything material has collapsed into one composition factor $\mu_e$ ; the rest is $\hbar$ , $c$ , $G$ , and the nucleon mass, sharpened by a single number from the $n=3$ Lane–Emden solution.

M_{\mathrm{Ch}} = \frac{\sqrt{3\pi}}{2}\,\omega_3^0\left(\frac{\hbar c}{G}\right)^{3/2}\frac{1}{(\mu_e m_H)^2} \approx 1.457\left(\frac{2}{\mu_e}\right)^2 M_\odot \qquad (1)

where: $\hbar$ — reduced Planck constant; $c$ — speed of light; $G$ — gravitational constant; $\mu_e$ — mean nucleons per electron ( $\approx2$ for He/C/O); $m_H$ — atomic mass unit ( $\approx$ nucleon mass, $1.66\times10^{-27}\ \mathrm{kg}$ ); $\omega_3^0=2.018$ — Lane–Emden mass constant $(-\xi^2\theta')_{\xi_1}$ for $n=3$ . Temperature and opacity do not appear — only constants and composition. The two equations of state bracketing the problem come from the same Fermi-gas integral at opposite ends of the velocity range. At low density the electrons are slow and pressure stiffens as the five-thirds power of density.

P_{\text{nr}} = \frac{1}{20}\left(\frac{3}{\pi}\right)^{2/3}\frac{h^2}{m_e}\left(\frac{\rho}{\mu_e m_H}\right)^{5/3} \qquad (2)

where: $h$ — Planck constant; $m_e$ — electron mass; $\rho$ — mass density; $\rho/(\mu_e m_H)=n_e$ — electron number density. The $5/3$ exponent makes low-mass white dwarfs stable and gives $R\propto M^{-1/3}$ . At high density the electrons move at nearly $c$ ; the integral loses a power of momentum and softens to the four-thirds power — the pivotal change in the whole subject.

P_{\text{ur}} = \frac{1}{8}\left(\frac{3}{\pi}\right)^{1/3}h c\left(\frac{\rho}{\mu_e m_H}\right)^{4/3} \qquad (3)

where: $c$ — speed of light, entering because $v\to c$ ; the coefficient is a pure number times $hc$ . Note what is absent: $m_e$ has dropped out. The ultra-relativistic pressure knows only $\hbar$ , $c$ , and the number density — which is exactly why the resulting mass limit is built from constants alone. The instability is cleanest in energy form: add the relativistic kinetic energy of $N_e$ electrons to the gravitational binding of mass $M$ , and both carry the same $1/R$ .

E_{\mathrm{tot}}(R) = \underbrace{A\,\frac{\hbar c\,N_e^{4/3}}{R}}_{\text{degeneracy}} - \underbrace{B\,\frac{GM^2}{R}}_{\text{gravity}} = \frac{C}{R} \qquad (4)

where: $N_e$ — number of electrons ( $\propto M$ ); $R$ — stellar radius; $A,B$ — order-unity constants; $C\equiv A\hbar cN_e^{4/3}-BGM^2$ . Because $M\propto N_e$ , the gravity term grows as $N_e^2$ , faster than degeneracy's $N_e^{4/3}$ ; at some $N_e$ the constant $C$ turns negative — that is the limit. For a precise value one solves the structure equation: hydrostatic equilibrium plus the polytropic law $P=K\rho^{1+1/n}$ gives the Lane–Emden equation, and the ultra-relativistic gas is the polytrope $n=3$ .

\frac{1}{\xi^2}\frac{\mathrm{d}}{\mathrm{d}\xi}\!\left(\xi^2\frac{\mathrm{d}\theta}{\mathrm{d}\xi}\right) = -\theta^{\,n},\qquad n=3 \qquad (5)

where: $\xi$ — dimensionless radius; $\theta$ — dimensionless variable with $\rho=\rho_c\theta^{\,n}$ ; the $n=3$ solution has its first zero at $\xi_1=6.897$ and $(-\xi^2\theta')_{\xi_1}=2.018$ . This equation is where the number $1.457$ comes from. Two more relations complete the picture: on the stable branch the star obeys a backwards mass–radius law, and its stability is set by one adiabatic exponent.

R \propto M^{-1/3},\qquad \Gamma_1 > \tfrac{4}{3}\ \Leftrightarrow\ \text{dynamical stability} \qquad (6)

where: $\Gamma_1=(\partial\ln P/\partial\ln\rho)_{\text{ad}}$ — adiabatic exponent. A non-relativistic gas has $\Gamma_1=5/3$ (safely stable); the ultra-relativistic gas sits at exactly the marginal $\Gamma_1=4/3$ , so any further softening (neutronisation, GR) tips it below $4/3$ and triggers collapse.

Derivation AWhy the equation of state softens from five-thirds to four-thirds

The goal is to show, from the single Fermi-gas pressure integral, that a degenerate electron gas obeys $P\propto n_e^{5/3}$ when its electrons are slow and $P\propto n_e^{4/3}$ when they are relativistic, and that this softening is a purely kinematic consequence of special relativity — not a thermal, gravitational, or many-body effect. That change of one-third in the exponent is the entire physical origin of the Chandrasekhar limit.

Step 1. At zero temperature every momentum state up to the Fermi momentum $p_F$ is filled and every state above is empty. Counting two spin states per phase-space cell, $n_e=\dfrac{8\pi}{3h^3}p_F^{\,3}$ , so $p_F\propto n_e^{1/3}$ . Why: the Fermi momentum is fixed by how densely the electrons are packed — compression forces them into higher-momentum states, the microscopic root of degeneracy pressure.

Step 2. Pressure is the momentum flux these electrons carry: $P=\dfrac{1}{3}\displaystyle\int_0^{p_F} p\,v(p)\,\dfrac{8\pi p^2}{h^3}\,\mathrm{d}p$ , where $v(p)$ is the electron speed. Why: the whole difference between the two regimes lives in the velocity law $v(p)$ — how fast a particle carries momentum depends on whether it is relativistic.

Step 3. Non-relativistic limit ( $p_F\ll m_e c$ ): $v=p/m_e$ , so the integrand goes as $p^4$ and $P\propto p_F^{\,5}\propto n_e^{5/3}$ . Restoring constants gives equation (2). Why: the extra factor of $p$ from $v=p/m_e$ pushes the integral to the fifth power of $p_F$ , and the cube-root link $p_F\propto n_e^{1/3}$ turns that into the five-thirds law that keeps low-mass dwarfs stable.

Step 4. Ultra-relativistic limit ( $p_F\gg m_e c$ ): now $v\to c$ is constant, the integrand goes only as $p^3$ , and $P=\dfrac{2\pi c}{3h^3}p_F^{\,4}\propto n_e^{4/3}$ , equation (3). Why: pinned at the speed of light, the electrons carry no extra momentum flux as they are squeezed, so the gas loses a full power of $p_F$ — and with it the ability to out-stiffen gravity.

Step 5. The electron mass $m_e$ has vanished from the ultra-relativistic coefficient, which depends only on $h$ , $c$ , and $n_e$ . Why: this is what makes $M_{\mathrm{Ch}}$ a number built from $\hbar$ , $c$ , $G$ , and the nucleon mass alone — the electron sets the density where relativity switches on, not the value of the limit.

Go deeper · Derivation B · 1 minOrder-of-magnitude Chandrasekhar mass from energy balance

Derivation BOrder-of-magnitude Chandrasekhar mass from energy balance

The aim here is to recover the scaling $M_{\mathrm{Ch}}\sim(\hbar c/G)^{3/2}/(\mu_e m_H)^2$ from dimensional estimates of the two competing energies, and to see why the limit is a statement about a critical number of particles rather than a radius or density. Treat the star as a uniform sphere of $N_e$ electrons, total mass $M\approx\mu_e m_H N_e$ and radius $R$ , and demand that the degeneracy and gravitational energies trade off — essentially the calculation Chandrasekhar did on the boat to England.

Step 1. The electron number density is $n_e\sim N_e/R^3$ , so $p_F\sim\hbar n_e^{1/3}\sim\hbar N_e^{1/3}/R$ . Why: the uncertainty principle in a volume $R^3$ fixes the typical momentum, which grows as the star shrinks — degeneracy energy fights compression.

Step 2. In the ultra-relativistic limit each electron carries energy $\sim p_F c$ , so $E_{\mathrm{deg}}\sim N_e p_F c\sim\hbar c\,N_e^{4/3}/R$ . Why: the kinetic budget scales as $1/R$ — the crucial power that matches gravity.

Step 3. The gravitational binding is $E_{\mathrm{grav}}\sim -GM^2/R\sim -G(\mu_e m_H N_e)^2/R$ . Why: gravity scales as $N_e^2$ , faster than degeneracy's $N_e^{4/3}$ , so beyond some $N_e$ gravity must dominate.

Step 4. Add them: $E_{\mathrm{tot}}\sim\dfrac{1}{R}\left(\hbar c\,N_e^{4/3}-G\mu_e^2 m_H^2 N_e^2\right)$ . Both terms scale as $1/R$ , so no minimum exists — the sign of the bracket decides: positive expands, negative collapses to $R\to0$ . Why: the shared $1/R$ dependence is why there is a sharp threshold rather than a gentle equilibrium — the radius drops out and only $N_e$ matters.

Step 5. Setting the bracket to zero gives the critical particle number $N_{e,\max}\sim\left(\dfrac{\hbar c}{G}\right)^{3/2}(\mu_e m_H)^{-3}\approx2\times10^{57}$ . Why: this is essentially the number of nucleons in a Chandrasekhar mass, fixed by the ratio of $\hbar c$ to the gravitational scale $G m_H^2$ .

Step 6. Convert to mass: $M_{\mathrm{Ch}}=\mu_e m_H N_{e,\max}\sim(\hbar c/G)^{3/2}/(\mu_e m_H)^2$ , of order a solar mass. Why: the crude estimate already lands within a factor of a few of $1.4\,M_\odot$ ; the missing order-unity number is what the Lane–Emden solution supplies.

Go deeper · Derivation C · 1 minExact mass from the n = 3 polytrope

Derivation CExact mass from the n = 3 polytrope

This derivation turns the scaling estimate into the precise coefficient $1.457\,(2/\mu_e)^2\,M_\odot$ by solving the real structure equation. The key fact is that the ultra-relativistic gas is a polytrope of index exactly $n=3$ , and for that single index the polytropic mass formula loses all dependence on central density — so the star can have only one mass. We assemble the general polytrope mass, specialise to $n=3$ , insert the ultra-relativistic pressure constant, and evaluate with tabulated Lane–Emden numbers.

Step 1. For any polytrope $P=K\rho^{(n+1)/n}$ , integrating the Lane–Emden solution gives the total mass $M=4\pi\left[\dfrac{(n+1)K}{4\pi G}\right]^{3/2}\rho_c^{(3-n)/2n}\,(-\xi^2\theta')_{\xi_1}$ , with $\theta'$ read at the surface $\xi_1$ . Why: this packages hydrostatic equilibrium into one closed expression — the bracket carries the pressure law, the Lane–Emden factor the geometry of the sphere.

Step 2. Set $n=3$ . The exponent on the central density is $(3-n)/2n=0$ , so $\rho_c$ drops out completely and the mass is a single fixed number. Why: this is the mathematical heart of the limit — the mass is independent of how compressed the star is, so there is exactly one allowed value, $M_{\mathrm{Ch}}$ .

Step 3. Read the pressure constant off equation (3): $K_2=\dfrac{hc}{8}\left(\dfrac{3}{\pi}\right)^{1/3}(\mu_e m_H)^{-4/3}$ , so the bracket becomes $M=4\pi\left(\dfrac{K_2}{\pi G}\right)^{3/2}(-\xi^2\theta')_{\xi_1}$ . Why: substituting the ultra-relativistic $K_2$ — the one with no electron mass — makes the answer depend only on $\hbar$ , $c$ , $G$ , and $\mu_e m_H$ .

Step 4. Insert the tabulated $n=3$ data, $\xi_1=6.89685$ and $(-\xi^2\theta')_{\xi_1}=2.01824$ . Why: these two numbers are the only geometry the answer needs, and are why $M_{\mathrm{Ch}}$ is a definite number rather than a scaling.

Step 5. Evaluating the constants gives $M_{\mathrm{Ch}}=1.457\,(2/\mu_e)^2\,M_\odot$ . Composition enters only through $(2/\mu_e)^2$ , so a He, C, or O core (all $\mu_e=2$ ) shares the same limit to within a fraction of a percent. Why: the derivation pins the order-unity coefficient, delivering the canonical $1.4\,M_\odot$ that anchors white-dwarf and supernova astrophysics.

Limiting cases

How the star behaves depends on how relativistic its electrons are, set by its mass and central density.

Case	Result	Meaning
Non-relativistic degeneracy	$P\propto\rho^{5/3}$ , $\ R\propto M^{-1/3}$	Low-mass white dwarfs; $n=3/2$ polytrope. A stable equilibrium exists at any mass, and heavier stars are smaller.
Ultra-relativistic degeneracy	$P\propto\rho^{4/3}$	High-mass limit; $n=3$ polytrope. The mass is fixed at $M_{\mathrm{Ch}}$ independent of central density — the instability threshold.
Composition dependence	$M_{\mathrm{Ch}}=1.457\,(2/\mu_e)^2\,M_\odot$	Pure hydrogen ( $\mu_e=1$ ) would give $\approx5.8\,M_\odot$ , but real degenerate cores are He/C/O with $\mu_e\approx2$ , giving $1.4\,M_\odot$ .
Beyond the ideal model	$M_{\mathrm{crit}}<M_{\mathrm{Ch}}$	Inverse beta decay (neutronisation) and GR binding push $\Gamma_1$ below $4/3$ , destabilising the star slightly before the ideal limit.

Common trapThe softening is special relativity, not general relativity or heat

The most common confusion is about what causes the limit. The softening to $P\propto\rho^{4/3}$ is a purely kinematic effect of special relativity on the electrons' momentum–energy relation: once $p_F\gtrsim m_e c$ , the energy grows as $pc$ instead of $p^2/2m_e$ and the Fermi integral loses a power. It is not a thermal effect, and it is not general relativity. GR corrections and inverse beta decay do exist, but they only lower the critical mass by a few percent; the ideal $M_{\mathrm{Ch}}\approx1.4\,M_\odot$ is a special-relativistic result. A second trap: the composition factor is $\mu_e$ , the number of nucleons per electron, not the mean molecular weight $\mu$ . For fully ionised He, C, and O, $\mu_e=A/Z\approx2$ , so the limit is nearly composition-blind among realistic cores.

Worked examples

1 · Building the Chandrasekhar mass from fundamental constants (fully worked)

Evaluate $M_{\mathrm{Ch}}$ from equation (1) for a carbon–oxygen core ( $\mu_e=2$ ), then read off the hydrogen and iron cases. Use $\hbar=1.055\times10^{-34}$ , $c=2.998\times10^8$ , $G=6.674\times10^{-11}$ (SI), $m_H=1.661\times10^{-27}\ \mathrm{kg}$ , $\omega_3^0=2.018$ , $M_\odot=1.989\times10^{30}\ \mathrm{kg}$ .

Step 1. Form the quantum-gravity scale $\dfrac{\hbar c}{G}=\dfrac{(1.055\times10^{-34})(2.998\times10^8)}{6.674\times10^{-11}}=4.74\times10^{-16}\ \mathrm{kg^2}$ , so $(\hbar c/G)^{3/2}=1.03\times10^{-23}\ \mathrm{kg^3}$ — the Planck mass cubed, the gravitational scale that controls the answer.

Step 2. Form the composition denominator $(\mu_e m_H)^2=(2\times1.661\times10^{-27})^2=(3.322\times10^{-27})^2=1.104\times10^{-53}\ \mathrm{kg^2}$ .

Step 3. Assemble the prefactor $\dfrac{\sqrt{3\pi}}{2}\,\omega_3^0=1.549\times2.018=3.098$ .

Step 4. Multiply: $M_{\mathrm{Ch}}=3.098\times\dfrac{1.03\times10^{-23}}{1.104\times10^{-53}}=2.89\times10^{30}\ \mathrm{kg}$ . Dividing by $M_\odot$ gives $\boxed{M_{\mathrm{Ch}}=1.456\,M_\odot}$ — the canonical result from first principles.

Step 5. Composition scaling ( $M_{\mathrm{Ch}}\propto\mu_e^{-2}$ ): pure hydrogen ( $\mu_e=1$ ) quadruples the limit to $\mathbf{5.8\,M_\odot}$ — never realised, since hydrogen fuses first; an iron core ( $\mu_e=56/26=2.15$ ) drops it to $1.456\times(2/2.15)^2=\mathbf{1.26\,M_\odot}$ . Realistic compositions span only $\approx1.26$ – $1.46\,M_\odot$ .

2 · When do the electrons turn relativistic? Sirius B on the edge

Find the density at which electrons become relativistic ( $p_F=m_e c$ ), and compare it to the mean density of Sirius B ( $M=1.02\,M_\odot$ , $R\approx5800\ \mathrm{km}$ ) to judge how close it sits to the limit. Use $\mu_e=2$ .

Step 1. Set $p_F=m_e c$ in $n_e=\dfrac{8\pi}{3h^3}p_F^{\,3}$ . With $m_e c=2.73\times10^{-22}\ \mathrm{kg\,m/s}$ , this gives $n_e=\dfrac{8\pi}{3(6.626\times10^{-34})^3}(2.73\times10^{-22})^3=5.9\times10^{35}\ \mathrm{m^{-3}}$ .

Step 2. Convert to mass density $\rho=\mu_e m_H n_e=2\times1.661\times10^{-27}\times5.9\times10^{35}=1.96\times10^{9}\ \mathrm{kg/m^3}$ , i.e. $\boxed{\rho_{\text{rel}}\approx2\times10^{6}\ \mathrm{g/cm^3}}$ . Below this the star is unconditionally stable; above it the softening begins.

Step 3. Sirius B's mean density: $\bar\rho=\dfrac{M}{\frac{4}{3}\pi R^3}=\dfrac{1.02\times1.989\times10^{30}}{\frac{4}{3}\pi(5.8\times10^6)^3}=2.4\times10^{9}\ \mathrm{kg/m^3}\approx\mathbf{2.4\times10^{6}\ g/cm^3}$ , with a central density several times higher.

Step 4. Interpret. Sirius B's mean density already exceeds the relativistic-onset density, so its core electrons are mildly relativistic — not a low-mass, non-relativistic dwarf. Consistently, $M/M_{\mathrm{Ch}}\approx1.02/1.44\approx\mathbf{0.71}$ . This is why it sits on the steepening part of the mass–radius curve: a slightly heavier white dwarf would be markedly smaller, and one at $1.44\,M_\odot$ could not exist at all.

Physics of the portrait

The Chandrasekhar mass is inseparable from the young man who found it and the older man who refused to believe it — one of the sharpest confrontations in twentieth-century astrophysics.

Origins · 1 minSubrahmanyan Chandrasekhar and how this was discovered

Subrahmanyan Chandrasekhar

1910–1995 · Nobel Prize in Physics 1983 · degenerate matter and the structure of stars

In the summer of 1930, a nineteen-year-old Subrahmanyan Chandrasekhar passed the steamer voyage from Madras to Cambridge combining R. H. Fowler's 1926 quantum theory of white dwarfs with special relativity. At the enormous densities involved, he realised, the electrons must move near the speed of light, softening the pressure and implying a maximum mass above which no white dwarf can exist. He published the limit in 1931. The idea met fierce public opposition from Sir Arthur Eddington, the most eminent astrophysicist of the age, who in 1935 rose at a Royal Astronomical Society meeting to call the relativistic degeneracy formula a reductio ad absurdum — "I think 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 and drove Chandrasekhar toward other fields. History vindicated the calculation completely: the limit underlies white-dwarf demographics, Type Ia supernovae, and the very existence of neutron stars and black holes as the alternative endpoints. Chandrasekhar shared the 1983 Nobel Prize, and NASA's Chandra X-ray Observatory carries his name — a lasting reminder that authority is no substitute for a correct equation of state.

History

1926

— Fowler explains white dwarfs as matter supported by electron degeneracy pressure

1930

— Chandrasekhar, aged 19, works out the relativistic mass limit en route to England

1931

— Chandrasekhar publishes the maximum white-dwarf mass in the Astrophysical Journal

1935

— Eddington publicly ridicules the limit at the Royal Astronomical Society

1983

— Chandrasekhar shares the Nobel Prize in Physics for stellar-structure theory

Going further

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., §16.3–16.4 · Shapiro & Teukolsky, Black Holes, White Dwarfs, and Neutron Stars, Ch. 2–4 (degenerate EOS and the exact limit) · Chandrasekhar, An Introduction to the Study of Stellar Structure (1939) · Kippenhahn, Weigert & Weiss, Stellar Structure and Evolution, 2nd ed., Ch. 15 · Padmanabhan, Theoretical Astrophysics, Vol. II, §5.

Provenance. The coefficient $M_{\mathrm{Ch}}=1.457\,(2/\mu_e)^2\,M_\odot$ and the $n=3$ Lane–Emden constants ( $\xi_1=6.89685$ , $(-\xi^2\theta')_{\xi_1}=2.01824$ ) follow Chandrasekhar (1931, ApJ 74, 81) and Carroll & Ostlie §16.4; the pressure coefficients in equations (2)–(3) are the standard Fermi-integral results in Shapiro & Teukolsky Ch. 2. A verifier should reconstruct $M_{\mathrm{Ch}}$ from $(\sqrt{3\pi}/2)\,\omega_3^0\,(\hbar c/G)^{3/2}/(\mu_e m_H)^2$ with $\mu_e=2$ and $m_H$ the atomic mass unit, confirming $1.46\,M_\odot$ ; the onset density $\rho_{\text{rel}}\approx2\times10^{6}\ \mathrm{g/cm^3}$ follows from $p_F=m_e c$ .

Prerequisites. Electron Degeneracy Pressure · White Dwarfs · Equation of State of Dense Matter. Leads to. Neutron Stars · The Tolman–Oppenheimer–Volkoff Equation · White Dwarf Cooling and Crystallisation.

Ask the tutor: Ask me to walk through why the electron mass cancels out of the ultra-relativistic pressure, or to show step by step how the same $\rho^{4/3}$ scaling turns Type Ia supernovae into standard candles.

Is the Chandrasekhar mass exactly 1.4 solar masses?

Not exactly, and the digits matter. The ideal special-relativistic value for a $\mu_e=2$ core is $1.457\,M_\odot$ ; the often-quoted " $1.4\,M_\odot$ " is a round number, and " $1.44\,M_\odot$ " folds in small corrections. Real progenitors are lowered a few percent by inverse beta decay and general-relativistic binding, and shifted by iron composition ( $\mu_e\approx2.15$ gives $1.26\,M_\odot$ ); temperature, Coulomb corrections, and rotation move it too. So $1.4\,M_\odot$ is the right number to remember, but the effective limit deciding a star's fate can differ by ten percent or more.

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

This backwards relation, $R\propto M^{-1/3}$ , is the fingerprint of degeneracy pressure and is genuinely counter-intuitive — for a rock, adding mass adds volume. Here, adding mass increases gravity, which compresses the star until the electrons are packed densely enough to push back; denser packing means a smaller star. The relation comes from balancing the non-relativistic pressure (equation 2) against gravity in a polytrope, and it explains why a Chandrasekhar-mass dwarf would formally have zero radius: no density is high enough for the softened relativistic pressure to hold gravity off.

What actually happens if you push a white dwarf past M_Ch?

It cannot become a slightly larger white dwarf — no cold degenerate equilibrium exists above the limit. One of two things happens. If the star is a carbon–oxygen dwarf accreting from a companion, the rising central density ignites carbon fusion before full collapse, and the thermonuclear runaway unbinds the entire star as a Type Ia supernova, leaving no remnant. If instead the core is already iron and cannot release fusion energy, it collapses: electrons are captured onto protons, pressure support vanishes, and the collapse halts only when nuclear-density neutron degeneracy takes over — forming a neutron star, or a black hole if even that fails.

Why does μ_e appear instead of the ordinary mean molecular weight μ?

Because the pressure holding the star up comes entirely from the electrons, not the ions. What matters is the electron number density $n_e=\rho/(\mu_e m_H)$ , where $\mu_e$ is the mass of gas per electron in nucleon masses. For a fully ionised nucleus with mass number $A$ and charge $Z$ , each supplies $Z$ electrons for $A$ nucleons, so $\mu_e=A/Z$ . Helium, carbon, and oxygen all have $A/Z=2$ , which is why the limit is nearly the same for any realistic core. Hydrogen would give $\mu_e=1$ and a much larger limit, but it fuses long before becoming degenerate on this scale.

How can a mass limit be built almost entirely from fundamental constants?

It is one of the most beautiful results in stellar physics. The ultra-relativistic pressure (equation 3) contains only $\hbar$ , $c$ , and the number density — the electron mass has cancelled out. Balancing that against gravity, which brings in $G$ and the nucleon mass, leaves $M_{\mathrm{Ch}}\sim(\hbar c/G)^{3/2}/(\mu_e m_H)^2$ . Since $(\hbar c/G)^{1/2}$ is the Planck mass, $M_{\mathrm{Ch}}$ is essentially the Planck mass cubed over the nucleon mass squared — a solar-scale mass conjured from the constants of quantum mechanics, relativity, and gravity, with composition the only tunable input.

Is the neutron-star maximum mass just the Chandrasekhar limit again?

The logic is the same but the numbers and physics differ. A neutron star is held up by neutron degeneracy pressure, so by the same argument it has a maximum mass. Naively rescaling the formula suggests a few solar masses, but the answer is dominated by two effects the white-dwarf case can ignore: general relativity is essential (the Tolman–Oppenheimer–Volkoff equation replaces Newtonian equilibrium), and the nuclear equation of state above saturation density is genuinely uncertain. The observed maximum lies near $2$ – $2.3\,M_\odot$ , and pinning it down is a major goal of gravitational-wave astronomy.

Does the white dwarf's temperature affect the limit?

Barely — and that is the point of degeneracy. The pressure is set by the Pauli-exclusion filling of momentum states, which persists even at zero temperature; thermal energy is a small correction on top of the huge Fermi energy. So unlike a main-sequence star, whose support depends on ongoing fusion, a white dwarf's structure and its limit are essentially temperature-independent — which is why white dwarfs can cool for billions of years without shrinking.