Codex @codex  0

Use $\mathcal{R}$ for the stellar gas constant, reserving $R$ for the stellar radius. The specific internal energy of a monatomic gas obeying the ideal gas law is $e=P/[(\gamma -1)\rho ]=3P/(2\rho )$. At fixed enclosed mass, the first law of thermodynamics gives the heat available to the outgoing luminosity as $ϵ=-{\partial }_{t}e-P{\partial }_t(1/\rho )$ when there is no stellar nuclear fusion heating. ThereforeThe second coefficient includes both the mass density dependence of $e$ and pressure work. This is gravitational energy generation in a homologously contracting ideal-gas star.

Write the characteristic scales asThen $\rho ={\rho }_{0}b$, $P=P_{0}p$, $T=T_{0}p/b$, $m=Mq$ and $L_r=Ll$. Since $d/dr=R^{-1}d/dx$, direct substitution in the stellar hydrostatic equation and mass conservation givesFor radiative diffusion in a star, dividing the temperature equation by $T_0/R$ givesIn particular the fourth power here is of $\mathcal{R}$, not of the radius.

Because $q(x)$ is fixed and $M$ is constant, a fixed mass label has fixed $x$. Thus stellar homology gives ${\partial }_t\rho =-3\rho \dot{R}/R$ and ${\partial }_{t}P=-4P\dot{R}/R$. The heating rate reduces to $ϵ=-3P\dot{R}/(2\rho R)$, positive during Kelvin-Helmholtz contraction. Substitution in $d{L}_r/dr=4\pi r^2\rho ϵ$ now yieldsThe dimensionless profiles fix $D$ and $E$. In a common stellar homology family with fixed stellar composition and opacity normalization, the expression for $D$ impliesIt is independent of radius and constant in time for a particular fixed-mass star. The formula for $E$ gives $\dot{R}=-2EL{R}^2/(3G{M}^2)$, which integrates toWhen the initial-radius term is negligible,This is the Kelvin-Helmholtz contraction time law in the prescribed stellar homology model. Infinite initial radius is a limiting idealization; a finite initial radius retains the first term.

At main sequence arrival, CNO cycle heating with fixed stellar composition has luminosityThe integral is dimensionless and fixed within the stellar homology family. The constant-opacity radiative scaling still gives $L\propto M^3$. Equating the stellar nuclear fusion and transported luminosities givesFinally the Kelvin-Helmholtz contraction law evaluated at $R_{\mathrm{MS}}$ gives $t_{\mathrm{MS}}\propto M^2/(L{R}_{\mathrm{MS}})\propto M^{-34/19}$. In a coeval star cluster, the mass just arriving at the main sequence consequently satisfiesThese results comprise constant-opacity homologous contraction to CNO ignition; their mass exponents compare fixed-composition models with common dimensionless profiles.

To derive the nonrelativistic electron-degeneracy pressure in a hydrogen-helium mixture, take fully ionized hydrogen and helium with hydrogen mass fraction $X$, the Electron number density isHere helium nuclei have mass $4m_p$ in the approximation used in the question. Filling the two spin states of each momentum mode up to the Fermi momentum $p_0$ givesThe isotropic momentum flux gives pressure $\frac{1}{3}\int pvdn$. For nonrelativistic Electrons $v=p/m_e$, so the electron degeneracy pressure isSubstitute the composition-dependent $n_e$ and simplify powers of two:Comparing $5/3=1+1/n$ gives the polytropic index $n=3/2$. The result assumes complete Electron degeneracy, $p_0\ll m_{e}c$, and negligible nuclear mass corrections.

Now adopt the stated approximate partially degenerate equation of state, keeping composition and hence $K,\mu$ fixed during contraction. Set $q=m/M$ and define the mass-averaged temperature ${\overline{T}}_M={\int }_0^{1}Tdq$. The stellar virial theorem from part (a) givesFor an index-$3/2$ stellar polytrope, the permitted gravitational-energy formula is $\Omega =-6G{M}^2/(7R)$. Also, using $\bar{\rho }=3M/(4\pi R^3)$,Thus, with $A=3K{\int }_0^1[3\rho /(4\pi \bar{\rho }){]}^{2/3}dq$,The normalized density $\rho /\bar{\rho }$ is a fixed function of $q$ for this homologous polytropic structure. Therefore its integral and $K$ fix $A$, which is independent of $M$ and $R$ for a common composition. The mean here is mass-weighted, rather than a volume average.

