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

Efficient convection keeps the interior at nearly uniform specific entropy. For the fully ionized monatomic gas, the heat capacity ratio is $5/3$, so the adiabatic process gives $P=K_{\rho }{\rho }^{5/3}$, with $K_{\rho }$ spatially constant. The ideal gas relation $P=\mathcal{R}\rho T/\mu$ then impliesThe constant is spatial, not temporal: loss of specific entropy during contraction changes it. This is a stellar polytrope with polytropic index $n=3/2$. We neglect radiation pressure and take the thin photosphere as an atmospheric matching layer rather than a significant fraction of the stellar mass.

For a fixed dimensionless profile of the Lane-Emden equation, the stellar homology relations have the formwhere $C_{\rho },C_P$ are positive dimensionless structure constants. More explicitly, write $\rho (r)={\rho }_{c}f(x)$ and $m(r)=Mh(x)$ with $x=r/R$ and fixed dimensionless profiles. Mass conservation gives $M=4\pi {\rho }_c{R}^3{\int }_0^{1}f(x)x^{2}dx$. Integrating the hydrostatic pressure support equation, with negligible surface pressure on the interior scale, gives $P_c=(GM{\rho }_c/R){\int }_0^{1}h(x)f(x)x^{-2}dx$. Both integrals are fixed positive numbers, proving the two central scalings. Substituting these relations into $K=P_c/T_c^{5/2}$ givesThus the same polytropic index fixes $K_0$ throughout this homologous sequence, while the dimensional stellar polytrope constant changes.

At the photosphere, whose optical depth is $2/3$, identify the temperature with the effective temperature $T_e$. The ideal gas relation and the interior adiabat give ${\rho }_{\mathrm{ph}}=\mu K{T}_e^{3/2}/\mathcal{R}$. With the prescribed opacity and stellar surface boundary condition,Eliminating $K$ therefore provesUsing the Stefan–Boltzmann law, with $\sigma =ac/4$, gives the luminosity

The gravitational energy of a stellar polytrope is $\Omega =-6G{M}^2/(7R)$ for $n=3/2$. The stellar virial theorem for a monatomic gas gives $2U+\Omega =0$, and hence the total energy is $E=U+\Omega =-3G{M}^2/(7R)$. With no appreciable stellar nuclear fusion or accretion, energy conservation givesCombining the two expressions for $L$ yieldsMultiplication by $-7R^{-9/2}/2$ and integration proves the contraction law with density-linear fourth-power stellar opacity:When the accumulated contraction term dominates the initial-radius term, $R\propto {\mu }^{-4/7}{M}^{2/7}(t-t_0{)}^{-2/7}$. Substitution into the central ideal gas relation giveswhich has the stated $t^{2/7}$ behavior when the time origin is negligible. This late-time approximation still requires the assumed fully convective star, fixed composition and opacity law to remain applicable.

For fully ionized hydrogen and helium with mass fractions $X$ and $Y=1-X$, the inverse mean molecular weight counts nuclei plus free electrons:Thus ${\mu }_1=8/13$ for $X_1=0.7$, ${\mu }_2=4/7$ for $X_2=0.8$, and ${\mu }_1/{\mu }_2=14/13$. Keep the same $M,{\kappa }_0$ and dimensionless structure constants when comparing these toy sequences. From part (a),Eliminating $R$ gives the composition shift of a convective pre-main-sequence track:Consequently the tracks in a Hertzsprung-Russell diagram are almost vertical, with $d\log L/d\log T_e=20$. At equal luminosity, $T_{e,1}/T_{e,2}=(14/13{)}^{2/5}\simeq 1.030$: the $X=0.7$ track lies to the hotter, left-hand side of the $X=0.8$ track. Since contraction decreases $R$, both $L$ and $T_e$ decrease in this particular opacity model. Both arrows therefore point downward and slightly rightward on the conventional hot-left diagram.

The plotted reference scales are arbitrary; the horizontal separation is the predicted $\Delta {\log }_{10}{T}_e=(2/5){\log }_{10}(14/13)$. They are schematic Hayashi tracks under the specified atmospheric opacity law, rather than calibrated stellar models.

A protostar leaves this contraction sequence when its assumptions fail. Increasing central stellar temperature makes stellar nuclear fusion significant: burning of deuterium can temporarily slow contraction, and sustained hydrogen burning brings the star to the main sequence. In a sufficiently massive pre-main-sequence star, a radiative core develops, the interior ceases to share one adiabat, and evolution can proceed along a hotter Henyey track. Changing atmospheric opacity, partial ionization, or substantial continuing accretion also changes the predicted path. For sufficiently low mass, electron degeneracy pressure instead limits the increase of central temperature, as in part (c).

Retain fixed dimensionless structure constants and fixed $M$. Write the assumed stellar homology relations as $P_c=C_{P}G{M}^2/R^4$ and ${\rho }_c=C_{\rho }M/R^3$. The equation of state giveswhere $A=C_P/C_{\rho }$ and $B=C_{\rho }^{2/3}$. The first term is the temperature required for ordinary thermal pressure support; the second subtracts the support supplied by nonrelativistic electron degeneracy pressure.

As a function of $u=1/R$, this is a concave quadratic polynomial. Differentiation givesand substitution provesThe scaling holds at fixed $K_{\mathrm{nr}}$ and homologous profile; changing the mean molecular weight per electron also changes $K_{\mathrm{nr}}$. At the maximum the ion thermal pressure and electron degeneracy pressure contribute equally. Subsequent contraction lets electron degeneracy pressure carry more of the load without raising the temperature; eventually the assumed positive-temperature branch approaches a cold, degenerate configuration.

A sufficiently low-mass object never reaches the temperature needed for sustained hydrogen burning, so it becomes a cooling brown dwarf rather than a hydrogen-burning star on the main sequence. The scaling explains a hydrogen-burning minimum mass; its numerical value additionally requires the composition, detailed structure and nuclear ignition conditions, which are not specified here.