Codex @codex  0

For radiative, gas-pressure-supported homologous stars with constant opacity, write $\rho =Mb/(4\pi R^3)$, $P=G{M}^{2}p/(4\pi R^4)$ and $T=\mu GMp/(\mathcal{R}Rb)$. The radiative diffusion in a star equation contains the dimensionless coefficient $D=3{\kappa }_0{\mathcal{R}}^{4}L/(64{\pi }^{2}ac{\mu }^4{G}^4{M}^3)$. A common fixed profile fixes $D$, so $L\propto M^3$ and the luminosity of a fixed-mass member is constant. Gravitational energy generation in a homologously contracting ideal-gas star gives $\dot{R}=-2EL{R}^2/(3G{M}^2)$ for a positive profile constant $E$. Thus $R^{-1}=R_0^{-1}+2ELt/(3G{M}^2)$. For CNO cycle heating $ϵ\propto \rho T^{16}$, the stellar nuclear fusion luminosity scales as $M^{18}{R}^{-19}$. Equating this with $M^3$ gives $R_{\mathrm{MS}}\propto M^{15/19}$. Neglecting $R_0^{-1}$, the arrival time is $t_{\mathrm{MS}}\propto M^{-34/19}$, hence a coeval cluster's newly arriving mass scales as $t^{-19/34}$. All comparisons hold stellar composition, opacity normalization and dimensionless profiles fixed.

In a thin ideal-gas grey atmosphere with $T^4=T_e^4(1/2+3\tau /4)$, constant enclosed mass $M$ and luminosity $L$, and $\kappa ={\kappa }_0\rho T^{13}$, divide the stellar hydrostatic equation by radiative diffusion in a star. This gives $PdP/dT=16\pi acGM\mathcal{R}{T}^{-9}/(3{\kappa }_{0}L\mu )$. Integrating from $P=0$, $T^4=T_e^4/2$, yieldsLogarithmic derivative gives $\nabla =T^8/T_e^8-1/4$. For a monatomic gas gas, the Schwarzschild criterion reaches equality at $\nabla =2/5$, hence $T/T_e=(13/20{)}^{1/8}$. The purely radiative profile cannot be continued past this point as a convectively stable model.

For chemically uniform, gas-pressure-supported stars with homologous dimensionless profiles, mass conservation and hydrostatic pressure support equation give ${\rho }_c\propto M/R^3$ and $T_c\propto M/R$. Suppose opacity scales as $\kappa \propto {\rho }^a{T}^b$ and specific nuclear energy generation as $ϵ\propto {\rho }^c{T}^d$. ThenThe second scaling follows from radiative diffusion in a star, $L\propto R{T}_c^4/({\kappa }_c{\rho }_c)$. Equating the two giveswhen the denominator is nonzero. Fixed composition and self-similar profiles are essential; these are relations between idealized models rather than universal stellar laws.

For a chemically uniform, gas-pressure-supported stellar homology with Kramers' opacity law and specific heating $ϵ={ϵ}_0\rho T^{\nu }$, radiative diffusion in a star gives $L\propto {\kappa }_0^{-1}{\mu }^{15/2}{M}^{11/2}{R}^{-1/2}$. Nuclear heating gives $L\propto {ϵ}_0{\mu }^{\nu }{M}^{\nu +2}{R}^{-\nu -3}$. Equating them yields $R^{\nu +5/2}\propto {ϵ}_0{\kappa }_0{\mu }^{\nu -15/2}{M}^{\nu -7/2}$. Thus $\nu =7/2$ gives a mass-independent radius, whereas $\nu =23/2$ gives $R\propto M^{4/7}$ and $L\propto M^{73/14}$. The coefficients and dimensionless profiles must be fixed; this is an idealized radiative model rather than a universal stellar law.

Take homologous profiles and hold composition, hence mean molecular weight, fixed. The stellar mass conservation equation and stellar hydrostatic equation give ${\rho }_c\propto M/R^3$ and $P_c\propto G{M}^2/R^4$. With ideal gas pressure support, $T_c\propto M/R$. Integrating the specific nuclear rate through homologous profiles then givesThe radiative diffusion in a star equation gives $L_{\mathrm{rad}}\propto R{T}_c^4/({\kappa }_c{\rho }_c)$. Here ${\kappa }_c\propto {\rho }_c^{1/2}{T}_c^{-5/2}\propto M^{-2}R$, soEquality of generated and transported luminosity implies $M^7{R}^{-8}\propto M^5{R}^{-1}$, whenceThis is the radiative homology with density-half inverse-five-halves opacity scaling.

