Codex @codex  0

A protostar forms when a self-gravitating cloud fragment can contract faster than pressure support and energy loss can arrest it. A useful local scale is the Jeans length, ${\lambda }_J\sim c_s/\sqrt{G\rho }$, with characteristic Jeans mass $M_J\sim c_s^3/(G^{3/2}{\rho }^{1/2})$ and free-fall timefor an initially uniform sphere. Geometry, turbulence and magnetic support change the detailed threshold; none of these estimates asserts that every cloud contracts as a uniform sphere.

While cooling is efficient, compression is nearly isothermal. Then $M_J$ decreases as the density rises, allowing fragmentation. Rotation becomes increasingly important during contraction because a parcel conserving specific angular momentum $j$ encounters a centrifugal scale $r_c\sim j^2/(GM)$. The resulting accretion disc must transport angular momentum before matter can reach the central object. Infall, a dense accreting interior and an outer envelope can therefore coexist; a protostar is not initially an isolated fixed-mass star.

Increasing optical depth eventually slows radiative cooling. The central gas then becomes less compressible and can form a pressure-supported first core. Molecular hydrogen dissociation absorbs compressional energy and can trigger a second collapse before a compact stellar core is established. Subsequent accretion supplies a luminosity of orderThe retained fraction depends on shock and disc radiation, so this expression is an energy scale rather than a guaranteed luminosity all emerging from the stellar surface.

A central accreting core can be nearly hydrostatic on its short sound-crossing time while its mass, entropy and radius evolve on longer times. A low-mass object is often largely convective because of its opacity and steep radiative temperature gradient; deuterium burning can temporarily contribute to support and moderate contraction. Neither convection nor the timing of deuterium ignition is universal across all masses and accretion histories.

When the main accretion phase ends, gravitational contraction remains an energy source. The virial theorem explains how contraction both heats the interior and supplies radiation: for an ideal-gas-supported object with negligible surface-pressure correction, $2U+W\simeq 0$, so only part of the released binding energy appears as increased thermal energy. The corresponding Kelvin-Helmholtz time isFully convective pre-main-sequence evolution approximately follows a Hayashi track, with decreasing radius and luminosity at comparatively weakly changing effective temperature. Objects developing radiative interiors can follow Henyey tracks. Eventually sustained core hydrogen burning balances the luminosity and the star reaches the main sequence. Massive objects can ignite hydrogen while substantial accretion continues, so “accretion ends, then hydrogen ignites” is not a universal ordering.

A useful structural distinction is therefore between an accretion-powered protostar and a largely fixed-mass pre-main-sequence contraction model. An infrared-bright envelope or disc traces obscured formation and reprocessing of radiation; a schematic evolutionary track should not replace the actual mass- and entropy-dependent stellar-structure problem.

Both neutron diffraction and X-ray diffraction use interference from atomic arrangements. Elastic scattering gives reciprocal-space peaks when the momentum transfer matches a reciprocal lattice vector. For a basis with sites ${\mathbf{r}}_j$, the structural amplitude is $F(\mathbf{Q})={\sum }_j{b}_j(\mathbf{Q})e^{i\mathbf{Q}\cdot {\mathbf{r}}_j}$, and peak intensities constrain positions, occupancies and thermal displacement as well as lattice symmetry. Powder averaging loses some directional information that a suitable single-crystal measurement retains.

X-rays couple primarily to the electron distribution. Their atomic form factors generally favor heavier atoms and decrease with increasing scattering angle; the wavelength can also be tuned near an absorption edge for element-sensitive anomalous contrast. Neutrons scatter from nuclei with scattering lengths that vary irregularly with element and isotope. This can make light atoms visible beside heavy ones, and isotope substitution can change contrast without changing the nominal chemical element. In particular, replacing hydrogen with deuterium changes the balance of coherent and incoherent neutron scattering.

Neutron magnetic scattering also couples to electronic magnetic moments, allowing magnetic order to produce peaks distinct from the nuclear structure. X-ray magnetic scattering is possible too, especially under suitable resonant conditions; it is incorrect to say X-rays can never probe magnetism. Neutron experiments often require larger samples and dedicated sources, whereas X-ray scattering is available in laboratories as well as synchrotrons. Absorption, multiple scattering, texture and instrumental resolution affect both methods. Their complementarity comes from different contrast weights and magnetic sensitivity, rather than from two identical pictures of the same density.

