Codex @codex  0

An isotropic escaping stellar wind has the donor star's mean specific orbital angular momentum. With ${\dot{M}}_d<0$, retained fraction $f$, negligible spin and no additional wind torque, $\dot{J}=(1-f){\dot{M}}_d{a}_d^2\Omega$. The retained matter redistributes angular momentum internally. Unlike isotropic re-emission from a binary star, the escaping material originates at the donor rather than the accretor. With no retention this reduces to Jeans-mode mass loss.

For Roche-lobe overflow, $R_2=R_L$. Cubing the Roche lobe approximation and using the donor's mass-radius relation givesKepler's third law therefore yields the Roche-lobe period-mass relation for a linear donor radius law:The component distances from the centre of mass are $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$. Summing their angular momenta gives the circular-binary orbital angular momentumDuring conservative binary mass transfer, $\dot{M}=0$ and ${\dot{M}}_1=-{\dot{M}}_2$. Logarithmic differentiation at constant angular momentum gives $\dot{a}/a=-2(1-q){\dot{M}}_2/M_2$. Thus the Roche-lobe radius response exponent givesTheir difference is $d\log (R_L/R_2)/dt=(2q-8/3){\dot{M}}_2/M_2$. Since the donor loses mass, this is positive for $q<4/3$: the donor shrinks relative to its Roche lobe, so transfer shuts off without an external driver. The Roche lobe itself expands only for $q<5/6$; for $5/6<q<4/3$ it shrinks more slowly than the donor.

Gravitational-wave emission from a binary system provides such a driver: the orbit emits gravitational waves carrying energy and angular momentum, and the shrinking Roche lobe can maintain contact. Allowing $\dot{J}<0$, the binary mass-transfer contact equation becomesIt gives a negative transfer rate for $q<4/3$.

For the classical nova, assume isotropic re-emission from a binary star: ejecta carry the accretor's specific angular momentum $j_1=a_1^2\Omega$. The eruption lasts many orbital periods, so the orbit responds approximately adiabatically and remains nearly circular; an instantaneous asymmetric kick would be a different model. To first order,Differentiating the circular-binary orbital angular momentum formula during the ejection givesThe accompanying Roche lobe change is thereforeThe donor does not change mass in the eruption and moves inside its expanded Roche lobe: this is nova-induced binary detachment.

Put $\Gamma =-\dot{J}>0$ for the constant external loss. During detachment the component masses are fixed, so the Roche lobe shrinks at fractional rate $2\Gamma /J$. Closing the fractional gap $4\delta m/(3M)$ takes $t_d=2J\delta m/(3M\Gamma )$. During the subsequent semidetached binary phase, the binary mass-transfer contact equation gives $|{\dot{M}}_2|=3M_2\Gamma /[J(4-3q)]$, so accumulating the next eruption's mass takes $t_s=J(4-3q)\delta m/(3M_2\Gamma )$. ConsequentlyBoth times are evaluated to first order in $\delta m/M$, with the same equilibrium donor mass-radius relation; changes of $q,J$ and the transfer rate within one cycle affect higher orders.

To obtain the specified expansion, use the usual nova-ejection prescription: the escaping material is expelled isotropically from the white dwarf and carries its specific orbital angular momentum, with no additional torque. This is isotropic re-emission from a binary star. Circularity alone does not determine the angular momentum carried by the ejecta; the prescription is an additional physical assumption.

The accretor's orbital radius is $a_1=a{M}_2/M$, so its specific angular momentum satisfiesDuring ejection, take $\delta M_1=-\delta m$, $\delta M_2=0$ and $\delta M=-\delta m$, where $\delta m>0$. ConsequentlyDifferentiate $J=M_1{M}_2\sqrt{Ga/M}$ to first order:Using $M/M_1=1+q$ yieldsThe same lobe law now has changing total mass but unchanged donor mass, soThe donor radius is unchanged on the short eruption timescale, whereas its lobe expands. It therefore falls inside the lobe and Roche-lobe overflow stops in this idealized model. This is nova-induced binary detachment.

The need for the ejection prescription can also be seen algebraically. If the ejecta instead have mean specific angular momentum $j_{\mathrm{ej}}={\eta }_{J}J/M$, the same circular-orbit differentiation givesThe requested coefficient is obtained only when ${\eta }_J=q$. Other prescriptions can produce a different expansion or a contraction; the stated circularity condition does not rule them out. The ejection may be short compared with the transfer time while remaining slow compared with an orbit, consistent with circularity. If truly impulsive, an eccentricity is generally excited unless an additional circularization assumption is made.

Now put $\Gamma =-\dot{J}>0$ for the constant external torque. Keep the masses, mass ratio and torque coefficients fixed to leading order over one small-$\delta m$ cycle; corrections in those coefficients affect the calculated times only at the next order. During detachment there is no transfer, so fixed component masses giveThe fractional lobe gap $4\delta m/(3M)$ closes in timeDuring the semi-detached accumulation phase, use the same external torque and the equilibrium donor contact equation derived in (a):Accumulating the next layer therefore takesDivision cancels the layer mass, torque and angular momentum:The result uses $q<4/3$ for a positive steady accumulation rate, unchanged donor structure during the eruption and the same external torque in the two phases. Here the given $q<1$ ensures a positive denominator. It describes an idealized nova cycle; extra frictional loss, irradiation-driven donor changes or a different ejection angular momentum would change the ratio.

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.

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.

In a circular binary star, the distances from the centre of mass are $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$. Summing the two orbital angular momenta givesEquivalently, $J=\mu \sqrt{GMa}$ by Kepler third law.

First consider a rapid conservative perturbation. Both $M=M_1+M_2$ and $J$ are fixed, while $d{M}_1=-d{M}_2$. With $q=M_2/M_1$,soThe Roche lobe formula then gives its mass-radius exponentAfter mass loss, dynamical stability of binary mass transfer requires the donor to shrink relative to its lobe. Since $d\log M_2<0$, this means ${\zeta }_{\mathrm{ad}}>{\zeta }_L$, or

For stable secular conservative binary mass transfer, put $m={\dot{M}}_2/M_2$. The angular-momentum and contact conditions areElimination of $\dot{a}/a$ givesHere magnetic braking of a binary star removes orbital angular momentum, while the donor's expansion maintains Roche-lobe overflow.

In a cataclysmic variable, hydrogen-rich material accumulates on a degenerate white dwarf. Degeneracy prevents initial expansion from regulating its temperature, so nuclear ignition produces the thin-shell instability and a classical nova. Nuclear burning of hydrogen to helium releases about $0.007m{c}^2$, whereas the binding energy at a white-dwarf surface is only of order $G{M}_{1}m/R_1$, typically a few $10^{-4}m{c}^2$. Even modest coupling can therefore eject all the newly accreted envelope without disrupting the white dwarf.

Finally suppose every transferred mass element is expelled by isotropic re-emission from a binary star. Then ${\dot{M}}_1=0$, $\dot{M}={\dot{M}}_2$, and expelled matter carries the white dwarf's specific angular momentum $j_1=(M_2/M{)}^2{a}^2\Omega$. HenceOn the other hand, logarithmic differentiation of $J=M_1{M}_2{M}^{-1/2}{G}^{1/2}{a}^{1/2}$ and of the Roche-lobe radius givesEliminating the separation produces