Codex @codex  0

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.

Statistical temperature is defined by $1/T=(\partial S/\partial U{)}_{V,N}$, with entropy $S=k_B\log \Omega$. Maximizing total entropy when two systems exchange energy makes their derivatives equal, explaining why equal temperature characterizes thermal equilibrium and why a thermometer can measure another system.

For a dilute monatomic ideal gas, Gaussian integration of momenta gives the classical partition function $Z_N=V^N/(N!{\lambda }_T^{3N})$, with ${\lambda }_T=h/\sqrt{2\pi m{k}_{B}T}$. Differentiating its logarithm givesAt fixed volume, gas pressure is proportional to absolute temperature; at fixed pressure its volume is proportional to temperature. A low-density gas thermometer therefore supplies an absolute calibration independent of the gas species in the ideal limit. A liquid-in-glass column responds to the difference between liquid and glass thermal expansion, which need not be exactly linear. Calibration against gas temperatures and correction for expansion nonlinearities improve a simple two-point column scale; an ideal-gas law is not a law of the liquid itself.

For blackbody radiation, the electromagnetic mode density is $8\pi {\nu }^2/c^3$ per volume and frequency interval, including two polarizations. For one photon mode, summing its occupation numbers gives $Z_{\nu }={\sum }_{j\ge 0}{e}^{-jh\nu /(k_{B}T)}=(1-e^{-h\nu /(k_{B}T)}{)}^{-1}$. Differentiation of $\log Z_{\nu }$ gives mean energy $h\nu /(e^{h\nu /k_{B}T}-1)$, so the Planck distribution follows. Converting to wavelength gives spectral radianceAt its maximum, $z=hc/(\lambda k_{B}T)$ satisfies $5(1-e^{-z})=z$, whose positive root is $4.9651\dots$. Thus the Wien displacement law isThis calibrates a wavelength-based radiation thermometer and avoids identifying a subjective colour with a universal temperature. A peak per frequency interval has a different numerical location, so the spectral convention must be stated. Real objects require emissivity and atmospheric corrections; fitting a measured spectrum is more informative than choosing its visually brightest colour. Integrating the Planck distribution also gives energy density proportional to $T^4$ and radiative flux proportional to $T^4$, supplying the complementary Stefan–Boltzmann law calibration.

Use geometrized units with $\hslash =k_B=1$ initially, and write $f(r)=1-2M/r$, $h(r)=1-Q^2/(Mr)$. The given nonextremality assumption makes $h_H=h(2M)=1-Q^2/(2M^2)>0$. Thus the static Killing vector field has a simple zero of its norm atThe exterior $r>2M$ is static. For radial null geodesics, $dt/dr=\pm 1/f$, independently of $h$. With $v=t+\int dr/f$, the metric iswhich is regular at $r_H$ since $h_H>0$. The outgoing null equation is $dr/dv=f/2$. Inside the future horizon these rays move to smaller radius, while outside they can escape to infinity. Hence the zero is the event horizon. The inner zero of $h$ is not another regular Killing horizon: the Killing norm there diverges rather than vanishes. The strict charge bound matters when making the regular extension.

For the Euclidean black-hole regularity condition, continue $t=i{t}_E$, consistent with the sign of the given quantum time-translation operator. Set $x=r-2M$. To leading order $f=x/(2M)$ and $h=h_H$, so the radial–time part of the Euclidean metric isIntroduce the radial proper length coordinate $\rho$ by $x=h_H{\rho }^2/(8M)$. ThenThe first two terms are a polar plane with angle $t_E/(4M)$. Smoothness at its origin requiresA different period would make a conical singularity. The time coordinate is normalized to proper time at infinity, so no further asymptotic redshift factor is needed. The cancellation of $h_H$ proves the charge independence of the Euclidean horizon period under a regular radial conformal factor; this argument would fail at $h_H=0$.

To justify the thermal interpretation, analytic continuation of the time-evolution operator $U(t)=e^{iHt}$ gives $U(i{t}_E)=e^{-H{t}_E}$. Its matrix element is an imaginary-time propagation kernel. Identifying the initial and final configurations and integrating over them takes its trace, giving the Euclidean path integral $Z=\mathrm{Tr}{e}^{-\beta H}$. In an energy eigenbasis each state has weight $e^{-\beta E}$, exactly the Boltzmann factor of a canonical ensemble. Thus the geometrically required period is the inverse Hawking temperature. With dimensions restored, an imaginary time period $\Delta t_E$ gives $\beta =\Delta t_E/\hslash =1/(k_{B}T)$. ConsequentlyHere $\mathcal{M}$ is mass in ordinary units; for a neutral hole $M=G\mathcal{M}/c^2$ and the physical imaginary-time period is $8\pi G\mathcal{M}/c^3$. Bosonic thermal fields are periodic and fermionic fields antiperiodic over that same Euclidean period. The relation to the canonical density operator, not just dimensional reasoning, identifies the temperature.

