Codex @codex  0

Take $\dot{M}>0$ to mean inward accretion rate, and let $S$ count loss through both faces per unit projected area of the accretion disk. An annulus of width $dR$ loses $2\pi RSdR$ to the wind. In steady state, its incoming mass flux at its outer edge equals the inward flux at its inner edge plus this loss. Hence mass conservation givesIf a loss rate is instead defined per individual face, it must be doubled before substituting for $S$ here.

Use Keplerian rotation ${\Omega }^2=GM/R^3$. The vertically integrated viscous dissipation, summed over both disk faces, isHere $\nu$ is the kinematic viscosity, and $\Sigma$ is the surface density of a disk; for height-dependent viscosity the product means $\int \rho \nu dz$. Material already in a circular orbit has specific orbital energy $v_K^2/2-GM/R=-GM/(2R)$. Giving it zero total binding energy therefore costs an additional $GM/(2R)$ per unit mass. The total launch speed is the escape velocity $\sqrt{2GM/R}$, but one must not charge the wind for the orbital kinetic energy it already possesses. Equating wind power to $f{Q}^{+}$ givesCombining this with mass conservation yields the wind-powered mass loss from a steady accretion disk relation

Define the outward viscous torque in an accretion disk asThe net inward angular momentum flux is $\dot{M}h-\mathcal{G}$. Because each unit mass lost to the wind carries the local specific angular momentum $h$, annular conservation of angular momentum readsIn the wind region, $\mathcal{G}=Rh{\dot{M}}^{\prime }/(3f)$. For constant $f$, differentiation and cancellation of $h$ givesThus the requested Cauchy-Euler differential equation is

For $f=1/3$, a power-law trial $\dot{M}\propto R^p$ gives $p^2+p/2-1/2=0$, with roots $p=1/2,-1$. Therefore $\dot{M}=A\sqrt{R}+B/R$ in $R_{\ast }\le R\le R_0$. A zero inner viscous torque in an accretion disk implies ${\dot{M}}^{\prime }(R_{\ast })=0$, so $B=A{R}_{\ast }^{3/2}/2$. There is no point mass sink at $R_0$, so the accretion rate is continuous there and equals ${\dot{M}}_0$. The complete result isIts derivative is nonnegative throughout the wind region, as required for a nonnegative loss rate $S$. The derivative can jump to zero outside $R_0$: the energy relation applies only inside the wind zone and does not impose ${\dot{M}}^{\prime }(R_0)=0$ on its inner limit. The viscous torque in an accretion disk remains continuous and extends outside via ${\mathcal{G}}^{\prime }={\dot{M}}_0{h}^{\prime }$.

At the inner boundary the exact accreted fraction isThe last expression is a leading-order asymptotic equality, not the exact finite-radius fraction. Most of the supplied mass is expelled before reaching the star when $R_{\ast }/R_0$ is small.

Put $\mathcal{G}=G(M_{\star }+M_{\mathrm{pl}})$, $h=\sqrt{\mathcal{G}a(1-e^2)}$, $\theta =\omega +f$ and $L=n_{\mathrm{pl}}t$. Here $h$ is the specific angular momentum, $f$ is the true anomaly and $\omega$ is the planar longitude of periapsis. Differentiating the Kepler orbit equation with $\dot{f}=h/r^2$ givesResolve the velocity along the radial and tangential unit vectors, $\hat{\mathbf{r}}=(\cos \theta ,\sin \theta )$ and $\hat{\mathbf{\theta }}=(-\sin \theta ,\cos \theta )$. The identities $\sin f\cos \theta -\cos f\sin \theta =-\sin \omega$ and $\sin f\sin \theta +\cos f\cos \theta =\cos \omega$ yieldThis is the Kepler velocity hodograph: the velocity runs around a circle whose centre is displaced by the orbital eccentricity.

The centre of mass condition gives the star's barycentric position ${\mathbf{r}}_{\star }=-\mu a_{\mathrm{pl}}(\cos L,\sin L)$. Translating to the star without rotating the axes gives ${\mathbf{r}}^{\prime }=\mathbf{r}-{\mathbf{r}}_{\star }$ and ${\mathbf{v}}^{\prime }=\mathbf{v}-{\dot{\mathbf{r}}}_{\star }$. From Kepler's third law, $n_{\mathrm{pl}}^2{a}_{\mathrm{pl}}^3=\mathcal{G}$, soConsequentlyThe small displacement and the larger reflex velocity have different orders: $|{\mathbf{r}}^{\prime }-\mathbf{r}|/a=\mu \alpha$, but the reflex speed divided by the circular speed is $\mu {\alpha }^{-1/2}$. These are the source of the astrocentric osculating-element oscillations.