For the mass-averaged temperature maximum in a partially degenerate polytrope, put $\alpha =2\mu GM/(7\mathcal{R})$ and $\delta =\mu A{M}^{2/3}/(3\mathcal{R})$, both positive. As a function of $s=1/R$ the temperature is the concave quadraticIt reaches its maximum when $R=2\delta /\alpha =7A{M}^{-1/3}/(3G)$, withAt large radius, contraction increases the temperature as in a thermally supported star. Once Electron degeneracy provides a sufficiently large fraction of the support, further contraction can lower the thermal temperature. The formal zero-temperature radius is $R_{\mathrm{deg}}=\delta /\alpha =R_{\max }/2$, with $R_{\mathrm{deg}}\propto M^{-1/3}$. The continuation to negative temperatures would be unphysical; the adopted sequence only describes its positive-temperature part.

The contraction-limited central temperature has the same mass scaling. Indeed, if the total polytropic pressure is $P=K_{\mathrm{tot}}{\rho }^{5/3}$, the approximate equation of state gives $T=(\mu /\mathcal{R})(K_{\mathrm{tot}}-K){\rho }^{2/3}$. The fixed density profile therefore fixes $T_c/{\overline{T}}_M$. A sufficiently small mass never reaches the central temperature required for sustained hydrogen burning, producing a hydrogen-burning minimum mass. Such an object becomes a cooling, degeneracy-supported brown dwarf instead of settling onto the main sequence. More massive contracting objects can ignite hydrogen before reaching this degeneracy-limited maximum. This conclusion uses the nonrelativistic, fixed-composition model; it does not require extrapolation into the relativistic regime.

The result is in the broad range expected for an old, relatively faint stellar population, but is high compared with many observed stellar mass-to-light ratios of a few solar units in optical bands. The comparison must specify the band and whether the measured mass includes dark matter: for example, early-type galaxy dynamical measurements find an $I$-band normalization near $3.8$ solar units at velocity dispersion $200\mathrm{km}{\mathrm{s}}^{-1}$, including a contribution from dark matter. That is an illustrative observational comparison, not a direct equality with a bolometric stellar calculation.

For an old stellar population, an upper surviving main sequence mass near $M_{\odot }$ is reasonable because more massive stars have shorter lifetimes; low-mass main sequence stars also contain much of the surviving stellar mass. However, red giants and other evolved stars can contribute a large fraction of the light despite their small numbers, so omitting them generally makes this calculated mass-to-light ratio too large. Omitting stellar remnants has the opposite effect on the inferred mass, and a single cubic mass-luminosity relation is only approximate across the stated range. Finally, a single $m^{-3}$ initial mass function down to $0.1M_{\odot }$ is steeper and more low-mass-heavy than commonly adopted distributions, which flatten at low mass. Ten solar units is a useful toy-model estimate, not a precision stellar-population prediction.

Write $x=m/M_{\odot }$ and absorb units into the initial mass function normalization, $dN=A{x}^{-3}dx$. The specified mass-luminosity relation is $L(x)=L_{\odot }{x}^3$, normalized to the Sun. Therefore the total stellar mass and luminosity areThe unknown normalization cancels. The mass-to-light ratio isLow-mass main sequence stars dominate the stellar mass, while equal intervals of $x$ contribute equal luminosity under these particular exponents. More generally the mass-to-light ratio of a power-law stellar population with these exponents on $a<x<b$ is $(M_{\odot }/L_{\odot })/(ab)$.

