Integrate the reduced vertical phase-space distribution function over $v_z$. The Gaussian velocity integral cancels its normalization, givingHere the constant vertical velocity dispersion is ${\sigma }_z$, and $\Phi (0)=0$ fixes the potential reference. Neglecting the radial derivatives in the Poisson equation for Newtonian gravity leavesWith $\varphi =\Phi /{\sigma }_z^2$, $\zeta =z/z_0$ and $z_0^2={\sigma }_z^2/(8\pi G{\rho }_0)$, this becomesThus the stated dimensionless equation follows with its factor of two.

For a direct solution, multiply by $2{\varphi }^{\prime }$ and integrate from the midplane. The initial conditions giveFor $\zeta >0$ take the increasing branch and put $u=e^{\varphi /2}$. Then $4(u^{\prime }{)}^2=u^2-1$, whose smoothly increasing solution consistent with the second-order equation is $u=\cosh (\zeta /2)$. Reflection about the midplane supplies the negative-$\zeta$ branch. ThereforeSubstitution verifies $2{\varphi }^{\prime \prime }={\mathrm{sech}}^2(\zeta /2)=e^{-\varphi }$ and both boundary conditions. The smooth second-order equation has a unique solution for those initial data, so no extra constant branch from the first-integral square root is admissible.

This is the self-gravitating isothermal slab. Its scale convention matters: if its profile is instead written $\rho ={\rho }_0{\mathrm{sech}}^2(z/H)$, then $H=2z_0={\sigma }_z/\sqrt{2\pi G{\rho }_0}$. Its surface density is $\Sigma =4{\rho }_0{z}_0$. The force approaches $|{\Phi }^{\prime }|={\sigma }_z^2/z_0=2\pi G\Sigma$, consistent with the field of a plane sheet. These checks fix both the scale-height and Poisson factors.

Magnetic flux concentration and flux expulsion describe redistribution of an existing field, not necessarily creation of net flux. For a weak field in a prescribed conducting flow, the resistive induction equation isIts ideal part combines transport, stretching and compression:Thus a converging flow can intensify a component, and stretching can increase magnetic energy even if the signed net flux is zero. At large magnetic Reynolds number $\mathrm{Rm}=UL/\eta$, diffusion is initially slow on the eddy scale, so magnetic flux freezing carries field lines with the flow. These statements do not imply that an ideal closed eddy can simply lose its initially enclosed flux.

The missing step is the production of small scales by differential circulation. In a two-dimensional incompressible eddy, the Cartesian magnetic flux function obeysDifferent streamlines have different circulation rates. An initially smooth field is wound into fine folds; diffusion and reconnection then smooth the variation of $A$ inside the eddy. As $A$ becomes nearly uniform there, its gradient and hence the interior magnetic field become small. The surviving flux lies mainly near the circulation boundary and adjacent separatrices, where stronger gradients remain. This is flux expulsion. With exactly zero diffusivity, field topology and the flux of material loops remain frozen, so a relaxed flux-free interior is not obtained by pure advection alone. Rigid rotation also fails to generate the differential winding responsible for this process.

The phase mixing in magnetic flux expulsion estimate makes the enhanced role of small diffusion explicit. If the angular-velocity gradient is of order $U/L^2$, a wound radial wavenumber grows as $k_r\sim Ut/L^2$. The accumulated damping exponent is of orderIt becomes order one atThis is much shorter than the large-scale diffusion time $(L/U)\mathrm{Rm}$ at large $\mathrm{Rm}$, but longer than one turnover time. It is a kinematic differential-shear estimate, not a universal time for every turbulent or rapidly evolving convection cell.

A related concentration mechanism occurs at converging flow boundaries. For a locally vertical bundle, conservation of flux gives $B\mathcal{A}\simeq \mathrm{constant}$, so reducing its cross-section increases $B$. In a simple converging boundary region with strain rate $U/L$, balancing strain and transverse diffusion givesThis is the thickness estimate for that local strain-diffusion model; it need not be the thickness or energy scaling of every closed-eddy flux-expulsion problem. Stronger field eventually exerts Lorentz force feedback and changes the flow. The equipartition magnetic field $B_{\mathrm{eq}}=\sqrt{{\mu }_0\rho }U$ is a useful indication that the kinematic assumption can fail, not a universal upper bound on concentrated fields.