Two stages of the thermal history. Big Bang nucleosynthesis is nuclear assembly in an almost fully ionized, radiation-dominated plasma. Around $t\sim 1\mathrm{s}$, at temperatures of order $1\mathrm{MeV}/k_B$, the weak neutron–proton conversion rate falls below the expansion rate. The neutron-to-proton abundance freezes out at order $1/5$ and then decreases by neutron decay. The deuterium bottleneck delays effective nuclear assembly until $k_{B}T\sim 0.07$–$0.1\mathrm{MeV}$, at ages of a few minutes: the high photon-to-baryon ratio permits photodissociation despite the deuteron binding energy exceeding the mean photon energy. Once deuterium survives, reactions rapidly put most remaining neutrons into helium-4. A neutron-to-proton ratio near $1/7$ gives a helium mass fraction $Y_p\simeq 2(n/p)/(1+n/p)\simeq 1/4$, with hydrogen dominant by number and small residual deuterium, helium-3 and lithium-7 abundances. Expansion and low density prevent substantial heavy-element production. Nuclear binding is involved; neutral atoms do not yet survive.

Photon decoupling occurs much later, near a few thousand kelvin and a few hundred thousand years, after matter has overtaken radiation in the expansion dynamics. Electrons and nuclei form neutral atoms during cosmological recombination. The resulting decline of free-electron density reduces Thomson opacity, and photons begin to free stream over cosmological distances. Those photons form the observed cosmic microwave background. Atomic binding and scattering opacity, rather than new nuclear reactions, control this stage. Helium forms neutral atoms earlier than hydrogen. Recombination, photon last scattering, baryon drag release and the later thermal decoupling of the gas are related but not identical events.

The temperature falls through expansion: after the relevant entropy releases have ended, $T_{\gamma }\propto a^{-1}$ and number densities scale as $a^{-3}$. Big Bang nucleosynthesis lies in the radiation-dominated universe; photon decoupling lies after matter–radiation equality, when a matter-dominated approximation is useful but radiation is not exactly negligible.

Photon density and baryon abundance. The cosmic microwave background has a black-body radiation spectrum obeying Planck's law:where $u_{\nu }$ is spectral energy density, not a photon-number density. The measured temperature is about $2.72548\mathrm{K}$.^[ref] The mean photon energy in a black body isUsing the supplied rounded energy density givesThe exact black-body number integral at $2.725\mathrm{K}$ gives $4.1\times 10^8{\mathrm{m}}^{-3}$; the difference reflects rounding of the stated energy density. Either estimate is adequate.

A numerical photon-to-baryon ratio additionally needs a baryon density, which is not printed in this question. Make that observational input explicit rather than infer baryons from radiation alone. For illustration take the standard rounded physical density ${\Omega }_b{h}_H^2\simeq 0.022$, with $H_0=100h_H\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1}$. This is consistent with the cosmological parameter fits available around the paper's date.^[ref] It givesCounting baryons, including those bound in helium, rather than counting nuclei then givesBoth populations dilute as $a^{-3}$ after electron–positron annihilation, so this ratio also applies at hydrogen recombination in the standard adiabatic history.

Estimating decoupling. The large photon excess explains why setting $k_{B}T$ equal to the Hydrogen binding energy ${\chi }_H=13.6\mathrm{eV}$ is wrong. Even a small ionizing tail can oppose neutral-atom formation. The crude exponential balance $e^{-{\chi }_H/(k_{B}T)}\sim \eta$ gives ${\chi }_H/(k_{B}T)\sim \ln (1/\eta )\sim 21$, or $T\sim 7\times 10^3\mathrm{K}$. This only establishes the scale and the large suppression below ${\chi }_H/k_B$; it neglects phase-space factors and the opacity criterion.