Low-mass main sequence stars live for much of the Milky Way's age, retaining representatives of many generations, while phase mixing disperses their birth associations around their orbits. In an approximately axisymmetric galactic disk, a sufficiently complete local census can therefore sample the mixed annular stellar population at the Sun's radius despite its small volume. This is a representative local-radius sample, not automatically an unbiased whole-disc census: different vertical scale heights, selection effects and radial migration require corrections when inferring integrated population fractions.

Keep the fixed stellar homology profile, so that ${\rho }_c=C_{\rho }M/R^3$ and $P_c=C_{P}G{M}^2/R^4$. Solving the partially degenerate equation of state for the central temperature giveswhere $A=C_P/C_{\rho }$ and $B=C_{\rho }^{2/3}$ are positive constants. The ion mean molecular weight is distinct from the total-particle mean molecular weight used in part (a).

Set $y=1/R$. The expression in parentheses is a concave quadratic in $y$. Its maximum occurs atSubstitution gives the contraction-limited central temperatureAt this point electron degeneracy pressure supplies half the required central pressure. Further contraction increasingly goes into the density-dependent degeneracy term rather than heating the ions. The formula applies along the positive-temperature part of the idealized homologous sequence; a negative value beyond it is not a physical temperature.

For sufficiently small mass, this maximum never reaches the temperature needed for sustained hydrogen burning. Such an object stays below the hydrogen-burning minimum mass and becomes a cooling brown dwarf, rather than an ordinary main sequence star. Some such objects burn their initial deuterium, but that finite reservoir does not provide lasting hydrogen burning.

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.

In stellar homology, the dimensionless radial profiles are the same after scaling radius, enclosed mass, pressure, temperature and luminosity. At corresponding radii let $r^{\prime }=\xi r$, $m_r^{\prime }=m{m}_r$, $T^{\prime }=tT$, $P^{\prime }=pP$ and $L_r^{\prime }=l{L}_r$. The printed $R^{\prime }=\xi r$ is understood as this local radial scaling, with surface radius $R^{\prime }=\xi R$.

Let $q={\rho }^{\prime }/\rho$ and keep composition, opacity coefficient and nuclear coefficient fixed first. Mass conservation gives $q=m/{\xi }^3$, and hydrostatic equilibrium gives $p=mq/\xi =m^2/{\xi }^4$. The ideal gas law then gives $t=p/q=m/\xi$. Scaling the stellar energy-generation rate equation and the stellar radiative temperature gradient gives respectivelyThe second follows from $dT/dr\propto -\kappa \rho L_r/(r^2{T}^3)$, not from energy production. Equating the two luminosity scalings givesThis assumes the denominator is nonzero and a consistent homologous family exists. If the mean molecular weight differs by $u={\mu }^{\prime }/\mu$, then $t=um/\xi$ andRatios of the opacity and energy-generation coefficients multiply the right-hand side. Thus fixed composition is a real restriction, not an automatic property of every stellar sequence.

Use the Stefan–Boltzmann law $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$ to locate the family on a Hertzsprung-Russell diagram. It implies $T_{\mathrm{eff}}\propto M^{(Y-2X)/4}$ andFor proton–proton chain burning take the usual local approximation $\eta =4$, with Kramers' opacity law $\lambda =1$, $\nu =-7/2$. ThenFor the CNO cycle with electron-scattering opacity, $\lambda =\nu =0$, soA conventional local choice $\eta =16$ gives $R\propto M^{15/19}$ and slope $76/9\simeq 8.44$. Choosing $\eta =17$ instead gives slope $60/7\simeq 8.57$. The source specifies no numerical nuclear exponents, and the effective exponent changes with temperature; the general expression is the unambiguous answer.

The Hertzsprung-Russell diagram places hotter stars to the left. Its branches rise toward higher luminosity and mass, and the CNO cycle/electron-scattering opacity branch has the steeper logarithmic slope. Their illustrative joining point and normalization are arbitrary because the proportional opacity and reaction laws do not specify absolute stellar scales. These are the fully radiative ideal-gas homology predictions, not exact observed main sequence relations: convection and increasing radiation support limit those assumptions in real stars.

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.

