Codex @codex  0

For a circular Kepler orbit, neglecting stellar spin, the orbital angular momentum of a binary star is $J=\mu a^2\Omega =\mu \sqrt{GMa}$, where $M$ is total mass, $\mu$ is the reduced mass, and $a$ is separation. Kepler's third law also gives $J=G^{2/3}{P}^{1/3}{M}_1{M}_2/[(2\pi {)}^{1/3}{M}^{1/3}]$. Distances from the center of mass are $a{M}_2/M$ and $a{M}_1/M$; adding their individual orbital angular momenta proves the formula.

Kepler's laws describe three connected features of an elliptic orbit: its shape, how fast it is traversed, and how its size determines its orbital period. Kepler's first law gives an ellipse with the attracting centre at one focus of a conic; Kepler's second law gives equal swept areas in equal times; Kepler's third law relates the squared orbital period to the cubed semi-major axis. Together they turn a drawing of a Kepler orbit into a prediction of where the bodies will be and when they will return. ^[ref]

Here the mathematical model is an isolated two-body problem in Newtonian gravity, with point masses $M,m>0$ and a bound, noncolliding elliptic orbit. Let $\mathbf{r}={\mathbf{r}}_m-{\mathbf{r}}_M$ be the relative displacement, $r=|\mathbf{r}|$ the separation, and $\mu =G(M+m)$ the gravitational parameter, where $G$ is the gravitational constant. Newton's law of universal gravitation then givesThe relative Kepler orbit has semi-major axis $a$ and orbital eccentricity $0\le e<1$. Both bodies move around their shared astronomical barycenter, the system's centre of mass. In its rest frame their displacements are ${\mathbf{r}}_m=M\mathbf{r}/(M+m)$ and ${\mathbf{r}}_M=-m\mathbf{r}/(M+m)$. Their individual semi-major axes therefore add to $a$; the relative ellipse's attracting focus of a conic represents zero separation, rather than a body fixed in space.

One example will connect all three laws. After nondimensionalization, choose $a=1$, $e=0.6$ and $\mu =1$. The semi-minor axis is $b=a\sqrt{1-e^2}=0.8$, while the periapsis distance and apoapsis distance areThus $a=(0.4+1.6)/2=1$ measures the ellipse's size even though the separation changes fourfold. The specific angular momentum will be $0.8$, the two endpoint speeds $2$ and $0.5$, and the orbital period $2\pi$. These are values in the chosen scaled units, not measurements of a particular planet.

Johannes Kepler published Kepler's first law and Kepler's second law in Astronomia Nova in 1609 and Kepler's third law in Harmonices Mundi in 1619. These publication dates are documented by ^[ref]. Their appeal is that three short statements organize both geometry and timing; Newtonian gravity supplies a common dynamical explanation.

As a heuristic approximation, inserting $r=6GM/c^2$ into the circular orbit form of Kepler's third law, with $M$ here denoting the total mass, givesThe dominant quadrupole gravitational-wave frequency is twice the orbital frequency, $f_{\mathrm{GW}}=2/T$.

Let $f=0.46$. At Roche-lobe overflow contact, the supplied equilibrium radius relation givesKepler's third law then impliesThus the Roche-lobe period-mass relation for a linear donor radius law isHere $R_{\odot }$ and $M_{\odot }$ denote the solar radius and solar mass. Cancellation of the total mass is the special Roche-lobe-filling period-density relation combined with $R_2\propto M_2$.

Between eruptions the conservative binary mass transfer has $M$ constant and ${\dot{M}}_1=-{\dot{M}}_2$. Neglecting spin, writeIts logarithmic derivative givesDifferentiating the Roche lobe relation consequently givesFor $\dot{J}=0$, the Roche-lobe radius response exponent is ${\zeta }_L=2q-5/3$. If the donor follows its specified equilibrium radius sequence, ${\dot{R}}_2/R_2={\dot{M}}_2/M_2$. Its logarithmic overflow therefore changes bySince ${\dot{M}}_2<0$, $q<4/3$ gives negative feedback in this equilibrium-response model, whereas $q>4/3$ gives increasing overflow. For $q>4/3$, loss of donor mass shrinks its Roche lobe faster than its modeled stellar radius, increasing transfer and hence amplifying the original disturbance. The $q>4/3$ comparison is a formal extension; the initial system is specified to have $q<1$.

