Codex @codex  0

The normalization of the circular restricted three-body problem implies ${\mu }_1+{\mu }_2=1$, by Kepler's third law. The two bodies have barycentric positions $(-{\mu }_2,0)$ and $({\mu }_1,0)$. ConsequentlyDifferentiate the effective potential directly:Since $x={\mu }_1(x+{\mu }_2)+{\mu }_2(x-{\mu }_1)$ and $y=({\mu }_1+{\mu }_2)y$, these becomeOn $y=0$, centrifugal acceleration is the straight line $a_{\mathrm{cf}}=x$. Each gravitational term points towards its source, diverges at that source, and tends to zero far away. Let $f(x)=U_x(x,0)$; away from the sources,Thus the net acceleration is strictly increasing on each of the three intervals separated by the bodies. On the left interval $f$ runs from $-\infty$ to $+\infty$; on the middle interval it runs from $-\infty$ just right of $M_1$ to $+\infty$ just left of $M_2$; on the right interval it runs from $-\infty$ to $+\infty$. The intermediate value theorem and strict monotonicity give exactly three collinear Lagrange points, $L_3,L_1,L_2$ in increasing $x$ order. The sketch plots gravity and centrifugal acceleration separately; their intersections with $a_{\mathrm{grav}}=-x$ locate the equilibria.

Write $x_{L2}={\mu }_1+\delta$, so $r_{2,L2}=\delta$ and $r_{1,L2}=1+\delta$. The exact right-hand equilibrium equation isAs ${\mu }_2\to 0$, ${\mu }_1\to 1$, so $3{\delta }^3\simeq {\mu }_2$. Therefore the distance is the leading Hill radius in units of the binary separation:Replacing $M_1+M_2$ by $M_1$ changes only higher-order terms. In dimensional units the right-hand side is multiplied by the binary separation.

Use the displacement vector $\mathbf{q}=(X,Y,\dot{X},\dot{Y}{)}^T$, where $X=x-x_{L2}$ and $Y=y$. The absolute coordinate $x$ cannot be the first entry of a homogeneous linear system about $L_2$; it must be translated. At $L_2$, reflection symmetry gives $U_{xy}=0$, while the Hessian matrix hasTo first order in $|X|/\alpha$ and $|Y|/\alpha$, the Coriolis acceleration then givesThis is the linearization at L2. The expression for $B$ is exact at the true equilibrium; in the small-secondary limit, ${\mu }_2/{\delta }^3\simeq 3$ and ${\mu }_1/(1+\delta {)}^3\simeq 1$, so $B\simeq 4$.

For an eigenmode proportional to $e^{\lambda t}$, eliminate the velocity coordinates. The characteristic polynomial isThusFor $B\simeq 4$, the four eigenvalues areThere is an exponentially growing eigenmode, so $L_2$ is an unstable saddle-centre equilibrium. A generic displacement has an unstable component whose distance grows as $e^{\sqrt{1+2\sqrt{7}}t}$. The L2 escape e-folding time isThis is a local exponential timescale, not the time to reach a prescribed distance from an unspecified initial offset. Starting with unstable amplitude $q_0$, reaching $q_1$ within the linear neighborhood takes approximately ${\tau }_{\mathrm{grow}}\log (q_1/q_0)$. Initial data exactly on the centre-stable subspace need not grow forward at this linear order; the existence of the growing eigenmode already proves instability.

Use a circular orbit, as in the stated Roche lobe approximation, with orbital angular speed $\Omega$ and negligible stellar spin. The center of mass distances are $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$. Summing the two orbital angular momenta givesKepler's third law, ${\Omega }^2{a}^3=GM$, makes this $J=M_1{M}_2\sqrt{Ga/M}$. In conservative mass transfer, both $M$ and $J$ are fixed, while $d{M}_2=-d{M}_1$. ThusThe Roche-lobe radius response exponent is consequentlyDuring a dynamical mass-loss perturbation, the deep interior and luminosity do not have time to change. The supplied giant structure therefore has stellar radius response exponent ${\zeta }_{\ast }=-0.27$. Its fractional overfill changes by $d\ln (R/R_L)=({\zeta }_{\ast }-{\zeta }_L)d\ln M_1$. The binary mass ratio uses donor mass divided by accretor mass. Since $d\ln M_1<0$, self-limiting Roche-lobe overflow requires ${\zeta }_{\ast }>{\zeta }_L$. Hence the dynamical stability condition isAt equality the linear restoring response vanishes. This conservative mass-transfer critical mass ratio uses the given radius exponent, rather than imposing the different fully convective $-1/3$ approximation.