A cataclysmic variable is a close semidetached binary in which a white dwarf accretes from a cool, usually low-mass donor star filling its Roche lobe. In an ordinary hydrogen-rich system the donor star is often near the lower main sequence. In a nonmagnetic cataclysmic variable (a disk-fed CV), the white dwarf's magnetic field is too weak to control the flow over the disk. This leaves a characteristic disk-fed geometry, summarized in the original schematic below. NASA's introduction to cataclysmic variables describes the basic components.

Roche-lobe overflow passes through the inner Lagrange point $L_1$. A nearly ballistic stream bends in the rotating frame and strikes the outer accretion disk, producing an accretion-disk stream-impact hot spot. Its retained angular momentum prevents direct radial infall. Viscous evolution of an accretion disk transports angular momentum outward while gas moves inward through a nearly Keplerian accretion disk. Close to the white dwarf, gas slows from orbital rotation toward the stellar rotation in an accretion-disk boundary layer. The disk, accretion-disk stream-impact hot spot, accretion-disk boundary layer, white dwarf and donor star all contribute to the spectrum and, where the orbital inclination permits, to the eclipses of an eclipsing binary. In particular, “nonmagnetic” describes the accretor's control of the flow; it does not imply that the donor star cannot sustain a magnetic field.

The gravitational power available from accretion is approximately $L_{\mathrm{acc}}=G{M}_{\mathrm{WD}}{\dot{M}}_{\mathrm{acc}}/R_{\mathrm{WD}}$. For a slowly rotating white dwarf and a thin steady Keplerian accretion disk, the specific energy changes from approximately zero far out to $-G{M}_{\mathrm{WD}}/(2R_{\mathrm{WD}})$ at the inner disk. Roughly half the available power is radiated by the accretion disk; the rest is released in the accretion-disk boundary layer as the orbital kinetic energy is dissipated. Stellar rotation and departures from a steady thin disk change this partition.

A classical nova has a different energy source. Transferred hydrogen accumulates on the white dwarf; compression heats the base of its envelope until temperature-sensitive hydrogen burning accelerates. Electron degeneracy pressure initially weakens the expansion response to heating, helping a thermonuclear runaway develop. The envelope subsequently expands and ejects material, producing a large optical outburst followed by a decline as the ejecta expand and residual burning ends. The white dwarf usually survives, so continued accretion can build another fuel layer. A rough recurrence scale is the ignition-envelope mass divided by the mean accretion rate; both this mass and the rate vary strongly among systems. The event is an envelope eruption, and the retained fraction is not automatically unity. Starrfield, Iliadis and Hix's nova calculations explains this nuclear mechanism. A classical nova is powered by unstable nuclear burning on the white dwarf.

A dwarf nova undergoes recurrent, shorter brightenings powered principally by enhanced gravitational accretion. The hydrogen-ionization disk instability creates cold, mostly neutral and hot, ionized branches of the disk's accretion-disk thermal S-curve, separated by unstable equilibria. In quiescence the cool disk stores matter because inward transport is slow. Once a critical surface density is reached, a heating transition puts the disk into a hotter, more state with higher effective viscosity: the inward accretion rate and luminosity rise and the disk drains. A cooling transition returns it to quiescence, completing the cycle. A persistent increase in the donor star's transfer rate is not required. Sufficiently high transfer rates can keep the disk on its hot stable branch, giving a nova-like variable rather than ordinary disk cycles. Lasota's disk-instability analysis and Hameury's disk-instability review develop this picture. A dwarf-nova outburst is a disk instability, not a white-dwarf thermonuclear explosion. The names classify mechanisms and need not identify permanently distinct binaries: a nova-producing binary can also possess an unstable disk between nuclear eruptions.

The cataclysmic-variable orbital-period distribution is not smooth. For ordinary hydrogen-rich cataclysmic variables, prominent features are a cataclysmic-variable period gap around two to three hours, a cataclysmic-variable period minimum near eighty minutes, and an accumulation near that minimum. These are population features rather than absolute exclusions. Selection effects matter: luminous high-accretion rate systems are easier to find than faint evolved systems. Helium-transferring binaries have a different period range and are not described by the hydrogen-rich minimum. Gänsicke and collaborators' period-minimum study documents the observed accumulation.

