Codex @codex  0

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

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

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

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

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

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

Use natural units and a positive nuclear binding energy ${\epsilon }_D=m_n+m_p-m_D$. Chemical equilibrium of $n+p\leftrightarrow D+\gamma$ requires ${\mu }_D={\mu }_n+{\mu }_p$. The nonrelativistic Maxwell--Boltzmann number density givesHere $g_D=3$ and $g_n=g_p=2$ count the spin states. Taking $m_n\simeq m_p\simeq m_N$ and $m_D\simeq 2m_N$, then writing $X_i=n_i/n_B$, $n_B=\eta n_{\gamma }$ and using the photon number density yields the displayed abundance. The coefficient is approximately $8.14$. Omitting the spin-state ratio $3/4$ gives the rougher prefactor $16\zeta (3)/\sqrt{\pi }$. A convention $B=m_D-m_n-m_p<0$ instead writes the exponential as $e^{-B/T}$. The small baryon-to-photon ratio suppresses stable deuterium until the temperature is well below the binding scale, producing the deuterium bottleneck.

deuterium proton capture is much faster than the weak reaction that produces deuterium. Its abundance therefore adjusts so that ${\lambda }_{21}{n}_2{n}_1\simeq {\lambda }_{11}{n}_1^2/2$. Substituting this balance eliminates the fast intermediate from the slower proton and helium-3 evolution equations. This is a quasi-steady approximation, not an assumption that hydrogen itself is in nuclear equilibrium.

An effective Proton–proton chain network can track protons, deuterium, helium-3, helium-4 and lithium-7 after eliminating a rapidly adjusting beryllium intermediate. Identical-reactant event rates include a factor one half. Stoichiometric changes must count both helium-4 products of lithium-7 proton capture. The weighted baryon mass density $n_1+2n_2+3n_3+4n_4+7n_7$ is conserved by the nuclear reactions at fixed volume.

In 1972, Stuart J. Freedman and John F. Clauser excited a calcium beam from an oven with a filtered deuterium arc lamp. A calcium emission cascade supplied paired photons; an optical lens system, optical filters, removable glass-plate polarizers and cooled photomultipliers selected their linear polarizations. They compared coincidence circuit rates at relative angles $22.5^{\circ }$ and $67.5^{\circ }$, normalized with polarizers removed. Alternating configurations controlled drift; delayed counts estimated accidental coincidences. The violation agreed with quantum mechanics, conditional on their explicit assumption that passage through a polarizer did not change a photon's detection probability. This was a static-analyzer test with a detection loophole, rather than the ideal four-output CHSH inequality statistic.^[ref]

Keep the same bond force constant $k$ and approximate the nuclear masses of ${}^1\mathrm{H}$, ${}^2\mathrm{H}$, and ${}^{35}\mathrm{Cl}$ by $1u$, $2u$, and $35u$, respectively. Here $u$ is the unified atomic mass unit, and $\mathrm{D}$ denotes deuterium. These are chosen mass number approximations, not measured atomic masses or natural-abundance averages. With the same chlorine isotope in both molecules,This is the isotope scaling of diatomic vibration prediction. Doubling the light partner's mass nearly doubles the reduced mass, lowering the frequency by about $28.3\%$. The factor differs from $1/\sqrt{2}$ because the chlorine partner also moves.

For the same small initial stretch $A$ and zero velocity, both models start with the same internal mechanical energy $k{A}^2/2$, but the DCl oscillation period is longer by $\sqrt{72/37}\approx 1.394972$. In the centre of mass frame, $u_{\mathrm{H}}=-35q/36$ and $u_{\mathrm{Cl}}=q/36$ for HCl; $u_{\mathrm{D}}=-35q/37$ and $u_{\mathrm{Cl}}=2q/37$ for DCl. These ratios explain which partner moves and are separate from the change in frequency.

The common-$k$ assumption belongs to the leading Born-Oppenheimer approximation: fixed nuclear positions and electric charges determine the electronic potential energy, while nuclear masses govern motion on that curve. Quantitative comparison requires actual atomic masses, departures from that approximation, and separation of anharmonicity and molecular rotation from the measured band. Even with common $k$, the quantum harmonic oscillator predicts a different zero-point energy for each isotope.^[ref]

An actual infrared spectroscopy route is to compare gaseous hydrogen chloride containing ordinary hydrogen with its deuterium counterpart. A documented teaching experiment uses a stopcock-equipped gas cell with potassium bromide windows and a Nicolet Avatar instrument for Fourier transform spectroscopy. A nitrogen-purged flask and drying tube supply the mixture; the sample also contains sulfur dioxide. An empty sample-chamber background precedes the filled-cell scan. The absorbance trace and peak-picked spectroscopic wavenumbers are the readout. The analysis assigns and indexes ${}^1{\mathrm{H}}^{35}\mathrm{Cl}$ and ${}^2{\mathrm{H}}^{35}\mathrm{Cl}$ transitions separately and propagates transition uncertainties. The instructions distinguish resolvable chlorine isotope peaks from adequate precision for the intended parameter analysis. ^[ref]

