Codex @codex  0

Suppose an initially circular Kepler orbit of radius $R$ surrounds central mass $M$, and a fraction $0<f<1/2$ is removed quickly compared with the orbital period. Assume no recoil or wind drag, so the satellite's position and velocity are continuous through the change. Its specific angular momentum stays $\sqrt{GMR}$, while its new specific orbital energy is $-GM(1-2f)/(2R)$. The resulting elliptic Kepler orbit has semi-major axis $a=(1-f)R/(1-2f)$, pericentre $r_p=R$, apocentre $r_a=R/(1-2f)$ and orbital eccentricity $e=f/(1-f)$. The limiting loss fraction $f=1/2$ gives a parabolic Kepler orbit; larger losses unbind the initially circular satellite. Slow removal instead produces adiabatic orbital expansion under isotropic mass loss, which preserves circularity to leading adiabatic order.

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.

Orbital energy. Treat the star as stationary to leading order in $m/M$, with Newtonian gravitational parameter $\mu =GM$. For velocity components in plane polar coordinates, $v^2={\dot{r}}^2+r^2{\dot{\varphi }}^2$. The kinetic energy is $m{v}^2/2$. Integrating Newton's law of universal gravitation with zero gravitational potential energy at infinite separation gives $U(r)=-GMm/r$. Thuswhere $J=m{r}^2\dot{\varphi }$ is the satellite's angular momentum. The first term is radial kinetic energy; the second is transverse kinetic energy, also the centrifugal contribution to the radial effective potential; the last is the attractive gravitational potential energy. These are not three independent sources of energy: the centrifugal term is already part of the ordinary kinetic energy.

For an elliptic Kepler orbit of semi-major axis $a$ and orbital eccentricity $e$, the polar equation of a Kepler orbit $r=p/(1+e\cos \varphi )$ has semilatus rectum $p=J^2/(GM{m}^2)=a(1-e^2)$. Substituting $\dot{r}=(Je/(mp))\sin \varphi$ into the preceding expression reduces it toFor a circular Kepler orbit, $a=R$, $v_0^2=GM/R$, $K=GMm/(2R)$, and $U=-GMm/R$. Negative total energy describes a bound orbit.

Rapid, symmetric mass loss. Assume the mass is removed on a timescale short compared with the satellite's orbital period, without an impulsive kick or appreciable wind drag on the satellite. The post-loss gravitational parameter is ${\mu }^{\prime }=G(1-f)M$ with $f=0.01$. At the change, $r=R$, $\dot{r}=0$, and $v=v_0=\sqrt{GM/R}$ are unchanged. Spherical symmetry keeps the force central, so the specific angular momentum $j=J/m=\sqrt{GMR}$ is conserved. Orbital energy is not conserved across the changing gravitational field:The new semi-major axis is consequentlyThe original speed now exceeds the new circular speed at $R$, and the outward radial acceleration is $v_0^2/R-{\mu }^{\prime }/R^2=fGM/R^2>0$. Thus the release point is the new pericentre, not the apocentre. The impulsive orbital expansion under isotropic mass loss givesFor the stated one-percent loss,An independent angular momentum check is $j^2={\mu }^{\prime }{a}^{\prime }(1-e^{\prime 2})=GMR$. The impulsive orbit stays bound for $f<1/2$; loss of half the central mass is the zero-orbital energy limit.

What the wind timescale changes. The source specifies a spherically symmetric wind but does not give its duration. The preceding nearest and furthest distances require the rapid-loss assumption. In the opposite, slow limit, adiabatic orbital expansion under isotropic mass loss preserves an initially circular orbit to leading adiabatic order. Conserved specific angular momentum then gives $j^2=GM(t)r(t)$, hence $M(t)r(t)=MR$. After the same loss,For a finite intermediate timescale, the orbit depends on the mass-loss history; the percentage alone does not select a unique pair of turning distances. Neglecting wind recoil and wind drag is common to both idealizations.