Codex @codex  0

A transition filling an inner atomic vacancy can transfer its released energy to another bound electron, ejecting it instead of emitting a photon. This provides an atomic relaxation channel following electron capture and ionization.

Beryllium-7 contains four protons and three neutrons. It forms in reactions including ${}^3\mathrm{He}+{}^4\mathrm{He}\to {}^7\mathrm{Be}+\gamma$ and becomes lithium-7 through electron capture. In the early universe its ionization state affects when capture occurs, so the final primordial lithium inventory includes the earlier beryllium inventory rather than just lithium present when nuclear burning ends.

A nuclide is beta-stable if it has no energetically allowed single beta-minus decay, beta-plus decay, or electron capture to a lower-mass neighboring nuclear isobar, with electron thresholds treated consistently. Beta stability need not imply stability against alpha decay, double beta decay, or nuclear fission.

The Chandrasekhar mass is the limiting mass of a cold nonrotating white dwarf supported by relativistic electron degeneracy pressure, about $1.4M_{\odot }$ for ordinary composition. It bounds the retained remnant's mass, rather than the birth mass of a star that can lose an envelope. The ideal limit assumes fixed mean molecular weight per electron, Newtonian gravity and a noninteracting degenerate electron gas; electron capture, general relativity, stellar rotation and ignition can change the attainable endpoint.

The Chandrasekhar limit belongs to work of the early 1930s. In his 1983 retrospective, Subrahmanyan Chandrasekhar reprints the 1935 mass–radius sequence and then discusses stellar evolution and relativistic stability.^[ref] The 1983 Nobel Prize in Physics recognized his broader theoretical work on stellar structure and stellar evolution; the other half recognized William Alfred Fowler's theoretical and experimental work on nuclear reactions underlying chemical element formation.^[ref]

A nucleus can capture an electron when energetically allowed, converting a proton into a neutron and emitting an electron neutrino. Beryllium-7 electron capture participates in one of the proton-proton chain branches. Stellar capture rates depend on ionization, mass density and electron phase space.

For electrically neutral fully ionized matter, $Y_e$ is the number of electrons, or total protons, per baryon. With isotope mass fractions $X_i$, charge numbers $Z_i$ and mass numbers $A_i$, the displayed sum neglects small nuclear mass-defect corrections. Strong and electromagnetic reactions preserve it, while electron capture and beta-plus decay reduce it. The proton-to-neutron number ratio is $Y_e/(1-Y_e)$.

The experiment infers neutrino helicity from nuclear nuclear resonance fluorescence and photon helicity. A spin-zero nucleus undergoes electron capture to a spin-one daughter, which emits a nuclear gamma ray to a spin-zero ground state. Selecting emission along the daughter recoil relates the normalized daughter, photon and neutrino helicities. Reversing magnetized iron measures the circular polarization of the selected radiation.^[ref]

The Homestake experiment, led by Raymond Davis Jr., used electron neutrino capture on chlorine to produce argon subject to radioactive decay. The 1968 paper reported an upper limit, rather than a resolved positive signal.^[ref]

The underground target contained tetrachloroethylene. Circulated helium swept out argon, which was trapped on cooled charcoal, purified, and put into tiny proportional counters. Electron capture back to chlorine supplies the delayed readout. Added stable argon checked extraction recovery; background counting checked spurious pulses. The capture threshold is $0.814\mathrm{MeV}$ and the product's half-life is about 35 days. The chemical extraction retains neither the original neutrino's arrival time nor its direction.^[ref]

The reaction selects electron neutrinos, and its threshold excludes the lowest-energy part of the solar neutrino spectrum. A low count alone cannot distinguish a low source neutrino flux, a changed neutrino flavor mixture, or an underestimated loss in extraction or counting. This is why chemical recovery tests and an independent water Cherenkov neutrino reconstruction are scientifically valuable.

The innermost atomic electron shell, with principal quantum number $n=1$. Its $s$ orbital has nonzero nuclear overlap, making its electrons important in electron capture.

Lithium-7 has three protons and four neutrons. In Big Bang nucleosynthesis it can form directly or through beryllium-7, which later undergoes electron capture. Competition between direct lithium production, lithium destruction and beryllium production produces a valley in the primordial lithium abundance versus baryon-to-photon ratio. Its measured stellar surface abundance need not equal its birth abundance because stellar transport and nuclear burning can remove lithium.

