Codex @codex  0

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.

Apply the preceding temperature profile at the envelope base. Introduce the dimensionless parameterThen $R/r_b-1=q$, so the geometrical thickness of a white-dwarf radiative envelope isThe smallness to be established is therefore $\mathcal{R}{T}_c\ll GM/R$, up to the numerical factor $17/4$.

That ordering follows from the cold, strongly electron-degenerate white dwarf model. Throughout the mass-bearing core, $k{T}_c\ll E_F$, where $E_F$ is the Fermi energy. For nonrelativistic electron degeneracy pressure,The stellar hydrostatic equation gives a characteristic core $P_{\mathrm{deg}}/\rho$ of order $GM/R$. Consequentlyfor fixed composition with order-one molecular-weight ratios, and $s\ll R$. Equivalently, the matched base mass density is much lower than a characteristic core mass density, since ${\rho }_b/{\rho }_{\mathrm{core}}\sim [\mathcal{R}{T}_c/(GM/R){]}^{3/2}$. The degeneracy ordering applies in the bulk core; at the matching point degeneracy is becoming marginal.

The inequalities $T_c\gg T_{\mathrm{eff}}$ and $k{T}_c\ll E_F$ express different facts. The first justifies ignoring surface-temperature terms, while the second makes the envelope geometrically thin. If the core were hot enough that $q$ were not small, the exact expression $s/R=q/(1+q)$ would not justify a thin-layer approximation.

Compare two cooling states of the same white dwarf, holding its mass, nearly temperature-independent degenerate radius and composition fixed. The matching of a white-dwarf envelope to a degenerate core gives $L\propto T_c^{7/2}$, while the geometrical thickness of a white-dwarf radiative envelope gives $s\propto T_c$ at leading thin-layer order. Eliminate $T_c$:Thus a luminosity drop makes the envelope thinner, with $s_2/s_1=(L_2/L_1{)}^{2/7}$. This is the white-dwarf envelope thickness-luminosity relation. Retaining finite thickness in the radiative-zero model gives the correctionwhich reduces to the boxed result for $q_1,q_2\ll 1$.

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

For the evolved-companion fraction, a red giant has another evolved component precisely for unordered $GG$ or $GW$ pairs. By conditional independence at fixed age, their combined slice probability is $g^2+2gw=9g^2$. ThereforeDivide by the system fraction containing at least one red giant to obtain the requested conditional fraction:A double-giant system is counted once in this denominator, and the factor two for a giant–white dwarf pair represents its two possible component orderings.

Use conditional independence at the common age $Y=y$. A system with a red giant and a second evolved component is either a pair of red giants, with conditional probability $g(y{)}^2$, or a red giant and a white dwarf, with conditional probability $2g(y)w(y)$. The factor two counts which component is the white dwarf; a pair of red giants is counted only once. Since $w=9g$, the unconditional probability of this event isThis event is contained in the event of at least one red giant. Therefore the requested conditional probability, among systems selected for containing a red giant, isIn particular the denominator is a fraction of systems, not the individual-star fraction $1/180$.

For this stellar population, let $m=M/M_{\odot }$ denote initial mass in units of the solar mass, and let $\tau =t/\mathrm{Gyr}$ denote present stellar age. Constant formation of equal numbers of stars per unit time makes the stellar age a uniform distribution on $[0,10]$ Gyr, so $Y=\tau /10$ is uniform on $[0,1]$. The normalized initial mass function has probability density function $f(m)=0.1m^{-2}$ for $m>0.1$, since ${\int }_{0.1}^{\infty }{m}^{-2}dm=10$. ConsequentlyThese formulae have the stated age and mass domains; outside them the relevant cumulative fractions saturate at zero or one. The probability integral transform makes $X$ uniform on $[0,1]$: $m=0.1/X$ gives $f(m)|dm/dX|=1$. The time-independent initial mass function and constant number formation rate give independence of $X$ and $Y$. All fractions here count objects, including white dwarfs, using the stipulated stellar evolution law.

