White Dwarf Stars and the Chandrasekhar Limit

A white dwarf is what most stars become: the cold, spent core left when a star below about $8\,M_\odot$ exhausts its nuclear fuel. Gravity crushes the core inward — and then something purely quantum-mechanical stops it. The electrons, forbidden by the Pauli exclusion principle from sharing states, are forced to high momenta and push back with a pressure that survives even at absolute zero. The result is an object with the mass of the Sun packed into the volume of the Earth, held up not by heat but by the statistical mechanics of a degenerate Fermi gas. No cleaner astrophysical laboratory for that gas exists. The whole drama is a contest of two scalings. Gravity's binding energy grows as $-GM^2/R$ . In the ordinary, non-relativistic regime the degenerate electrons supply a kinetic energy $\propto +\hbar^2 N^{5/3}/(m_e R^2)$ — falling off more slowly than gravity as the star shrinks, so an equilibrium radius always exists and the star is stable. But squeeze harder and the Fermi momentum climbs past $m_e c$ ; the electron energy crosses over to the ultra-relativistic form $\propto +\hbar c\,N^{4/3}/R$ , which shares the same $1/R$ dependence as gravity. Now stability no longer depends on radius at all — only on the coefficient of $1/R$ , that is, on the total mass. Chandrasekhar's 1931 insight was that this relativistic softening implies a hard ceiling: above a critical mass, gravity's coefficient wins at every radius and no equilibrium can exist. Solving hydrostatic equilibrium with the relativistic equation of state (the $n=3$ Lane–Emden polytrope) gives $M_{\mathrm{Ch}}\approx1.44\,M_\odot$ for a helium/carbon/oxygen composition. Remarkably, the answer is built only from $\hbar$ , $c$ , $G$ and the proton mass — the electron mass cancels. It is a quantum-statistical bound on how much gravitating matter one star can hold, and the trigger for neutron-star formation and Type Ia supernovae.

Key takeaways

A white dwarf is the cold remnant of a low-mass star, held up not by thermal pressure but by the zero-temperature degeneracy pressure of its electrons, forced to high momenta by the Pauli exclusion principle.
In the non-relativistic regime $P\propto\rho^{5/3}$ gives stable stars of finite radius, with the inverted mass–radius law $R\propto M^{-1/3}$ : heavier white dwarfs are smaller.
As the mass grows the central electrons become ultra-relativistic and the equation of state softens to $P\propto\rho^{4/3}$ — the marginal $\Gamma=4/3$ case, where pressure and gravity scale identically with radius.
The $\rho^{4/3}$ polytrope has a unique mass independent of central density: the Chandrasekhar mass $M_{\mathrm{Ch}}\approx1.44\,M_\odot$ (for $\mu_e=2$ ), above which no degenerate equilibrium exists.
The limit is $M_{\mathrm{Ch}}=\tfrac{\sqrt{3\pi}}{2}\omega_3(\hbar c/G)^{3/2}(\mu_e m_p)^{-2}$ — quantum mechanics, relativity and gravity in one number; the electron mass drops out.
Above $M_{\mathrm{Ch}}$ : collapse to a neutron star or black hole; near it, accreting white dwarfs ignite as Type Ia supernovae, the standard candles that revealed cosmic acceleration.
Key idea: degeneracy pressure scales as $P\propto\rho^{5/3}$ (non-relativistic) but softens to $\rho^{4/3}$ once electrons turn relativistic; because the $\rho^{4/3}$ pressure and gravity both scale as the same power of radius, equilibrium stops depending on radius and fixes a single maximum mass $M_{\mathrm{Ch}}\approx1.44\,M_\odot$ .

The core insight

Degeneracy pressure scales with density as $P\propto\rho^{5/3}$ (non-relativistic) but softens to $P\propto\rho^{4/3}$ when the electrons become ultra-relativistic; since the $\rho^{4/3}$ pressure and gravity both scale as the same power of radius, equilibrium ceases to depend on radius and instead fixes a single maximum mass $M_{\mathrm{Ch}}\approx1.44\,M_\odot$ .

