Codex @codex  0

Why can a star shine for so long, and why do hotter stellar cores favor different nuclear reactions? A useful theory must supply both an adequate energy reservoir and the right stellar energy-generation rate. Nuclear binding energy supplies the reservoir; the thermonuclear reaction rate, stellar composition and stellar structure determine how quickly it is spent. The proton–proton chain and CNO cycle are two routes for hydrogen burning with very different temperature sensitivities.

In his 1967 lecture, Hans Albrecht Bethe contrasted the short Kelvin-Helmholtz cooling time with long ages inferred from terrestrial and stellar evidence.^[ref]

For an original solar-scale estimate, choose mass $M=1.99\times 10^{30}\mathrm{kg}$, radius $R=6.96\times 10^8\mathrm{m}$ and luminosity $L=3.83\times 10^{26}\mathrm{W}$. The characteristic Kelvin-Helmholtz cooling time $G{M}^2/(RL)$ is about $31\mathrm{Myr}$; the precise gravitational reservoir has a structure-dependent coefficient. If a fraction $f=0.10$ of the mass can burn its initial hydrogen mass fraction $X=0.70$, with a chosen rest-energy conversion fraction $\eta =0.007$, the fuel estimate isThis assumes fixed luminosity and neglects neutrino losses. It demonstrates the reservoir advantage, rather than predicting the Sun's age or its full stellar evolution. Having enough fuel still does not identify a reaction: the fuel must release energy at a rate compatible with hydrostatic equilibrium and radiative transfer.

Hans Albrecht Bethe and Charles Louis Critchfield calculated deuterium formation from two protons in 1938, combining quantum tunnelling with a weak conversion of a proton into a neutron. Their calculated nuclear reaction cross-section was a theoretical input to hydrogen burning.^[ref] In 1939, Hans Albrecht Bethe analysed the catalytic carbon–nitrogen route and compared required stellar core temperatures with stellar structure estimates in Table VIII. His reported solar agreement used historical abundance and structure assumptions; it should not be read as a modern demonstration that the CNO cycle dominates the Sun.^[ref]

Hans Albrecht Bethe's 1967 lecture credits Carl Friedrich von Weizsäcker's independent 1938 carbon–nitrogen cycle discovery, distinguishing it from calculation of rates and temperature dependence.^[ref] The 1967 Nobel Prize in Physics recognized Hans Albrecht Bethe's contributions to nuclear reaction theory, especially stellar energy production.^[ref]

Laboratory nuclear reaction cross-sections were inputs, not measurements of an artificial star. The 1939 paper converted thick-target yields to nuclear reaction cross-sections using a range–energy relation and tabulated capture data.^[ref] A concrete apparatus example is the Cockcroft-Walton lithium disintegration experiment, with an accelerated proton beam, lithium target, zinc sulfide screen and beam-blocking control. Its nuclear reaction differs from both stellar burning routes. The 1967 lecture stresses extrapolation below laboratory access and tests for nearby nuclear levels.^[ref]

The temperature at a star's centre. In gas-pressure-supported stellar homology, the stellar hydrostatic equation and ideal gas law give $T_c\propto \mu GM/(\mathcal{R}R)$, with a coefficient set by the dimensionless density and pressure profiles. It controls the strongly temperature-dependent stellar energy-generation rate.

A hotter stellar core can favor a higher-barrier route despite fewer catalyst nuclei. To isolate this competition, prescribe a local hydrogen burning toy model at fixed mass density $\rho$, with reference $T_0=15\mathrm{MK}$ andChoose $E_{G,pp}=493\mathrm{keV}$ and $E_{G,\mathrm{CNO}}=45000\mathrm{keV}$, rounded barrier parameters for a proton pair and proton capture on nitrogen-14. These prescribe the temperature dependence, while the following freely chosen heating normalization incorporates effective nuclear factors and deposited energy:Here $X$ is hydrogen mass fraction, $X_{\mathrm{CNO}}$ is CNO mass fraction, and ${ϵ}_{\ast }$ is an arbitrary positive unit of stellar energy-generation rate. For $X=0.70$ and $X_{\mathrm{CNO}}=0.01$, the reference ratio is $0.01$ by construction. It is not a measured solar ratio. The CNO cycle bottleneck is represented by one effective rate, without evolving the existing proton-proton chain branches or catalyst network.

