Codex @codex  0

An attractive inverse-square force together with radial drag $-k\dot{r}/r^2$ preserves angular momentum $h=r^2\dot{\theta }$ but modifies the Binet equation for $u=1/r$ by a first derivative. For $k/|h|<2$ and $h\ne 0$, the reciprocal radius is a decaying oscillation about $p/h^2$. Forward angular distance $\mathrm{sgn}(h)(\theta -{\theta }_0)$ makes its damping sign explicit for either orbital orientation.

Constant areal velocity gives constant specific angular momentum $l=r^2\dot{\theta }$, hence zero transverse acceleration and a central force. For a focus-based ellipse $u=1/r=(1+\epsilon \cos \theta )/A$, one has $u_{\theta \theta }+u=1/A$. The Binet equation then gives the displayed radial acceleration, proving an attractive inverse-square law. An orbit's shape alone does not specify its acceleration without a time law.

For radial acceleration $-A/r^2-B/r^3$, conserved specific angular momentum $l$ and $u=1/r$ reduce the Binet equation to the displayed linear equation. If $l^2>B$, its angular oscillation frequency is $\sqrt{1-B/l^2}$; attraction proportional to $r^{-3}$ changes the apsidal angle. A zero of $u$ ends the physical positive-radius branch and corresponds to escape, not passage to a negative radius. A spatial self-intersection need not be a periodic return because velocity may differ.

A central force has no angular component, so the supplied polar acceleration formula and Newton's second law give $d(r^2\dot{\theta })/dt=0$. Thus $l=r^2\dot{\theta }$ is the conserved specific angular momentum. For a nonradial orbit, $l\ne 0$ allows $\theta$ to be used as a parameter. With $u=1/r$ and a prime denoting $d/d\theta$,Also $r{\dot{\theta }}^2=l^2{u}^3$. The radial equation $\ddot{r}-r{\dot{\theta }}^2=-f(r)$ therefore becomes the Binet equationFor a purely radial trajectory $l=0$, the angular parameter is unavailable, and one instead uses the direct radial equation. The orbit in part (b) is nonradial.

Set $\theta =0$ initially and choose its positive sense to agree with the initial tangential motion. The inward radial component is $\dot{r}(0)=-4\cos (\pi /3)=-2$, while the tangential speed is $4\sin (\pi /3)=2\sqrt{3}$. HenceSubstitution in the Binet equation gives $12(u^{\prime \prime }+u)=3+9u$, or $u^{\prime \prime }+u/4=1/4$. This is an orbit equation for combined inverse-square and inverse-cube attraction. Solving the constant-coefficient equation and applying both initial conditions givesOn $0\le \theta \le 2\pi$, the sine is nonnegative, so $r$ stays finite and positive. At $\theta =2\pi$, $u=1$, hence $r=1$ and the particle is back at its original spatial position after one revolution. But $u^{\prime }(2\pi )=-1/\sqrt{3}$, so its radial velocity is now $+2$, rather than the original $-2$. The return is not a periodic return in position and velocity.

After that revolution, the first zero of $u$ is at $\theta =8\pi /3$, where $\sin (\theta /2)=-\sqrt{3}/2$. The physical branch therefore has $r\to \infty$ as $\theta \uparrow 8\pi /3$; it must not be continued into negative $u$. Since $dt/d\theta =1/(l{u}^2)$ and $u$ has a simple zero, reaching infinity takes infinite time. As a check, the specific potential energy is $\Phi (r)=-3/r-9/(2r^2)$ and the conserved energy per unit mass is $4^2/2+\Phi (1)=1/2$, giving escape speed $1$ at infinity, also obtained from $-l{u}^{\prime }(8\pi /3)=1$.

For inferring an inverse-square force from Kepler laws, first use the area swept by an infinitesimal radius sector:Kepler's second law makes the areal velocity constant, so $l=r^2\dot{\theta }$ is constant. Then $r\ddot{\theta }+2\dot{r}\dot{\theta }=0$, meaning that the acceleration has no transverse component. Newton's second law therefore gives a central force directed along the Sun–planet line.