Where this lives today

The Chandrasekhar mass is one of the most-used numbers in modern astrophysics. A carbon–oxygen white dwarf accreting from a companion (or merging with another) is pushed toward $M_{\mathrm{Ch}}$ ; because degenerate pressure is temperature-independent, the carbon fusion that ignites near the limit runs away instead of self-regulating, detonating the whole star. Since this happens at a nearly universal mass, Type Ia supernovae have a standardizable peak luminosity — they are the standard candles that revealed the accelerating expansion of the universe (1998 Nobel-winning discovery). The precise masses and radii of nearby white dwarfs (Sirius B, 40 Eridani B) from astrometry and gravitational redshift confirm the $R\propto M^{-1/3}$ law, and no white dwarf above $M_{\mathrm{Ch}}$ has ever been found in equilibrium — a clean empirical test of a Pauli-driven bound written in fundamental constants.

Origins · 1 minChandrasekhar limit and how this was discovered

Chandrasekhar limit

maximum mass of a white dwarf, ≈1.44 M⊙ · the quantum-statistical ceiling on gravitating matter

The Chandrasekhar limit is the maximum mass a white dwarf can have while remaining supported by electron degeneracy pressure. Computed by Subrahmanyan Chandrasekhar in 1930 — famously during his voyage from Madras to Cambridge, aged nineteen — it emerges from the relativistic softening of the electron equation of state to $\Gamma=4/3$ . Its value $M_{\mathrm{Ch}}=\tfrac{\sqrt{3\pi}}{2}\omega_3(\hbar c/G)^{3/2}(\mu_e m_p)^{-2}\approx1.44\,M_\odot$ combines quantum mechanics, special relativity and gravity into a single stellar mass, and it marks the boundary between white dwarfs and the more compact remnants — neutron stars and black holes.

History

1926

— R. H. Fowler explains white dwarfs as degenerate electron gas, applying the brand-new Fermi–Dirac statistics.

1930

— Subrahmanyan Chandrasekhar, aged nineteen, computes the relativistic maximum mass on his voyage to Cambridge.

1930

— E. C. Stoner independently derives a relativistic limiting mass (the name Chandrasekhar's, the priority partly Stoner's).

1931

— Chandrasekhar publishes 'The Maximum Mass of Ideal White Dwarfs' in the Astrophysical Journal.

1935

— Eddington publicly ridicules the limit at the Royal Astronomical Society; he is wrong, and the limit stands.

1983

— Chandrasekhar shares the Nobel Prize in Physics for the theory of stellar structure and evolution.

The physical picture

Squeeze electrons into a smaller volume and exclusion forces them into higher-momentum states; the outward push that results exists even at $T=0$ . It is the macroscopic face of fermions resist crowding, and a white dwarf is that resistance frozen into the structure of a star. The pressure law was worked out in Degeneracy Pressure of Fermions, which owns the result $P_0=\tfrac25 n\varepsilon_F$ and the $n^{5/3}\to n^{4/3}$ softening; here we take that pressure and let gravity be the wall. Why a maximum mass? Think of the total energy as a function of radius. Gravity contributes $-A/R$ . The non-relativistic kinetic energy contributes $+B/R^2$ , and the sum $-A/R+B/R^2$ has a genuine minimum — a stable star at finite $R$ . But the relativistic kinetic energy contributes $+B'/R$ , the same shape as gravity. Now the total is $(B'-A)/R$ : if $B'>A$ the star simply expands back to the non-relativistic branch; if $B'<A$ it plunges to $R\to0$ with no minimum at all. The borderline $A=B'$ has no $R$ in it — it is a condition on mass alone. That is the Chandrasekhar mass. The hero figure below is the mass–radius relation. On the low-mass branch $R\propto M^{-1/3}$ : pile on mass and the star shrinks. Follow the curve up and the radius plunges toward zero at a single mass — the vertical wall at $M_{\mathrm{Ch}}$ . No white dwarf can exist to the right of that wall.

