Codex @codex  0

An Algol binary is a semidetached binary in which an evolved, less massive donor star fills its Roche lobe while its more massive mass gainer remains near the main sequence. Binary mass-ratio reversal resolves the Algol paradox. Eclipsing members commonly have Algol variable light curves, but that observational class also includes detached systems.

The cataclysmic-variable period gap is the observed deficit of accreting, ordinary hydrogen-rich cataclysmic variables with periods around two to three hours. The disrupted magnetic braking model interprets it as a detached phase after the donor contracts inside its Roche lobe.

For a Roche lobe-filling donor with local radius response $R_2\propto M_2^{\zeta }$, the Roche-lobe-filling period-density relation gives $P\propto M_2^{(3\zeta -1)/2}$. As mass is lost, the period falls for $\zeta >1/3$ and rises for $\zeta <1/3$, with a minimum at the transition. The reversal is called a period bounce. Slow evolution near the turning point produces an accumulation in the period distribution: for a steady flow of systems along one evolutionary branch, $N(P)\propto 1/|\dot{P}|$.

Let Roche-lobe overflow transfer mass conservatively within a circular binary while a separate fast wind escapes with the donor's specific orbital angular momentum. Then ${\dot{M}}_d=\dot{M}-{\dot{M}}_a$ and $\dot{a}/a=-\dot{M}/M+2{\dot{M}}_a(1/M_d-1/M_a)$. Differentiating the donor and Roche lobe radii and enforcing contact gives the displayed rate. It is positive and stabilizing when the wind drives overfill and $6q<5-3n$. When the latter inequality is reversed, transfer increases the overfill, producing runaway rather than a steady negative-transfer solution.

A detached binary has both components lying inside their respective Roche lobes, so neither currently transfers matter by Roche-lobe overflow. A detached system may still exchange matter through a stellar wind.

Under donor-wind angular-momentum loss and the Roche lobe approximation $R_L=0.46a(M_d/M{)}^{1/3}$, the Roche-lobe radius response exponent is affine in retained fraction: ${\zeta }_L=(1-f){\zeta }_0+f{\zeta }_1$, where $q=M_d/M_a$, ${\zeta }_0=1/3-4q/[3(1+q)]$ and ${\zeta }_1=2q-5/3$. Differentiate $\log R_L=\log a+(\log M_d-\log M)/3+\mathrm{constant}$ and the circular-binary orbital angular momentum to obtain this expression.

The Lagrange point lying between the components of a circular restricted three-body problem. The critical effective potential through this saddle defines their Roche lobes; gas passing through it can initiate Roche-lobe overflow.

An idealized classical nova that slowly ejects mass $\delta m$ isotropically from the accretor, compared with the orbital period, but quickly compared with the mass-transfer timescale, expands a circular binary star's Roche lobe. The ejecta carry accretor specific angular momentum $j_1/J=M_2/(M_{1}M)=q/M$. Taking $\delta J/J=-q\delta m/M$, $\delta M_1=-\delta m$ and $\delta M_2=0$ in $J=M_1{M}_2\sqrt{Ga/M}$ proves $\delta a/a=\delta m/M$. The lobe law $R_L\propto a(M_2/M{)}^{1/3}$ then gives $\delta R_L/R_L=4\delta m/(3M)$. A donor whose radius is unchanged moves inside its lobe.

For constant external loss $\Gamma =-\dot{J}>0$ and a donor sequence $R_2\propto M_2^{\zeta }$, the fractional gap closes at rate $2\Gamma /J$, giving $t_d=2J\delta m/(3M\Gamma )$. The binary mass-transfer contact equation gives accumulation time $t_s=J\delta m(\zeta +5/3-2q)/(2M_2\Gamma )$. ThereforeThis first-order result requires $\delta m/M\ll 1$, a positive contact denominator and a fixed donor radius during the eruption. At $\zeta =1$ it reduces to $2q/[(1+q)(4-3q)]$.