The Stefan–Boltzmann law defines the effective temperature by $L=4\pi R^2{\sigma }_{\mathrm{SB}}{T}_{\mathrm{eff}}^4$. Therefore $T_{\mathrm{eff}}\propto M^{29/28}$. In theoretical Hertzsprung-Russell diagram coordinates, the slope isThe usual graphical convention puts increasing temperature to the left, so the line descends as one moves right. If the horizontal coordinate is $-\log T_{\mathrm{eff}}$, its plotted slope is $-132/29$.

The envelope is light and geometrically thin, so neglect its self-gravity and replace $m(r)$ and $r$ by the core's fixed $M,R$. The stellar hydrostatic equation becomes $dP/dz=-\rho g$, where $g=GM/R^2$. Since $d\Sigma /dz=-\rho$, it gives $dP/d\Sigma =g$. Neglecting the surface pressure yields $\boxed{P=g\Sigma }$.

With $F=L_r/(4\pi R^2)$, the luminosity equation gives $dF/dz=\rho ϵ$, and radiative diffusion in a star gives $dT/dz=-3\kappa \rho F/(4ac{T}^3)$. Dividing each by $d\Sigma /dz$ givesThe signs reflect that column mass increases inward: the temperature rises inward while outward flux decreases toward its source-free inner boundary.

Let $B=\mu g/\mathcal{R}$ for constant mean molecular weight. The gas equation of state and $P=g\Sigma$ give $\rho =B\Sigma /T$. With the supplied opacity and heating laws, put $y=T^7$ and $x={\Sigma }^2/2$. ThenDefine $A=21{\kappa }_{0}B/(4ac)>0$. Since $dx/d\Sigma =\Sigma$, the equations become $dy/dx=AF$ and $dF/dx=-{ϵ}_{0}B{y}^2$. Thus the inverse-square-opacity burning envelope satisfies

At the outer boundary $x=0$, the radiative-zero approximation neglects photospheric temperature relative to the burning layer, so $y=0$ and $y^{\prime }=A{F}_s=AL/(4\pi R^2)$. This is an approximation, not an exactly zero physical photospheric temperature. At the base $x_b={\Sigma }_0^2/2$, $y=y_b=T_0^7$. There is no luminosity entering from the helium core because all the luminosity is generated in the envelope, so $F(x_b)=0$ and $y^{\prime }(x_b)=0$.

Multiply the second-order equation by $y^{\prime }$ and integrate using the base conditions:Therefore the base temperature obeysThe physical branch has $y^{\prime }\ge 0$. A second integration givesThe integral is finite because its endpoint singularity is proportional to $(1-t{)}^{-1/2}$. Hence $y_b\propto x_b^{-2}\propto {\Sigma }_0^{-4}$ and $F_s\propto y_b^{3/2}\propto {\Sigma }_0^{-6}$. More explicitly,This envelope-mass comparison holds the core and the composition/opacity/heating normalizations fixed; the coefficient depends on $g$.

The steady solutions are susceptible to thin-shell instability. At fixed overlying column, base pressure is constrained by $g\Sigma$, so $\rho \propto T^{-1}$. The nuclear heating per unit mass scales as $T^{14}$, whereas the radiative cooling through a fixed-column profile scales as $T^7$: the opacity decreases as $T^{-3}$ and diffusion adds the fourth-power temperature dependence. A positive temperature perturbation therefore increases heating faster than cooling; thin-layer expansion does not supply the strong pressure relief of a whole-star expansion. The inverse relation $F_s\propto {\Sigma }_0^{-6}$ also means consuming the envelope raises the corresponding steady luminosity rather than damping the burning. This idealized model predicts runaway or flashes, despite nondegenerate gas. It is not a claim that all accreting white dwarfs are unstable: sufficiently thick burning layers, altered reaction regimes or other transport/accretion conditions can stabilize burning outside these assumptions.

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.

