Codex @codex  0

An Algol binary is a semidetached binary in which an evolved, less massive donor star fills its Roche lobe while its more massive mass gainer remains near the main sequence. Binary mass-ratio reversal resolves the Algol paradox. Eclipsing members commonly have Algol variable light curves, but that observational class also includes detached systems.

The Algol paradox is the apparently inconsistent evolutionary ordering in an Algol binary: the less massive star is more evolved, whereas higher-mass isolated main sequence stars normally exhaust core fuel sooner. Binary mass-ratio reversal resolves it by identifying the present donor as the initially more massive star.

Case A mass transfer starts while the donor star is on the main sequence, with core hydrogen burning. Its onset precedes giant-envelope formation; donor structure and the binary mass ratio still determine stability.

The Hertzsprung gap is a sparsely populated region of the Hertzsprung-Russell diagram traversed rapidly by many stars evolving away from the main sequence toward the red giant region.

The main-sequence turnoff is the location in a coeval population's Hertzsprung-Russell diagram where stars are leaving the main sequence after exhausting their core hydrogen. More massive stars generally have larger luminosity per unit mass and shorter nuclear lifetimes, so the most massive stars still on the main sequence constrain the population's age. In a simple model, equating the stellar nuclear timescale $t\sim qM/L(M)$ to the common age determines a turnoff mass; detailed calibration requires composition and stellar models.

A mass-gaining star accretes material from a companion. Its changing mass and composition can alter its main sequence location and cause stellar rejuvenation.

The radiative family is a stylized approximation to gas-pressure-supported main sequence stars with mostly radiative interiors and moderately temperature-sensitive hydrogen burning, such as the roughly solar-mass regime. Its opacity exponents are local approximations, not a universal composition-independent opacity law. The fully convective family is most relevant to the low-mass main sequence, below the transition to a radiative core. Its steep surface temperature dependence resembles the role of negative hydrogen ion opacity in setting a Hayashi track; real surface opacity is not exactly the supplied density-linear tenth-power law everywhere.

At higher mass, the CNO cycle replaces the proton–proton chain as the dominant source of stellar nuclear fusion and has a much larger local temperature exponent. This promotes a convective core. Hot highly ionized matter increasingly has electron-scattering opacity, rather than the specified density/temperature-dependent absorption. Radiation pressure becomes important in the equation of state, changing $T_c\propto M/R$ and eventually making the luminosity relation approach an Eddington luminosity scaling. Composition gradients and strong mass loss also break strict homology.

At lower mass and effective temperature, partial ionization and thermal dissociation change the heat capacity and adiabatic index; molecules and eventually condensates alter the surface opacity. Increasing electron degeneracy pressure invalidates pure ideal-gas support. Below the hydrogen-burning minimum mass, sustained hydrogen-burning equilibrium is absent, so brown dwarfs cool instead of following a zero-age hydrogen-burning sequence. These changes explain why neither derived slope should be extrapolated indefinitely.

For a circular binary with negligible spin, the two orbital radii about the center of mass are $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$. Summing their orbital angular momenta givesUsing Kepler's third law, ${\Omega }^2=GM/a^3$, also gives $J=M_1{M}_2\sqrt{Ga/M}$. Under conservative binary mass transfer, $M$ and $J$ stay fixed, so $a\propto (M_1{M}_2{)}^{-2}$. Since $d\log M_2/d\log M_1=-q$, the Roche-lobe radius response exponent isFor the supplied giant response, luminosity stays fixed during the mass-loss perturbation and ${\zeta }_{\ast }=d\log R/d\log M_1=-0.27$. The change of the logarithmic overfill is $({\zeta }_{\ast }-{\zeta }_L)d\log M_1$. Because the donor loses mass, this reduces overfill and stabilizes transfer precisely when ${\zeta }_{\ast }>{\zeta }_L$. Thus the stated model requiresThis is a response criterion within the given radius law; a different adiabatic donor response would change the dynamical threshold.

The initially more massive star evolves first. If it waits until it has the supplied giant structure while still having $q>1$, conservative overflow is unstable by this criterion. If Roche-lobe overflow instead begins while the donor is on the main sequence or shortly after central hydrogen exhaustion, its largely radiative-envelope response can allow transfer before a deep convective giant envelope develops. Mass loss and accretion can then reverse the mass ratio before it reaches the giant state. This is the binary mass-ratio reversal route to an Algol binary, through Case A mass transfer or early case B mass transfer.