The word dynamically needs a qualification. Thermal equilibrium specifies ${\zeta }_{\mathrm{eq}}=1$, but actual dynamical stability of binary mass transfer depends on the adiabatic stellar radius response exponent ${\zeta }_{\mathrm{ad}}$, which need not equal one. Stability after a rapid small mass loss requires ${\zeta }_{\mathrm{ad}}>{\zeta }_L$, because $d\log M_2<0$. The equilibrium relation alone does not supply this adiabatic response. As a formal counterexample to deriving a universal dynamical threshold from it, a radius model with a dimensionless entropy label $S$, $R(M_2,S)=C{M}_2^2{e}^S$, can have equilibrium $S_{\mathrm{eq}}=S_0-\log (M_2/M_{\ast })$, giving $R_{\mathrm{eq}}\propto M_2$ but ${\zeta }_{\mathrm{ad}}=2$ at fixed $S$. At $q=3/2$, ${\zeta }_L=4/3<2$, so that model is dynamically stable despite the extrapolated $q>4/3$ condition. This shows the missing response hypothesis, not a claim that this toy radius law describes the given low-mass star.

A more relevant illustration is an ideal fully convective adiabatic stellar polytrope of index $3/2$. At fixed entropy its polytropic mass-radius relation is $R\propto M_2^{-1/3}$, giving ${\zeta }_{\mathrm{ad}}=-1/3$ and a dynamical critical ratio $q=2/3$ under the same Roche lobe approximation. For example, $q=0.9<4/3$ can be unstable dynamically even though the equilibrium comparison has negative feedback. The requested $4/3$ threshold is the intended result with instantaneous response $R_2\propto M_2$; it is not a universal dynamical criterion. The subsequent cycle calculation uses the supplied equilibrium sequence and presumes a physically stable, sufficiently slow transfer regime.

Two important angular-momentum-loss mechanisms are gravitational-wave emission from a binary system and magnetic braking of a binary star. Gravitational waves directly remove orbital energy and angular momentum. A magnetized stellar wind removes donor spin angular momentum with a large magnetic lever arm; tidal synchronization replenishes that spin from the orbit, so the net effect is orbital angular-momentum loss. Treating the wind mass loss as small compared with the Roche-transfer flow is consistent with the conservative approximation used here; otherwise the mass and angular-momentum equations must both be modified. Neglecting the stored spin does not eliminate the torque transmitted through it.

Such losses shrink the Roche lobe and offset the donor's tendency to move inside it. Maintaining contact on the equilibrium sequence requires ${\dot{R}}_L/R_L={\dot{M}}_2/M_2$, hence the binary mass-transfer contact equation becomesFor $q<4/3$ and $\dot{J}<0$, this drives ${\dot{M}}_2<0$ and maintains the modeled semi-detached binary.

Now take $\delta m>0$ for the mass expelled from the accretor in one classical nova. The white dwarf loses $\delta M_1=-\delta m$, while the donor has $\delta M_2=0$ and the total mass change is $\delta M=-\delta m$. Isotropic ejection in the accretor rest frame supplies no mean kick and carries the accretor's specific orbital angular momentum. Its orbital radius is $a_1=a{M}_2/M$, soThis is the angular-momentum budget for isotropic re-emission from a binary star, without an additional spin or frictional torque on the ejecta. The general logarithmic change of $J=M_1{M}_2\sqrt{Ga/M}$ isUsing $M/M_1=1+q$ gives the first-order nova changesThe donor mass and its assumed equilibrium radius are unchanged during the eruption. It is therefore left inside the enlarged Roche lobe, and Roche-lobe overflow ceases in this sharp-contact model. An extended atmosphere, irradiation-induced donor expansion or additional ejecta torques would require a different treatment.

The stipulated circular orbit is an important idealization. Isotropic loss that is slow compared with the orbital period can preserve circularity. If it is truly impulsive compared with that period, an initially circular point-mass orbit instead acquires eccentricity. With unchanged velocity and position, its new specific orbital energy is ${\epsilon }^{\prime }=GM/(2a)-G(M-\delta m)/a=-G(M-2\delta m)/(2a)$ and its squared specific angular momentum stays $h^2=GMa$. The Kepler orbit energy and eccentricity relations then give $a^{\prime }=a(M-\delta m)/(M-2\delta m)$ and $e^{\prime }=\delta m/(M-\delta m)$. The latter still has the same first-order separation increase. The question assumes circularity; its first-order circular budget is consistent with adiabatic orbital mass loss or with subsequent circularization that changes the separation only at second order. Work symbolically with $\delta m$, since the printed numerical layer mass has no explicit unit.

