Past exam of the mathematics course of the University of Cambridge / 2015 / iii / Paper 58 (Structure and Evolution of Stars)

Write $\mathcal{R}=k/(\mu H)$, where $H$ is the atomic mass unit used in the mean molecular weight convention, and let $\kappa$ be the Rosseland mean opacity per unit mass. For a static spherical star with negligible thermal storage, the stellar structure equations areHere $ϵ$ is the net stellar energy-generation rate per unit mass; $a$ is the radiation energy-density constant and $c$ the speed of light. The last differential equation follows from radiative diffusion, $F=-(c/3\kappa \rho )d(a{T}^4)/dr$, with $L=4\pi r^{2}F$. There is no gravothermal storage term in the specified stellar energy balance equation. Neglecting radiation pressure in force balance does not remove radiation as an energy carrier. Regular central conditions are $m(0)=L(0)=0$; take negligible external pressure, $P(R)=0$, for the analytic hydrostatic equilibrium model below.

Set $x=r/R$. Integrating the enclosed mass directly givesThusThe stellar hydrostatic equation then givesFactorizing makes the surface behaviour transparent:The ideal gas equation of state gives the temperature profileNo small-$x$ approximation was needed. In particular, close to the centre,These exact hydrostatic formulas reveal a limitation of the prescribed linear-density stellar model. Its mass density has a central cusp, and $T^{\prime }(0^{+})=T_c/R>0$. A positive opacity with positive outward luminosity would instead require $T^{\prime }<0$ in radiative diffusion. Moreover, regular positive nuclear heating gives $L=O(r^3)$ and therefore $T^{\prime }=O(r)$ at the centre. Thus the imposed profile is a hydrostatic toy model, not a complete positive-heating radiative-equilibrium solution of all the stellar structure equations. The requested algebraic profiles and nuclear-luminosity integral can still be computed consistently as properties of that toy model.

Evaluate radiation pressure on the gas-supported hydrostatic equilibrium solution. At the centre,The factors $R^{-4}$ cancel between $T_c^4$ and $P_c$, so at fixed mean molecular weight the ratio depends only on mass. Equivalently it is $(M/M_{\mathrm{eq}}{)}^2$, withFor $H\simeq 1.6605\times 10^{-27}\mathrm{kg}$ this is approximately $14.6{\mu }^{-2}{M}_{\odot }$, or $40.5M_{\odot }$ if $\mu =0.6$. If $H$ denotes hydrogen-atom mass instead, use the corresponding mean molecular weight convention consistently; the coefficient is about $14.4{\mu }^{-2}{M}_{\odot }$. Equality marks the breakdown of the gas-pressure approximation. It is a diagnostic based on the gas-only model, not an exact mixed-pressure stellar solution.

Since $ϵ$ is per unit mass, each shell contributes $dL=ϵdm=4\pi r^2\rho ϵdr$. Both factors of mass density must be retained. Substituting the linear-density stellar model and its temperature profile givesWriting $I_{\nu }(z)$ for this incomplete integral, with $I_{\nu }=I_{\nu }(1)$, the luminosity isThis supplies both the local $L_r$ and the total luminosity associated with the printed integral extending to one. For fixed composition and fixed $\nu$, the total scales as $M^{2+\nu }{R}^{-3-\nu }$. The physical nuclear-heating exponents are positive; mathematically the surface integral converges for $\nu >-3$. As a simple normalization check, $I_0=1/30$, hence $L(R)=6{ϵ}_0{M}^2/(5\pi R^3)$ for temperature-independent heating. This is the generated luminosity; the radiative-transport incompatibility noted above is not removed by integrating the heating.

For a uniform-density stellar model, $\rho =3M/(4\pi R^3)$ and $m(r)=4\pi \rho r^3/3$. With zero surface pressure, integrate the stellar hydrostatic equation:Using the ideal gas equation of state, $T(r)=P(r)/(\mathcal{R}\rho )$. ThereforeEquivalently $P_c=3G{M}^2/(8\pi R^4)$ and $T_c=GM/(2\mathcal{R}R)$. The constant mass density is imposed as a structural approximation; a perfect gas can realize it hydrostatically by allowing its temperature and entropy to vary with radius.

