Past exam of the mathematics course of the University of Cambridge / 2003 / iii / Paper 63 / 4

A stellar polytrope replaces a detailed equation of state and thermal structure by $P=K{\rho }^{1+1/n}$ for $n>0$, with spatially constant $K$ and polytropic index $n$. Together with spherical hydrostatic equilibrium and mass conservation, it gives a tractable family of stellar density profiles. Introduce the Lane-Emden variables for a stellar polytropeCombining $dP/dr=-Gm\rho /r^2$ with $dm/dr=4\pi r^2\rho$ eliminates $m$ to give $r^{-2}d[r^2{\rho }^{-1}(dP/dr)]/dr=-4\pi G\rho$. Substitution of the pressure law and the dimensionless variables gives the Lane-Emden equationA regular finite-radius model for $0\le n<5$ ends at the first zero ${\xi }_1$. Its radius and mass areThese formulas supply central-to-mean density ratios and stellar homology scalings without solving a full evolutionary model. Eliminating ${\rho }_c$, for $0<n<5$ and $n\ne 3$, gives the polytropic mass-radius relation

An efficient-convection monatomic gas has $n=3/2$. The same index describes a cold, nonrelativistic electron degeneracy pressure law; at fixed composition it gives $R\propto M^{-1/3}$ for a white dwarf. Ultrarelativistic electron degeneracy pressure has $n=3$, for which the mass is independent of ${\rho }_c$ and scales as $(K/G{)}^{3/2}$. This is the origin of the Chandrasekhar limit. A radiation-dominated, approximately constant-entropy interior also approaches $P\propto {\rho }^{4/3}$ and $n=3$. The polytrope of index zero is an incompressible comparison model, useful for illustrating the role of central concentration.

The gravitational energy of a stellar polytrope, $\Omega =-3G{M}^2/[(5-n)R]$, and the stellar virial theorem make polytropes useful for contraction-energy estimates. They also provide structural models for perturbation and stability calculations. A polytrope is a structural approximation, not a complete stellar evolution calculation. One must still specify stellar nuclear fusion, opacity, energy transport, composition and boundary conditions to determine the luminosity and history. Moreover, the structural exponent $1+1/n$ need not equal the adiabatic exponent relevant to a rapid perturbation. Real stars can have different effective indices in their cores and envelopes, as well as ionization zones or varying electron degeneracy pressure.

The first argument is the energy budget. The gravitational energy reservoir is of order $G{M}^2/R$, so the Kelvin-Helmholtz cooling time is $t_{\mathrm{KH}}\sim G{M}^2/(RL)$. For the Sun, $M\simeq 2\times 10^{30}\mathrm{kg}$, $R\simeq 7\times 10^8\mathrm{m}$ and $L\simeq 3.8\times 10^{26}\mathrm{W}$ give $t_{\mathrm{KH}}\sim 3\times 10^7$ years. This is much shorter than the approximately $4.6\times 10^9$-year solar-system age. Chemical binding energies are only of order electronvolts per atom and provide an even smaller reservoir. Thus neither ordinary chemical energy nor continuing gravitational contraction can supply the long-term solar luminosity.

Stellar nuclear fusion provides an adequate reservoir. Conversion of four hydrogen atoms into one helium-4 atom releases about $26.7\mathrm{MeV}$, corresponding to a mass defect of about $0.007$ of the initial mass; some energy leaves as neutrinos. If a fraction $f$ of stellar mass with hydrogen fraction $X$ is available for core burning, the energy scale is $0.007fXM{c}^2$. Taking $f\sim 0.1$ and $X\sim 0.7$ supplies roughly $9\times 10^{43}\mathrm{J}$, enough for several billion years at the solar luminosity, consistent with a stellar nuclear timescale of order $10^{10}$ years. The core temperature and density inferred from stellar structure permit the proton–proton chain; hotter stellar cores favor the CNO cycle.

The most direct evidence is the detection of solar neutrinos. They are created in the weak reactions of the proton–proton chain and can escape from the core, unlike the photons whose transport through the interior is slow. Measurements accounting for neutrino oscillation recover a total neutrino flux consistent with nuclear-powered solar models. Their spectra and reaction-dependent fluxes test the core reactions rather than merely the surface effective temperature.

