Equation of State for Stellar Matter

A star is a self-gravitating ball of gas that survives by winning, moment to moment, a standoff between gravity pulling inward and pressure pushing out. The equations that govern that balance — hydrostatic equilibrium, mass continuity, energy generation and transport — carry one more unknown than there are equations. The missing relation is the equation of state: the thermodynamic law that returns the pressure $P$ once you specify the local density $\rho$ , temperature $T$ , and composition. Without it the structure equations cannot be closed and no stellar model can be built. What makes stellar matter interesting is that no single formula works everywhere. In a main-sequence core the gas is a nearly ideal, fully ionized plasma and $P$ tracks $\rho T$ . In a massive star's centre the photon gas adds a $\tfrac13 aT^4$ term that rivals it and softens the medium. And in a white dwarf the electrons are squeezed so hard that the Pauli principle — not heat — supplies the pressure, which then depends on density alone and, in the relativistic limit, so weakly that a maximum mass appears from nowhere. This article assembles those pieces into one $P(\rho,T,\text{composition})$ , derives the exponents that decide stability, and shows how the balance among gas, radiation, and degeneracy maps onto a track through the density–temperature plane that every star follows from birth to remnant.

Key takeaways

The equation of state (EOS) closes the equations of stellar structure by giving pressure $P$ as a function of density $\rho$ , temperature $T$ , and composition.
In most interiors three terms add: an ideal gas of ions and electrons, radiation pressure $P_{\mathrm{rad}}=\tfrac13 aT^4$ , and — at high density — electron degeneracy pressure.
The ideal-gas term is $P_{\mathrm{gas}}=\rho k_B T/(\mu m_H)$ ; the mean molecular weight $\mu$ counts free particles per unit mass.
When electrons degenerate, $P$ stops depending on $T$ : $P\propto\rho^{5/3}$ non-relativistically, $P\propto\rho^{4/3}$ ultra-relativistically — the physics of white dwarfs and the Chandrasekhar mass.
Radiation pressure grows as $T^4$ and softens the EOS toward the marginal value $\Gamma_1\to 4/3$ , dominating in the most massive stars.
Key idea: pressure in a star is the sum of gas, radiation, and degeneracy contributions; which term dominates — set by where $(\rho,T)$ falls in the density–temperature plane — controls the star's structure, stability, and fate.

The core insight

Pressure in a star is the sum of gas, radiation, and degeneracy contributions; which term dominates — set by where $(\rho,T)$ falls in the density–temperature plane — controls the star's structure, stability, and evolutionary fate.

Where this lives today

The stellar EOS is live, precision physics. Modern codes such as MESA interpolate tabulated equations of state — the OPAL and SCVH tables and the more recent FreeEOS — that fold in partial ionization, Coulomb corrections, and degeneracy across the whole plane a star traverses. Helioseismology has made this a measurement: the Sun's oscillation frequencies, sensitive to $\Gamma_1$ through the sound speed $c_s^2=\Gamma_1 P/\rho$ , pin down the solar EOS to better than a part in a thousand. At the high-density frontier the Chandrasekhar mass remains the organizing number of stellar death — it decides which stars leave white dwarfs and which collapse, and the standard-candle uniformity of Type Ia supernovae rests on white dwarfs approaching that one critical mass. Open questions remain: the carbon ignition density, Coulomb crystallization in the white-dwarf "cooling delay" now dated in globular clusters, and the sub-Chandrasekhar detonation channel all turn on the EOS at densities no laboratory can reach.

The physical picture

Pressure is momentum flux: each particle crossing an imaginary surface carries momentum, and the pressure is the rate of that transport per unit area. That single fact explains why $P$ depends both on how many particles there are and on how fast they move — whether "fast" means thermal speed set by temperature or the Fermi speed forced on electrons by quantum crowding. Everything specific about stellar matter is a statement about which effect sets the speeds. In an ordinary gas, heat sets the speed: warm it and $P$ rises, cool it and it deflates. Degenerate matter breaks that link — pack electrons densely enough and the Pauli principle forces them into high-momentum states even at zero temperature, so their pressure survives cooling to $T\to0$ . A white dwarf is literally a stone held up by quantum mechanics rather than heat. The cleanest way to hold this in mind is the density–temperature plane, which divides into four provinces — ideal gas, radiation-dominated, non-relativistic degenerate, and relativistic degenerate — mapped in the first figure. Every star is a moving point here, and its whole life is a track wandering among these provinces.

Figure

