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

Introduce the characteristic scalesThen $\rho ={\rho }_{\ast }b$, $P=P_{\ast }p$ and $T=T_{\ast }t$. The gas-only equation of state gives the important dimensionless identity $p=bt$. The PDF's CNO cycle approximation has temperature power 13; the converted TeX's power 1.3 is an OCR error.

Since $d/dr=R^{-1}d/dx$, the mass and hydrostatic pressure support equation immediately becomeThus $dq/dx=x^{2}b$ and $dp/dx=-bq/x^2$. The stellar structure equations for nuclear heating givesoFor radiative diffusion in a star, $\kappa \rho ={\kappa }_0{\rho }^2{T}^{-3}$, and thereforeThis gives the radiative gradient coefficientIn an efficiently convective region, the monatomic adiabatic temperature gradient is $2/5$. Hence $dt/dx=(2/5)(t/p)(dp/dx)=-2q/(5x^2)$. If the radiative gradient is shallower, radiation suffices; if it is steeper than the adiabatic gradient, the Schwarzschild criterion selects convection and the shallower adiabatic profile. Thus

The regular central conditions are $q(0)=l(0)=0$ with finite positive $p(0),t(0)$. In particular $q\sim b(0)x^3/3$ and $l\sim Eb(0{)}^{2}t(0{)}^{13}{x}^3/3$, so the pressure and temperature derivatives tend to zero at the centre. The normalized surface has $q(1)=l(1)=1$. In the idealized zero-pressure, zero-temperature surface used for stellar homology, $p(1)=t(1)=0$; a physical photosphere supplies small nonzero surface values that are neglected in these bulk homology relations.

The common dimensionless profiles and boundary conditions keep $D,E$ constant. With fixed ${\kappa }_0,{ϵ}_0$, $D$ constant yields $L{\mu }^{-7}{M}^{-5}=\mathrm{constant}$. Substitution into the nuclear coefficient then gives $R^{16}\propto {\mu }^6{M}^{10}$. Consequently the CNO homology with density-dependent inverse-cubic opacity isThe Stefan–Boltzmann law gives $T_e^4\propto L/R^2\propto {\mu }^{25/4}{M}^{15/4}$. Eliminating mass, rather than dropping composition dependence, gives the Hertzsprung-Russell diagram relationThe additive constant depends on the fixed physical and structural coefficients; logarithms refer to fixed choices of units.

For physical evolution time $s$, distinct from the dimensionless temperature $t(x)$, hold stellar mass fixed and put $u=3+5X$, $u_0=3+5X_0$. The mean molecular weight approximation and the homology relation give $L=L_0(u_0/u{)}^7$. Converting hydrogen mass to helium releases energy at rate $L$, so fuel conservation is $M{E}_{H}dX/ds=-L$. ThereforeIntegration with the initial composition proves the homogeneous fuel-depletion luminosity feedback lawIt is valid while hydrogen remains and the homology assumptions hold. Fuel exhaustion occurs when $u=3$, before the formal denominator vanishes; the apparent infinite luminosity is not a physical endpoint of this stellar model.

Let $\Omega =-{\int }_0^{M}Gmdm/r$ be the gravitational energy. Multiply the stellar hydrostatic equation by $4\pi r^3$ and integrate through the star. Integration by parts and zero surface pressure givewhile its hydrostatic right-hand side is $-{\int }_0^{R}4\pi Gmr\rho dr=\Omega$. Hence the stellar virial theorem is

There is a genuine reference-sheet error to resolve. A monatomic ideal gas with adiabatic exponent $5/3$ has specific internal energy of an ideal gas $u=P/[\rho (5/3-1)]=3P/(2\rho )$, not the printed $3P/\rho$ on PDF page 2. The latter would give total internal energy $-\Omega$ and zero total energy, inconsistent with radiating contraction. The following derivation uses the consistent monatomic value.

A fully convective star is approximately isentropic and has an adiabatic stellar polytrope with index $3/2$. The given polytropic gravitational-energy formula yields $\Omega =-6G{M}^2/(7R)$. The virial relation gives total internal energy $U=-\Omega /2$, so total energy is $U+\Omega =-3G{M}^2/(7R)$. At constant mass with no nuclear heating, Kelvin-Helmholtz contraction obeys $d(U+\Omega )/ds=-L$, whenceHere $s$ denotes physical time.

