Codex @codex  0

Interpret projection from very far away as the asymptotic incoming state at $r\to \infty$. For $v>0$ and impact parameter $b>0$, conservation of energy and conservation of angular momentum give $E=v^2/2$ and $|h|=bv$. At the closest approach the radial velocity vanishes. With $k=-\kappa <0$, the central-force radial turning point equation isOnly the positive root is physical. Put $q=\kappa /v^2$. The pericentre distance and the purely tangential speed there areThus $p<b$ and $u>v$: attraction bends the path inward and increases the speed. With positive energy, the Kepler orbit is a hyperbola. Measured from the pericentre direction, its polar coordinates satisfy $r=h^2/[\kappa (1+e\cos \varphi )]$ with $e=\sqrt{1+b^2{v}^4/{\kappa }^2}>1$.

The hyperbola shown has a nonzero impact parameter. For $b=0$, the attractive radial orbit instead reaches the singular origin; there is no regular turning point with a finite closest-approach speed.

For constant-density incompressible inviscid flow with conservative body force, Euler equations for an inviscid fluid give ${\partial }_t\mathbf{u}+(\mathbf{u}\cdot \nabla )\mathbf{u}=-\nabla (p/\rho +\Phi )$. Using $(\mathbf{u}\cdot \nabla )\mathbf{u}=\nabla (|\mathbf{u}{|}^2/2)-\mathbf{u}\times \mathbf{\omega }$, take the curl to obtain ${\partial }_t\mathbf{\omega }=\nabla \times (\mathbf{u}\times \mathbf{\omega })$. Since both divergences vanish, the vector identity expands this toWithout incompressibility or the constant-density/barotropic hypothesis, extra compression or baroclinic terms are present; the printed equation uses these flow assumptions.

Let $\kappa (t)=\dot{h}/h$, with $h>0$. Axisymmetric mass conservation is $r^{-1}{\partial }_r(r{u}_r)+{\partial }_z{u}_z=0$. If $u_r$ is independent of $z$, then $u_z$ is affine in $z$. The impermeable moving-plate conditions $u_z(r,\pm h,t)=\pm \dot{h}$ force $u_z=\kappa z$, and integrating the radial equation gives $u_r=-\kappa r/2+C(t)/r$. Regularity at the axis excludes $C(t)/r$. The meridional part of the axisymmetric inviscid flow between moving parallel plates is thereforeFor no swirl, $u_{\theta }=0$ completes the field. Continuity alone does not determine azimuthal motion; this must be specified separately. More generally, prescribe a smooth regular initial profile $U(r)$ independent of $z$. The azimuthal part of Euler equations for an inviscid fluid isso conservation of $r{u}_{\theta }$ along radial particle paths gives $u_{\theta }(r,t)=\lambda U(\lambda r)$, where $\lambda =\sqrt{h/h_0}$. Together with the meridional components above, this gives a family of complete fields, with a suitable pressure determined by the radial and vertical equations. The field is therefore not unique without rotational data.

With the stated solid-body rotation, $u_{\theta }=\Omega (t)r$ and $\mathbf{\omega }=2\Omega {\mathbf{e}}_z$. The vorticity equation gives $2\dot{\Omega }=2\Omega {\partial }_z{u}_z=2\Omega \kappa$, henceFor a fluid particle, $\dot{r}=-\kappa r/2$, so $r^{2}h$ is constant and $r^2\Omega$ is constant. This is conservation of angular momentum about the axis: radial spreading when the plates approach decreases the angular velocity. The meridional and rotating fields satisfy the inviscid equations with an appropriate pressure; no tangential no-slip condition is imposed on inviscid plates.