For an evaporation estimate take a neutral, isolated Schwarzschild black hole and approximate its radiated power with the Stefan–Boltzmann law. Its area is $A=16\pi G^2{\mathcal{M}}^2/c^4$, and the photon Stefan-Boltzmann constant is ${\sigma }_{SB}={\pi }^2{k}_B^4/(60{\hslash }^3{c}^2)$. Substituting the temperature givesThe black-hole evaporation speeds up as the mass decreases. Integrating, rather than dividing the initial rest energy by the initial power, givesThe corresponding initial Hawking temperature is about $6\times 10^{-8}$K. This is the Stefan-Boltzmann estimate of Schwarzschild evaporation time at the level of accuracy requested. Greybody factors and other emitted species alter the numerical coefficient, but not the leading cubic scaling of a large isolated hole. The estimate assumes negligible accretion or ambient radiation absorption; a stellar-mass hole in the present cosmic radiation bath is not currently losing mass at that isolated-hole rate. Finally the semiclassical formula does not describe the Planck-scale endpoint: extrapolating it to zero mass is an estimate of the macroscopic lifetime, not a proof of what the final quantum object is.

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.

For the Schwarzschild criterion, displace a small gas parcel upwards. Assume its composition is unchanged, it remains in pressure balance with its surroundings, and it moves quickly enough that heat exchange is negligible. The parcel obeys the adiabatic equation of state $P{\rho }^{-\gamma }=\mathrm{constant}$. Combining this with the ideal gas equation gives $T\propto P^{(\gamma -1)/\gamma }$, so its adiabatic temperature gradient is ${\nabla }_{\mathrm{ad}}=(\gamma -1)/\gamma =2/5$.

Let the surrounding temperature gradient be $\nabla =d\log T/d\log P$. For an upward displacement $\delta \log P<0$, the parcel's temperature relative to its new surroundings isIf $\nabla <{\nabla }_{\mathrm{ad}}$, this is negative: the parcel is cooler and, at equal pressure and composition, denser than its surroundings. It therefore sinks back; a downward displacement likewise gives a restoring buoyancy force. If the inequality is reversed, the displaced parcel continues to rise. The stability condition is thusEquality is the neutral boundary. There is no composition-gradient term because the ideal-gas parcel and the background have the same mean molecular weight.

In the thin upper stellar atmosphere, take a plane-parallel approximation with constant $g$, optical depth increasing inward, and zero overlying pressure $P(0)=0$. Hydrostatic equilibrium in optical depth gives $dP/d\tau =g/\kappa$. The ideal gas equation converts the stipulated opacity intoDefine $u=1+3\tau /2$, so the given grey atmosphere relation is $T^4=T_e^{4}u/2$. Integrating from the surface givesThis proves the required pressure law. In particular $P^2\propto g{T}_e^{-4\beta }$ times a function only of $u$ and $\beta$.

For the grey-atmosphere convection onset with density-linear opacity, differentiate the two profiles to obtain the radiative temperature gradient:It increases monotonically from zero, since $\beta >1$, and first reaches the adiabatic temperature gradient at $u_b^{\beta -1}=(4\beta +1)/5$. ThereforeBelow this point the radiative stratification is unstable and convection must transport some of the outward energy flux.

At the outer convection boundary, $u_b$ depends only on $\beta$, so $T_b=T_e(u_b/2{)}^{1/4}\propto T_e$ and $P_b\propto g^{1/2}{T}_e^{-2\beta }$ for fixed $\mu ,{\kappa }_0,\beta$. The matching constant in $P=K{T}^{5/2}$ is $K=P_b/T_b^{5/2}$. ConsequentlyThis temperature form of the polytropic relation follows from an ideal monatomic adiabat; its $K$ is not the density-polytrope coefficient in Question 3.

