Past exam of the mathematics course of the University of Cambridge / 2012 / iii / Paper 63 / 1

Efficient convection makes the interior follow a nearly adiabatic process and maintains uniform specific entropy. For the fully ionized monatomic ideal gas, the specific-heat ratio is $5/3$, so the adiabatic temperature gradient gives $d\log T/d\log P=2/5$. Integrating at a fixed time givesThe coefficient is spatially constant because the composition and specific entropy are uniform. It can change as the star loses energy. Equivalently the gas is an adiabatic stellar polytrope with $P=K_{\rho }{\rho }^{5/3}$ and polytropic index $n=3/2$; $K_{\rho }$ and the pressure-temperature coefficient $K$ are different quantities.

This fixes the dimensionless Lane-Emden equation profile throughout a stellar homology sequence. Writing $\rho (r)={\rho }_{c}f(r/R)$ in the mass integral gives ${\rho }_c=C_{\rho }M/R^3$, with a fixed positive shape constant $C_{\rho }$. Integrating the stellar hydrostatic equation through the same profile, with negligible surface pressure compared with central pressure, gives $P_c=C_{P}G{M}^2/R^4$. HenceThese are homology scalings of central stellar pressure. Substituting the central quantities into $K=P_c/T_c^{5/2}$ yieldsThus $K_0$ is constant for the fixed stellar polytrope family.

At the approximate photosphere, take $T=T_e$ and use the ideal gas relation $\rho =\mu P/(\mathcal{R}T)$. The prescribed opacity givesThe stellar surface boundary condition and surface gravity of a star therefore implyThis is the temperature-four special case of convective contraction with power-law opacity.

Using the Stefan–Boltzmann law, the luminosity isFor $n=3/2$, the supplied gravitational-energy formula is $\Omega =-6G{M}^2/(7R)$. The stellar virial theorem for a monatomic gas gives $2U+\Omega =0$, so the total energy is $E=U+\Omega =-3G{M}^2/(7R)$. During Kelvin-Helmholtz contraction, with constant mass and no significant nuclear source, $L=-dE/dt$. ThusMultiplying by $-(7/2)R^{-9/2}$ and integrating givesThe sign makes the radius decrease as time increases.

When the accumulated contraction term dominates the initial-radius term, this yields $R\propto {\mu }^{-4/7}{M}^{2/7}(t-t_0{)}^{-2/7}$. Combining it with the central ideal gas temperature gives the late-time central-temperature scalingReplacing $t-t_0$ by $t$ gives the stated late-time form. The asymptotic step also requires a time long enough to forget the initial radius, not merely a choice of zero for the age.

The helium mass fraction is $Y=1-X-Z$. For the two compositions, $Y_1=0.28$ and $Y_2=0.18$. The mean molecular weight is thereforeAt fixed mass and with the same prescribed opacity coefficient, part (a) gives $T_e\propto {\mu }^{1/2}{R}^{1/8}$ and $L\propto {\mu }^2{R}^{5/2}$. Eliminating radius givesThis is the composition shift of a convective pre-main-sequence track. The higher-helium mass fraction case $X_1=0.7$ is hotter at a given luminosity, lying to the left of the $X_2=0.8$ case in a Hertzsprung-Russell diagram. Its fixed-luminosity temperature ratio is $(349/324{)}^{2/5}\simeq 1.030$.

Both stars contract, so both $L$ and $T_e$ decrease in this particular opacity model. Their nearly vertical Hayashi tracks therefore run downward and slightly to the right, with arrows in that direction:

The plotted luminosity and effective temperature are normalized to arbitrary common references; no opacity normalization or absolute age was specified. The comparison holds the given opacity coefficient fixed.

These paths cease to apply when their assumptions fail. Deuterium-burning energy-generation rate can temporarily slow contraction; sustained hydrogen burning ends the purely gravitational phase. Formation of a radiative core breaks the fully convective adiabat and can lead to a Henyey track toward hotter effective temperature. Partial ionization, a changing atmospheric opacity, and electron degeneracy pressure alter the equation of state or boundary relation. Accretion luminosity, changing mass, and magnetic inhibition of convection can also modify a young star's path.

Keep the fixed stellar homology profile, so that ${\rho }_c=C_{\rho }M/R^3$ and $P_c=C_{P}G{M}^2/R^4$. Solving the partially degenerate equation of state for the central temperature giveswhere $A=C_P/C_{\rho }$ and $B=C_{\rho }^{2/3}$ are positive constants. The ion mean molecular weight is distinct from the total-particle mean molecular weight used in part (a).

Set $y=1/R$. The expression in parentheses is a concave quadratic in $y$. Its maximum occurs atSubstitution gives the contraction-limited central temperatureAt this point electron degeneracy pressure supplies half the required central pressure. Further contraction increasingly goes into the density-dependent degeneracy term rather than heating the ions. The formula applies along the positive-temperature part of the idealized homologous sequence; a negative value beyond it is not a physical temperature.

For sufficiently small mass, this maximum never reaches the temperature needed for sustained hydrogen burning. Such an object stays below the hydrogen-burning minimum mass and becomes a cooling brown dwarf, rather than an ordinary main sequence star. Some such objects burn their initial deuterium, but that finite reservoir does not provide lasting hydrogen burning.