The Roche-lobe-filling period-density relation makes the orbital period a measure of donor structure. Combining $R_d\simeq 0.46a(M_d/M{)}^{1/3}$ with Kepler's third law givesAs the donor star loses mass, its stellar radius response exponent $\zeta =d\log R_d/d\log M_d$ impliesA donor with $\zeta >1/3$ evolves toward shorter orbital periods. When its effective response falls below $1/3$, continued mass loss instead lengthens the orbital period: this is the cataclysmic-variable period bounce. A very low-mass donor may be substellar and increasingly affected by electron degeneracy pressure; the ideal degenerate scaling $R_d\propto M_d^{-1/3}$ illustrates the reversal. The precise cataclysmic-variable period minimum depends on thermal disequilibrium and the strength of orbital angular momentum loss. Near a turning point $|\dot{P}|$ is small; for an approximately steady evolutionary flow of systems, the number per period interval scales as $N(P)\propto 1/|\dot{P}|$, explaining the accumulation. Knigge, Baraffe and Patterson's donor-based evolutionary study relates the donor star's response to these period features.

Long-term transfer is driven mainly by losses of orbital angular momentum, rather than by disk outbursts. Magnetic braking of a binary star removes the cool donor star's spin through a magnetized stellar wind. Tidal synchronization makes the orbit replenish that spin, so the wind extracts orbital angular momentum. This is usually the dominant standard driving mechanism above the cataclysmic-variable period gap. Gravitational-wave emission from a binary system supplies a baseline loss, especially important below the gap. For a weak-field, slowly moving circular orbit, the circular gravitational-wave inspiral givesBecause ordinary transfer from the lighter donor star tends to expand its Roche lobe if orbital angular momentum is conserved, an external loss is needed to sustain contact. In the conservative contact approximation, the binary mass-transfer contact equation readswhere the stellar radius response exponent must match the evolutionary timescale and additional donor expansion has been neglected. The positive denominator on the stable branch makes ${\dot{J}}_{\mathrm{ext}}<0$ drive ${\dot{M}}_d<0$. Nova ejecta can introduce additional nonconservative binary mass transfer. Knigge's evolutionary discussion describes the standard loss mechanisms and their limitations.

In the disrupted magnetic braking model, relatively rapid mass loss above the cataclysmic-variable period gap keeps the donor star inflated relative to stellar thermal equilibrium. When the donor approaches the fully convective star transition, the model postulates a substantial reduction in magnetic braking of a binary star. The donor star can contract within its Roche lobe, suppressing Roche-lobe overflow near the upper edge of the gap. Gravitational-wave emission from a binary system continues to shrink the detached binary; near the lower edge the Roche lobe again reaches the donor radius and transfer resumes. No mass transfer is needed during the detached crossing. The density relation predicts a radius ratio $(3/2{)}^{2/3}\simeq 1.31$ between contact at three and two hours if the masses remain nearly fixed, illustrating the required inflation before detachment. The torque reduction is a model ingredient, not a claim that all fully convective stars lose their magnetic fields. Zorotovic and collaborators' detached-binary study tests the predicted detached population in the gap.

The standard formation channel starts with an initially wider binary star containing two main sequence stars. The initially more massive component evolves first and becomes a giant. Unstable Roche-lobe overflow can engulf the companion in a common envelope. Drag causes inward orbital motion, releasing two-body orbital energy and transferring angular momentum to the envelope. If the envelope is expelled before merger, a close detached binary survives, containing the exposed core, which becomes a white dwarf, and the lower-mass companion. This explains how a binary becomes much tighter than the giant progenitor's radius would have allowed. The outcome depends on envelope binding and the efficiency of energy deposition, summarized approximately by the common-envelope energy formalism; ejection is not guaranteed. Ivanova and collaborators' common-envelope analysis discusses the relevant physics and uncertainties.

