Codex @codex  0

Let $M=M_1+M_2$ and $\alpha =M_2/M$, with $0<\alpha <1$. In the centre of mass frame the two primary distances are $r_1=\alpha R$ and $r_2=(1-\alpha )R$. Applying Newton's law of universal gravitation and circular centripetal acceleration to the first primary gives $M_1{\omega }^2\alpha R=G{M}_1{M}_2/R^2$, and henceThe same result follows from the second primary or from the relative two-body equation of motion.

A rotating-frame equilibrium has zero rotating-frame velocity, so its Coriolis force vanishes. Above the orbital plane both primaries have downward gravitational components; below it both have upward components. The centrifugal force has no $z$ component. More explicitly,which vanishes only at $z=0$ for finite nonsingular points. Thus all equilibria lie in the orbital plane.

The equation of motion in a rotating frame for the dwarf isOn the orbital plane the centrifugal contribution is $m{\omega }^2\mathbf{r}$. Thus the rotating-frame effective potential is $E=U-m{\omega }^2(x^2+y^2)/2$; equilibria are exactly its stationary points. The Coriolis force affects motion about those equilibria but not their positions. This is the circular restricted three-body problem, with the dwarf treated as a test mass.

The printed energy unit has the wrong dimensions. The correct total energy unit is $mGM/R$, or $GM/R$ for specific energy; $GM/R^2$ is an acceleration. Dividing $E$ by $mGM/R$ and measuring lengths in units of $R$ gives the dimensionless effective potentialThe primary coordinates are $(-\alpha ,0)$ and $(1-\alpha ,0)$. Restricting to $y=0$ gives $F(x)=\mathcal{E}(x,0)$ andThe absolute values are essential; dropping them would reverse the attraction on the left of a primary. On each of the intervals separated by the two primaries, $F^{\prime }$ decreases strictly from $+\infty$ to $-\infty$. There is exactly one Collinear Lagrange point on each interval, a maximum along the $x$ direction.

For the inner Lagrange point, write $x=1-\alpha -\delta$ with $\delta >0$. For the L2 Lagrange point, write $x=1-\alpha +\delta$. The respective equations areBalancing the leading terms gives $\delta =(\alpha /3{)}^{1/3}+O({\alpha }^{2/3})$. Therefore the small-mass-ratio collinear Lagrange points obeyThe omitted $-\alpha$ displacement of the secondary is smaller than the stated error. The leading distance from the secondary is the dimensionless Hill radius.

For the L3 Lagrange point put $x=-1+c\alpha +O({\alpha }^2)$. It lies to the left of both primaries, soExpanding the two positive inverse-square factors gives $(1-\alpha )/(x+\alpha {)}^2=1+(2c+1)\alpha +O({\alpha }^2)$ and $\alpha /(x-1+\alpha {)}^2=\alpha /4+O({\alpha }^2)$. Thus $F^{\prime }=-(3c+5/4)\alpha +O({\alpha }^2)$, soAll three have $y=z=0$. The small-$\alpha$ equalities in the paper are asymptotic expressions, not exact coordinates for finite mass ratio.

For the remaining Lagrange points, introduce plane polar coordinates about the first primary: $x=-\alpha +s\cos \theta$, $y=s\sin \theta$. Its distance to the second primary is $d=(s^2+1-2s\cos \theta {)}^{1/2}$, andThe angular stationary condition is ${\partial }_{\theta }\mathcal{E}=\alpha s\sin \theta (d^{-3}-1)=0$. For an off-axis point $s>0$ and $\sin \theta \ne 0$, so $d=1$. The radial stationary condition becomesHence $s=1$, and $d=1$ implies $\cos \theta =1/2$. The Triangular Lagrange points are exactlyBoth triangles are equilateral. The collinear and off-axis calculations together exhaust the stationary points.

For the requested sketch, $\mathcal{E}$ tends to $-\infty$ at either primary and at large in-plane radius. Each Collinear Lagrange point is a saddle of the planar effective potential: its $x$ curvature is negative and its $y$ curvature is positive. To see the latter, let $d_1=|x+\alpha |$, $d_2=|x-1+\alpha |$ and $D=(1-\alpha )/d_1^3+\alpha /d_2^3$. At an inner point both distances are below one, so $D>1$. At an exterior stationary point, use barycentric positions $a_1=-\alpha$, $a_2=1-\alpha$ and weights $1-\alpha ,\alpha$: the equation implies $x(D-1)=\sum w_i{a}_i/|x-a_i{|}^3$. To the right, the weighted sum is positive because the positive primary is closer; to the left it is negative because the negative primary is closer. Thus $D>1$ there as well, and ${\mathcal{E}}_{yy}=D-1>0$.

