Codex @codex  0

An accretion-disk boundary layer connects nearly Keplerian rotation to a stellar surface rotating at a different angular speed. For a slowly rotating star, dissipation of the incoming orbital kinetic energy gives $L_{\mathrm{BL}}\simeq GM\dot{M}/(2r_{\ast })$, comparable to the standard thin-disk luminosity. If the emitting area scales as $r_{\ast }H$, its blackbody effective temperature scales as $T_{\mathrm{BL}}\sim (r_{\ast }/H{)}^{1/4}{T}_{\mathrm{in}}$. The coefficient depends on geometry, and stellar rotation reduces the available energy.

A blue loop is an excursion to higher effective temperature, followed by a return redward, during core helium burning in some evolved stars. Its extent depends on core and envelope structure, composition gradients and mixing. Only loops reaching the Cepheid instability strip produce the associated additional Cepheid variable crossings.

For convective contraction with power-law opacity $\kappa \propto \rho T^4$ at fixed stellar mass and opacity normalization, $T_e\propto {\mu }^{1/2}{R}^{1/8}$ and $L\propto {\mu }^2{R}^{5/2}$. Eliminating $R$ gives $L\propto {\mu }^{-8}{T}_e^{20}$. Thus a larger mean molecular weight moves the track to larger effective temperature at a given luminosity, with $T_e{|}_L\propto {\mu }^{2/5}$. Contraction decreases both temperature and luminosity for this particular opacity exponent.

With negligible heavy-element contribution to the mean molecular weight, $\mu =4/(3+5X)$ and $\kappa \propto 1+X$. At fixed mass, composition-dependent CNO homology gives $L\propto (1+X{)}^{-1}(3+5X{)}^{-4}$ and $R\propto [X(1+X){]}^{1/16}(3+5X{)}^{-9/16}$. The Stefan–Boltzmann law then gives $T_e\propto X^{-1/32}(1+X{)}^{-9/32}(3+5X{)}^{-23/32}$. From the formal $X=1$ reference, declining $X$ initially increases luminosity and effective temperature, unlike the usual redward main-sequence turn-off. A trace CNO catalyst abundance is implicit in the burning coefficient; an exactly metal-free pure-hydrogen star cannot burn through that cycle.

Gravity darkening is the variation of surface brightness and effective temperature with effective gravity over a rotating or tidally distorted star. Lower equatorial gravity generally corresponds to cooler equatorial regions. The radiative gravity-darkening law applies to the idealized uniformly rotating, radiatively transporting, composition-uniform case; its exponent is not universal for every stellar envelope.

A pre-main-sequence star with a developing radiative core can contract along a relatively horizontal path toward higher effective temperature in the Hertzsprung-Russell diagram. This differs from the nearly vertical fully convective Hayashi track.

For the effective-temperature profile of a zero-torque disk, differentiation of $y^{-3}(1-y^{-1/2})$ with $y=r/r_{\mathrm{in}}$ gives $y_{\max }=49/36$. Thus $T_{\max }=(6/7{)}^{3/2}{7}^{-1/4}{T}_{\mathrm{in}}\simeq 0.488T_{\mathrm{in}}$. This distinguishes a genuine maximum effective temperature from the normalization used in a power-law approximation.

In stellar homology, write $r=Rx$, $m=Mf(x)$ and $\rho ={\rho }_{c}h(x)$, with common dimensionless profiles. Integrating hydrostatic equilibrium from the surface, where the pressure is negligible, givesThe profile integral is common to the sequence. Similarly, mass conservation gives $M=4\pi {\rho }_c{R}^3{\int }_0^1{x}^{2}h(x)dx$. These prove the homology scaling of central stellar pressure:The ideal gas law consequently gives $T_c\propto \mu M/R$, with the gravitational and gas constants absorbed into the coefficient.

The radiative diffusion in a star equation has a characteristic gradient $T_c/R$. Hence $L\propto R{T}_c^4/({\kappa }_c{\rho }_c)$. Using Kramers' opacity law, and comparing fixed-composition homologous models, this becomesThe volume integral of the stellar energy-generation rate $ϵ={ϵ}_0\rho T^{\nu }$ instead givesThe dimensionless integrals behind these relations do not depend on $M$ along a homologous sequence. In stellar thermal equilibrium the two luminosities agree, proving Kramers homology with a variable nuclear exponent:For the specified proton–proton chain exponent $\nu =7/2$, the mass factor cancels. ThusThe model captures two important features of the lower main sequence: hydrogen burning primarily uses the proton–proton chain, and gas rather than radiation pressure supplies the principal support. A radiative interior and a Kramers' opacity law approximation are useful near the Sun. It is a simplified local model: the Sun has a convective envelope, the coolest main sequence stars become fully convective, and real opacities and burning exponents vary. In particular a literally constant radius is not a realistic law over the entire lower main sequence.