Figure

Mass–radius relation of white dwarfs — radius falls as R ∝ M^(−1/3) then plunges to zero at the Chandrasekhar wall M_Ch ≈ 1.44 M_⊙

Strikingly, the electron mass drops out of the final answer. On the branch that sets the limit the electrons are ultra-relativistic, $\varepsilon\approx pc$ , with no $m_e$ in sight; the pressure $P\propto\hbar c\,n_e^{4/3}$ carries only $\hbar$ , $c$ and density. The balance against gravity then involves only $\hbar$ , $c$ , $G$ and the mass carried per electron, $\mu_e m_p$ . The electron mass sets where relativity begins, not how much mass a star can hold.

Figure

Energy versus radius — the non-relativistic gas has a stable minimum, but the ultra-relativistic gas shares gravity's 1/R and has no minimum above M_Ch

The core equations

Everything starts with counting electrons in the filled Fermi sphere. At $T=0$ every momentum state up to $p_F$ is occupied once per spin, so the number density fixes the Fermi momentum. This is the same phase-space count used throughout the degenerate-gas articles.

n_e=\frac{p_F^3}{3\pi^2\hbar^3},\qquad p_F=\hbar(3\pi^2 n_e)^{1/3},\qquad \rho=\mu_e m_p\, n_e \qquad (1)

where: $n_e$ — electron number density; $p_F$ — Fermi momentum; $\rho$ — mass density (carried by the ions, not the electrons); $\mu_e$ — number of nucleons per electron ( $=2$ for He/C/O); $m_p$ — nucleon mass. The electrons supply all the pressure but almost none of the mass — the split that makes $\mu_e$ appear. The pressure of the gas comes from the momentum flux of the fast electrons at the Fermi surface. In full generality the $T=0$ pressure is a single integral over the filled sphere with the relativistic velocity $v=pc^2/\varepsilon$ .

P=\frac{1}{3\pi^2\hbar^3}\int_0^{p_F}\frac{p^4 c^2}{\sqrt{p^2c^2+m_e^2c^4}}\,\mathrm{d}p \qquad (2)

where: $P$ — electron degeneracy pressure; the factor $\tfrac13$ inside is the isotropy average of $p_x v_x$ ; the square root is the relativistic energy $\varepsilon(p)$ . Read it as one-third the flux of momentum carried by electrons streaming across any surface — a purely kinematic pressure that does not vanish at $T=0$ . Evaluating this integral in the two limits $p_F\ll m_e c$ and $p_F\gg m_e c$ gives the two polytropes that decide the star's fate. The stiff one supports any mass; the soft one does not.

P_{\mathrm{NR}}=\frac{(3\pi^2)^{2/3}}{5}\frac{\hbar^2}{m_e}n_e^{5/3}\ \ (\Gamma=\tfrac53),\qquad P_{\mathrm{UR}}=\frac{(3\pi^2)^{1/3}}{4}\hbar c\,n_e^{4/3}\ \ (\Gamma=\tfrac43) \qquad (3)

where: $\Gamma$ — polytropic exponent, $P\propto\rho^\Gamma$ . The non-relativistic $\tfrac53$ is stiff — compress the gas and it pushes back hard; the ultra-relativistic $\tfrac43$ is soft — the same compression buys much less pressure, because a fast electron gains only $\mathrm{d}\varepsilon\approx c\,\mathrm{d}p$ from extra momentum, not $p\,\mathrm{d}p/m_e$ . The star's structure follows from hydrostatic equilibrium: at each radius the pressure gradient must support the weight of the material above.

\frac{\mathrm{d}P}{\mathrm{d}r}=-\frac{G\,m(r)\,\rho(r)}{r^2},\qquad \frac{\mathrm{d}m}{\mathrm{d}r}=4\pi r^2\rho \qquad (4)