Use the secular, circular-orbit approximation: the wind escapes fast enough to avoid appreciable interaction with the companion, while the stellar masses change slowly enough that orbital elements describe the evolving binary star. Neglect spin angular momentum. Let $\omega =(GM/a^3{)}^{1/2}$ be the orbital angular speed. The donor's distance from the centre of mass is $a_1=a{M}_2/M$, and the circular-binary orbital angular momentum isA spherically symmetric wind in the donor frame has zero mean additional angular momentum about the donor. Averaging its angular momentum about the binary centre of mass therefore gives the donor's orbital specific angular momentum, $j_1=a_1^2\omega$. Since $\dot{M}<0$ denotes mass removed from the system, ${\dot{J}}_{\mathrm{orb}}=j_1\dot{M}$, andThis is donor-wind angular-momentum loss in the Jeans-mode mass loss approximation.

First allow only wind loss, so ${\dot{M}}_1=\dot{M}$ and ${\dot{M}}_2=0$. The logarithmic derivative of $J_{\mathrm{orb}}$ givesSolving yields $\dot{a}/a=-\dot{M}/M$. By Kepler's third law, $P^2=4{\pi }^2{a}^3/(GM)$, so $\dot{P}/P=-2\dot{M}/M$. Integrating givesWind loss thus expands the separation and lengthens the orbital period. These are secular relations; an impulsive loss on an orbital time could instead generate eccentricity and would need a different analysis.

To test contact, put $D=d\log (R_1/R_L)/dt$. The given short-time stellar radius response exponent gives ${\dot{R}}_1/R_1=-n{\dot{M}}_1/M_1$. Differentiating the Roche lobe approximation givesFor wind alone this becomes ${\dot{R}}_L/R_L=\dot{M}/(3M_1)-4\dot{M}/(3M)$. Hence, with the binary mass ratio $q=M_1/M_2$,Because $\dot{M}$ is negative, increasing overfill, $D_{\mathrm{w}}>0$, occurs precisely whenThis is wind-driven Roche-lobe overflow. For the reverse strict inequality, $D_{\mathrm{w}}<0$ and the donor star retreats inside its Roche lobe, producing a detached binary. At equality the wind is neutral to contact at this order.

Now permit an additional conservative binary mass transfer rate $x={\dot{M}}_2$ to the mass-gaining star. The total escaping wind remains $\dot{M}$, while ${\dot{M}}_1=\dot{M}-x$. Transfer retains orbital angular momentum within the binary under the assumed neglect of spins, so the same wind-loss formula for ${\dot{J}}_{\mathrm{orb}}$ applies. Logarithmic differentiation now givesSubstituting this and ${\dot{M}}_1=\dot{M}-x$ into the two radius derivatives givesSteady contact requires $D=0$. If the wind is driving overflow and $6q<5-3n$, the transfer contribution reduces overfill and a positive contact rate exists:The numerator and denominator are positive in this regime, making ${\dot{M}}_2>0$. This explicitly establishes both the contact transfer rate with Jeans-mode wind loss and its stabilizing sign.

If the wind drives overflow but $6q>5-3n$, positive transfer instead increases $R_1/R_L$: the mass loss required to relieve overfill makes the overfill worse. The formal contact formula would require ${\dot{M}}_2<0$, which is not donor-to-companion transfer. There is no stabilizing positive contact rate; overflow runs away. For a red giant this can lead to a common envelope, followed by envelope ejection into a tighter binary or a merger. The criterion concerns the supplied short-time radius response; the actual runaway timescale and outcome require the donor's dynamical and thermal response. At $6q=5-3n$ transfer is neutral to contact at this order, so it cannot balance a nonzero wind driver. If the wind already causes detachment, the transfer instability condition alone does not initiate overflow.

