Neutron Stars and Nuclear Degeneracy

A neutron star is the collapsed core left behind when a massive star ends its life in a supernova. Once the density climbs high enough, electron capture $p+e^-\to n+\nu_e$ converts nearly all the protons and electrons into neutrons, and what remains is a colossal Fermi sea of neutrons squeezed to nuclear density. Because the neutron Fermi temperature $T_F\sim10^{12}\,$ K towers over the interior temperature $\sim10^{8}\,$ K, the gas is, to superb approximation, a fully degenerate ideal Fermi gas at $T=0$ . The pressure it exerts — degeneracy pressure, a direct consequence of the Pauli exclusion principle — is all that holds the star up against its own crushing weight. The mechanics is the same contest as in a white dwarf, one mass scale up. Gravity lowers the energy as $-GM^2/R$ . The degenerate neutrons resist: confining $N$ fermions to radius $R$ forces them to a Fermi momentum $p_F\propto N^{1/3}/R$ , and in the non-relativistic regime this kinetic energy scales as $1/R^2$ . Since gravity only goes as $1/R$ , the total energy has a genuine minimum — a stable star with the inverted mass–radius law $R\propto M^{-1/3}$ : heavier neutron stars are smaller. The whole calculation is the self-gravitating twin of the white dwarf, with the neutron replacing the electron. The drama again lives in relativity. As the mass grows, $p_F$ rises until $p_F c\gtrsim m_n c^2$ and the neutrons turn ultra-relativistic; their kinetic energy then scales as $1/R$ , identical to gravity, the radius cancels, and equilibrium is possible only at one special mass — the neutron-star analogue of the Chandrasekhar limit. But here honesty demands two caveats. First, these objects are so compact that general relativity is essential: the Tolman–Oppenheimer–Volkoff equation replaces Newtonian hydrostatics and lowers the free-gas maximum to $\approx0.7\,M_\odot$ . Second, real nuclear matter is strongly interacting, and repulsive nuclear forces stiffen the equation of state and push the observed maximum back up to $\sim2$ – $2.3\,M_\odot$ . The ideal-gas calculation does not give the exact ceiling — but it gives the correct scale (a $\sim10$ km, $\sim M_\odot$ ball of degenerate neutrons) and the essential physics: Pauli pressure versus gravity.

Key takeaways

A neutron star is supported against gravitational collapse by the degeneracy pressure of a nearly ideal, fully degenerate neutron Fermi gas at $T=0$ — exactly as electron degeneracy supports a white dwarf, with $m_e\to m_n$ .
Balancing the gravitational energy $-\alpha GM^2/R$ against the Fermi kinetic energy and minimizing gives the equilibrium radius $R\propto M^{-1/3}$ in the non-relativistic regime: more mass, smaller star.
When the neutrons become ultra-relativistic the kinetic energy scales as $1/R$ just like gravity, the radius drops out, and equilibrium survives only at a single maximum mass.
Because $m_n\gg m_e$ , the gas can be compressed to nuclear density ( $\sim10^{18}\,\mathrm{kg\,m^{-3}}$ ) and the star is $\sim500$ times smaller than a white dwarf — a radius of only $\sim10$ km.
The ideal Newtonian limiting mass is $\sim5.7\,M_\odot$ ; general relativity (TOV) lowers the free-gas maximum to $\approx0.7\,M_\odot$ , while nuclear repulsion raises the real maximum to the observed $\sim2$ – $2.3\,M_\odot$ .
Typical numbers: $M\sim1.4\,M_\odot$ , $R\sim10$ km, central density $\sim10^{18}\,\mathrm{kg\,m^{-3}}$ , neutron $E_F\sim$ tens of MeV; degenerate because $T_F\sim10^{12}\,$ K dwarfs $T\sim10^{8}\,$ K.
Key idea: a neutron star is a self-gravitating, fully degenerate ideal Fermi gas — equate gravitational and Fermi kinetic energies. Non-relativistically the balance gives a stable $R\propto M^{-1/3}$ star; ultra-relativistically the radius cancels and a maximum (Chandrasekhar-type) mass appears, though GR and nuclear forces set the true value.