where: $r$ — radial coordinate; $m(r)$ — mass enclosed within $r$ ; $G$ — Newton's constant. With a polytropic $P(\rho)=K\rho^\Gamma$ these two equations collapse into the single dimensionless Lane–Emden equation, whose $n=3$ solution delivers the limit.

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

where: $M_{\mathrm{Ch}}$ — limiting white-dwarf mass; $\omega_3=-\xi^2\theta'|_{\xi_1}\approx2.018$ — the $n=3$ Lane–Emden eigenvalue; $M_\odot=1.989\times10^{30}$ kg. Every constant in it is fundamental except the composition $\mu_e$ ; the electron mass is gone. For $\mu_e=2$ this is $\approx1.44\,M_\odot$ .

Limiting cases

The equation of state has two universal regimes set by the Fermi momentum relative to $m_e c$ :

Case	Result	Meaning
Non-relativistic ( $p_F\ll m_e c$ )	$P\propto\rho^{5/3}$ , polytrope $n=\tfrac32$	Radius decreases with mass, $R\propto M^{-1/3}$ . Always stable; low-mass, low-density white dwarfs.
Ultra-relativistic ( $p_F\gg m_e c$ )	$P\propto\rho^{4/3}$ , polytrope $n=3$	Pressure and gravity scale identically with $R$ ; mass becomes independent of central density and approaches $M_{\mathrm{Ch}}$ . Marginal stability.
GR + neutronization corrections	true max $<M_{\mathrm{Ch}}$	Before the ideal limit, electron capture (neutronization) and general-relativistic instability soften the gas further and begin collapse slightly below $1.44\,M_\odot$ .

Derivation AExistence of a maximum mass from an energy argument

The goal is to see, with a back-of-envelope variational estimate, why the ultra-relativistic gas cannot support arbitrary mass and to locate the limit up to order-unity factors. We model the star as $N$ electrons (and $\mu_e N$ nucleons) in a uniform sphere of radius $R$ , write its total energy as kinetic plus gravitational, and study how the two terms compete as $R$ changes. The whole point is which power of $R$ each term carries.

Step 1. Fermi momentum from confinement. Packing $N$ electrons into volume $\sim R^3$ gives $n_e\sim N/R^3$ , so $p_F\sim\hbar n_e^{1/3}\sim\hbar N^{1/3}/R$ . Why: exclusion sets a momentum scale from the inter-particle spacing $n_e^{-1/3}$ — the tighter the box, the faster the electrons.

Step 2. Ultra-relativistic kinetic energy. Each fast electron carries $\varepsilon\approx p_F c$ , so $E_{\mathrm{kin}}\sim N p_F c\sim \hbar c\,N^{4/3}/R$ . Why: with the rest mass negligible the energy is linear in momentum, and $p_F\propto 1/R$ drags the total into a $1/R$ law — the same shape as gravity.

Step 3. Gravitational energy. With $M\approx\mu_e m_p N$ (nucleons dominate the mass), $E_{\mathrm{grav}}\sim-G M^2/R\sim -G(\mu_e m_p)^2 N^2/R$ . Why: the self-binding of a uniform sphere is $\propto M^2/R$ ; both energies now live on the same $1/R$ curve.

Step 4. The instability condition. The total is $E\sim\tfrac1R\big[\hbar c\,N^{4/3}-G(\mu_e m_p)^2 N^2\big]$ . If the bracket is positive, $E\to0^+$ as $R\to\infty$ and the star drifts back to the stable non-relativistic branch; if negative, $E\to-\infty$ as $R\to0$ — runaway collapse. Why: with no competing power of $R$ , only the sign of the bracket decides the fate — radius has dropped out of the question entirely.