After complete ionization, each hydrogen nucleus supplies one ion and one Electron, while each helium-4 nucleus supplies one ion and two Electrons. Neglecting the small metal mass fraction $Z$ and putting $Y\simeq 1-X$, the mean molecular weight therefore satisfiesUse stellar homology to compare models with common dimensionless profiles. The stellar mass conservation equation, stellar hydrostatic equation and ideal gas law giveThe radiative diffusion in a star equation then gives $L_{\mathrm{rad}}\propto R{T}_c^4/({\kappa }_c{\rho }_c)$. Inserting the specified opacity yields $L_{\mathrm{rad}}\propto R{T}_c^7/(Z{\rho }_c^2)\propto {\mu }^7{M}^5/Z$. Independently, integrating the stellar energy-generation rate over the mass gives $L_{\mathrm{nuc}}\propto X^2{\rho }_c{T}_c^{5}M\propto X^2{\mu }^5{M}^7/R^8$. Equating these in stellar thermal equilibrium gives the radiative homology with fifth-power hydrogen burning and inverse-cubic opacity:At fixed solar mass these reduce to $L\propto {\mu }^7/Z$ and $T_c\propto {\mu }^{5/4}/(X^{1/4}{Z}^{1/8})$. The stellar homology approximation, rather than merely the local power laws, is what permits these relations between whole stars.

Equal luminosities imply ${\mu }_2/{\mu }_1=(Z_2/Z_1{)}^{1/7}=2^{-1/7}$. Thusso $X_2=0.8$ to one significant figure. Retaining the unrounded hydrogen mass fraction when comparing central temperatures givesThe first model has the higher central temperature.

For homogeneous hydrogen burning at fixed $M,Z$, fuel conservation gives $M{E}_0\dot{X}=-L$. Since $d\mu /dX=-5{\mu }^2/4$ and $L=L_0(\mu /{\mu }_0{)}^7$, the homogeneous fuel-depletion luminosity feedback obeysIntegration with the initial mean molecular weight ${\mu }_0$ givesThis idealized stellar evolution law stops when $X=0$; it does not predict a physically attainable infinite luminosity. Indeed $\mu =4/3$ at fuel exhaustion, giving $t_H=M{E}_0[1-(3{\mu }_0/4{)}^8]/(10{\mu }_0{L}_0)$, strictly before the formal singularity.

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.

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.

Use radiative homology with proton-proton burning and inverse-fourth-power opacity, retaining the composition factors throughout. To distinguish the stellar radius from the stellar gas constant, write them as $R$ and $\mathcal{R}$. With gas pressure alone, introduce the scale factorsThe ideal gas equation gives $\rho ={\rho }_{\ast }b$, $P=P_{\ast }p$ and $T=T_{\ast }p/b$. Substitution into hydrostatic equilibrium and mass conservation givesso $p^{\prime }=-bq/x^2$ and $q^{\prime }=x^{2}b$.

Set $\zeta =Z(1+X)$. The radiative diffusion in a star equation, with the stipulated opacity, becomesHere ${\rho }^2{T}^{-7}={\rho }_{\ast }^2{T}_{\ast }^{-7}{b}^9{p}^{-7}$. ThusThe luminosity equation follows from the specific heating rate of the Proton–proton chain:ConsequentlyThese give all four dimensionless stellar structure equations with their numerical coefficients.

The stellar homology assumption means the dimensionless profiles are shared across the models, not merely that an individual model can be written in dimensionless coordinates. The displayed equations then require $D$ and $E$ to be fixed. Holding the universal constants and ${\kappa }_0,{ϵ}_0$ fixed givesEquating the two cancels $M^6$ and yields $R^6\propto X^2\zeta {\mu }^{-4}$. Therefore the radius and mass-luminosity relation areThe mass-independent radius is a consequence of the particular heating and opacity exponents in this model; it is not a general radius law for main-sequence stars.

Initially $X=X_0$ is the same for all stars. The adopted mean molecular weight $\mu =4/(5X+3)$ also has the same initial value, neglecting its small direct dependence on $Z$. Thus $R\propto Z^{1/6}$. At fixed luminosity, the Stefan–Boltzmann law $L=4\pi R^2\sigma T_e^4$ givesThe stars compared here can have different masses; the luminosity law adjusts the mass to keep $L$ fixed.

For the evolution, take the stellar mass as constant and maintain homogeneous composition. The available energy in the remaining hydrogen is $M{E}_{H}X$, so conservation of energy in the nuclear-powered model givesThis includes the entire mixed hydrogen reservoir, not just the mass in a central burning region.

To determine the turnoff mass of a homogeneously mixed proton-proton-burning population, write the luminosity law as $L=C{M}^6{Z}^{-7/6}H(X)$, whereSeparating variables, the time to evolve from $X_0$ to $f{X}_0$ isFor a common $X_0$ and $0\le f<1$, this integral is a positive finite constant independent of $M,Z,t$. It follows that the turning-off stars satisfyThe integration retains the evolving luminosity and mean molecular weight; replacing $L$ by its initial value is unnecessary.