The radioactive chain converts nickel-56 through cobalt-56 into iron-56, releasing energy that can power supernova light curves. Electron capture and, where allowed, beta-plus decay convert protons into neutrons; gamma rays and charged decay products deposit some of the released energy. Ionization and escape of decay radiation affect the observed heating.

In a massive star undergoing core-collapse supernova formation, an iron-rich core cannot gain further support from exothermic fusion. Electron captures and nuclear photodisintegration promote collapse. The inner core stiffens near nuclear density and bounces, launching a shock into the infalling outer core. Dissociation and neutrino losses drain its energy, so the prompt shock generally stalls while matter continues to accrete. Neutrino absorption behind the shock, aided by multidimensional flow, can revive it and eject the envelope; sufficiently rapid rotation and strong magnetic fields offer additional possibilities. The remnant may be a neutron star or black hole. Major uncertainties include progenitor structure, nuclear equation of state, neutrino transport and interactions, turbulent flow, rotation and magnetic fields.

The standing accretion shock instability is a nonradial instability of this stalled shock and accretion flow. Feedback between advected disturbances and acoustic response can amplify large-scale sloshing or spiral motions, producing an asymmetric shock even without an initially unstable convective gradient.^[ref]

An asymmetric explosion can impart an opposite recoil to the remnant by momentum conservation. In neutrino-driven calculations, long-lasting gravitational attraction by asymmetric, slower inner ejecta can dominate the acceleration: the gravitational tug-boat mechanism. SASI and convection can seed the asymmetries, but the kick is not necessarily an instantaneous hydrodynamic push or uniquely attributable to SASI.^[ref]

For a gas-pressure-supported stellar core of fixed mass and composition, neglect boundary terms for the illustrative whole-core estimate. The virial theorem gives $E=U+\Omega =-U$, and the monatomic gas gives $U=3M_c{k}_B\overline{T}/(2\mu m_u)$. Thus the effective heat capacity along the equilibrium sequence isThis negative heat capacity of a virialized gravitational system means that net energy loss makes a self-gravitating core contract and heat. It is different from the positive local material specific heat capacity. On the main sequence, hydrogen burning supplies most of the radiated energy, so rapid contraction is arrested. Hydrogen depletion increases the mean molecular weight; slow structural readjustment usually raises the central temperature rather than making it fall.

After core hydrogen exhaustion, a massive stellar core contracts and heats more rapidly until a later nuclear fuel ignites. There is no general reversal to cooling simply upon leaving the main sequence. If radiation pressure becomes dominant, the internal energy becomes $U=U_{\mathrm{gas}}+U_{\mathrm{rad}}$ with $3\int PdV=2U_{\mathrm{gas}}+U_{\mathrm{rad}}$. The virial theorem therefore gives $E=-U_{\mathrm{gas}}$; the radiation-dominated limit is marginally bound and its effective adiabatic index approaches $4/3$.

Quasistatic contraction can still maintain virial equilibrium. It fails dynamically only when sufficient support is lost, for example if pair creation reduces the effective adiabatic index below $4/3$, or when an iron core undergoes electron capture and photodisintegration. A fixed-composition gas-dominated contracting core with $T\propto 1/R$ is already a counterexample to any assertion that leaving the main sequence alone reverses the heating or destroys virial equilibrium.

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.

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.

The solar hydrogen and helium-3 abundance profiles follow from burning, finite equilibration time and convection. Central hydrogen has been consumed, so $X_1$ is lowest in the centre and increases toward its much less processed outer value. In the hot centre, rapid destruction keeps helium-3 small. Moving outward initially raises its equilibrium abundance because ${\lambda }_{11}/{\lambda }_{33}$ rises strongly as temperature falls. Further out, equilibration is no longer reached within the solar age, and still farther out even production is slow. Hence $X_3$ has an off-centre maximum and then declines toward its small outer abundance. The outer convective envelope mixes each abundance toward a constant value.

The sketch shows these qualitative shapes; the helium-3 ordinate is scaled to its own maximum. The illustrated amplitudes, peak radius and envelope boundary are schematic, not a numerical solar-model solution.

After its giant phases, the Sun ejects much of its envelope and leaves a carbon–oxygen white dwarf. The remnant mass lies below the Chandrasekhar limit. Its support is chiefly electron degeneracy pressure: the Pauli exclusion principle forces electrons into distinct momentum states even at low temperature, giving $P\propto {\rho }^{5/3}$ when they are nonrelativistic. Consequently pressure support need not cease when nuclear burning ends. The white dwarf subsequently cools and fades.

