Equation of State of Dense Matter

Squeeze a gram of matter to the density inside a neutron star and every familiar structure dissolves. Atoms are crushed until their electrons form a single degenerate sea; nuclei swell, grow neutron-rich, and finally dissolve into a uniform fluid of neutrons laced with a few percent protons; and above nuclear saturation the strong force — repulsive at short range — takes over as the dominant source of pressure. The equation of state (EOS) is the one function that summarizes all of this: the cold, zero-temperature relation $P=P(\varepsilon)$ between pressure and energy density obeyed at each depth. Because thermal energies $k_BT$ inside even a young neutron star are minuscule beside Fermi energies of hundreds of MeV, temperature drops out and the EOS collapses to a single curve — a one-parameter map from density to pressure. That curve is the crucial microphysical input to stellar structure. Feed it into the Tolman–Oppenheimer–Volkoff equation of general-relativistic hydrostatic equilibrium and it fixes everything macroscopic: the radius at each mass, the maximum mass before collapse to a black hole, the moment of inertia, the tidal deformability. The map from EOS to the mass–radius relation $M(R)$ is essentially one-to-one, turning the whole enterprise into an inverse problem — measure masses and radii of real stars, and read the strong interaction back out at densities no laboratory can reach. Two-solar-mass pulsars demand a stiff EOS from below; NICER radii and the GW170817 tidal bound squeeze it from the side and above. Between these anchors lies the least-understood regime of nuclear physics.

Key takeaways

The equation of state of dense matter is the cold relation $P(\varepsilon)$ spanning densities from $\sim10^{6}\ \mathrm{g\,cm^{-3}}$ in white-dwarf material to several times nuclear saturation $\rho_0\simeq2.7\times10^{14}\ \mathrm{g\,cm^{-3}}$ ( $n_0\simeq0.16\ \mathrm{fm^{-3}}$ ) in a neutron-star core.
Fed into the Tolman–Oppenheimer–Volkoff equation it maps one-to-one onto the mass–radius relation $M(R)$ and fixes the maximum mass; a stiffer EOS gives larger radii and higher $M_\mathrm{max}$ .
In degenerate regimes the EOS is polytropic — non-relativistic gives $P\propto\rho^{5/3}$ ( $\Gamma=5/3$ ), ultra-relativistic gives $P\propto\rho^{4/3}$ ( $\Gamma=4/3$ ), and $\Gamma=4/3$ is exactly the threshold of gravitational instability.
Above $\rho_0$ the pressure is set by the strong force — symmetry energy, three-body forces, possibly hyperons or deconfined quarks — and cannot be computed from first principles with certainty.
Observations bracket the curve: $2\,M_\odot$ pulsars force $M_\mathrm{max}\gtrsim2\,M_\odot$ , NICER gives $R\approx11$ – $13\ \mathrm{km}$ , and GW170817's bound $\Lambda_{1.4}\lesssim800$ disfavours the stiffest models.
Key idea: a single microphysical function $P(\varepsilon)$ — the cold, dense equation of state — is the sole input beyond gravity that fixes the entire mass–radius relation and maximum mass of neutron stars through the TOV equation, so measuring $M$ and $R$ inverts to constrain the strong interaction at supranuclear density.

The core insight

A single microphysical function $P(\varepsilon)$ — the cold, dense equation of state — is the sole input beyond gravity that determines the entire mass–radius relation and maximum mass of neutron stars via the Tolman–Oppenheimer–Volkoff equation; measuring $M$ and $R$ inverts to constrain the strong interaction at supranuclear density.

Where this lives today

The dense-matter EOS is a frontier because three independent channels now converge on it. Radio pulsar timing weighs the heaviest stars: PSR J0348+0432 ( $2.01\,M_\odot$ ) and PSR J0740+6620 ( $2.08\,M_\odot$ ) set a hard floor, ruling out every EOS too soft to reach $2\,M_\odot$ and creating the "hyperon puzzle" — hyperons in the core soften the EOS and push $M_\mathrm{max}$ below the observed value, so their interactions must be strongly repulsive or their onset suppressed. NASA's NICER X-ray telescope reads mass and radius simultaneously by modelling the pulse profiles of hot magnetic-polar-cap spots, delivering $R\approx11$ – $13\ \mathrm{km}$ for both a $1.4\,M_\odot$ and a $2.08\,M_\odot$ star and pinning $P(\varepsilon)$ at a few times $\rho_0$ . Gravitational waves add the tidal deformability: GW170817's bound $\Lambda_{1.4}\lesssim800$ excludes the stiffest, largest-radius models. In parallel, terrestrial nuclear experiments — neutron-skin measurements like PREX-II and heavy-ion collisions — constrain the symmetry-energy slope $L$ that controls the radius, tying a laboratory nucleus to a star. Ahead lie next-generation detectors (Einstein Telescope, Cosmic Explorer) sensitive to post-merger oscillations, and the open question of whether a deconfined quark-matter core exists inside the most massive neutron stars.