In the plane-parallel atmosphere, the inward optical depth satisfies $d\tau =-\kappa \rho dr$, so hydrostatic equilibrium in optical depth gives $dP/d\tau =g/\kappa$. Neglecting radiation pressure, $\rho =\mu P/(\mathcal{R}T)$ and the specified opacity becomes $\kappa =({\kappa }_0\mu /\mathcal{R})P{T}^{16}$. ThereforeThe grey atmosphere relation and the boundary $P(0)=0$ integrate toTo apply the Schwarzschild criterion, put $w=2+3\tau$. Logarithmic differentiation givesThis increases inward, starting from zero. Marginal convection has ${\nabla }_{\mathrm{rad}}=2/5$, so $w_b^3=136/5$. Since $T_b^{12}/T_e^{12}=w_b^3/64$, the grey-atmosphere convection onset with seventeenth-power opacity isThe purely radiative atmosphere is used only up to this boundary.

The convective relation $P=K{T}^{5/2}$ is spatially uniform in its entropy coefficient $K$ at each instant; $K$ can change as the star radiates and contracts. At the convection boundary $w_b$ is a fixed number, so $P_b\propto g^{1/2}{T}_e^{-8}$ at fixed composition and opacity normalization. ThusOn the other hand, the central scalings for the $n=3/2$ stellar polytrope give $P_c\propto G{M}^2/R^4$ and $T_c=\mu P_c/(\mathcal{R}{\rho }_c)\propto \mu GM/(\mathcal{R}R)$. Hence, at fixed composition,Equating the interior and atmosphere expressions yields $T_e^{21/2}\propto M{R}^{1/2}$. Applying the Stefan–Boltzmann law then proves the Hayashi contraction with seventeenth-power opacity relation

At fixed mass write $L=A{M}^{8/21}{R}^{\alpha }$, with constant $A$ and $\alpha =46/21$. The contraction equation is $\dot{R}=-C{R}^{\alpha +2}$, where $C=7A{M}^{8/21}/(3G{M}^2)$. Direct integration givesMultiplying through by $R^{\alpha +1}$ and substituting the current luminosity givesFor the stated very large initial radius, or after that initial-radius term becomes negligible,This is the luminosity-age relation within the contraction model, rather than an assertion that a physical star begins at literally infinite radius.

For an isotropic degenerate electron gas, pressure is one third of the momentum flux. In the ultrarelativistic limit each occupied mode has speed $c$, so the supplied momentum-space state count givesA fully ionized helium nucleus provides two electrons for approximately $4m_p$ of mass, so $n_e=\rho /(2m_p)$. Eliminating the Fermi momentum $p_0=(3h^3\rho /(16\pi m_p){)}^{1/3}$ therefore givesThe two-spin-state degeneracy is already included in the supplied state count; no further factor 2 is needed.

If electron degeneracy pressure supplies the star's pressure, this is an index-three stellar polytrope. Its central pressure and density obeyThe radius cancels, leaving $M^{2/3}=K(3a_3/(4\pi ){)}^{4/3}/(c_{3}G)$. Thus the Chandrasekhar mass from central polytropic constants isThe dependence on microscopic constants reflects the balance of fermionic pressure and gravity. The absence of radius dependence is special to the limiting $4/3$ pressure exponent. Ordinary partially relativistic white dwarfs can have a range of masses below this Chandrasekhar limit; the formula does not force every white dwarf to have this mass. As a degenerate stellar core grows toward the limiting mass, compression cannot provide an arbitrarily larger supporting pressure relative to gravity. Above it the ideal cold, nonrotating electron-supported configuration cannot remain in equilibrium; nuclear ignition, collapse or other changes in the support mechanism must intervene. Which outcome occurs depends on the core's composition and evolutionary state.

Outside the thin burning shell, the enclosed mass and luminosity are approximately $M_c$ and $L$. Divide the hydrostatic pressure support equation by radiative diffusion in a star with constant opacity:In the constant-opacity radiative-zero giant envelope approximation, the outer pressure and temperature constants are negligible. Integration givesThe radiation pressure is $a{T}^4/3$, so the gas pressure fraction $\beta$ is constant and satisfies $(1-\beta )\Lambda =a/3$. Thus alsoThe physical branch has positive gas pressure, hence $L<4\pi cG{M}_c/{\kappa }_0$, the constant-opacity Eddington luminosity for the core mass.