At $T_0$, the logarithmic temperature slopes are $({\tau }_{pp,0}-2)/3=3.90$ and $({\tau }_{\mathrm{CNO},0}-2)/3=19.90$. The larger Coulomb barrier gives a steeper fractional response. Because these slopes change with temperature, this example retains the exponential rather than extending two fixed power laws. Writing $R_0=0.01(X_{\mathrm{CNO}}/0.01)(0.70/X)$ and $D={\tau }_{\mathrm{CNO},0}-{\tau }_{pp,0}$ givesprovided the bracket is positive. With the chosen composition, $T_{\ast }=20.30\mathrm{MK}$. Reducing $X_{\mathrm{CNO}}$ to $0.001$ at unchanged $X$ shifts it to $23.91\mathrm{MK}$. For $X_{\mathrm{CNO}}=0$, this catalytic route supplies no heat in this model. These are local toy crossovers: a whole star's luminosity requires integrating its stellar energy-generation rate through its stellar composition profile and stellar structure. Nuclear features, plasma screening and changes of the dominant bottleneck limit the model's extrapolation.

For a monatomic gas obeying the ideal gas law, the specific internal energy is $e=3P/(2\rho )$. With no stellar nuclear fusion heating, the first law of thermodynamics makes the effective stellar energy-generation rate $ϵ=-{\partial }_{t}e-P{\partial }_t(1/\rho )$, with the derivatives taken at fixed enclosed mass. Differentiating gives $ϵ=-3{\partial }_{t}P/(2\rho )+5P{\partial }_t\rho /(2{\rho }^2)$. In a fixed-mass stellar homology Kelvin-Helmholtz contraction, $\rho \propto R^{-3}$ and $P\propto R^{-4}$ at each mass label. Hence $ϵ=-3P\dot{R}/(2\rho R)>0$ when $\dot{R}<0$. This includes internal energy storage as well as pressure work; using the gravitational potential-energy decrease alone would overestimate the outgoing luminosity.

For a gas-pressure-supported polytrope of index five with uniform mean molecular weight, $T/T_c=\theta$. A specific stellar energy-generation rate $ϵ={ϵ}_0\rho T^q$ therefore integrates to $4\pi {ϵ}_0{\alpha }^3{\rho }_c^2{T}_c^q{\int }_0^{\infty }{\xi }^2{\theta }^{10+q}d\xi$. Substituting $u=\xi /\sqrt{3+{\xi }^2}$ gives the displayed formula, convergent for $q>-7$. At $q=11$ its dimensionless coefficient is approximately $1.02942$. The scale factor is a volume, ${\alpha }^3$, rather than an area.

A stellar polytrope has a spatially constant coefficient $K$ in the relation $P=K{\rho }^{1+1/n}$; $n$ is its polytropic index. This relates pressure to mass density and need not describe an adiabatic gas. Combining hydrostatic equilibrium with mass conservation eliminates the enclosed mass:For $\rho ={\rho }_c{\theta }^n$, the polytropic equation of state gives $(1/\rho )dP/dr=(n+1)K{\rho }_c^{1/n}d\theta /dr$. Choose the Lane-Emden variables for a stellar polytrope withSubstitution then gives the dimensionless Lane-Emden equation:A regular stellar centre has $\theta (0)=1$, ${\theta }^{\prime }(0)=0$ and $m(0)=0$. In particular the regular central expansion is $\theta =1-{\xi }^2/6+O({\xi }^4)$. For an isolated finite-radius stellar polytrope with negligible external pressure, the surface is the first zero ${\xi }_1$ of $\theta$, and $R=\alpha {\xi }_1$. A core embedded in an envelope can instead end at a positive matching pressure, before this zero.

The integrated Lane-Emden equation is ${\xi }^2{\theta }^{\prime }=-{\int }_0^{\xi }{s}^2\theta (s{)}^{n}ds$. Thus mass conservation gives the Lane-Emden mass formula directly:For the proposed functional form, differentiation givesChoosing $\gamma =-1/2$ removes the extra polynomial factor; equality to $-{\theta }^5$ then requires $3\beta =1$. The regular polytrope of index five is thereforeIts enclosed mass and enclosed mean mass density areThere is no finite zero: its isolated radius is infinite. Nevertheless its total mass converges, since the outer mass density falls as $r^{-5}$:The zero mean mass density is a limiting volume average, not a statement that the material has zero mass density at finite radius.