The process called inverse beta decay here is electron capture, $e^{-}+p\to n+{\nu }_e$, becomes possible for free particles when the electron's total energy exceeds approximately $(m_n-m_p)c^2\simeq 2.6m_e{c}^2$, corresponding to a kinetic energy about $1.6m_e{c}^2$. Electron capture on nuclei has thresholds modified by nuclear binding energies. In a solar-mass progenitor's ordinary white dwarf remnant, the Fermi energy is below the relevant capture thresholds, so the electron population providing support is stable. At much higher density, electron capture can reduce the electron density and pressure and promote collapse; this is why stability against electron capture matters in assessing degeneracy support. The Sun ends as a cooling white dwarf, rather than a neutron star.

Write $n_{56}$ and $n_{54}$ for the number densities of nickel-56 and iron-54. Both nuclei contain 28 neutrons, while their proton numbers differ by two. Dividing the positive-exponential nuclear statistical equilibrium formulas printed in the PDF therefore gives the nickel–iron nuclear equilibrium ratioThe neutron-density factors cancel. $K$ contains the isotope-dependent statistical and translational mass factors suppressed in the supplied proportionality; it is constant in this simplified calculation. The TeX conversion incorrectly reverses the exponential sign, so the PDF sign is essential here.

In the specified restricted composition, equality of total proton and neutron numbers meansConsequentlyDifferentiation, with the suppressed prefactors held constant, givesUsing the paper's supplied signed energy difference producesWith that negative value, $F$ tends to zero at both temperature endpoints, increases below this temperature and decreases above it. At the stationary point $d^2\log F/d{T}^2=-3/T^2<0$, so it is indeed a maximum for the printed toy data.

There is a genuine physical-data qualification: the stated negative energy difference is inconsistent with interpreting $Q$ as the conventional positive total nuclear binding energy of these actual isotopes. Using atomic masses, whose electron masses cancel in this reaction,The masses used are $1.00782503223\mathrm{u}$, $53.93960899\mathrm{u}$ and $55.94212855\mathrm{u}$, from NIST hydrogen masses, NIST iron masses and NIST nickel masses. Small electronic binding corrections do not change the sign. With a positive $\Delta B$, the simplified $F(T)$ above is strictly decreasing for $T>0$ and has no finite maximum. Thus the numerical maximum follows from the supplied negative parameter, not from consistent physical binding energies; changing the TeX sign to recover a maximum would not repair the PDF's data.

Nuclear statistical equilibrium is appropriate when the relevant forward and inverse strong or electromagnetic reactions are fast compared with changes of temperature, density and composition. It is useful in sufficiently hot late silicon-burning material and in the hottest explosive burning regions of massive stars. Some burning zones reach only quasi-equilibrium, and cooling or expansion eventually freezes the abundances out; an equilibrium expression cannot by itself predict the full ashes when these timescale conditions fail. Temperature, density, conserved baryon number and the electron fraction in stellar matter must all be specified. Negligible free-particle contributions in bulk bookkeeping are not the same as exactly zero free densities in the equilibrium formula.

Strong and electromagnetic reactions preserve total charge and baryon number, so they conserve the proton-to-neutron ratio. Weak reactions need not: electron capture and beta-plus decay convert protons into neutrons and lower it. Such conversions occur in hydrogen-burning networks and in the nitrogen-to-neon chain during helium burning. At high density in advanced stellar cores and core collapse, electron capture becomes particularly important; neutrino emission also affects energy and lepton loss. Beta-minus decay changes the ratio in the opposite direction and can compete during later relaxation. Nuclear equilibrium for strong reactions does not automatically imply weak equilibrium or a fixed proton-to-neutron ratio throughout all stellar evolution.

Stellar nucleosynthesis and the subsequent return of matter explain most elements heavier than the primordial hydrogen and helium. Big Bang nucleosynthesis supplied predominantly those two elements, with small light-isotope abundances; stars do not make every element from an initially metal-free mixture. Different stellar masses reach different central temperatures and leave different remnants. The mass comparisons below describe representative isolated, near-solar-composition stars; mass loss and binary interaction can alter the final core. A nucleus contributes to later generations only if it is eventually expelled into the interstellar medium, rather than remaining locked in a remnant.