Assemble the star from concentric shells. A shell of mass $dm$ at radius $r$ has interaction energy $-Gm(r)dm/r$ with the already assembled interior. This counts each gravitational pair once. Hence the Newtonian gravitational potential energy isUsing $M=4\pi \rho R^3/3$ givesThe negative sign expresses gravitational binding. A factor of one half should not be inserted again: the shell assembly already avoids double counting.

For a spherical hydrostatic star with negligible surface pressure, the stellar virial theorem isFor a monatomic nonrelativistic perfect gas, $U=\frac{3}{2}\int PdV$, so $2U+\Omega =0$ and the total energy is $E=U+\Omega =\Omega /2$. Half the released binding energy raises the internal energy; only the other half can be radiated. With fixed mass, no nuclear supply and no external work, conservation of energy givesFor the uniform-density stellar model, this becomes the contraction-luminosity evolution equationIt relates luminosity to the contraction rate; it does not by itself prescribe $L(t)$. Given $L(t)$,For constant luminosity, $R(t)=R_0/[1+t/t_0]$ with $t_0=3G{M}^2/(10R_{0}L)$. This describes quasi-static Kelvin-Helmholtz contraction, not dynamical free fall. If the gas has a different constant specific-heat ratio ${\gamma }_g$, the appropriate stellar virial theorem is $3({\gamma }_g-1)U+\Omega =0$, so the radiated fraction of binding release is $(3{\gamma }_g-4)/[3({\gamma }_g-1)]$ rather than universally one half.

The conventional Kelvin-Helmholtz cooling time is of order $G{M}^2/(RL)$. For the monatomic uniform-density stellar model, its precise accessible total-energy reservoir is $|E|=3G{M}^2/(10R)$, soUsing $M_{\odot }\simeq 1.9885\times 10^{30}\mathrm{kg}$, $R_{\odot }\simeq 6.957\times 10^8\mathrm{m}$ and $L_{\odot }\simeq 3.828\times 10^{26}\mathrm{W}$, the usual order-of-magnitude normalization without the structural factor is $t_{\mathrm{KH}}\simeq 3.1\times 10^7\mathrm{yr}$. The full gravitational binding reservoir $|\Omega |/L$ is about $1.9\times 10^7\mathrm{yr}$, but that overcounts the radiatable energy in virial equilibrium.

The Sun and the Solar system are approximately $4.6\times 10^9\mathrm{yr}$ old, much older than any of these contraction estimates. Gravitational contraction cannot sustain the Sun's long-lived present power. The dominant long-term source is stellar nuclear fusion. Contraction can still supply transient luminosity, especially in a pre-main-sequence star. The timescale comparison rules out a contraction-only explanation over the observed age; it does not assert that all gravitational energy release is absent.

For a stellar polytrope of index $n>0$, write $P=K{\rho }^{1+1/n}$ and combine hydrostatic equilibrium with mass conservation to obtainIntroduce Lane-Emden variables for a stellar polytrope, $\rho ={\rho }_c{\theta }^n$, $P=P_c{\theta }^{n+1}$ and $r=\alpha \xi$, whereSince ${\rho }^{-1}dP/dr=(n+1)P_c{\theta }^{\prime }/{\rho }_c\alpha$, the mechanical equation reduces to the Lane-Emden equationA constant-density interior corresponds to the formal polytrope of index zero, with $\rho ={\rho }_c{\theta }^0={\rho }_c$ where $\theta >0$. At $n=0$, regularity first gives ${\xi }^2{\theta }^{\prime }=-{\xi }^3/3$, and a second integration givesHere $P=P_c\theta =P_c(1-r^2/R^2)$ and ${\alpha }^2=P_c/(4\pi G{\rho }_c^2)$, reproducing $P_c=2\pi G{\rho }_c^2{R}^2/3$. The mass density jumps from its constant interior value to zero at the surface, while pressure vanishes continuously. The relation $P=K{\rho }^{1+1/n}$ is singular at $n=0$; the regular dimensionless pressure/mass density formulation defines this incompressible structural limit. It does not mean that the gas's perturbative stellar adiabatic exponent is infinite: the hydrostatic ideal-gas toy model and its adiabatic response are distinct choices.