The density–temperature plane and its four regimes, separated by T ∝ ρ²ᐟ³ for the onset of degeneracy, ρ ∝ T³ for the gas–radiation balance, and the vertical line ρ ≈ μ_e × 10⁶ g cm⁻³ where electrons turn relativistic; the Sun sits in the ideal-gas wedge, a white dwarf in the degenerate corner

The number that decides whether a star is stable is the stiffness of its EOS — how steeply the pressure climbs when you squeeze the gas, measured by $\Gamma_1=(\partial\ln P/\partial\ln\rho)_s$ . Stiff matter ( $\Gamma_1>4/3$ ) fights back hard and keeps the star in equilibrium; soft matter, pushed toward $\Gamma_1\to 4/3$ by radiation or relativistic degeneracy, sits on a knife-edge where a small squeeze barely raises $P$ and the star is one nudge from collapse. The second figure shows this directly, plotting the pressure of cold degenerate matter against density and marking where the exponent slides from $5/3$ down to the dangerous $4/3$ .

Figure

Pressure of cold degenerate electrons versus density on log scales: the curve rises with slope 5/3 in the non-relativistic regime, P ≈ 1.0×10¹³ (ρ/μ_e)⁵ᐟ³, then softens to slope 4/3 above ρ ≈ μ_e × 10⁶ g cm⁻³ where electrons turn relativistic; the amber ideal-gas isotherm at T = 10⁷ K sits far below, so degeneracy dominates thermal pressure here

The core equations

In the bulk of a normal star the pressure is the sum of an ideal-gas term from the ions and electrons and a radiation term from the photon field. This is the workhorse EOS for main-sequence and giant interiors, and the two terms scale so differently with $\rho$ and $T$ that their competition alone organizes much of stellar structure.

P = P_{\mathrm{gas}} + P_{\mathrm{rad}} = \dfrac{\rho k_B T}{\mu m_H} + \dfrac{1}{3}aT^4 \qquad (1)

where: $P$ — total pressure; $\rho$ — mass density; $T$ — temperature; $k_B$ — Boltzmann constant, $1.381\times10^{-16}\,\mathrm{erg\,K^{-1}}$ ; $m_H$ — hydrogen mass, $1.6726\times10^{-24}\,\mathrm{g}$ ; $\mu$ — mean molecular weight (mass per free particle in units of $m_H$ ); $a$ — radiation constant, $7.566\times10^{-15}\,\mathrm{erg\,cm^{-3}\,K^{-4}}$ . Reading it: the gas term is linear in $\rho$ and $T$ ; the radiation term ignores $\rho$ and explodes as $T^4$ . The gas term hides its microphysics inside $\mu$ , which counts how many free particles each gram of matter supplies once the atoms are ionized. Because ions and electrons contribute separately, $\mu$ splits into an ion part and an electron part that add as reciprocals — a small algebra that determines the pressure a plasma of given composition can muster.

\dfrac{1}{\mu} = \dfrac{1}{\mu_{\mathrm{ion}}} + \dfrac{1}{\mu_e},\qquad \dfrac{1}{\mu_{\mathrm{ion}}}\approx X + \dfrac{Y}{4} + \dfrac{Z}{\langle A\rangle},\qquad \dfrac{1}{\mu_e}\approx X + \dfrac{Y}{2} + \dfrac{Z}{2} \qquad (2)

where: $X$ — hydrogen mass fraction; $Y$ — helium mass fraction; $Z$ — metal mass fraction ( $X+Y+Z=1$ ); $\mu_{\mathrm{ion}}$ — mean molecular weight of ions; $\mu_e$ — mean molecular weight per free electron; $\langle A\rangle$ — mean atomic mass of the metals. Reading it: hydrogen is generous (two particles per proton mass), He and metals stingy, so $\mu$ climbs as H is processed into heavier nuclei. The balance between the two terms in equation (1) is captured by one dimensionless number, Eddington's $\beta$ , the fraction of the total pressure carried by gas — the natural variable in which the adiabatic exponent is expressed.

\beta \equiv \dfrac{P_{\mathrm{gas}}}{P},\qquad 1-\beta = \dfrac{P_{\mathrm{rad}}}{P} \qquad (3)

where: $\beta$ — gas fraction of the pressure; $\beta\to1$ for low-mass stars where gas dominates, $\beta\to0$ for very massive stars where radiation takes over. Reading it: $\beta$ is essentially fixed along the main sequence by the stellar mass through Eddington's quartic, $1-\beta\propto M^2\mu^4\beta^4$ . When the electrons stop being a classical gas, equation (1) fails and the honest expression is the full Fermi–Dirac integral for the electron pressure. Every regime the electrons can occupy — classical, non-relativistic degenerate, relativistic degenerate — is a limit of this one integral.