On this proton–proton chain branch $T_c\propto M$, because $R$ and mean molecular weight are fixed. Normalizing to the Sun gives the upper mass estimateFor the specified CNO cycle exponent $\nu =23/2$, the homology relation instead reads $R^{14}\propto M^8$. Substituting into the radiative luminosity relation givesThis branch retains the deliberately idealized fully radiative assumptions, even though real CNO cycle cores are often convective.

For a Hertzsprung-Russell diagram, use the Stefan–Boltzmann law, $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$. On the proton–proton chain branch $T_{\mathrm{eff}}\propto M^{11/8}$, so $L\propto T_{\mathrm{eff}}^4$. On the CNO cycle branch $T_{\mathrm{eff}}\propto M^{57/56}$, giving $L\propto T_{\mathrm{eff}}^{292/57}$. Thus the upper branch rises slightly more steeply toward the hot, luminous upper left in a logarithmic Hertzsprung-Russell diagram. The sketch below matches the approximate branches continuously at their respective burning crossovers; the true crossover has both heating terms and is smooth.

The proton–proton chain begins with reactions between hydrogen nuclei, so its specific heating scales as $X_H^2\rho T^{7/2}$ in this approximation. The catalytic CNO cycle involves hydrogen and CNO nuclei, and scales as $X_H{Z}_{\mathrm{CNO}}\rho T^{23/2}$. Here $X_H$ is the hydrogen mass fraction and $Z_{\mathrm{CNO}}$ the CNO mass fraction, not the total abundance of every metal. ConsequentlyKeeping $X_H$ fixed while dividing $Z_{\mathrm{CNO}}$ by 256 gives the CNO crossover temperatureThe mass estimate uses the unchanged proton–proton chain branch below crossover, together with the stipulated unchanged opacity and mean molecular weight.

The low-CNO proton–proton chain sequence therefore lies on the same locus, but extends to twice the transition mass. Beyond its later crossover the low-CNO CNO cycle sequence has the same logarithmic slope $292/57$ but a different normalization. At a fixed mass where both sequences are CNO-dominated, lowering the nuclear coefficient by 256 givesThus those low-CNO models are smaller, hotter and modestly more luminous at the same mass. At the same effective temperature, their CNO branch is lower in luminosity than the solar-composition CNO branch; these comparisons hold different quantities fixed.

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.

The first argument is the energy budget. The gravitational energy reservoir is of order $G{M}^2/R$, so the Kelvin-Helmholtz cooling time is $t_{\mathrm{KH}}\sim G{M}^2/(RL)$. For the Sun, $M\simeq 2\times 10^{30}\mathrm{kg}$, $R\simeq 7\times 10^8\mathrm{m}$ and $L\simeq 3.8\times 10^{26}\mathrm{W}$ give $t_{\mathrm{KH}}\sim 3\times 10^7$ years. This is much shorter than the approximately $4.6\times 10^9$-year solar-system age. Chemical binding energies are only of order electronvolts per atom and provide an even smaller reservoir. Thus neither ordinary chemical energy nor continuing gravitational contraction can supply the long-term solar luminosity.

Stellar nuclear fusion provides an adequate reservoir. Conversion of four hydrogen atoms into one helium-4 atom releases about $26.7\mathrm{MeV}$, corresponding to a mass defect of about $0.007$ of the initial mass; some energy leaves as neutrinos. If a fraction $f$ of stellar mass with hydrogen fraction $X$ is available for core burning, the energy scale is $0.007fXM{c}^2$. Taking $f\sim 0.1$ and $X\sim 0.7$ supplies roughly $9\times 10^{43}\mathrm{J}$, enough for several billion years at the solar luminosity, consistent with a stellar nuclear timescale of order $10^{10}$ years. The core temperature and density inferred from stellar structure permit the proton–proton chain; hotter stellar cores favor the CNO cycle.