At the Sun's photosphere, hot upflows form solar granules, and their spreading horizontal motion transports weak vertical field toward cooler intergranular lanes and downflows. A thin-layer illustration, neglecting vertical transport terms, is ${\partial }_t{B}_z+{\nabla }_h\cdot ({\mathbf{u}}_h{B}_z)=0$. It gives $D_h{B}_z/Dt=-B_z{\nabla }_h\cdot {\mathbf{u}}_h$: horizontal convergence intensifies the vertical flux density. This surface transport is related to the redistribution by convective eddies, but is not identical to the closed-streamline mechanism in the preceding planar example. Larger-scale supergranulation also gathers magnetic elements toward its boundaries, helping organize the solar magnetic network.

Concentration can be further intensified by compressible thermal dynamics. In the thin-tube picture of convective collapse, an unstable downward flow drains gas, reduces internal pressure and allows contraction at approximately conserved flux. Lateral balance $p_i+B^2/(2{\mu }_0)\simeq p_e$ then permits field stronger than the local kinetic equipartition estimate. This evacuation mechanism is excluded by a strict Boussinesq approximation. It is not a description of every evolving surface concentration: gas pressure, dynamic stresses and magnetic curvature can all matter when thin-tube equilibrium fails. The resulting strong, intermittent magnetic elements are associated with intergranular structure; their radiative appearance also depends on thermal structure and opacity. The organized fields of sunspots require a broader account of flux supply and magnetic feedback, rather than attribution to weak-field expulsion alone.

Thus convection redistributes and concentrates flux, finite diffusion enables classical expulsion, and nonlinear magnetic and thermal feedback determines the strong-field outcome. An expulsion argument alone neither supplies the original solar field nor establishes a self-sustaining solar dynamo.

Let $g>0$ be the magnitude of the downward gravity and write ${\alpha }_T$ for the thermal expansion coefficient, to distinguish it from the dynamo coefficient of the previous question. The conductive state has $\mathbf{u}=0$, $\mathbf{B}=B_0\hat{\mathbf{z}}$ and $T_0(z)=T_r-\Delta Tz/d$. Its gravitational term is absorbed in a background hydrostatic pressure.

For the dimensional two-dimensional fields, choose the stream function and Cartesian magnetic flux function conventionsThese automatically give zero divergence. Define the Jacobian determinant bracket $J(f,h)=f_x{h}_z-f_z{h}_x$, so $\mathbf{u}\cdot \nabla h=J(\psi ,h)$. The full magnetic flux function, including the imposed field, must be retained in the magnetic force; using the perturbation alone would lose its linear restoring term.

The scalar vorticity is $(\nabla \times \mathbf{u}{)}_y=-{\nabla }^2\psi$. Also $(\nabla \times \mathbf{B}{)}_y=-{\nabla }^{2}A$, and henceThe buoyancy perturbation is $g{\alpha }_T(T-T_0)\hat{\mathbf{z}}$, whose curl in the $y$ direction is $-g{\alpha }_T(T-T_0{)}_x$. Taking the curl of the momentum equation and multiplying by minus one eliminates pressure and therefore givesWe have restored ${\mu }_0$ in this dimensional expression. The printed magnetic force uses units in which ${\mu }_0=1$, or equivalently a correspondingly rescaled field.

Use length scale $d$, time scale $d^2/\kappa$, velocity scale $\kappa /d$, stream-function scale $\kappa$, magnetic flux-function scale $B_{0}d$, and temperature perturbation scale $\Delta T$. Thus the dimensionless total flux function isHere $x,z$ now denote coordinates divided by $d$, and the same $T$ notation in the second identity denotes the physical temperature. The dimensionless parameters for thermal-diffusion scaling of planar magnetoconvection areThese are the Rayleigh number, Chandrasekhar number, Prandtl number and magnetic-to-thermal diffusivity ratio. The nondimensional magnetic-force coefficient is $B_0^2{d}^2/({\mu }_0{\rho }_r{\kappa }^2)=\sigma \zeta Q$, so omitting the factor $\zeta$ would use a different magnetic parameter.