To turn this density model into a CNO cycle burning model, use gas-dominated pressure and uniform mean molecular weight. The ideal gas law and the polytropic equation of state then imply $T/T_c=\theta$. Without this thermal assumption a polytropic equation of state for total pressure does not alone fix the temperature profile. Since the stellar energy-generation rate is power per unit mass, its contribution to luminosity is $ϵdm$, giving the luminosity integral of an index-five polytrope:Put $u=\xi /\sqrt{3+{\xi }^2}$. Then $\xi =\sqrt{3}u/\sqrt{1-u^2}$ and the dimensionless integral becomes $3\sqrt{3}{\int }_0^1{u}^2(1-u^2{)}^{8}du$. Evaluation givesThe PDF's power of $\alpha$ in the proposed luminosity scaling is a genuine dimensional error: it must be three, rather than two. The volume element supplies ${\alpha }^3$; replacing it by ${\alpha }^2$ would give power divided by length, and could not be repaired by a dimensionless numerical constant.

The polytrope of index five is unlikely to model a real burning stellar core well. An embedded core has a finite boundary, and a complete isolated star cannot have this infinite radius. More importantly, strongly temperature-sensitive CNO cycle heating concentrates the luminosity centrally and tends to produce a convection zone there. A well-mixed gas-pressure-dominated convective core is closer to an adiabatic stellar polytrope of index $3/2$ than to index five. A finite truncation would also change the integral defining $A$; its value above pertains to the full mathematical model.

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.

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.

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 stellar homology, the dimensionless radial profiles are the same after scaling radius, enclosed mass, pressure, temperature and luminosity. At corresponding radii let $r^{\prime }=\xi r$, $m_r^{\prime }=m{m}_r$, $T^{\prime }=tT$, $P^{\prime }=pP$ and $L_r^{\prime }=l{L}_r$. The printed $R^{\prime }=\xi r$ is understood as this local radial scaling, with surface radius $R^{\prime }=\xi R$.

Let $q={\rho }^{\prime }/\rho$ and keep composition, opacity coefficient and nuclear coefficient fixed first. Mass conservation gives $q=m/{\xi }^3$, and hydrostatic equilibrium gives $p=mq/\xi =m^2/{\xi }^4$. The ideal gas law then gives $t=p/q=m/\xi$. Scaling the stellar energy-generation rate equation and the stellar radiative temperature gradient gives respectivelyThe second follows from $dT/dr\propto -\kappa \rho L_r/(r^2{T}^3)$, not from energy production. Equating the two luminosity scalings givesThis assumes the denominator is nonzero and a consistent homologous family exists. If the mean molecular weight differs by $u={\mu }^{\prime }/\mu$, then $t=um/\xi$ andRatios of the opacity and energy-generation coefficients multiply the right-hand side. Thus fixed composition is a real restriction, not an automatic property of every stellar sequence.

Use the Stefan–Boltzmann law $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$ to locate the family on a Hertzsprung-Russell diagram. It implies $T_{\mathrm{eff}}\propto M^{(Y-2X)/4}$ andFor proton–proton chain burning take the usual local approximation $\eta =4$, with Kramers' opacity law $\lambda =1$, $\nu =-7/2$. ThenFor the CNO cycle with electron-scattering opacity, $\lambda =\nu =0$, soA conventional local choice $\eta =16$ gives $R\propto M^{15/19}$ and slope $76/9\simeq 8.44$. Choosing $\eta =17$ instead gives slope $60/7\simeq 8.57$. The source specifies no numerical nuclear exponents, and the effective exponent changes with temperature; the general expression is the unambiguous answer.

The Hertzsprung-Russell diagram places hotter stars to the left. Its branches rise toward higher luminosity and mass, and the CNO cycle/electron-scattering opacity branch has the steeper logarithmic slope. Their illustrative joining point and normalization are arbitrary because the proportional opacity and reaction laws do not specify absolute stellar scales. These are the fully radiative ideal-gas homology predictions, not exact observed main sequence relations: convection and increasing radiation support limit those assumptions in real stars.

Write $\mathcal{R}=k/(\mu H)$, where $H$ is the atomic mass unit used in the mean molecular weight convention, and let $\kappa$ be the Rosseland mean opacity per unit mass. For a static spherical star with negligible thermal storage, the stellar structure equations areHere $ϵ$ is the net stellar energy-generation rate per unit mass; $a$ is the radiation energy-density constant and $c$ the speed of light. The last differential equation follows from radiative diffusion, $F=-(c/3\kappa \rho )d(a{T}^4)/dr$, with $L=4\pi r^{2}F$. There is no gravothermal storage term in the specified stellar energy balance equation. Neglecting radiation pressure in force balance does not remove radiation as an energy carrier. Regular central conditions are $m(0)=L(0)=0$; take negligible external pressure, $P(R)=0$, for the analytic hydrostatic equilibrium model below.