Subsequent magnetic braking of a binary star and gravitational-wave emission from a binary system shrink the detached binary until its companion fills its Roche lobe. Stable Roche-lobe overflow, for a suitable binary mass ratio and stellar radius response exponent, then creates a cataclysmic variable. Its later nuclear eruptions, disk cycles and secular orbital period evolution occur on different timescales. The formation sequence is a wide binary, envelope ejection, a close detached white-dwarf binary, and angular-momentum-driven contact.

For a slowly evolving spherical star, use a Lagrangian enclosed mass coordinate $m$. Let $u$ be specific thermal internal energy and $D/Dt$ follow a mass element. The stellar energy balance equation isThe outward luminosity is energy transported through the mass shell, ${ϵ}_{\mathrm{nuc}}$ is nuclear rest-mass energy released per unit mass and time, and ${ϵ}_{\nu }$ is escaping neutrino power per unit mass. Define the gravothermal stellar energy generation by ${ϵ}_{\mathrm{grav}}=-Du/Dt-PD(1/\rho )/Dt$. At fixed composition it is $-TDS/Dt$ by the first law of thermodynamics. During composition changes the full internal-energy derivative must use the appropriate equation of state, with nuclear rest-mass changes counted only once. If a tabulated nuclear rate already subtracts reaction-neutrino energy, that loss must not be subtracted a second time in ${ϵ}_{\nu }$.

In a nearly steady main sequence star, stellar nuclear fusion supplies most of the luminosity. Before stable hydrogen ignition, Kelvin-Helmholtz contraction releases gravitational energy. In a gas-supported hydrostatic star the stellar virial theorem gives $U\simeq -\Omega /2$, so about half the change in gravitational binding energy heats the star and half is available for radiation. The associated timescale is $G{M}^2/(RL)$. After nuclear burning ends, a white dwarf can shine by losing stored thermal energy; crystallization releases latent heat, and composition separation can add gravitational energy. These are heat and gravitational reservoirs rather than sustained hydrogen fusion.

Stellar nuclear fusion is possible because lighter nuclei can combine into products with larger nuclear binding energy. The energy per reaction is $Q=(m_{\mathrm{initial}}-m_{\mathrm{final}})c^2$, including the relevant particles consistently. Although thermal energies are below the Coulomb barrier, quantum tunnelling permits fusion. The thermonuclear reaction rate averages a cross-section over the thermal relative-speed distribution:For nonresonant charged-particle reactions the energy weighting contains $\exp [-E/(kT)-\sqrt{E_G/E}]$; the compromise between the thermal tail and penetration produces the Gamow peak. A resonance can greatly enhance the rate. Density, abundances and temperature therefore all affect the stellar energy-generation rate, and power laws such as $T^4$ or $T^{16}$ are local approximations.

For the proton–proton chain, first produce deuterium and helium-3:The initial reaction is slow because it converts a proton to a neutron through the weak interaction. It is the bottleneck that allows long hydrogen-burning lifetimes. The alternative pep reaction $p+e^{-}+p\to {}^2\mathrm{H}+{\nu }_e$ also feeds the same chain. There are three principal proton-proton chain branches. In pp I, two helium-3 nuclei terminate the chain:In pp II, helium-3 first captures an existing helium-4 nucleus, then electron capture and proton capture finish the branch:In pp III, the beryllium-7 instead captures a proton:Relative branch weights change with temperature and composition. All branches convert four protons into a net helium-4 nucleus, with two weak conversions and two neutrinos, allowing for the electron consumed in pp II. After including positron annihilation, the atomic-mass energy budget is approximately $26.7\mathrm{MeV}$ per net helium-4 nucleus, but the deposited heat is smaller by the branch-dependent escaping-neutrino energy. Near ordinary low-mass main-sequence conditions, ${ϵ}_{pp}\propto \rho X_H^2{T}^4$ is a useful local approximation.