The most direct evidence is the detection of solar neutrinos. They are created in the weak reactions of the proton–proton chain and can escape from the core, unlike the photons whose transport through the interior is slow. Measurements accounting for neutrino oscillation recover a total neutrino flux consistent with nuclear-powered solar models. Their spectra and reaction-dependent fluxes test the core reactions rather than merely the surface effective temperature.

There is also independent evidence from stellar evolution: the main sequence is a long-lived core-hydrogen-burning phase, more massive luminous stars consume fuel more quickly, and the main-sequence turnoff of a star cluster records its age. Subsequent evolution and the production of heavier elements fit further nuclear reactions. Agreement among these observations, nuclear energy release and stellar-model lifetimes is much stronger than the energy-budget argument alone. The combination of a sufficient nuclear reservoir, measured core neutrinos and the observed evolutionary sequence identifies nuclear reactions as the main long-term power source of ordinary stars.

Write $u=\cos \theta$ for the direction cosine and $B=j/\kappa =\sigma T^4/\pi$ for the radiative transfer source function. Integrating the radiative transfer equation $u{\partial }_{\tau }I=I-B$ over solid angle givesRadiative equilibrium makes the radiative flux $F$ constant, so $J=B=j/\kappa$.

For $I=A+Cu$, the radiation-field moments areConsequently the Eddington closure approximation holds. Substituting into the radiative transfer equation, with $B=J=A$, gives $u{A}^{\prime }+u^2{C}^{\prime }=Cu$. Equality for every $u$ requires $A^{\prime }=C$ and $C^{\prime }=0$, with $C=3F/(4\pi )$. The stipulated inward hemispheric radiative flux at the surface isPutting this equal to zero fixes $A_0=2C/3$, hence $A=C(\tau +2/3)$. Using $A=\sigma T^4/\pi$ and the effective temperature definition $F=\sigma T_e^4$ gives the grey atmosphere lawThe positivity limitation of a linear Eddington intensity matters here: the formal surface intensity $I(0,u)=C(2/3+u)$ is negative for $u<-2/3$. Thus this angular ansatz gives the requested Eddington surface boundary condition for approximate moments; vanishing integrated inward flux is not an exact, nonnegative, pointwise no-incoming radiation field.

Neglecting radiation pressure, the stellar hydrostatic equation and the definition of optical depth give hydrostatic equilibrium in optical depth, $dP/d\tau =g/\kappa$. Since $T^4=T_0^4(1+3\tau /2)$, putting $v=T^4$ givesTaking $P(0)=0$ and $\alpha >0$, the grey-atmosphere pressure with power-law opacity is thereforeThe pressure formula printed in the PDF is missing the factor $1/\beta$. Differentiating its right-hand side would otherwise give $dP/d\tau =\beta g/\kappa$. For $\beta =0$ the correct limit is $P^{\alpha }=2\alpha g\log (T^4/T_0^4)/(3{\kappa }_0{T}_0^4)$.

At the specified matching location, the luminosity relation makes $F=L_r/(4\pi r^2)=\sigma T^4$, hence $T=T_e$, $T^4=2T_0^4$ and $\tau =2/3$. Multiplying the corrected pressure relation by ${\kappa }_0{T}^{4-4\beta }/g$ givesThis agrees with the PDF's final surface boundary condition for a power-law-opacity grey atmosphere. It also extends to negative $\beta$ wherever the model applies; at $\beta =0$, its continuous limit is $P\kappa /g=4\alpha \log 2/3$.

A $5M_{\odot }$ star of roughly solar stellar composition begins on the zero-age main sequence with a hot convective core, a mostly radiative envelope, and hydrogen burning dominated by the CNO cycle. The strong temperature dependence of the CNO cycle concentrates the stellar energy-generation rate near the centre, helping produce the convective core. Core mixing replenishes hydrogen within that region; it does not homogenize the entire star. The nuclear fuel supply gives a stellar nuclear timescale of order $10^8$ years. For a numerical comparison, Table 2 of Ekström et al.'s solar-metallicity models gives core hydrogen burning lifetimes of about $88$ and $109$ million years for nonrotating and rotating $5M_{\odot }$ models, respectively.