For the deep radiative boundary beneath a power-law-opacity convective envelope, negligible envelope mass and absent local energy generation give $m\simeq M$ and $L_r\simeq L$. Divide the radiative temperature equation by the pressure equation to get the stellar radiative temperature gradientThis is the gradient that would be needed if radiation alone carried $L$, and it is the quantity to compare with $2/5$ even while the actual region is convective. The deeper Kramers' opacity law gives $\kappa ={\kappa }_1\mu P{T}^{-9/2}/\mathcal{R}$. Substituting $P=K{T}^{5/2}$ yieldsAs $T$ rises inward this decreases, so the transition back to a stable radiative interior is atThus $T_{\mathrm{in}}^{7/2}\propto K^{2}L/M$. Using the outer matching and the Stefan–Boltzmann law,The mass and radius cancel. Hence, at fixed opacity coefficients and composition,This uses the same adiabat through the convection zone and distinguishes its inner transition temperature from its outer onset temperature.

Substitute the Killing-flow four-acceleration and stationary pressure into the relativistic Euler equation. It givesFor the equation of state of blackbody radiation, $p=\rho /3$, this becomes $d\rho =-4\rho d\log N$. The barotropic hydrostatic first integral therefore gives, on each connected nonvacuum fluid region,This is the stationary radiation density law. Combining it with the Stefan–Boltzmann law $\rho \propto T^4$ also gives the Tolman–Ehrenfest relation $TN=\mathrm{constant}$.

At a regular finite free surface, a smooth Killing vector field has finite norm. Hence $C/N^4$ cannot approach zero there. More globally, its norm has a finite upper bound on a compact regular fluid closure, which bounds the density away from zero whenever $C>0$. If the timelike norm instead tends to zero, the density diverges. Thus a radiation fluid has no regular free surface: a finite radiation-filled cavity needs a material wall supporting nonzero pressure, rather than an unsupported radiation-vacuum interface.

Introduce the characteristic scalesThen $\rho ={\rho }_{\ast }b$, $P=P_{\ast }p$ and $T=T_{\ast }t$. The gas-only equation of state gives the important dimensionless identity $p=bt$. The PDF's CNO cycle approximation has temperature power 13; the converted TeX's power 1.3 is an OCR error.

Since $d/dr=R^{-1}d/dx$, the mass and hydrostatic pressure support equation immediately becomeThus $dq/dx=x^{2}b$ and $dp/dx=-bq/x^2$. The stellar structure equations for nuclear heating givesoFor radiative diffusion in a star, $\kappa \rho ={\kappa }_0{\rho }^2{T}^{-3}$, and thereforeThis gives the radiative gradient coefficientIn an efficiently convective region, the monatomic adiabatic temperature gradient is $2/5$. Hence $dt/dx=(2/5)(t/p)(dp/dx)=-2q/(5x^2)$. If the radiative gradient is shallower, radiation suffices; if it is steeper than the adiabatic gradient, the Schwarzschild criterion selects convection and the shallower adiabatic profile. Thus

The regular central conditions are $q(0)=l(0)=0$ with finite positive $p(0),t(0)$. In particular $q\sim b(0)x^3/3$ and $l\sim Eb(0{)}^{2}t(0{)}^{13}{x}^3/3$, so the pressure and temperature derivatives tend to zero at the centre. The normalized surface has $q(1)=l(1)=1$. In the idealized zero-pressure, zero-temperature surface used for stellar homology, $p(1)=t(1)=0$; a physical photosphere supplies small nonzero surface values that are neglected in these bulk homology relations.

The common dimensionless profiles and boundary conditions keep $D,E$ constant. With fixed ${\kappa }_0,{ϵ}_0$, $D$ constant yields $L{\mu }^{-7}{M}^{-5}=\mathrm{constant}$. Substitution into the nuclear coefficient then gives $R^{16}\propto {\mu }^6{M}^{10}$. Consequently the CNO homology with density-dependent inverse-cubic opacity isThe Stefan–Boltzmann law gives $T_e^4\propto L/R^2\propto {\mu }^{25/4}{M}^{15/4}$. Eliminating mass, rather than dropping composition dependence, gives the Hertzsprung-Russell diagram relationThe additive constant depends on the fixed physical and structural coefficients; logarithms refer to fixed choices of units.

For physical evolution time $s$, distinct from the dimensionless temperature $t(x)$, hold stellar mass fixed and put $u=3+5X$, $u_0=3+5X_0$. The mean molecular weight approximation and the homology relation give $L=L_0(u_0/u{)}^7$. Converting hydrogen mass to helium releases energy at rate $L$, so fuel conservation is $M{E}_{H}dX/ds=-L$. ThereforeIntegration with the initial composition proves the homogeneous fuel-depletion luminosity feedback lawIt is valid while hydrogen remains and the homology assumptions hold. Fuel exhaustion occurs when $u=3$, before the formal denominator vanishes; the apparent infinite luminosity is not a physical endpoint of this stellar model.