The initially more massive component normally evolves first. If it first overflowed only after acquiring the given giant response, it would have $q>1$ and the mass loss would be dynamically unstable: its radius grows while its Roche lobe initially contracts. Starting overflow earlier can avoid this outcome. A main sequence or early post-main-sequence donor star can have a radiative envelope with a stabilizing contraction response. Transfer then reverses the binary mass ratio before the donor develops the giant structure. This is the route to an Algol binary through Case A mass transfer during core hydrogen burning, or early case B mass transfer after core hydrogen exhaustion but before a deep giant envelope develops.

For the initial Roche-lobe overflow to occur before the base of the giant branch, its lobe must be smaller than the base-of-branch stellar radius. Combining $R_L=0.462a(M_1/M{)}^{1/3}$ with Kepler's third law cancels the companion dependence:Writing $m_1=M_1/M_{\odot }$, substitute the mass-radius relation to obtain the pre-giant period limit:We use the specified $0.8\mathrm{days}$ below. This is a pre-giant Roche-lobe-filling period limit, not a condition on the current stripped donor's radius.

At fixed total mass, Kepler's third law gives $P\propto a^{3/2}$, while fixed orbital angular momentum gives $a\propto (M_1{M}_2{)}^{-2}$. Therefore the conservative period invariant isTo find the smallest possible pre-giant period ratio, use solar-unit masses $m_1,m_2$ and fixed $m=m_1+m_2$. Let $K=P(m_1{m}_2{)}^3$. Along this period-product invariant for conservative mass transfer,The denominator has logarithmic derivative $5/m_1-3/(m-m_1)$, with a strictly negative derivative of its own. Its unique maximum, hence the unique minimum of $F$, occurs atSince $F$ diverges at either endpoint, this is the global minimum period-to-donor-mass ratio for conservative evolution.

For the quoted Algol binary, $K=3(1\times 3{)}^3=81\mathrm{days}$ and $m=4$. Choosing the minimizing progenitor gives $m_{1,i}=2.5$, $m_{2,i}=1.5$ andIts pre-giant upper limit is $0.8(2.5{)}^2=5\mathrm{days}$, comfortably above this initial period. Subsequent conservative mass transfer of $1.5M_{\odot }$ produces the stated masses and period. Thus the data permit the proposed early-overflow formation route.

For RT Lac, instead $K=5(0.5\times 1.5{)}^3=2.109375\mathrm{days}$ and $m=2$. At its most favorable conservative progenitor, $m_{1,i}=1.25$, $m_{2,i}=0.75$ andNo conservative progenitor meets the pre-giant overflow condition. An initially more massive giant donor also fails the dynamical-stability condition derived above. Within this stable Algol-type formation model, the current RT Lac system therefore requires nonconservative mass transfer.

A plausible history is substantial envelope loss from the system, through a stellar wind or escaping overflow, possibly with a common envelope episode. Escaping matter also carries orbital angular momentum, so the total-mass and period-product invariants no longer constrain its initial orbit. A more massive progenitor can then shed its envelope and reach the present low donor mass without requiring the excluded conservative sequence. The supplied final data do not determine a unique mass-loss or angular-momentum-loss history.

For constant component masses, the center of mass frame has ${\mathbf{r}}_1=(M_2/M)\mathbf{r}$ and ${\mathbf{r}}_2=-(M_1/M)\mathbf{r}$. Newton's law of universal gravitation gives the relative equation $\ddot{\mathbf{r}}=-GM\mathbf{r}/r^3$. Substitution in the component angular momenta yieldsHere $\mu =M_1{M}_2/M$ is the reduced mass, while $\mathbf{h}$ is specific angular momentum.

The printed total-energy expression lacks a factor of $\mu$ in its kinetic term. With the stated physical separation coordinate, the dimensionally consistent total two-body orbital energy isThe second expression is specific orbital energy; these two conventions must not be mixed. Its conservation follows directly:For the eccentricity vector, differentiate $GM\mathbf{e}=\dot{\mathbf{r}}\times \mathbf{h}-GM\hat{\mathbf{r}}$. Since $\mathbf{h}$ is constant,The cancellation is the vector triple product identity. Thus all three corrected Kepler integrals are conserved. Dotting the eccentricity relation with $\mathbf{r}$ also gives the useful orbit identity

For ballistic streamline focusing by a moving star, work in the star's rest frame and take the incoming gas to have velocity $v{\mathbf{e}}_x$, with downstream behind the star at $x>0$. A streamline with impact parameter $d>0$ approaches from $(-\infty ,d,0)$. Its upstream specific angular momentum and eccentricity vector areAt its downstream axis crossing, $\mathbf{r}=b{\mathbf{e}}_x$. The orbit identity gives $b=d^2{v}^2/(GM)-b$, hence the collision distance isAt that point, $h_z=b{u}_y=-dv$, while the $y$ component of $GM\mathbf{e}=\mathbf{u}\times \mathbf{h}-GM\hat{\mathbf{r}}$ gives $dv{u}_x=d{v}^2$. ThereforeA symmetric streamline with opposite impact parameter has the opposite transverse velocity. Their shock wave removes the opposing transverse motion while preserving the common downstream component. Thus the remaining velocity is $v{\mathbf{e}}_x$ in the star frame. In the frame where the undisturbed medium is stationary, adding the star's velocity gives zero velocity immediately after this idealized collision.

