Polytropic Models of Stars

Real stars obey four coupled differential equations — hydrostatic equilibrium, mass continuity, energy generation, and energy transport — closed by an equation of state and constitutive laws for opacity and nuclear burning, and in general nothing short of a computer will integrate them. Polytropic models cut this Gordian knot with a single audacious postulate: that pressure is a fixed power of density, $P = K\rho^{1+1/n}$ , so that the mechanical structure of a star closes off from its thermal structure entirely. The assumption is not a convenience plucked from the air — it is exact for a degenerate electron gas (both non-relativistic and ultra-relativistic), exact for a region stirred to uniform entropy by convection, and very nearly exact in the Eddington standard model where radiation and gas pressure keep a fixed ratio. The payoff is spectacular. Once $P$ is a known function of $\rho$ alone, hydrostatic balance and mass conservation fuse into one second-order ordinary differential equation — the Lane–Emden equation — for a dimensionless density $\theta$ versus a dimensionless radius $\xi$ . Every physical star of a given index $n$ is the same dimensionless solution, stretched by two scale factors, so solving it once describes a whole one-parameter family: white dwarfs, fully convective red dwarfs, isothermal cloud cores, and the interiors that seed modern evolution codes. The single number $n$ carries deep meaning. Small $n$ is stiff, incompressible matter; $n=1.5$ is non-relativistic degeneracy and adiabatic $\gamma=5/3$ gas; $n=3$ is ultra-relativistic degeneracy and $\gamma=4/3$ , and it delivers the most famous result in stellar physics — a mass fixed by fundamental constants alone, the Chandrasekhar limit near $1.4\,M_\odot$ above which no white dwarf can stand. Cross $n=5$ and the star becomes infinite; beyond it the mass diverges. Polytropes remain the first tool any theorist reaches for when a self-gravitating fluid must be understood by hand.

Key takeaways

A polytrope replaces a star's full thermal physics with one power law, $P=K\rho^{1+1/n}$ , decoupling the mechanical structure from temperature, luminosity, and opacity.
Hydrostatic equilibrium plus mass continuity then collapse into the dimensionless Lane–Emden equation, whose universal solution $\theta(\xi)$ fixes density, pressure, and mass; two scale factors turn it into a real star.
Closed-form solutions exist only for $n=0$ ( $\xi_1=\sqrt6$ ), $n=1$ ( $\theta=\sin\xi/\xi$ , $\xi_1=\pi$ ), and $n=5$ (infinite radius, finite mass); all others are integrated numerically.
The decisive indices are $n=1.5$ (fully convective stars and non-relativistic degeneracy, giving $MR^3=$ const) and $n=3$ (the Eddington model and ultra-relativistic degeneracy, giving the Chandrasekhar limiting mass).
Polytropes reproduce white-dwarf mass–radius relations and give central pressures and temperatures good to a factor of about two, without ever solving the energy equations.
Key idea: by assuming $P = K\rho^{1+1/n}$ the pressure decouples from temperature, so hydrostatic equilibrium plus mass continuity reduce to a single dimensionless equation — the Lane–Emden equation — whose solution characterizes an entire family of stars through just the index $n$ and two scale factors.

Where this lives today

Polytropes are working infrastructure, not a historical curiosity. Modern stellar-evolution codes such as MESA seed their pre-main-sequence models with $n=1.5$ polytropes, because a fully convective protostar has nearly uniform entropy and so is genuinely polytropic. In white-dwarf astrophysics the Chandrasekhar mass is the trigger of Type Ia supernovae: a carbon–oxygen white dwarf accreting toward $\sim1.4\,M_\odot$ ignites explosively, and because that threshold is nearly universal, Type Ia events are the standardizable candles that revealed the accelerating Universe. Neutron-star modelling still leans on piecewise-polytropic equations of state above nuclear density, and mass–radius constraints from NICER X-ray timing and gravitational-wave tidal deformabilities (notably GW170817) are expressed as bounds on polytropic stiffness. On the smallest scales the $n\to\infty$ isothermal sphere underlies the Bonnor–Ebert analysis of cloud cores. From the ignition of the Sun to the standard candles of cosmology, the same one-parameter family keeps reappearing.

The physical picture