At lobe filling, the lobe approximation and Kepler's law give the pre-giant Roche-lobe-filling period limitIt increases with the donor's radius. Requiring $R_L<R_{\mathrm{BGB}}=1.68R_{\odot }(M_1/M_{\odot }{)}^{5/3}$ therefore yieldsThe companion mass cancels in this lobe approximation.

Since $P\propto a^{3/2}$ at fixed total mass, the separation scaling gives $\boxed{P(M_1{M}_2{)}^3=\text{constant}}$. Use solar-unit masses $m_i=M_i/M_{\odot }$ and $m=m_1+m_2$, and let $K=P(m_1{m}_2{)}^3$. ThenThe derivative vanishes at $m_1=5m/8$, and its derivative is $5/m_1^2+3/(m-m_1{)}^2>0$, while the ratio diverges at both endpoints. Hence this is the global minimum, the conservative pre-giant period-ratio minimum.

For the stated Algol parameters, $m=4$ and $K=3(1\cdot 3{)}^3=81$ days. At the minimum $m_1=2.5,m_2=1.5$, the period is $81/(2.5\cdot 1.5{)}^3=1.536$ days and $P/m_1^2=0.24576$ days, well below $P_0=0.8$ days. This is a possible initially more-massive donor configuration with overflow before the giant branch. The current $q=1/3$ is also below the giant stability threshold. Thus Algol is compatible with the proposed conservative early-overflow route, although this necessary geometric test alone does not prove its full evolutionary history.

For RT Lac, $m=2.4$ and $K=5.1(0.9\cdot 1.5{)}^3$. Its minimum occurs at $m_1=1.5,m_2=0.9$, where the mass product equals its current value, so $P=5.1$ days andNo conservative antecedent on this orbit can begin overflow before the giant branch. Beginning with an initially more massive giant instead encounters the instability just derived. RT Lac cannot follow the stated coeval, conservative Algol-formation model. A plausible alternative is rapid initial transfer with substantial mass and angular-momentum loss, via an outflow, strong wind or envelope-ejection episode. Such nonconservative binary mass transfer changes both total mass and $J$ and removes the period-product constraint. Its present mass ratio $0.9/1.5=0.6$ permits stable current giant transfer. The given data do not determine a unique mass-loss history or establish a specific ejection mechanism.

Nuclear burning explains both the large available energy and its release under stellar-core conditions. Converting hydrogen into helium releases approximately $0.007$ of the fuel rest-mass energy. The net proton–proton chain reaction can be written $4p+2e^{-}\to {}^4\mathrm{He}+2{\nu }_e$, with about $26.7$ MeV available per helium nucleus, including the neutrino energy loss. The CNO cycle achieves the same net hydrogen-to-helium conversion with catalysts. Measured nuclear reaction rates and stellar-core temperatures give rates capable of supplying the observed luminosity.

The energy budget distinguishes stellar nuclear fusion from alternatives. Kelvin-Helmholtz contraction has a characteristic Kelvin-Helmholtz cooling time $t_{\mathrm{KH}}\sim G{M}^2/(RL)$, only tens of millions of years for a solar-type star. Geological ages and stellar-population ages are much longer. Chemical binding energy per unit mass is smaller still. Nuclear fuel in a substantial fraction of the core gives a stellar nuclear timescale of billions of years, explaining long-lived main-sequence stars.

Main sequence luminosities, the correlation of lifetime with mass, and cluster turnoffs are reproduced by models in which central hydrogen is consumed and helium accumulates. Spectroscopic abundances and the appearance of helium-, carbon- and heavier-element-rich evolved material agree with sequential nuclear processing and stellar dredge-up. These observations support the reaction network and evolutionary consequences, rather than independently measuring every reaction in every star.