The physical picture

Picture matter held up not by heat but by the Pauli principle: confine fermions to a smaller volume and they are forced into ever higher momentum states, so a quantum kinetic pressure resists compression even at $T=0$ . The single most important quantity is the stiffness — how steeply $P$ climbs with $\varepsilon$ . A stiff EOS resists compression and builds large, massive stars; a soft EOS yields small stars with a low mass ceiling. Compact-star structure is a contest between this degeneracy pressure (augmented above nuclear density by the repulsive strong force) and self-gravity. As you descend into a compact star the composition changes in distinct layers, each governed by different physics. The first figure is that ladder: from electron-degenerate white-dwarf matter, through the crust where nuclei neutronize and free neutrons "drip" out at $\rho_\mathrm{drip}\approx4\times10^{11}\ \mathrm{g\,cm^{-3}}$ , down to the uniform strong-interaction fluid of the core above $\rho_0\simeq2.7\times10^{14}\ \mathrm{g\,cm^{-3}}$ .

Figure

The density ladder of compact-star matter — pressure passes from degenerate electrons, through neutron drip at ρ ≈ 4×10¹¹ g cm⁻³, to the strong-force core above nuclear saturation ρ₀ ≈ 2.7×10¹⁴ g cm⁻³ (n₀ ≈ 0.16 fm⁻³)

The physics that decides a star's fate lives at the crossover from non-relativistic to ultra-relativistic degeneracy. For $p_F\ll mc$ the pressure scales as $P\propto\rho^{5/3}$ and gravity can always be balanced; for $p_F\gg mc$ it softens to $P\propto\rho^{4/3}$ , where pressure and gravity scale identically with radius, equilibrium turns marginal, and a maximum mass appears. Above $\rho_0$ the strong force intervenes: nucleons repel hard at short range, adding pressure far beyond the free-gas value — which is why real neutron stars reach $2\,M_\odot$ where a free gas caps near $0.7\,M_\odot$ . The payoff is the map from EOS to observables, the mass–radius relation of the second figure.

Figure

The mass–radius relation is a one-to-one image of the EOS: a stiff EOS (blue) yields large radii and reaches beyond 2 M⊙, while a soft EOS (red) is ruled out because it cannot support the observed 2 M⊙ pulsars; NICER radii cluster near 11–13 km

The core equations

The starting point is the pressure of a completely degenerate fermion gas — particles filling every momentum state up to the Fermi surface. Pressure is one-third the flux of momentum $p$ carried at velocity $v$ , integrated over the filled Fermi sphere, with two spin states per cell $h^3$ .

P = \frac{1}{3}\,n\,\langle p\,v\rangle = \frac{8\pi}{3h^3}\int_0^{p_F} p^3\,v(p)\,\mathrm{d}p \qquad (1)

where: $P$ — pressure; $n$ — number density of the fermions; $p$ — particle momentum with velocity $v=\partial E/\partial p$ ; $p_F$ — Fermi momentum, fixed by $n=(8\pi/3h^3)\,p_F^3$ for spin- $\tfrac12$ ; $h$ — Planck's constant. Reading it: at $T=0$ the only free parameter is $p_F$ (equivalently $n$ ), so pressure is a pure function of density. When Fermi momenta are small compared with $mc$ the dispersion is non-relativistic, $E\approx p^2/2m$ , and the integral gives a $\Gamma=5/3$ polytrope. For electron-degenerate matter the mass density is carried by the nucleons, hence $\mu_e m_H$ .