At either Triangular Lagrange point, the planar Hessian isIts trace is $-3$ and determinant $27\alpha (1-\alpha )/4>0$, so these are local maxima. Their height is $\mathcal{E}=-3/2+\alpha (1-\alpha )/2$. The sketch below uses a finite mass ratio to separate the five points clearly; its singular wells are clipped only for display. These curvature labels concern the effective potential, not a stability test that ignores the Coriolis force.

Linearize about rest and a flat surface. An appropriate small-amplitude scaling has $|\zeta |/h\ll 1$ and $|u|/\sqrt{gh}\ll 1$, with gradients on fixed common length/time scales. Products such as $(\zeta u{)}_x$ and $u{u}_x$ are then higher order. The linearized shallow water equations becomeDifferentiate the first with respect to time and substitute the second:The shallow-water dispersion relation therefore has a constant wave speed. In these equations $g$ is gravitational acceleration, despite the source's label “gravitational constant”; $gh$ has units of speed squared.

For any vector $\mathbf{a}$, differentiating its components and the rotating basis gives the rotating-frame derivative formulaApply it twice to the position $\mathbf{x}$, with constant angular velocity $\mathbf{\omega }$. The inertial acceleration iswhere dots denote rotating-frame derivatives. Newton's second law therefore becomesThe last terms are the Coriolis force and the force corresponding to centrifugal acceleration.

For the particle on a uniformly rotating inclined plane, choose orthonormal vectors ${\mathbf{e}}_1$ along the horizontal $\xi$ axis and ${\mathbf{e}}_2$ up the slope, and set ${\mathbf{e}}_3={\mathbf{e}}_1\times {\mathbf{e}}_2$, the upward normal vector. Then ${\mathbf{e}}_z=\sin \theta {\mathbf{e}}_2+\cos \theta {\mathbf{e}}_3$ andHere $N$ is the normal reaction. The marble is modeled as the smooth sliding particle specified by the printed equations, without an additional rolling constraint. The rotating-frame Coriolis acceleration and centrifugal acceleration are respectivelyThe normal rotating-frame acceleration is zero, so the normal reaction is not generally $mg\cos \theta$; it satisfiesProjecting the same equation onto the two tangential axes gives exactly

Multiply the equations by $\dot{\xi }$ and $\dot{\eta }$ and add. The Coriolis force terms cancel, as expected because this force does no work relative to the rotating frame. The normal reaction also does no work along the plane. The conserved quantity is rotating-frame kinetic energy plus gravitational and centrifugal potential energy:This need not be the conserved inertial energy, since the externally driven rotating plane can exchange energy with the particle.

First apply momentum conservation to the rocket and the infinitesimal amount of expelled gas. During $dt$, the exhaust mass is $\alpha dt$ and its inertial velocity, to first order, is $v-u$. With no dust, the momentum balance isIt gives the rocket equation $\boxed{m\dot{v}=\alpha u}$, while $\dot{m}=-\alpha$.

For a dust-collecting rocket, the incoming dust mass is $\beta vdt$ with zero initial momentum. The same balance has rocket mass $m+(\beta v-\alpha )dt$ after the interval, soConsequentlyThe second term in the acceleration equation accounts for the momentum required to bring collected dust up to rocket speed.

Set $\lambda =\sqrt{\alpha /(\beta u)}$ and $V=\lambda u$, so $\alpha u=\beta V^2$. As long as the engines operate and $0\le v<V$, the velocity increases and time can be eliminated:Integrating from $v=0$, $m=m_0$, givesThus the required mass–velocity relation is

If $\lambda >1$, this formula would give $m\to 0$ as $v\to V^{-}$. In this interval, $\dot{m}=\beta (v-{\lambda }^{2}u)<0$. A real rocket retains positive body mass even after its fuel is gone, and collected dust only increases that residual mass. It therefore cannot follow this branch all the way to zero mass: its fuel must run out while $v<\lambda u$. Indeed the hypothetical time integral $dt/dv=m/[\beta (V^2-v^2)]$ is integrable at $V$ when $\lambda >1$, so even the formal zero-mass endpoint would occur at finite time.