P_e = \dfrac{8\pi}{3h^3}\int_0^{\infty} \dfrac{p^3\,v(p)}{e^{(E-\mu_F)/k_BT}+1}\,\mathrm{d}p \qquad (4)

where: $P_e$ — electron partial pressure; $h$ — Planck constant; $p$ — electron momentum; $v(p)$ — electron speed, $=pc^2/E$ relativistically; $E$ — electron energy including rest mass; $\mu_F$ — chemical potential (the Fermi level). Reading it: the Fermi–Dirac factor $\to$ Maxwell–Boltzmann when sparsely filled (giving $P_e=n_e k_B T$ ) and $\to$ a filled step at $T\to0$ (degeneracy pressure). The frontier between the non-relativistic and relativistic degenerate regimes is set by when the Fermi energy reaches the electron rest energy. Beyond this density the electrons move at nearly the speed of light and the EOS softens decisively.

E_F/(m_ec^2) \sim 1 \;\Leftrightarrow\; \rho_{\mathrm{rel}} \approx \mu_e \times 10^{6}\,\mathrm{g\,cm^{-3}} \qquad (5)

where: $E_F$ — Fermi energy of the electrons; $\rho_{\mathrm{rel}}$ — density at which they turn relativistic. Reading it: for a CO white dwarf ( $\mu_e=2$ ) this is $\approx2\times10^{6}\,\mathrm{g\,cm^{-3}}$ , inside the range of massive white dwarfs — which is why the Chandrasekhar limit is reachable. Beyond $P$ itself, the EOS supplies the thermodynamic derivatives that govern stability and convection. Chandrasekhar's three adiabatic exponents measure how $P$ , $T$ , and $\rho$ respond to one another under an adiabatic squeeze.

\Gamma_1=\left(\dfrac{\partial\ln P}{\partial\ln\rho}\right)_s,\quad \dfrac{\Gamma_2-1}{\Gamma_2}=\left(\dfrac{\partial\ln T}{\partial\ln P}\right)_s=\nabla_{\mathrm{ad}},\quad \Gamma_3-1=\left(\dfrac{\partial\ln T}{\partial\ln\rho}\right)_s \qquad (6)

where: $\Gamma_1,\Gamma_2,\Gamma_3$ — Chandrasekhar's adiabatic exponents, all equal to $5/3$ for a monatomic ideal gas; $\nabla_{\mathrm{ad}}$ — adiabatic temperature gradient, $=0.4$ for an ideal gas; $s$ — specific entropy, held fixed. Reading it: $\Gamma_1$ enters dynamical stability, $\nabla_{\mathrm{ad}}$ enters the Schwarzschild and Ledoux convection criteria. It is often convenient to isolate the two ways pressure responds to state, at fixed $T$ and at fixed $\rho$ . These exponents feed the stability analysis and diagnose which physics is in charge.

\chi_\rho=\left(\dfrac{\partial\ln P}{\partial\ln\rho}\right)_T,\qquad \chi_T=\left(\dfrac{\partial\ln P}{\partial\ln T}\right)_\rho \qquad (7)

where: $\chi_\rho$ — isothermal density exponent, $=1$ for an ideal gas; $\chi_T$ — temperature exponent, $=1$ for gas, $=4$ for pure radiation, $=0$ for fully degenerate matter. Reading it: a fully degenerate gas has $\chi_T=0$ — pressure is blind to temperature — which is the mathematical signature of quantum support.

Derivation ANon-relativistic electron degeneracy pressure

The goal is to obtain the cold, temperature-independent pressure $P\propto\rho^{5/3}$ that holds up a white dwarf, from nothing but the Pauli principle and phase-space counting. At $T=0$ every state below the Fermi momentum $p_F$ is occupied and every state above it empty; we count how many electrons that Fermi sphere holds, find $p_F$ as a function of density, and compute the pressure as the momentum flux those electrons carry — arriving at a pressure that depends on density alone and survives to zero temperature.

Step 1. Count states. Phase space is divided into cells of volume $h^3$ , each holding two electrons (spin up and down), so a Fermi sphere of radius $p_F$ holds $n_e=\tfrac{2}{h^3}\cdot\tfrac43\pi p_F^3$ . Why: the factor $2/h^3$ is the Pauli principle made quantitative — it forces electrons to high momenta once space runs short.