A more controlled elementary Saha-Thomson decoupling estimate uses the hydrogen-only Saha equation. With $X_e=n_e/n_H$ and ground-state statistical weights combined in the usual way,The electron translational phase-space factor is significant inside the logarithm. This equilibrium relation places the sharp neutralization interval around $T\sim (3\text{–}4)\times 10^3\mathrm{K}$. Combine it with the approximate Thomson decoupling rate criterionFor $Y_p\simeq 0.25$, $n_{b0}\simeq 0.25{\mathrm{m}}^{-3}$, $H_0\simeq 68\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1}$, and ${\Omega }_m\simeq 0.31$, solving these two approximate relations gives $T\simeq 3.1\times 10^3\mathrm{K}$, with equilibrium $X_e$ of order $10^{-2}$. Thus the appropriate estimate isThis is not a precision recombination calculation. Hydrogen-only equilibrium omits the kinetic bottleneck in escaping Lyman photons, two-photon processes, the detailed helium history and freeze-out; $\Gamma /H\sim 1$ is not identical to the maximum of the last-scattering visibility function. These approximations must be improved to predict an accurate free-electron fraction or a percent-level last-scattering redshift. The high photon-to-baryon ratio alone does not fix that detailed history.

The electrons are highly nonrelativistic: $k_{B}T/(m_e{c}^2)\sim 5\times 10^{-7}$, and their thermal root-mean-square speed is only about $10^{-3}c$. Photon coupling does not make their kinetic energies relativistic.

Age–temperature relation. During a nearly flat pressureless matter-dominated phase, with radiation, curvature and vacuum energy neglected in the expansion equation,Here the Big Bang time is chosen as $t=0$. Conserved photon entropy gives $T\propto a^{-1}$, soTo normalize this relation, use the matter-era expansion rate expressed in present-day measured parameters, rather than assume the late universe is still matter dominated:The same illustrative parameters give $t_{\mathrm{dec}}\sim 4.5\text{–}5\times 10^5\mathrm{yr}$ for $T_{\mathrm{dec}}\sim 3.0\text{–}3.1\times 10^3\mathrm{K}$. Radiation was still a meaningful correction; retaining it reduces the age, and a more complete recombination calculation gives the familiar scale of roughly $4\times 10^5\mathrm{yr}$.^[ref] A numerical age cannot follow from a proportionality $t\propto T^{-3/2}$ alone: an expansion or density normalization is essential. The small amount of dark energy at decoupling can be neglected even though it affects the present age.

As the temperature falls, electrons and protons form neutral hydrogen through cosmological recombination. The free electron fraction drops sharply, and so does the scattering rate ${\Gamma }_{\gamma }=n_e{\sigma }_{T}c$. When this rate becomes comparable to or smaller than the Hubble parameter, the plasma becomes transparent and photon decoupling occurs. The Thomson scattering cross-section itself does not abruptly vanish. The rapid change is in the abundance of free targets.

The Hydrogen binding energy is $13.6\mathrm{eV}$, but recombination occurs at a much smaller thermal energy, of order $0.3\mathrm{eV}$, or a few thousand kelvin. There are roughly a billion photons per baryon, so the high-energy tail can keep hydrogen ionized long after $k_{B}T$ drops below the binding energy. The falling exponential ionizing tail then makes the transition comparatively sharp. Last scattering follows when sufficiently few electrons remain free.^[ref]

To estimate a horizon, specify an expansion model as well as the redshift. The proper particle horizon at emission isAs a transparent analytic estimate, approximate a spatially flat matter-dominated universe by $a(t)=(t/t_{\mathrm{now}}{)}^{2/3}$ throughout. ThenTaking the illustrative present value $H_0=70\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1}$, this model has $t_{\mathrm{now}}=2/(3H_0)\simeq 9.3\times 10^9\mathrm{yr}$. At $z_{\ast }=10^3$,This is a proper horizon radius at emission; its corresponding present comoving length is about $270\mathrm{Mpc}$.

For the same model the comoving distance to the Cosmic microwave background last-scattering surface and its angular diameter distance areTherefore the angular extent associated with one radial horizon length isThis is an order-one-degree estimate. Including the earlier radiation domination and late dark energy changes the coefficient; it does not make the horizon span the sky. A horizon patch diameter has twice this angle, whereas the acoustic scale uses a sound horizon and differs again. The redshift by itself, without an expansion model or normalization, is insufficient for a unique horizon distance.