Use the separation vector $\mathbf{r}$ from the star to the comet, and let $m_{\mathrm{red}}=m{M}_{\star }/(m+M_{\star })$. The two-body problem separates into uniform centre of mass motion andTaking a cross product with $\mathbf{r}$ proves conservation of angular momentum; taking the scalar product with $\dot{\mathbf{r}}$ proves conservation of energy. The relative constants areThe total relative angular momentum and energy are $m_{\mathrm{red}}\mathbf{h}$ and $m_{\mathrm{red}}\epsilon$, respectively. The conserved eccentricity vector is $\mathbf{e}=\dot{\mathbf{r}}\times \mathbf{h}/\mu -\hat{\mathbf{r}}$. Its scalar product with $\hat{\mathbf{r}}$ gives $r=h^2/[\mu (1+e\cos f)]$, with $f$ the true anomaly. The periapsis and apoapsis distances are therefore $h^2/[\mu (1+e)]$ and $h^2/[\mu (1-e)]$. Their sum is twice the semi-major axis, soAt either apsis $v=h/r$; substituting an apsidal radius gives $\epsilon =-\mu /(2a)$. Thus the specific orbital energy implies the vis-viva equationThese are exact relative-coordinate formulas for an elliptic orbit, $0\le e<1$. In the barycentric frame the comet's speed is $M_{\star }/(M_{\star }+m)$ times the relative speed; the printed identification with the comet's stellar orbital speed uses $M_{\star }\gg m$.

For Roche-lobe overflow in the stipulated stellar thermal equilibrium, $R_2=R_L$. Using the mass-radius relation and Kepler's third law,The total mass cancels, givingHere $R_{\odot }$ is the solar radius. This is the Roche-lobe-filling period-density relation specialized to $R_2\propto M_2$.

Neglecting spin, the angular momentum of a circular orbit is $J=M_1{M}_2\sqrt{Ga/M}$. Between classical novae the transfer conserves total mass, so $\dot{M}=0$ and ${\dot{M}}_1=-{\dot{M}}_2$. Its logarithmic derivative givesIn particular, for conservation of angular momentum,The source PDF has $\dot{J}=0$ here; the TeX's $J=0$ is a transcription defect. For $q>4/3$, the Roche-lobe radius response exponent exceeds the equilibrium stellar radius response exponent $1$. Since ${\dot{M}}_2<0$, the Roche lobe then shrinks faster than the donor star; the overfill increases and drives more transfer, giving positive feedback. Equality at $q=4/3$ is marginal in this linear response test.

Two qualifications matter. First, $q>4/3$ is a formal extrapolation outside the stated $q<1$ range, and the supplied Roche lobe approximation need not remain accurate there. Second, true dynamical stability of binary mass transfer compares ${\zeta }_L$ with the adiabatic stellar radius response exponent ${\zeta }_{\mathrm{ad}}$, not ${\zeta }_{\mathrm{eq}}=1$. The stellar thermal equilibrium argument supplies the displayed feedback threshold under the imposed radius law, rather than a necessary and sufficient physical dynamical threshold. For example, the polytropic mass-radius relation for a fully convective adiabatic stellar polytrope with index $3/2$ gives ${\zeta }_{\mathrm{ad}}=-1/3$; within this same Roche lobe approximation its dynamical threshold is $q=2/3$. Thus $q<4/3$ by itself does not guarantee dynamical stability.

Gravitational-wave emission from a binary system removes orbital energy and angular momentum. Magnetic braking of a binary star provides another important loss: a magnetized stellar wind carries away donor spin, while tidal locking couples that spin to the orbit. Tides alone redistribute angular momentum and are not an external sink. Matter expelled from the system can also carry orbital angular momentum, although appreciable continuous mass loss would require nonconservative modifications of the preceding equations. Setting ${\dot{R}}_L/R_L={\dot{R}}_2/R_2$ gives the binary mass-transfer contact equationFor $q<4/3$, a negative external $\dot{J}$ therefore sustains negative ${\dot{M}}_2$ along the stipulated stellar thermal equilibrium contact sequence, provided the system is otherwise stable and can remain thermally relaxed.

During the classical nova, the ejecta carry the white dwarf's specific angular momentum, not zero angular momentum. With ${\Omega }^2=GM/a^3$ and the accretor's distance from the centre of mass $a_1=a{M}_2/M$, isotropic escape givesThe slow ejection relative to the orbital period permits the adiabatic circular orbit approximation; the donor star does not transfer appreciable additional mass during it. Take $\delta M_1=-\delta m$, $\delta M_2=0$ and $\delta M=-\delta m$. Differentiating $J=M_1{M}_2\sqrt{Ga/M}$ givesSince $M/M_1=1+q$, this yieldsThese are first-order formulae with errors of order $(\delta m/M{)}^2$. The orbital widening agrees with Jeans-mode mass loss. The donor star's radius is unchanged to this order by the assumed ejection, so its Roche lobe expands away from it, causing nova-induced binary detachment. This idealization neglects changes in the donor from irradiation or interaction with ejecta; those are not supplied in the model.