Step 2. Invert for the Fermi momentum, $p_F=(3h^3 n_e/8\pi)^{1/3}=\hbar(3\pi^2 n_e)^{1/3}$ , with $n_e=\rho/(\mu_e m_H)$ . Why: this ties the quantum quantity $p_F$ to the astrophysical variable $\rho$ , so the pressure will emerge as a function of density.

Step 3. Write the pressure as one-third of the momentum-flux integral over the filled sphere, $P_e=\tfrac{8\pi}{3h^3}\int_0^{p_F}p^3 v\,\mathrm{d}p$ , with $v=p/m_e$ . Why: pressure is momentum flux, and the one-third is the isotropic average over directions of motion.

Step 4. Evaluate: with $v=p/m_e$ the integrand is $p^4/m_e$ and $\int_0^{p_F}p^4\,\mathrm{d}p=p_F^5/5$ , so $P_e=\tfrac{8\pi}{15 m_e h^3}p_F^5$ . Why: the fifth power of $p_F$ is the origin of the steep $\rho^{5/3}$ dependence.

Step 5. Substitute $p_F(n_e)$ to collapse the constants into $P_e=\dfrac{(3\pi^2)^{2/3}}{5}\dfrac{\hbar^2}{m_e}\left(\dfrac{\rho}{\mu_e m_H}\right)^{5/3}$ . Why: the result carries no $T$ , so a white dwarf can cool toward $T\to0$ without losing support — the defining property of degenerate matter. *The coefficient evaluates to $K_1=(3\pi^2)^{2/3}\hbar^2/(5m_e)(\mu_e m_H)^{-5/3}=1.0\times10^{13}\,\mu_e^{-5/3}$ cgs, and the corresponding polytropic index is $n=3/2$ .*

Go deeper · Derivation B · 1 minThe adiabatic exponent of a gas plus radiation mixture

Derivation BThe adiabatic exponent of a gas plus radiation mixture

The goal is to show how the stiffness $\Gamma_1$ of a mixed gas-and-radiation medium interpolates between the safe ideal-gas value $5/3$ and the marginal value $4/3$ , as a function of the single parameter $\beta=P_{\mathrm{gas}}/P$ . This is the quantitative statement of why very massive stars, in which radiation carries much of the pressure, live near the edge of dynamical instability — we set up the two-component pressure, find how each responds to compression and heating, and assemble the exponent.

Step 1. Set up the mixture: $P=P_{\mathrm{gas}}+P_{\mathrm{rad}}$ with $P_{\mathrm{gas}}=\rho\Re T/\mu$ (monatomic, $\gamma=5/3$ ) and $P_{\mathrm{rad}}=\tfrac13 aT^4$ , and $\beta=P_{\mathrm{gas}}/P$ . Why: the two terms depend on $\rho$ and $T$ so differently that the total's response to a squeeze is a weighted blend, and $\beta$ is the weight.

Step 2. Read off the isothermal exponents: compressing at fixed $T$ moves only the gas, so $\chi_\rho=\beta$ ; heating at fixed $\rho$ moves gas as $T^1$ and radiation as $T^4$ , so $\chi_T=4-3\beta$ . Why: the radiation term's factor of four is what softens the response as $\beta$ falls.

Step 3. Assemble the exponent from the specific heats: $\Gamma_1=\beta+\dfrac{(4-3\beta)^2(\gamma-1)}{\beta+12(\gamma-1)(1-\beta)}$ with $\gamma=\tfrac53$ . Why: this closed form turns the composition of the pressure into a stability diagnostic in one line.

Step 4. Check limits: $\beta=1$ gives $\Gamma_1=5/3$ ; $\beta\to0$ gives $\Gamma_1\to4/3$ ; every mixture lies between, and $\nabla_{\mathrm{ad}}$ falls from $0.4$ toward $0.25$ as radiation takes over. Why: this softening is why radiation-heavy envelopes convect readily and the most luminous stars flirt with collapse. *As a check on the formula: setting $\beta=1$ returns $\Gamma_1=5/3$ and taking $\beta\to0$ returns $\Gamma_1\to4/3$ , confirming the two limiting cases.*

Go deeper · Derivation C · 1 minThe Chandrasekhar mass from the relativistic equation of state

Derivation CThe Chandrasekhar mass from the relativistic equation of state