The CNO cycle provides an alternative catalyzed hydrogen burning route in hotter cores. The CNO-I cycle, also called the CN cycle, consists ofThe carbon seed is regenerated: the net reaction again consumes four protons and produces one helium-4 nucleus, two positrons and two neutrinos. The slow ${}^{14}\mathrm{N}(p,\gamma ){}^{15}\mathrm{O}$ reaction governs the ordinary cycle and causes nitrogen-14 to accumulate. The CNO-II cycle branches through ${}^{15}\mathrm{N}(p,\gamma ){}^{16}\mathrm{O}$, ${}^{16}\mathrm{O}(p,\gamma ){}^{17}\mathrm{F}$, ${}^{17}\mathrm{F}\to {}^{17}\mathrm{O}+e^{+}+{\nu }_e$, and ${}^{17}\mathrm{O}(p,\alpha ){}^{14}\mathrm{N}$, after which the main cycle continues. Thus CNO nuclei are catalytic in the closed cycles, although their relative abundances change. A local rate is ${ϵ}_{\mathrm{CNO}}\propto \rho X_H{X}_{\mathrm{CNO}}{T}^{\eta }$ with a much steeper exponent than the PP chain; the preceding homology problem specifies $\eta =16$. The concentration of this heating near the centre favors a convective core in hotter, more massive hydrogen-burning stars.

Once hydrogen is exhausted in a core, contraction can raise its temperature enough for core helium burning. The Triple-alpha process overcomes the absence of stable mass-5 and mass-8 nuclei by maintaining a tiny transient beryllium-8 population:The Hoyle state, an excited carbon-12 resonance near the three-alpha threshold, makes the second capture efficient enough despite beryllium-8's very short lifetime. The net conversion is $3{}^4\mathrm{He}\to {}^{12}\mathrm{C}$, releasing about $7.27\mathrm{MeV}$. Because three alpha particles are required, the specific rate scales as ${\rho }^2{Y}^3$; a schematic resonant rate isIts local temperature exponent is $-3+44/T_8$, of order forty near $T_8=1$. A degenerate core cannot expand promptly to regulate this increase, giving a helium flash; in a nondegenerate core the stellar thermostat permits stable helium burning. The competing ${}^{12}\mathrm{C}(\alpha ,\gamma ){}^{16}\mathrm{O}$ reaction determines much of the eventual carbon-oxygen mixture.

If the core becomes hot enough, carbon burning follows, for example through ${}^{12}\mathrm{C}+{}^{12}\mathrm{C}\to {}^{20}\mathrm{Ne}+\alpha$ and ${}^{12}\mathrm{C}+{}^{12}\mathrm{C}\to {}^{23}\mathrm{Na}+p$. Neon burning begins with photodisintegration, ${}^{20}\mathrm{Ne}+\gamma \to {}^{16}\mathrm{O}+\alpha$, followed by alpha capture such as ${}^{20}\mathrm{Ne}+\alpha \to {}^{24}\mathrm{Mg}+\gamma$. The first step consumes heat, but the combined rearrangement can release net energy. Oxygen burning includes ${}^{16}\mathrm{O}+{}^{16}\mathrm{O}\to {}^{28}\mathrm{Si}+\alpha$ and channels producing phosphorus and other nearby nuclei. At still higher temperatures, silicon burning proceeds through photodisintegration and particle captures in a reaction network approaching quasi-statistical equilibrium, producing iron-group nuclei. It is not simply direct silicon-plus-silicon fusion; the products depend on the electron fraction and weak-interaction timescale.

The increasing nuclear binding energy per nucleon supplies progressively less energy as products approach the iron group. Further fusion past that region cannot provide sustained net heat to support a core. Massive stars can therefore end with an unstable iron-group core; lower-mass stars do not attain all these burning stages and instead leave remnants such as carbon-oxygen white dwarfs. Heavy-element neutron captures can synthesize nuclei beyond the iron group without being the principal hydrostatic power source.

