Stellar Remnants: White Dwarfs, Neutron Stars & Black Holes

In many astrophysical regimes—such as completely convective stars, degenerate white dwarfs, and neutron stars—the pressure depends solely on the local density according to a power-law polytropic equation of state:

where $K$ is the polytropic constant, $\gamma = 1 + 1/n$ is the polytropic exponent, and $n$ is the polytropic index ($n = \frac{1}{\gamma - 1}$).

Derivation of the Lane-Emden Equation

Combine the hydrostatic equation $\frac{dP}{dr} = -\frac{G M(r)\rho}{r^2}$ with mass conservation $\frac{dM(r)}{dr} = 4\pi r^2 \rho$ by dividing by $\rho$, multiplying by $r^2$, and differentiating with respect to $r$:

Define dimensionless variables for density and radius anchored to the central density $\rho_c$:

Substituting into the combined hydrostatic equation yields the celebrated Lane-Emden equation of index $n$:

subject to central boundary conditions $\theta(0) = 1$ (since $\rho(0) = \rho_c$) and $\theta'(0) = 0$ (zero gravitational force at the center).

Analytical Solutions

Exact closed-form analytical solutions exist for three integer indices:

• $n = 0$ (Uniform incompressible sphere, $\rho = \text{const}$): $$\theta(\xi) = 1 - \frac{\xi^2}{6}, \quad \xi_1 = \sqrt{6} \approx 2.449$$
• $n = 1$ (Linear gas): $$\theta(\xi) = \frac{\sin\xi}{\xi}, \quad \xi_1 = \pi \approx 3.142$$
• $n = 5$ (Infinitely extended sphere): $$\theta(\xi) = \left(1 + \frac{\xi^2}{3}\right)^{-1/2}, \quad \xi_1 = \infty$$

For all other indices (notably $n = 3/2$ for non-relativistic degeneracy and $n = 3$ for relativistic degeneracy), the equation is solved numerically. The surface occurs at the first zero $\xi = \xi_1$ where $\theta(\xi_1) = 0$.

In a white dwarf, stellar matter is compressed to extreme densities ($\rho \sim 10^5 - 10^9\text{ g cm}^{-3}$). Under these conditions, atoms are completely pressure-ionized into bare nuclei surrounded by a sea of free electrons. The average separation between electrons becomes comparable to their quantum mechanical de Broglie wavelengths, and Fermi-Dirac quantum statistics completely dominates thermal pressure.

Fermi Energy & Fermi Momentum

Electrons are spin-$1/2$ fermions subject to the Pauli Exclusion Principle: each quantum phase-space cell $h^3 = (2\pi\hbar)^3$ accommodates at most 2 electrons (spin up and spin down). In the zero-temperature limit ($T \to 0$, valid since $k T \ll E_F$ even at $10^7\text{ K}$), electrons pack the Fermi sphere up to a maximum momentum $p_F$:

For fully ionized matter with mass fraction of hydrogen $X$, helium $Y$, and metals $Z$, the electron number density is $n_e = \frac{\rho}{\mu_e m_u}$, where the mean molecular weight per electron is $\mu_e = \frac{2}{1+X} \approx 2.0$ for pure helium, carbon, or oxygen.

1. Non-Relativistic Degeneracy Pressure ($p_F \ll m_e c$)

When electron velocities are non-relativistic ($v = p/m_e$), the pressure is computed from the momentum flux tensor $P = \frac{1}{3}\int_0^{p_F} p v \, n(p) dp$:

This corresponds exactly to a polytrope of index $n = 3/2$ ($\gamma = 5/3$). Crucially, $P_e$ is completely independent of temperature $T$! A white dwarf does not shrink as it radiates heat and cools.

2. Ultra-Relativistic Degeneracy Pressure ($p_F \gg m_e c$)

At extreme densities ($\rho \gtrsim 10^6\text{ g cm}^{-3}$), the Fermi momentum approaches and exceeds $m_e c$. Electrons move at speeds near the speed of light ($v \approx c$):

The polytropic index shifts to $n = 3$ ($\gamma = 4/3$). The equation of state 'softens' dramatically from $\rho^{5/3}$ to $\rho^{4/3}$, which will destabilize the star.

In 1930, on his voyage from India to England at the age of nineteen, Subrahmanyan Chandrasekhar applied relativistic Fermi-Dirac statistics to stellar polytropes, making the revolutionary discovery that there exists an absolute upper mass limit beyond which an electron-degenerate star cannot exist.

