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

Use $\mathcal{R}$ for the stellar gas constant, reserving $R$ for the stellar radius. The specific internal energy of a monatomic gas obeying the ideal gas law is $e=P/[(\gamma -1)\rho ]=3P/(2\rho )$. At fixed enclosed mass, the first law of thermodynamics gives the heat available to the outgoing luminosity as $ϵ=-{\partial }_{t}e-P{\partial }_t(1/\rho )$ when there is no stellar nuclear fusion heating. ThereforeThe second coefficient includes both the mass density dependence of $e$ and pressure work. This is gravitational energy generation in a homologously contracting ideal-gas star.

Write the characteristic scales asThen $\rho ={\rho }_{0}b$, $P=P_{0}p$, $T=T_{0}p/b$, $m=Mq$ and $L_r=Ll$. Since $d/dr=R^{-1}d/dx$, direct substitution in the stellar hydrostatic equation and mass conservation givesFor radiative diffusion in a star, dividing the temperature equation by $T_0/R$ givesIn particular the fourth power here is of $\mathcal{R}$, not of the radius.

Because $q(x)$ is fixed and $M$ is constant, a fixed mass label has fixed $x$. Thus stellar homology gives ${\partial }_t\rho =-3\rho \dot{R}/R$ and ${\partial }_{t}P=-4P\dot{R}/R$. The heating rate reduces to $ϵ=-3P\dot{R}/(2\rho R)$, positive during Kelvin-Helmholtz contraction. Substitution in $d{L}_r/dr=4\pi r^2\rho ϵ$ now yieldsThe dimensionless profiles fix $D$ and $E$. In a common stellar homology family with fixed stellar composition and opacity normalization, the expression for $D$ impliesIt is independent of radius and constant in time for a particular fixed-mass star. The formula for $E$ gives $\dot{R}=-2EL{R}^2/(3G{M}^2)$, which integrates toWhen the initial-radius term is negligible,This is the Kelvin-Helmholtz contraction time law in the prescribed stellar homology model. Infinite initial radius is a limiting idealization; a finite initial radius retains the first term.

At main sequence arrival, CNO cycle heating with fixed stellar composition has luminosityThe integral is dimensionless and fixed within the stellar homology family. The constant-opacity radiative scaling still gives $L\propto M^3$. Equating the stellar nuclear fusion and transported luminosities givesFinally the Kelvin-Helmholtz contraction law evaluated at $R_{\mathrm{MS}}$ gives $t_{\mathrm{MS}}\propto M^2/(L{R}_{\mathrm{MS}})\propto M^{-34/19}$. In a coeval star cluster, the mass just arriving at the main sequence consequently satisfiesThese results comprise constant-opacity homologous contraction to CNO ignition; their mass exponents compare fixed-composition models with common dimensionless profiles.

Consider a small outward displacement $\delta r$ of a fluid element. It adjusts to the ambient pressure while retaining its specific entropy and stellar composition, so an adiabatic process gives $d\log {\rho }_{\mathrm{parcel}}=d\log P/\gamma$. Its excess mass density over the new surroundings is, to first order,Restoring buoyancy requires this bracket to be positive. With $\gamma =5/3$, the Schwarzschild criterion is thereforeEquality is neutral stability. For uniform mean molecular weight, the ideal gas law gives $d\log \rho =d\log P-d\log T$. Since hydrostatic equilibrium pressure decreases outward, the inequality is equivalent to $\nabla =d\log T/d\log P<{\nabla }_{\mathrm{ad}}=2/5$. Dividing the radiative diffusion in a star and hydrostatic equilibrium stellar structure equations givesThis is the stellar radiative temperature gradient test for stellar convective stability.

In the thin upper stellar atmosphere take $m\simeq M$, $L_r\simeq L$ and use gas pressure dominance. The prescribed opacity becomes $\kappa ={\kappa }_0\mu P{T}^{12}/\mathcal{R}$. Dividing the two structure equations then yieldsAt zero optical depth, the atmospheric boundary has $T_s^4=T_e^4/2$ and negligible gas pressure. Integrating from that boundary, with $T_s^{-8}=4T_e^{-8}$, givesTo locate the onset, put $y=T_e^8/T^8$. Logarithmic derivative gives $d\log P/d\log T=4y/(4-y)$, henceThis starts at zero at the outer boundary and increases inward. It first reaches $2/5$ atThe corresponding optical depth is ${\tau }_b=\frac{4}{3}(\sqrt{13/20}-\frac{1}{2})>0$, so the onset lies inside the stellar atmosphere. This is grey-atmosphere convection onset with thirteenth-power opacity.