The goal is to explain why the ultra-relativistic degenerate EOS, $P\propto\rho^{4/3}$ , does not merely soften a star but imposes an absolute ceiling on the mass it can support — the Chandrasekhar limit. The structural fact is that a $\gamma=4/3$ EOS is a polytrope of index $n=3$ , whose mass is fixed entirely by physical constants with the central density scaling out. That independence of central density is the fingerprint of a unique critical mass; above it, no equilibrium exists.

Step 1. Recognize the polytrope: in the limit $v\to c$ the Fermi-integral pressure becomes $P=K_2\rho^{4/3}=K_2\rho^{1+1/n}$ with $n=3$ and $K_2=\dfrac{(3\pi^2)^{1/3}\hbar c}{4}\left(\dfrac{1}{\mu_e m_H}\right)^{4/3}$ . Why: casting the EOS as a polytrope lets us borrow the Lane–Emden solution, and $n=3$ is what makes the mass special.

Step 2. Write the polytropic mass: for $n=3$ , $M=4\pi\left(\dfrac{K_2}{\pi G}\right)^{3/2}(-\xi^2\theta')_{\xi_1}$ with the surface constant $(-\xi^2\theta')_{\xi_1}=2.018$ . Why: the central density has dropped out of $M$ entirely — the star has this one mass at any central density, so equilibrium exists only here.

Step 3. Insert $K_2$ and the constants: $M_{\mathrm{Ch}}=\dfrac{5.83}{\mu_e^2}\,M_\odot$ . Why: the mass depends only on fundamental constants and $\mu_e$ — a universal ceiling, not a property of any one star.