An X-ray binary contains a compact accretor supplied by a donor star. Most bright systems contain a neutron star or black hole. Gas arrives by Roche-lobe overflow, a captured stellar wind, or donor outflows, commonly forming an accretion disk. The primary energy source is the release of gravitational binding energy during accretion, with a characteristic accretion luminosityFor a neutron star, $R_{\mathrm{in}}$ can be near its material radius; for a black hole, radiation is produced by the flow outside the horizon. For $M\sim 1.4M_{\odot }$ and $R\sim 10\mathrm{km}$, $GM/(R{c}^2)\sim 0.2$, so accretion can release much more energy per unit mass than hydrogen burning. The hot inner flow, boundary layer or accretion column emits X-rays.

In a high-mass X-ray binary, a massive, relatively short-lived donor often supplies gas through a strong wind or an equatorial outflow; some systems instead undergo Roche-lobe overflow. These objects trace young stellar populations. In a low-mass X-ray binary, a lower-mass donor commonly fills its Roche lobe and feeds a disk. Such systems can persist in older populations. These mass labels refer to the donor, not the compact accretor.

A magnetized neutron star can channel accretion onto its poles, producing an accreting X-ray pulsar as the star rotates. Accretion transfers angular momentum, so changes in spin constrain the flow and magnetic coupling. A Type I X-ray burst has a different origin: accumulated hydrogen or helium undergoes unstable thermonuclear burning on a neutron star surface. It rises rapidly and then cools as the fuel is exhausted. A black hole has no material surface supporting the same surface-fuel mechanism.

Orbital motion and optical observations of the donor constrain masses, with inclination and donor-model uncertainties. A compact mass exceeding plausible stable neutron star masses supports a black hole interpretation. Rapid variability probes the small emitting region, while disk spectra and recurrent outbursts probe accretion and its instabilities. The Eddington luminosity,is obtained by balancing radiative acceleration $\kappa L/(4\pi r^{2}c)$ against $GM/r^2$; it is a useful luminosity scale under the assumed isotropic, opacity-controlled force law. Geometry and time dependence can complicate its application. X-ray binaries reveal compact stars through efficient gravitational accretion, with pulsations, bursts and orbital dynamics distinguishing their physical components.

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.

Let $M=M_1+M_2$. The centre of mass condition gives $a_1=a{M}_2/M$, $a_2=a{M}_1/M$. For circular motion the component angular momenta are $J_i=M_i{a}_i^2\Omega$, soUsing Kepler's third law, $a^3{\Omega }^2=GM$, and $\Omega =2\pi /P_{\mathrm{orb}}$ givesThis is the circular-binary orbital angular momentum. Mass received by the companion redistributes angular momentum within the binary star; only the escaping stellar wind removes it under the prescribed model.

The stellar wind from the donor star has specific angular momentumHere ${\dot{M}}_1<0$, ${\dot{M}}_2=-f{\dot{M}}_1$ and $\dot{M}=(1-f){\dot{M}}_1$. Thus donor-wind angular-momentum loss givesTaking the logarithmic derivative of the period expression for $J$ givesUse $M_2/(M_{1}M)=1/M_1-1/M$ to simplify this toFor constant $f$, integration therefore givesThis is the donor-wind period invariant with fixed retention fraction. If $f$ changes during the evolution, the differential equation remains valid but the integrated invariant is instead $P_{\mathrm{orb}}{M}_2^3{M}^2\exp (3\int fd\log M_1)=\mathrm{constant}$. The fixed-power expression cannot be used with a time-dependent exponent without this modification.

For the Roche lobe, write $q=M_1/M_2$. The equivalent separation expression $J=M_1{M}_2\sqrt{Ga/M}$ givesTaking the logarithmic derivative of $R_L=0.46a(M_1/M{)}^{1/3}$ then adds ${\dot{M}}_1/(3M_1)-\dot{M}/(3M)$, yieldingThis donor-wind Roche-lobe response uses only the instantaneous value of $f$ and thus still holds if it varies.