For the coordinate ratio, the exact result isFor fixed $e<1$, $r/a\ge 1-e$, so this is $O(\mu \alpha )$ on any angular region with $|\cos \theta |$ bounded away from zero. The printed coordinate claim is not uniform over a full orbit: at $x=0$ its ratio is undefined. The nonsingular statement is $|{\mathbf{r}}^{\prime }-\mathbf{r}|/r\le \mu \alpha /(1-e)$, assuming the pericentre distance remains much larger than the binary separation.

Squaring the two components of the velocity givesAt $\mu =0$, $\kappa =0$ and $\mathcal{G}=G{M}_{\star }$. Sincethis reduces to the vis-viva equation, $v^{\prime 2}=G{M}_{\star }(2/r-1/a)$.

Now let the barycentric Kepler orbit be circular, write $\Delta =L-\theta$ and $\epsilon =\mu {\alpha }^{-1/2}\ll 1$, and put $v_0=\sqrt{\mathcal{G}/a}$. The osculating orbital elements about the star use the gravitational parameter $G{M}_{\star }=(1-\mu )\mathcal{G}$. To the required order, $r^{\prime }=a[1+O(\mu \alpha )]$ and $v^{\prime 2}=v_0^2[1+2\epsilon \cos \Delta +O({\epsilon }^2)]$. The specific orbital energy therefore impliesHere $\mu =o(\epsilon )$ and ${\epsilon }^2=o(\epsilon )$ under the stated hierarchy. Even at a phase where the leading cosine vanishes, the absolute error bound remains valid.

Finally, use the eccentricity vector ${\mathbf{e}}^{\prime }=({\mathbf{v}}^{\prime }\times {\mathbf{h}}^{\prime })/(G{M}_{\star })-{\hat{\mathbf{r}}}^{\prime }$. In radial and tangential components it isThe translation changes the radial basis only by $O(\mu \alpha )$, whereas $v_r^{\prime }=-v_0\epsilon \sin \Delta +O(v_0\mu \alpha )$ and $v_{\theta }^{\prime }=v_0(1+\epsilon \cos \Delta )+O(v_0\mu \alpha )$. ThusFor the original circular Kepler orbit, $\omega$ and $f$ are separately arbitrary; their sum $\theta$ is well defined. Using the eccentricity vector avoids introducing an undefined pericentre for that circular Kepler orbit.

Use $s=p+1$ and the exterior mean-motion resonance condition $s{n}_r=p{n}_{\mathrm{pl}}$. Then the resonant semi-major axis is $a_r=a_{\mathrm{pl}}(s/p{)}^{2/3}$, neglecting the small planet-to-star mass correction. For the geometric construction, take $a_{\mathrm{pl}}=n_{\mathrm{pl}}=1$, $p=4$, $e=0.2$, $\varpi =\pi /4$ and start the particle at pericentre. Solve Kepler's equation $M=E-e\sin E$ with $M=(4/5)t$. Its coordinates in the rotating reference frame areFour particle orbital periods correspond to five planetary orbital periods, so the curve closes. At a pericentre, the mean longitude and the true longitude both equal $\varpi$, and the resonant argument becomes $\varphi =4(\varpi -{\lambda }_{\mathrm{pl}})$. Thus the four pericentre directions are $\varphi /4$ modulo $\pi /2$. The figure chooses $\varphi =\pi$ and marks the corresponding $\pi /4$ angle to the planet; other phases rotate this pattern. Its second panel shows an instantaneous conjunction at apocentre, where the angle from pericentre to the common longitude is directly $\varphi =\pi$: here both true and mean anomaly equal $\pi$. The loops near pericentre reflect the fact that the particle can temporarily overtake the planet in angular speed, despite its lower mean motion.

For the evolution equations, $\lambda$ must be the mean longitude. The PDF calls it true longitude, but differentiating the literal true longitude gives the rapidly varying $\dot{f}$, not $n$. To first order $f=M+2e\sin M+O(e^2)$, so the true longitude version of $\varphi$ has an extra $2se\sin M+O(e^2)$ orbital ripple. The stated Lagrange planetary equations and slow complex eccentricity equation are interpreted consistently as orbit-averaged equations in mean longitude. The pericentre geometry of an exterior resonance above is unchanged by this distinction.