After central hydrogen exhaustion, an inert helium core contracts and heats through Kelvin-Helmholtz contraction, while a hydrogen-burning shell supplies the growing luminosity. The envelope expands and cools: the track crosses the Hertzsprung gap toward the red giant region. Its deepening convective envelope causes first dredge-up, bringing material altered by hydrogen burning to the surface. The rapid structural crossing has a thermal timescale, roughly $10^5$–$10^6$ years for a star of this mass, rather than the full stellar nuclear timescale; an estimate follows from the Kelvin-Helmholtz cooling time $G{M}^2/(RL)$ with the evolving radius and luminosity.

The helium core reaches ignition before becoming strongly supported by electron degeneracy pressure. Thus a normal $5M_{\odot }$ star starts core helium burning without the helium flash characteristic of lower-mass stars. The Triple-alpha process produces carbon-12, and carbon-12 alpha capture produces oxygen-16; these reactions build a carbon-oxygen core. There is a convective core and continued hydrogen-burning shell activity. Core helium burning lasts order $10^7$ years: the same model table gives approximately $19$ and $18$ million years. The changing core and envelope structure can produce a blue loop in the Hertzsprung-Russell diagram; its extent depends on stellar composition, mixing and stellar rotation.

After central helium exhaustion, the carbon-oxygen core contracts and becomes increasingly supported by electron degeneracy pressure. A helium-burning shell forms above it, beneath the hydrogen-burning shell. The star returns to a cool, luminous Asymptotic giant branch configuration. On this early branch, helium-burning shell activity is an important energy source; the overlying hydrogen-burning shell can weaken during readjustment and resume later. At this mass, a deepening convective envelope can cause second dredge-up. This shell-supported evolution to the first AGB thermal pulse takes order $10^6$ years, with model-dependent duration. As the helium-burning shell becomes sufficiently thin, the thin-shell instability ends steady shell burning and marks the onset of AGB thermal pulses.

The original sketch shows the direction of stellar evolution, with higher effective temperature to the left. It illustrates a possible blue loop; its coordinates and loop shape are schematic, rather than a computed stellar model. The stages before the first AGB thermal pulse are the requested track; the later remnant evolution is described below.

In the usual final evolution, stellar winds and an AGB superwind remove the envelope before the carbon-oxygen core can ignite substantial carbon burning. Once shell fuel is exhausted, the exposed core contracts and heats, potentially ionizing the expelled gas as a planetary nebula, and finally becomes a carbon-oxygen white dwarf that loses its stored thermal energy by cooling. The relevant competition is between the AGB envelope-loss and core-growth timescales, not just the total initial mass.

Much higher late-stage mass loss strips the envelope sooner, shortens the Asymptotic giant branch phase, and allows less core growth, usually producing a smaller carbon-oxygen white dwarf. Much lower mass loss prolongs shell burning and permits greater core growth. In an extreme hypothetical case, retention of the envelope could let a degenerate carbon-oxygen core approach the Chandrasekhar mass and undergo explosive carbon burning: the proposed type 1.5 supernova channel. This is not an inevitable fate for a $5M_{\odot }$ star. Doherty et al.'s evolutionary calculations emphasize the uncertainty of that channel and find that their adopted mass loss removes carbon-oxygen envelopes before such ignition. If carbon burning instead proceeds without disruption, the possible further evolution depends on the resulting core and continued envelope retention; lowering mass loss alone does not establish an electron-capture supernova outcome.

Mechanical balance and equipotentials. Let $\Phi$ be the Newtonian gravitational potential and $\varpi$ the cylindrical distance from the rotation axis. For rigid angular velocity $\Omega$, a fluid element at rest in the rotating frame has no Coriolis acceleration. The centrifugal acceleration is ${\Omega }^2\varpi \hat{\varpi }$, so hydrostatic equilibrium becomesUsing the Poisson equation for gravity and ${\nabla }^2(x^2+y^2)=4$ gives

On any regular connected equipotential surface, every tangent displacement satisfies $d\varphi =0$ and therefore $dP=0$. Locally $P=P(\varphi )$, with $dP/d\varphi =-\rho$. This makes $\rho$ a function of $\varphi$ as well. Equivalently, taking the curl of the hydrostatic equation gives $\nabla \rho \times \nabla \varphi =0$. Thus $P,\rho$ and ${\nabla }^2\varphi$ are constant on each surface. The last conclusion follows from the modified Poisson equation; it does not constrain the magnitude of the first derivative $|\nabla \varphi |$ to be constant there.