A red giant hasThus the mass boundaries are $m=9/\tau$ and $m=10/\tau$ for positive $\tau$. Equivalently, for $0.9<m<1$ the red giant region runs from $9/m$ to the age cap $10$, and for $m\ge 1$ it runs from $9/m$ to $10/m$. There are no red giants with $m\le 0.9$. Boundaries have zero probability and their endpoint convention does not affect the fractions. In the uniform $(X,Y)$ square the red giant region is $9X<Y<10X$, with triangle vertices $(0,0)$, $(1/10,1)$ and $(1/9,1)$. The white dwarf region is the triangle $0<X<Y/10$.

At fixed $Y=y$, define the conditional probabilities of a red giant, white dwarf and main sequence star by $g(y)$, $w(y)$ and $s(y)$. Their interval widths areIntegration over the uniform distribution of age gives the individual-star fractionsThe systems form a coeval binary population: the two binary star components have the same age. Their masses are independent, so $(X_1,X_2,Y)$ has uniform probability density function on $[0,1{]}^3$ and their states have conditional independence given $Y$. Unconditional independence of their states would be incorrect: older systems make both evolved states more likely. The law of total probability now givesThe subtraction removes the double counting of systems containing two red giants.

For Roche-lobe overflow in the stipulated stellar thermal equilibrium, $R_2=R_L$. Using the mass-radius relation and Kepler's third law,The total mass cancels, givingHere $R_{\odot }$ is the solar radius. This is the Roche-lobe-filling period-density relation specialized to $R_2\propto M_2$.

Neglecting spin, the angular momentum of a circular orbit is $J=M_1{M}_2\sqrt{Ga/M}$. Between classical novae the transfer conserves total mass, so $\dot{M}=0$ and ${\dot{M}}_1=-{\dot{M}}_2$. Its logarithmic derivative givesIn particular, for conservation of angular momentum,The source PDF has $\dot{J}=0$ here; the TeX's $J=0$ is a transcription defect. For $q>4/3$, the Roche-lobe radius response exponent exceeds the equilibrium stellar radius response exponent $1$. Since ${\dot{M}}_2<0$, the Roche lobe then shrinks faster than the donor star; the overfill increases and drives more transfer, giving positive feedback. Equality at $q=4/3$ is marginal in this linear response test.

Two qualifications matter. First, $q>4/3$ is a formal extrapolation outside the stated $q<1$ range, and the supplied Roche lobe approximation need not remain accurate there. Second, true dynamical stability of binary mass transfer compares ${\zeta }_L$ with the adiabatic stellar radius response exponent ${\zeta }_{\mathrm{ad}}$, not ${\zeta }_{\mathrm{eq}}=1$. The stellar thermal equilibrium argument supplies the displayed feedback threshold under the imposed radius law, rather than a necessary and sufficient physical dynamical threshold. For example, the polytropic mass-radius relation for a fully convective adiabatic stellar polytrope with index $3/2$ gives ${\zeta }_{\mathrm{ad}}=-1/3$; within this same Roche lobe approximation its dynamical threshold is $q=2/3$. Thus $q<4/3$ by itself does not guarantee dynamical stability.

Gravitational-wave emission from a binary system removes orbital energy and angular momentum. Magnetic braking of a binary star provides another important loss: a magnetized stellar wind carries away donor spin, while tidal locking couples that spin to the orbit. Tides alone redistribute angular momentum and are not an external sink. Matter expelled from the system can also carry orbital angular momentum, although appreciable continuous mass loss would require nonconservative modifications of the preceding equations. Setting ${\dot{R}}_L/R_L={\dot{R}}_2/R_2$ gives the binary mass-transfer contact equationFor $q<4/3$, a negative external $\dot{J}$ therefore sustains negative ${\dot{M}}_2$ along the stipulated stellar thermal equilibrium contact sequence, provided the system is otherwise stable and can remain thermally relaxed.

During the classical nova, the ejecta carry the white dwarf's specific angular momentum, not zero angular momentum. With ${\Omega }^2=GM/a^3$ and the accretor's distance from the centre of mass $a_1=a{M}_2/M$, isotropic escape givesThe slow ejection relative to the orbital period permits the adiabatic circular orbit approximation; the donor star does not transfer appreciable additional mass during it. Take $\delta M_1=-\delta m$, $\delta M_2=0$ and $\delta M=-\delta m$. Differentiating $J=M_1{M}_2\sqrt{Ga/M}$ givesSince $M/M_1=1+q$, this yieldsThese are first-order formulae with errors of order $(\delta m/M{)}^2$. The orbital widening agrees with Jeans-mode mass loss. The donor star's radius is unchanged to this order by the assumed ejection, so its Roche lobe expands away from it, causing nova-induced binary detachment. This idealization neglects changes in the donor from irradiation or interaction with ejecta; those are not supplied in the model.