The horizon problem is that the observed cosmic microwave background has nearly the same mean temperature across angular separations vastly greater than this scale. In the purely decelerating Hot Big Bang model, sufficiently separated last-scattering regions had neither time to exchange signals nor a common post-Big-Bang causal region that could establish their observed thermal uniformity. The puzzle concerns the missing causal explanation of those initial correlations; it is not an inconsistency of specifying uniform initial conditions mathematically. An earlier period of cosmic inflation can enlarge a formerly causally connected region beyond the later observable horizon.

Carl Wilhelm Scheele and Joseph Priestley independently prepared and investigated the gas now called oxygen; Antoine Lavoisier was central to interpreting it within a new account of combustion and chemical composition. Assigning a unique discoverer requires deciding whether discovery means first preparation, first public report, recognition of a distinct kind of gas, or a successful theoretical interpretation. The historical controversy is largely about why one of these achievements should be privileged over the others.

Carl Wilhelm Scheele's investigations, generally dated to the early 1770s, preceded publication of Chemical Treatise on Air and Fire by Carl Wilhelm Scheele in 1777. His work distinguished the component that sustained burning from the remainder of ordinary air. He obtained the active component through several chemical preparations, including heating suitable substances, and interpreted the results in terms of phlogiston. Delayed publication matters to public priority, but it does not show that the experimental work had not occurred. Equally, a publication date alone does not establish the precise date of every earlier trial.^[ref]

Joseph Priestley's 1775 account is unusually revealing because it records hesitation and revision. On 1 August 1774 he heated red mercury(II) oxide with a large burning lens and collected the released gas. A candle burned particularly strongly in it; later trials involving animals and comparison with nitric oxide helped him distinguish its properties from ordinary air. He initially lacked a settled understanding of what he had obtained. Calling the product dephlogisticated air eventually made it intelligible within his own framework: it had an unusually great capacity to accept the phlogiston supposedly released in burning and animal breathing.^[ref]

The modern chemical equation,is a useful reconstruction, but it must not be passed off as Joseph Priestley's own explanation. In modern terms, heating decomposes mercury(II) oxide into mercury and molecular oxygen. The eighteenth-century experimenter instead worked with such categories as calces and differently modified airs. The material operation can be recognized across the change in theory while the meaning of its product remains historically contested.

Antoine Lavoisier made a different contribution. His programme tracked what entered and left substances in calcination and combustion, placing gases within quantitative material balances. In the mercury experiment described in Elements of Chemistry by Antoine Lavoisier, a retort communicated with an enclosed supply of air over mercury. Prolonged heating formed a red calx and diminished the enclosed air. Stronger heating of the recovered calx released an especially combustion-supporting gas; adding it back to the residual air restored the properties of the original mixture. The persuasive feature was the linked decomposition and reconstruction, not just a single spectacular flame.^[ref]

That sequence supported the claim that air was a mixture and that calcination involved combining a metal with part of it. A metal could gain mass because material entered it from the air, without requiring an invisible substance to leave with anomalous properties. Conservation of mass made the weighing scale a constraint on the explanation. However, Antoine Lavoisier's entire system was not modern chemistry. He treated caloric as an explanatory principle and associated the newly named oxygen with the formation of all chemical acids. The latter universal claim is false: hydrochloric acid is an chemical acid whose solute, hydrogen chloride, contains hydrogen and chlorine but no oxygen. This concrete counterexample prevents us from defining discovery as possession of today's complete theory.^[ref]

There are consequently several legitimate allocations of credit. An investigator concerned with earliest preparation may emphasize Carl Wilhelm Scheele; one concerned with the well-documented 1774 experiment and its rapid communication may emphasize Joseph Priestley; one concerned with the explanatory reconstruction of combustion may emphasize Antoine Lavoisier. The last achievement depended on experimental results that the others had produced. Correct interpretation was neither a substitute for preparing the material nor an achievement wholly independent of communication between investigators.