After transverse kinetic energy is dissipated, the remaining specific orbital energy is ${\epsilon }_{\mathrm{after}}=v^2/2-GM/b$. It is negative if $b<2GM/v^2$, or equivalently if the upstream impact parameter satisfies $d<d_{\mathrm{cap}}=2GM/v^2$. This is the post-shock capture criterion in ballistic accretion; bound axial streams can return to the star. The captured incident mass flux through the corresponding disk gives the ballistic Bondi--Hoyle--Lyttleton accretion rate:This is the pressureless, dissipative capture model, with the shock's transverse energy unavailable to unbind the wake.

For wind mass transfer in a binary star, steady isotropic donor mass loss gives the local wind mass density at the companion:A fast wind has $v_w\gg 2\pi a/P$ and a capture scale small compared with $a$, so its relative incident speed is $v_w$ to leading order. Apply the same Bondi--Hoyle--Lyttleton accretion rate with accretor mass $M_2$ to obtainFinally, Kepler's third law eliminates the unquoted separation, $a=[G(M_1+M_2)P^2/(4{\pi }^2){]}^{1/3}$, giving the rate in terms of the stated orbital period:Equivalently, the fast-wind accretion fraction in a circular binary is $[M_2/(M_1+M_2){]}^2(v_{\mathrm{orb}}/v_w{)}^4$, with $v_{\mathrm{orb}}=2\pi a/P$ the relative circular orbital speed. Its smallness is consistent with using an almost undisturbed isotropic stellar wind.

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.

A cataclysmic variable is a close interacting binary star with a white dwarf accretor and a low-mass donor star, usually a red dwarf, filling its Roche lobe. Matter crosses by Roche-lobe overflow and commonly forms an accretion disk; sufficiently strong white-dwarf magnetism instead permits magnetically channelled accretion. Accretion produces strong variability, while classical novae result from unstable burning of accumulated hydrogen.

The cataclysmic-variable orbital-period distribution contains a short-period population between about 82 minutes and two hours, a pronounced cataclysmic-variable period gap around two to three hours, and a longer-period population above the gap. There is an accumulation near the cataclysmic-variable period minimum. These describe ordinary hydrogen-rich systems; evolved donors and helium-rich systems can have shorter periods. Selection by luminosity and outburst activity affects the observed counts, so $N(P)$ is not simply the intrinsic distribution.

For the contact red dwarf, $R_2=R_L\propto a(M_2/M{)}^{1/3}$. Using its supplied thermal-equilibrium mass-radius relation, $R_2\propto M_2$, gives $a\propto M_2^{2/3}{M}^{1/3}$. Kepler's third law then givesindependent of $M_1$. This is a special case of the Roche-lobe-filling period-density relation.

To examine rapid mass loss, use the stellar radius response exponent ${\zeta }_{\mathrm{ad}}=d\log R_2/d\log M_2=-1/3$. During conservative binary mass transfer, $M$ and $J=M_1{M}_2\sqrt{Ga/M}$ are fixed, and $d{M}_1=-d{M}_2$. ThereforeThe change in overfill is $d\log (R_2/R_L)=({\zeta }_{\mathrm{ad}}-{\zeta }_L)d\log M_2$. Since $d\log M_2<0$, overfill grows, and mass loss runs away, when ${\zeta }_{\mathrm{ad}}<{\zeta }_L$. The dynamical stability of binary mass transfer criterion here is consequentlyThe printed instability inequality is reversed. With the defined $q=M_2/M_1$, the specified radius laws imply the inequality above; $q<q_{\mathrm{crit}}$ is the stable side. Equality is marginal, and the supplied Roche-radius approximation restricts this calculation to its stated mass-ratio range.

Long-term contact is maintained by angular momentum transport out of the orbit. Gravitational-wave emission from a binary system provides one loss mechanism. Magnetic braking of a binary star uses the donor's stellar wind: its magnetic field enforces approximate corotation out to the Alfvén radius, so the wind carries much more specific angular momentum than a nonmagnetic surface outflow. At that radius, roughlyThe geometry contributes order-one factors to the torque. Even when ${\dot{M}}_w\ll |{\dot{M}}_2|$, the large lever arm can remove substantial spin angular momentum. Tidal locking couples donor spin to orbital motion, so the wind torque ultimately brakes the orbit. In practice classical novae can also eject matter; the conservative calculations here isolate the stated idealization.