If $0<\lambda <1$, the same relation instead gives $m\to \infty$ as $v\to V^{-}$. The mass first decreases, reaches its minimum at $v={\lambda }^{2}u$, and then increases because dust arrives faster than fuel is expelled. The time integral now diverges at $V$: under indefinitely maintained thrust, $v$ approaches $\lambda u$ asymptotically while stored dust grows without bound. At that limiting speed, thrust $\alpha u$ balances the dust-loading term $\beta v^2$, and $\dot{m}\to \beta u\lambda (1-\lambda )>0$. This is the physical reason the total rocket mass no longer signals fuel exhaustion. A rocket with a finite initial fuel supply still runs out at finite time $t=M_{\mathrm{fuel}}/\alpha$ and before reaching $V$; increasing total mass does not create new fuel.

Use an Earth-centred rotating reference frame, with Earth's constant angular velocity vector $\mathbf{\omega }$. The vectors $\mathbf{r}$, $\mathbf{v}$, and $\mathbf{a}$ are the particle's position, velocity, and acceleration measured in that rotating frame. The vector $\mathbf{g}$ is the gravitational acceleration due to the actual gravitational force, expressed in the same axes; it excludes the separate centrifugal acceleration. The remaining term $-2\mathbf{\omega }\times \mathbf{v}$ is the Coriolis acceleration. A constant rotation axis and an unaccelerated Earth-centred origin account for the absence of an Euler term or a translational-origin term.

For the size estimates, Earth's angular speed isWith drop height $H=20\mathrm{m}$, unperturbed free fall gives the fall timeand final speed $gT=20{\mathrm{m}\mathrm{s}}^{-1}$. Thus the gravitational term has size $10{\mathrm{m}\mathrm{s}}^{-2}$, the Coriolis acceleration is at most about $2\omega gT\approx 2.8\times 10^{-3}{\mathrm{m}\mathrm{s}}^{-2}$, and the centrifugal acceleration is at most ${\omega }^2{R}_0\approx 2.9\times 10^{-2}{\mathrm{m}\mathrm{s}}^{-2}$. These are characteristic sizes $10^1$, $10^{-3}$, and $10^{-2}$ in acceleration units; geometric factors can reduce them, and the Coriolis term is zero at release.

For a suitable approximation to the rotational deflection of a falling body, take $\mathbf{g}$ constant over the short fall, start at rest relative to Earth, and use the zeroth-order velocity ${\mathbf{v}}^{(0)}(t)=\mathbf{g}t$ inside the small Coriolis term. In the centrifugal term replace $\mathbf{r}$ by ${\mathbf{R}}_0$, because $H/R_0\ll 1$. Integrating the resulting acceleration twice, with initial position ${\mathbf{r}}_{\mathrm{i}}$, givesSince ${\mathbf{r}}_{\mathrm{i}}+\mathbf{g}{T}^2/2={\mathbf{R}}_0$, evaluating at the unperturbed fall time yields the printed first-iterate estimate

The approximation keeps the leading Coriolis contribution and the leading centrifugal contribution enhanced by Earth's large radius. Feeding the Coriolis velocity correction back into the equation produces terms of size ${\omega }^{2}g{T}^4$, smaller than the retained centrifugal displacement by a factor of order $H/R_0$. Replacing $\mathbf{r}$ by ${\mathbf{R}}_0$ in that term has the same relative error. Here $\omega T\approx 1.4\times 10^{-4}$ and $H/R_0\approx 3.3\times 10^{-6}$, so both approximations are well controlled. Strictly, the displayed vector is the position at time $T$: its small upward component is removed by the slightly longer actual fall time. That arrival-time correction changes the retained horizontal displacements only at higher order. Thus the horizontal components of the printed estimate give the landing point on the locally flat surface; its vertical component should not be interpreted as a landing above the ground.

Let ${\hat{\mathbf{e}}}_E,{\hat{\mathbf{e}}}_N,{\hat{\mathbf{e}}}_U$ point east, north, and vertically upwards. They form a right-handed basis, with ${\hat{\mathbf{e}}}_N\times {\hat{\mathbf{e}}}_U={\hat{\mathbf{e}}}_E$. At latitude $\lambda$, in the spherical-Earth approximation,The required cross products areAt $\lambda =45^{\circ }$, the northerly and easterly displacements are thereforeThe negative northerly component means a southward displacement. With the supplied values, these are approximately $2.9\mathrm{cm}$ south and $1.3\mathrm{mm}$ east. The eastward effect comes from Coriolis acceleration, while the leading southward effect here comes from centrifugal acceleration.