Picture a ball of gas holding itself up against its own gravity by pressure alone. How centrally concentrated it becomes depends on one thing: how stiffly the gas pushes back when squeezed. If pressure rises steeply with density (small $n$ , large $\gamma$ ), the matter resists compression and the star is only mildly peaked. If pressure rises weakly (large $n$ , $\gamma$ toward $4/3$ ), the star piles a towering central density under a wispy envelope and teeters near instability. The index $n$ is a dial for how compressible the matter is; turning it from $0$ toward $5$ takes you from a nearly uniform sphere to a monstrously condensed one. The Lane–Emden variable $\theta$ measures how far the density has fallen from center to surface, through $\rho=\rho_c\theta^{\,n}$ , with $\theta=1$ at the center and $\theta=0$ at the surface, while $\xi$ measures distance in units of a natural length set by $K$ and $\rho_c$ . Solve the equation once and you have the shape of the star; the two scale factors stretch it to any size and mass. The first figure shows how violently the central concentration climbs with $n$ .

Figure

The Lane–Emden family θ(ξ) for n = 0, 1, 3: larger n means a more centrally condensed star. Each curve falls from θ = 1 at the center to its first zero ξ₁ (the surface) at ξ₁ = 2.449, π = 3.142, 6.897, and the central-to-mean density ratio ρc/ρ̄ climbs from 1 to 3.29 to 54.2

Why should this matter for real stars? Because the two most important indices are handed to us by microphysics, not chosen for convenience. A degenerate electron gas — the stuff of white dwarfs — obeys $P\propto\rho^{5/3}$ when the electrons are slow ( $n=1.5$ ) and $P\propto\rho^{4/3}$ when they are relativistic ( $n=3$ ). In the relativistic case $K$ is built from $h$ , $c$ , and the electron mass alone, so the total mass ceases to depend on how hard you squeeze it: an $n=3$ polytrope pinned to a single mass. That value is the Chandrasekhar mass, and it is why a white dwarf cannot exceed about $1.4\,M_\odot$ before collapsing. The second figure shows the payoff: the white-dwarf mass–radius relation bending toward that wall.

Figure

White-dwarf mass–radius relation from polytropes: non-relativistic degeneracy (n = 1.5) gives R ∝ M^(−1/3), so heavier white dwarfs are smaller; relativistic corrections (toward n = 3) drive the radius to zero at the Chandrasekhar mass M_Ch ≈ 1.44 M⊙ for μe = 2. Sirius B sits near 1.0 M⊙, 0.88 Earth radii

The dividing index $n=5$ is the frontier of bounded configurations. For $0\le n<5$ the density reaches zero at a finite $\xi_1$ — a genuine surface. At $n=5$ the density tail stretches to infinity even though, remarkably, the total mass stays finite; beyond $n>5$ the mass diverges and no such sphere can exist. Real stars live in $0\le n<5$ , with the interesting physics clustered around $n=1.5$ and $n=3$ .

The core equations

The whole edifice rests on one closure assumption: the polytropic relation makes pressure a pure function of density, so temperature never enters the mechanical structure.

P = K\,\rho^{1+1/n} = K\,\rho^{\gamma} \qquad (1)

where: $P$ — pressure; $\rho$ — mass density; $K$ — polytropic constant, fixed by the microphysics of the gas; $n$ — polytropic index; $\gamma$ — polytropic exponent, $\gamma=1+1/n$ . A small $n$ (large $\gamma$ ) is stiff, incompressible matter; a large $n$ (small $\gamma$ ) is soft, highly compressible matter. Feed this into hydrostatic equilibrium and mass continuity, strip out every dimensional constant, and the mechanical problem collapses to one universal equation for the shape function $\theta(\xi)$ .

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

where: $\theta$ — dimensionless density variable, $\rho=\rho_c\theta^{\,n}$ ; $\xi$ — dimensionless radius, $r=\alpha\xi$ ; $n$ — polytropic index. The right-hand side is self-gravity; the left is the divergence of the pressure gradient. Derivation A builds this line by line. Being second order, the equation needs two conditions, both imposed at the center where regularity demands a smooth, gradient-free density peak.

\theta(0)=1,\qquad \left.\frac{\mathrm{d}\theta}{\mathrm{d}\xi}\right|_{\xi=0}=0 \qquad (3)