The core insight

A neutron star is a self-gravitating, fully degenerate ideal Fermi gas: equate gravitational and Fermi kinetic energies. Non-relativistically the balance gives a stable $R\propto M^{-1/3}$ star; in the ultra-relativistic limit the radius cancels and a maximum (Chandrasekhar-type) mass appears — with the precise ceiling set by general relativity and nuclear interactions.

Where this lives today

Integrating the TOV equation with a realistic nuclear equation of state predicts a mass–radius curve — $R\approx10$ – $12$ km nearly independent of mass over a broad range, terminating at a maximum. The discovery of $\approx2\,M_\odot$ pulsars (e.g. PSR J0740+6620) rules out the softest equations of state and directly constrains the physics of dense matter, with the degeneracy argument here providing the conceptual backbone. Cooling of young neutron stars over $10^5$ – $10^6$ years is governed by the tiny thermal excitation near the Fermi surface (a Sommerfeld heat capacity $\propto T/T_F$ , the same result as for metals). And in a binary neutron-star merger, the tidal deformability — set by the stiffness of the degenerate-plus-nuclear equation of state — imprints on the gravitational-wave signal (GW170817), turning the abstract pressure–density relation of Fermi statistics into a directly measured astrophysical quantity.

Origins · 1 minFritz Zwicky and how this was discovered

Fritz Zwicky

1898–1974 · the astronomer who, with Walter Baade, coined 'neutron star'

Fritz Zwicky was a Swiss-American astrophysicist of famously combative brilliance. In 1934, barely two years after Chadwick discovered the neutron, Zwicky and Walter Baade published a pair of papers proposing that supernovae mark the transition of an ordinary star into a neutron star — an object made almost entirely of degenerate neutrons packed to nuclear density. They coined the term and predicted the object's compactness decades before pulsars were found (1967). Zwicky also inferred the existence of dark matter from galaxy-cluster velocities (1933) and catalogued thousands of supernovae. The neutron star, held up by exactly the Pauli degeneracy pressure treated in this article, was his most spectacular prediction.

History

1932

— James Chadwick discovers the neutron, making the idea of a neutron-composed star conceivable.

1934

— Walter Baade and Fritz Zwicky propose neutron stars as supernova remnants and coin the term.

1939

— Oppenheimer and Volkoff derive the TOV equation and a free-gas maximum mass of about 0.7 solar masses.

1967

— Jocelyn Bell Burnell and Antony Hewish discover pulsars, the first observed neutron stars.

2017

— GW170817, a binary neutron-star merger, is detected in gravitational waves and light, constraining the equation of state.

The physical picture

Think of the neutrons as a quantum gas in a box, but now the walls are gravity. The Pauli exclusion principle forbids two neutrons from sharing a state, so packing $N$ of them into radius $R$ forces them into ever-higher momentum states; the resulting outward push is the degeneracy pressure. The star sits precisely where that push balances its own weight. This is the same machinery that gives Degeneracy Pressure of Fermions its $P\propto n^{5/3}$ law and Fermi Energy, Momentum and Fermi Temperature their scales — recalled here by name and applied with $m_n$ in place of $m_e$ . The size comparison is the punchline. A white dwarf, supported by electron degeneracy, is roughly Earth-sized ( $\sim10^4$ km). A neutron star, supported by neutron degeneracy, is a thousand times smaller — about 10 km — because the heavier neutron ( $m_n\approx1839\,m_e$ ) can be squeezed to nuclear density before relativistic softening caps the mass. Both hold up a mass near the Sun's. The figure below places all three on a logarithmic size axis.

Figure

Size comparison on a log radius scale — the Sun (~10⁵ km), a white dwarf (~10⁴ km, Earth-sized), and a neutron star (~10 km), all near one solar mass