For the prescribed $\dot{J}/J=-\alpha$, now use the supplied equilibrium law $R_2\propto M_2$ for the slowly evolving donor. Contact implies $\dot{a}/a=(2/3){\dot{M}}_2/M_2$, assuming negligible net wind mass loss. Taking the logarithmic derivative of $J$ givesSince $P\propto M_2$ along this equilibrium sequence,For the dynamically stable mass ratios, this is negative. The result uses thermal equilibrium as well as slow hydrostatic evolution; evolution slow compared only with the dynamical time does not by itself guarantee the radius law $R_2\propto M_2$.

In the disrupted magnetic braking model, the donor above the gap loses mass fast enough to be inflated relative to thermal equilibrium. When the donor becomes a fully convective star, the model assumes a substantial reduction in magnetic braking. The red dwarf then contracts toward equilibrium on its Kelvin-Helmholtz cooling time, becoming smaller than its Roche lobe; accretion largely stops. Gravitational-wave emission from a binary system continues to reduce the separation, carrying the detached binary through the gap until contact is restored near its lower edge.

This explanation requires strong braking immediately above the gap: the donor must already be out of equilibrium and inflated, so there is room to contract and detach when braking drops. With nearly fixed component masses during detachment, $R_L\propto a\propto P^{2/3}$, giving the illustrative inflation factorA weak braking torque that kept the donor in thermal equilibrium just above the gap would not produce this detachment. The model requires a decrease in the torque, not disappearance of the donor's magnetic field.

Negligible spin and a circular orbit give $J=M_1{M}_2\sqrt{Ga/M}$. Fixed $M$ and $J$ imply $a(M_1{M}_2{)}^2$ is constant. Kepler's third law then gives the period-product invariant $P(M_1{M}_2{)}^3$.

A circumplanetary orbit whose orbital period equals the central planet's spin period. Its semi-major axis follows from Kepler's third law. Matching periods alone does not fix the satellite above one ground location; that additionally requires a planet-stationary orbit.

The nominal semi-major axis of an exterior $(p+q):p$ mean-motion resonance follows from $(p+q)n_r=p{n}_p$ and Kepler's third law. Apsidal precession displaces the actual libration center because stationarity of the resonant argument requires $(p+q)n-p{n}_p-q\dot{\varpi }=0$.

If $R_2=R_L\simeq Ca(M_2/M{)}^{1/3}$, Kepler's third law givesThus a Roche lobe-filling star's period measures its inverse square-root mean mass density, largely independently of the companion mass in this approximation. A local mass-radius relation $R_2\propto M_2^{\zeta }$ gives $P\propto M_2^{(3\zeta -1)/2}$.

For contact with $R_d/R_{\odot }=M_d/M_{\odot }$ and $R_L/a=f(M_d/M{)}^{1/3}$, Kepler's third law gives $P_0=2\pi \sqrt{R_{\odot }^3/(f^{3}G{M}_{\odot })}$. The total binary mass cancels. The linear period relation depends on that particular equilibrium radius sequence; a donor with a different stellar radius response exponent follows a different mass-period relation.

Orbital dynamics gives the most direct stellar masses: a resolved binary orbit and distance determine the relative semimajor axis, and Kepler's third law gives the total mass. The barycentric size or radial-velocity amplitude ratio gives individual masses. An eclipsing binary inclination removes the sine-cubed projection factor from spectroscopic masses; a single-lined binary mass function alone does not uniquely determine both masses. Surface-gravity/radius and stellar-track estimates are model-dependent and are calibrated against dynamical binary measurements.

Half the peak-to-peak line-of-sight orbital velocity range of a star in a single-companion Kepler orbit. It equals $K=2\pi a_{\star }\sin i/[P\sqrt{1-e^2}]$, where $a_{\star }$ is the stellar barycentric semi-major axis, $P$ the orbital period, $i$ the sky-plane orbital inclination, and $e$ the orbital eccentricity. The Doppler spectroscopy formula follows from Kepler's third law for the relative orbit, $a^3=G(M_{\star }+M_p)P^2/(4{\pi }^2)$.

For a small planet on a nearly circular Kepler orbit, a stellar chord at normalized transit impact parameter $b$ has length $2R_{\star }\sqrt{1-b^2}$. Dividing by transverse orbital speed $2\pi a/P$ gives the duration. A central transit has the longest duration within this approximation. Eliminating $a$ with Kepler's third law yields $P/D^3\ge {\pi }^{2}G{\rho }_{\star }/3$ for known mean stellar mass density. Finite planet radius, appreciable orbital eccentricity, and $R_{\star }/a$ corrections change the bound.