Finally put ${\mathcal{L}}_J=-\dot{J}>0$, constant over both intervals as specified, and let $J,M,q$ denote their leading-order values over one small cycle. While detached, masses and donor radius are fixed, so ${\dot{R}}_L/R_L=2\dot{J}/J=-2{\mathcal{L}}_J/J$. The fractional gap $4\delta m/(3M)$ closes afterOnce contact resumes, the contact equation gives $-{\dot{M}}_2=3M_2{\mathcal{L}}_J/[J(4-3q)]$. Accumulation of the next layer takesDividing yields the detached-to-semidetached time ratioBoth times are first-order in $\delta m/M$. Corrections to either time are second-order in that small mass fraction, so the displayed ratio has relative corrections of order $\delta m/M$. The positive denominator, negligible eruption duration relative to the cycle, fixed donor radius during detachment and common external loss rate are all essential to this idealized nova-induced binary detachment cycle.

Synchronous radius and its limits. Kepler's third law for a circumplanetary orbit givesHere $G$ is the gravitational constant. A planet-synchronous orbit repeats after one spin period. A planet-stationary orbit additionally requires a circular orbit, zero orbital inclination to the equator and prograde motion. Then an antenna fixed on the ground points continuously at the same satellite; merely matching the orbital period does not give this property.

Write the planetary radius as $R_p=(3M/(4\pi \rho ){)}^{1/3}$. Requiring the semi-major axis to exceed $R_p$ gives $P^2>3\pi /(G\rho )$. The stellar tidal force restricts the circumplanetary orbit to the Hill sphere, of radius $r_H=a(M/(3M_{\star }){)}^{1/3}$. ThusThese are the surface and Hill-radius constraints in the idealized spherical, small-$M/M_{\star }$ model. Long-term prograde stability generally requires a radius appreciably inside the Hill sphere; for nonzero orbital eccentricity, the surface constraint applies to the periapsis, not just to the semi-major axis.

Collision time and the launch population. A phase-mixed isotropic swarm occupies a shell of radial scale $e{a}_s$. Its number density scales as $N_s/(e{a}_s^3)$, its relative speed as $a_s/P$, and its geometric collision cross-section as $r^2$. Consequently the total collision rate scales as $N_s^2{r}^2/(Pe{a}_s^2)$. To give the numerical normalization used here, adopt an effective shell volume $V=4\pi e{a}_s^3$. At a fixed position the velocity directions are uniformly distributed in the tangent plane, so their mean relative speed is $4v_K/\pi =8a_s/P$. For an unordered pair, the collision cross-section is $4\pi r^2$, givingThus the mean interval between collisions is $A/N_s^2$ in this large-population kinetic model. The effective shell width fixes an order-one coefficient: small orbital eccentricity and random planes alone do not specify a unique radial probability density. Phase mixing, negligible gravitational focusing, $r\ll e{a}_s$, and uncorrelated encounters are implicit in this estimate. Exactly identical orbital periods with perfectly fixed phases do not themselves produce a memoryless collision process.

For nearly planet-stationary orbits with $e\ll I\ll 1$, the swarm volume is smaller by a factor of order $I$, while the relative speed is smaller by the same factor, since vertical motion of scale $I{v}_K$ dominates the eccentric motion. The two changes cancel in the rate $n\sigma v$: there is no parametric factor $I$ in the collision time within the same phase-mixed kinetic approximation. Numerical factors and phase correlations can differ. This is not an argument that bringing all satellites into one nearly circular plane makes their phases random.

With $N_s=Rt$, an Inhomogeneous Poisson process has cumulative hazard functionThe survival function of the first collision is $\exp [-H(t)]$. Setting the expected number of collisions to one givesThis is a characteristic first-event population, with probability $1-e^{-1}$ of an earlier event. It is not a median: the median has an additional factor $(\log 2{)}^{1/3}$.