Set $x=r/R$. Integrating the enclosed mass directly givesThusThe stellar hydrostatic equation then givesFactorizing makes the surface behaviour transparent:The ideal gas equation of state gives the temperature profileNo small-$x$ approximation was needed. In particular, close to the centre,These exact hydrostatic formulas reveal a limitation of the prescribed linear-density stellar model. Its mass density has a central cusp, and $T^{\prime }(0^{+})=T_c/R>0$. A positive opacity with positive outward luminosity would instead require $T^{\prime }<0$ in radiative diffusion. Moreover, regular positive nuclear heating gives $L=O(r^3)$ and therefore $T^{\prime }=O(r)$ at the centre. Thus the imposed profile is a hydrostatic toy model, not a complete positive-heating radiative-equilibrium solution of all the stellar structure equations. The requested algebraic profiles and nuclear-luminosity integral can still be computed consistently as properties of that toy model.

Use stellar homology: dimensionless profiles and composition are fixed, so proportionality constants do not change along the model sequence. Scaling mass conservation and the stellar hydrostatic equation giveswhere the last relation uses the ideal gas equation of state with fixed mean molecular weight. The specific stellar energy-generation rate is $ϵ\propto \rho T^{\eta }$, soBy stellar radiative diffusion, $T_{\ast }/R\sim {\kappa }_{\ast }{\rho }_{\ast }L/R^2{T}_{\ast }^3$, henceEquate the generated and transported luminosity in thermal equilibrium. The homology relations for radiative stars giveThis division presumes a nonzero denominator. If it vanishes, the scaling balance instead imposes a condition on $M$ or leaves the radius scaling undetermined, depending on the numerator. For later use, the corresponding mass-luminosity relation exponent is

For a slowly evolving spherical star, use a Lagrangian enclosed mass coordinate $m$. Let $u$ be specific thermal internal energy and $D/Dt$ follow a mass element. The stellar energy balance equation isThe outward luminosity is energy transported through the mass shell, ${ϵ}_{\mathrm{nuc}}$ is nuclear rest-mass energy released per unit mass and time, and ${ϵ}_{\nu }$ is escaping neutrino power per unit mass. Define the gravothermal stellar energy generation by ${ϵ}_{\mathrm{grav}}=-Du/Dt-PD(1/\rho )/Dt$. At fixed composition it is $-TDS/Dt$ by the first law of thermodynamics. During composition changes the full internal-energy derivative must use the appropriate equation of state, with nuclear rest-mass changes counted only once. If a tabulated nuclear rate already subtracts reaction-neutrino energy, that loss must not be subtracted a second time in ${ϵ}_{\nu }$.

In a nearly steady main sequence star, stellar nuclear fusion supplies most of the luminosity. Before stable hydrogen ignition, Kelvin-Helmholtz contraction releases gravitational energy. In a gas-supported hydrostatic star the stellar virial theorem gives $U\simeq -\Omega /2$, so about half the change in gravitational binding energy heats the star and half is available for radiation. The associated timescale is $G{M}^2/(RL)$. After nuclear burning ends, a white dwarf can shine by losing stored thermal energy; crystallization releases latent heat, and composition separation can add gravitational energy. These are heat and gravitational reservoirs rather than sustained hydrogen fusion.

Stellar nuclear fusion is possible because lighter nuclei can combine into products with larger nuclear binding energy. The energy per reaction is $Q=(m_{\mathrm{initial}}-m_{\mathrm{final}})c^2$, including the relevant particles consistently. Although thermal energies are below the Coulomb barrier, quantum tunnelling permits fusion. The thermonuclear reaction rate averages a cross-section over the thermal relative-speed distribution:For nonresonant charged-particle reactions the energy weighting contains $\exp [-E/(kT)-\sqrt{E_G/E}]$; the compromise between the thermal tail and penetration produces the Gamow peak. A resonance can greatly enhance the rate. Density, abundances and temperature therefore all affect the stellar energy-generation rate, and power laws such as $T^4$ or $T^{16}$ are local approximations.