where: $\theta(0)=1$ — density equals $\rho_c$ at the center; $\theta'(0)=0$ — zero gradient at the center, because the enclosed mass vanishes there and so must the gravitational force. Integrating outward from these two numbers generates the whole star. The dimensional information the equation discarded is carried by a single length scale $\alpha$ : multiply $\xi$ by $\alpha$ to recover physical distance, and the surface sits at the first zero $\xi_1$ of $\theta$ .

\alpha=\left[\frac{(n+1)K}{4\pi G}\,\rho_c^{\frac{1-n}{n}}\right]^{1/2},\qquad R=\alpha\,\xi_1 \qquad (4)

where: $\alpha$ — physical length for $\xi=1$ ; $\rho_c$ — central density; $G$ — gravitational constant; $R$ — stellar radius; $\xi_1$ — first zero of $\theta$ , the surface. All the star's individuality lives in $K$ and $\rho_c$ ; $\xi_1$ is tabulated once per index. The mass follows from integrating $4\pi r^2\rho$ , but the Lane–Emden equation turns the volume integral into a boundary term, so the mass depends only on the surface value of $\xi^2\theta'$ .

M=4\pi\alpha^3\rho_c\left(-\,\xi^2\frac{\mathrm{d}\theta}{\mathrm{d}\xi}\right)_{\xi_1} \qquad (5)