Which population collides next? Immediately after the first disruption let $S=N_{s0}-2$ and $F=2x^{-3}$. The latter follows from mass conservation for equal-density spherical pieces. Using the same unordered-pair counting and geometric collision cross-sections as above, define $b=(1+x{)}^2/2$; thenHere $S^2,F^2$ stand for the usual large-population approximations to $S(S-1),F(F-1)$. The factor two distinguishing identical and different species is essential. The probability that the next collision is fragment–fragment isIt is the most likely type if $S<4/[x(1+x{)}^2]$ and $S<2/x^2$. It has probability greater than one half ifFor $x\ll 1$, a strongly fragment-dominated next event therefore requires $Sx\ll 4$; the fragment–satellite comparison is more restrictive than the satellite–satellite comparison. If $S=0$, the next event must be fragment–fragment, provided fragments remain.

Population equations and the normalization discrepancy. Every satellite–satellite event produces $2x^{-3}$ fragments and destroys two satellites. Every satellite–fragment event produces a net $x^{-3}-1$ fragments and destroys one satellite; every fragment–fragment event destroys two fragments. Consistent collision counting therefore givesFor $x\ll 1$, $b\simeq 1/2$ and $x^{-3}-1\simeq x^{-3}$. The last two terms in $A{\dot{N}}_f$ are then twice those in the printed equation. This is a genuine factor-of-two inconsistency: equal-size fragment pairs have a collision cross-section smaller by $x^2$, so their event rate must be $x^2{N}_f^2/A$ if the satellite event rate is $N_s^2/A$; destroying both fragments necessarily gives the sink $-2x^2{N}_f^2/A$.

If the printed equation is taken as a prescribed approximate rate model instead, its implicit event rates are $N_s^2/A$, $N_f{N}_s/(4A)$ and $x^2{N}_f^2/(2A)$. Under precisely that mixed normalization, its corresponding satellite equation isThe printed source term also neglects the consumed fragment in a satellite–fragment event, a legitimate relative $O(x^3)$ approximation. That approximation does not repair the pair-counting discrepancy.

The ensuing cascade. This is a two-size fragmentation cascade. First use the printed approximate model, dropping satellite–satellite events as requested. Put $S=N_s$, $F=N_f$, $k=4x^2$ and take $x<1/2$. ThenIntegrating this linear differential equation, with $F_0=2x^{-3}$ and $S_0=N_{s0}-2>0$, givesInitially the collisional cascade grows if $S_0>8x^2$, with approximate early exponential growth time $4A{x}^3/S_0$ when the satellite population is nearly fixed. The fragment population reaches its maximum atAfterwards satellites are depleted and fragment–fragment losses dominate. For $S_0>0$ the continuum solution has $S\to 0$, $F\to 0$, withWhen no satellites remain initially, $F(t)=F_0/[1+x^2{F}_{0}t/A]$ directly. Small integer populations eventually invalidate these deterministic differential equations.

The consistently counted model has, to leading order in $x$, exactly the same curve $F(S)$ and peak, but both retained time derivatives are twice as large. Its growth time is $2A{x}^3/S_0$, and $F(t)\sim A/(2x^{2}t)$. Keeping $(1+x{)}^2$ and the consumed fragment replaces $k$ by $4x^2/(1+x{)}^2$ and $x^{-3}$ in the curve by $x^{-3}-1$. This explicitly separates the physical collision bookkeeping from the printed normalization while giving the evolution under both conventions.

Write $R_p$ for the physical planet radius and $R_c=\eta R_{H1}$ for the debris-cloud radius; the source uses $R_1$ for both even though their stated size orderings are opposite. A central transit duration of a small planet on a nearly circular Kepler orbit is $D_p\simeq 2R_{\star }/(n_1{a}_1)=P_1{R}_{\star }/(\pi a_1)$. Kepler's third law, with $M_{\star }=4\pi {\rho }_{\star }{R}_{\star }^3/3$, givesA noncentral planetary exoplanet transit is shorter, so within the small-radius and nearly circular approximations, a duration requiringis too long for the planet alone and favors extended debris. For a central cloud transit with $R_c\gg R_{\star }$, the path length is approximately $2R_c$. Using the Hill radius $R_{H1}=a_1(M_1/(3M_{\star }){)}^{1/3}$ gives $D_1/P_1\simeq R_c/(\pi a_1)$, hence the Hill-cloud transit mass estimateThis equality assumes a central cloud chord. An orbit described merely as nearly edge-on need not have central cloud crossings. If its projected cloud impact distance is $s$, the measured half-chord is $\sqrt{R_c^2-s^2}$, so the boxed mass is multiplied by $[1-(s/R_c{)}^2{]}^{-3/2}$. It is a lower estimate when that geometry is unknown.