For uniform composition, an equation of state $P=P(\rho ,T)$ determines $T$ uniquely at fixed $P,\rho$ in the usual stellar regime. For the gas-plus-radiation equation, the derivative with respect to $T$ is positive. Hence $T=T(\varphi )$, and the opacity and radiative conductivity are also functions of $\varphi$.

Surface flux and the requested sketch. Define $f(\varphi )=-\chi dT/d\varphi$, positive for the outward radiative flux. ThenOn an outer equipotential surface, $f$ is constant and the outward normal is parallel to $\nabla \varphi$. Its effective gravity is $g=|\nabla \varphi |$, so $F_n=fg$. The Stefan–Boltzmann law $F_n={\sigma }_{\mathrm{SB}}{T}_e^4$ gives the radiative gravity-darkening law,A rapidly rotating star is oblate. Centrifugal support reduces the effective gravity near the equator, which is therefore cooler and dimmer; the poles are hotter. The original sketch below uses a subcritical Roche model of a uniformly rotating star to show this geometry and the corresponding temperature variation. The material temperature on an interior equipotential is not the same quantity as the emergent effective temperature defined by surface flux.

Why purely local radiative balance fails. Differentiating the flux givesOn one equipotential, $f,f^{\prime },{\nabla }^2\varphi$ and the proposed source $\rho ϵ(\rho ,T)$ are all constant, but $|\nabla \varphi {|}^2$ generally varies with latitude. Unless $f^{\prime }=0$ or there is another special cancellation, no single value of $\rho ϵ$ can match this flux divergence at all points on the surface. This is the radiative-equilibrium obstruction in a rotating barotropic star. Slow stellar meridional circulation can transport heat to compensate; it is distinct from abandoning hydrostatic balance at leading order.

Integrated energy balance with circulation. The steady continuity equation is $\nabla \cdot (\rho \mathbf{v})=0$. Combining the entropy–enthalpy relation with hydrostatic balance givesIntegrate the energy equation over $V$ and apply the divergence theorem:Uniform composition makes $h=h(\rho ,T)$ constant on the enclosing equipotential, and $\varphi$ is constant there by definition. Pull $(h+\varphi )$ out of the last integral; the remaining net mass flux is zero by steady continuity. Thus the equipotential luminosity conservation with stellar circulation isThis does not require $\mathbf{v}\cdot \hat{n}=0$ pointwise on $S$: inward and outward circulation may cancel in the net mass flux.

Finally, $f$ is constant on $S$ and the divergence theorem applied to the effective potential givesTherefore $-\chi dT/d\varphi =L/(4\pi Gm-2{\Omega }^{2}V)$. Since $dP/d\varphi =-\rho$,The mass $m$ is the mass enclosed by this surface, not necessarily the total stellar mass. The denominator is positive for the regular hydrostatic surfaces used here; near loss of centrifugal confinement the assumed stellar model must be reconsidered. The derivation uses slow circulation, retaining hydrostatic/equipotential thermodynamics to leading order while allowing a nonzero advective heat flux.

Use stellar homology: dimensionless density, temperature and luminosity profiles are taken to have the same shape, so the dimensional scale relations follow from the structure equations. Mass conservation and hydrostatic support giveThe ideal gas equation then gives $T_c\propto \mu M/R$, with the dimensional constants held fixed. The specific stellar energy-generation rate has the prescribed hydrogen factor and temperature exponent, soIndependently, radiative diffusion in a star gives $T_c/R\propto \kappa {\rho }_{c}L/(R^2{T}_c^3)$, henceEquating the nuclear supply and radiative transport establishes composition-dependent CNO homology:At one fixed composition, these reduce to $R\propto M^{3/4}$ and $L\propto M^3$ as required.

The Stefan–Boltzmann law gives $T_e^4\propto L/R^2\propto M^{3/2}$, so $T_e\propto M^{3/8}$. Thus the logarithmic Hertzsprung-Russell diagram slope along this zero-age main sequence isIncreasing mass moves upward and leftward on the conventional diagram, whose temperature axis is reversed.

For composition changes, neglect the small heavy-element contribution to the mean molecular weight and set $Y\simeq 1-X$. Then ${\mu }^{-1}=(3+5X)/4$. The fixed trace catalyst abundance is absorbed into ${ϵ}_0$; this is not an assertion that an exactly metal-free star can undergo CNO burning. Substitution into the homology laws givesFor a retained, nonnegligible fixed metallicity $Z$, the corresponding mean-molecular-weight factor is $3+5X-Z$ instead. The displayed requested formula uses the trace-metal approximation.