Use the printed stellar gas-pressure fraction, $\beta =P_g/P$, with $P=P_g+P_{\mathrm{rad}}$, and keep composition fixed. For a monatomic perfect gas plus equilibrium blackbody radiation, the specific heats of a monatomic gas-radiation mixture follow from its specific internal energy and specific enthalpy:At fixed total pressure, differentiating the equation of state givesIn particular, $\beta$ is not held constant during that derivative. The specific heat capacity at constant pressure is $c_P=(\partial h/\partial T{)}_P$, soSince $a{T}^3/(3\rho )=\mathcal{R}(1-\beta )/\beta$, this simplifies toThe pure-gas limit is $5\mathcal{R}/2$. For a gas with unspecified molecular degrees of freedom, replace $3\mathcal{R}/2$ in $u$ by its gas heat capacity $c_{V,g}$: then $c_P=c_{V,g}+\mathcal{R}+4\mathcal{R}(1-\beta )(4+\beta )/{\beta }^2$. The boxed formula uses the conventional monatomic stellar-gas interpretation.

The stellar adiabatic exponents are fixed-composition, constant-specific entropy derivatives:The chain rule gives ${\Gamma }_1({\Gamma }_2-1)/{\Gamma }_2={\Gamma }_3-1$. Define the pressure derivatives ${\chi }_{\rho }=(\partial \log P/\partial \log \rho {)}_T$ and ${\chi }_T=(\partial \log P/\partial \log T{)}_{\rho }$. The first law of thermodynamics, with $du=Tds+Pd\rho /{\rho }^2$, givesCombining this with $d\log P={\chi }_{\rho }d\log \rho +{\chi }_{T}d\log T$ on an adiabat givesThus the specific-heat ratio is not generally equal to the three stellar adiabatic exponents.

For the mixture, ${\chi }_{\rho }=\beta$ and ${\chi }_T=4-3\beta$. The specific heat capacity at constant volume, obtained by differentiating $u$ at fixed mass density, isTherefore the adiabatic exponents of a monatomic gas-radiation mixture and its specific-heat ratio areIf a relation involving only the exponents is wanted, eliminate $\beta$ from ${\Gamma }_1=\beta +(4-3\beta )({\Gamma }_3-1)$:For a pure monatomic perfect gas, $\beta =1$ and $\gamma ={\Gamma }_1={\Gamma }_2={\Gamma }_3=5/3$. For any calorically perfect gas with constant heat capacities, the same equality holds with its own ${\gamma }_g$. In a genuine gas-radiation mixture, $0<\beta <1$, $\gamma ={\Gamma }_1/\beta$ while the exponents are given separately above.

In the radiation limit, ${\Gamma }_1={\Gamma }_2={\Gamma }_3=4/3$ and the adiabatic temperature gradient is $1/4$. This follows independently from photon entropy: a comoving volume $V$ has $S\propto V{T}^3$, so an adiabat obeys $T\propto V^{-1/3}$ and $P\propto V^{-4/3}$. However, $c_P/c_V$ is singular in the pure-radiation limit, not $4/3$. This is the pure-radiation constant-pressure heat-capacity singularity. The equation $P=a{T}^4/3$ fixes temperature whenever pressure is fixed, so an ordinary constant-pressure temperature derivative is not available; along the mixture limit $\gamma \to \infty$ and $\beta \gamma \to 4/3$. For photons alone, mass-specific quantities additionally require a material mass label. The often quoted radiation index $4/3$ is its pressure-density adiabatic exponent, not a finite constant-pressure/constant-volume heat-capacity ratio.

Let $g=Gm/r^2>0$, $H_P=-dr/d\log P>0$, $\nabla =d\log T/d\log P$, and ${\nabla }_{\mu }=d\log \mu /d\log P$. Consider a small radial fluid displacement $\xi$ whose sound-crossing time is short enough to keep its pressure equal to its surroundings. Its composition is frozen and its heat exchange negligible. This is the local stellar convective stability test.