Historians disagree because their accounts embody different pictures of scientific change. A simple accumulation story can make discovery look like identifying one additional item in a ready-made list of chemical elements. A conceptual-change story emphasizes that the list and its criteria were themselves changing: oxygen was not just another variety of modified air in the new chemistry. National priority narratives and the dramatic reputation of the Chemical Revolution can further encourage a single heroic winner. Yet replacing that winner with a collective label is insufficient unless the different experimental and theoretical contributions are explained.

The most defensible answer distributes credit while specifying the achievement. Carl Wilhelm Scheele and Joseph Priestley deserve independent experimental discovery credit; Antoine Lavoisier deserves central credit for the chemical interpretation and reorganization that gave oxygen its place in a new system. Discovery was a sequence of preparation, discrimination, communication and reinterpretation, rather than one instantaneous encounter with a theory-free fact.

The specific intensity of a blackbody is direction-independent over its outward hemisphere. Its outward flux per unit frequency is thereforeMultiplying by the stellar surface area gives the spectral luminosityUsing the Planck law, the hydrogen-ionizing photon production rate follows by dividing each spectral contribution by its photon energy and integrating above the threshold:With $x=h\nu /(k{T}_s)$ and $y=I_H/(k{T}_s)$, this becomes

For the initial hydrogen-only photon balance, assume an ionization-bounded nebula at steady photoionization equilibrium, negligible dust absorption and ionizing-photon leakage, negligible collisional ionization, and local reabsorption of photons from direct ground-state recombinations. A ground-state recombination then causes another photoionization and produces no net removal of a proton. Net removal is counted by Case B recombination, whose coefficient is $\beta ={\sum }_{n\ge 2}{\alpha }_n$. Integrating the local balance givesThe stellar temperature $T_s$ sets the photon supply; the gas temperature sets the recombination coefficients. They are not generally equal. For uniform hydrogen and constant gas temperature, the Strömgren radius iswith $N_e=N_H$ in fully ionized pure hydrogen.

For the hydrogen-helium mixture, define the blackbody ionizing photon functionThe three relevant thresholds are approximately $13.6$, $24.6$ and $54.4\mathrm{eV}$, for ${\mathrm{H}}^0$, ${\mathrm{He}}^0$ and ${\mathrm{He}}^{+}$ respectively. Higher stellar temperature increases the relative supply of the harder photons. These cumulative photon rates cannot be assigned independently to the three species: a photon above a helium threshold can instead ionize hydrogen, and helium recombination radiation can also ionize hydrogen. The absorption cross-sections and the diffuse radiation determine the effective allocation.

A useful sharp-front model of nested hydrogen-helium ionization zones has an inner ${\mathrm{He}}^{++}$ sphere of radius $R_2$, a ${\mathrm{He}}^{+}$ shell ending at $R_1$, and an outer hydrogen-ionized region ending at $R_H$, with $R_2\le R_1\le R_H$. Let $n_H,n_{He}$ be uniform nuclear number densities. The electron number densities in these three regions are $n_H+2n_{He}$, $n_H+n_{He}$ and $n_H$. For constant gas temperature, integrating the recombination rates givesHere the effective budgets mean net species-specific ionizations after eliminating local ground-state recycling. Their calculation must include photon competition and any interspecies diffuse contribution; the helium coefficients label the recombining ion's emitted spectrum. Once these budgets are specified, solve the third equation for $R_2^3$, the second for $R_1^3-R_2^3$, and the first for $R_H^3$.

For an isolated helium front with approximately common electron number density, the useful scaling isThus the small helium abundance partly compensates for the smaller hard-photon supply. A cooler central star produces little ${\mathrm{He}}^{++}$ and can have a substantially smaller helium-ionized zone. A sufficiently hot star can ionize helium throughout most of the hydrogen region. If independent radius estimates place a helium front beyond the hydrogen front, the independent photon allocations are inconsistent: the fronts and transfer must be treated together. The radii depend on photon hardness, abundance, electron number density and recombination rates, not on thresholds alone. The sharp-front equations are an approximation; they do not assert three exact independently conserved photon budgets for a real mixed nebula.