P = K_\mathrm{NR}\,\rho^{5/3},\qquad K_\mathrm{NR}=\frac{h^2}{20\,m_e}\left(\frac{3}{\pi}\right)^{2/3}\frac{1}{(\mu_e m_H)^{5/3}} \qquad (2)

where: $K_\mathrm{NR}$ — non-relativistic polytropic constant; $m_e$ — electron mass (for electron-degenerate matter); $\mu_e$ — mean molecular weight per electron; $m_H$ — hydrogen mass; $\rho$ — mass density. Reading it: the pressure carries $1/m_e$ , which is why the light electrons, not the heavy ions, supply it. Push the same gas to $p_F\gg mc$ and the particles move at essentially $c$ ; the dispersion becomes $E\approx pc$ and the polytrope softens to $\Gamma=4/3$ — the marginal case, where the constant no longer contains the particle mass.

P = K_\mathrm{UR}\,\rho^{4/3},\qquad K_\mathrm{UR}=\frac{hc}{8}\left(\frac{3}{\pi}\right)^{1/3}\frac{1}{(\mu_e m_H)^{4/3}} \qquad (3)

where: $K_\mathrm{UR}$ — ultra-relativistic polytropic constant; $c$ — speed of light. Reading it: because $K_\mathrm{UR}$ is fixed by fundamental constants alone, the mass it supports is a fixed number — the origin of the Chandrasekhar-type mass ceiling. Real cores are not free gases; the strong force adds an interaction energy. At $T=0$ the entire EOS follows from one function — the energy per baryon — through a thermodynamic identity requiring only $\varepsilon(n)$ .

P = n^2\,\frac{\mathrm{d}(\varepsilon/n)}{\mathrm{d}n} = n\,\frac{\mathrm{d}\varepsilon}{\mathrm{d}n} - \varepsilon \qquad (4)

where: $\varepsilon$ — total energy density (rest mass + kinetic + interaction); $n$ — baryon number density; $\varepsilon/n$ — energy per baryon. Reading it: nuclear theory need only predict how the energy per baryon rises with compression — differentiate and the pressure, hence the star, follows. The dominant unknown in $\varepsilon(n)$ is the cost of making matter neutron-rich, the nuclear symmetry energy. Expanding the energy per nucleon about symmetric matter in the neutron excess $(1-2x)$ isolates it.

\frac{E}{A}(n,x)\simeq E_0(n) + S(n)\,(1-2x)^2,\qquad x=\frac{n_p}{n} \qquad (5)

where: $E_0(n)$ — energy per nucleon of symmetric matter; $S(n)$ — symmetry energy, $\approx30$ – $34\ \mathrm{MeV}$ at saturation; $x$ — proton fraction ( $x\to0$ for pure neutron matter); $n$ — baryon density. Reading it: the density slope $L=3n_0\,\mathrm{d}S/\mathrm{d}n|_{n_0}$ sets the pressure of neutron matter near saturation, hence the radius of a $1.4\,M_\odot$ star. No EOS is free to be arbitrarily stiff: the adiabatic sound speed cannot exceed light, which caps $\mathrm{d}P/\mathrm{d}\varepsilon$ and hence the maximum mass.

c_s^2 = \frac{\mathrm{d}P}{\mathrm{d}\varepsilon}\,c^2 \le c^2 \qquad (6)

where: $c_s$ — adiabatic sound speed; $\mathrm{d}P/\mathrm{d}\varepsilon$ — local EOS stiffness (dimensionless in $c=1$ units). Reading it: the stiffest causal EOS gives the absolute mass bound $M_\mathrm{max}\lesssim3\,M_\odot$ . Finally, the observable a gravitational-wave detector reads directly is the tidal deformability, which depends on the EOS through the radius and the Love number and scales as the fifth power of the inverse compactness.

\Lambda = \frac{2}{3}\,k_2\left(\frac{c^2 R}{G M}\right)^5 \qquad (7)

where: $\Lambda$ — dimensionless tidal deformability (bounded by GW170817); $k_2$ — second tidal Love number, itself EOS-dependent; $R,\,M$ — stellar radius and gravitational mass. Reading it: the fifth-power dependence on $R$ makes $\Lambda$ an exquisitely sharp probe of the radius, hence of the EOS at a few times $\rho_0$ .

Derivation APressure of a completely degenerate relativistic fermion gas