Step 5. The critical mass. Setting the bracket to zero, $N_{\max}\sim(\hbar c/G(\mu_e m_p)^2)^{3/2}$ , hence $M_{\max}\sim\mu_e m_p N_{\max}\sim(\hbar c/G)^{3/2}/(\mu_e m_p)^2$ . Why: this is the exact scaling of the Chandrasekhar mass; only the dimensionless Lane–Emden factor $\tfrac{\sqrt{3\pi}}{2}\omega_3\approx3.1$ separates the estimate from the precise number. *Verified: the scaling $M_{\max}\propto(\hbar c/G)^{3/2}(\mu_e m_p)^{-2}$ and the electron mass dropping out are checked symbolically, and the numerical prefactor and $1.44\,M_\odot$ result in verify_white-dwarf-chandrasekhar-limit.py.*

Go deeper · Derivation B · 1 minThe non-relativistic mass–radius relation R ∝ M to the minus one third

Derivation BThe non-relativistic mass–radius relation R ∝ M to the minus one third

The goal is the famous inverted mass–radius law for low-mass white dwarfs, where the electrons are non-relativistic ( $\Gamma=\tfrac53$ , polytrope $n=\tfrac32$ ). We use the standard polytrope relations that express radius and mass as powers of the central density $\rho_c$ , then eliminate $\rho_c$ between them. Watch how the exponents conspire to give a negative power of mass — the counter-intuitive signature of degeneracy support.

Step 1. Polytrope radius. For index $n$ , $R\propto\rho_c^{(1-n)/(2n)}$ . With $n=\tfrac32$ the exponent is $(1-\tfrac32)/3=-\tfrac16$ , so $R\propto\rho_c^{-1/6}$ . Why: a denser centre means a more compact star; the mild $-\tfrac16$ power comes straight from the Lane–Emden scaling.

Step 2. Polytrope mass. Likewise $M\propto\rho_c^{(3-n)/(2n)}$ ; with $n=\tfrac32$ the exponent is $(3-\tfrac32)/3=\tfrac12$ , so $\rho_c\propto M^2$ . Why: heavier stars are built from denser cores — the mass grows as the square root of central density, so $\rho_c$ grows as $M^2$ .

Step 3. Eliminate the central density. Substituting, $R\propto\rho_c^{-1/6}\propto(M^2)^{-1/6}=M^{-1/3}$ . Why: the two power laws combine so that more massive white dwarfs are physically smaller — the exact opposite of an ordinary self-gravitating gas sphere, and a direct, observable fingerprint of degeneracy.

Step 4. Approach to the limit. As $M\to M_{\mathrm{Ch}}$ the central density diverges and $R\to0$ ; in reality the relativistic softening ( $n\to3$ ) takes over and drives $R\to0$ precisely at $M_{\mathrm{Ch}}$ . Why: the non-relativistic law already points toward a vanishing radius; relativity supplies the hard wall where it actually happens. *Verified: the $R\propto M^{-1/3}$ and $\rho_c\propto M^2$ scalings are confirmed symbolically from the polytrope exponents in verify_white-dwarf-chandrasekhar-limit.py; Sirius B ( $M\approx1.0\,M_\odot$ , $R\approx0.008\,R_\odot$ ) sits on this relation.*

Go deeper · Derivation C · 1 minEvaluating the Chandrasekhar mass numerically

Derivation CEvaluating the Chandrasekhar mass numerically

The goal is to turn equation (5) into an actual number, $1.44\,M_\odot$ , by inserting the fundamental constants and the Lane–Emden eigenvalue. This is exactly the step where the AllQFT base article slipped a numerical error, so we do the grouping carefully and check each factor. The strategy: assemble the constant combination $(\hbar c/G)^{3/2}/m_p^2$ , multiply by the dimensionless prefactor, and read the result in solar masses.