Let $\Gamma =-\dot{J}>0$ be the fixed external loss rate. During the detached binary phase both component masses are fixed, and ${\dot{R}}_L/R_L=2\dot{J}/J=-2\Gamma /J$. Closing the fractional gap $4\delta m/(3M)$ takesIn the subsequent semidetached binary, the binary mass-transfer contact equation gives $-{\dot{M}}_2=3M_2\Gamma /[J(4-3q)]$. Accumulating a fresh layer of mass $\delta m$ therefore takesThroughout these first-order estimates $J,M,q$ can be evaluated at the start of the cycle: their fractional changes during it are of order $\delta m/M$. The limit $\Gamma =0$ is excluded, since no finite reconnection time follows without an external shrinkage mechanism.

Since $R=a{\xi }_1$, the scaling in part (b) givesSubstitute this into $M=4\pi (R/{\xi }_1{)}^3{\rho }_c{Y}_n$. After collecting powers,whereThis is the polytropic mass-radius relation.

An $n=3$ polytrope models a white dwarf whose electrons are an ultrarelativistic degenerate gas. For $n=3$, the radius exponent vanishes and the mass is independent of both $R$ and ${\rho }_c$. This fixed limiting value is the Chandrasekhar mass: increasing the central density can no longer support a larger mass, so a cold nonrotating white dwarf above the limit cannot remain in hydrostatic equilibrium.

For the nonrelativistic white dwarf, combine the hydrostatic pressure support equation and enclosed mass equation:Eliminating $m$ yieldsIntroduce the Lane-Emden variables for a stellar polytrope,The resulting Lane-Emden equation isThe white dwarf surface is at the first zero ${\xi }_1$ of $\theta$, with $R=a{\xi }_1$.

A cataclysmic variable is a close interacting binary star with a white dwarf accretor and a low-mass donor star, usually a red dwarf, filling its Roche lobe. Matter crosses by Roche-lobe overflow and commonly forms an accretion disk; sufficiently strong white-dwarf magnetism instead permits magnetically channelled accretion. Accretion produces strong variability, while classical novae result from unstable burning of accumulated hydrogen.

The cataclysmic-variable orbital-period distribution contains a short-period population between about 82 minutes and two hours, a pronounced cataclysmic-variable period gap around two to three hours, and a longer-period population above the gap. There is an accumulation near the cataclysmic-variable period minimum. These describe ordinary hydrogen-rich systems; evolved donors and helium-rich systems can have shorter periods. Selection by luminosity and outburst activity affects the observed counts, so $N(P)$ is not simply the intrinsic distribution.

For the contact red dwarf, $R_2=R_L\propto a(M_2/M{)}^{1/3}$. Using its supplied thermal-equilibrium mass-radius relation, $R_2\propto M_2$, gives $a\propto M_2^{2/3}{M}^{1/3}$. Kepler's third law then givesindependent of $M_1$. This is a special case of the Roche-lobe-filling period-density relation.

To examine rapid mass loss, use the stellar radius response exponent ${\zeta }_{\mathrm{ad}}=d\log R_2/d\log M_2=-1/3$. During conservative binary mass transfer, $M$ and $J=M_1{M}_2\sqrt{Ga/M}$ are fixed, and $d{M}_1=-d{M}_2$. ThereforeThe change in overfill is $d\log (R_2/R_L)=({\zeta }_{\mathrm{ad}}-{\zeta }_L)d\log M_2$. Since $d\log M_2<0$, overfill grows, and mass loss runs away, when ${\zeta }_{\mathrm{ad}}<{\zeta }_L$. The dynamical stability of binary mass transfer criterion here is consequentlyThe printed instability inequality is reversed. With the defined $q=M_2/M_1$, the specified radius laws imply the inequality above; $q<q_{\mathrm{crit}}$ is the stable side. Equality is marginal, and the supplied Roche-radius approximation restricts this calculation to its stated mass-ratio range.