The goal is to obtain the exact closed-form pressure of a cold fermion gas that interpolates between the two polytropic limits, so a single function $f(x)$ contains both the $\Gamma=5/3$ white-dwarf regime and the $\Gamma=4/3$ marginally stable regime. We integrate the momentum-flux definition (1) over the filled Fermi sphere with the relativistic dispersion, then extract the low- and high-density limits and show they reproduce equations (2) and (3). This is Chandrasekhar's classical result and the microphysical spine of both white-dwarf and free-neutron-star structure.

Step 1. Define the dimensionless Fermi momentum $x\equiv p_F/(mc)$ . Counting two spin states per cell $h^3$ inside a sphere of radius $p_F$ gives the number density $n=\dfrac{8\pi}{3h^3}p_F^3$ , so $n\propto x^3$ and the whole EOS is a function of $x$ alone. Why: at $T=0$ the Fermi surface is sharp, so one parameter — the radius of the filled sphere — carries all the thermodynamics.

Step 2. Insert the relativistic velocity $v=pc^2/E$ with $E=\sqrt{(pc)^2+(mc^2)^2}$ into (1): $P=\dfrac{8\pi}{3h^3}\displaystyle\int_0^{p_F}\dfrac{p^4c^2}{\sqrt{(pc)^2+(mc^2)^2}}\,\mathrm{d}p$ . Why: the only physics beyond counting states is the correct dispersion relation, which smoothly connects the $E\approx p^2/2m$ and $E\approx pc$ regimes.

Step 3. Substitute $u=p/(mc)$ so the integral runs from $0$ to $x$ ; it evaluates in closed form to $P=\dfrac{\pi m^4c^5}{3h^3}\,f(x)$ with $f(x)=x(2x^2-3)\sqrt{x^2+1}+3\,\mathrm{arcsinh}\,x$ . Why: the antiderivative is elementary once the dispersion is rationalized, giving a single analytic curve valid at all densities.

Step 4. Take $x\ll1$ : expanding gives $f(x)\to\tfrac{8}{5}x^5$ , so $P\propto p_F^5\propto n^{5/3}\propto\rho^{5/3}$ — exactly the $\Gamma=5/3$ polytrope of equation (2). Why: at low density the particles are non-relativistic, $E\approx p^2/2m$ , and the leading term must reproduce the ordinary degenerate-gas result.

Step 5. Take $x\gg1$ : now $f(x)\to2x^4$ , so $P\propto p_F^4\propto n^{4/3}\propto\rho^{4/3}$ — the $\Gamma=4/3$ polytrope of equation (3), the marginally stable limit. Why: at high density every particle moves at nearly $c$ , pressure loses its dependence on $m$ , and stiffness drops to the critical value that admits a maximum mass.

Go deeper · Derivation B · 1 minPolytropic mass–radius scaling and the origin of the maximum mass

Derivation BPolytropic mass–radius scaling and the origin of the maximum mass

The goal is to show, from hydrostatic equilibrium alone, why the ultra-relativistic index $\Gamma=4/3$ produces a mass independent of central density — a single critical mass — whereas any stiffer EOS supports a family of stable stars. We combine Newtonian hydrostatic balance with a polytrope $P=K\rho^\Gamma$ by dimensional analysis, solve how radius and mass scale with central density, and read off the exponent that flips sign at $\Gamma=4/3$ . This is the root of the Chandrasekhar limit and, with GR corrections, of every neutron-star maximum mass.

Step 1. For $P=K\rho^\Gamma$ , hydrostatic balance $\mathrm{d}P/\mathrm{d}r=-Gm\rho/r^2$ estimated with $M\sim\rho R^3$ gives $K\rho^\Gamma/R\sim G(\rho R^3)\rho/R^2$ , i.e. $K\rho^{\Gamma-2}\sim G\rho R^2$ . Why: dimensional balance of the pressure gradient against gravity captures the scaling without solving the full Lane–Emden equation.

Step 2. Solve for the radius: $R^2\sim(K/G)\,\rho^{\Gamma-2}$ , so $R\propto\rho_c^{(\Gamma-2)/2}$ where $\rho_c$ is the central density. Why: this already shows a stiffer gas ( $\Gamma>2$ ) swells with compression while a softer one shrinks — the geometric response to squeezing.