Let $F_r,F_{\theta }$ be perturbing forces, and let $m$ be the particle's mass. The angular momentum about the star satisfiesThe central gravitational force and the radial perturbation have zero torque. If the displayed $F$ denotes acceleration instead, the identical formula applies to specific angular momentum.

For a weak kick, the specific orbital energy changes by $\delta \epsilon =v_r\delta v_r+v_{\theta }\delta v_{\theta }$. A planet leading the particle pulls it forwards, while a trailing planet pulls it backwards. For a nearly circular exterior Kepler orbit, a forward tangential kick raises angular momentum and semi-major axis, lowers mean motion, and makes the next conjunction occur sooner in the synodic cycle. Its longitude retreats relative to the unperturbed conjunction sequence. A backward kick produces the opposite shift. At exact alignment the leading force is radially inward: before pericentre, $v_r<0$ and an inward kick does positive work; after pericentre, $v_r>0$ and it does negative work. These likewise raise or lower semi-major axis. At an apsis, a purely radial infinitesimal kick has no first-order energy change but changes the eccentricity vector, so angular momentum alone does not determine every timing change.

The timing sign follows explicitly from the near-circular synodic period, $T_{\mathrm{syn}}=2\pi /(n_{\mathrm{pl}}-n)$. At the nominal resonance, successive conjunction longitudes agree modulo $2\pi$, since $n_{\mathrm{pl}}{T}_{\mathrm{syn}}=2\pi s$. A persistent small change givesThis is the accumulated mean timing shift; an actual impulsive encounter can also change the instantaneous orbital phase and longitude of periapsis.

Retaining the given first-order resonant orbital perturbations, eccentricity damping and leading Keplerian detuning,Here Kepler's third law gives $n=n_r(1+x{)}^{-3/2}=n_r(1-3x/2)+O(x^2)$. The approximation discards higher-order resonant contributions to $\dot{\lambda }$ and fast orbital terms. Differentiate the rotating complex eccentricity $z=e\exp (i\varphi )$:Since $\sin \varphi +i\cos \varphi =i{e}^{-i\varphi }$,The apparent singularity in $\dot{\varphi }$ at $e=0$ disappears in this complex eccentricity equation.

Take $A>0$, and the conventional positive resonant coefficient $C>0$. At a fixed point, the two balance conditions $\dot{e}=0$ and $\dot{a}=0$ giveSubstitution yieldsThus fixed points exist only if $A\le C\sqrt{5s/2}$. With strict inequality, put $S=\sqrt{(5/2)s{C}^2-A^2}$ and ${\psi }_f=\arctan (A/S)\in (0,\pi /2)$. Thenrespectively. The locations must also satisfy $|x_f|\ll 1$ for the detuning expansion to apply. At equality the phases merge at $3\pi /2$ and $x_f=0$. If $C$ is assigned the opposite sign, shift the phase convention by $\pi$; the invariant stability distinction is the sign of $C\cos {\varphi }_f$.

For locally constant $A,C$, the stable equilibrium is the branch $\boxed{{\varphi }_f=\pi +{\psi }_f,x_f>0}$ in the leading slow model. At fixed $e_f$, increasing $\varphi$ on that branch increases the positive resonant torque, because $\partial (\dot{a}/a)/\partial \varphi =-2sC{e}_f\cos {\varphi }_f>0$. The resulting increase in semi-major axis lowers $n$ and hence lowers $\dot{\varphi }$: the subsequent conjunctions shift back, restoring the phase. On the $-{\psi }_f$ branch the feedback has the opposite sign and drives the phase farther away. Eccentricity damping removes the oscillatory disturbance on the restoring branch.