Using the gas part of the equation of state, define $H=\mu \beta \Lambda /\mathcal{R}$, so $\rho =H{T}^3$. Substitution into hydrostatic balance givesConsequently, in the same deep-envelope radiative-zero approximation,Precisely, a finite outer matching point $(R_s,T_s)$ gives $T=T_s+\mu \beta G{M}_c(1/r-1/R_s)/(4\mathcal{R})$. The displayed $1/r$ profile neglects that additive outer-boundary constant; it is appropriate near a compact core when the envelope extends to much smaller temperatures.

To estimate the thin-shell luminosity, write $K_T=\mu \beta G{M}_c/(4\mathcal{R})$, so $T=K_T/r$. The hydrogen-layer heating law givesThis steep radial dependence strongly concentrates the emission near $R_c$. Continuing the outer profile across the thin emitting layer and taking a distant upper limit givesThe same radiative coefficient implies $H=16\pi ac{K}_T/(3{\kappa }_{0}L)$. Hence the high-power shell burning in a constant-opacity giant envelope satisfiesTaking the cube root gives the requested form,Thus $C$ contains only the specified opacity and heating normalizations and physical constants. The explicit coefficient uses the outer radiative profile as a thin-shell estimate; resolving the variation of $L_r$ within the emitting layer can change that coefficient. Also $\beta$ depends on $L$, so the radiation-pressure-inclusive relation is implicit, not a pure explicit eighth-power core-mass law with independently fixed $\beta$.

Put $M=M_1+M_2$. In a circular orbit, the distances from the centre of mass are $r_1=a{M}_2/M$ and $r_2=a{M}_1/M$. Summing the components' angular momentum givesUse Kepler's third law ${\Omega }^2{a}^3=GM$ and $\Omega =2\pi /P_{\mathrm{orb}}$ to obtain the circular-binary orbital angular momentumThe numerator contains orbital period, not the erroneous $M_{\mathrm{orb}}$ token in the converted TeX.

The total mass changes as $\dot{M}=(1-f){\dot{M}}_1$. With the prescribed wind carrying specific angular momentum $\lambda J/M$, the angular-momentum loss law is $\dot{J}/J=\lambda \dot{M}/M$. For constant $\lambda$, integration gives $J\propto M^{\lambda }$. Cubing the orbital formula then proves the binary period invariant for a constant wind angular-momentum factor:This integration does not require $f$ to be constant, although its instantaneous value enters the following response coefficients. In the conservative case $f=1$, both $M$ and $J$ are constant and the usual period-product invariant for conservative mass transfer is recovered.

Kepler's law and the given Roche lobe approximation imply $R_L\propto M_1^{1/3}{P}_{\mathrm{orb}}^{2/3}$, since the total-mass factor cancels. Using the period law therefore givesBecause ${\dot{M}}_2=-f{\dot{M}}_1$, the instantaneous Roche-lobe radius response exponent isFor $x=\log (R/R_L)$, subtracting this response from the prescribed stellar radius evolution givesThe overflow-rate model requires $f>0$ and is used for small positive overfill. It gives ${\dot{M}}_1/M_1=-x/(f{t}_{\mathrm{dyn}})$. Hence the nuclear-driven Roche-overfill relaxation equation is

Let the bracket be $A(q)$. On the short response time, freeze the slowly changing masses and coefficients. If $A>0$, the local solution isThe equilibrium donor mass-loss rate is ${\dot{M}}_1/M_1=-1/(A{t}_{\mathrm{nuc}})$. Thus the required stabilizing mass-ratio condition isWhen $A$ is of order one and $t_{\mathrm{dyn}}\ll t_{\mathrm{nuc}}$, the overfill is small and relaxation is fast, while mass transfer proceeds on a nuclear timescale. A positive $A$ arbitrarily close to zero does not by itself ensure this timescale separation or small overfill.

For an explicit threshold put $B=n+5/3$ and $C_w=2\delta (1-f)/3$. The inequality isIf $B>0$ and $f>0$, its positive root givesFor $f=1$ this reduces to $q<(n+5/3)/2$, the conservative mass-transfer critical mass ratio. For other parameter signs the unsimplified inequality is the general condition.

If $A<0$, a larger overfill causes further growth rather than restoration, on the short scale $f{t}_{\mathrm{dyn}}/|A|$. Transfer runs away towards a dynamical episode, potentially a common envelope or merger; the small-overfill model then fails. At $A=0$ the overfill grows linearly under nuclear expansion and there is no steady regulated contact. These conclusions use the stipulated donor-radius response; an actual dynamical episode requires its adiabatic response and hydrodynamics.