A separate published teaching dataset reports a Nicolet Nexus 670 instrument, a $10\mathrm{cm}$ gas cell at $20\mathrm{Torr}$, and chlorine isotope peaks resolved to baseline. These apparatus details describe that dataset, not the preceding experiment. ^[ref]

The inference connects several models. Fit correctly assigned transitions to separate molecular rotation from molecular vibration, then compare the inferred vibrational scales using isotope scaling of diatomic vibration. Agreement supports the approximate separation between a nearly unchanged chemical bond stiffness and a changed nuclear mass. It does not make a real molecule an exact classical spring: anharmonicity shifts the fundamental transition away from the local harmonic oscillator scale, and coupling to molecular rotation shifts individual lines.

For an independent test, calculate the reduced masses from the specified isotopes and compare the inferred values $k=\mu (2\pi c{\tilde{\nu }}_{\mathrm{harm}}{)}^2$. Agreement within uncertainty tests the common-curvature assumption; choosing separate $k$ values just to fit each frequency would not test it. The documented background addresses instrument and ambient contributions, but does not by itself subtract window absorption, identify contaminating species, or calibrate the spectroscopic wavenumber scale. Those require separate controls. Finite spectral resolution, blended peaks, incorrect isotope assignments, and departure from the harmonic approximation near a potential minimum limit the comparison. No numerical measured line positions are assumed here.

Heavy water replaces both ordinary hydrogen nuclei of water with deuterium.

With $A={\lambda }_{33}$, $B={\lambda }_{34}{n}_4$ and $C={\lambda }_{11}{n}_1^2/2$, helium-3 evolves as ${\dot{n}}_3=C-A{n}_3^2-B{n}_3$ after deuterium elimination. The positive root is its unique attracting equilibrium for a fixed background. In the pp-I limit it reduces to $n_1\sqrt{{\lambda }_{11}/(2{\lambda }_{33})}$, making equilibrium abundance strongly temperature dependent.

A hydrogen bond acceptor supplies the electron-rich atom or group $Y$ interacting with the attached hydrogen of a hydrogen bond donor in $X-\mathrm{H}\cdots Y$. An oxygen lone pair in water is the familiar example. Acceptance describes the association, not necessarily acquisition of a proton. A single molecule, such as water or ammonia, can offer both hydrogen bond donor and hydrogen bond acceptor sites. See Glasfeld: water and nitrogen-containing accepting sites.

A lone pair drawn in a Lewis structure is a useful clue, not a universal test. At an ordinary neutral amide nitrogen, the nominal lone pair participates in chemical resonance with the carbonyl group, so this nitrogen is ordinarily excluded from lists of hydrogen bond acceptors; the carbonyl group's oxygen remains an accepting site. Protonation changes the test again: ${\mathrm{NH}}_4^{+}$ has four N–H covalent bonds and no nitrogen lone pair, whereas ammonia has three and one. These distinctions follow the actual chemical form, and accepting ability also depends on the partner and whether it can approach the site. See McMurry: amide resonance and protonation and Glasfeld: limitations of nitrogen accepting sites.

Conversely, an accepting group need not supply an atom-centred lone pair. The following phenol–benzene experiment probes an O–H donor interacting with the ring's distributed $\pi$ electron density.

Saggu, Levinson and Boxer (2011) used infrared spectroscopy to compare $200\mathrm{mmol}{\mathrm{L}}^{-1}$ phenol in different solvents. The O–H stretching peak was at $3611\mathrm{c}{\mathrm{m}}^{-1}$ in ${\mathrm{CCl}}_4$ and $3555.9\mathrm{c}{\mathrm{m}}^{-1}$ in benzene. To distinguish association from a general solvent effect, they diluted the ring-containing chemical solutions with ${\mathrm{CCl}}_4$: two peaks appeared, with the free peak nearly fixed while the associated peak depended on the partner. Replacing the O–H hydrogen with deuterium reproduced the ordering of shifts for different partners. Thus the perturbed O–H group identifies the hydrogen bond donor; the ring supplies the hydrogen bond acceptor. A peak shift alone does not specify a unique molecular geometry. ^[ref]

The roles make recognition between molecules intelligible: exposed hydrogen bond donors need compatible hydrogen bond acceptors in reachable positions. A count of sites is a starting description, not a prediction of binding strength. In water, forming a contact between two partners can also replace their existing hydrogen bonds to water; the Gibbs free energy of association concerns the whole change. See Glasfeld, pp. 15–16: competition with solvent interactions.