The resistive induction equation becomes $A_t+J(\psi ,A)=\zeta {\nabla }^{2}A$. Since $J(\psi ,x)=-{\psi }_z$, its perturbation equation has source $+{\psi }_z$. Similarly, advection of the conductive temperature gradient gives $J(\psi ,-z)=-{\psi }_x$, producing the source $+{\psi }_x$ in the perturbation heat equation. The resulting nonlinear equations are thereforeThis derives both forcing signs from the background fields, rather than selecting them by a stream-function convention after the fact.

For the specified roll mode put $K=\pi /\lambda$ and $p_L=K^2+{\pi }^2$, the positive eigenvalue of minus the Laplacian on each spatial factor. Writeand take $\psi =(p_L/K)a{f}_s$, $\Theta =b{f}_T$, $\chi =\lambda c{f}_A$, with amplitudes depending on $\tau =p_{L}t$. These factors satisfy the printed boundary conditions: sine factors make the normal velocity and $B_x=-{\chi }_z$ vanish at the required faces; the second derivatives of the velocity factors give the stress-free conditions; $\Theta$ vanishes at the top and bottom and ${\Theta }_x$ vanishes at the sides.

Discard quadratic perturbation terms. Since ${\nabla }^{2}A={\nabla }^2\chi$, the magnetic force linearizes to $J(x,{\nabla }^2\chi )=({\nabla }^2\chi {)}_z$. The needed derivatives areThe heat and induction equations give $b^{\prime }=a-b$ and $c^{\prime }=a-\zeta c$. In the vorticity equation, its left side is $-(p_L^3/K)a^{\prime }{f}_s$ and its three right-side contributions are $\sigma [(p_L^3/K)a-RKb+\zeta Q{p}_L\pi \lambda c]f_s$. Consequently, withthe three-amplitude vertical-field magnetoconvection system isPrimes now mean $d/d\tau$. They are genuine time derivatives, as the original PDF shows.

For its linear stability, seek amplitudes proportional to $e^{s\tau }$. Eliminating $b,c$ or taking the determinant giveswhereAssume the physical range $\sigma >0$, $\zeta >0$, $q\ge 0$. The Routh-Hurwitz stability criterion gives strict decay exactly when $A_2>0$, $A_3>0$ and $A_1{A}_2>A_3$. At $r=0$ these hold. As $r$ increases, the first loss of stability can occur through a zero eigenvalue or a conjugate imaginary pair.

The stationary threshold isAt this value, $P(s)=s(s^2+A_{1}s+C)$ withIf $C>0$, the other two roots have negative real parts and the zero root is simple. Its derivative with respect to $r$ is $ds/dr=\sigma \zeta /C>0$. This is stationary onset; the reflection symmetry of the roll permits a generic pitchfork bifurcation. The linear calculation does not determine whether the nonlinear branch is supercritical or subcritical.

For oscillatory onset set $s=i\omega$, $\omega \ne 0$. Separating real and imaginary parts gives ${\omega }^2=A_2$ and $A_1{A}_2=A_3$. Solving the latter for $r$ yieldsA genuine Hopf bifurcation onset requires ${\omega }_H^2>0$, soAt this threshold $P(s)=(s+A_1)(s^2+{\omega }_H^2)$, giving the third root $-A_1<0$. The pair crosses transversely: differentiating $P(s,r)=0$ at $i{\omega }_H$ gives

The first instability of the conductive state is therefore stationary at $r_s$ if $\zeta \ge 1$, or if $0<\zeta <1$ and $q<q_{\ast }$. It is oscillatory at $r_H<r_s$ if $0<\zeta <1$ and $q>q_{\ast }$. At $q=q_{\ast }$, the two thresholds meet and $P(s)=s^2(s+A_1)$: the frequency is zero and neither a simple pitchfork nor a nonzero-frequency Hopf test applies. The oscillatory marginality in vertical-field magnetoconvection formula must not be used where its frequency squared is negative.