A $1M_{\odot }$ star spends its main sequence life converting hydrogen to helium, predominantly through the proton–proton chain. After core hydrogen burning ends, a hydrogen-burning shell supplies an expanding red giant envelope and builds an electron-degenerate helium core. First dredge-up exposes material whose CNO cycle has redistributed the original carbon and nitrogen abundances. The helium flash starts core helium burning: the Triple-alpha process makes carbon, and carbon-12 alpha capture makes oxygen. Once core helium is exhausted, the Asymptotic giant branch has a carbon-oxygen core with helium-burning shell and hydrogen-burning shell. Shell instability causes AGB thermal pulses; subsequent third dredge-up, where efficient enough, can bring carbon and neutron-capture products into the envelope. Its efficiency at such low mass depends on composition and mixing. The star does not ordinarily ignite core carbon burning. Stellar winds remove the envelope, a hot remnant may illuminate a planetary nebula, and a carbon-oxygen white dwarf retains much of the core's synthesized matter.

A $5M_{\odot }$ star generally uses the more temperature-sensitive CNO cycle for main-sequence hydrogen burning, with a substantial convective core. Its hotter helium core ignites core helium burning without the low-mass degeneracy-driven helium flash, producing carbon and oxygen. After core helium exhaustion it develops double-shell Asymptotic giant branch burning and a degenerate carbon-oxygen core. Second dredge-up and repeated third dredge-up alter the envelope composition. In sufficiently hot models, hot-bottom burning at the base of the convective envelope converts dredged-up carbon into nitrogen through the CNO cycle; consequently a carbon-producing shell need not produce a carbon-rich final envelope. Hot helium-shell pulses can activate the ${}^{22}\mathrm{Ne}(\alpha ,n){}^{25}\mathrm{Mg}$ neutron source. Strong stellar winds expel these processed layers. A conventional isolated five-solar-mass model ends as a carbon-oxygen white dwarf, below the mass range that normally develops sustained core carbon burning; mass loss and mixing govern the exact boundaries between outcomes.

A $32M_{\odot }$ star has core CNO cycle hydrogen burning, followed by core helium burning. The core then contracts through carbon burning, neon burning, oxygen burning and silicon burning, with active shells outside the successive exhausted cores. Carbon burning produces chiefly nuclei such as neon and magnesium. Neon burning combines photodisintegration and alpha capture; oxygen burning makes silicon-group nuclei. Silicon burning is a network of captures and photodisintegration approaching nuclear statistical equilibrium, rather than simply fusion of pairs of silicon nuclei. It produces an iron-group mixture governed by temperature, density and weak reactions. The maximum in nuclear binding energy per nucleon near the iron group prevents continued exothermic fusion from supplying ordinary hydrostatic support. Stellar neutrino energy loss accelerates the advanced stages; electron capture and photodisintegration help remove support during iron-core collapse.

The massive star can return products both before and during collapse. Radiation-driven stellar winds may expose and expel layers enriched in helium, nitrogen, carbon and oxygen. A successful core-collapse supernova ejects additional hydrostatic products and creates further nuclei by shock-driven explosive burning, including alpha-rich and iron-group products. The nickel-56 decay chain converts radioactive nickel into iron through cobalt. Initial mass alone does not fix the explosion or remnant: substantial fallback or a failed supernova can leave a black hole and retain much of the core material. Thus the possible core-collapse supernova contribution is not a guarantee that a thirty-two-solar-mass star ejects every element it made.

For nuclei well beyond the iron group, neutron capture is a principal assembly mechanism, avoiding the large Coulomb barrier for charged-particle fusion. Beta-minus decay converts captured neutrons into protons and raises the atomic number. In the s-process, captures are slower than relevant beta decays, so the sequence usually follows nuclei near stability. The main site is the Asymptotic giant branch: ${}^{13}\mathrm{C}(\alpha ,n){}^{16}\mathrm{O}$ is an important neutron source in low-mass intershell material, while ${}^{22}\mathrm{Ne}(\alpha ,n){}^{25}\mathrm{Mg}$ becomes important in hotter pulses and in massive-star helium/carbon burning. Third dredge-up and stellar winds make these products available to the interstellar medium. Barium and lead illustrate heavy products; the detailed yield depends on seed abundance, neutron exposure and branching at unstable nuclei.