Step 3. Combine with $M\sim\rho_c R^3$ to get $M\propto\rho_c\,\rho_c^{3(\Gamma-2)/2}=\rho_c^{(3\Gamma-4)/2}$ . Why: the mass scaling is the single diagnostic that decides stability, and its exponent $(3\Gamma-4)/2$ is controlled entirely by $\Gamma$ .

Step 4. Read the exponent. For $\Gamma>4/3$ it is positive, so raising $\rho_c$ adds supportable mass and stable equilibria exist at every mass. For $\Gamma=4/3$ it vanishes: $M$ is fixed independent of $\rho_c$ , a unique critical mass $M_\mathrm{crit}$ . Why: when mass no longer depends on central density, there is exactly one mass the configuration can hold — the definition of a mass limit.

Step 5. For $\Gamma<4/3$ the exponent is negative: compressing the star lowers the mass it can support, so a perturbation runs away and the star collapses. Why: a decreasing $M(\rho_c)$ is the hallmark of gravitational instability, and it is why an ultra-relativistic gas — which sits exactly at $\Gamma=4/3$ — is only marginally stable.

Go deeper · Derivation C · 1 minNeutron drip and the cold pressure from the energy per baryon

Derivation CNeutron drip and the cold pressure from the energy per baryon

The goal is twofold: to justify equation (4) — that the entire cold EOS follows from the energy per baryon $\varepsilon/n$ — and to use it to locate the neutron-drip density where nuclei first shed free neutrons and the inner crust begins. We start from the first law at $T=0$ , convert a volume derivative into a density derivative, then apply beta equilibrium to find where the neutron chemical potential reaches the continuum, fixing the outer/inner-crust boundary at $\rho_\mathrm{drip}\approx4\times10^{11}\ \mathrm{g\,cm^{-3}}$ .

Step 1. At $T=0$ the first law reads $P=-\mathrm{d}E/\mathrm{d}V$ at fixed baryon number $N$ . Writing $E=N(\varepsilon/n)$ and $V=N/n$ gives $\mathrm{d}V=-N\,\mathrm{d}n/n^2$ , hence $P=n^2\,\mathrm{d}(\varepsilon/n)/\mathrm{d}n$ . Why: this reduces the whole thermodynamics to one modelled function, the energy per baryon, and rearranges to the compact form $P=n\,\mathrm{d}\varepsilon/\mathrm{d}n-\varepsilon$ .

Step 2. As density rises, the electron Fermi energy grows until electron capture $p+e^-\to n+\nu_e$ becomes favourable, neutronizing nuclei and driving the proton fraction down. Why: matter minimizes its energy under beta equilibrium, and at high electron pressure the cheapest way to add a baryon is as a neutron, not a proton plus electron.

Step 3. The neutron chemical potential $\mu_n$ climbs as nuclei grow neutron-rich. Free neutrons appear once $\mu_n$ reaches the free-neutron rest mass — equivalently, when the neutron separation energy of the equilibrium nucleus falls to zero. Why: below that point a neutron is bound to the lattice nucleus; at it, adding a neutron costs nothing, so neutrons spill into a continuum coexisting with the nuclei.

Step 4. Detailed nuclear mass tables place this at $\rho_\mathrm{drip}\approx4\times10^{11}\ \mathrm{g\,cm^{-3}}$ , the boundary between the outer crust (nuclei + electrons) and the inner crust (nuclei + electrons + a free-neutron superfluid). Why: the drip density is where a new pressure-bearing component — dripped neutrons — switches on, softening then, above $\rho_0$ , stiffening the EOS as the strong force takes hold.

Limiting cases

In clean limits the EOS collapses to a polytrope $P=K\rho^\Gamma$ or to a causal bound, and these brackets frame the complicated intermediate behaviour that nuclear theory must supply:

Case	Result	Meaning
Non-relativistic degeneracy	$P\propto\rho^{5/3},\ \Gamma=5/3$	$p_F\ll mc$ . White-dwarf electrons and crust neutrons. Always stable, since $\Gamma>4/3$ .
Ultra-relativistic degeneracy	$P\propto\rho^{4/3},\ \Gamma=4/3$	$p_F\gg mc$ . Marginal stability — pressure and gravity scale alike, giving a Chandrasekhar-type mass limit. The softest self-supporting EOS.
Causal / maximally stiff	$P=\varepsilon-\varepsilon_0,\ c_s=c$	$\mathrm{d}P/\mathrm{d}\varepsilon=c^2$ everywhere. Stiffest EOS allowed by causality; sets the absolute bound $M_\mathrm{max}\lesssim3\,M_\odot$ .
Free relativistic neutron gas	$M_\mathrm{max}\approx0.7\,M_\odot,\ R\approx9.6\ \mathrm{km}$	Oppenheimer–Volkoff (1939) with no interactions — far below observation, proving nuclear forces are indispensable.

Worked examples

1 · Are the neutrons relativistic at nuclear saturation? (confidence builder)

Estimate the dimensionless Fermi momentum $x=p_F/(m_nc)$ for pure neutron matter at nuclear saturation density $n_0=0.16\ \mathrm{fm^{-3}}$ , and say which polytropic regime applies.

Step 1. For a single spin- $\tfrac12$ species the Fermi momentum is $p_F=\hbar(3\pi^2 n)^{1/3}$ . Evaluate the argument: $3\pi^2\times0.16=4.74\ \mathrm{fm^{-3}}$ , whose cube root is $1.68\ \mathrm{fm^{-1}}$ , so $p_F=1.68\,\hbar\ \mathrm{fm^{-1}}$ .

Step 2. Convert to energy with $\hbar c=197.3\ \mathrm{MeV\,fm}$ : $p_Fc=197.3\times1.68=\mathbf{331\ MeV}$ . The neutron rest energy is $m_nc^2=939.6\ \mathrm{MeV}$ .

Step 3. Hence $x=p_Fc/(m_nc^2)=331/939.6=\mathbf{0.35}$ . Since $x<1$ , neutrons at saturation are only mildly relativistic and the EOS sits nearer the $\Gamma=5/3$ regime than $\Gamma=4/3$ — but by a few times $\rho_0$ in the core, $x\gtrsim1$ and relativistic softening is well underway. At $\rho_0$ , $x\approx0.35$ : mildly relativistic, closer to the stiff $\Gamma=5/3$ limit.

2 · Tidal deformability of a canonical neutron star (fully worked)

A $1.4\,M_\odot$ neutron star has radius $R=12\ \mathrm{km}$ and tidal Love number $k_2=0.09$ . Compute its dimensionless tidal deformability $\Lambda$ from equation (7) and compare with the GW170817 bound $\Lambda_{1.4}\lesssim800$ .

Step 1. Form the compactness using the convenient constant $GM_\odot/c^2=1.477\ \mathrm{km}$ : $\dfrac{GM}{c^2R}=\dfrac{1.477\times1.4}{12}=\dfrac{2.068}{12}=\mathbf{0.172}$ . The inverse compactness is $c^2R/GM=1/0.172=\mathbf{5.80}$ .

Step 2. Raise it to the fifth power: $5.80^5=5.80^2\times5.80^2\times5.80=33.64\times33.64\times5.80\approx\mathbf{6560}$ . The steep fifth-power dependence is why $\Lambda$ is such a sharp radius probe.

Step 3. Assemble equation (7): $\Lambda=\tfrac{2}{3}\times0.09\times6560=0.06\times6560=\mathbf{394}$ .

Step 4. Compare: $\Lambda\approx394$ sits comfortably below the GW170817 bound of $\approx800$ , so a $12\ \mathrm{km}$ radius is allowed. Because $\Lambda\propto R^5$ , a stiffer star with $R=14\ \mathrm{km}$ would give $\approx394\times(14/12)^5\approx700$ , pressing against the bound. $\Lambda\approx390$ , consistent with GW170817; a radius much above $13$ km would violate it — how a single event excludes the stiffest EOS.

Physics of the portrait

The equation of state of dense matter is where nuclear physics and general relativity were first forced to meet — a curve assembled over ninety years from degeneracy theory, the neutron's discovery, and the first computations of collapsed-star structure.

Origins · 1 minequation of state of dense matter and how this was discovered

equation of state of dense matter

the cold relation P(ε) from ρ ≈ 10⁶ to > 10¹⁵ g cm⁻³ · dense-matter and neutron-star physics