Let the interior of the fully convective star have $P=K_T{T}^{5/2}$, with $K_T$ spatially constant. At its boundary $T_b/T_e$ is fixed, so the atmospheric pressure formula givesThe Stefan–Boltzmann law, $L=4\pi R^2\sigma T_e^4$ with $\sigma =ac/4$, consequently impliesThe ideal gas equation transforms the interior relation intoFor completeness, the mass-radius scaling follows by taking $\rho ={\rho }_c{\theta }^{3/2}$ and $r=\alpha \xi$ with ${\alpha }^2=5K_{\rho }{\rho }_c^{-1/3}/(8\pi G)$. Hydrostatic equilibrium and mass conservation reduce to the Lane-Emden equation $({\xi }^2{\theta }^{\prime }{)}^{\prime }=-{\xi }^2{\theta }^{3/2}$, with $\theta (0)=1$, ${\theta }^{\prime }(0)=0$. Its fixed dimensionless profile gives $R\propto \alpha$ and $M\propto {\alpha }^3{\rho }_c$. Eliminating ${\rho }_c$ gives $R\propto K_{\rho }{G}^{-1}{M}^{-1/3}$, or $K_{\rho }\propto G{M}^{1/3}R$. Thus $K_T^2\propto M^{-1}{R}^{-3}$ at fixed stellar composition. Equating this with the atmospheric expression gives $L^{17/4}\propto M^2{R}^{19/2}$, and henceThis Hayashi relation with thirteenth-power opacity concerns the interior of a fully convective star beneath a thin stellar atmosphere. The coefficient $K_T$ may vary between stars; spatially constant specific entropy does not mean equal specific entropy across the whole sequence.

A fully ionized helium nucleus supplies two electrons for approximately $4m_p$ of mass, so $n_e=\rho /(2m_p)$. Filling all electron momentum states up to the Fermi momentum givesFor a nonrelativistic electron, the speed is $p/m_e$. The three mean squared Cartesian components of momentum are equal by rotational symmetry and sum to $p^2$, so each is $p^2/3$. This is the isotropic tensor integral. Thus the momentum flux, or electron degeneracy pressure, isSubstituting the expression for $p_0$ and simplifying the powers of two givesHere $h$ is the Planck constant. This is the pure-helium specialization of nonrelativistic electron-degeneracy pressure in a hydrogen-helium mixture.

For a cold helium white dwarf, neglect ion thermal pressure and use this electron degeneracy pressure as the support against Newtonian gravity. Combining hydrostatic equilibrium with mass conservation givesSet $\rho ={\rho }_c{\theta }^{3/2}$ and $r=\alpha \xi$, where ${\alpha }^2=5K{\rho }_c^{-1/3}/(8\pi G)$. Direct substitution gives the Lane-Emden equationIts dimensionless solution is independent of ${\rho }_c$. Let ${\xi }_1$ be its first zero and ${\omega }_1=-{\xi }_1^2{\theta }^{\prime }({\xi }_1)={\int }_0^{{\xi }_1}{\xi }^2{\theta }^{3/2}d\xi$. Integrating mass conservation givesEliminating ${\rho }_c$ yields the polytropic mass-radius relationThus a more massive nonrelativistic white dwarf is smaller at fixed stellar composition. The assumptions require $p_0\ll m_{e}c$ and $k_{B}T\ll p_0^2/(2m_e)$ in the region providing most support; the law is not extended into the relativistic or thermally supported regimes.

Outside the emitting layer, take $m\simeq M_c$, $L_r\simeq L$ and $P=\mathcal{R}\rho T/\mu$. The Kramers' opacity law then reads $\kappa =({\kappa }_0\mu /\mathcal{R})P{T}^{-9/2}$. Division of the radiative diffusion in a star equation by the stellar hydrostatic equation givesWith the outer pressure constant of integration neglected, integration yieldsSubstitute this into hydrostatic equilibrium:Define $B=4\mu G{M}_c/(17\mathcal{R})$. The general antiderivative is $T=B/r+T_0$. The Kramers radiative-zero envelope around a stellar core sets $T_0=0$, givingThis is exact for the idealized boundary $T\to 0$ as $r\to \infty$. With a finite outer radius $R_s$ and temperature $T_s$, instead $T=T_s+B(1/r-1/R_s)$; the displayed profile is the deep-envelope approximation when $R_s\gg R_c$ and $|T_s-B/R_s|\ll B/R_c$. Small outer pressure alone does not eliminate this temperature constant of integration.

The specific hydrogen burning rate becomesConsequentlyIts local radial e-folding length is $R_c/20$, so burning is concentrated within a small fraction of the stellar core radius. The volume heating is even steeper: $\rho ϵ\propto C^2{T}^{93/4}$. Using the exterior envelope profile to estimate the thin shell's luminosity givesA finite but very extended upper limit changes the factor by $1-(R_c/R_s{)}^{81/4}$. The large exponent makes this correction small and further verifies the hydrogen burning thin-shell approximation. Since $C^2\propto M_c/L$ and $B\propto M_c$, the integral implies $L^2\propto M_c^{97/4}{R}_c^{-81/4}$. HenceThis is the fixed-radius-core shell-burning luminosity relation. The proportionality coefficient is a thin-shell estimate: luminosity rises from its value beneath the shell to $L$ through the emitting layer, so treating it as constant there uses the exterior profile rather than resolving the burning region. The scaling holds while its dimensionless shell structure and stellar composition remain fixed.