For conversion back to physical parameters, $R_{s,H}={\lambda }^2{p}_L^3{r}_{s,H}/{\pi }^2$, while the dimensional Hopf frequency is $(\kappa p_L/d^2){\omega }_H$. These are onset conditions for the specified spatial mode. Identifying the globally first mode would additionally require comparison over the allowed spatial modes. Nonlinear saturation and the stability of the bifurcating finite-amplitude states are not determined by this linear three-amplitude system alone.

Choose the poloidal magnetic field convention ${\mathbf{B}}_p=\nabla \times (A\hat{\mathbf{y}})=(-A_z,0,A_x)$. The term $\alpha B$ regenerates this poloidal flux function from the toroidal magnetic field through the alpha effect. It represents a component of the mean-field electromotive force produced by correlated velocity and magnetic fluctuations. The shear $G=\pi V/d$ stretches the poloidal vertical component $B_{p,z}=A_x$ into toroidal field, giving $G{A}_x$: this is the Omega effect. Both $\eta {\nabla }^2$ terms diffuse and damp the field. In a turbulent mean-field model, $\eta$ can denote an effective diffusivity including turbulent transport. The alpha-Omega dynamo approximation retains the shear-dominated toroidal source, and the kinematic dynamo assumption prescribes the velocity and coefficients without their Lorentz force feedback.

For a horizontal wavenumber $k$ and vertical wavenumber $\pi /d$, set $K^2=k^2+{\pi }^2/d^2$ and seek amplitudes proportional to $e^{st}$. The two modal equations areA nonzero mode requires the dispersion relationFor real nonzero $\alpha Gk$, the root with positive real square-root part isThe other root has strictly negative real part. For $k=0$ there is no oscillatory feedback and both roots have real part $-\eta {\pi }^2/d^2$; the possible triangular-system prefactor cannot yield exponential growth.

Let $h=|k|d/\pi >0$ and $D=\alpha V{d}^2/({\pi }^2{\eta }^2)$. The growing-branch exponent satisfiesConsequently the plane-layer alpha-Omega dynamo threshold at fixed $h$ isTo find its least value, differentiate $D_c(h)=2(h^{-1}+2h+h^3)$:The unique positive solution is $h^2=1/3$. The threshold diverges at both ends $h\downarrow 0$ and $h\to \infty$, so this is its global minimum. HenceWith the usual positive-$D$ convention, the onset value is $D_c=32/(3\sqrt{3})$ and exponentially growing waves exist for $D>D_c$. Exactly at the threshold the real part is zero, so it is a neutral dynamo wave: the quoted minimum is the onset threshold, or infimum for strict growth, not a growing solution at equality. If signed $D$ is allowed, the growth condition concerns $|D|$; there is no smallest signed negative value among all growing cases.

The real field is obtained by taking the real part of a Fourier mode. Its phase is $kx+(\Im s_{+})t+\pi z/d$, so the horizontal phase speed is $-(\Im s_{+})/k$. Reversing the sign of $\alpha V$ reverses propagation, while leaving the growth threshold unchanged. Reversing $k$ gives the conjugate real-field representation.

A density deficit in a flux tube is not by itself an instability calculation. For example, a straight magnetized tube in otherwise field-free gas at the same temperature has lateral balance $p_i+B_i^2/(2{\mu }_0)=p_e$. The perfect gas law then gives ${\rho }_i<{\rho }_e$, so the tube rises unless some other force balances its weight deficit. This configuration simply lacks magnetostatic equilibrium. A magnetic buoyancy instability, in contrast, starts from a force-balanced atmosphere and asks whether a small allowed displacement grows. A balanced horizontal magnetic layer can be stable or unstable; the presence of magnetic pressure alone does not decide this.

