Past exam of the mathematics course of the University of Cambridge / 2015 / iii / Paper 65 / 1

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.