Step 4. Read the physics: for a CO white dwarf $\mu_e=2$ , so $M_{\mathrm{Ch}}=5.83/4\approx1.46\,M_\odot$ (canonically quoted $\approx1.4\,M_\odot$ ). Above it the soft $\gamma=4/3$ EOS cannot stiffen fast enough, and the star collapses toward a neutron star or ignites as a Type Ia supernova. Why: this one number gates the fate of every white dwarf that approaches it. *The Lane–Emden surface constant $(-\xi^2\theta')_{\xi_1}=2.018$ for $n=3$ is what yields the numerical prefactor $5.83\,\mu_e^{-2}M_\odot$ .*

Limiting cases

The general EOS reduces to clean, memorable forms depending on which physics is in charge:

Case	Result	Meaning
Ideal non-degenerate gas	$P = \dfrac{\rho k_B T}{\mu m_H}$	Valid when the de Broglie wavelength is far below the spacing and $k_BT$ exceeds the Fermi energy; the workhorse for main-sequence interiors.
Radiation-dominated	$P \approx \tfrac{1}{3}aT^4$	Takes over in very massive stars where $\beta\to0$ ; drives $\Gamma_1\to4/3$ and sets the Eddington luminosity.
Non-relativistic degenerate	$P = K_1\rho^{5/3},\; K_1=\dfrac{(3\pi^2)^{2/3}\hbar^2}{5 m_e}\left(\dfrac{1}{\mu_e m_H}\right)^{5/3}$	Supports low-mass white dwarfs and He/C cores; independent of $T$ .
Ultra-relativistic degenerate	$P = K_2\rho^{4/3},\; K_2=\dfrac{(3\pi^2)^{1/3}\hbar c}{4}\left(\dfrac{1}{\mu_e m_H}\right)^{4/3}$	Massive white-dwarf cores; the $\rho^{4/3}$ softness admits no equilibrium above $M_{\mathrm{Ch}}$ .

Common trapThe ideal-gas law is not universal, and μ is not a constant

The most common error in stellar-interior estimates is to use $P=\rho k_B T/(\mu m_H)$ everywhere. It fails in two opposite directions. Where electrons are degenerate — white-dwarf interiors, red-giant cores — the true pressure comes from Fermi momenta, can exceed the ideal estimate by orders of magnitude, and is nearly $T$ -independent; the ideal law there lets a cooling core collapse in your model when the real star is stable. Where radiation matters — massive-star centres — the ideal law understates $P$ and misses the softening of $\Gamma_1$ toward $4/3$ that governs stability. And even in its proper domain $\mu$ is not fixed: it depends on ionization and composition ( $\mu_{\mathrm{ion}}\neq\mu_e$ ) and drifts upward as a star fuses hydrogen into heavier elements.

Worked examples

1 · Mean molecular weight of solar-composition plasma (confidence builder)

Take a fully ionized plasma with the Sun's composition, $X=0.70$ , $Y=0.28$ , $Z=0.02$ , and mean metal mass $\langle A\rangle\approx16$ . Find $\mu_{\mathrm{ion}}$ , $\mu_e$ , and $\mu$ .

Step 1. Ion part from equation (2): $1/\mu_{\mathrm{ion}}=X+Y/4+Z/\langle A\rangle=0.70+0.070+0.00125=0.771$ , so $\mu_{\mathrm{ion}}=\mathbf{1.30}$ . A helium nucleus is four mass units but one particle, which is why He and metals depress the count.

Step 2. Electron part: $1/\mu_e=X+Y/2+Z/2=0.70+0.140+0.010=0.850$ , so $\mu_e=\mathbf{1.18}$ . H liberates one electron per mass unit; He and metals roughly one per two.

Step 3. Combine: $1/\mu=0.771+0.850=1.621$ , giving $\mu=\mathbf{0.62}$ — the familiar $\mu\approx0.6$ for solar gas. As the core burns H to He, $\mu$ rises toward $\approx0.84$ and the gas pressure at fixed $\rho,T$ drops.

2 · Central pressure of the Sun: gas versus radiation (fully worked)

Using solar-centre values $\rho_c\approx150\,\mathrm{g\,cm^{-3}}$ , $T_c\approx1.5\times10^{7}\,\mathrm{K}$ , and the evolved-core $\mu\approx0.84$ , estimate $P_{\mathrm{gas}}$ , $P_{\mathrm{rad}}$ , and $\beta$ .

Step 1. Gas pressure from equation (1), with $\Re=k_B/m_u=8.314\times10^{7}\,\mathrm{erg\,g^{-1}\,K^{-1}}$ : $P_{\mathrm{gas}}=\rho_c\Re T_c/\mu=(150)(8.314\times10^{7})(1.5\times10^{7})/0.84=\mathbf{2.2\times10^{17}\,dyn\,cm^{-2}}$ , right on the value from detailed solar models.

Step 2. Radiation pressure: $P_{\mathrm{rad}}=\tfrac13 aT_c^4=\tfrac13(7.566\times10^{-15})(1.5\times10^{7})^4=\mathbf{1.3\times10^{14}\,dyn\,cm^{-2}}$ — the $T^4$ is steep, but at $1.5\times10^7\,\mathrm{K}$ still tiny.

Step 3. Gas fraction: $1-\beta=P_{\mathrm{rad}}/P\approx\mathbf{5.8\times10^{-4}}$ , so $\beta\approx0.9994$ — radiation supplies under 0.1% of the Sun's central pressure, a firmly gas-dominated, $\Gamma_1\approx5/3$ star. A $30\,M_\odot$ star has $\beta$ nearer $0.7$ , and the most massive stars push $\beta$ below $0.5$ .

3 · Degeneracy versus thermal pressure in a white dwarf (fully worked)

A carbon white-dwarf interior has $\rho=10^{6}\,\mathrm{g\,cm^{-3}}$ , $\mu_e=2$ , and $T\approx10^{7}\,\mathrm{K}$ . Compare the electron degeneracy pressure with the pressure the same electrons would exert as an ideal thermal gas, and decide whether the star cares about its temperature.

Step 1. Degeneracy pressure from the non-relativistic limit, $P_{\mathrm{deg}}=K_1(\rho/\mu_e)^{5/3}$ with $K_1=1.0\times10^{13}$ cgs: $(\rho/\mu_e)=5\times10^{5}$ , so $(\rho/\mu_e)^{5/3}=3.15\times10^{9}$ and $P_{\mathrm{deg}}=\mathbf{3.1\times10^{22}\,dyn\,cm^{-2}}$ . This density is only mildly relativistic — the onset (equation 5) is $\rho_{\mathrm{rel}}\approx2\times10^{6}\,\mathrm{g\,cm^{-3}}$ — so the $5/3$ law is a good first estimate.

Step 2. Ideal thermal pressure of the same plasma (fully ionized carbon, $\mu=12/7\approx1.71$ ): $P_{\mathrm{th}}=\rho\Re T/\mu=(10^{6})(8.314\times10^{7})(10^{7})/1.71=\mathbf{4.8\times10^{20}\,dyn\,cm^{-2}}$ .

Step 3. Compare: $P_{\mathrm{deg}}/P_{\mathrm{th}}\approx\mathbf{65}$ . Degeneracy outweighs thermal pressure by about 65×, so the star is held up by quantum crowding and is almost blind to temperature. That is why a white dwarf can cool for billions of years and never collapse: shed the entire $4.8\times10^{20}$ thermal term and you have barely dented the $3.1\times10^{22}$ that supports it.

Physics of the portrait

The equation of state for stellar matter is the biography of one insight: that quantum statistics, not heat, holds up dead stars — and that this support has a limit. Both halves belong to one man.

Origins · 1 minSubrahmanyan Chandrasekhar and how this was discovered

Subrahmanyan Chandrasekhar

1910–1995, Nobel Prize in Physics 1983 · degenerate stellar matter and the limiting mass of white dwarfs

Subrahmanyan Chandrasekhar was nineteen and aboard a ship from Madras to Cambridge in 1930 when he combined Fowler's idea that white dwarfs are degenerate electron gases with special relativity, and found something no one expected: a maximum mass. Ralph Fowler had shown in 1926 that the Pauli principle supplies the pressure holding up a white dwarf even at zero temperature — but Fowler used the non-relativistic gas, which supports any mass. Chandrasekhar noticed that in a dense enough star the electrons must move near the speed of light, softening the EOS from $P\propto\rho^{5/3}$ to $P\propto\rho^{4/3}$ , and that a $\rho^{4/3}$ gas is an $n=3$ polytrope whose mass is fixed by fundamental constants alone at about $1.4\,M_\odot$ ; above it, degeneracy pressure cannot win. Eddington ridiculed the result publicly in 1935, insisting "there should be a law of Nature to prevent a star from behaving in this absurd way." Chandrasekhar was right; the limit now bears his name, gates the formation of neutron stars and black holes, and underwrites the Type Ia supernovae used to measure cosmic acceleration. He also systematized the thermodynamics of stellar matter — the exponents $\Gamma_1,\Gamma_2,\Gamma_3$ still carry his notation — and shared the 1983 Nobel Prize for work begun on that voyage.

History

1926

— Eddington's Internal Constitution of the Stars formalizes radiation pressure and the quartic for β

1926

— Fowler shows that white-dwarf matter is a degenerate electron gas held up by the Pauli principle

1930

— Chandrasekhar derives the relativistic degeneracy softening en route to Cambridge

1931

— Chandrasekhar publishes the maximum mass for white dwarfs, about 1.4 solar masses

1939

— Chandrasekhar's monograph systematizes the adiabatic exponents and the stellar equation of state

Going further

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., ch. 10 and ch. 16 (stellar interiors and degenerate remnants) · Kippenhahn, Weigert & Weiss, Stellar Structure and Evolution, 2nd ed., ch. 13–16 · Clayton, Principles of Stellar Evolution and Nucleosynthesis, ch. 2 · Chandrasekhar, An Introduction to the Study of Stellar Structure (1939), ch. X–XI · Hansen, Kawaler & Trimble, Stellar Interiors, 2nd ed., ch. 3. For broader context, the standard field references are Dodelson (Modern Cosmology), Mo, van den Bosch & White (Galaxy Formation and Evolution), and Binney & Tremaine (Galactic Dynamics). Key numerical anchors used throughout this article: the degeneracy coefficient $K_1=(3\pi^2)^{2/3}\hbar^2/(5m_e)(\mu_e m_H)^{-5/3}=1.0\times10^{13}\,\mu_e^{-5/3}$ cgs with $\rho^{5/3}$ scaling; radiation constant $a=7.566\times10^{-15}$ and gas constant $\Re=8.314\times10^{7}$ cgs; Chandrasekhar mass $5.83\,\mu_e^{-2}M_\odot\to1.46\,M_\odot$ for $\mu_e=2$ , with $n=3$ Lane–Emden constant $2.018$ (canonical $\approx1.4\,M_\odot$ ); $\mu=0.5$ for pure ionized H and $\mu\approx0.62$ for solar composition; solar-centre $P\approx2.2\times10^{17}\,\mathrm{dyn\,cm^{-2}}$ , $\beta\approx0.9994$ , and white-dwarf degeneracy-to-thermal ratio $P_{\mathrm{deg}}/P_{\mathrm{th}}\approx65$ ; and the $\Gamma_1(\beta)$ limits running from $5/3$ down to $4/3$ .

Prerequisites. Hydrostatic Equilibrium in Stars · Mass Continuity Equation · Radiation Pressure in Stellar Interiors. Leads to. Equations of Stellar Structure · Adiabatic Temperature Gradient · Schwarzschild and Ledoux Criteria for Convection · Virial Theorem for Stars.

Ask the tutor: Ask me to walk you through where any real star — the Sun's core, a $30\,M_\odot$ supergiant, or a cooling white dwarf — sits in the density–temperature plane and which term of the equation of state is holding it up.

Why does the pressure become independent of temperature when matter degenerates?

Because the pressure then comes from momentum, not heat. In a classical gas the speeds are thermal, so cooling drops the pressure. In degenerate matter the electrons are packed so densely that the Pauli principle forces them into high-momentum states regardless of temperature — the Fermi momentum $p_F=\hbar(3\pi^2 n_e)^{1/3}$ depends only on density, so lowering $T$ toward zero barely changes $p_F$ or the pressure. That is why a white dwarf, radiating its residual heat with no fusion, does not collapse as it cools: its support was never thermal.

What is the mean molecular weight μ, and why is it not just the mass of a hydrogen atom?

$\mu$ is the average mass, in units of $m_H$ , carried by each free particle — and "free particle" is the whole point. Ionizing a hydrogen atom gives two free particles from one mass unit, so hydrogen contributes $\mu\approx0.5$ ; helium gives one nucleus plus two electrons, three particles from four mass units, so $\mu\approx1.33$ . Gas pressure depends on the number of particles, not their weight, so hydrogen-rich plasma has more particles per gram and more pressure at given $\rho$ and $T$ . As a star fuses H into He, $\mu$ climbs and the star adjusts.

Is the Chandrasekhar mass exactly 1.4 solar masses?

Not exactly, and the spread in quoted values is real physics. The ideal polytrope gives $M_{\mathrm{Ch}}=5.83\,\mu_e^{-2}M_\odot$ , which for $\mu_e=2$ is $1.46\,M_\odot$ . The canonical " $1.4$ " includes corrections the polytrope leaves out: finite temperature, electrostatic Coulomb interactions (which lower the pressure), general relativity, and composition — bringing real CO white dwarfs to an effective limit near $1.38$ – $1.40\,M_\odot$ . The $\mu_e^{-2}$ scaling means an iron white dwarf would have a lower limit. So " $1.4\,M_\odot$ " is a rounded, correction-included value of the leading-order $5.83\,\mu_e^{-2}M_\odot$ .

Why does radiation pressure make massive stars unstable?

Because it softens the EOS toward the marginal value $\Gamma_1=4/3$ . Dynamical stability requires the pressure-weighted average of $\Gamma_1$ to exceed $4/3$ ; at exactly $4/3$ the restoring force against uniform compression vanishes to leading order and gravity is no longer opposed. A pure monatomic gas has the comfortable $\Gamma_1=5/3$ , but radiation contributes with an effective exponent nearer $4/3$ , so as the radiation fraction $1-\beta$ grows in the most massive stars, $\Gamma_1$ is dragged toward the critical value. Such stars sit near the edge of instability, prone to violent pulsations and mass loss — one reason there appears to be an upper limit to stellar masses.

Where does the Fermi–Dirac integral reduce to the simple ideal-gas law?

In the non-degenerate limit, when the states are sparsely occupied — that is, $e^{(E-\mu_F)/k_BT}\gg1$ for the occupied states. Then the $+1$ in the Fermi–Dirac denominator is negligible, the occupation becomes the Maxwell–Boltzmann factor, and the integral collapses to $P_e=n_e k_BT$ , the ideal-gas law. Physically this is where the thermal de Broglie wavelength is much smaller than the interparticle spacing, so particles rarely contend for the same state. Main-sequence interiors live here; white-dwarf interiors need the full integral.

How do the adiabatic exponents connect the EOS to convection?

Through the adiabatic temperature gradient $\nabla_{\mathrm{ad}}=(\Gamma_2-1)/\Gamma_2$ , which the EOS fixes. A blob displaced upward cools along an adiabat at the rate $\nabla_{\mathrm{ad}}$ ; if the star's actual gradient is steeper it stays hotter and less dense than its surroundings and keeps rising — the Schwarzschild criterion $\nabla_{\mathrm{rad}}>\nabla_{\mathrm{ad}}$ (with a composition term in the Ledoux version). For an ideal gas $\nabla_{\mathrm{ad}}=0.4$ , but in partial-ionization zones the energy spent ionizing rather than heating drives $\Gamma_1$ and $\nabla_{\mathrm{ad}}$ down sharply — which makes the outer layers of cool stars convect and powers the pulsations of Cepheids and RR Lyrae stars.