Long-term contact is maintained by angular momentum transport out of the orbit. Gravitational-wave emission from a binary system provides one loss mechanism. Magnetic braking of a binary star uses the donor's stellar wind: its magnetic field enforces approximate corotation out to the Alfvén radius, so the wind carries much more specific angular momentum than a nonmagnetic surface outflow. At that radius, roughlyThe geometry contributes order-one factors to the torque. Even when ${\dot{M}}_w\ll |{\dot{M}}_2|$, the large lever arm can remove substantial spin angular momentum. Tidal locking couples donor spin to orbital motion, so the wind torque ultimately brakes the orbit. In practice classical novae can also eject matter; the conservative calculations here isolate the stated idealization.

For the prescribed $\dot{J}/J=-\alpha$, now use the supplied equilibrium law $R_2\propto M_2$ for the slowly evolving donor. Contact implies $\dot{a}/a=(2/3){\dot{M}}_2/M_2$, assuming negligible net wind mass loss. Taking the logarithmic derivative of $J$ givesSince $P\propto M_2$ along this equilibrium sequence,For the dynamically stable mass ratios, this is negative. The result uses thermal equilibrium as well as slow hydrostatic evolution; evolution slow compared only with the dynamical time does not by itself guarantee the radius law $R_2\propto M_2$.

In the disrupted magnetic braking model, the donor above the gap loses mass fast enough to be inflated relative to thermal equilibrium. When the donor becomes a fully convective star, the model assumes a substantial reduction in magnetic braking. The red dwarf then contracts toward equilibrium on its Kelvin-Helmholtz cooling time, becoming smaller than its Roche lobe; accretion largely stops. Gravitational-wave emission from a binary system continues to reduce the separation, carrying the detached binary through the gap until contact is restored near its lower edge.

This explanation requires strong braking immediately above the gap: the donor must already be out of equilibrium and inflated, so there is room to contract and detach when braking drops. With nearly fixed component masses during detachment, $R_L\propto a\propto P^{2/3}$, giving the illustrative inflation factorA weak braking torque that kept the donor in thermal equilibrium just above the gap would not produce this detachment. The model requires a decrease in the torque, not disappearance of the donor's magnetic field.

The exact initial-mass boundaries depend on metallicity, convective overshooting, rotation, and mass loss, but stars near the upper intermediate-mass range, roughly $6.5$--$9M_{\odot }$, can become Super-AGB stars. They ignite carbon but do not immediately ignite neon hydrostatically throughout the core.

Their final fate is set by competition between shell-driven core growth and envelope loss. If a stellar wind removes the envelope first, the remnant is an oxygen-neon, or oxygen-neon-magnesium, white dwarf. If the degenerate core grows toward the Chandrasekhar mass, electron captures on magnesium and neon reduce the electron pressure and can trigger an electron-capture supernova, leaving a neutron star. Slightly higher-mass stars can ignite further fuels and proceed to ordinary iron-core collapse.

Let $a_{\mathrm{ns}}$ be the neutron star's distance from the centre of mass. In a circular orbit,Combining this with Kepler third law,gives the binary mass functionSince $\sin i\le 1$ and $(M_{\mathrm{ns}}+M_{\mathrm{wd}}{)}^2>M_{\mathrm{wd}}^2$,

When a low-mass red giant undergoes stable Roche-lobe overflow, the neutron star receives matter with substantial specific angular momentum, usually through an accretion disk. The accretion torque spins it up and weakens its external magnetic field, producing a recycled pulsar with a millisecond period. Removal of the donor's envelope exposes its helium core as a low-mass white dwarf.

At detachment, the donor mass and core mass both become $M_{\mathrm{wd}}$. Its luminosity is fixed by the red-giant core-mass--luminosity relation, and the supplied radius law therefore makes its final radius $R_L$ a function only of $M_{\mathrm{wd}}$. Cubing the Roche-lobe relation givesEliminating $a$ with Kepler third law yieldsThis is the white-dwarf mass--orbital-period relation: dependence on the neutron-star mass cancels.

For randomly oriented binaries, the isotropic binary inclination distribution isThe observed period gives $M_{\mathrm{wd}}$ from the preceding relation. If $M_{\mathrm{ns}}=1.35M_{\odot }$ is adopted, the measured mass function then givesso the inclination is inferred without astrometry.

Two plausible causes of an apparent shortage of edge-on systems are:

Scatter or bias in the core-mass--period relation can reinforce the second effect.