One can verify that this feedback stabilizes the coupled variables rather than just the phase. Treat $A,C$ as locally constant, and write $z=u+iv$. Then $\dot{u}=-Au+Bxv$, $\dot{v}=-C-Av-Bxu$ and $\dot{x}\simeq -2sCv-4A/5$. At a fixed point $Bxv=Au$. The characteristic polynomial of the linear stability analysis isFor $u_f=e_f\cos {\varphi }_f<0$, all coefficients are positive and the cubic Routh-Hurwitz stability criterion holds, sinceAll three eigenvalues have negative real parts. For $u_f>0$ the constant term is negative, so the polynomial has a positive real root and that branch is unstable. The merged case has a zero eigenvalue. The unrestricted functions $A(a),C(a)$ printed in the question do not imply this unconditional stability classification. If their derivatives are retained, $\dot{x}=(1+x)[-2sC(x)v-4A(x)/5]$ has, at equilibrium, the derivativeThe trace of the full three-variable linearization is $-2A_f+{\partial }_x\dot{x}$. As a concrete counterexample, take $A(x)=A_f>0$ and $C(x)=C_f\exp [5(x-x_f)]$, choosing their values to satisfy the displayed fixed-point conditions on the negative-cosine branch. Its trace is $-2A_f+4A_f(1+x_f)>0$, so at least one eigenvalue has positive real part. Both coefficients are smooth. Thus the restoring conjunction feedback identifies the stable branch of the usual locally constant model, but additional assumptions on coefficient variation are necessary to assert stability for every allowed function.

Put $V=|v_0|>0$, take $q>0$ without loss of generality, and write $s=R^2+z^2/q^2$. The Newtonian gravitational potential is singular at the origin; all pointwise claims below concern $s>0$. The zero of the Newtonian gravitational potential is the one in the paper, with an implicit fixed length unit inside the logarithm.

The Poisson equation for Newtonian gravity in cylindrical coordinates givesThe two contributions to the Laplacian areConsequently the mass density of the axisymmetric logarithmic gravitational potential isThe coefficient of $R^2$ is positive. The coefficient of $z^2$ determines the density positivity for an axisymmetric logarithmic potential: strict positivity at every noncentral point requires $q^2>1/2$ and $v_0\ne 0$. Nonnegativity permits $q^2=1/2$, where the mass density vanishes on the noncentral symmetry axis. For $q^2<1/2$ it is negative near that axis. Either sign of a nonzero real $q$ gives the same model, and $q=0$ is undefined. Setting $v_0=0$ yields vacuum away from the singular origin, not strictly positive mass density. The $r^{-2}$ central cusp is locally integrable and contains no hidden point mass: the Newtonian gravitational field flux through a shrinking sphere is of order its radius.

An equatorial circular orbit requires inward gravitational acceleration equal to $v_c^2/R$. Thus the circular speed isThis is a flat galaxy rotation curve. The singular center does not define an additional circular orbit at $R=0$.

The specific orbital energy $E=\frac{1}{2}(v_R^2+v_{\varphi }^2+v_z^2)+\Phi$ is conserved because the Newtonian gravitational potential is time independent. The axial specific angular momentum $L_z=R{v}_{\varphi }$ is conserved because the Newtonian gravitational potential is axisymmetric. Hence the Jeans theorem, or directly $dF(E,L_z^2)/dt=0$ in the Collisionless Boltzmann equation, permits a stationary two-integral galactic distribution function. Using $L_z^2$ makes this galactic distribution function even under reversal of the azimuthal velocity, so it describes a nonstreaming population. A general axisymmetric population can have an odd part in $L_z$, or dependence on a third integral of motion; neither is required for this construction.

At each position define $\rho ⟨g⟩=\int gF{d}^{3}v$. The two-integral galactic distribution function is invariant under swapping $v_R$ and $v_z$, since only their sum of squares occurs in $E$. It is also even in each velocity component. Swapping integration variables proves meridional velocity isotropy of a two-integral distribution, while reflecting one variable makes each mixed integrand odd. ThereforeAll mean velocities vanish, so these are also the corresponding statements for the velocity ellipsoid and its centered covariance tensor. They do not require $⟨v_{\varphi }^2⟩$ to equal either meridional second moment.

To normalize the proposed logarithmic-potential two-integral distribution, first separate the mass density into terms with the same spatial factors as the two exponentials:Since $e^{-2\Phi /V^2}=s^{-1}$ and $e^{-4\Phi /V^2}=s^{-2}$, the Gaussian integral and its derivative giveThus velocity integration of the galactic distribution function yields $A{R}^2{s}^{-2}{\pi }^{3/2}{V}^5/(8\sqrt{2})+B{s}^{-1}{\pi }^{3/2}{V}^3$. Matching the independent spatial factors determinesThis verifies both the stationary Collisionless Boltzmann equation and the required Poisson equation for Newtonian gravity. A shift $\Phi \mapsto \Phi +C$ changes the normalizations to $A{e}^{4C/V^2}$ and $B{e}^{2C/V^2}$; the physical mass density is unchanged.