Step 1. Group the fundamental constants — and fix the base error. Compute $\hbar c/G=(1.0546\times10^{-34})(2.998\times10^{8})/(6.674\times10^{-11})=4.74\times10^{-16}\ \mathrm{kg^2}$ . Why: this is a mass-squared; its square root $\sqrt{\hbar c/G}=2.18\times10^{-8}$ kg is the Planck mass. The base article wrongly wrote $\hbar c/G\approx3.16\times10^{-3}\ \mathrm{kg^2}$ — off by thirteen orders of magnitude; had you used it, $M_{\mathrm{Ch}}$ would come out as $\sim10^{19}\,M_\odot$ . Always sanity-check this combination against the Planck mass.

Step 2. Raise to the $3/2$ power and divide by $m_p^2$ . $(\hbar c/G)^{3/2}=1.031\times10^{-23}\ \mathrm{kg^3}$ , and dividing by $m_p^2=2.80\times10^{-54}\ \mathrm{kg^2}$ gives $3.69\times10^{30}\ \mathrm{kg}=1.85\,M_\odot$ . Why: this is the raw Planck-mass-cubed over nucleon-mass-squared scale — already of order a solar mass, before any structure factor.

Step 3. Apply the dimensionless prefactor. $\tfrac{\sqrt{3\pi}}{2}\omega_3=\tfrac12(3.070)(2.018)=3.10$ . Why: $\sqrt{3\pi}=3.070$ ; the $n=3$ Lane–Emden eigenvalue $\omega_3\approx2.018$ carries the actual stellar structure — the density profile is not uniform, and this factor accounts for it.