For consistency, $M_c/R_c$ must provide a hot hydrogen burning shell while the stellar core remains inert to core helium burning; the radiative envelope mass must be much less than $M_c$, and $R_c$ must be well inside an extended radiative envelope. At the shell base the gas must remain a nondegenerate gas and dominated by gas pressure, in particular $a{T}_b^4/3\ll C{T}_b^{17/4}$ with $T_b=B/R_c$. The assumed stellar core sequence must genuinely have nearly fixed $R_c$ over the mass interval considered. A cold helium stellar core supported by a nonrelativistic degenerate electron gas instead has the mass-radius law from part (a), so substituting that law would describe a different model and change the luminosity exponent.

Let $M=M_1+M_2$. The centre of mass condition gives $a_1=a{M}_2/M$, $a_2=a{M}_1/M$. For circular motion the component angular momenta are $J_i=M_i{a}_i^2\Omega$, soUsing Kepler's third law, $a^3{\Omega }^2=GM$, and $\Omega =2\pi /P_{\mathrm{orb}}$ givesThis is the circular-binary orbital angular momentum. Mass received by the companion redistributes angular momentum within the binary star; only the escaping stellar wind removes it under the prescribed model.

The stellar wind from the donor star has specific angular momentumHere ${\dot{M}}_1<0$, ${\dot{M}}_2=-f{\dot{M}}_1$ and $\dot{M}=(1-f){\dot{M}}_1$. Thus donor-wind angular-momentum loss givesTaking the logarithmic derivative of the period expression for $J$ givesUse $M_2/(M_{1}M)=1/M_1-1/M$ to simplify this toFor constant $f$, integration therefore givesThis is the donor-wind period invariant with fixed retention fraction. If $f$ changes during the evolution, the differential equation remains valid but the integrated invariant is instead $P_{\mathrm{orb}}{M}_2^3{M}^2\exp (3\int fd\log M_1)=\mathrm{constant}$. The fixed-power expression cannot be used with a time-dependent exponent without this modification.

For the Roche lobe, write $q=M_1/M_2$. The equivalent separation expression $J=M_1{M}_2\sqrt{Ga/M}$ givesTaking the logarithmic derivative of $R_L=0.46a(M_1/M{)}^{1/3}$ then adds ${\dot{M}}_1/(3M_1)-\dot{M}/(3M)$, yieldingThis donor-wind Roche-lobe response uses only the instantaneous value of $f$ and thus still holds if it varies.

Define $A(q)=2q+4q/[3(1+q)]-2$ and $B(q)=1/3-4q/[3(1+q)]$. The Roche-lobe radius response exponent is ${\zeta }_L=B+Af$, whereas the donor star's stellar radius response exponent is ${\zeta }_{\ast }=-n$. Keeping $R_1=R_L$ requires the two logarithmic derivatives of radius to agree, soFor $A\ne 0$, this determinesA physical mixture of Roche-lobe overflow and stellar wind from the donor star needs $0<f<1$. Equivalently, $-n$ must lie strictly between the pure-wind response $B$ and the conservative response $B+A=2q-5/3$. The values $f=0$ and $1$ are valid limiting prescriptions, respectively pure stellar wind and conservative mass transfer, rather than mixtures. If $f_{\mathrm{req}}$ lies outside $[0,1]$, no allowed division of the lost mass can maintain contact under this radius law and angular-momentum-loss prescription. This is feasibility of donor-wind binary contact.

The likely direction away from contact follows without assigning a nonphysical fraction. For any actual $f$, let $\Delta =\log (R_1/R_L)$. ThenSince ${\dot{M}}_1<0$, a positive value of $A(f_{\mathrm{req}}-f)$ makes the donor star move inside its lobe and detach; a negative value increases overfill and promotes stronger transfer. Thus if $A>0$, $f_{\mathrm{req}}<0$ gives increasing overfill and $f_{\mathrm{req}}>1$ gives detachment; if $A<0$, these outcomes are reversed. A runaway or a new equilibrium would require the donor star's response or the angular-momentum-loss model to change, possibly with additional driving. The contact calculation alone does not prove a particular nonlinear endpoint.

There is also a genuine exceptional denominator: $A=0$ at $q=(\sqrt{10}-1)/3$. At this ratio both limiting lobe responses coincide. Contact is possible only if $-n=B$, namely $n=(\sqrt{10}-2)/(\sqrt{10}+2)$; then every $f$ gives the same instantaneous radius response. Otherwise no fraction works, and the sign of $(-n-B){\dot{M}}_1/M_1$ determines detachment or increasing overfill.