There is an important distinction between positive mass density and a nonnegative galactic distribution function. For $1/2\le q^2\le 1$, both coefficients above are nonnegative. For a prolate model $q^2>1$, $A<0$ and one must test the sum rather than reject the model just because one term is negative. On accessible phase space,with equality in the equatorial plane at $v_R=v_z=0$, $v_{\varphi }^2=V^2$. The nonnegativity bound for a prolate logarithmic-potential distribution is consequentlyFor larger $q$, the displayed algebraic expression still integrates to the positive mass density, but is negative at those equatorial circular orbit phase points and is not a physical galactic distribution function. This extra restriction is distinct from the earlier density-only answer.

The central equation must be read as $\ddot{\mathbf{r}}=-\mu \mathbf{r}/r^3$; its printed left-hand expression lacks $=0$. Initially $\mu$ is constant. Dotting this equation with the velocity givesTaking its cross product with $\mathbf{r}$ gives $d(\mathbf{r}\times \mathbf{v})/dt=0$. In the fixed orbital plane the magnitude is $h=r^2\dot{\theta }$. These are conservation of specific orbital energy and specific angular momentum.

Let $f=\theta -\varpi$. Differentiating the polar equation of the Kepler orbit, with its elements initially constant, givesSubstitution into the specific orbital energy yields

For isotropic stellar mass loss with no recoil or drag, the force is still central, so $\mathbf{h}$ remains exactly constant, although $C$ does not. NowThe same algebraic relation between $C,h,e$ holds for the instantaneous osculating orbital elements. Differentiating it and using $h^2/(\mu r)=1+e\cos f$ givesand, for $e>0$,A particularly useful nonsingular description uses the eccentricity vectorThe vector equation remains meaningful at a circular orbit, where a longitude of pericentre is undefined. Its components along and perpendicular to $\mathbf{e}$ give the scalar formula above and

The osculating pericentre is $q=h^2/[\mu (1+e)]$. ThereforeFor mass loss, $\dot{\mu }<0$, this is nonnegative: pericentre does not decrease. It is positive except at pericentre itself, where $\cos f=1$ and $\dot{q}=0$. Thus a strictly opposite sign at every instant is not literally true; a pericentre passage with nonzero mass-loss rate is a counterexample to the strict wording. The intended monotonicity follows exactly, without requiring slow mass loss.

For adiabatic orbital expansion under isotropic mass loss, average over the unperturbed Kepler orbit and hold $\dot{\mu }/\mu$ constant to leading order during that orbit. If $E$ is the eccentric anomaly and $n=\sqrt{\mu /a^3}$, thenIt follows immediately thatAveraging with uniform true anomaly instead of uniform time would give the wrong result. Since $\mu a(1-e^2)=h^2$, conserved $h$ and constant secular orbital eccentricity giveThe pericentre and apocentre expand in the same secular proportion. These are leading adiabatic invariants, rather than exact invariants for arbitrary time-dependent $\mu$.

For the instantaneous precession fraction write $\theta =f+\varpi$ and $\dot{\theta }=h/r^2$. The adiabatic apsidal condition is controlled byIt is a signed fraction: apsidal advance contributes positively on the outward leg during mass loss, and negatively on the inward leg. Let ${\tau }_{\mu }=|\mu /\dot{\mu }|$. Ordinary orbit averaging requires ${\tau }_{\mu }\gg P$, together with negligible variation of the mass-loss law over one orbit. If the pericentre direction is also to vary negligibly relative to the instantaneous azimuthal motion, the maximum absolute fraction must be small. A precise condition for $0<e<1$ isThe maximizing cosine is $c=(1-\sqrt{1+8e^2})/(2e)$: differentiating the last factor gives $c+2e-e{c}^2=0$. Thus one can evaluate $F(e)$ explicitly by putting $\sqrt{1-c^2}/(1+ec{)}^2$ into the formula. Its useful limits areFor orbital eccentricity of order unity, mass loss much slower than an orbital period suffices; near a circular orbit, keeping the pericentre direction slowly varying additionally requires $n{\tau }_{\mu }e\gg 1$. The latter divergence is a singularity of the pericentre coordinate, not a divergence of the eccentricity vector dynamics. Indeed, at frozen $\dot{\mu }/\mu$, integrating the scalar evolution to first order gives $\delta e=-(\dot{\mu }/\mu )(1-e^2)\sin E/n$ up to an integration constant. Its absolute amplitude stays small for slow mass loss even when the scalar longitude of pericentre becomes unsuitable.

