Past exam of the mathematics course of the University of Cambridge / 2015 / iii / Paper 65 (Binary Stars)

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 the center of mass frame, the velocities are ${\mathbf{v}}_1=(M_2/M)\mathbf{v}$ and ${\mathbf{v}}_2=-(M_1/M)\mathbf{v}$. Hence the kinetic energy is $T=\frac{1}{2}\mu v^2$, while the Newtonian gravitational potential energy is $U=-G{M}_1{M}_2/r=-GM\mu /r$. The two-body orbital energy is thereforeThe relative equation of Newtonian gravity is $\ddot{\mathbf{r}}=-GM\mathbf{r}/r^3$. It is a central force, so conservation of angular momentum confines the motion to a plane and preserves the specific angular momentum $h=|\mathbf{r}\times \mathbf{v}|=r^2\dot{\theta }$.

For completeness, derive the polar equation of a Kepler orbit. Let $w=1/r$ and use primes for $\theta$ derivatives. Then $\dot{r}=-h{w}^{\prime }$ and $\ddot{r}-r{\dot{\theta }}^2=-h^2{w}^2(w^{\prime \prime }+w)$. The radial equation gives the Binet equationChoose the angular origin at closest approach. Its solution is $w=(GM/h^2)(1+e\cos \theta )$, soFor the ellipse, the orbital eccentricity satisfies $0\le e<1$. Its extreme separations are $r_p=l/(1+e)$ and $r_a=l/(1-e)$. The semi-major axis is half their sum, giving $a=l/(1-e^2)$ and thus $l=a(1-e^2)$. The true anomaly $\theta$ runs through a full $2\pi$; the endpoints are the same position, and for a circular Kepler orbit the angular origin is arbitrary.

The radial and transverse relative speeds are $\dot{r}=(GM/h)e\sin \theta$ and $r\dot{\theta }=h/r$. Substituting these into the specific orbital energy givesTherefore the conserved orbital energy isThe negative specific orbital energy and fixed angular momentum characterize the bound Kepler orbit; neither should be assumed unchanged through a mass-ejecting explosion.

Initially the circular Kepler orbit has $r=a$ and $v^2=GM/a$. Treat the supernova kick in a binary star as impulsive: the relative position does not change, the companion's velocity is unchanged during the impulse, and the new neutron star receives $\mathbf{u}$ with $|\mathbf{u}|=\alpha v$. Thus ${\mathbf{v}}^{\prime }=\mathbf{v}+\mathbf{u}$, and with $\psi$ the angle between these vectors before the kick,Write $M^{\prime }=M_1^{\prime }+M_2$, $\beta =M^{\prime }/M$ and ${\mu }^{\prime }=M_1^{\prime }{M}_2/M^{\prime }$. The post-explosion specific orbital energy isUsing ${\epsilon }^{\prime }=-G{M}^{\prime }/(2a^{\prime })$ givesHere $a^{\prime }>0$ for a bound Kepler orbit. For a hyperbolic Kepler orbit, this formula uses the signed energy parameter $a^{\prime }<0$, rather than the positive geometric magnitude of the hyperbola's semi-major axis. At the parabolic Kepler orbit boundary, $1/a^{\prime }=0$.

The kick-direction binary survival criterion is bound if $S<2\beta$ and unbound with positive asymptotic speed if $S>2\beta$. To prove the printed sufficient disruption condition, choose a perpendicular kick, $\cos \psi =0$: then $S=1+{\alpha }^2$, so $M^{\prime }<\frac{1}{2}(1+{\alpha }^2)M$ produces positive specific orbital energy. This is sufficient, not the sharp existence threshold. Sincethe complete kick-direction binary survival criterion isThe second inequality follows by taking a kick directly opposite to $\mathbf{v}$. All directions remain bound if $2\beta >(1+\alpha {)}^2$, while every direction has positive escape energy if $2\beta <(1-\alpha {)}^2$. For $\alpha >0$, the marginal direction obeys $\cos \psi =(2\beta -1-{\alpha }^2)/(2\alpha )$ when the right-hand side lies in $[-1,1]$. Equalities give marginal parabolic Kepler orbits for the relevant extreme direction. With zero kick, losing more than half the original total mass unbinds the circular Kepler orbit.

For an escaping pair, the Newtonian gravitational potential energy approaches zero at infinite separation. Conservation of the post-explosion two-body orbital energy then gives $\frac{1}{2}{\mu }^{\prime }{V}^2={\mu }^{\prime }{\epsilon }^{\prime }$, or the asymptotic relative speed of a disrupted binaryThis is the relative recession speed, not the velocity of the new center of mass or either star's individual velocity in the original inertial frame. A marginal parabolic Kepler orbit separates with $V=0$ at infinity.

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.