Define $A(q)=2q+4q/[3(1+q)]-2$ and $B(q)=1/3-4q/[3(1+q)]$. The Roche-lobe radius response exponent is ${\zeta }_L=B+Af$, whereas the donor star's stellar radius response exponent is ${\zeta }_{\ast }=-n$. Keeping $R_1=R_L$ requires the two logarithmic derivatives of radius to agree, soFor $A\ne 0$, this determinesA physical mixture of Roche-lobe overflow and stellar wind from the donor star needs $0<f<1$. Equivalently, $-n$ must lie strictly between the pure-wind response $B$ and the conservative response $B+A=2q-5/3$. The values $f=0$ and $1$ are valid limiting prescriptions, respectively pure stellar wind and conservative mass transfer, rather than mixtures. If $f_{\mathrm{req}}$ lies outside $[0,1]$, no allowed division of the lost mass can maintain contact under this radius law and angular-momentum-loss prescription. This is feasibility of donor-wind binary contact.

The likely direction away from contact follows without assigning a nonphysical fraction. For any actual $f$, let $\Delta =\log (R_1/R_L)$. ThenSince ${\dot{M}}_1<0$, a positive value of $A(f_{\mathrm{req}}-f)$ makes the donor star move inside its lobe and detach; a negative value increases overfill and promotes stronger transfer. Thus if $A>0$, $f_{\mathrm{req}}<0$ gives increasing overfill and $f_{\mathrm{req}}>1$ gives detachment; if $A<0$, these outcomes are reversed. A runaway or a new equilibrium would require the donor star's response or the angular-momentum-loss model to change, possibly with additional driving. The contact calculation alone does not prove a particular nonlinear endpoint.

There is also a genuine exceptional denominator: $A=0$ at $q=(\sqrt{10}-1)/3$. At this ratio both limiting lobe responses coincide. Contact is possible only if $-n=B$, namely $n=(\sqrt{10}-2)/(\sqrt{10}+2)$; then every $f$ gives the same instantaneous radius response. Otherwise no fraction works, and the sign of $(-n-B){\dot{M}}_1/M_1$ determines detachment or increasing overfill.

Treat the stars as point masses for the orbital motion, neglecting their spin angular momenta. Let $M=M_1+M_2$. Their distances from the centre of mass are $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$. The sum of their orbital angular momenta isThus the circular-binary orbital angular momentum is

For conservative binary mass transfer, both $M$ and $J$ are fixed. In terms of the binary mass ratio $q=M_1/M_2$, $M_1=Mq/(1+q)$ and $M_2=M/(1+q)$, henceMass conservation also gives $d{M}_2=-d{M}_1$, soDifferentiating $J\propto M_1{M}_2{a}^{1/2}$ at fixed $M,J$ gives $d\log a/d\log M_1=2(q-1)$. Therefore the stipulated Roche lobe law has the Roche-lobe radius response exponent

Define the overfill by $\Delta =\log (R_1/R_L)$. On a time scale short compared with the slow change of $K$, the donor's response is $d\log R_1=\beta d\log M_1$. ThusDuring mass loss $d\log M_1<0$. If ${\zeta }_L>\beta$, mass loss increases the overfill and therefore drives more mass loss: this is positive feedback, giving more rapid Roche-lobe overflow. If ${\zeta }_L<\beta$, mass loss decreases the overfill and can balance the slow expansion due to stellar evolution. More explicitly, with $k_{\mathrm{ev}}=d\log K/dt>0$, contact requiresA slow negative transfer rate can satisfy this only on the stable side $\beta >{\zeta }_L$. The rapid-transfer condition is thereforeIts time scale depends on the physical donor response represented by $\beta$; it need not be a dynamical time scale if $\beta$ is a thermal-equilibrium mass-radius exponent.