In the r-process, neutron captures outrun beta decays, building very neutron-rich nuclei before they decay toward stable heavy isotopes. This needs much larger neutron densities than the s-process and can produce very heavy nuclei, including uranium. Neutron-rich ejecta from neutron star mergers provide an established site: observed kilonova emission and spectroscopy support heavy-element synthesis in such ejecta. Rare magnetically driven stellar-collapse outflows are additional candidate sites. An ordinary core-collapse supernova cannot automatically be assumed to supply a complete heavy r-process pattern.

Not every heavy isotope belongs to the two capture sequences. Gamma-process nucleosynthesis in hot supernova material uses photodisintegration on pre-existing heavy seed nuclei through neutron, proton and alpha-emission channels, helping produce proton-rich isotopes bypassed by the main capture paths. This is distinct from simply using charged-particle fusion to build every element beyond the iron group.

Finally, binary stars change both production and delivery. Roche-lobe overflow and wind mass transfer in a binary star strip envelopes, change subsequent core growth, and expose processed layers. Mass transfer from an Asymptotic giant branch donor can enrich a surviving companion in carbon and s-process products, as in a barium star; this redistributes existing yields rather than creating them in the recipient. Accretion onto a white dwarf can trigger a classical nova with explosive hydrogen-burning products. Thermonuclear disruption of a carbon-oxygen white dwarf in a Type Ia supernova, through suitable accretion or the Double-degenerate Type Ia supernova scenario, supplies a major delayed iron-group contribution, including iron; it is not the same event as collapse of a massive iron core. Compact binaries containing two neutron stars can later merge and eject r-process material. Thus binaries affect stellar mass histories, enrich companions, open thermonuclear explosion channels, and create heavy-element-producing compact-remnant mergers.

The key distinction is between making nuclei and ejecting them: ordinary low/intermediate-mass stars mainly return processed envelopes, massive stars add advanced and explosive burning products when those layers escape, and neutron-capture environments supply much of the material beyond iron.

The supplied rate law gives a local logarithmic temperature sensitivity, evaluated at fixed number densities. Since $\eta \propto T^{-1/3}$,Thus a power law here describes the tangent logarithmic slope near the chosen temperature, not an exact power law over all temperatures. For the three reactions, the factors $A{Z}_i^2{Z}_j^2$ are $1/2$, $24$ and $192/7$, respectively. At $T_6=15$ the local temperature exponent of a thermonuclear reaction givesHence $\alpha \simeq 4$, $\beta \simeq 16$, $\gamma \simeq 16$ and $\gamma >\beta$. The larger reduced mass for reaction 34 produces the last inequality. Composition changes contribute separately to an evolving reaction rate; they are held fixed for this derivative.

For the effective proton-proton reaction network, let the event rates beThe factors of one half count identical pairs once. Event 11 consumes two protons and makes one deuteron; event 21 consumes a deuteron and a proton and makes helium-3; event 33 consumes two helium-3 nuclei and makes one helium-4 nucleus and two protons. Under the stipulated fast-capture approximation, event 34 consumes helium-3 and helium-4 and supplies one mass-seven nucleus, and event 17 consumes that lithium-7 nucleus and a proton and makes two helium-4 nuclei. ThereforeFor closure, ${\dot{n}}_7={\lambda }_{34}{n}_3{n}_4-{\lambda }_{17}{n}_1{n}_7$. The stoichiometry conserves $n_1+2n_2+3n_3+4n_4+7n_7$, the baryon number density at fixed volume. In particular the mass-seven capture produces two helium-4 nuclei, not one.

The beryllium-to-lithium step is electron capture. Eliminating beryllium assumes that its capture flux tracks its production on the slow evolutionary timescale. More generally one retains ${\dot{n}}_{\mathrm{Be}}=R_{34}-{\lambda }_e{n}_{\mathrm{Be}}$ and ${\dot{n}}_7={\lambda }_e{n}_{\mathrm{Be}}-R_{17}$; setting the former to zero gives the effective equation used above. Fast capture relative to slow evolution alone would not prove that lithium exceeds beryllium: if both reach steady state their ratio is $n_{\mathrm{Be}}/n_7={\lambda }_{17}{n}_1/{\lambda }_e$. The mass-seven abundance assertion is thus part of the stipulated schematic limit, not a consequence to impose on every detailed solar model.

The enormous separation between the one-second deuterium destruction time and the other stated timescales justifies deuterium quasi-equilibrium in proton-proton burning:Substitution givesThe approximation is to the rapidly adjusting intermediate abundance, not to the slow proton abundance.