Assume a smooth equilibrium, with downward gravitational acceleration $g>0$, and write $\mathcal{R}$ for the gas constant per unit mass. The ambient total pressure and force balance areWe use ideal magnetic flux freezing on the displacement time scale and rapid lateral pressure equilibration. A bodily displacement perpendicular to a straight horizontal tube does not change its length or bend its field. If its area is $\mathcal{A}$, conservation of magnetic flux and mass per unit length gives $B\mathcal{A}=\mathrm{constant}$ and $\rho \mathcal{A}=\mathrm{constant}$. Thus the tube preservesEfficient heat exchange fixes its temperature to the ambient temperature at its new height, not to its original temperature. This is the thermal assumption that distinguishes the requested criterion from the adiabatic one.

Take a tube initially at $z_0$ and let $f_0=B_0(z_0)/\rho (z_0)$. After a small upward displacement $\xi$, let ${\rho }_i(z)$ be the density it would have in lateral balance at $z=z_0+\xi$. Its preserved $f_0$ and imposed local temperature giveAt the initial position ${\rho }_i=\rho$. Differentiate this relation at $z_0$:For the ambient atmosphere the same differentiation must also differentiate $f(z)=B_0(z)/\rho (z)$. It givesSubtracting isolates the actual displaced density contrast:The pressure-gradient force on the tube balances the weight of the displaced ambient gas, so its leading vertical buoyancy acceleration isThe coefficient multiplying $\xi$ is positive precisely when the tube becomes lighter on upward displacement. Thus the isothermal interchange criterion for magnetic buoyancy isWhere $B_0>0$, this is simply $d(B_0/\rho )/dz<0$, or $d\log B_0/dz<d\log \rho /dz$. The squared form is independent of the field's choice of direction. Equality is neutral at this order, and the reversed inequality is restoring for these displacements. A layer with an interval satisfying the strict inequality admits local growing interchange disturbances under the assumed pressure/thermal equilibration. The calculation does not address field-bending undular modes, finite heat-transfer times or magnetic diffusion; those can change the criterion. In particular, the localized interchange criterion for a magnetized atmosphere with an adiabatic gas response is not interchangeable with this fast-thermal-exchange result.

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.

In stellar homology, write $r=Rx$, $m=Mf(x)$ and $\rho ={\rho }_{c}h(x)$, with common dimensionless profiles. Integrating hydrostatic equilibrium from the surface, where the pressure is negligible, givesThe profile integral is common to the sequence. Similarly, mass conservation gives $M=4\pi {\rho }_c{R}^3{\int }_0^1{x}^{2}h(x)dx$. These prove the homology scaling of central stellar pressure:The ideal gas law consequently gives $T_c\propto \mu M/R$, with the gravitational and gas constants absorbed into the coefficient.

The radiative diffusion in a star equation has a characteristic gradient $T_c/R$. Hence $L\propto R{T}_c^4/({\kappa }_c{\rho }_c)$. Using Kramers' opacity law, and comparing fixed-composition homologous models, this becomesThe volume integral of the stellar energy-generation rate $ϵ={ϵ}_0\rho T^{\nu }$ instead givesThe dimensionless integrals behind these relations do not depend on $M$ along a homologous sequence. In stellar thermal equilibrium the two luminosities agree, proving Kramers homology with a variable nuclear exponent:For the specified proton–proton chain exponent $\nu =7/2$, the mass factor cancels. ThusThe model captures two important features of the lower main sequence: hydrogen burning primarily uses the proton–proton chain, and gas rather than radiation pressure supplies the principal support. A radiative interior and a Kramers' opacity law approximation are useful near the Sun. It is a simplified local model: the Sun has a convective envelope, the coolest main sequence stars become fully convective, and real opacities and burning exponents vary. In particular a literally constant radius is not a realistic law over the entire lower main sequence.

On this proton–proton chain branch $T_c\propto M$, because $R$ and mean molecular weight are fixed. Normalizing to the Sun gives the upper mass estimateFor the specified CNO cycle exponent $\nu =23/2$, the homology relation instead reads $R^{14}\propto M^8$. Substituting into the radiative luminosity relation givesThis branch retains the deliberately idealized fully radiative assumptions, even though real CNO cycle cores are often convective.