For the dynamical calculation use ${\mu }_2=M_2/M_{\star }$, $A=G{M}_2/a_2$, and $\mathcal{R}=Af(\alpha )e_1\cos {\varphi }_1$. The derivatives in the supplied Lagrange planetary equations hold the other osculating orbital elements fixed:Since $A/(n_1{a}_1^2)=n_1{\mu }_2\alpha$, the physical mean-longitude equation for a first-order resonant term isThe source's epoch notation needs care when $n_1$ varies. Literally differentiating ${\lambda }_1=n_1(t)t+{ϵ}_1(t)$ adds $t{\dot{n}}_1$; one cannot retain that definition and also identify the printed ${\dot{ϵ}}_1$ with ${\dot{\lambda }}_1-n_1$. A consistent version writes ${\lambda }_1={\int }_0^t{n}_1(s)ds+{ϵ}_{\mathrm{int}}$ and assigns the displayed correction to ${\dot{ϵ}}_{\mathrm{int}}$. Equivalently, the canonical mean-longitude equation for a first-order resonant term follows from the disturbing function directly. The result for $g$ is the intended physical equation under that convention.

Away from resonant-argument libration, neglect changes in $e_1,\alpha ,{\varpi }_1$ on one slow cycle to this order. With $n_0=2\pi /P_1$ the reference mean motion, the signed slow frequency isThis is the near-resonant transit-timing superperiod. Choose ${\varphi }_1(0)=0$. Differentiating $n_1\propto a_1^{-3/2}$ with the semi-major axis equation givesLet $C=3j{\mu }_2{n}_0^2\alpha f{e}_1$ and $n_s=n_1(0)$. Direct integration giveswith zero initial mean longitude. The constant and linear terms are absorbed into the fitted transit ephemeris: the measured mean frequency is $\bar{n}=n_s-C/\omega$. This distinction avoids mistaking an arbitrary initial frequency for the long-term fitted period.

For $|\Delta |\ll 1$, the double-integrated frequency term dominates the direct $g$ term provided $g/f$ stays finite. Relative to the fitted linear ephemeris, ${\lambda }_1(t_k)=2\pi k$ givesHere $P_1$ denotes the fitted mean period at the retained order. The signed-frequency convention fixes the phase sign. Replacing $\omega$ by $2\pi /P_t>0$ changes the sign of the signed $A_t$ when needed; a nonnegative amplitude hasUsing a signed $P_t=2\pi /\omega$ instead gives the source's signed proportionality $A_t/P_t\simeq 3{\mu }_2\alpha f{e}_1/(2\pi \Delta )$. The leading implicit transit equation can equally use $P_{1}k$ inside the sine, since the difference is higher order in the perturbation. This is the eccentricity-enhanced near-resonant transit-timing variation; its derivation assumes small longitude oscillations, ${\mu }_2{e}_1\ll {\Delta }^2$ at fixed resonance coefficients, and the effectively fixed eccentricity used in this truncated model. Dominance over forced-eccentricity timing terms typically also requires $e_1\gg |\Delta |$ at fixed $j$. A full transit calculation also includes forced orbital eccentricity and the conversion from mean longitude to the actual sky-crossing angle. Those are lower order in this eccentricity-enhanced regime, rather than identically absent; broader near-resonant formulas are derived by Lithwick, Xie and Wu.

Inside a mean-motion resonance, ${\varphi }_1$ undergoes resonant-argument libration and its evolution must be solved together with $a_1,e_1,{\varpi }_1$. The constant circulation frequency $\omega$ and the associated ${\Delta }^{-2}$ expansion then fail. At approximately fixed eccentricity, the resonant resonant pendulum approximation has libration frequency of order $j{n}_0\sqrt{{\mu }_2\alpha |f|e_1}$, up to numerical factors. The transit-timing variation follows that libration period and depends on its libration amplitude of a resonant argument. Exactly at a resonant equilibrium there need be no libration signal; near a separatrix the period grows, and eccentricity dynamics can modify the simple estimate. Resonant amplitudes remain controlled by the bounded resonant motion rather than diverging as the circulating detuning tends to zero.