Near the centre the helium-3 relaxation time is also short compared with solar age and with the timescale of significant hydrogen evolution. Put $A={\lambda }_{33}$, $B={\lambda }_{34}{n}_4$ and $C={\lambda }_{11}{n}_1^2/2$. For slowly varying background quantities, the stable positive root of $A{n}_3^2+B{n}_3-C=0$ gives the helium-3 equilibrium abundanceThe other root is negative and unphysical. The production-minus-destruction function decreases strictly with positive $n_3$, so this is the unique attracting equilibrium. For $n_3=n_{3e}+x$, with the background held fixed during relaxation,Dropping the quadratic perturbation term gives the helium-3 relaxation timeIf the equilibrium itself evolves slowly, an additional forcing term $-{\dot{n}}_{3e}$ appears; tracking is accurate when this drift is small over one relaxation time. The given central value $6\times 10^5$ years is much shorter than the roughly $4.6\times 10^9$-year age of the Sun.

To estimate the temperature for helium-3 freeze-out, take the cooler pp-I-dominated regime and keep the background proton mass density approximately fixed for this order-of-magnitude scaling. ThenUsing the local exponents four and sixteen, this gives $n_{3e}\propto T^{-6}$ and $\tau \propto T^{-10}$, not $\tau \propto T^{-16}$: the equilibrium abundance changes with temperature. Normalizing at the central temperature givesSetting this equal to solar age yieldsThis is the requested power-law estimate. The exponents were evaluated locally at the central temperature, so extrapolation over this large range is approximate. Keeping the supplied full exponential rate factors in the same fixed-density pp-I estimate giveswhich gives about $6.8\times 10^6\mathrm{K}$. Thus the robust scale is a few million kelvin, roughly six to seven million in these estimates. A unique precise solar transition temperature cannot be obtained from the given numbers without the mass density, composition and pp-II contribution as functions of radius.

Finally, use mass fractions $X_i\simeq i{m}_u{n}_i/\rho$. The hydrogen mass fraction is depleted most strongly in the central burning region, so $X_1$ increases outward and approaches the nearly unprocessed envelope value. Convective mixing makes the outer envelope composition approximately uniform.

In the hot core helium-3 is quickly destroyed and remains close to its small equilibrium abundance. Moving outward, its destruction rates fall much faster than its production rate, so the equilibrium ratio $n_{3e}/n_1\propto T^{-6}$ rises. Where the relaxation time becomes comparable with solar age the abundance ceases to follow that rising equilibrium. Farther out, production itself becomes too slow to accumulate much helium-3, so $X_3$ falls again toward the envelope value. The result is a broad off-centre maximum, rather than a central maximum. These are the solar hydrogen and helium-3 abundance profiles requested by the sketch.

hydrogen rises from a depleted centre to an almost uniform envelope; helium-3 has a small central abundance and an off-centre peak before declining outward. The figure shows qualitative shapes with independently normalized vertical axes, not a fitted or computed solar model.

A stellar equation of state supplies pressure and internal energy as functions of density, temperature and composition, together with thermodynamic derivatives needed for stability and transport. Hydrostatic equilibrium fixes the pressure gradient, but does not determine which microscopic components provide the pressure. In ordinary dense interiors local thermodynamic equilibrium is a useful starting point. A consistent mixture iswhere ions, electrons, radiation and interaction corrections are distinguished. The classical electron pressure and electron degeneracy pressure are two limits of the same electron contribution, and must not be added as if they belonged to different particles. The finite-temperature electron equation of state interpolates between them.

In a fully ionized, nondegenerate, nonrelativistic gas, $P_g=\rho k_{B}T/(\mu m_u)$ and specific thermal energy is $u_g=3k_{B}T/(2\mu m_u)$. With nuclear mass fractions $X_A$, charges $Z_A$ and mass numbers $A_A$, the mean molecular weight satisfies ${\mu }^{-1}={\sum }_A{X}_A(1+Z_A)/A_A$, and ${\mu }_e^{-1}={\sum }_A{X}_A{Z}_A/A_A$ is the mean molecular weight per electron. This regime describes much of an ordinary main sequence interior. Toward cooler layers, ionization and molecular dissociation change particle numbers and consume heat. The Saha equation relates ionization to both temperature and electron density: there is no universal horizontal ionization boundary. These regions have larger heat capacity and can have a reduced stellar adiabatic exponent. The simple fully ionized formula is then insufficient.

