Past exam of the mathematics course of the University of Cambridge / 2021 / iii / Paper 322 (Binary Stars)

Let $a_{\mathrm{ns}}$ be the neutron star's distance from the centre of mass. In a circular orbit,Combining this with Kepler third law,gives the binary mass functionSince $\sin i\le 1$ and $(M_{\mathrm{ns}}+M_{\mathrm{wd}}{)}^2>M_{\mathrm{wd}}^2$,

When a low-mass red giant undergoes stable Roche-lobe overflow, the neutron star receives matter with substantial specific angular momentum, usually through an accretion disk. The accretion torque spins it up and weakens its external magnetic field, producing a recycled pulsar with a millisecond period. Removal of the donor's envelope exposes its helium core as a low-mass white dwarf.

At detachment, the donor mass and core mass both become $M_{\mathrm{wd}}$. Its luminosity is fixed by the red-giant core-mass--luminosity relation, and the supplied radius law therefore makes its final radius $R_L$ a function only of $M_{\mathrm{wd}}$. Cubing the Roche-lobe relation givesEliminating $a$ with Kepler third law yieldsThis is the white-dwarf mass--orbital-period relation: dependence on the neutron-star mass cancels.

For randomly oriented binaries, the isotropic binary inclination distribution isThe observed period gives $M_{\mathrm{wd}}$ from the preceding relation. If $M_{\mathrm{ns}}=1.35M_{\odot }$ is adopted, the measured mass function then givesso the inclination is inferred without astrometry.

Two plausible causes of an apparent shortage of edge-on systems are:

Scatter or bias in the core-mass--period relation can reinforce the second effect.

The component distances from the centre of mass are $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$, where $M=M_1+M_2$. Their orbital angular momenta add to

Let $W>0$ be the isotropic wind-loss magnitude, so ${\dot{M}}_1=-W$ and ${\dot{M}}_2=0$. The wind carries star 1's specific orbital angular momentum $j_1=a_1^2\Omega$, and henceUsing Kepler third law to write $J=M_1{M}_2\sqrt{Ga/M}$ and differentiating showsSince $P^2\propto a^3/M$, it follows thatThis is Jeans-mode mass loss.

Now allow transfer to star 2 at rate ${\dot{M}}_2>0$ while the wind continues. ThenThe same conservation of angular momentum calculation givesThe donor response $R\propto M_1^{-n}$ iswhile

Before contact, put ${\dot{M}}_2=0$. The wind drives the star farther into its Roche lobe when $\dot{\log }(R/R_L)>0$. With $q=M_1/M_2$, this condition reduces toIf the inequality is reversed, the Roche lobe expands relative to the donor, so wind loss detaches the star and no wind-driven Roche-lobe transfer is sustained; later nuclear expansion may restore contact.

During stable contact, impose $\dot{R}/R={\dot{R}}_L/R_L$ and solve for the transfer rate. Straightforward algebra givesIn the paper's signed notation $W=-\dot{M}$, this is exactlyIfthe stationary response has the wrong sign: transfer enlarges the overfill rather than removing it. Dynamical stability of binary mass transfer is lost, leading to runaway transfer and usually a common envelope or merger.

An intrinsic S-type star is a thermally pulsing AGB star whose third dredge-up exposes carbon and products of the slow neutron-capture process. An extrinsic S-type star has similar surface pollution but no current internal source; it received the material from a former AGB companion. A barium star is the warmer main-sequence or giant counterpart, recognized particularly through strong barium and other slow-neutron-capture elements.

This interpretation predicts a white dwarf companion, the remnant of the former AGB donor. Barium and extrinsic S-type stars are indeed binaries, often with periods long enough that the donor could not have undergone ordinary Roche-lobe overflow. Their abundance patterns, white-dwarf companions, and wide or eccentric orbits therefore point to wind mass transfer in a binary star.

The Bondi–Hoyle accretion estimate treats the companion as moving through a locally uniform wind with relative speed $v_{\mathrm{rel}}$, sound speed $c_s$, and density ${\rho }_w$. Gravity focuses gas from the accretion radiusgivingFor a roughly spherical donor wind, ${\rho }_w\simeq {\dot{M}}_w/(4\pi a^2{v}_w)$ and $v_{\mathrm{rel}}$ combines wind and orbital velocities. This supplies an order-of-magnitude accreted fraction; wind acceleration, density gradients, orbital deflection, and Wind Roche-lobe overflow can substantially change it.

For initial masses $5M_{\odot }$ and $0.7M_{\odot }$, increasing the initial period gives three broad channels:

The transition periods are set by the primary's maximum radius relative to its Roche lobe, and their exact values depend on mass-transfer efficiency, wind speed, and common-envelope energy formalism.