1. The Non-Relativistic Mass-Radius Relation ($n = 3/2$)

For any polytrope of index $n$, the mass and radius scale with central density $\rho_c$ as:

For non-relativistic degenerate electrons ($n = 3/2$):

Substituting the physical constants for $\mu_e = 2$:

More massive white dwarfs are physically smaller and denser! Adding mass compresses the star to higher densities to generate the requisite degeneracy pressure.

2. Exact Derivation of the Chandrasekhar Limit ($n = 3$)

As the mass $M$ increases, the radius shrinks, driving the central density $\rho_c$ up until electrons become ultra-relativistic, transitioning the star into an $n = 3$ polytrope ($P = K_2 \rho^{4/3}$). For $n = 3$:

The mass of an $n=3$ polytrope is a single unique eigenvalue:

From the numerical integration of the Lane-Emden equation for $n=3$, the boundary values are $\xi_1 = 6.89685$ and $\omega_3 = -\xi_1^2 \theta'(\xi_1) = 2.01824$. Substituting $K_2 = \frac{(3\pi^2)^{1/3} \hbar c}{4 (\mu_e m_u)^{4/3}}$:

Inserting fundamental physical constants ($h, c, G, m_u$):

For carbon-oxygen white dwarfs with $\mu_e = 2.00$:

If a white dwarf accretes matter from a binary companion pushing its mass beyond $1.44 M_\odot$, relativistic electron degeneracy pressure is mathematically incapable of supporting it. The star either undergoes runaway thermonuclear detonation (Type Ia Supernova) or collapses into a neutron star.

When an iron stellar core collapses beyond the Chandrasekhar limit, electrons and protons are crushed together via inverse beta decay: $p + e^- \to n + \nu_e$. The resulting remnant is a neutron star—a gigantic macroscopic atomic nucleus of mass $M \sim 1.4 - 2.0 M_\odot$ compressed into a radius of merely $R \sim 10-12\text{ km}$ with core densities $\rho \sim 3-8 \times 10^{14}\text{ g cm}^{-3}$.

The Need for General Relativity: Compactness Parameter

The gravitational compactness parameter $\Xi$ measures the severity of spacetime curvature:

For the Earth, $\Xi \sim 10^{-9}$; for the Sun, $\Xi \sim 4 \times 10^{-6}$. For a $1.4 M_\odot, 11\text{ km}$ neutron star:

Spacetime curvature is massive ($R \approx 2.6 R_s$). Newtonian gravity fails completely, requiring Einstein's General Relativity.

Derivation of the Tolman-Oppenheimer-Volkoff (TOV) Equation

Starting with the spherically symmetric static Schwarzschild metric $ds^2 = -e^{2\Phi(r)} c^2 dt^2 + e^{2\Lambda(r)} dr^2 + r^2 d\Omega^2$ and solving Einstein's field equations $G^\mu_\nu = \frac{8\pi G}{c^4} T^\mu_\nu$ with a perfect fluid stress-energy tensor $T^\mu_\nu = \text{diag}(-\rho c^2, P, P, P)$ yields the Tolman-Oppenheimer-Volkoff (TOV) equation of relativistic hydrostatic equilibrium (1939):

Notice the three relativistic correction bracket terms:

1. $\left[1 + \frac{P}{\rho c^2}\right]$: In General Relativity, pressure itself possesses effective mass-energy ($E = m c^2$), contributing directly to gravitational attraction! High pressure intensifies gravity.
2. $\left[1 + \frac{4\pi r^3 P}{M(r) c^2}\right]$: Accounts for the enclosed pressure work contributing to the total active gravitational mass.
3. $\left[1 - \frac{2 G M(r)}{r c^2}\right]^{-1}$: The metric spatial curvature correction ($e^{2\Lambda}$), amplifying the local gravitational gradient as $r$ approaches the Schwarzschild radius.

The Tolman-Oppenheimer-Volkoff Mass Limit

Because pressure self-gravitates in General Relativity, increasing the central density eventually increases inward gravitational pull faster than outward degenerate pressure. Even with the strong repulsive nuclear force (meson exchange between neutrons), there exists a rigorous upper limit for neutron stars:

Any core remnant exceeding $M_{\text{TOV}}$ (such as the primary object in binary merger GW170817) cannot stabilize and collapses into a black hole.

Discovered in 1967 by Jocelyn Bell Burnell and Antony Hewish as periodic radio beeps ($P = 1.337\text{ s}$), pulsars were identified by Thomas Gold and Franco Pacini as rapidly rotating, highly magnetized neutron stars.