Step 4. Collect and set the composition. Combining all factors gives $M_{\mathrm{Ch}}=5.83\,M_\odot/\mu_e^2$ ; with $\mu_e=2$ (hydrogen-free He/C/O), $M_{\mathrm{Ch}}=1.44\,M_\odot$ . Why: one electron per two nucleons gives $\mu_e=2$ ; small Coulomb, GR and electron-capture corrections nudge the true value slightly below, to the canonical $1.4\,M_\odot$ . *Verified: $\hbar c/G=4.74\times10^{-16}\ \mathrm{kg^2}$ (NOT the base's $3.16\times10^{-3}$ ), the prefactor $3.10$ , the coefficient $5.83\,M_\odot/\mu_e^2$ , and $M_{\mathrm{Ch}}=1.435\,M_\odot$ at $\mu_e=2$ are all recomputed in verify_white-dwarf-chandrasekhar-limit.py.*

Common trapThe Chandrasekhar limit is where the central electrons 'become relativistic'

Not quite — and the distinction matters. The limit is not the density at which the central electrons cross $p_F=m_e c$ ; that crossover happens at a modest $\rho\sim2\times10^6\ \mathrm{g\,cm^{-3}}$ , well inside many stable white dwarfs. The limit is the asymptotic mass of the $n=3$ polytrope, reached as the central density tends to infinity. Two more traps: the electrons provide the pressure but the ions provide the mass, which is why $\mu_e$ (nucleons per electron) appears squared; and the star is not literally at $T=0$ — it still shines — but $T\ll T_F$ makes thermal pressure negligible to better than a percent.

Worked examples

Example 1 — The crossover density where relativity switches on.

Step 1. Electrons turn relativistic when $p_F\approx m_e c=(9.109\times10^{-31})(2.998\times10^{8})=2.73\times10^{-22}\ \mathrm{kg\,m/s}$ . From $p_F=\hbar(3\pi^2 n_e)^{1/3}$ , invert: $3\pi^2 n_e=(p_F/\hbar)^3=(2.59\times10^{12})^3=1.74\times10^{37}\ \mathrm{m^{-3}}$ , so $n_e=\mathbf{5.9\times10^{29}\ cm^{-3}}$ (i.e. $5.9\times10^{35}\ \mathrm{m^{-3}}$ ).

Step 2. Convert to mass density with $\mu_e=2$ : $\rho=\mu_e m_p n_e=2\times(1.673\times10^{-27})\times(5.9\times10^{35})=\mathbf{2.0\times10^{9}\ kg\,m^{-3}}=2\times10^6\ \mathrm{g\,cm^{-3}}$ — roughly two tonnes per cubic centimetre, a typical white-dwarf core density.

Step 3. Below this the electrons are non-relativistic ( $P\propto\rho^{5/3}$ ); above it the law softens toward $P\propto\rho^{4/3}$ . Massive white dwarfs, being denser, live on the soft branch — which is exactly why they approach the mass limit. The most massive observed white dwarfs sit at $\rho\sim10^9$ – $10^{10}\ \mathrm{g\,cm^{-3}}$ , firmly ultra-relativistic.

Example 2 — The Chandrasekhar mass, digit by digit, with the base error avoided.

Step 1. Form $\hbar c/G=(1.0546\times10^{-34})(2.998\times10^{8})/(6.674\times10^{-11})=\mathbf{4.74\times10^{-16}\ kg^2}$ . Sanity check: $\sqrt{4.74\times10^{-16}}=2.18\times10^{-8}$ kg $=$ Planck mass. ✓ (The base's $3.16\times10^{-3}$ fails this check instantly.)

Step 2. Raise and divide: $(\hbar c/G)^{3/2}=1.031\times10^{-23}\ \mathrm{kg^3}$ ; with $m_p^2=(1.673\times10^{-27})^2=2.80\times10^{-54}\ \mathrm{kg^2}$ , the ratio is $3.69\times10^{30}\ \mathrm{kg}=\mathbf{1.85\,M_\odot}$ .

Step 3. Multiply by $\tfrac{\sqrt{3\pi}}{2}\omega_3=3.10$ and divide by $\mu_e^2=4$ : $M_{\mathrm{Ch}}=1.85\times3.10/4\,M_\odot=\mathbf{1.43\,M_\odot}$ . The standard quoted value is $1.4$ – $1.44\,M_\odot$ ; small physical corrections lower it slightly. Note the result used $\hbar$ , $c$ , $G$ , $m_p$ and $\mu_e$ only — no electron mass anywhere.

Going further

References. Pathria & Beale, Statistical Mechanics, 3rd ed., §8.5 (Theory of white dwarf stars) · Chandrasekhar, An Introduction to the Study of Stellar Structure (1939), Ch. XI · Kardar, Statistical Physics of Particles, problems on the degenerate Fermi gas · Landau & Lifshitz, Statistical Physics Part 1, §106 · Shapiro & Teukolsky, Black Holes, White Dwarfs and Neutron Stars (1983), Ch. 3–4 · Chandrasekhar, 'The Maximum Mass of Ideal White Dwarfs', Astrophys. J. 74, 81 (1931).

Provenance. Equation of state, hydrostatic equilibrium and Lane–Emden treatment follow Pathria & Beale §8.5, Chandrasekhar (1939) Ch. XI, and Shapiro & Teukolsky Ch. 3–4. This article recalls the pressure law $P_0=\tfrac25 n\varepsilon_F$ and the $n^{5/3}\to n^{4/3}$ softening from Degeneracy Pressure of Fermions by name, and the relativistic equation of state from Relativistic Degenerate Fermi Gas, and owns the mass–radius relation and the Chandrasekhar mass. Base-error note: the AllQFT base article's numerical Chandrasekhar derivation stated $\hbar c/G\approx3.16\times10^{-3}\ \mathrm{kg^2}$ , wrong by $\sim10^{13}$ ; the correct value is $4.74\times10^{-16}\ \mathrm{kg^2}$ (square root $=$ Planck mass), fixed here and machine-verified. The final $1.44\,M_\odot$ is unaffected because it uses the correct grouping.

Builds on: Degenerate Fermi Gas at Zero Temperature · Degeneracy Pressure of Fermions · Fermi Energy, Momentum and Fermi Temperature · Density of States for Free Particles Leads to: Neutron Stars and the Tolman–Oppenheimer–Volkoff Limit · Type Ia Supernovae as Standard Candles · Polytropes and the Lane–Emden Equation

If the pressure is finite at T = 0, how can the white dwarf still cool and shine?

Two separate things are happening. The mechanical support is the $T=0$ degeneracy pressure, which is essentially temperature-independent because $T\ll T_F$ . The luminosity comes from the leftover thermal energy stored in the ions (and a thin non-degenerate surface layer), which leaks out slowly over billions of years. The star cools without changing its structure appreciably — the ions lose heat while the electron gas keeps holding the star up. So 'cold' means $T\ll T_F$ , not literally zero, and the two processes are almost decoupled.

Why does μ_e (nucleons per electron) enter the Chandrasekhar mass squared?

It appears once from the pressure and once from the mass. The pressure is set by the electron density $n_e$ , but the mass density is $\rho=\mu_e m_p n_e$ — so for a given mass density you get fewer electrons (hence less pressure) when $\mu_e$ is larger. Tracing this through the balance of $P\propto n_e^{4/3}$ against gravity $\propto\rho^{4/3}=(\mu_e m_p n_e)^{4/3}$ , plus the mass carried per electron, leaves $M_{\mathrm{Ch}}\propto\mu_e^{-2}$ . Physically: more nucleons per electron means less quantum pressure per unit weight, so the star can hold less mass.

Is the Chandrasekhar mass exactly 1.44 M⊙, or does it depend on the star?

The ideal-gas value is $5.83\,M_\odot/\mu_e^2$ , so it depends on composition through $\mu_e$ : for hydrogen-free matter (He, C, O, or even $^{56}$ Fe) $\mu_e=2$ and you get $1.44\,M_\odot$ ; pure hydrogen ( $\mu_e=1$ ) would give $\sim5.8\,M_\odot$ , but hydrogen-rich cores never form white dwarfs. Real white dwarfs cluster near $1.4\,M_\odot$ because they are C/O remnants. Small corrections — Coulomb interactions between ions and electrons, general relativity, and electron capture (neutronization) — all lower the practical limit a little below $1.44\,M_\odot$ .

How does this relate to a neutron star — do they have the same kind of limit?

Yes, structurally. Replace the electrons by neutrons (with the mass now carried by the neutrons themselves) and repeat the identical energy-balance argument; you again find a maximum mass where the ultra-relativistic $\Gamma=4/3$ pressure loses to gravity. The scale is set by $(\hbar c/G)^{3/2}/m_n^2$ instead of $/(\mu_e m_p)^2$ . But neutron stars are far more compact, so general relativity is essential and the ideal-gas estimate is only a guide — the details are the subject of Neutron Stars and Nuclear Degeneracy.

Why does a Type Ia supernova have a standard brightness?

Because it detonates at a nearly universal mass. An accreting C/O white dwarf approaches $M_{\mathrm{Ch}}$ ; near the limit the central density and temperature rise until carbon ignites. Crucially, the pressure is degenerate and temperature-independent, so the burning does not expand and cool the star to throttle itself — it runs away and consumes the whole $\sim1.4\,M_\odot$ . Since the fuel mass is always about the same, the energy released and hence the peak luminosity is standardizable, letting astronomers use Type Ia supernovae as cosmic distance markers.

The base article said ħc/G ≈ 3.16×10⁻³ kg². How would I have caught that mistake?

By dimensional and order-of-magnitude sanity checks. First, $\hbar c/G$ has units of mass squared, and its square root should be the Planck mass, $\sqrt{\hbar c/G}=2.18\times10^{-8}$ kg — a famous number. If you square that you get $4.74\times10^{-16}\ \mathrm{kg^2}$ , not $3.16\times10^{-3}$ . Second, propagate the wrong value: $(3.16\times10^{-3})^{3/2}/m_p^2\sim10^{49}\ \mathrm{kg}\sim10^{19}\,M_\odot$ , an absurd stellar mass. Whenever a combination of fundamental constants appears, check it against a known landmark (here the Planck mass) before trusting downstream arithmetic.