Write $M=M_1+M_2$, $a=a_1+a_2$, and $\mu =M_1{M}_2/M$, the reduced mass. In the center of mass frame, $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$. Both components of the circular orbit have the same angular velocity $\Omega$. Thus their total orbital angular momentum isUsing Kepler's third law, ${\Omega }^2=GM/a^3$, and the orbital period $P=2\pi /\Omega$, gives the circular-binary orbital angular momentum

Take ${\dot{M}}_1<0$, so ${\dot{M}}_2=-f{\dot{M}}_1$ and $\dot{M}=(1-f){\dot{M}}_1$. A parcel in an isotropic stellar wind has, on average, the donor star's orbital velocity. Its wind velocity relative to the donor star averages to zero, so its mean specific orbital angular momentum about the center of mass is $j_1=a_1^2\Omega$. The escaping mass per unit time is $-(1-f){\dot{M}}_1$. With negligible stellar spin, no additional wind torque, and internal redistribution of the angular momentum of retained matter, donor-wind angular-momentum loss therefore givesIsotropy is in the donor star's frame: it does not make the escaping orbital angular momentum vanish. This is also different from isotropic re-emission from a binary star, where matter escapes from the accretor.

Logarithmically differentiate the expression for $J$ along a slowly evolving sequence of circular orbits:Substitution of the donor-wind angular-momentum loss and mass rates givesFor constant $f$, the last expression is $-3fd\log M_1/dt-3d\log M_2/dt-2d\log M/dt$. Hence the period invariant for constant-fraction donor-wind mass loss isFor a time-dependent $f$, the differential equation still holds, but this integrated power law does not. At $f=1$ it reduces to the period-product invariant for conservative mass transfer; at $f=0$, $M_2$ is constant and $P{M}^2$ is constant, as in Jeans-mode mass loss.

Set the binary mass ratio $q=M_1/M_2$. Since $a^3=GM{P}^2/(4{\pi }^2)$, the Kepler third law and the preceding orbital period derivative implyThe specified approximation to the Roche lobe gives $\log R_L=\log (0.46)+\log a+\frac{1}{3}(\log M_1-\log M)$. Consequently the donor-wind Roche-lobe response isThis is a local Roche-lobe radius response exponent, so it remains valid instantaneously even if $f$ varies.

The stellar radius response exponent of $R_1\propto M_1^{-n}$ is ${\zeta }_{\ast }=-n$. Maintaining Roche-lobe overflow in exact contact requires ${\dot{R}}_1/R_1={\dot{R}}_L/R_L$, yieldingFor a precise feasibility of donor-wind binary contact test, put ${\zeta }_0=\frac{1}{3}-\frac{4q}{3(1+q)}$ and ${\zeta }_1=2q-\frac{5}{3}$. The Roche-lobe radius response exponent is ${\zeta }_L=(1-f){\zeta }_0+f{\zeta }_1$. Thus a permitted contact fraction exists exactly whenUnless $A={\zeta }_1-{\zeta }_0=2q+4q/[3(1+q)]-2$ vanishes, the required fraction is $f=(-n-{\zeta }_0)/A$. At $q=(\sqrt{10}-1)/3$, both limiting Roche-lobe radius response exponents coincide: contact is possible for any $f$ only if $n=5/3-2q$, and for no $f$ otherwise.

The sign of the overfilling change resolves the failure of contact:If ${\zeta }_{\ast }>\max ({\zeta }_0,{\zeta }_1)$, mass loss makes the donor star underfill its Roche lobe: the system detaches and contact-driven transfer stops. If ${\zeta }_{\ast }<\min ({\zeta }_0,{\zeta }_1)$, mass loss increases the overfilling: transfer is destabilized, and rapid transfer or a common envelope may result. Calling this a failure of dynamical stability of binary mass transfer specifically requires $-n$ to be the adiabatic stellar radius response exponent; a thermal or equilibrium response concerns a different timescale. Additional angular momentum losses or intrinsic stellar expansion can change these outcomes by changing the contact equation.

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.