Solar neutrinos provide direct information from the burning core because they escape without the long radiative diffusion delay of photons. Their measured energies and fluxes test weak reactions in the hydrogen-burning network. The historical deficit of electron neutrinos did not require abandoning nuclear power; flavor conversion changes the measured electron component. Later flavor-sensitive measurements confirmed the total flux expected from solar nuclear burning, as described in the SNO collaboration's first-results report. Thus the evidence combines an adequate energy source, quantitatively successful stellar models and a direct nuclear by-product. Cooling remnants and brief gravitationally powered evolutionary phases are exceptions to the statement that a star's instantaneous luminosity is nuclear.

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.

A constant number star formation rate over $10\mathrm{Gyr}$ produces a uniform present-age distribution. The birth-age distribution under constant star formation therefore givesOutside this age interval the cumulative distribution function is zero or one. The census here retains stellar remnants as objects, as required by the population model.

Write $m=M/M_{\odot }$. The initial mass function is proportional to $m^{-3}$ on $m\ge 0.2$. Normalizing the tail integral gives the power-law initial-mass-function tail coordinate:For lower thresholds $X=1$. Because the birth mass and formation time have independence, the probability integral transform makes $(X,Y)$ uniform on the unit square. $X$ is an upper-tail coordinate, so larger masses have smaller $X$.

A red giant of current age $t$ lies in the stellar lifetime regions in a mass-age diagramThere are no red giants below $m=\sqrt{0.8}$ in this model. For $\sqrt{0.8}<m<1$ the upper boundary is the population's age limit, while for $m\ge 1$ both lifetime boundaries are present. A white dwarf lies above $t/\mathrm{Gyr}=10/m^2$.

In uniform probability coordinates the main-sequence boundary is $Y=20X$ and the end of the giant phase is $Y=25X$. The giant region is therefore $20X\le Y<25X$, clipped to the unit square: it is the triangle with vertices $(0,0),(1/25,1),(1/20,1)$. Its area, and hence the individual giant fraction, isThe white dwarf region $Y\ge 25X$ is the triangle with vertices $(0,0),(0,1),(1/25,1)$, givingThe remaining $97.5\%$ are on the main sequence. Boundaries have zero probability and do not affect these counts. Treating every post-giant object as a white dwarf is a stipulated toy-model assumption; real high-mass stars have other stellar remnants.

For the binary calculations, use a coeval binary population: components share their birth time even though their birth masses are independently drawn and their subsequent evolution is independent. Hence $(X_1,X_2,Y)$ is uniform in the unit cube. At a fixed age coordinate $Y=y$, a component is a giant when $y/25<X<y/20$, and a white dwarf when $0<X<y/25$. Define the conditional state fractionsComponents have conditional independence given $y$; averaging over the common age must be done after forming the conditional pair probabilities. The continuous-birth binary state fractions below are consequently volumes in this cube, not products of the marginal single-star fractions.

The elements have several complementary origins: primordial light nuclei, stellar and explosive fusion products, and nuclei built by neutron capture. Chemical enrichment also requires their ejection, rather than merely their production inside a star. Nucleosynthesis and transport into the interstellar medium must therefore be considered together.

Hydrogen and much of the helium are primordial. During Big Bang nucleosynthesis, a cooling, expanding plasma assembled surviving neutrons and protons into deuterium, helium-3 and especially helium-4, with small abundances of lithium-7 or its beryllium-7 progenitor. Unstable nuclei at mass numbers five and eight, together with the rapidly falling density and short available time, obstructed sustained production of heavier nuclei. The carbon, oxygen and heavier material in later generations consequently require other sites. Cosmic rays fragment heavier interstellar nuclei and contribute importantly to the light elements lithium, beryllium and boron, which ordinary stellar burning tends to destroy.

Burning stages depend on initial stellar mass. Stellar nuclear fusion can release energy when the products have greater nuclear binding energy than the reactants. The resulting stellar thermostat links nuclear burning to hydrostatic equilibrium: when fuel exhaustion lowers heating, contraction raises the central temperature until another fuel can ignite, if the core becomes hot enough. Higher nuclear charges bring larger Coulomb barriers, so later charged-particle burning generally requires progressively higher temperatures. Whether the core reaches them depends on its mass, electron degeneracy pressure, mixing and envelope loss.

An approximately solar-composition $1M_{\odot }$ star first spends most of its life on the main sequence, converting hydrogen into helium mainly through the proton–proton chain. The net conversion is four protons into one helium-4 nucleus, with positrons and neutrinos; the small CNO cycle contribution redistributes CNO catalyst abundances but is not a net source of carbon at this stage. Exhaustion of central hydrogen leaves a helium-rich core. A surrounding hydrogen-burning shell continues to add helium, while the core contracts and the envelope expands into a red giant.