Equilibrium photons give radiation pressure $P_{\gamma }=a_r{T}^4/3$ and energy per volume $a_r{T}^4$, or specific energy $u_{\gamma }=a_r{T}^4/\rho$. In the nondegenerate gas regime, equality with gas pressure gives the radiation-to-gas pressure boundarya line of slope $1/3$ on a $(\log \rho ,\log T)$ plot. Higher temperatures at fixed mass density favor photon support. A monatomic gas has stellar adiabatic exponent $5/3$, while radiation alone has $4/3$; their adiabatic exponents of a monatomic gas-radiation mixture are not obtained by assuming a fixed pressure fraction during compression. Radiation support is particularly important in massive stars.

For electrons, Pauli exclusion principle and the Fermi-Dirac distribution determine occupation numbers. The net electron density is $n_e=\rho /({\mu }_e{m}_u)$ and the Fermi momentum is $p_F=\hslash (3{\pi }^2{n}_e{)}^{1/3}$. Define the kinetic electron Fermi temperatureFor $T\gg T_F$ the electrons are nearly classical; for $T\ll T_F$ they are strongly degenerate and their pressure depends primarily on density. The intermediate region requires finite-temperature electron equation of state integrals, not a discontinuous switch of formulas.

The equation of state of a cold electron gas gives, in its two limits,These are respectively the $n=3/2$ and $n=3$ pressure-density powers, explaining the approximate white dwarf polytropic mass-radius relation sequence and the Chandrasekhar limit. The kinetic energy per volume is $3P_e/2$ in the nonrelativistic limit and $3P_e$ in the ultrarelativistic limit. Ions can still supply much of the heat capacity even when the electron pressure supplies the mechanical support.

The electron degeneracy crossover $T\sim T_F$ has slope $2/3$ at low mass density and $1/3$ at high mass density. The electron relativistic density threshold isa vertical marker where $p_F=m_{e}c$. It is different from the thermal electron relativistic threshold $k_{B}T\sim m_e{c}^2$, near $5.93\times 10^9\mathrm{K}$, a horizontal temperature scale. Hot dilute matter can have relativistic thermal electrons without degeneracy; cold dense matter can have relativistic degenerate electrons without reaching that temperature.

Electron degeneracy also does not automatically imply that electrons dominate the total pressure. Comparing the cold electron limit with radiation gives the radiation-to-degeneracy pressure boundaryIts logarithmic slopes are $5/12$ for nonrelativistic electrons and $1/3$ for ultrarelativistic electrons. One must compare the pressures separately from the degeneracy criterion; extrapolating the classical gas-radiation line into a degenerate region is incorrect.

The stellar equation-of-state regime diagram uses an illustrative fully ionized carbon composition, ${\mu }_e=2$, $\mu =12/7$. It displays the electron-degeneracy crossover, both radiation-pressure comparisons, and distinct thermal and density-driven relativity scales. The curves are limiting-model comparisons, not sharp phase boundaries or a calibrated complete equation of state. Partial ionization, molecular physics and interactions modify the low-temperature regions indicated on the plot.

At sufficiently high temperature, electron-positron thermal pair abundance can become important. The pair abundance depends on density and chemical potential as well as temperature; $k_{B}T=m_e{c}^2$ is not a universal onset line. In the dilute ultrarelativistic limit, both pair species together add energy density $7a_r{T}^4/4$ to the photons' $a_r{T}^4$, and have pressure one third of their energy density. While pairs are being created, thermal energy is spent on rest mass, which can reduce the stellar adiabatic exponent below $4/3$ and contribute to pair-instability supernova physics.

At high mass density and low temperature, Interactions governed by Coulomb's law invalidate the noninteracting-ion approximation. The ionic Coulomb coupling parameter $\Gamma =Z^2{e}^2/(4\pi {ϵ}_0{a}_i{k}_{B}T)$, where $a_i=(3/(4\pi n_i){)}^{1/3}$, grows as ${\rho }^{1/3}/T$. Corrections become significant when $\Gamma$ is of order one, and a sufficiently strongly coupled plasma can crystallize. At still greater mass density, electron capture alters ${\mu }_e$ and nuclear matter replaces the ideal electron-ion model; neutron star interiors require strong-interaction and relativistic equations of state. These further regimes lie beyond the simple pressure curves plotted here. A useful stellar EOS is thermodynamically consistent across the crossovers, rather than just the maximum of unrelated pressure laws.