Why neutrons rather than electrons? In a white dwarf the electrons provide the pressure. As density rises, the electron Fermi energy exceeds the neutron–proton mass difference, electron capture $p+e^-\to n+\nu$ proceeds, and neutrons take over. Because the neutron is so much heavier, the gas can reach nuclear density before its momenta become relativistic — which is exactly why a neutron star is so tiny and so dense.

Figure

Energy versus radius — the non-relativistic neutron gas has a stable minimum (a real star); the ultra-relativistic gas shares gravity's 1/R and admits equilibrium only at one mass

The core equations

The neutron Fermi energy follows from the same phase-space count as any degenerate Fermi gas, with spin degeneracy $g=2$ . It is the kinetic energy of the highest filled state at $T=0$ .

E_F=\frac{\hbar^2}{2m_n}\left(3\pi^2 n\right)^{2/3},\qquad p_F=\hbar(3\pi^2 n)^{1/3},\qquad n=\frac{N}{V} \qquad (1)

where: $E_F$ — Fermi energy; $m_n$ — neutron mass, $1.675\times10^{-27}\,\mathrm{kg}$ ; $n$ — neutron number density; $p_F$ — Fermi momentum. Read it as: the density fixes how big the Fermi sphere is, and $E_F$ is the energy at its surface. At nuclear density $E_F$ is tens of MeV. Whether relativity matters is decided by a single dimensionless ratio comparing the Fermi momentum to the neutron rest-mass scale $m_n c$ .

\varepsilon(p)=\sqrt{p^2c^2+m_n^2c^4},\qquad x_F\equiv\frac{p_F}{m_n c} \qquad (2)

where: $\varepsilon(p)$ — single-particle energy including rest mass; $x_F$ — relativity parameter, $x_F\ll1$ non-relativistic, $x_F\gg1$ ultra-relativistic. The full relativistic degenerate-gas equation of state (its Chandrasekhar functions) is owned by Relativistic Degenerate Fermi Gas; here we use its two limits. The degenerate pressure comes from the momentum flux of the fast neutrons, exactly as in the electron case, with $m_n$ in the dispersion.

P=\frac{1}{3\pi^2\hbar^3}\int_0^{p_F}\frac{p^4 c^2}{\sqrt{p^2c^2+m_n^2c^4}}\,\mathrm{d}p\ \xrightarrow{x_F\ll1}\ \frac{(3\pi^2)^{2/3}}{5}\frac{\hbar^2}{m_n}n^{5/3}\ \xrightarrow{x_F\gg1}\ \frac{(3\pi^2)^{1/3}}{4}\hbar c\,n^{4/3} \qquad (3)

where: the non-relativistic limit gives the stiff $\Gamma=\tfrac53$ polytrope (stable), the ultra-relativistic limit the soft $\Gamma=\tfrac43$ polytrope (marginal). It is the crossover between these two that both permits a stable star at low mass and forbids one above a critical mass. Newtonian hydrostatic equilibrium is the leading description, but for a neutron star the relativistic corrections are decisive, so the true structure equation is the Tolman–Oppenheimer–Volkoff (TOV) equation.

\frac{\mathrm dP}{\mathrm dr}=-\frac{G m\rho}{r^2}\left(1+\frac{P}{\rho c^2}\right)\!\left(1+\frac{4\pi r^3 P}{m c^2}\right)\!\left(1-\frac{2Gm}{rc^2}\right)^{-1} \qquad (4)

where: $P/\rho c^2$ — pressure-to-energy-density ratio; $2Gm/rc^2$ — compactness, approaching 1 near the Schwarzschild radius. Each of the three factors is $>1$ and strengthens effective gravity, so the TOV maximum mass is smaller than the Newtonian one — for a free neutron gas, $\approx0.7\,M_\odot$ . When all three factors $\to1$ (weak field, $P\ll\rho c^2$ ) it reduces to Newtonian hydrostatics.

Limiting cases

The behaviour is governed by $x_F=p_F/(m_n c)$ and by $T/T_F$ :