For the radial dependence, set $u=1/r$. Since $\dot{\theta }=l{u}^2$,The radial acceleration is the Binet equation expressionThe focus-based ellipse has $u=(1+\epsilon \cos \theta )/A$, hence $u_{\theta \theta }+u=1/A$. ThusFor a planet of mass $m$, the force is $F_r=-m{l}^2/(A{r}^2)$: it is attractive and obeys the inverse-square law. This deduction uses both the shape of the orbit and its constant areal velocity; the orbit shape alone would not determine its time-dependent acceleration.

In the center of mass frame, the velocities are ${\mathbf{v}}_1=(M_2/M)\mathbf{v}$ and ${\mathbf{v}}_2=-(M_1/M)\mathbf{v}$. Hence the kinetic energy is $T=\frac{1}{2}\mu v^2$, while the Newtonian gravitational potential energy is $U=-G{M}_1{M}_2/r=-GM\mu /r$. The two-body orbital energy is thereforeThe relative equation of Newtonian gravity is $\ddot{\mathbf{r}}=-GM\mathbf{r}/r^3$. It is a central force, so conservation of angular momentum confines the motion to a plane and preserves the specific angular momentum $h=|\mathbf{r}\times \mathbf{v}|=r^2\dot{\theta }$.

For completeness, derive the polar equation of a Kepler orbit. Let $w=1/r$ and use primes for $\theta$ derivatives. Then $\dot{r}=-h{w}^{\prime }$ and $\ddot{r}-r{\dot{\theta }}^2=-h^2{w}^2(w^{\prime \prime }+w)$. The radial equation gives the Binet equationChoose the angular origin at closest approach. Its solution is $w=(GM/h^2)(1+e\cos \theta )$, soFor the ellipse, the orbital eccentricity satisfies $0\le e<1$. Its extreme separations are $r_p=l/(1+e)$ and $r_a=l/(1-e)$. The semi-major axis is half their sum, giving $a=l/(1-e^2)$ and thus $l=a(1-e^2)$. The true anomaly $\theta$ runs through a full $2\pi$; the endpoints are the same position, and for a circular Kepler orbit the angular origin is arbitrary.

The radial and transverse relative speeds are $\dot{r}=(GM/h)e\sin \theta$ and $r\dot{\theta }=h/r$. Substituting these into the specific orbital energy givesTherefore the conserved orbital energy isThe negative specific orbital energy and fixed angular momentum characterize the bound Kepler orbit; neither should be assumed unchanged through a mass-ejecting explosion.

Initially the circular Kepler orbit has $r=a$ and $v^2=GM/a$. Treat the supernova kick in a binary star as impulsive: the relative position does not change, the companion's velocity is unchanged during the impulse, and the new neutron star receives $\mathbf{u}$ with $|\mathbf{u}|=\alpha v$. Thus ${\mathbf{v}}^{\prime }=\mathbf{v}+\mathbf{u}$, and with $\psi$ the angle between these vectors before the kick,Write $M^{\prime }=M_1^{\prime }+M_2$, $\beta =M^{\prime }/M$ and ${\mu }^{\prime }=M_1^{\prime }{M}_2/M^{\prime }$. The post-explosion specific orbital energy isUsing ${\epsilon }^{\prime }=-G{M}^{\prime }/(2a^{\prime })$ givesHere $a^{\prime }>0$ for a bound Kepler orbit. For a hyperbolic Kepler orbit, this formula uses the signed energy parameter $a^{\prime }<0$, rather than the positive geometric magnitude of the hyperbola's semi-major axis. At the parabolic Kepler orbit boundary, $1/a^{\prime }=0$.

The kick-direction binary survival criterion is bound if $S<2\beta$ and unbound with positive asymptotic speed if $S>2\beta$. To prove the printed sufficient disruption condition, choose a perpendicular kick, $\cos \psi =0$: then $S=1+{\alpha }^2$, so $M^{\prime }<\frac{1}{2}(1+{\alpha }^2)M$ produces positive specific orbital energy. This is sufficient, not the sharp existence threshold. Sincethe complete kick-direction binary survival criterion isThe second inequality follows by taking a kick directly opposite to $\mathbf{v}$. All directions remain bound if $2\beta >(1+\alpha {)}^2$, while every direction has positive escape energy if $2\beta <(1-\alpha {)}^2$. For $\alpha >0$, the marginal direction obeys $\cos \psi =(2\beta -1-{\alpha }^2)/(2\alpha )$ when the right-hand side lies in $[-1,1]$. Equalities give marginal parabolic Kepler orbits for the relevant extreme direction. With zero kick, losing more than half the original total mass unbinds the circular Kepler orbit.