Let $\Gamma =-\dot{J}>0$ be the fixed external loss rate. During the detached binary phase both component masses are fixed, and ${\dot{R}}_L/R_L=2\dot{J}/J=-2\Gamma /J$. Closing the fractional gap $4\delta m/(3M)$ takesIn the subsequent semidetached binary, the binary mass-transfer contact equation gives $-{\dot{M}}_2=3M_2\Gamma /[J(4-3q)]$. Accumulating a fresh layer of mass $\delta m$ therefore takesThroughout these first-order estimates $J,M,q$ can be evaluated at the start of the cycle: their fractional changes during it are of order $\delta m/M$. The limit $\Gamma =0$ is excluded, since no finite reconnection time follows without an external shrinkage mechanism.

Set $\mathbf{r}={\mathbf{x}}_1-{\mathbf{x}}_2$, $M=M_1+M_2$ and $\mu =M_1{M}_2/M$, the reduced mass. Newton's law of universal gravitation gives ${\ddot{\mathbf{x}}}_1=-G{M}_2\mathbf{r}/r^3$ and ${\ddot{\mathbf{x}}}_2=G{M}_1\mathbf{r}/r^3$. Subtraction establishesFor fixed masses, the orbital energy and angular momentum areTaking derivatives gives $\dot{E}=\mu \mathbf{v}\cdot (\ddot{\mathbf{r}}+GM\mathbf{r}/r^3)=0$ and $\dot{\mathbf{J}}=\mu \mathbf{r}\times \ddot{\mathbf{r}}=0$. These are conservation of energy and conservation of angular momentum in the unperturbed Kepler orbit. To see directly that the orbital eccentricity is fixed, put $\mathbf{h}=\mathbf{r}\times \mathbf{v}$ and $\mathbf{n}=\mathbf{r}/r$. The eccentricity vector isSince $\dot{\mathbf{h}}=0$, the cross product identity $\mathbf{r}\times (\mathbf{r}\times \mathbf{v})=\mathbf{r}(\mathbf{r}\cdot \mathbf{v})-r^2\mathbf{v}$ gives $\ddot{\mathbf{r}}\times \mathbf{h}=GM\dot{\mathbf{n}}$, so $\dot{\mathbf{e}}=0$. Its magnitude $e$ is the orbital eccentricity. Equivalently,Thus the semi-major axis and orbital eccentricity of a bound noncollision Kepler orbit are constant, including $e=0$.

Take the centre of mass as origin. The two positions are ${\mathbf{x}}_1=(M_2/M)\mathbf{r}$ and ${\mathbf{x}}_2=-(M_1/M)\mathbf{r}$. For $x\gg r$, a Taylor expansion givesThe mass dipole vanishes because ${\sum }_A{M}_A{\mathbf{x}}_A=0$, and the second mass moment tensor is ${\sum }_A{M}_A{x}_{Ai}{x}_{Aj}=\mu r_i{r}_j$. The gravitational quadrupole potential of a point-mass binary is consequentlyBoth printed tensors are trace-free. Contracting them with the Kronecker delta summed over three spatial dimensions givesThis proves the requested expression through second order. In particular the printed convention is $q_{ij}=\frac{3}{2}{Q}_{ij}$, where $Q_{ij}=\mu (r_i{r}_j-r^2{\delta }_{ij}/3)$ is the standard mass quadrupole moment; the radiation coefficient must be adjusted accordingly. No mass dipole term may be retained in centre of mass coordinates.

To evaluate the quadrupole formula without assuming a circular orbit, put $k=GM$, $u=\dot{r}=\mathbf{n}\cdot \mathbf{v}$ and ${\mathbf{v}}_{\perp }=\mathbf{v}-u\mathbf{n}$. Differentiating the gravitational acceleration gives the kinematic jerkFor $S_{ij}=r_i{r}_j$, the product rule then givesInserting these into the printed mass quadrupole moment convention yieldsThe two tensors in brackets are orthogonal in their tensor contraction, because $\mathbf{n}\cdot {\mathbf{v}}_{\perp }=0$. Their squared norms are $6$ and $2v_{\perp }^2$, respectively. HenceThe instantaneous quadrupole luminosity of a Kepler binary is thereforeThe radial derivative $\dot{r}$ is not the vector speed. The bracket equals $v_{\perp }^2+u^2/12$, so the radiated luminosity is nonnegative for every instantaneous velocity.