The low-mass helium core becomes supported substantially by electron degeneracy pressure. At roughly $10^8\mathrm{K}$, the Triple-alpha process ignites, first making transient beryllium-8 and then carbon-12 through the Hoyle state. Because degeneracy initially weakens the usual expansion response, this ignition is a helium flash. It heats and lifts the degeneracy of the core, leading to a stable phase of core helium burning. The Triple-alpha process makes carbon; competing carbon-12 alpha capture makes oxygen:After central helium exhaustion, the star develops a carbon-oxygen core, a helium-burning shell and a more external hydrogen-burning shell. On the Asymptotic giant branch, recurrent thermal pulses of an asymptotic-giant-branch star arise from unstable helium-shell burning. Stellar dredge-up can bring processed material toward the surface, with its efficiency depending on mass and composition. The eventual stellar wind removes the envelope; the hot exposed core can illuminate a planetary nebula, and the remnant becomes a carbon-oxygen white dwarf. A usual isolated $1M_{\odot }$ star does not ignite carbon burning. In particular, an element synthesized in its core need not all escape: much of the carbon and oxygen remains locked in the white dwarf.

An approximately $5M_{\odot }$ star burns central hydrogen predominantly through the CNO cycle, usually in a convective core. The CNO nuclei act as catalysts; the slow nitrogen-14 proton-capture step tends to accumulate nitrogen-14 while converting hydrogen into helium. After central fuel exhaustion, a hydrogen-burning shell surrounds the contracting helium core. Unlike the $1M_{\odot }$ case, core helium burning normally begins before strong electron degeneracy, without a classical helium flash. The same Triple-alpha process and carbon-12 alpha capture create a carbon-oxygen core.

The $5M_{\odot }$ star then enters the Asymptotic giant branch, with hydrogen and helium burning in shells. Repeated thermal pulses and stellar dredge-up can enrich its envelope with freshly made carbon and s-process nuclei. Hot-bottom burning can occur when the base of the convective envelope is sufficiently hot: the CNO cycle processes dredged-up carbon into nitrogen, so a massive AGB star need not become carbon-rich despite making carbon below its envelope. Helium-shell reactions on neon-22 can also provide neutrons for the s-process. For an ordinary $5M_{\odot }$ model, strong envelope loss generally precedes carbon ignition, leaving a carbon-oxygen white dwarf. The exact threshold for carbon ignition and the size of these yields depend on metallicity, mass loss and mixing; they are not universal integers fixed solely by initial mass.

A $20M_{\odot }$ star reaches much more advanced stages. It starts with convective-core CNO cycle hydrogen burning, followed by core helium burning. After central helium exhaustion, contraction can ignite carbon burning at temperatures approaching $10^9\mathrm{K}$. Carbon-fusion channels make products including neon, sodium and magnesium. Neon burning then uses photodisintegration of neon-20 to release oxygen-16 and an alpha particle, followed by capture of that alpha particle on another neon-20 nucleus; the combined network releases energy and produces magnesium-24. Oxygen burning makes silicon-group nuclei through channels such asAt still higher temperature, silicon burning rearranges nuclei by photodisintegration and captures of protons, neutrons and alpha particles. It is not principally direct fusion of two silicon nuclei. Networks approach quasi-equilibrium and then nuclear statistical equilibrium, producing an iron-group core whose precise mixture depends on the electron fraction in stellar matter. Outer shells retain different stages of burning, giving an approximately layered structure. Increasing stellar neutrino energy loss makes advanced burning stages much shorter than hydrogen burning.

Near the iron group, nuclear binding energy per nucleon is close to its maximum. Further fusion cannot provide the sustained energy source required to support the core. As the core grows, electron capture reduces electron-pressure support and photodisintegration absorbs energy, assisting gravitational collapse. A successful core-collapse supernova ejects some of the envelope and processed shells, and leaves a neutron star or black hole. The shock also causes explosive burning, changing the pre-existing abundances. A $20M_{\odot }$ initial mass does not uniquely fix the explosion or remnant: envelope loss, rotation, composition and fallback all matter, and some collapses fail to eject much material.