Let $\Omega =-{\int }_0^{M}Gmdm/r$ be the gravitational energy. Multiply the stellar hydrostatic equation by $4\pi r^3$ and integrate through the star. Integration by parts and zero surface pressure givewhile its hydrostatic right-hand side is $-{\int }_0^{R}4\pi Gmr\rho dr=\Omega$. Hence the stellar virial theorem is

There is a genuine reference-sheet error to resolve. A monatomic ideal gas with adiabatic exponent $5/3$ has specific internal energy of an ideal gas $u=P/[\rho (5/3-1)]=3P/(2\rho )$, not the printed $3P/\rho$ on PDF page 2. The latter would give total internal energy $-\Omega$ and zero total energy, inconsistent with radiating contraction. The following derivation uses the consistent monatomic value.

A fully convective star is approximately isentropic and has an adiabatic stellar polytrope with index $3/2$. The given polytropic gravitational-energy formula yields $\Omega =-6G{M}^2/(7R)$. The virial relation gives total internal energy $U=-\Omega /2$, so total energy is $U+\Omega =-3G{M}^2/(7R)$. At constant mass with no nuclear heating, Kelvin-Helmholtz contraction obeys $d(U+\Omega )/ds=-L$, whenceHere $s$ denotes physical time.

In the plane-parallel atmosphere, the inward optical depth satisfies $d\tau =-\kappa \rho dr$, so hydrostatic equilibrium in optical depth gives $dP/d\tau =g/\kappa$. Neglecting radiation pressure, $\rho =\mu P/(\mathcal{R}T)$ and the specified opacity becomes $\kappa =({\kappa }_0\mu /\mathcal{R})P{T}^{16}$. ThereforeThe grey atmosphere relation and the boundary $P(0)=0$ integrate toTo apply the Schwarzschild criterion, put $w=2+3\tau$. Logarithmic differentiation givesThis increases inward, starting from zero. Marginal convection has ${\nabla }_{\mathrm{rad}}=2/5$, so $w_b^3=136/5$. Since $T_b^{12}/T_e^{12}=w_b^3/64$, the grey-atmosphere convection onset with seventeenth-power opacity isThe purely radiative atmosphere is used only up to this boundary.

The convective relation $P=K{T}^{5/2}$ is spatially uniform in its entropy coefficient $K$ at each instant; $K$ can change as the star radiates and contracts. At the convection boundary $w_b$ is a fixed number, so $P_b\propto g^{1/2}{T}_e^{-8}$ at fixed composition and opacity normalization. ThusOn the other hand, the central scalings for the $n=3/2$ stellar polytrope give $P_c\propto G{M}^2/R^4$ and $T_c=\mu P_c/(\mathcal{R}{\rho }_c)\propto \mu GM/(\mathcal{R}R)$. Hence, at fixed composition,Equating the interior and atmosphere expressions yields $T_e^{21/2}\propto M{R}^{1/2}$. Applying the Stefan–Boltzmann law then proves the Hayashi contraction with seventeenth-power opacity relation

At fixed mass write $L=A{M}^{8/21}{R}^{\alpha }$, with constant $A$ and $\alpha =46/21$. The contraction equation is $\dot{R}=-C{R}^{\alpha +2}$, where $C=7A{M}^{8/21}/(3G{M}^2)$. Direct integration givesMultiplying through by $R^{\alpha +1}$ and substituting the current luminosity givesFor the stated very large initial radius, or after that initial-radius term becomes negligible,This is the luminosity-age relation within the contraction model, rather than an assertion that a physical star begins at literally infinite radius.

Write the photon energy density as $ϵ=a_R{T}^4$, where $a_R=4\sigma /c$ by the Stefan–Boltzmann law. The photon gas has pressure $P=ϵ/3$ and zero photon chemical potential. The first law of thermodynamics therefore givesDividing by $T$ and integrating gives the entropy density of a photon gas:The integration constant is fixed to zero at zero temperature, consistently with the Third law of thermodynamics.

The inequality $\Gamma \gg H$ says interactions act much faster than the expansion time, maintaining local thermal equilibrium. For adiabatic reversible expansion with no additional source of entropy, cosmological entropy conservation fixes the entropy in a comoving volume. Since that volume scales as $a^3$, and the relativistic degrees of freedom are fixed here, $T^3{a}^3$ is constant. ThusDuring radiation domination, the Friedmann equation gives $H^2\propto \rho \propto T^4$, and hence $\boxed{H\propto T^2}$. Fast interactions alone would not exclude an additional entropy-producing process; changing relativistic species likewise changes the simple temperature scaling.

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.