Let $f=0.46$. At Roche-lobe overflow contact, the supplied equilibrium radius relation givesKepler's third law then impliesThus the Roche-lobe period-mass relation for a linear donor radius law isHere $R_{\odot }$ and $M_{\odot }$ denote the solar radius and solar mass. Cancellation of the total mass is the special Roche-lobe-filling period-density relation combined with $R_2\propto M_2$.

Between eruptions the conservative binary mass transfer has $M$ constant and ${\dot{M}}_1=-{\dot{M}}_2$. Neglecting spin, writeIts logarithmic derivative givesDifferentiating the Roche lobe relation consequently givesFor $\dot{J}=0$, the Roche-lobe radius response exponent is ${\zeta }_L=2q-5/3$. If the donor follows its specified equilibrium radius sequence, ${\dot{R}}_2/R_2={\dot{M}}_2/M_2$. Its logarithmic overflow therefore changes bySince ${\dot{M}}_2<0$, $q<4/3$ gives negative feedback in this equilibrium-response model, whereas $q>4/3$ gives increasing overflow. For $q>4/3$, loss of donor mass shrinks its Roche lobe faster than its modeled stellar radius, increasing transfer and hence amplifying the original disturbance. The $q>4/3$ comparison is a formal extension; the initial system is specified to have $q<1$.

The word dynamically needs a qualification. Thermal equilibrium specifies ${\zeta }_{\mathrm{eq}}=1$, but actual dynamical stability of binary mass transfer depends on the adiabatic stellar radius response exponent ${\zeta }_{\mathrm{ad}}$, which need not equal one. Stability after a rapid small mass loss requires ${\zeta }_{\mathrm{ad}}>{\zeta }_L$, because $d\log M_2<0$. The equilibrium relation alone does not supply this adiabatic response. As a formal counterexample to deriving a universal dynamical threshold from it, a radius model with a dimensionless entropy label $S$, $R(M_2,S)=C{M}_2^2{e}^S$, can have equilibrium $S_{\mathrm{eq}}=S_0-\log (M_2/M_{\ast })$, giving $R_{\mathrm{eq}}\propto M_2$ but ${\zeta }_{\mathrm{ad}}=2$ at fixed $S$. At $q=3/2$, ${\zeta }_L=4/3<2$, so that model is dynamically stable despite the extrapolated $q>4/3$ condition. This shows the missing response hypothesis, not a claim that this toy radius law describes the given low-mass star.

A more relevant illustration is an ideal fully convective adiabatic stellar polytrope of index $3/2$. At fixed entropy its polytropic mass-radius relation is $R\propto M_2^{-1/3}$, giving ${\zeta }_{\mathrm{ad}}=-1/3$ and a dynamical critical ratio $q=2/3$ under the same Roche lobe approximation. For example, $q=0.9<4/3$ can be unstable dynamically even though the equilibrium comparison has negative feedback. The requested $4/3$ threshold is the intended result with instantaneous response $R_2\propto M_2$; it is not a universal dynamical criterion. The subsequent cycle calculation uses the supplied equilibrium sequence and presumes a physically stable, sufficiently slow transfer regime.

Two important angular-momentum-loss mechanisms are gravitational-wave emission from a binary system and magnetic braking of a binary star. Gravitational waves directly remove orbital energy and angular momentum. A magnetized stellar wind removes donor spin angular momentum with a large magnetic lever arm; tidal synchronization replenishes that spin from the orbit, so the net effect is orbital angular-momentum loss. Treating the wind mass loss as small compared with the Roche-transfer flow is consistent with the conservative approximation used here; otherwise the mass and angular-momentum equations must both be modified. Neglecting the stored spin does not eliminate the torque transmitted through it.

Such losses shrink the Roche lobe and offset the donor's tendency to move inside it. Maintaining contact on the equilibrium sequence requires ${\dot{R}}_L/R_L={\dot{M}}_2/M_2$, hence the binary mass-transfer contact equation becomesFor $q<4/3$ and $\dot{J}<0$, this drives ${\dot{M}}_2<0$ and maintains the modeled semi-detached binary.