Ejection determines the enrichment of the interstellar medium. Low- and intermediate-mass stars return helium, CNO-processed matter, dredged-up carbon and s-process products through stellar winds and envelope expulsion. Massive-star winds can expose hydrogen- or helium-burning products before collapse. Successful core-collapse supernovae eject substantial oxygen and other alpha-capture products, together with some iron-group matter formed in hydrostatic and explosive burning. Material that falls back or remains in a stellar remnant is withheld from the interstellar medium. Subsequent mixing, cooling and star formation distribute the escaped material into new stars. The distinction between production and ejection is why an internal burning inventory alone is not a stellar yield.

Binary stars provide an additional major iron-group source. In a Type Ia supernova, a carbon-oxygen white dwarf in a binary undergoes thermonuclear disruption after an appropriate accretion or merger history. Explosive carbon and oxygen burning makes intermediate-mass and iron-group nuclei. Much initially synthesized nickel-56 decays through cobalt-56 into iron-56; this radioactive decay also powers the supernova light curve. Such explosions have delay times and yields different from those of massive-star collapse, and help explain why iron enrichment need not track oxygen enrichment exactly.

Neutron capture builds much of the material beyond iron. Charged-particle reactions encounter increasing Coulomb barriers as nuclear charge grows, but a neutral neutron does not encounter that barrier. Neutron capture raises mass number, and subsequent beta-minus decay can turn a neutron into a proton and raise atomic number. The relative capture and decay times select different paths through the nuclear chart.

In the slow neutron-capture process, most unstable intermediates decay before the next capture, so the path stays near stable nuclei. In Asymptotic giant branch stars, important neutron sources areThe carbon-13 source operates in suitable partially mixed layers, while the hotter helium-shell pulses can activate the neon-22 source. Captures on pre-existing seed nuclei produce elements such as strontium and barium and, under suitable conditions, lead. Stellar dredge-up and stellar winds return these products to the interstellar medium. Massive stars also have a weak s-process, notably using neon-22 during helium- and carbon-burning phases; it need not produce the same abundance range as the main AGB contribution.

In the rapid neutron-capture process, very high neutron densities permit many captures before typical beta-minus decays. The path moves to neutron-rich nuclei, which later decay toward stability, producing heavy nuclei including lanthanides and actinides. Neutron-rich ejecta from neutron star mergers provide an established site: a merger can eject dynamical material and disk outflows, and radioactive decay of the newly made nuclei powers a kilonova. Special neutron-rich outflows associated with some massive-star collapses are additional candidate sites. It would be unjustified to assume that every ordinary core-collapse supernova produces the full heavy r-process pattern; the required neutron richness and expansion history are restrictive.

The two neutron-capture processes do not make every isotope beyond iron. Some proton-rich isotopes bypass their main paths. Gamma-process nucleosynthesis, for example, uses photodisintegration of pre-existing heavy seed nuclei in hot explosive material. Other proton-rich reaction paths can matter in suitable environments. Thus “heavier than iron” describes a family of production mechanisms rather than a single further ordinary fusion stage.

Binaries change both stellar evolution and the delivery of elements. Roche-lobe overflow, wind transfer and mergers can strip a donor, rejuvenate a mass-gaining star, alter core and envelope masses, and change the timing and character of an eventual explosion. Binary stripping can remove the hydrogen envelope of a massive core-collapse supernova progenitor. Accretion onto a white dwarf can cause a classical nova, ejecting hydrogen-burning products and particular isotopes while ordinarily leaving the white dwarf intact; this differs from the destructive Type Ia supernova. Compact binary stars can lose orbital energy through gravitational radiation and merge, providing the neutron star merger route to heavy r-process material. Consequently chemical enrichment cannot be inferred entirely from isolated-star tracks, even when a simplified population calculation treats binary components as evolving independently.

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

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

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

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

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

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

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.