Case	Result	Meaning
Non-relativistic ( $x_F\ll1$ )	$P\propto n^{5/3}$ , $R\propto M^{-1/3}$	Light neutron stars (and the bulk of a white dwarf). Stable: pressure ( $1/R^2$ ) beats gravity ( $1/R$ ). More mass, smaller star.
Ultra-relativistic ( $x_F\gg1$ )	$P\propto n^{4/3}$ , $M\to M_{\max}=\mathrm{const}$	Approached as $M\to M_{\max}$ . Pressure and gravity scale identically with $R$ ; equilibrium fixes a unique limiting mass.
Degenerate limit ( $T\ll T_F$ )	$T\sim10^{8}\,\mathrm{K}\ll T_F\sim10^{12}\,\mathrm{K}$	Thermal corrections $\mathcal O((T/T_F)^2)\sim10^{-8}$ ; the $T=0$ equation of state is essentially exact for structure.
GR + nuclear-force corrections	$M_{\max}^{\text{free,GR}}\approx0.7\,M_\odot\ \to\ M_{\max}^{\text{real}}\sim2\,M_\odot$	TOV (GR) lowers the free-gas limit; nuclear repulsion stiffens the equation of state and raises the observed maximum. The ideal estimate is an order-of-magnitude guide.

Derivation AEquilibrium radius from energy minimization, non-relativistic

The goal is to derive the neutron star's radius by balancing the Fermi kinetic energy of the degenerate gas against its own gravitational self-energy and minimizing over $R$ . We model the star as $N$ neutrons of total mass $M=N m_n$ in a uniform sphere of radius $R$ , write the total energy as a sum of a $1/R^2$ kinetic term and a $-1/R$ gravitational term, and locate the minimum. The competition of these two powers of $R$ is the entire content.

Step 1. Kinetic energy of the Fermi sea. The average energy of a $T=0$ Fermi gas is $\tfrac35 E_F$ , so $E_{\mathrm{kin}}=\tfrac35 N E_F=\tfrac35 N\tfrac{\hbar^2}{2m_n}(3\pi^2 n)^{2/3}$ . With $n=3N/(4\pi R^3)$ this is $\propto N^{5/3}/R^2$ . Why: the $\tfrac35$ is the purely geometric ratio of $\int p^4\mathrm dp$ to $E_F\int p^2\mathrm dp$ over the sphere, and $n^{2/3}\propto R^{-2}$ carries the radius dependence.

Step 2. Gravitational energy. The self-energy of a uniform sphere is $E_{\mathrm{grav}}=-\tfrac35 GM^2/R$ . Why: the $\tfrac35$ here is the uniform-density coefficient of the gravitational binding integral; it scales as $1/R$ , one power weaker than the kinetic term.

Step 3. Total energy. Writing $E(R)=A/R^2-B/R$ with $A\propto N^{5/3}$ and $B=\tfrac35 GM^2\propto M^2$ , the two terms compete: kinetic dominates at small $R$ (it blows up faster), gravity at large $R$ . Why: because the powers of $R$ differ ( $-2$ versus $-1$ ), the sum must turn over — there is a genuine minimum, unlike the ultra-relativistic case where both are $1/R$ .

Step 4. Minimize. Setting $\mathrm dE/\mathrm dR=-2A/R^3+B/R^2=0$ gives $R=2A/B$ , a true minimum since $E\to+\infty$ as $R\to0$ . Substituting $A\propto M^{5/3}$ and $B\propto M^2$ yields $R\propto M^{5/3}/M^2=M^{-1/3}$ . Why: the negative exponent is the signature of degeneracy support — adding mass strengthens gravity ( $\propto M^2$ ) faster than the Pauli term, so the star shrinks. *Verified: $E''>0$ at $R=2A/B$ (genuine minimum) and the $R\propto M^{-1/3}$ scaling are confirmed symbolically, and $R\approx11$ km for $M=1.4\,M_\odot$ numerically, in verify_neutron-stars-degeneracy.py.*

Go deeper · Derivation B · 1 minExistence of a limiting mass, ultra-relativistic

Derivation BExistence of a limiting mass, ultra-relativistic