For gas-pressure-dominated ideal gas matter, $\rho \propto \mu P/T$. The ambient gradient is thereforeThe displaced element conserves specific entropy and mean molecular weight, giving $(d\log \rho /dr{)}_{\mathrm{parcel}}=-(1-{\nabla }_{\mathrm{ad}})/H_P$, where ${\nabla }_{\mathrm{ad}}=({\Gamma }_2-1)/{\Gamma }_2$. The parcel-minus-environment mass density difference isIts buoyancy acceleration is minus $g$ times this difference. ThusPositive stellar buoyancy frequency squared gives a restoring force. The gas-pressure-dominated Ledoux criterion is thereforeEquality is marginal in this ideal adiabatic test. For monatomic gas ${\nabla }_{\mathrm{ad}}=2/5$. An inward increase of mean molecular weight has ${\nabla }_{\mu }>0$ and stabilizes the layer; uniform composition recovers the Schwarzschild criterion. For a radiative layer one substitutes the stellar radiative temperature gradient, ${\nabla }_{\mathrm{rad}}=3\kappa LP/(16\pi acGm{T}^4)$.

More generally define $\delta =-(\partial \log \rho /\partial \log T{)}_{P,\mu }$ and $\phi =(\partial \log \rho /\partial \log \mu {)}_{P,T}$. The same displacement argument gives $N^2=(g/H_P)[\delta ({\nabla }_{\mathrm{ad}}-\nabla )+\phi {\nabla }_{\mu }]$ and stability for $\nabla <{\nabla }_{\mathrm{ad}}+(\phi /\delta ){\nabla }_{\mu }$. Both coefficients equal one in the gas-dominated ideal-gas limit used above. This is local dynamical stability against convection; it is distinct from global radial collapse and from instabilities requiring heat or composition diffusion.

Take a static, nonrotating, nonmagnetic spherical equilibrium with $\nabla P=-\rho \nabla \Phi$ and ${\nabla }^2\Phi =4\pi G\rho$. Use a fluid displacement $\mathbf{\xi }(\mathbf{r})e^{-i\omega t}$. Write $p_1,{\rho }_1,{\varphi }_1$ for Eulerian fluid perturbations. The Lagrangian pressure perturbation and corresponding mass density change obey $\Delta P=p_1+\mathbf{\xi }\cdot \nabla P$, $\Delta \rho ={\rho }_1+\mathbf{\xi }\cdot \nabla \rho$. Linearizing mass conservation, the Euler equations for an inviscid fluid and Poisson equation for Newtonian gravity givesThere is no equilibrium acceleration to multiply a perturbed mass density. The adiabatic process condition, with composition carried by the parcel, closes the system:These are the complete linear adiabatic stellar oscillation equations, including the perturbation of self-gravity. Neglect of heat exchange is appropriate when oscillation periods are short compared with relevant thermal relaxation times; it does not determine nonadiabatic excitation or damping.

For the pressure and buoyancy modes, separate angular dependence using spherical harmonics:The horizontal amplitude ${\xi }_h$ has dimensions of length. Angular differentiation gives $\nabla \cdot \mathbf{\xi }=[r^{-2}(r^2{\xi }_r{)}^{\prime }-\ell (\ell +1){\xi }_h/r]Y_{\ell }^m$, and tangential momentum gives ${\omega }^{2}r{\xi }_h=p_{\ell }/\rho +{\varphi }_{\ell }$. Define the adiabatic sound speed, stellar buoyancy frequency and Lamb frequency byThe adiabatic mass density relation becomes ${\rho }_{\ell }/\rho =p_{\ell }/(\rho c_s^2)+(N^2/g){\xi }_r$. Substitution produces a radial form of the full oscillation equations for nonzero $\omega$:The apparent $N^2/g$ factor is evaluated through its defining gradient at the centre rather than by dividing two zeros.