At fixed mass, differentiating with respect to $X$ at the formal reference $X=1$ givesConsequentlyand the initial fully mixed hydrogen-burning evolutionary track has slopeHydrogen depletion has $dX<0$, so both luminosity and effective temperature initially increase. The star therefore moves upward and leftward, on a shallower logarithmic track than the mass sequence. The figure uses dimensionless luminosity and temperature relative to its initial reference star; the track is computed for declining $X$ at constant mass.

Ordinary cluster main-sequence turn-offs instead evolve toward brighter, cooler, more extended configurations and eventually the giant region. Normal stars are therefore not mixed throughout their entire mass during hydrogen burning: the core composition changes while the envelope retains hydrogen, and composition and entropy gradients develop. Core exhaustion and subsequent shell burning then produce very different structural evolution. Complete mixing is a special evolutionary scenario, not the standard explanation of ordinary giant formation.

These hydrogen-burning scalings must not be extrapolated literally to $X=0$. Their vanishing nuclear coefficient would make the inferred radius tend to zero and the effective temperature singular; in reality hydrogen exhaustion requires a different structural and energy-generation model. A homogeneous remnant may be helium-rich, but its later structure is not described by the same CNO law. Also, the cover-sheet radiation-pressure expression has a temperature-power typo: the physical term is $a{T}^4/3$. Radiation pressure was neglected in this derivation, so that typo does not enter the calculation.

In the radiative layer, take the enclosed mass to be $M_1$, the transported luminosity to be $L$, and the mean molecular weight to be constant. The envelope's own gravity is neglected. Divide the hydrostatic and radiative temperature equations, and then eliminate density using the gas equation of state:Put $p=n+1$ and $b=4+m+n$. For $p,b\ne 0$, integration givesDefining $C^p=pB/b$ and $K=K_P/C^p$ yieldsIn the relevant positive-exponent case, writing $K=T_0^b$ produces the stated form whenever this integration constant is nonnegative. The exceptional values excluded above lead to logarithms: $p=0$ gives $\log P=B{T}^b/b+\mathrm{const}$, $b=0$ gives $P^p=pB\log T+\mathrm{const}$, and both zero give $\log P=B\log T+\mathrm{const}$.

For radiative-envelope convection matching, compute the logarithmic gradient of the pressure relation:At the marginal boundary, the Schwarzschild criterion equates this with the ideal-gas adiabatic gradient ${\nabla }_{\mathrm{ad}}=(\gamma -1)/\gamma$. HenceThe positive real $T_0$ notation requires a nonnegative bracket; a general signed $K$ is the safer integration constant otherwise. For the later case $n=1,m=3$ and a monatomic gas, this bracket is $3/5$, so $T_0$ is of the same order as $T_b$.

Deep below the convective boundary, with the temperature-power term dominating $K$, the pressure law becomes $P\simeq C{T}^{b/p}$. Differentiating this approximation and using hydrostatic balance givesThe general temperature is a term proportional to $1/r$ plus a boundary-dependent constant. Well inside the boundary radius, where that term dominates, it reduces toThe inverse-radius profile is an inner approximation, not an assertion that the boundary term vanishes everywhere.

Now take $n=1,m=3$, so $p=2$, $b=8$. Then $P\simeq C{T}^4$, $\rho \simeq (\mu C/\mathcal{R})T^3$, and define $A_T=\mu G{M}_1/(4\mathcal{R})$ so that $T\simeq A_T/r$. The luminosity supplied by the thin burning layer isBecause the integral is strongly concentrated near $R_1$ and $r_b\gg R_1$, the outer term is negligible. ThusThe factor $\rho ϵ$, rather than merely $ϵ$, is essential since the specified energy-generation rate is per unit mass. This integration uses the deep radiative profile as a leading thin-shell approximation.

Crucially, $C$ is not independent of luminosity. Its defining relation givesSubstitution into the shell integral therefore yields $L^2\propto M_1^{17}{R}_1^{-13}$. The degenerate-core relation gives $R_1\propto M_1^{-1/3}$, soThis degenerate-core shell-burning mass-luminosity relation is specific to the adopted opacity, burning exponent and structural approximations.