Conservation of Magnetic Flux & Magnetic Field Strength

During core collapse, the progenitor's magnetic field is trapped within the highly conducting collapsing plasma. By Alfvén's theorem of flux freezing, total magnetic flux is conserved: $\Phi_B = B R^2 = \text{const}$. If a progenitor of radius $R_0 = 10^6\text{ km}$ and field $B_0 = 100\text{ G}$ collapses to $R = 10\text{ km}$:

Neutron stars possess the strongest magnetic fields in the cosmos (up to $10^{14}-10^{15}\text{ G}$ in magnetars).

Magnetic Dipole Radiation & Spin-Down Power

A neutron star spinning with angular frequency $\Omega = 2\pi/P$ with magnetic dipole moment $\vec{m}$ ($m = B_p R^3$) inclined at magnetic obliquity angle $\alpha$ relative to its rotation axis emits magnetic dipole radiation into space at a rate governed by electrodynamics:

This radiant power is extracted directly from the star's rotational kinetic energy $E_{\text{rot}} = \frac{1}{2} I \Omega^2$ (where $I \approx 0.4 M R^2 \sim 10^{45}\text{ g cm}^2$ is the moment of inertia):

Equating $\dot{E}_{\text{rot}} = \dot{E}_{\text{dipole}}$:

Characteristic Age & Surface Magnetic Field Determination

From the observable period $P$ and its measured spin-down rate $\dot{P} = dP/dt > 0$:

For the Crab Pulsar ($P = 0.033\text{ s}, \dot{P} = 4.2 \times 10^{-13}\text{ s/s}$), $\tau_c \approx 1240\text{ years}$ (closely matching the historical supernova of 1054 CE) and $B_p \approx 3.8 \times 10^{12}\text{ G}$.

When the mass of a collapsing stellar core exceeds the TOV limit ($M > 2.3 M_\odot$), no known physical force can oppose gravity. The core undergoes complete gravitational collapse to a spacetime singularity, cloaked behind an event horizon.

The Schwarzschild Spacetime Metric

Karl Schwarzschild (1916) found the exact static, spherically symmetric vacuum solution to Einstein's Field Equations:

The metric coefficient $g_{00} = -(1 - R_s/r)$ vanishes and $g_{rr} = (1 - R_s/r)^{-1}$ diverges at the Schwarzschild radius:

For a stellar black hole of $10 M_\odot$, $R_s \approx 29.5\text{ km}$.

Gravitational Redshift & Time Dilation

A clock at rest at coordinate radius $r$ ticks at proper time $d\tau = \sqrt{-g_{00}} dt = \sqrt{1 - R_s/r} dt$. Light emitted at radius $r_e$ with frequency $\nu_e$ is received by an observer at infinity with frequency $\nu_\infty$:

As an infalling source reaches the event horizon ($r_e \to R_s$), $z \to \infty$: emitted photons red-shift to infinite wavelength, and the object appears frozen and dimmed to blackness.

Geodesics, Photon Sphere & The ISCO

Particle orbits in the equatorial plane ($\theta = \pi/2$) obey the effective potential equation $\frac{1}{2}\left(\frac{dr}{d\tau}\right)^2 + V_{\text{eff}}(r) = \frac{\mathcal{E}^2 - c^2}{2}$:

Notice the general relativistic term $-\frac{G M L^2}{c^2 r^3}$, which pulls particles inward at small radii, absent in Newtonian physics.

• The Photon Sphere ($r_{\text{ph}} = 1.5 R_s = 3GM/c^2$): The radius at which photons can execute unstable circular orbits. Light emitted tangentially here circles the black hole.
• Innermost Stable Circular Orbit (ISCO, $r_{\text{ISCO}} = 3 R_s = 6GM/c^2$): Inside $6GM/c^2$, circular orbits of massive particles become dynamically unstable ($d^2 V_{\text{eff}}/dr^2 < 0$). Matter in an accretion disk spirals into the horizon.

Accretion Efficiency & X-Ray Binaries

The binding energy of a particle at the ISCO of a Schwarzschild black hole is:

Accretion onto a non-rotating black hole converts $\eta \approx 5.7\%$ of rest-mass energy into radiation (for a maximally rotating Kerr black hole, $\eta \approx 42.3\%$!), dwarfing nuclear fusion ($\eta \approx 0.7\%$). This immense efficiency powers Galactic X-ray binaries (Cygnus X-1) and quasars.