The idea that cold matter resists compression by quantum pressure alone reached its first crisis in 1931, when the nineteen-year-old Subrahmanyan Chandrasekhar showed that special relativity softens the electron gas to $\Gamma=4/3$ and imposes a maximum mass — a result Eddington famously refused to accept. The discovery of the neutron by James Chadwick in 1932 opened a denser chapter: Lev Landau speculated about stars made of a single gigantic nucleus, and Baade and Zwicky proposed in 1934 that neutron stars form in supernovae. The first genuine dense-matter EOS calculation was Oppenheimer and Volkoff's 1939 integration of the relativistic hydrostatic equation — now the Tolman–Oppenheimer–Volkoff equation — for a free neutron gas, giving a discouraging $M_\mathrm{max}\approx0.7\,M_\odot$ and proving the strong force could not be neglected. The crust EOS was placed on firm footing by Baym, Pethick and Sutherland in 1971, and the modern programme — connecting laboratory nuclear data (symmetry energy, its slope $L$ , three-body forces) to neutron-star observables — was crystallized by Lattimer and Prakash. What began as an idealized Fermi gas is now a precision inverse problem: $2\,M_\odot$ pulsars, NICER's mass–radius points, and GW170817's tidal bound together carve the allowed $P(\varepsilon)$ curve to a band, closing the loop Chandrasekhar and Oppenheimer opened.

History

1931

— Chandrasekhar derives the relativistic degeneracy mass limit for white dwarfs

1932

— Chadwick discovers the neutron; Landau speculates on collapsed dense stars

1939

— Oppenheimer and Volkoff compute the maximum mass of a free-neutron star

1971

— Baym, Pethick and Sutherland construct the crust equation of state

2017

— GW170817 places the first tidal-deformability bound on the equation of state

Going further

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., Ch. 16 (Degenerate Remnants) and Ch. 17 (Neutron Stars) · Shapiro & Teukolsky, Black Holes, White Dwarfs, and Neutron Stars (1983), Ch. 2, 8–9 · Glendenning, Compact Stars (2000) · Haensel, Potekhin & Yakovlev, Neutron Stars 1: Equation of State and Structure (2007) · Lattimer & Prakash, "Neutron Star Observations: Prognosis for Equation of State Constraints", Phys. Rep. 442 (2007) 109 · Oppenheimer & Volkoff, Phys. Rev. 55 (1939) 374. The polytropic constants $K_\mathrm{NR}, K_\mathrm{UR}$ and the exact degeneracy function $f(x)=x(2x^2-3)\sqrt{x^2+1}+3\,\mathrm{arcsinh}\,x$ are Chandrasekhar's classical results; their limiting forms are $f\to\tfrac{8}{5}x^5$ for $x\ll1$ and $f\to2x^4$ for $x\gg1$ , recovering the $\Gamma=5/3$ and $\Gamma=4/3$ polytropes respectively, and the mass scaling is $M\propto\rho_c^{(3\Gamma-4)/2}$ . Nuclear saturation $\rho_0\simeq2.7\times10^{14}\ \mathrm{g\,cm^{-3}}$ ( $n_0\simeq0.16\ \mathrm{fm^{-3}}$ ), neutron-drip $\approx4\times10^{11}\ \mathrm{g\,cm^{-3}}$ , symmetry energy $\approx30$ – $34\ \mathrm{MeV}$ , and the free-gas $M_\mathrm{max}\approx0.7\,M_\odot$ (Oppenheimer & Volkoff 1939) are standard values. The worked results $x\approx0.35$ at $n_0$ and $\Lambda\approx390$ for a $1.4\,M_\odot$ , $12\ \mathrm{km}$ , $k_2=0.09$ star follow from the equations above; pulsar masses, NICER radii ( $11$ – $13\ \mathrm{km}$ ), and the GW170817 bound $\Lambda_{1.4}\lesssim800$ are drawn from the observational literature.

Prerequisites. Electron Degeneracy Pressure · Chandrasekhar Mass Limit · Tolman–Oppenheimer–Volkoff Equation · Neutron Stars. Leads to. Neutron Star Structure and Interior · Neutron Stars · Pulsar Glitches and Neutron Star Superfluidity · Magnetars.

Ask the tutor: Ask me to trace a pulsar-mass or NICER-radius measurement back to a constraint on the $P(\varepsilon)$ curve, or to work the degeneracy-pressure integral to the exact $f(x)$ by hand.