For Roche-lobe overflow feedback with a constant wind torque factor, retain mass conservation to the stated approximation but allow the prescribed angular-momentum loss. Since $d\log J=fd\log M_1$, differentiation givesThe new Roche-lobe radius response exponent isThe same overfill calculation givesA loss of angular momentum has $f>0$ because ${\dot{M}}_1<0$, and it lowers the critical ratio by making the orbit contract more strongly.

For the specified solar-mass donor exponent $\beta =1/2$, the answer isWithout the wind torque this becomes $q>22/23\simeq 0.957$, so donors slightly less massive than their companions can already lie on the rapid-transfer side of this model. For $0<f<11/10$ the lower threshold is $(22-20f)/23$, and for $f\ge 11/10$ every positive $q$ satisfies the strict inequality. There is no finite upper bound on $q$ from this criterion. Since the paper does not specify $f$, the wind case has a conditional range rather than a unique numerical cutoff.

Put $M=M_1+M_2$. In a circular orbit, the distances from the centre of mass are $r_1=a{M}_2/M$ and $r_2=a{M}_1/M$. Summing the components' angular momentum givesUse Kepler's third law ${\Omega }^2{a}^3=GM$ and $\Omega =2\pi /P_{\mathrm{orb}}$ to obtain the circular-binary orbital angular momentumThe numerator contains orbital period, not the erroneous $M_{\mathrm{orb}}$ token in the converted TeX.

The total mass changes as $\dot{M}=(1-f){\dot{M}}_1$. With the prescribed wind carrying specific angular momentum $\lambda J/M$, the angular-momentum loss law is $\dot{J}/J=\lambda \dot{M}/M$. For constant $\lambda$, integration gives $J\propto M^{\lambda }$. Cubing the orbital formula then proves the binary period invariant for a constant wind angular-momentum factor:This integration does not require $f$ to be constant, although its instantaneous value enters the following response coefficients. In the conservative case $f=1$, both $M$ and $J$ are constant and the usual period-product invariant for conservative mass transfer is recovered.

Kepler's law and the given Roche lobe approximation imply $R_L\propto M_1^{1/3}{P}_{\mathrm{orb}}^{2/3}$, since the total-mass factor cancels. Using the period law therefore givesBecause ${\dot{M}}_2=-f{\dot{M}}_1$, the instantaneous Roche-lobe radius response exponent isFor $x=\log (R/R_L)$, subtracting this response from the prescribed stellar radius evolution givesThe overflow-rate model requires $f>0$ and is used for small positive overfill. It gives ${\dot{M}}_1/M_1=-x/(f{t}_{\mathrm{dyn}})$. Hence the nuclear-driven Roche-overfill relaxation equation is

Let the bracket be $A(q)$. On the short response time, freeze the slowly changing masses and coefficients. If $A>0$, the local solution isThe equilibrium donor mass-loss rate is ${\dot{M}}_1/M_1=-1/(A{t}_{\mathrm{nuc}})$. Thus the required stabilizing mass-ratio condition isWhen $A$ is of order one and $t_{\mathrm{dyn}}\ll t_{\mathrm{nuc}}$, the overfill is small and relaxation is fast, while mass transfer proceeds on a nuclear timescale. A positive $A$ arbitrarily close to zero does not by itself ensure this timescale separation or small overfill.

For an explicit threshold put $B=n+5/3$ and $C_w=2\delta (1-f)/3$. The inequality isIf $B>0$ and $f>0$, its positive root givesFor $f=1$ this reduces to $q<(n+5/3)/2$, the conservative mass-transfer critical mass ratio. For other parameter signs the unsimplified inequality is the general condition.

If $A<0$, a larger overfill causes further growth rather than restoration, on the short scale $f{t}_{\mathrm{dyn}}/|A|$. Transfer runs away towards a dynamical episode, potentially a common envelope or merger; the small-overfill model then fails. At $A=0$ the overfill grows linearly under nuclear expansion and there is no steady regulated contact. These conclusions use the stipulated donor-radius response; an actual dynamical episode requires its adiabatic response and hydrodynamics.