The goal is to show that when the neutrons become ultra-relativistic, a maximum mass appears — the neutron-star analogue of the Chandrasekhar limit. We repeat the energy argument with the relativistic dispersion $\varepsilon\approx p c$ , note that now both the kinetic and gravitational energies scale as $1/R$ , and see that the radius drops out, leaving a condition on mass alone. The strategy is identical to the white-dwarf case with $m_n$ replacing $\mu_e m_p$ .

Step 1. Ultra-relativistic kinetic energy. With $\varepsilon\approx p_F c$ the average energy is $\tfrac34 E_F$ and $E_{\mathrm{kin}}\approx\tfrac34 N\hbar c(3\pi^2 n)^{1/3}\propto \hbar c\,N^{4/3}/R$ . Why: linear dispersion makes the energy scale as $p_F\propto N^{1/3}/R$ , so the total carries a single power of $1/R$ — the same as gravity.

Step 2. Total energy. Now $E(R)=\tfrac1R(a\,\hbar c\,N^{4/3}-b\,G m_n^2 N^2)$ with $M=N m_n$ . Why: both terms share the $1/R$ factor, so the sign of the bracket — not $R$ — decides whether the star expands or collapses.

Step 3. Stability criterion. If $a\,\hbar c\,N^{4/3}>b\,G m_n^2 N^2$ the bracket is positive, $E$ falls as $R$ grows, and the gas expands until it becomes non-relativistic and stabilizes; if the inequality reverses, $E\to-\infty$ as $R\to0$ — unbounded collapse. The threshold is $N<N_{\max}\sim(\hbar c/G m_n^2)^{3/2}$ . Why: setting the bracket to zero gives a critical particle number built only from fundamental constants and the neutron mass.

Step 4. Limiting mass. Thus $M_{\max}=N_{\max}m_n\sim M_{\mathrm{Pl}}^3/m_n^2$ with $M_{\mathrm{Pl}}=\sqrt{\hbar c/G}$ . The careful Lane–Emden ( $n=3$ ) coefficient gives $\approx5.7\,M_\odot$ for an ideal Newtonian gas — but GR (TOV) lowers the free-gas value to $\approx0.7\,M_\odot$ . Why: remarkably the maximum mass is set by the Planck mass cubed over the neutron mass squared; the ideal-gas number is only a guide because GR and nuclear forces both matter at these compactnesses. *Verified: $N_{\max}\propto(\hbar c/G m_n^2)^{3/2}$ symbolically, $M_{\mathrm{Pl}}=2.18\times10^{-8}$ kg, $M_{\mathrm{Pl}}^3/m_n^2=1.85\,M_\odot$ , and the ideal Newtonian $\approx5.7\,M_\odot$ are recomputed in verify_neutron-stars-degeneracy.py; the TOV free-gas $0.7\,M_\odot$ is the literature value (Oppenheimer & Volkoff 1939).*

Common trapNewtonian gravity gives the exact maximum mass of a neutron star

It does not — and this is the central honesty of the topic. The ideal-gas, Newtonian calculation gives $\sim5.7\,M_\odot$ , which is wrong as a physical maximum. General relativity is essential at these compactnesses ( $2GM/Rc^2\sim0.3$ – $0.5$ ): the TOV equation strengthens effective gravity and lowers the free-gas maximum to $\approx0.7\,M_\odot$ . Then repulsive nuclear interactions stiffen the real equation of state and raise the observed maximum to $\sim2$ – $2.3\,M_\odot$ . So the Newtonian degeneracy argument is trustworthy only for the scale — a $\sim10$ km, $\sim M_\odot$ object — and for the existence of a limit, not its precise value. Also: never use the non-relativistic equation of state near the maximum mass, since the limit is entirely a relativistic effect.

Worked examples

Example 1 — Why the neutron gas is degenerate: the Fermi temperature.