For a slowly evolving circular orbit, there is no preferred orbital phase at which to excite a persistent eccentricity vector. More precisely, the rotating mass quadrupole moment emits at twice the orbital frequency, with angular harmonic number two; its gravitational-wave energy and angular-momentum balance is $\dot{E}=\Omega \dot{J}$. The circular Kepler orbit sequence has $J=\mu \sqrt{GMa}$ and $dE/dJ=\Omega$, so this loss is tangent to the sequence. At leading adiabatic order it preserves zero secular orbital eccentricity. This is a circular gravitational-wave inspiral, with a small radial drift rather than an exactly fixed-radius Newtonian circle. The energy flux alone would not establish circularity without this symmetry and angular momentum balance.

To leading radiation order use $r=a$, $u=0$ and ${\mathbf{v}}^2=GM/a$ in the loss formula. Combining it with $E=-GM\mu /(2a)$ givesKeeping the fixed masses and leading quadrupole formula, integration givesThis formal coalescence time displays the strong $a_0^4$ dependence: gravitational-wave emission from a binary system matters far more for close binary stars than for wide ones. The orbital period decreases as $a^{3/2}$ and the gravitational-wave frequency increases, producing a chirp. A detached binary can be driven into Roche-lobe overflow; further evolution then depends on mass-transfer stability. A sufficiently close compact pair may merge. The weak-field, slow-motion, point-mass approximation ceases to apply before literal $a=0$; finite stellar radii or strong relativistic effects determine the final interaction.

Use $\dot{M}>0$ for inward mass accretion rate, and $T_m>0$ for magnetic torque adding angular momentum to the disk. Combining the continuity equation with conservation of angular momentum identifies the inward mass flux as $(dh/dr{)}^{-1}d\mathcal{G}/dr$. In a steady Keplerian accretion disk, this is the supplied constant $\dot{M}$, soThe boundary condition is $\mathcal{G}(r_{\mathrm{in}})=T_m$. Since $d\Omega /dr=-3\Omega /(2r)$, the viscous torque in an accretion disk is $\mathcal{G}=3\pi \bar{\nu }\Sigma h$. Hence the steady Keplerian accretion disk with an inner torque hasIn particular, $\Sigma (r_{\mathrm{in}})=\lambda \dot{M}/(3\pi \bar{\nu })$. For a physical steady solution with this boundary condition, $\lambda \ge 0$.

For $0\le \lambda \ll 1$, the surface density of a disk rises from a small inner value toward $\dot{M}/(3\pi \bar{\nu })$: this is the usual Keplerian accretion disk with an approximately zero-torque inner boundary condition. For $\lambda \gg 1$, the inner surface density of a disk is large and decreases approximately as $r^{-1/2}$, giving a torque-dominated accretion disk. This approximation applies where $r/r_{\mathrm{in}}\ll (\lambda -1{)}^2$; sufficiently far out, any fixed finite $\lambda$ again approaches the same constant asymptote.

The constant net inward angular-momentum flux isFor $\lambda >1$, net angular momentum moves outward. The strong-torque profile is therefore decretion disk-like in its stress-dominated structure, but the specified mass flux remains inward: a genuine net decretion disk requires an outward mass flux, not merely large $\lambda$.

At $r_{\mathrm{in}}=r_m$, the Maxwell stress tensor estimate givesUsing the magnetospheric truncation radius scaling in part (b), all dimensional parameters cancel:If the displayed $4\pi$ is retained and the scaling radius is assigned unit coefficient, the estimate is $1/(4\pi )$; retaining $\Sigma \simeq 2H\rho$ and exact free-fall speed throughout gives $\sqrt{2}/(4\pi )$. Those prefactors are not controlled by the scaling argument. In particular, it does not generate a parametrically large $\lambda \gg 1$.