There is also independent evidence from stellar evolution: the main sequence is a long-lived core-hydrogen-burning phase, more massive luminous stars consume fuel more quickly, and the main-sequence turnoff of a star cluster records its age. Subsequent evolution and the production of heavier elements fit further nuclear reactions. Agreement among these observations, nuclear energy release and stellar-model lifetimes is much stronger than the energy-budget argument alone. The combination of a sufficient nuclear reservoir, measured core neutrinos and the observed evolutionary sequence identifies nuclear reactions as the main long-term power source of ordinary stars.

A Type Ia supernova is a thermonuclear disruption associated with a carbon-oxygen white dwarf in a binary system. Electron degeneracy pressure supports the dense star, and initially depends much less on temperature than an ideal gas pressure does. Once carbon burning heats degenerate material, the immediate pressure response is insufficient to produce the usual expansion-and-cooling thermostat. The strongly temperature-sensitive nuclear reactions can therefore undergo a thermonuclear runaway.

One channel involves accretion toward the Chandrasekhar limit, about $1.4$ solar masses for a typical electron composition. Other proposed channels involve merging white dwarfs or an initial helium-layer detonation that triggers carbon burning below that mass. Thus the Chandrasekhar limit motivates an important model, but it is not an assumption that every Type Ia supernova has exactly the same progenitor mass or ignition history. Burning carbon and oxygen into more tightly bound nuclei releases enough nuclear binding energy to overcome the stellar gravitational binding energy. An ordinary successful event ejects the star rather than leaving the compact core characteristic of a core-collapse supernova.

The ejecta contain iron-group material, including radioactive nickel-56, and intermediate-mass elements from incomplete burning. The radioactively powered Type Ia supernova light curve is powered largely by the nickel-56 to cobalt-56 to iron-56 decay chain: its energetic products heat the ejecta, and radiative diffusion delays the emerging radiation. Near maximum light, ordinary Type Ia supernovae show prominent silicon absorption and lack the hydrogen lines defining Type II events. The spectral classification and the thermonuclear mechanism provide complementary descriptions.

The similar progenitors and radioactive power make Type Ia supernovae useful standard candles, but their peak luminosities are not identical. The relation between light curve width and peak brightness, together with color corrections, allows standardization for distance measurements. Extinction, composition and explosion diversity still introduce systematic uncertainties. A Type Ia event is a white-dwarf thermonuclear explosion with a radioactively powered light curve; its corrected luminosity makes it a valuable distance indicator.

An X-ray binary contains a compact accretor supplied by a donor star. Most bright systems contain a neutron star or black hole. Gas arrives by Roche-lobe overflow, a captured stellar wind, or donor outflows, commonly forming an accretion disk. The primary energy source is the release of gravitational binding energy during accretion, with a characteristic accretion luminosityFor a neutron star, $R_{\mathrm{in}}$ can be near its material radius; for a black hole, radiation is produced by the flow outside the horizon. For $M\sim 1.4M_{\odot }$ and $R\sim 10\mathrm{km}$, $GM/(R{c}^2)\sim 0.2$, so accretion can release much more energy per unit mass than hydrogen burning. The hot inner flow, boundary layer or accretion column emits X-rays.

In a high-mass X-ray binary, a massive, relatively short-lived donor often supplies gas through a strong wind or an equatorial outflow; some systems instead undergo Roche-lobe overflow. These objects trace young stellar populations. In a low-mass X-ray binary, a lower-mass donor commonly fills its Roche lobe and feeds a disk. Such systems can persist in older populations. These mass labels refer to the donor, not the compact accretor.

A magnetized neutron star can channel accretion onto its poles, producing an accreting X-ray pulsar as the star rotates. Accretion transfers angular momentum, so changes in spin constrain the flow and magnetic coupling. A Type I X-ray burst has a different origin: accumulated hydrogen or helium undergoes unstable thermonuclear burning on a neutron star surface. It rises rapidly and then cools as the fuel is exhausted. A black hole has no material surface supporting the same surface-fuel mechanism.

Orbital motion and optical observations of the donor constrain masses, with inclination and donor-model uncertainties. A compact mass exceeding plausible stable neutron star masses supports a black hole interpretation. Rapid variability probes the small emitting region, while disk spectra and recurrent outbursts probe accretion and its instabilities. The Eddington luminosity,is obtained by balancing radiative acceleration $\kappa L/(4\pi r^{2}c)$ against $GM/r^2$; it is a useful luminosity scale under the assumed isotropic, opacity-controlled force law. Geometry and time dependence can complicate its application. X-ray binaries reveal compact stars through efficient gravitational accretion, with pulsations, bursts and orbital dynamics distinguishing their physical components.