Step 1. Take a representative interior density $\rho\sim4\times10^{17}\,\mathrm{kg\,m^{-3}}$ (a few times nuclear saturation). The neutron number density is $n=\rho/m_n=(4\times10^{17})/(1.675\times10^{-27})=\mathbf{2.4\times10^{44}\ m^{-3}}$ — a few neutrons per cubic femtometre.

Step 2. The Fermi momentum is $p_F=\hbar(3\pi^2 n)^{1/3}$ . With $3\pi^2 n=7.1\times10^{45}$ , this gives $p_F c\approx60$ – $80$ MeV-scale energies; computing the relativistic Fermi energy $E_F=\sqrt{(p_F c)^2+(m_n c^2)^2}-m_n c^2\approx\mathbf{74\ MeV}$ (neutrons are mildly relativistic at this density).

Step 3. The Fermi temperature is $T_F=E_F/k_B=(74\times1.602\times10^{-13})/(1.381\times10^{-23})=\mathbf{8.5\times10^{11}\ K}\sim10^{12}\,$ K. Compared to the interior $T\sim10^{8}\,$ K, $T/T_F\sim10^{-4}$ , so thermal corrections are $\mathcal O((T/T_F)^2)\sim\mathbf{10^{-8}}$ — utterly negligible. The $T=0$ Fermi-gas equation of state is justified to eight decimal places.

Example 2 — The radius of a solar-mass neutron star.

Step 1. Using the non-relativistic estimate $R=(9\pi/4)^{2/3}\hbar^2/(G m_n^{8/3}M^{1/3})$ for $M=1.4\,M_\odot=2.78\times10^{30}\,$ kg: the prefactor $(9\pi/4)^{2/3}=3.66$ , and $\hbar^2/(G m_n^{8/3})=(1.055\times10^{-34})^2/[(6.674\times10^{-11})(1.675\times10^{-27})^{8/3}]$ .

Step 2. Evaluating, $m_n^{8/3}=(1.675\times10^{-27})^{8/3}=3.96\times10^{-72}$ , so the ratio is $\approx1.11\times10^{-68}/(2.64\times10^{-82})\approx4.2\times10^{13}$ , and dividing by $M^{1/3}=(2.78\times10^{30})^{1/3}=1.4\times10^{10}$ gives $R\approx3.66\times4.2\times10^{13}/1.4\times10^{10}\approx\mathbf{1.1\times10^{4}\ m=11\ km}$ .

Step 3. So a solar-mass neutron star is about 10 km across — a thousand times smaller than a white dwarf of the same mass, because $R\propto1/m^{...}$ shrinks with the much larger neutron mass. Real neutron-star radii from X-ray and gravitational-wave observations are $R\approx10$ – $12$ km, in gratifying agreement with this order-of-magnitude degeneracy estimate.

Going further

References. Pathria & Beale, Statistical Mechanics, 3rd ed., Ch. 8 (Ideal Fermi Systems; white-dwarf/neutron-star application) · Kardar, Statistical Physics of Particles, Ch. 7 · Reif, Fundamentals of Statistical and Thermal Physics, Ch. 9 · Shapiro & Teukolsky, Black Holes, White Dwarfs and Neutron Stars (1983) · Oppenheimer & Volkoff, 'On Massive Neutron Cores', Phys. Rev. 55, 374 (1939).

Provenance. Equations are standard Fermi-gas and stellar-structure results. This article recalls the pressure law and Fermi scales from Degeneracy Pressure of Fermions, Fermi Energy, Momentum and Fermi Temperature, and the white-dwarf mechanism from White Dwarf Stars and the Chandrasekhar Limit and Relativistic Degenerate Fermi Gas by name, and owns the neutron-degeneracy scaling (the $m_e\to m_n$ radius shrink to $\sim10$ km) and the honest GR/nuclear-force caveat. Scope note: the Newtonian ideal-gas maximum ( $\sim5.7\,M_\odot$ ) is not the physical ceiling — GR (TOV) gives $\approx0.7\,M_\odot$ for a free gas and realistic nuclear forces $\sim2$ – $2.3\,M_\odot$ ; the Newtonian argument is reliable only for the scale and the existence of a limit, machine-verified in verify_neutron-stars-degeneracy.py.