Now take $\delta m>0$ for the mass expelled from the accretor in one classical nova. The white dwarf loses $\delta M_1=-\delta m$, while the donor has $\delta M_2=0$ and the total mass change is $\delta M=-\delta m$. Isotropic ejection in the accretor rest frame supplies no mean kick and carries the accretor's specific orbital angular momentum. Its orbital radius is $a_1=a{M}_2/M$, soThis is the angular-momentum budget for isotropic re-emission from a binary star, without an additional spin or frictional torque on the ejecta. The general logarithmic change of $J=M_1{M}_2\sqrt{Ga/M}$ isUsing $M/M_1=1+q$ gives the first-order nova changesThe donor mass and its assumed equilibrium radius are unchanged during the eruption. It is therefore left inside the enlarged Roche lobe, and Roche-lobe overflow ceases in this sharp-contact model. An extended atmosphere, irradiation-induced donor expansion or additional ejecta torques would require a different treatment.

The stipulated circular orbit is an important idealization. Isotropic loss that is slow compared with the orbital period can preserve circularity. If it is truly impulsive compared with that period, an initially circular point-mass orbit instead acquires eccentricity. With unchanged velocity and position, its new specific orbital energy is ${\epsilon }^{\prime }=GM/(2a)-G(M-\delta m)/a=-G(M-2\delta m)/(2a)$ and its squared specific angular momentum stays $h^2=GMa$. The Kepler orbit energy and eccentricity relations then give $a^{\prime }=a(M-\delta m)/(M-2\delta m)$ and $e^{\prime }=\delta m/(M-\delta m)$. The latter still has the same first-order separation increase. The question assumes circularity; its first-order circular budget is consistent with adiabatic orbital mass loss or with subsequent circularization that changes the separation only at second order. Work symbolically with $\delta m$, since the printed numerical layer mass has no explicit unit.

Finally put ${\mathcal{L}}_J=-\dot{J}>0$, constant over both intervals as specified, and let $J,M,q$ denote their leading-order values over one small cycle. While detached, masses and donor radius are fixed, so ${\dot{R}}_L/R_L=2\dot{J}/J=-2{\mathcal{L}}_J/J$. The fractional gap $4\delta m/(3M)$ closes afterOnce contact resumes, the contact equation gives $-{\dot{M}}_2=3M_2{\mathcal{L}}_J/[J(4-3q)]$. Accumulation of the next layer takesDividing yields the detached-to-semidetached time ratioBoth times are first-order in $\delta m/M$. Corrections to either time are second-order in that small mass fraction, so the displayed ratio has relative corrections of order $\delta m/M$. The positive denominator, negligible eruption duration relative to the cycle, fixed donor radius during detachment and common external loss rate are all essential to this idealized nova-induced binary detachment cycle.

Write $S(R)$ for the genuine surface density of a disk, and $\rho (R,z)=S(R)\delta (z)$ for its three-dimensional mass density. The delta function in the printed surface density formula belongs to $\rho$, not to $S$. Introduce a reference length $R_{\ast }$ to make logarithms dimensionless. The intended model is the infinite, self-gravitating, scale-free Mestel disk, with no extra source or imposed gravitational field.

The circular speed satisfies $v_c^2=R{\partial }_R\Phi (R,0)$. Thus the flat galaxy rotation curve givesA reflection-symmetric harmonic function with this midplane boundary value isVerify this continuation using the Poisson equation for Newtonian gravity. For $z>0$, set $r=\sqrt{R^2+z^2}$; thenHence ${\nabla }^2\Phi =0$ off the astrophysical disk, and reflection gives the lower-half-space solution. Across the astrophysical disk the derivative jumps by ${\Phi }_z(0^{+})-{\Phi }_z(0^{-})=2v_0^2/R$. The distributional Poisson equation for Newtonian gravity therefore givesAt the origin the enclosed disk mass tends to zero linearly with radius, so there is no additional central point mass. This verifies the Mestel disk potential-density pair. The astrophysical disk has infinite total mass and a logarithmic Newtonian gravitational potential, so it is not an isolated finite-mass model with Newtonian gravitational potential zero at infinity. The usual scale-free boundary condition is important: the midplane rotation curve alone would also permit an added term $a|z|$, representing an extra uniform sheet without changing the radial circular force. That contribution is excluded in the intended Mestel disk model.