In a steady state, mass conservation makes $\mathcal{F}$ independent of radius. Taking the accretion rate $\dot{M}>0$ for inward flow gives $\mathcal{F}=-\dot{M}$. The conservation of angular momentum equation then makes $-\dot{M}h+\mathcal{G}$ constant. The zero-torque inner boundary condition fixes this constant to $-\dot{M}{h}_{\mathrm{in}}$, soThe viscous torque in an accretion disk is $\mathcal{G}=-2\pi r^3\bar{\nu }\Sigma {\Omega }^{\prime }=2\pi qh\bar{\nu }\Sigma$. Combining these expressions gives the steady accretion disk with arbitrary rotation lawFor $x=r/R_g$, the Paczyński-Wiita circular orbit formulas giveConsequentlyIn particular, $f(6)=0$ and $f\to 1$ far from the black hole, recovering the outer Keplerian accretion disk.

Inside the innermost stable circular orbit, the gas enters the plunging region of a black-hole accretion disk. Its rapid inward motion leaves little time for stresses to exchange angular momentum, motivating the zero-torque inner boundary condition. This is a thin-disc approximation; a strong magnetic stress could change it.

Matter supplied from very large radius has negligible specific orbital energy, while matter crossing the inner edge carries ${\epsilon }_{\mathrm{in}}=-\eta c^2$. With zero inner torque, no energy is supplied by a stress at that edge. If the heat released outside it escapes by radiative transfer, rather than being lost through inward advection, the integrated conservation of energy balance gives

Let $\mathbf{r}={\mathbf{r}}_2-{\mathbf{r}}_1$ and $\mathbf{v}=\dot{\mathbf{r}}$. Relative to the centre of mass, the two positions are ${\mathbf{r}}_1=-(M_2/M)\mathbf{r}$ and ${\mathbf{r}}_2=(M_1/M)\mathbf{r}$. Their velocities have the same mass factors. Adding their kinetic energies and angular momenta gives, with reduced mass $\mu =M_1{M}_2/M$,These are the internal orbital quantities, excluding any uniform centre-of-mass motion. Newton's law of universal gravitation givesFor constant masses, differentiating $E$ makes the gravitational work cancel the derivative of the potential energy, so $\dot{E}=0$. Also $\dot{\mathbf{h}}=\mathbf{r}\times \ddot{\mathbf{r}}=0$, since the force is central; hence conservation of angular momentum gives $\dot{\mathbf{J}}=0$.

The eccentricity vector, the normalized attractive Laplace-Runge-Lenz vector, isUsing $\dot{\hat{\mathbf{r}}}=\mathbf{v}/r-\mathbf{r}(\mathbf{r}\cdot \mathbf{v})/r^3$ andshows that its derivative also vanishes. Adding a relative perturbing acceleration $\mathbf{f}$ changes these cancellations only by the perturbing terms. The perturbed Kepler-orbit conservation laws are therefore

For the distant Newtonian gravitational potential, expand each source term about the binary's centre of mass:The dipole vanishes because ${\sum }_i{M}_i{\mathbf{r}}_i=0$. With $\mathbf{a}=\mathbf{r}$, the quadratic mass moment is ${\sum }_i{M}_i{\mathbf{r}}_i{\mathbf{r}}_i=\mu \mathbf{aa}$, so the quadrupole potential of a circular binary isOn the inner binary's fast timescale, the distant companion's position is approximately fixed. Averaging over the inner circular orbit in the common plane gives $⟨{\cos }^2\theta ⟩=1/2$ and henceThis averaged potential function is independent of both time and azimuth. Consequently the companion's energy and axial angular momentum are conserved in the averaged problem. The exact motion has small fast variations; the statement is a secular approximation for a hierarchical system away from an orbital resonance.

For a circular outer orbit, force balance givesso the circumbinary orbital-period correction isConservative binary mass transfer leaves $M$ unchanged, so it does not change the leading monopole force. It changes $\mu$ and $a$, and therefore the quadrupole correction. If the inner orbital angular momentum is conserved, $J_{\mathrm{in}}=\mu \sqrt{GMa}$ gives $a\propto {\mu }^{-2}$ and $Q=\mu a^2\propto {\mu }^{-3}$. The correction is thus largest as the binary becomes strongly unequal in mass and widens, within the hierarchy $a\ll x$.