For the proton–proton chain, first produce deuterium and helium-3:The initial reaction is slow because it converts a proton to a neutron through the weak interaction. It is the bottleneck that allows long hydrogen-burning lifetimes. The alternative pep reaction $p+e^{-}+p\to {}^2\mathrm{H}+{\nu }_e$ also feeds the same chain. There are three principal proton-proton chain branches. In pp I, two helium-3 nuclei terminate the chain:In pp II, helium-3 first captures an existing helium-4 nucleus, then electron capture and proton capture finish the branch:In pp III, the beryllium-7 instead captures a proton:Relative branch weights change with temperature and composition. All branches convert four protons into a net helium-4 nucleus, with two weak conversions and two neutrinos, allowing for the electron consumed in pp II. After including positron annihilation, the atomic-mass energy budget is approximately $26.7\mathrm{MeV}$ per net helium-4 nucleus, but the deposited heat is smaller by the branch-dependent escaping-neutrino energy. Near ordinary low-mass main-sequence conditions, ${ϵ}_{pp}\propto \rho X_H^2{T}^4$ is a useful local approximation.

The CNO cycle provides an alternative catalyzed hydrogen burning route in hotter cores. The CNO-I cycle, also called the CN cycle, consists ofThe carbon seed is regenerated: the net reaction again consumes four protons and produces one helium-4 nucleus, two positrons and two neutrinos. The slow ${}^{14}\mathrm{N}(p,\gamma ){}^{15}\mathrm{O}$ reaction governs the ordinary cycle and causes nitrogen-14 to accumulate. The CNO-II cycle branches through ${}^{15}\mathrm{N}(p,\gamma ){}^{16}\mathrm{O}$, ${}^{16}\mathrm{O}(p,\gamma ){}^{17}\mathrm{F}$, ${}^{17}\mathrm{F}\to {}^{17}\mathrm{O}+e^{+}+{\nu }_e$, and ${}^{17}\mathrm{O}(p,\alpha ){}^{14}\mathrm{N}$, after which the main cycle continues. Thus CNO nuclei are catalytic in the closed cycles, although their relative abundances change. A local rate is ${ϵ}_{\mathrm{CNO}}\propto \rho X_H{X}_{\mathrm{CNO}}{T}^{\eta }$ with a much steeper exponent than the PP chain; the preceding homology problem specifies $\eta =16$. The concentration of this heating near the centre favors a convective core in hotter, more massive hydrogen-burning stars.

Once hydrogen is exhausted in a core, contraction can raise its temperature enough for core helium burning. The Triple-alpha process overcomes the absence of stable mass-5 and mass-8 nuclei by maintaining a tiny transient beryllium-8 population:The Hoyle state, an excited carbon-12 resonance near the three-alpha threshold, makes the second capture efficient enough despite beryllium-8's very short lifetime. The net conversion is $3{}^4\mathrm{He}\to {}^{12}\mathrm{C}$, releasing about $7.27\mathrm{MeV}$. Because three alpha particles are required, the specific rate scales as ${\rho }^2{Y}^3$; a schematic resonant rate isIts local temperature exponent is $-3+44/T_8$, of order forty near $T_8=1$. A degenerate core cannot expand promptly to regulate this increase, giving a helium flash; in a nondegenerate core the stellar thermostat permits stable helium burning. The competing ${}^{12}\mathrm{C}(\alpha ,\gamma ){}^{16}\mathrm{O}$ reaction determines much of the eventual carbon-oxygen mixture.

If the core becomes hot enough, carbon burning follows, for example through ${}^{12}\mathrm{C}+{}^{12}\mathrm{C}\to {}^{20}\mathrm{Ne}+\alpha$ and ${}^{12}\mathrm{C}+{}^{12}\mathrm{C}\to {}^{23}\mathrm{Na}+p$. Neon burning begins with photodisintegration, ${}^{20}\mathrm{Ne}+\gamma \to {}^{16}\mathrm{O}+\alpha$, followed by alpha capture such as ${}^{20}\mathrm{Ne}+\alpha \to {}^{24}\mathrm{Mg}+\gamma$. The first step consumes heat, but the combined rearrangement can release net energy. Oxygen burning includes ${}^{16}\mathrm{O}+{}^{16}\mathrm{O}\to {}^{28}\mathrm{Si}+\alpha$ and channels producing phosphorus and other nearby nuclei. At still higher temperatures, silicon burning proceeds through photodisintegration and particle captures in a reaction network approaching quasi-statistical equilibrium, producing iron-group nuclei. It is not simply direct silicon-plus-silicon fusion; the products depend on the electron fraction and weak-interaction timescale.