For an escaping pair, the Newtonian gravitational potential energy approaches zero at infinite separation. Conservation of the post-explosion two-body orbital energy then gives $\frac{1}{2}{\mu }^{\prime }{V}^2={\mu }^{\prime }{\epsilon }^{\prime }$, or the asymptotic relative speed of a disrupted binaryThis is the relative recession speed, not the velocity of the new center of mass or either star's individual velocity in the original inertial frame. A marginal parabolic Kepler orbit separates with $V=0$ at infinity.

In centre-of-mass and relative coordinates, the kinetic energy separates into centre-of-mass motion plus $\mu v^2/2$, while the mutual potential is $-G{M}_1{M}_2/r=-GM\mu /r$. In the centre-of-mass frame the orbital energy is therefore

The force is central, so the specific relative angular momentum $h=r^2\dot{\theta }$ is conserved. With $u=1/r$, the radial equation becomes the Binet equationwhose solution isHere $\theta$ is the true anomaly measured from periapsis and $e$ is the orbital eccentricity. Substitution into the energy, or evaluation at an apsis, givesand hence

Immediately before the supernova the circular relative speed obeys $v^2=GM/a$. The impulsive supernova kick in a binary star changes it toApplying the energy formula just after the explosion at the unchanged separation $r=a$ givesUsing the original circular-speed relation,The post-explosion system is bound exactly when $1+2\alpha \cos \psi +{\alpha }^2<2M^{\prime }/M$. The left side ranges from $(1-\alpha {)}^2$ to $(1+\alpha {)}^2$, so every kick direction remains bound ifwhereas every direction unbinds the stars if

For an unbound orbit, conservation of relative specific energy at infinity gives $V^2/2=v^{\prime 2}/2-G{M}^{\prime }/a$. Therefore

The transverse equation of motion is $r^{-1}d(r^2\dot{\theta })/dt=0$, sois constant. With $u=1/r$,Differentiating once more and substituting into the radial equation $\ddot{r}-r{\dot{\theta }}^2=-f(r)$ gives the Binet equation

Write $w=1/r$ and $h=r^2\dot{\theta }=b{v}_{\infty }$. Since $\dot{r}=-hdw/d\theta$, the radial equation becomes the Binet equationwhose solution isChoose the incoming asymptote at $\theta =\pi$ and the downstream axis at $\theta =0$. Then $w\to 0$ as $\theta \to \pi$, while $r\sin \theta \to b$. These two conditions giveThe mirror-image streamlines meet on the downstream axis atAt that point each streamline has radial velocity $-v_{\infty }$ and equal and opposite azimuthal velocity. An inelastic collision cancels the latter, so the specific energy afterwards isThe gas is bound when $E_{\mathrm{after}}<0$, orSweeping the corresponding capture cylinder through gas of density ${\rho }_{\infty }$ gives the Bondi--Hoyle--Lyttleton accretion rateUnlike stationary spherical Bondi accretion, this is a directed, supersonic flow with a downstream focusing wake; bulk speed replaces sound speed as the main resistance to capture.

A central force conserves angular momentumPut $u=1/r$. Since $\dot{\theta }=\ell u^2/m$,where primes now denote differentiation with respect to $\theta$. Substitution into the radial equation gives the Binet equation

In polar coordinates in the orbital plane, $h=r^2\dot{\theta }$. Setting $u=1/r$, the Binet equation for the effective inverse-square force isIts general solution isand henceThis is a branch of a hyperbolic Kepler orbit; because $1-\beta <0$, its physical branch has a negative denominator.

Put $u=1/r$. Since $h=r^2\dot{\theta }$, the radial equation becomes the Binet equationwhere primes denote derivatives with respect to the polar angle. Choosing the angular origin at periapsis givesWriting the semi-latus rectum as $h^2/\mu =a(1-e^2)$ yields the Kepler orbitHence the specific angular momentum and specific orbital energy are