One must also allow the outer radius to respond. For slow axisymmetric evolution its specific angular momentum $h_3$ remains constant, while its energy need not: the averaged potential now depends on time. A nearly circular outer orbit has $h_3^2=GMx+3GQ/(4x)$ and $P_3=2\pi x^2/h_3$. At fixed $M,h_3$, differentiation givesThus mass exchange inside the binary can produce a small measurable outer-period change even without total mass loss. Nonconservative mass loss additionally changes the leading monopole term.

Work in the rest frame of the hole and neglect stellar velocities compared with $v$. In addition to $v\gg {\sigma }_{\ast }$, assume the galactic encounter speed is nonrelativistic, $v\ll c$. Treat a star as a test particle with specific angular momentum $\ell =bv$. At its pericentre $r_{\min }$, its radial velocity vanishes. Conservation of energy and conservation of angular momentum giveThe second term is gravitational focusing. For the rough Newtonian estimate in the question, assume capture when the pericentre reaches the Schwarzschild radius, $r_S=2G{M}_{\mathrm{BH}}/c^2$. It follows thatThe approximation keeps the focusing term because $v\ll c$.

Let ${\rho }_{\ast }=m{n}_{\ast }$ be the stellar mass density. Multiplying the geometric collision cross-section by the incoming mass flux gives the collisionless black-hole capture rateFor a supersonic gas with negligible sound speed, Bondi--Hoyle--Lyttleton accretion has capture radius $r_a=2G{M}_{\mathrm{BH}}/v^2$ and rateConsequentlyFor comparable densities and $v=200\mathrm{km}{\mathrm{s}}^{-1}$, the suppression is approximately $4\times 10^{-7}$. Gas can collide, dissipate orbital energy, and be captured from a region of size $GM/v^2$; collisionless particles generally leave on an unbound orbit unless their angular momentum is small enough for direct capture. Chandrasekhar dynamical friction can slow the hole by transferring energy to a stellar wake, but does not mean that the stars producing that wake are swallowed.

There is an important accuracy limitation in the requested numerical coefficient. Newtonian orbital mechanics evaluated at a relativistic event horizon cannot give the exact capture boundary. In Schwarzschild spacetime, a slowly incoming massive particle has dimensionless energy $\mathcal{E}\simeq 1$ and timelike geodesic effective potentialThe double turning point satisfying $V=1$ and $dV/dr=0$ occurs at $r=4r_g$ and $\ell =4G{M}_{\mathrm{BH}}/c$. The exact low-speed boundary is therefore $b_{\mathrm{cap}}\simeq 4G{M}_{\mathrm{BH}}/(cv)$, twice the rough value; it makes the rate ratio $4({\rho }_{\ast }/{\rho }_g)(v/c{)}^2$. The key suppression scaling is unchanged.

Direct capture of ordinary stars or dark matter is therefore an inefficient route to rapid growth of a supermassive black hole under these conditions. Tidal disruption events can capture stellar debris without swallowing an intact point-like star, and special dense environments can modify the supply. These qualifications do not make collisionless capture equivalent to dissipative gas accretion.

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.

Take $+x$ along the upper ship's initial velocity, and take $+y$ from the lower ship to the upper ship at closest approach. The upper and lower ships then have initial velocities $v_1\hat{\mathbf{x}}$ and $-v_2\hat{\mathbf{x}}$. Their equal masses put each at distance $l/2$ from the centre of mass, soBy momentum conservation, the centre of mass keeps moving with constant velocityIn that frame the ships have opposite velocities of magnitude $u=(v_1+v_2)/2$. Their initial angular momentum about the centre of mass is nonzero even though their trajectories were straight:The tether exerts central internal forces, so conservation of angular momentum givesEach ship follows uniform circular motion of radius $l/2$ around the translating centre of mass, clockwise viewed from $+z$. More explicitly, put the initial centre at the origin and defineThe kinetic-energy decomposition about the centre of mass verifies conservation of mechanical energy directly:At closest approach the relative radial velocity is already zero, so connecting the tether requires no dissipative radial impulse. Subsequently its tension is perpendicular to each ship's velocity relative to the centre of mass and does no net work on the pair at fixed $l$.