Removing some envelope mass while retaining the hydrogen-burning shell and the deep-envelope regime leaves the leading luminosity approximately unchanged at fixed core mass. The envelope can instead change the radius and effective temperature. This independence is not valid for arbitrarily complete stripping: removing the fuel or the pressure-confined shell structure quenches or modifies shell burning, and the assumed giant model then fails.

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.

Let $\Psi$ be the gravitational potential, with ${\nabla }^2\Psi =4\pi G\rho$, and let ${\varpi }^2=x^2+y^2$ measure distance from the rotation axis. Uniform solid-body rotation contributes the centrifugal potential $-{\Omega }^2{\varpi }^2/2$. In the corotating frame, static force balance therefore gives the uniformly rotating stellar hydrostatic equilibrium equationsThe final coefficient follows because ${\nabla }^2(x^2+y^2)=4$.

A tangent displacement along a regular connected equipotential surface has zero dot product with $\nabla \varphi$, and hence with $\nabla P$. Thus $P=P(\varphi )$ locally. Taking the curl of the hydrostatic equilibrium equation gives $\nabla \rho \times \nabla \varphi =0$, so $\rho =\rho (\varphi )$ on such surfaces as well. The Poisson equation then shows that ${\nabla }^2\varphi$ is constant on each equipotential surface.

The magnitude $g=|\nabla \varphi |$ need not be constant. For instance $\varphi =A(x^2+y^2)+B{z}^2$, with positive $A\ne B$, has constant Laplacian but different pole and equator gradient magnitudes on an ellipsoidal equipotential. In a rotating star, flattening and centrifugal acceleration likewise make the effective gravity latitude dependent.

With uniform composition, an equation of state determines $T$ from $P$ and $\rho$; for the supplied ideal gas plus radiation pressure law its right side increases monotonically with $T$. Hence $T=T(\varphi )$, and the radiative conductivity $\chi (\rho ,T)$ is also constant on an equipotential surface. Put $f(\varphi )=-\chi dT/d\varphi$. Then $\mathbf{F}=f(\varphi )\nabla \varphi$. On the stellar surface, below mass shedding, the outward normal is $\nabla \varphi /g$, soThis is the radiative gravity-darkening law. The effective temperature measures emergent flux; it is not the same quantity as the interior material temperature that was constant on an equipotential.

A rapidly rotating star is oblate, with reduced effective gravity and cooler effective temperature around the equator. The poles are hottest. The original sketch uses a Roche model of a uniformly rotating star only to illustrate the shape and surface-flux variation:

The obstruction to purely local radiative equilibrium is visible by taking the divergence:On an equipotential surface, $\rho ϵ$, $f$, $f^{\prime }$ and the second term are constant, but $g^2$ generally varies. Unless $f^{\prime }=0$ or another special circumstance removes the mismatch, $\nabla \cdot \mathbf{F}=\rho ϵ$ cannot hold at every latitude. This is the radiative-equilibrium obstruction in a rotating barotropic star; heat advection by stellar meridional circulation resolves the missing energy transport in the leading hydrostatic treatment.

For the steady circulation, mass conservation gives $\nabla \cdot (\rho \mathbf{v})=0$. The specific enthalpy identity and hydrostatic equilibrium giveConsequentlyBoth $h$ and $\varphi$ are constant on the bounding equipotential surface. By the divergence theorem and steady continuity equation,Integrating the given thermal energy equation therefore gives equipotential luminosity conservation with stellar circulation:Individual circulation streams can cross the surface; only their net steady mass flux has to vanish.

Since $f$ is constant on that surface, integrating the Poisson equation givesThus $dT/d\varphi =-L/[\chi (4\pi Gm-2{\Omega }^{2}V)]$, while $dP/d\varphi =-\rho$. Dividing these derivatives and inserting $\chi =4ac{T}^3/(3\kappa \rho )$ yields the radiative gradient in a uniformly rotating star:The bracket must be positive for the regular outward-gravity configuration being used.

For an interior equipotential, both mass factors must be the enclosed mass $m$. The printed final bracket uses $M$, but its derivation gives $m$; at the outer surface $m=M$ and the printed bracket is recovered. For a direct check, a uniform-density model has $m=\rho V$, so the true correction ${\Omega }^{2}V/(2\pi Gm)={\Omega }^2/(2\pi G\rho )$ is independent of the chosen enclosed volume. Replacing $m$ by a fixed total mass would incorrectly make it volume dependent.