For a nonuniform nebula, an observed recombination line with effective recombination coefficient ${\alpha }_{\ell }^{\mathrm{eff}}$ has luminosityIf gas temperature and density permit a known approximately constant ratio $\beta /{\alpha }_{\ell }^{\mathrm{eff}}$, the integrated line givesThe unknown density distribution cancels because the same recombination integral appears in both quantities. For a stellar continuum measurement at a nonionizing frequency ${\nu }_o$,When the inferred recombination budget traces the full appropriate ionizing-photon supply, define ${\Phi }_{E_0}(T)={\int }_{E_0/h}^{\infty }{B}_{\nu }(T)/(h\nu )d\nu$ and obtain the Zanstra method equationBoth $R_s$ and the distance $d$ cancel. Hydrogen and helium recombination lines probe different hardness thresholds, so their ratios and the stellar continuum constrain $T_s$ and test the assumed spectrum and trapping. For a mixed nebula, use the corresponding transfer-corrected photon budgets. temperature gradients require emission-weighted coefficients; dust, leakage, density-bounded helium zones or departures from a blackbody can produce discrepant hydrogen and helium Zanstra estimates. Line fluxes alone without an appropriate budget or continuum normalization need not uniquely determine $T_s$.

Finally, locate the derived $T_s$ and luminosity on the Hertzsprung-Russell diagram. The central stars are hot exposed cores following post-asymptotic giant branch evolution: envelope loss creates the planetary nebula, contraction heats the core, and continued shell burning can support an approximately constant luminosity as the star moves leftward. Later the luminosity falls and the remnant approaches the white dwarf cooling sequence. Comparing central-star positions with stellar evolution tracks constrains core masses and evolutionary ages; comparison with nebular expansion ages tests whether the core heats rapidly enough to ionize the expelled envelope. The central-star locus traces the transition from an envelope-losing giant to a compact stellar remnant.

For instantaneous complete hydrogen reionization, take $n_{e0}={\Omega }_b{\rho }_{\mathrm{crit},0}/(m_p{c}^2)=3{\Omega }_b{H}_0^2/(8\pi G{m}_p)$ and $n_e(z)=n_{e0}(1+z{)}^3$. Helium or incomplete ionization would multiply this by the appropriate free-electron count per baryon. In the matter-plus-curvature model,The reionization optical depth in an open matter universe is thereforeSet $u=\sqrt{1+{\Omega }_{m}z}$ and ${\beta }^2=1-{\Omega }_m$. The integral is $2{\Omega }_m^{-2}[u^3/3-{\beta }^{2}u{]}_1^{u_{\mathrm{re}}}$, givingAt high redshift where matter dominates over curvature, this approaches $\kappa \simeq {\Omega }_b{C}_H(1+z_{\mathrm{re}}{)}^{3/2}/\sqrt{{\Omega }_m}$. The exact expression should be retained when ${\Omega }_m{z}_{\mathrm{re}}$ is not large.

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.

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).

The event rate counts unordered pairs of reacting nuclei. For unlike species there are $n_i{n}_j$ pairs per unit-volume normalization; for identical species there are approximately $n_i^2/2$ pairs, since counting ordered pairs counts each collision twice. The Kronecker delta therefore supplies exactly the required divisor $1+{\delta }_{ij}$. It is the event rate, not the destruction rate of identical nuclei: one identical-pair reaction destroys two nuclei.

The temperature dependence combines the high-energy tail of the Maxwell-Boltzmann distribution with quantum tunnelling through the Coulomb barrier. The principal exponential factor in the energy integral is $\exp [-E/(k_{B}T)-b/\sqrt{E}]$. Its stationary point, the Gamow peak, has $E_0\propto T^{2/3}$ and exponent $\eta =3E_0/(k_{B}T)\propto T^{-1/3}$. Integrating around this peak gives a prefactor proportional to $T^{-2/3}$, or ${\eta }^2$, when the nuclear factor varies slowly. Thus the supplied thermonuclear reaction rate is locally proportional to ${\eta }^2{e}^{-\eta }$. The slow weak interaction in the proton-proton reaction also makes its overall coefficient very small; its local thermal slope follows the supplied formula.