Why is the equation of state called "cold" when neutron stars form at 10¹¹ K?

Because what matters is temperature relative to the Fermi energy, not its absolute value. In a neutron-star core the Fermi energy of the degenerate neutrons is hundreds of MeV, equivalent to $\sim10^{12}\ \mathrm{K}$ , while even a hot young star cools below $10^{9}$ – $10^{10}\ \mathrm{K}$ within seconds to years. The ratio $k_BT/E_F$ is tiny, so thermal excitations occupy only a thin sliver near the Fermi surface and the matter behaves as if $T=0$ : the EOS is a single curve $P(\varepsilon)$ . The exception is the first seconds of a supernova or a merger remnant, where finite-temperature corrections matter and a full $P(\varepsilon,T,Y_e)$ table is needed.

What exactly do "stiff" and "soft" mean for an equation of state?

Stiffness is how steeply pressure rises with energy density, $\mathrm{d}P/\mathrm{d}\varepsilon$ . A stiff EOS produces a lot of pressure for a given density, resists compression, and supports more mass in a larger star; a soft EOS gives less pressure, compresses easily, and yields smaller stars with a lower maximum mass. In polytropic language a larger index $\Gamma$ is stiffer, with $\Gamma=4/3$ the marginal value below which no stable star exists. Stiffness is not one number for the whole star — the EOS can be soft at one density and stiff at another — which is why different observations (radius, maximum mass, tidal deformability) constrain the curve at different densities.

Why can't we just compute the equation of state from QCD?

In principle the strong interaction — quantum chromodynamics — governs everything above $\rho_0$ , but it is not solvable in the relevant regime. Lattice QCD suffers a "sign problem" at finite baryon density that makes direct simulation intractable, while perturbative QCD works only at densities far above those in neutron-star cores. In between — the few-times-saturation regime that actually sets neutron-star structure — one must use effective theories: chiral effective field theory anchored to nucleon–nucleon scattering data up to about $2\rho_0$ , then phenomenological or interpolated models beyond. This gap is precisely why astronomical measurements are indispensable inputs, not just consistency checks.

What is the "hyperon puzzle"?

As density climbs in the core, it can become energetically favourable to convert some neutrons into heavier strange baryons — hyperons such as the Λ — because doing so relieves the enormous neutron Fermi pressure. But adding a species softens the EOS: the freed-up Fermi energy goes into rest mass rather than pressure, lowering the maximum mass, often below the observed $2\,M_\odot$ . That is the puzzle — naive hyperonic models are ruled out by the very existence of the heaviest pulsars. Resolutions invoke strongly repulsive hyperon–nucleon and three-body forces that keep the EOS stiff, push the hyperon onset to densities never reached, or a quark-matter transition that intervenes first. It remains open.

How do NICER and GW170817 constrain different parts of the curve?

They probe overlapping but distinct density ranges. NICER measures the radius of a star of known mass, and the radius of a $\sim1.4\,M_\odot$ star is set mainly by the pressure at one-to-two times saturation — it is most sensitive to the symmetry-energy slope $L$ . GW170817's tidal deformability also responds to the radius (scaling as roughly $R^5$ ), constraining a similar density range but with a sharp fifth-power lever that excludes large radii. The two-solar-mass pulsars, by contrast, constrain the high-density stiffness that sets $M_\mathrm{max}$ . Together — floor from pulsars, radius from NICER and GW170817 — they bracket the EOS across the full supranuclear range.

Is there an absolute maximum mass for a neutron star?

Yes, set by causality. The sound speed cannot exceed light, which caps $\mathrm{d}P/\mathrm{d}\varepsilon\le c^2$ and hence how stiff the EOS can be. Combining the stiffest causal EOS ( $P=\varepsilon-\varepsilon_0$ ) with the TOV equation, anchored at low density to known nuclear physics, gives an absolute bound of roughly $3\,M_\odot$ for a non-rotating star; rapid rotation can raise it by perhaps 20 percent. Any compact object heavier than this cannot be supported by matter and must be a black hole — which is why the $\sim2.6\,M_\odot$ object in GW190814 sits in an ambiguous "mass gap". The precise value depends on the true EOS, so tightening it is a payoff of the whole programme.