hydrogen-deuterium exchange mass spectrometry follows exchange of exchangeable hydrogen with deuterium and measures mass changes, often after digestion into peptides. Protection can reflect hydrogen bonds, burial or restricted structural opening; altered exchange is not uniquely diagnostic of a binding interface.

Near a stable equilibrium bond length $r_e$, a Taylor expansion of a diatomic molecule's potential energy gives $U(r)\simeq U(r_e)+k(r-r_e{)}^2/2$, with $k=U^{\prime \prime }(r_e)>0$. Treat the internal stretch as a harmonic oscillator of reduced mass, initially omitting molecular rotation. If replacing an isotope leaves this local potential energy curvature unchanged, thenHere $c$ is the speed of light and ${\tilde{\nu }}_{\mathrm{harm}}$ a spectroscopic wavenumber; use consistent length units for $c$ and ${\tilde{\nu }}_{\mathrm{harm}}$. The quantum harmonic oscillator gives adjacent energy eigenvalues separated by $h\nu$, with $h$ the Planck constant. Consequently a vibrational photon absorption can probe the same frequency scale, subject to electric-dipole selection rules.

For an explicitly approximate isotope example, assign the atomic masses of ${}^1\mathrm{H}$, ${}^2\mathrm{H}$, ${}^{35}\mathrm{Cl}$, and ${}^{37}\mathrm{Cl}$ the mass-number values $1u,2u,35u,37u$, respectively, with $u$ the unified atomic mass unit. These are stipulated approximations to individual atomic masses, not measured masses of molecules or natural-abundance averages. They giveThe symbol $\mathrm{D}$ means deuterium, ${}^2\mathrm{H}$. Doubling the light partner's approximate atomic mass therefore produces a large predicted shift, whereas changing the already heavy chlorine partner by two mass units produces a much smaller one. A quantitative comparison with infrared spectroscopy requires actual isotope-specific atomic masses, fitted band parameters, and uncertainties, rather than identifying any arbitrary peak with ${\tilde{\nu }}_{\mathrm{harm}}$.

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.

Measure the cluster's baryonic mass, dominated by its intracluster medium with a smaller stellar contribution, and divide it by an independently inferred total mass. If this cluster baryon fraction $f_b$ is representative of the cosmic baryon fraction, then $f_b={\Omega }_b/{\Omega }_m$. Thus the cluster baryon-fraction estimate of matter density isBig Bang nucleosynthesis constrains the baryon-to-photon ratio and hence ${\Omega }_b{h}^2$ through primordial light-element abundances. A specified $h$ converts this to ${\Omega }_b$; the cluster mass and gas-mass estimates also have distance-dependent $h$ factors, which must be used consistently. This method estimates the total nonrelativistic matter density, including baryonic and dark matter. It does not add radiation or vacuum energy to that inferred ${\Omega }_m$; in a matter-dominated interpretation the question's $\Omega$ is this matter parameter.

Primordial deuterium as a baryon-density indicator works because increasing baryon density makes the destruction of deuterium into heavier nuclei more efficient during Big Bang nucleosynthesis. Smaller inferred primordial D/H therefore implies larger ${\Omega }_b$, with the other nucleosynthesis assumptions fixed. For a fixed measured cluster baryon fraction, a smaller primordial deuterium abundance raises the inferred matter density. A low abundance caused instead by stellar processing is not a lower primordial abundance; depletion, gas ejection and the representativeness of clusters must also be assessed before identifying the cluster fraction with the cosmic one.

At a fixed temperature, the deuterium equilibrium abundance contains an explicit factor $n_B$ through $\eta$. Two nucleons have to encounter each other to make one bound nucleus. Their capture rate per volume is proportional to $n_n{n}_p$, whereas nuclear photodisintegration competes through the very numerous thermal photons. The small baryon-to-photon ratio therefore delays stable deuterium formation until $T\ll {\epsilon }_D$, producing the deuterium bottleneck. Changing the cosmological baryon density changes both its release time and the subsequent efficiency of nuclear burning. The residual primordial deuterium abundance is particularly sensitive to these rates: on the usual nucleosynthesis branch, increasing $n_B$ burns more deuterium into heavier nuclei and lowers the remaining D/H. The positive dependence of the instantaneous equilibrium expression and the negative dependence of the final relic D/H are different statements.

By contrast, once nuclear reactions are efficient, nearly every surviving neutron ends in helium-4. Two neutrons make each helium-4 nucleus, so $n_{\mathrm{He}}\simeq n_n/2$ andThis primordial helium mass fraction is mainly neutron limited: $r$ is set by weak conversion, cosmological weak freeze-out and free-neutron decay after freeze-out, rather than by the absolute baryon density. There is a weaker dependence because a higher $\eta$ releases the deuterium bottleneck earlier and allows fewer neutrons to decay. Thus $Y_P$ is comparatively insensitive, not exactly independent, of $n_B$. At such exceptionally small baryon densities that nuclei never form efficiently, the neutron-limited approximation would fail.