Regularity at the centre excludes singular solutions. At a free surface the Lagrangian pressure perturbation vanishes, $\Delta P=0$; outside the star the gravitational perturbation decays as $r^{-\ell -1}$. For a model with mass density tending to zero at its surface, continuity of ${\varphi }_{\ell }$ and its radial derivative gives ${\varphi }_{\ell }^{\prime }(R)=-(\ell +1){\varphi }_{\ell }(R)/R$. If the equilibrium mass density jumps to vacuum, include the displaced-surface mass sheet: the outward-minus-inward derivative jump is $4\pi G\rho (R){\xi }_r(R)$, so the interior condition is ${\varphi }_{\ell }^{\prime }(R)=-(\ell +1){\varphi }_{\ell }(R)/R-4\pi G\rho (R){\xi }_r(R)$. An atmospheric boundary condition can replace the ideal free surface.

These conditions make ${\omega }^2$ an eigenvalue, not an arbitrary local sound frequency. With conservative boundary conditions the adiabatic operator is self-adjoint, giving real ${\omega }^2$; negative values describe instability. The frequencies depend on $P(r),\rho (r),{\Gamma }_1(r)$, self-gravity, stratification and boundaries, as well as angular degree $\ell$ and radial order. In a spherical nonrotating star they are degenerate in $m$. For homologous equilibrium structures, their scale isthus the typical oscillation time measures inverse square root of mean mass density, while individual frequencies probe the interior adiabatic sound speed and stellar buoyancy frequency profiles.

For radial modes, write ${\xi }_r=r\eta$. Eliminating the pressure and gravitational perturbations gives the radial stellar pulsation equationThis is a Sturm-Liouville problem. Multiplication by $\eta$ and integration, with vanishing boundary terms, gives its Rayleigh quotientFor constant ${\Gamma }_1>4/3$, $P^{\prime }<0$ makes both numerator contributions nonnegative. At ${\Gamma }_1=4/3$, a homologous displacement is neutral; for constant ${\Gamma }_1<4/3$ the same trial displacement makes the quotient negative. For the uniform-density stellar model, $\eta =$ constant is an exact mode and ${\omega }^2=(3{\Gamma }_1-4)GM/R^3$. With varying ${\Gamma }_1$, the integral criterion, rather than a universal pointwise threshold, controls radial stability.

The Cowling approximation neglects ${\varphi }_{\ell }$ while retaining the equilibrium gravitational field. It is useful for short-wavelength modes, but is not needed for the full derivation above. In a locally slowly varying region, take both remaining amplitudes proportional to $e^{i\int k_{r}dr}$ and retain the leading derivative terms. ThenEliminating either amplitude gives the acoustic-gravity propagation relationPositive $k_r^2$ is oscillatory propagation; negative $k_r^2$ means an evanescent wave. The high-frequency branch, ${\omega }^2>S_{\ell }^2,N^2$, describes stellar acoustic modes, restored chiefly by compressibility and pressure. At frequencies well above $N$, this gives ${\omega }^2\simeq c_s^2[k_r^2+\ell (\ell +1)/r^2]$. The low-frequency propagating branch in stable stratification, ${\omega }^2<N^2,S_{\ell }^2$, instead describes stellar gravity modes, restored by buoyancy. At $\ell =0$ there is no such nonradial gravity-wave cavity.

A stellar acoustic mode is trapped between an inner turning point near $\omega =S_{\ell }$ and an outer reflecting region. Low-$\ell$ modes penetrate deeply; radial modes reach the centre. Higher-degree modes turn farther out. Standing waves require the WKB quantization conditionwhere the phase $\alpha$ depends on the central or turning-point behaviour and surface reflection. For high radial order and small degree, the leading acoustic travel-time result is the large frequency separationThe $\ell /2$ term is the leading central angular phase shift; smaller frequency separations depend on detailed interior gradients. The near-surface mass density stratification sets the acoustic cutoff frequency. In a plane-parallel isothermal atmosphere with mass density scale height $H_{\rho }$, ${\omega }_{\mathrm{ac}}\simeq c_s/(2H_{\rho })$. Modes below this cutoff can reflect and form a cavity; waves sufficiently above it escape and need an outgoing-wave boundary condition. Nonadiabatic damping, driving and rotation alter real-star mode properties, but the adiabatic frequency problem isolates their dependence on the equilibrium structure.