where: $M$ — total stellar mass; $(-\xi^2\theta')_{\xi_1}$ — dimensionless surface constant, equal to $4.899$ ( $n=0$ ), $\pi$ ( $n=1$ ), $2.714$ ( $n=1.5$ ), $2.018$ ( $n=3$ ), $\sqrt3$ ( $n=5$ ). Derivation B evaluates this and reads off the Chandrasekhar mass. Finally, the polytrope sharpens the crude hydrostatic pressure estimate into an exact scaling with a coefficient set by the index alone.

P_c = W_n\,\frac{G M^2}{R^4},\qquad W_n=\frac{1}{4\pi(n+1)\,[\theta'(\xi_1)]^2} \qquad (6)

where: $P_c$ — central pressure; $W_n$ — dimensionless coefficient, e.g. $W_3\approx11.05$ . The naive guess $P_c\sim GM^2/R^4$ understates the true coefficient by an order of magnitude. Derivation C assembles this and converts it to a central temperature.

Derivation AFrom two structure equations to one dimensionless equation

The goal here is to show how the two mechanical equations of stellar structure — hydrostatic equilibrium and mass continuity — fuse, under the polytropic closure $P=K\rho^{1+1/n}$ , into the one universal Lane–Emden equation whose solution depends on the index $n$ and nothing else. We eliminate the enclosed mass to reach a Poisson-like form, insert the polytrope, and choose a length scale $\alpha$ that erases every dimensional constant.

Step 1. Start from the two mechanical equations of a static spherical star: hydrostatic equilibrium $\mathrm{d}P/\mathrm{d}r=-Gm\rho/r^2$ , and mass continuity $\mathrm{d}m/\mathrm{d}r=4\pi r^2\rho$ . Multiply the first by $r^2/\rho$ , differentiate with respect to $r$ , and substitute the second to eliminate $m(r)$ , giving the Poisson form $\tfrac{1}{r^2}\tfrac{\mathrm{d}}{\mathrm{d}r}\!\big(\tfrac{r^2}{\rho}\tfrac{\mathrm{d}P}{\mathrm{d}r}\big)=-4\pi G\rho$ . Why: this single equation now contains only $P$ and $\rho$ , so the moment a relation between them is supplied the problem is closed.

Step 2. Impose the polytrope. Write $\rho=\rho_c\theta^{\,n}$ , which makes $\theta=1$ at the center where $\rho=\rho_c$ , and then $P=K\rho^{1+1/n}=K\rho_c^{1+1/n}\theta^{\,n+1}$ ; the leading coefficient is a constant carried outside all derivatives. Why: raising density to the power $n$ and pressure to the power $n+1$ is exactly what makes the powers of $\theta$ collapse cleanly in the next step — this is the payoff of defining $\theta$ through $\rho=\rho_c\theta^n$ rather than through $\rho$ directly.

Step 3. Compute the pressure-gradient term: $\tfrac{1}{\rho}\tfrac{\mathrm{d}P}{\mathrm{d}r}=(n+1)K\rho_c^{1/n}\,\tfrac{\mathrm{d}\theta}{\mathrm{d}r}$ , because the factor $\theta$ from differentiating pressure divides against the $\theta^{\,n}$ in the density, leaving no leftover powers of $\theta$ outside the derivative. Why: this is the algebraic heart of the reduction — the polytropic exponent is precisely the value that linearizes the pressure term in $\theta$ .

Step 4. Nondimensionalize. Introduce $r=\alpha\xi$ with $\alpha^2=(n+1)K\rho_c^{1/n-1}/(4\pi G)$ , chosen so that the constants $(n+1)$ , $K$ , $\rho_c$ , and $4\pi G$ all cancel. What remains is $\tfrac{1}{\xi^2}\tfrac{\mathrm{d}}{\mathrm{d}\xi}\!\big(\xi^2\tfrac{\mathrm{d}\theta}{\mathrm{d}\xi}\big)=-\theta^{\,n}$ , a pure equation with no reference to any star's absolute size, mass, or composition. As a check, $n=1$ gives $\theta=\sin\xi/\xi$ ( $\xi_1=\pi$ ) and $n=0$ gives $\theta=1-\xi^2/6$ ( $\xi_1=\sqrt6$ ), both satisfying the boundary conditions. Why: the entire physical individuality of a star has been packed into $\alpha$ and $\rho_c$ , so one numerical integration serves every polytrope of index $n$ at once.

Go deeper · Derivation B · 1 minThe mass integral and the birth of the Chandrasekhar limit

Derivation BThe mass integral and the birth of the Chandrasekhar limit

The goal is to evaluate the total mass of a polytrope and expose why the index $n=3$ produces a mass fixed entirely by fundamental constants — the most consequential fact in white-dwarf physics. We turn the mass integral into a surface term using the Lane–Emden equation, express radius and mass as powers of the central density, and eliminate that density to reach a mass–radius relation that breaks down at $n=3$ .

Step 1. Write $M=\int_0^R 4\pi r^2\rho\,\mathrm{d}r=4\pi\alpha^3\rho_c\int_0^{\xi_1}\xi^2\theta^{\,n}\,\mathrm{d}\xi$ using $r=\alpha\xi$ and $\rho=\rho_c\theta^{\,n}$ . The Lane–Emden equation lets you replace $\xi^2\theta^{\,n}=-\tfrac{\mathrm{d}}{\mathrm{d}\xi}(\xi^2\theta')$ , an exact derivative, so the integral evaluates at the surface: $M=4\pi\alpha^3\rho_c\,(-\xi^2\theta')_{\xi_1}$ . Why: the mass never requires the interior profile, only a single tabulated boundary number.

Step 2. Expose the density scaling. From $\alpha\propto\rho_c^{(1-n)/(2n)}$ follows $R=\alpha\xi_1\propto\rho_c^{(1-n)/(2n)}$ and $M\propto\alpha^3\rho_c\propto\rho_c^{(3-n)/(2n)}$ ; eliminating $\rho_c$ gives the mass–radius relation $R\propto M^{(1-n)/(3-n)}$ . Why: this one exponent governs every degenerate-star mass–radius curve, and its blow-up as $n$ passes the critical values is the whole story.

Step 3. Read off the non-relativistic case. For $P\propto\rho^{5/3}$ , i.e. $n=1.5$ , $R\propto M^{(1-1.5)/(3-1.5)}=M^{-1/3}$ , equivalently $MR^3=\text{const}$ : more massive white dwarfs are physically smaller, a direct fingerprint of degeneracy pressure. Why: this shrinking-with-mass behaviour is the observational signature that a star is held up by degenerate electrons, not thermal gas.

Step 4. Take the relativistic limit. As electrons become ultra-relativistic, $P\propto\rho^{4/3}$ , i.e. $n=3$ , and the mass exponent $(3-n)/(2n)$ vanishes: $M\propto\rho_c^{0}$ — the mass no longer depends on central density. The radius exponent $1/(3-n)$ diverges, so $R$ can take any value while $M$ is pinned to a single number set by $K$ (built from $h$ , $c$ , and the electron mass) through $M\propto K^{3/2}$ . Why: this is the Chandrasekhar mass — a maximum white-dwarf mass emerging not from a balance of forces but from the pure kinematics of relativistic degeneracy, with value $\approx1.44\,M_\odot$ (worked below).

Go deeper · Derivation C · 1 minCentral pressure, temperature, and the true hydrostatic coefficient

Derivation CCentral pressure, temperature, and the true hydrostatic coefficient

The goal is to turn the polytrope's scale factors into a sharp statement about the stellar center — the pressure and temperature that drive nuclear burning — and to show that the honest numerical coefficient far exceeds the order-of-magnitude guess $GM^2/R^4$ . We evaluate the central pressure, substitute the scale factors so that the central density cancels, and convert pressure and density into a central temperature via the ideal-gas law.

Step 1. Since $\theta=1$ at the center, $P_c=K\rho_c^{1+1/n}$ . Substituting $K$ , $\alpha$ , $M$ , and $R$ , every appearance of $\rho_c$ cancels, leaving $P_c=W_n\,GM^2/R^4$ with the pure number $W_n=1/\{4\pi(n+1)[\theta'(\xi_1)]^2\}$ . Why: the cancellation of $\rho_c$ signals that central pressure is fixed by global mass and radius alone — interior detail enters only through the slope $\theta'(\xi_1)$ .

Step 2. Evaluate the coefficient for $n=3$ : the numerical solution gives $\xi_1=6.897$ , $-\xi_1^2\theta'(\xi_1)=2.018$ , hence $\theta'(\xi_1)=-0.04243$ and $W_3=1/\{4\pi(4)(0.04243)^2\}=11.05$ — about eleven times the naive dimensional estimate. Why: forgetting this factor of eleven underestimates the central pressure by an order of magnitude, badly misjudging whether a star can ignite hydrogen.

Step 3. Convert to a central temperature. For a core described by the ideal gas law $P=\rho k_B T/(\mu m_H)$ , invert at the center to get $T_c=(\mu m_H/k_B)(P_c/\rho_c)$ , where $\mu$ is the mean molecular weight and $m_H$ the atomic mass unit. Pairing the polytropic $P_c$ and $\rho_c$ yields a central temperature good to a factor of about two against detailed models. Why: this is exactly the estimate that tells you whether the core is hot enough for the proton–proton chain or the CNO cycle, and it comes for free once the polytrope is solved.

Limiting cases

The index $n$ sets the entire character of the model; a handful of special values are either exactly solvable or physically privileged.

Case	Result	Meaning
$n=0$ (incompressible)	$\theta = 1 - \tfrac{1}{6}\xi^2,\ \ \xi_1=\sqrt6\approx2.449$	Uniform density $\rho=\rho_c$ ; pressure is a parabola. Crude but exact, and $\rho_c/\bar\rho=1$ .
$n=1$	$\theta = \dfrac{\sin\xi}{\xi},\ \ \xi_1=\pi$	Radius $R\propto\sqrt K$ is independent of central density and mass — a curiosity exploited in neutron-star toy models; $\rho_c/\bar\rho=\pi^2/3=3.29$ .
$n=1.5$ (adiabatic $\gamma=5/3$ )	$P\propto\rho^{5/3},\ \ \xi_1\approx3.654,\ \ \rho_c/\bar\rho=5.99$	Non-relativistic degenerate electrons and fully convective stars; delivers the white-dwarf $MR^3=\text{const}$ relation.
$n=3$ (relativistic $\gamma=4/3$ )	$P\propto\rho^{4/3},\ \ \xi_1\approx6.897,\ \ \rho_c/\bar\rho=54.2$	Eddington standard model and ultra-relativistic degeneracy; mass fixed by $K$ alone — the Chandrasekhar mass.
$n=5$	$\theta = \left(1+\tfrac{1}{3}\xi^2\right)^{-1/2},\ \ \xi_1=\infty$	Infinite radius but finite mass; the exact borderline of bounded configurations. For $n>5$ the mass diverges.
$n\to\infty$ (isothermal)	$P\propto\rho,\ \ \rho\propto r^{-2}$ far out	Isothermal self-gravitating sphere; unbounded mass, the model behind Bonnor–Ebert cloud cores.

Common trapA polytrope is not the same thing as an adiabatic gas element

The polytropic relation $P=K\rho^{1+1/n}$ is a global constraint on the whole star at once, whereas the adiabatic relation $P\propto\rho^{\gamma_{\rm ad}}$ describes how a single gas element changes under compression. Their exponents can coincide numerically — an $n=1.5$ polytrope and an adiabatic monatomic gas both have $\gamma=5/3$ — but numerical agreement is not physical equivalence. A polytrope is only self-consistent when some real microphysics actually enforces the assumed $P(\rho)$ throughout: electron degeneracy, uniform-entropy convection, or a fixed radiation-to-gas pressure ratio. Assuming a polytrope automatically satisfies energy generation and transport is the classic error; those equations decouple only in the mechanical problem and must be checked separately.

Worked examples

1 · The Sun as an n = 3 polytrope (confidence builder)

Model the Sun as an $n=3$ polytrope with $M=1.989\times10^{30}\ \mathrm{kg}$ and $R=6.96\times10^{8}\ \mathrm{m}$ . Estimate its central pressure, central density, and central temperature, and compare to a detailed solar model.

Step 1. Central pressure from equation (6) with $W_3=11.05$ . The dimensional factor is $GM^2/R^4=(6.674\times10^{-11})(1.989\times10^{30})^2/(6.96\times10^{8})^4=1.125\times10^{15}\ \mathrm{Pa}$ , so $P_c=11.05\times1.125\times10^{15}=\mathbf{1.24\times10^{16}\ Pa}$ . A detailed Standard Solar Model gives $2.3\times10^{16}\ \mathrm{Pa}$ — within a factor of two.

Step 2. Central density from the central-to-mean ratio. The mean density is $\bar\rho=3M/(4\pi R^3)=1.41\times10^{3}\ \mathrm{kg\,m^{-3}}$ , and for $n=3$ the tabulated ratio is $\rho_c/\bar\rho=54.2$ , so $\rho_c=\mathbf{7.6\times10^{4}\ kg\,m^{-3}}$ . The detailed model has $1.5\times10^{5}\ \mathrm{kg\,m^{-3}}$ — again within a factor of two, with the extreme central concentration correctly captured.

Step 3. Central temperature from the ideal-gas inversion $T_c=(\mu m_H/k_B)(P_c/\rho_c)$ with $\mu=0.62$ (ionized solar composition), $m_H=1.66\times10^{-27}\ \mathrm{kg}$ , $k_B=1.38\times10^{-23}\ \mathrm{J\,K^{-1}}$ gives $T_c=\mathbf{1.2\times10^{7}\ K}$ , against the true value $1.57\times10^{7}\ \mathrm{K}$ . A single polytrope, with no energy equations, reproduces the Sun's core conditions to within a factor of two — enough to confirm hydrogen fusion is possible.

2 · The Chandrasekhar mass, worked in full

Derive the maximum mass of a white dwarf by treating it as an $n=3$ polytrope supported by ultra-relativistic degenerate electrons, with mean molecular weight per electron $\mu_e=2$ .

Step 1. Fix the polytropic constant. Ultra-relativistic degenerate electrons give $P=K\rho^{4/3}$ with $K=\tfrac18\left(\tfrac{3}{\pi}\right)^{1/3}\dfrac{hc}{(\mu_e m_H)^{4/3}}$ . Inserting $h=6.626\times10^{-34}$ , $c=2.998\times10^{8}$ , $m_H=1.673\times10^{-27}$ (SI), $\mu_e=2$ gives $K=4.89\times10^{9}$ . This $K$ depends only on fundamental constants and composition — no free thermal parameter.

Step 2. Insert into the polytropic mass. For $n=3$ the density dependence drops out, leaving $M=4\pi\,(-\xi^2\theta')_{\xi_1}\left[\dfrac{(n+1)K}{4\pi G}\right]^{3/2}$ with $(-\xi^2\theta')_{\xi_1}=2.018$ . The bracket is $(4)(4.89\times10^{9})/[4\pi(6.674\times10^{-11})]=2.33\times10^{19}$ , whose $3/2$ power is $1.13\times10^{29}$ .

Step 3. Assemble the mass: $M=4\pi(2.018)(1.13\times10^{29})=2.85\times10^{30}\ \mathrm{kg}$ , and dividing by $M_\odot=1.989\times10^{30}\ \mathrm{kg}$ gives $M_{\rm Ch}=\mathbf{1.44\,M_\odot}$ . In closed form this is the celebrated $M_{\rm Ch}\approx5.83\,\mu_e^{-2}\,M_\odot$ ( $1.46\,M_\odot$ for $\mu_e=2$ ; the difference is rounding). The mass is set by $h$ , $c$ , $G$ , and the electron count alone — no star heavier than about $1.4\,M_\odot$ can hold itself up on electron degeneracy, the wall observed white dwarfs never cross.

Physics of the portrait

The polytrope is where nineteenth-century thermodynamics of gas spheres grew into the twentieth-century theory of stellar structure. Its equation carries two names, Lane and Emden, and its most famous solution carries a third, Chandrasekhar — but it began with a self-taught American patent examiner asking what temperature the Sun must have at its center.

Origins · 1 minJonathan Homer Lane and how this was discovered

Jonathan Homer Lane

1819–1880, American physicist and inventor · the first mathematical model of a star as a gas sphere in equilibrium

Jonathan Homer Lane was an American physicist, patent examiner, and inventor who, in an 1870 paper "On the Theoretical Temperature of the Sun," became the first to treat a star as a gaseous sphere held in equilibrium between gravity and pressure and to write down the equation now bearing his name. Working alone and against the prevailing belief that the Sun was a cooling liquid or solid, Lane assumed the solar material behaved as a perfect gas and derived its governing equation, deducing the astonishing conclusion — decades ahead of its time — that a contracting gaseous star must actually grow hotter. His analysis was systematized by Robert Emden, whose 1907 treatise Gaskugeln ("Gas Spheres") tabulated the solutions in the form still used, so the equation is universally called Lane–Emden. Arthur Eddington built his 1926 standard model on the $n=3$ polytrope, and Subrahmanyan Chandrasekhar, computing the mass of a relativistically degenerate $n=3$ sphere on a 1930 voyage from India to England, found it could not exceed about $1.4\,M_\odot$ — the limit that won him the 1983 Nobel Prize and explained why massive stars must end as neutron stars or black holes. The humble power law Lane wrote down runs like a thread from the temperature of the Sun to the death of stars.

History

1870

— Lane derives the equation for a self-gravitating gas sphere and infers a contracting gas star heats up

1878

— Ritter develops the theory of polytropic gas spheres and adiabatic stability

1907

— Emden publishes Gaskugeln, systematizing and tabulating the Lane–Emden solutions

1926

— Eddington builds the standard model on the n = 3 polytrope in The Internal Constitution of the Stars

1931

— Chandrasekhar publishes the limiting mass of a relativistically degenerate n = 3 configuration

1939

— Chandrasekhar's An Introduction to the Study of Stellar Structure fixes polytropes as textbook canon

Going further

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., §10.5 (Polytropic Models and the Lane–Emden Equation) and §16.4 (white-dwarf mass–radius relation) · Chandrasekhar, An Introduction to the Study of Stellar Structure, ch. IV · Kippenhahn, Weigert & Weiss, Stellar Structure and Evolution, 2nd ed., ch. 19 · Hansen, Kawaler & Trimble, Stellar Interiors, 2nd ed., §7.2 · Binney & Tremaine, Galactic Dynamics, 2nd ed., §4.3 (polytropic stellar systems). The tabulated eigenvalues used throughout this article are: $n=0$ gives $\xi_1=\sqrt6=2.449$ , $(-\xi^2\theta')_{\xi_1}=2\sqrt6=4.899$ ; $n=1$ gives $\xi_1=\pi$ , $(-\xi^2\theta')_{\xi_1}=\pi$ ; $n=1.5$ gives $\xi_1=3.654$ , $2.714$ , $\rho_c/\bar\rho=5.99$ ; $n=3$ gives $\xi_1=6.897$ , $2.018$ , $\rho_c/\bar\rho=54.2$ ; $n=5$ gives $\xi_1=\infty$ , $(-\xi^2\theta')\to\sqrt3$ . The central-pressure coefficient $W_3=11.05$ and the Chandrasekhar mass $M_{\rm Ch}=1.44\,M_\odot$ ( $\mu_e=2$ ) follow directly from these eigenvalues and standard constants.

Prerequisites. Hydrostatic Equilibrium in Stars · Mass Continuity Equation · Equation of State for Stellar Matter · Equations of Stellar Structure. Leads to. Radiation Pressure in Stellar Interiors · Convective Energy Transport · Adiabatic Temperature Gradient.

Ask the tutor: Ask me to walk through the Lane–Emden numerical integration for any index $n$ , or to show exactly how the $n=3$ mass loses its dependence on central density to become the Chandrasekhar limit.

What is the difference between the polytropic index n and the exponent γ?

They are two names for the same information, related by $\gamma=1+1/n$ . The exponent $\gamma$ is the power of density in the pressure law $P\propto\rho^{\gamma}$ , read straight off an equation of state — $5/3$ for a non-relativistic gas, $4/3$ for a relativistic one. The index $n$ is what appears in the Lane–Emden equation and its tabulated eigenvalues. Small $n$ means large $\gamma$ (stiff matter); large $n$ means $\gamma$ approaching $4/3$ (soft matter, extreme central concentration). Whenever a source quotes one, convert with $n=1/(\gamma-1)$ .

Does a polytrope tell you anything about a star's temperature or luminosity?

Not directly, and that is the whole point. The polytropic assumption closes the mechanical structure — density, pressure, mass, radius — without ever invoking temperature, energy generation, or opacity. Once $P(r)$ and $\rho(r)$ are solved you can afterward attach a temperature via the ideal gas law $T=(\mu m_H/k_B)P/\rho$ , as we did for the Sun's core. But the luminosity and the detailed temperature gradient require the energy equations, which the polytrope deliberately leaves untouched — a strength for quick estimates, a limitation for anything demanding a real thermal profile.

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

Because degeneracy pressure, unlike thermal pressure, does not care about temperature — it comes from the Pauli exclusion principle packing electrons into momentum states. Adding mass strengthens gravity and compresses the star, raising the density and hence the degeneracy pressure until balance returns; but that pressure grows less steeply with compression than gravity demands, so the equilibrium radius shrinks. The $n=1.5$ polytrope makes this precise: $R\propto M^{-1/3}$ , i.e. $MR^3=$ constant. This shrinking-with-mass behaviour is the observational fingerprint of a degenerate star, exactly opposite to ordinary gas spheres.

What happens physically at the Chandrasekhar mass?

As a white dwarf approaches $\sim1.4\,M_\odot$ , its central electrons turn ultra-relativistic and the equation of state softens from $P\propto\rho^{5/3}$ toward $P\propto\rho^{4/3}$ — the index climbs toward $n=3$ . At exactly $\gamma=4/3$ pressure grows with compression at precisely the rate gravity demands, so equilibrium becomes indifferent to central density and the mass is pinned to one value regardless of radius. Above that mass no static balance exists: gravity always wins and the star collapses, either igniting carbon burning (a Type Ia supernova) or proceeding to a neutron star.

Why do polytropes with n greater than 5 have no finite mass?

The mass scales as $M\propto\rho_c^{(3-n)/(2n)}$ , and the density profile stretches farther out as $n$ grows. At $n=5$ the profile $\theta=(1+\xi^2/3)^{-1/2}$ never reaches zero, so the star is formally infinite in radius, yet the density falls fast enough that the total mass stays finite. Push past $n=5$ and the tail decays too slowly: the mass integral diverges. So $n=5$ is the exact boundary between bounded and unbounded self-gravitating polytropic spheres, and physical stars require $n<5$ .

How accurate is a polytrope compared to a real stellar model?

For mechanical quantities, remarkably good. Modelling the Sun as an $n=3$ polytrope reproduces central pressure, density, and temperature to within a factor of about two, and the white-dwarf mass–radius relation and Chandrasekhar mass come out essentially exact because degenerate matter is genuinely polytropic. Where polytropes fail is anything driven by the energy equations — luminosities, temperature gradients, convective boundaries, and evolutionary timescales all need the full coupled system. Rule of thumb: trust a polytrope for structure and scaling, verify numerically for anything thermal.