For a Hertzsprung-Russell diagram, use the Stefan–Boltzmann law, $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$. On the proton–proton chain branch $T_{\mathrm{eff}}\propto M^{11/8}$, so $L\propto T_{\mathrm{eff}}^4$. On the CNO cycle branch $T_{\mathrm{eff}}\propto M^{57/56}$, giving $L\propto T_{\mathrm{eff}}^{292/57}$. Thus the upper branch rises slightly more steeply toward the hot, luminous upper left in a logarithmic Hertzsprung-Russell diagram. The sketch below matches the approximate branches continuously at their respective burning crossovers; the true crossover has both heating terms and is smooth.

The proton–proton chain begins with reactions between hydrogen nuclei, so its specific heating scales as $X_H^2\rho T^{7/2}$ in this approximation. The catalytic CNO cycle involves hydrogen and CNO nuclei, and scales as $X_H{Z}_{\mathrm{CNO}}\rho T^{23/2}$. Here $X_H$ is the hydrogen mass fraction and $Z_{\mathrm{CNO}}$ the CNO mass fraction, not the total abundance of every metal. ConsequentlyKeeping $X_H$ fixed while dividing $Z_{\mathrm{CNO}}$ by 256 gives the CNO crossover temperatureThe mass estimate uses the unchanged proton–proton chain branch below crossover, together with the stipulated unchanged opacity and mean molecular weight.

The low-CNO proton–proton chain sequence therefore lies on the same locus, but extends to twice the transition mass. Beyond its later crossover the low-CNO CNO cycle sequence has the same logarithmic slope $292/57$ but a different normalization. At a fixed mass where both sequences are CNO-dominated, lowering the nuclear coefficient by 256 givesThus those low-CNO models are smaller, hotter and modestly more luminous at the same mass. At the same effective temperature, their CNO branch is lower in luminosity than the solar-composition CNO branch; these comparisons hold different quantities fixed.

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

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

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

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

Work at fixed enclosed mass $m$, so the perturbations are fluid Lagrangian perturbations. The equilibrium acceleration vanishes. Varying the factor $r^{-4}$ in the gravitational force contributes $4Gm\delta r/(4\pi r^5)$, while varying the factor multiplying the acceleration gives no first-order term because the background acceleration is zero. With $\delta r=r\xi e^{-i\omega t}$, the linearized momentum equation isConservation of the mass of each displaced shell gives its fractional density change as minus the fractional volume change. ThusThe adiabatic relation is $\delta p=\gamma p\delta \rho /\rho =-\gamma p\chi$. Use $dm/dr=4\pi r^2\rho$ and the stellar hydrostatic equation $p^{\prime }=-Gm\rho /r^2$ to turn the mass-coordinate momentum equation into $(\delta p{)}^{\prime }-{\omega }^2\rho r\xi +4p^{\prime }\xi =0$. It follows that the radial stellar pulsation equation isThe stellar adiabatic exponent may vary with radius, so it must remain inside the derivative.

Expand the derivative and multiply by $r^3$. The self-adjoint differential equation becomesHence $f(r)=\gamma p{r}^4$ and $g(r)=r^3[(3\gamma -4)p{]}^{\prime }+{\omega }^2\rho r^4$ in the requested form. Multiplying by a real eigenfunction and integrating gives the boundary term $[\gamma p{r}^4\xi {\xi }^{\prime }{]}_0^R$. It vanishes for a regular centre and bounded admissible surface displacement and derivative, since $p(R)=0$. Therefore the weighted stellar pulsation Rayleigh quotient isReal coefficients and the symmetric endpoint conditions give the Sturm-Liouville problem; the Rayleigh-Ritz variational principle identifies its minimum with the lowest squared frequency. For complex amplitudes the corresponding quadratic forms use modulus squares, and the real eigenfunctions may be chosen for this calculation.