Builds on: Degenerate Fermi Gas at Zero Temperature · Fermi Energy, Momentum and Fermi Temperature · Degeneracy Pressure of Fermions · Density of States for Free Particles Leads to: Sommerfeld Expansion · Electronic Heat Capacity of Metals

Why can't electron degeneracy support a neutron star like it supports a white dwarf?

As density rises, the electron Fermi energy climbs above the neutron–proton mass difference ( $\approx1.3$ MeV), and it becomes energetically favourable for electrons to be captured on protons, $p+e^-\to n+\nu_e$ . This removes the very electrons that were providing the pressure. The matter neutronizes, and support must come from neutron degeneracy instead. Because the neutron is $\sim1839$ times heavier, its degeneracy pressure allows compression to nuclear density — which is why the resulting object is so much smaller and denser than a white dwarf.

If it's called a 'neutron' star, are there really no protons or electrons at all?

Not exactly — it's dominated by neutrons but not pure. A small equilibrium fraction of protons and electrons remains, set by requiring that neutron decay and electron capture balance (beta equilibrium): typically a few percent protons, with an equal number of electrons for charge neutrality. This 'npe' matter, plus possibly muons and exotic phases deeper in, is what a realistic equation of state must describe. The ideal-neutron-gas model in this article keeps only the dominant neutron Fermi sea, which is why it captures the scale but not the fine detail.

How can the same R ∝ M^(−1/3) law describe both white dwarfs and neutron stars?

Because both are self-gravitating non-relativistic degenerate Fermi gases; only the particle mass and what carries the mass density differ. Minimizing $E=A/R^2-B/R$ always gives $R\propto M^{-1/3}$ , with the overall coefficient depending on the fermion mass $m$ (through the kinetic term). For a white dwarf the electrons give the pressure while nucleons give the mass; for a neutron star the neutrons do both. The much larger neutron mass shrinks the coefficient, taking the radius from $\sim10^4$ km down to $\sim10$ km.

Does the TOV equation always lower the maximum mass compared to Newtonian gravity?

For a given equation of state, yes — the three relativistic factors in the TOV equation are each greater than one and all strengthen effective gravity, so the star can support less mass than Newtonian hydrostatics would predict. That is why the free-gas maximum drops from the Newtonian $\sim5.7\,M_\odot$ to the TOV $\approx0.7\,M_\odot$ . The observed $\sim2\,M_\odot$ maximum is higher only because real nuclear matter is much stiffer than a free gas; with the same stiff equation of state, GR still lowers the ceiling relative to the Newtonian result.

Why does the maximum mass depend on the Planck mass?

Because it comes from balancing a quantum pressure (built from $\hbar$ and $c$ ) against gravity (carrying $G$ ). The critical particle number is $N_{\max}\sim(\hbar c/G m_n^2)^{3/2}$ , so $M_{\max}=N_{\max}m_n\sim(\hbar c/G)^{3/2}/m_n^2=M_{\mathrm{Pl}}^3/m_n^2$ with $M_{\mathrm{Pl}}=\sqrt{\hbar c/G}$ . It is a beautiful fact that a macroscopic stellar mass is fixed by the Planck mass (a quantum-gravity scale) and the neutron mass alone — the same structure as the Chandrasekhar mass with $m_n$ in place of $\mu_e m_p$ .

How do observations actually pin down the neutron-star equation of state?

Several independent handles now exist. Precise pulsar masses (radio timing of binaries) give the maximum mass a viable equation of state must reach — the $\sim2\,M_\odot$ pulsars are decisive. NASA's NICER X-ray telescope measures radii to $\sim10$ % by modelling hot spots on rotating neutron stars. And gravitational waves from mergers like GW170817 measure the tidal deformability, which depends on how stiff the star is. Together they are narrowing the pressure–density relation of matter beyond nuclear density — the ultimate empirical test of the degeneracy-plus-interactions picture sketched here.