Use $X_i=n_i/n_B$, so $X_D$ counts deuterium nuclei per baryon, not baryons bound into deuterium. For the radiative capture reaction, chemical equilibrium and the zero photon chemical potential give ${\mu }_D={\mu }_n+{\mu }_p$. Dividing the three nonrelativistic Maxwell--Boltzmann number densities removes all chemical potentials:Here $g_D=3$ and $g_n=g_p=2$. The binding convention in this PDF is a negative mass defect, $B_D=m_D-m_n-m_p<0$. The positive nuclear binding energy is ${\epsilon }_D=-B_D$, so the abundance is enhanced by $e^{-B_D/T}$, not by $e^{B_D/T}$.

Since $X_D/(X_n{X}_p)=n_B{n}_D/(n_n{n}_p)$, insert the baryon-to-photon ratio $n_B=\eta n_{\gamma }$ and the photon number density $n_{\gamma }=2\zeta (3)T^3/{\pi }^2$. This givesAt equal nucleon mass $m_N$ in the prefactors, $m_D\simeq 2m_N$, the deuterium equilibrium abundance simplifies toThe dimensionless prefactor is about $8.14$. These expressions describe thermal equilibrium; residual primordial deuterium after nuclear reactions cease must instead be obtained from the reaction kinetics.

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.

At neutrino decoupling, the interaction rate is comparable with the Hubble parameter. Equating the supplied rates givesHere the Planck mass is the unreduced $m_P=G^{-1/2}\simeq 1.22\times 10^{19}\mathrm{GeV}$ in natural units, $G_F\simeq 10^{-5}{\mathrm{GeV}}^{-2}$, and the effective number of relativistic energy degrees of freedom is $g_{\ast }=10.75$. The supplied rate gives an order-of-magnitude neutrino decoupling estimate; a precise neutron-proton conversion temperature uses its own weak-rate coefficients.

The chemical equilibrium condition for $n+{\nu }_e\leftrightarrow p+e^{-}$ is ${\mu }_n+{\mu }_{{\nu }_e}={\mu }_p+{\mu }_e$. With negligible lepton chemical potentials, ${\mu }_n={\mu }_p$. The nonrelativistic Maxwell--Boltzmann number density therefore givesbecause $g_n=g_p=2$ and the small mass difference can be neglected in the prefactor. This is the neutron-proton chemical equilibrium relation.

For the deuterium equilibrium abundance, ${\mu }_D={\mu }_n+{\mu }_p$ because photons have zero chemical potential. The sign convention in the paper is important: $B=m_D-m_n-m_p<0$, while the positive nuclear binding energy is ${\Delta }_D=-B$. Using spin multiplicities $g_D=3$, $g_n=g_p=2$ givesWith $n_B=\eta n_{\gamma }$, $m_n\simeq m_p\simeq m_N$, $m_D\simeq 2m_N$ and photon number density $n_{\gamma }=2\zeta (3)T^3/{\pi }^2$, this becomesHere $X_D=n_D/n_B$ is a number fraction; the fraction of baryons bound in deuterium is $2X_D$. The formula estimates the onset of the deuterium bottleneck opening while depletion of free neutrons and protons remains small. Because $\eta \ll 1$, stable deuterium becomes abundant only far below the positive nuclear binding energy, not simply at $T={\Delta }_D$.

Let $T_f$ denote cosmological weak freeze-out of the neutron-to-proton ratio, $r_f=e^{-Q/T_f}$, and let $t_N$ mark efficient nuclear burning. Free-neutron decay after freeze-out givesAlmost every surviving neutron enters helium-4, whose atomic nucleus contains two neutrons and four baryons. Thus neutron-limited helium synthesis givesA longer neutron half-life leaves more neutrons available, increasing the primordial helium mass fraction. If the longer lifetime reflects weaker weak interactions, the accompanying earlier cosmological weak freeze-out can reinforce that effect.

During radiation domination, $H\propto \sqrt{G{g}_{\ast }}{T}^2$ whereas the weak conversion rate is proportional to $T^5$. Hence $T_f\propto (G{g}_{\ast }{)}^{1/6}$. Increasing the effective number of relativistic energy degrees of freedom raises $T_f$, increasing $r_f$, and shortens the interval for free-neutron decay after freeze-out. Both increase the primordial helium mass fraction while nuclear burning remains efficient. The same reasoning gives the primordial helium response to a varying gravitational constant: larger $G$ over this interval produces more helium, with fixed microscopic rates and the approximate Friedmann equation unchanged. Smaller $G$ has the opposite effect. A particular varying-$G$ model must specify its background expansion and time history; the simple trend must not be extrapolated to changes so large that nuclear burning itself fails.

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.

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.