At fixed number densities, define the local temperature exponent of a thermonuclear reaction by a logarithmic derivative. Since $d\eta /d\log T=-\eta /3$,For proton-proton reactions, $A=1/2$ and $Z_1=1$. At $T_6=15$, ${\eta }_{11}=13.67$ and $\alpha =3.89$. For two helium-3 nuclei, $A=3/2$ and $Z_3^2{Z}_3^2=16$, giving ${\eta }_{33}=49.68$ and $\beta =15.89$. For helium-3 with helium-4, $A=12/7$ and the same charge product gives ${\eta }_{34}=51.95$ and $\gamma =16.65$. HenceThese are local powers near the specified temperature, not constant exponents over an arbitrarily large range. The greater reduced nuclear mass in the $34$ reaction increases its tunnelling exponent and its thermal sensitivity.

For the effective proton-proton reaction network, denote the event rates byThe short-lived beryllium intermediate can be eliminated, so production through $R_{34}$ feeds $n_7$, representing predominantly lithium-7. This rapid conversion is electron capture. Count the nuclei consumed and produced in each reaction: $11$ consumes two protons and makes one deuterium; $21$ consumes a proton and a deuterium and makes one helium-3; $33$ consumes two helium-3 and makes one helium-4 and two protons; $34$ consumes one of each helium isotope; $17$ consumes a proton and lithium-7 and makes two helium-4. It follows thatFor closure, ${\dot{n}}_7={\lambda }_{34}{n}_3{n}_4-{\lambda }_{17}{n}_1{n}_7$. These nuclear terms conserve $n_1+2n_2+3n_3+4n_4+7n_7$, providing an independent stoichiometric check. They describe a fixed local volume; compression or mixing would add transport terms.

The roughly one-second destruction time of deuterium is negligible compared with the proton-burning and helium-3 evolution times. Its abundance rapidly adjusts to the slowly evolving production rate. Set ${\dot{n}}_2\simeq 0$, giving the deuterium quasi-equilibrium in proton-proton burningSubstitution, without setting ${\dot{n}}_1=0$, givesThe helium-3 relaxation time of about $6\times 10^5$ years is likewise short compared with the central proton-burning time and the age of the Sun. On this relaxation time the proton and helium-4 abundances, and the background temperature, can be treated as fixed. Put $A={\lambda }_{33}$, $B={\lambda }_{34}{n}_4$, $C={\lambda }_{11}{n}_1^2/2$. Then ${\dot{n}}_3=C-A{n}_3^2-B{n}_3$, and its nonnegative equilibrium isThe other root is negative and unphysical. Write $n_3=n_{3e}+x$ and subtract the equilibrium equation. Exactly, for a fixed background,Dropping the quadratic displacement givesThus small displacements decay exponentially. A slowly moving equilibrium introduces an additional $-{\dot{n}}_{3e}$ term, negligible to leading order when $\tau$ is short compared with background evolution.

To estimate the temperature for helium-3 freeze-out, use the cooler pp-I-dominated region and neglect $B$ compared with $2A{n}_{3e}$. Then $n_{3e}\simeq n_1\sqrt{{\lambda }_{11}/(2{\lambda }_{33})}$ and ${\tau }^{-1}\simeq n_1\sqrt{2{\lambda }_{11}{\lambda }_{33}}$. With the local powers just obtained and fixed proton number density, $n_{3e}\propto T^{-6}$ and $\tau \propto T^{-10}$. Adopting $t_{\odot }\simeq 4.6\times 10^9\mathrm{yr}$ givesThis is an order-of-magnitude extrapolation of local exponents. Keeping the supplied exponential factors gives the helium-3 freeze-out estimate with Gamow factorswhich yields about $6.8\times 10^6\mathrm{K}$ with the same normalization and fixed-density approximation. Thus the robust estimate is several million kelvin. The supplied data do not fix the local density, composition or pp-II fraction away from the centre, so they cannot determine an exact transition temperature or radius. In particular, using $\tau \propto T^{-16}$ alone would omit the temperature dependence of the equilibrium helium-3 abundance.