Finally, stellar neutrino energy loss competes with every heat source. Reaction neutrinos accompany the proton–proton and CNO chains. Hot or dense material also emits neutrino pairs through $e^{-}+e^{+}\to \nu +\bar{\nu }$, plasmon decay, $e^{-}+\gamma \to e^{-}+\nu +\bar{\nu }$, and electron-ion bremsstrahlung. These neutrinos usually escape, unlike the photons whose transport is diffusive; neutrino trapping requires much more extreme collapse conditions. Late burning has smaller fuel-energy reservoirs and strong neutrino cooling, so its duration is much shorter than core hydrogen burning. The surface luminosity is the escaping photon power, $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$, not the sum of photons and unobserved neutrino power. Radiative diffusion in a star, convection and electron conduction redistribute energy, while the energy-balance equation distinguishes true sources, losses and storage. Accretion or mass loss adds boundary energy and mechanical work when present, rather than changing the nuclear reaction budget.

Let $T=10\mathrm{Gyr}$ and let $t$ denote present stellar age, rather than time measured from the galaxy's birth. With a fixed initial mass function, constant star formation rate also gives a constant number of births per unit time, as in the birth-age distribution under constant star formation. Counting the supplied stellar evolution model's remnants as part of the population, ages have a uniform distribution on $[0,T]$. Integrating the constant birth rate therefore givesOutside this interval the cumulative fraction is zero or one as appropriate.

Write $m=M/M_{\odot }$ for the birth-mass parameter and $M_{\min }=0.2M_{\odot }$. The normalization cancels when integrating the initial mass function:Thus the mass-tail fraction isIt is one at and below the lower cutoff. The normalization constant has no effect on any number fraction.

For evolutionary-state selection by stellar lifetimes, a red giant has completed its main sequence lifetime but not its giant lifetime. Its stellar lifetime regions in a mass-age diagram are bounded byThe lower boundary meets $t=T$ at $m=\sqrt{0.8}$ and the upper boundary meets it at $m=1$. A white dwarf lies above the upper lifetime curve, within the same age interval.

The probability integral transform, applied to the decreasing mass-tail coordinate, makes $X$ uniform on $(0,1)$: $\mathrm{Pr}(X<x)=x$. The independent birth age gives $Y=t/T$ uniform on $(0,1)$, so number fractions are areas in a unit square. Since $m^{-2}=25X$, the red giant region becomes $20X<Y<25X$, a triangle with vertices $(0,0)$, $(1/25,1)$ and $(1/20,1)$. Its area is $\frac{1}{2}(1/20-1/25)=1/200$. The white dwarf region $Y>25X$ is a triangle of area $\frac{1}{2}(1/25)=1/50$. Hence the present individual-star fractions areThe rest are on the main sequence under the given toy lifetime model.

For the binary stars, draw two independent mass-tail coordinates $X_1,X_2$, but only one shared age coordinate $Y$: components born together are coeval. The coeval binary population therefore occupies a uniform unit cube $(X_1,X_2,Y)$. At a fixed age $Y=y$, the state probabilities for either component areTheir evolutionary states have conditional independence given $y$. They generally do not have unconditional independence, because sharing an age correlates their states. The binary fractions below are slice areas averaged over $0<y<1$.

Use a circular orbit, as in the stated Roche lobe approximation, with orbital angular speed $\Omega$ and negligible stellar spin. The center of mass distances are $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$. Summing the two orbital angular momenta givesKepler's third law, ${\Omega }^2{a}^3=GM$, makes this $J=M_1{M}_2\sqrt{Ga/M}$. In conservative mass transfer, both $M$ and $J$ are fixed, while $d{M}_2=-d{M}_1$. ThusThe Roche-lobe radius response exponent is consequentlyDuring a dynamical mass-loss perturbation, the deep interior and luminosity do not have time to change. The supplied giant structure therefore has stellar radius response exponent ${\zeta }_{\ast }=-0.27$. Its fractional overfill changes by $d\ln (R/R_L)=({\zeta }_{\ast }-{\zeta }_L)d\ln M_1$. The binary mass ratio uses donor mass divided by accretor mass. Since $d\ln M_1<0$, self-limiting Roche-lobe overflow requires ${\zeta }_{\ast }>{\zeta }_L$. Hence the dynamical stability condition isAt equality the linear restoring response vanishes. This conservative mass-transfer critical mass ratio uses the given radius exponent, rather than imposing the different fully convective $-1/3$ approximation.