A solar-type main sequence star has a radiative interior and an outer convection zone. Rotating convection and differential rotation continually regenerate its magnetic field by dynamo action. The field is therefore not simply an initially stored field undergoing passive magnetic diffusion. In an alpha-Omega dynamo, the Omega effect winds a poloidal magnetic field into a toroidal magnetic field, while the alpha effect associated with helical motions can regenerate poloidal flux. In a solar dynamo, the Babcock-Leighton mechanism provides another poloidal source through the emergence and dispersal of tilted active regions. Shear near the tachocline and transport through the convection zone help organize the cycle; no single dynamo model is needed to claim that rotation and convection are essential ingredients.

Rotation acts on convective motions through the Coriolis force, altering their correlations and helicity. A useful activity measure is the stellar activity Rossby number, ${\mathrm{Ro}}_{\ast }=P_{\mathrm{rot}}/{\tau }_c$, where ${\tau }_c$ is a convective turnover time. At a given stellar structure, faster rotation ordinarily increases rotational influence and produces greater magnetic activity in the unsaturated regime. Young solar-type main sequence stars commonly rotate rapidly and show strong sunspots, flares, and chromospheric and coronal emission, although their initial rotation rates have substantial scatter. Very rapid rotators enter a saturated regime, so activity does not increase without limit as rotation increases. The saturation mechanism is not uniquely fixed by the correlation itself. The observed unsaturated correlation and saturated regime are documented in primary rotation–activity measurements.

The reverse interaction is wind-driven magnetic braking of a solar-type star. An ionized stellar wind is coupled to the star's open magnetic field. Matter carries angular momentum, and the Maxwell stress tensor also transmits a torque. Magnetic coupling provides an effective lever arm comparable to the Alfvén radius, rather than only the stellar radius. Writing ${\dot{M}}_w>0$ for the outward mass-loss rate, a schematic angular-momentum equation iswhere a geometry factor can be absorbed into an effective $R_A$. This lever-arm interpretation follows from the angular-momentum invariant of a magnetized wind; it is the central result of the original magnetized solar-wind calculation. Magnetic braking can be effective despite a small mass-loss rate, because $R_A$ can greatly exceed the stellar radius.

One can derive the familiar rotation-age scaling in a simple unsaturated model. For a radial open field, take $B_r(R_A)\simeq B_{\ast }(R_{\ast }/R_A{)}^2$, ${\dot{M}}_w\simeq 4\pi R_A^2{\rho }_A{v}_w$, and $v_w^2\simeq B_r(R_A{)}^2/({\mu }_0{\rho }_A)$. Eliminating ${\rho }_A$ givesIf the unsaturated large-scale field satisfies $B_{\ast }\propto \Omega$ and the other factors vary slowly, this gives $\dot{J}=-K{\Omega }^3$ with positive K. Approximating I as constant on this part of the main sequence, integration yieldsThus $\boxed{\Omega \propto t^{-1/2},P_{\mathrm{rot}}\propto t^{1/2}}$ in the late-time part of this approximation. This Skumanich rotation law captures the classic observed secular decline of rotation and magnetic activity; its observational basis is the original rotation and chromospheric-emission analysis. It is a scaling model, not an exact trajectory for every star. With a saturated field, the same simplified wind model instead gives a torque approximately proportional to $-\Omega$, so initially rapid rotation can decay more nearly exponentially until the unsaturated regime is reached.

Over the main sequence lifetime, wind loss therefore generally makes the star rotate more slowly, increases its stellar activity Rossby number, and reduces its unsaturated mean activity. Later changes of moment of inertia, wind properties and open-field topology alter the simple power law. Internal transport of angular momentum can also feed the braked envelope from the radiative interior, so surface rotation is not automatically a measure of the entire star's angular momentum. In a star with weak braking, this secular history differs from the ideal cubic-torque model; an age law alone is not a fundamental dynamo equation.

The stellar rotation-activity feedback also acts on shorter time scales. A growing magnetic field exerts a Lorentz force that changes differential rotation and convective motions, extracting kinetic energy and limiting further growth through dynamo quenching. Conversely, varying shear and poloidal regeneration change the magnetic cycle. The present solar-like regime can sustain repeated solar cycles, with polarity reversals and activity modulation superposed on the much slower rotational evolution. Rotation promotes dynamo activity; dynamo-generated fields brake rotation through the wind and react back on the motions that generate them. This coupled feedback explains the broad progression from young, rapidly rotating, strongly active stars toward older, more slowly rotating, less active solar-type stars, while allowing cycles, irregular modulation and different initial histories.