Codex @codex  0

Take an illustrative elliptic Kepler orbit with semi-major axis $a=1$, orbital eccentricity $e=0.6$ and gravitational parameter $\mu =1$, in consistent dimensionless units. Its semi-minor axis is $b=a\sqrt{1-e^2}=0.8$, mean motion is $n=\sqrt{\mu /a^3}=1$, and specific angular momentum magnitude is $h=\sqrt{\mu a(1-e^2)}=0.8$. About the attracting focus,The two sectors below each span $\Delta t=0.5$, giving $\Delta A=0.2$. Near periapsis the short radius sweeps a wide angle; near apoapsis the long radius sweeps a narrow angle. Constant areal velocity therefore does not require constant angular velocity or constant speed.

To compute the sectors without mistaking equal angles for equal times, parametrize the ellipse by eccentric anomaly $E$:This follows by integrating $\frac{1}{2}(xdy-ydx)$. Kepler's equation gives $E-e\sin E=M=n(t-t_{\mathrm{p}})$, with mean anomaly $M$ measured from periapsis time $t_{\mathrm{p}}$. Consequently $\Delta A=abn\Delta t/2=h\Delta t/2$. The displayed sectors use mean anomaly intervals $[-0.25,0.25]$ and $[\pi -0.25,\pi +0.25]$, rather than equal eccentric anomaly increments.

For an elliptic Kepler orbit, Kepler's equation relates mean anomaly $M$ and eccentric anomaly $E$ by

The relative displacement of a Kepler orbit sweeps equal areas about the attracting focus of a conic in equal times. A short interval near periapsis covers a wide angle at small separation; the same interval near apoapsis covers a narrow angle at large separation. The swept areas are equal even though the travelled distances are different. In particular, the sectors have curved ellipse boundaries: triangles formed by joining their endpoints only approximate the swept areas for short intervals.

Let $h=|\mathbf{r}\times \dot{\mathbf{r}}|$ be the relative specific angular momentum. The areal velocity isA central force exerts zero torque about its centre, so conservation of angular momentum makes this rate constant. This explanation also applies to other central forces; the ellipse shape specifically comes from Newtonian gravity.

In the $a=1$, $e=0.6$, $\mu =1$ example, specific angular momentum is $h=\sqrt{\mu a(1-e^2)}=\sqrt{0.64}=0.8$, giving areal velocity $0.4$. At periapsis and apoapsis, the velocity is perpendicular to the relative displacement, so $h=rv$. The endpoint speeds are consequentlyThe fourfold speed difference is consistent with the fourfold separation difference. Away from these endpoints, $h/r$ gives only the transverse component of velocity; the full speed follows from the vis-viva equation.

To actually place equal-time markers, use Kepler's equation. Write mean anomaly as $\mathcal{M}=nt$, with $t=0$ at periapsis and mean motion $n=\sqrt{\mu /a^3}=1$. The different symbol $\mathcal{M}$ avoids confusing this angle with the mass $M$. Solve for eccentric anomaly $E$ in radians, then compute the relative position:This timing procedure is explained in NASA educational archive, How Orbital Motion is Calculated.

For the first interval $0\le t\le \pi /6$, the final eccentric anomaly is $E\approx 1.041495$. For the equal interval $\pi \le t\le 7\pi /6$ starting at apoapsis, it ends at $E\approx 3.471065$. Both swept areas are $\Delta A=(h/2)\Delta t=0.4\pi /6=\pi /15\approx 0.209440$, one twelfth of the entire ellipse's area $0.8\pi$. The diagram computes its markers from Kepler's equation, making visible how uniform time steps become uneven spatial steps.

For a quantity on a periodic Kepler orbit, its orbit average is $⟨x⟩=P^{-1}{\int }_0^{P}x(t)dt$. Equal increments of true anomaly are not equal increments of time. Eccentric anomaly gives the weight $dt=(1-e\cos E)dE/n$, which yields $⟨\cos f⟩=-e$. Slow orbital perturbations can be averaged using elements held fixed over a single orbital period.

Write ${\mu }_g=\gamma m$ for the gravitational parameter of the relative Kepler orbit, and use specific orbital energy and momenta per unit reduced mass. The Hamiltonian isIf physical momenta rather than these specific momenta are used, multiply the actions by the reduced mass. This only introduces a constant factor into the phase space measure and leaves the eccentricity fractions unchanged. The formulas in the question adopt the specific-action convention.

Here is an explicit local canonical transformation. A separated complete integral of the Hamilton-Jacobi equation isChoose branches along an orbit and take $p_j=\partial S/\partial q_j$, $l=\partial S/\partial L$, $g=\partial S/\partial G$, $h=\partial S/\partial H$. Thenwhose exterior derivative proves preservation of the canonical symplectic form. Thus these variables are canonical, and the absolute Jacobian determinant is one. The angular separation gives $G=|\mathbf{r}\times \mathbf{v}|$ and $H=G\cos i$, including the sign of the angular momentum along the reference axis.

For clarity, the radial action can be evaluated rather than assumed. Let $r=a(1-e\cos u)$, with $u$ the eccentric anomaly. Then $p_r=\sqrt{{\mu }_g/a}e\sin u/(1-e\cos u)$, andIn the integral, use $e^2{\sin }^{2}u/(1-e\cos u)=1+e\cos u-(1-e^2)/(1-e\cos u)$ and ${\int }_0^{2\pi }(1-e\cos u{)}^{-1}du=2\pi /\sqrt{1-e^2}$. The polar action is $J_{\theta }=G-|H|$, while the signed azimuthal action is $J_{\varphi }=H$. The separated angular coordinates have the usual Delaunay variables interpretation: $g$ is the argument of periapsis, and $h$ is the longitude of ascending node. In particular,with a suitable origin at periapsis, so $l$ is the mean anomaly. We have therefore recoverednot merely assigned names to the new coordinates.

For nondegenerate bound Kepler orbits, the ranges areThe boundaries describe circular, coplanar or radial degeneracies, where some angles cease to be defined; they have zero measure in the present calculation. The actions are canonical coordinates with nested bounds, rather than three independently unrestricted Cartesian coordinates. The square of $H$ alone would miss the distinction between prograde and retrograde inclinations.

By the canonical transformation, the population measure becomesAt fixed $L$, the condition $e^{\prime }<e$ is $L\sqrt{1-e^2}<G<L$. Integrating all three angles and then $H$ givesConsequently the Delaunay phase-volume proof of the thermal eccentricity distribution yieldsThe PDF's $e^2$ dependence must be interpreted as a cumulative distribution function, not an eccentricity probability density function or the value of the original six-dimensional phase-space distribution function. A finite value of $F(E)$ at every energy also does not suffice to make the total population finite: $F=1$ on all bound energies gives the divergent integral ${\int }_0^{\infty }{L}^{2}dL$. The fractions are well defined at each fixed energy, over a finite allowed energy range, or for an energy distribution with a finite total integral.

A Boltzmann distribution is therefore unnecessary. Any normalizable energy-only distribution has the same eccentricity marginal, although such a distribution is a collisionless steady state by Jeans theorem. The marginal alone does not even prove stationarity: for $0<|\eta |<1$,is a positive time-dependent solution of the Collisionless Boltzmann equation and has exactly the same eccentricity marginal after integration over $l$. An observed $e^2$ cumulative law establishes neither thermal relaxation nor the age of the system or Universe.

For a bound Kepler orbit, the semimajor axis is half the long diameter of the ellipse; if $r_p,r_{mathrma}$ are its apsidal distances, $a=(r_p+r_{mathrma})/2$. The orbital eccentricity is $e=(r_{mathrma}-r_p)/(r_{mathrma}+r_p)$, with $0\le e<1$ and semiminor axis $b=a\sqrt{1-e^2}$. The gravitating centre is a focus, displaced by $ae$ from the ellipse's geometric centre.

Choose a reference plane, a positive normal and a fixed reference direction in that plane. The ascending node is where the orbit crosses toward positive normal displacement. The longitude of ascending node $\Omega$ is the reference-plane angle from the fixed direction to that node. The argument of periapsis $\omega$ is the angle from the ascending node to periapsis, measured in the orbital plane in the direction of motion. The orbital inclination $i$ is the angle between the reference normal and the orbital angular-momentum normal, with $0\le i\le \pi$. The diagram separates the ellipse's shape from these orientation angles.

Let $\mu =GM$ and use actions per unit test-particle mass. In the prograde spherical-action chart underlying the stated Hamiltonian, $J_{\varphi }=L_z$ and $J_{\theta }=L-L_z$, so $J_{\theta }+J_{\varphi }=L$. The radial action then satisfies $J_r+L=\sqrt{\mu a}$, consistently with $H_K=-\mu /(2a)$. DefineIf $(w_r,w_{\theta },w_{\varphi })$ are the old angles, a type-two generating function isDifferentiation with respect to the old angles recovers $J_r,J_{\theta },J_{\varphi }$, and differentiation with respect to the new actions gives the canonical transformationIts triangular integer action matrix has determinant one, so it preserves the angle periodicities. In the new coordinates,These are the Delaunay variables:With the usual choices of angle origins, $w_1=\ell$ is mean anomaly, $w_2=\omega$ is argument of periapsis, and $w_3=\Omega$ is longitude of ascending node. The generator $\mathcal{G}=L$ rotates the orbit within its plane, whereas $\mathcal{H}=L_z$ rotates that plane about the reference normal; this explains the angular identification. The fast angle obeys Kepler's equation $\ell =\mathcal{E}-e\sin \mathcal{E}$, where $\mathcal{E}$ is eccentric anomaly.

For retrograde orbits, signed $J_{\varphi }=L_z$ gives $J_{\theta }=L-|L_z|$ and the spherical Hamiltonian uses $J_r+J_{\theta }+|J_{\varphi }|$. The stated plus-sign form is therefore a prograde chart, not a global signed-action formula. The Delaunay variables still use $\mathcal{H}=L_z$; on the retrograde branch the generating function is $(\mathcal{L}-\mathcal{G})w_r+(\mathcal{G}+\mathcal{H})w_{\theta }+\mathcal{H}{w}_{\varphi }$. At $e=0$, periapsis is undefined, and at $i=0$ or $\pi$ the node is undefined: these are angle-chart degeneracies, not physical singularities.

For the perturbed planet orbit, the radial and vertical frequencies mean small oscillations about a circular equatorial orbit of cylindrical radius $R$. Write $r=\sqrt{R^2+z^2}$ and retain the stated quadrupole potential. On the equator,At fixed axial specific angular momentum $L_z$, the radial effective potential is $\varphi (R,0)+L_z^2/(2R^2)$. Its second derivative at circular balance $L_z=R^2{\Omega }_c$ gives the radial epicyclic frequencyReflection symmetry decouples vertical motion to first order. Expanding the potential in $z$ givesso the vertical epicyclic frequency isThe circular orbit is radially stable only for $R^2>3{\alpha }^2$. For $\alpha \ne 0$ in that domain, $\nu >{\Omega }_c>\kappa$. Vertical node crossings repeat at frequency $\nu$ while azimuth advances at ${\Omega }_c$, so the line of nodes has nodal precession rate ${\Omega }_c-\nu <0$. Radial periapsis passages repeat at $\kappa$, giving apsidal precession rate ${\Omega }_c-\kappa >0$. To leading quadrupole order, with $n=\sqrt{\mu /R^3}$,This is oblate-quadrupole epicyclic precession. Here $\varpi$ denotes the nearly equatorial longitude of periapsis, not $\omega$ alone.

The precession signs above have a small-eccentricity, small-inclination domain. They are not a universal claim for every inclined Kepler ellipse. To see the distinction with the same Delaunay variables, average the perturbation over an unperturbed orbit. Since $z/r=\sin i\sin (\omega +f)$ and $dt=r^{2}df/\mathcal{G}$, the term linear in $e\cos f$ integrates to zero and the ${\sin }^2(\omega +f)$ average in the $r^{-3}$ integral is $1/2$. ThusDifferentiating at fixed other actions yields the first-order secular ratesTheir nearly equatorial sum reproduces the prograde $\dot{\varpi }$. For a polar orbit, however, the node is fixed and $\dot{\omega }<0$. This supplies a concrete counterexample to extending the requested forward periapsis precession beyond the intended epicyclic regime.

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.

The printed line-density statement is incorrect for a phase-mixed orbit. A steady line density on a Kepler orbit is inversely proportional to speed: for a cross-sectional-area current $J_{\sigma }$, the area per unit arc length is $d\sigma /ds=J_{\sigma }/v$. Equivalently, phase mixing gives $d\sigma =({\sigma }_{\mathrm{tot}}/P)dt$, uniform in mean anomaly, where $P=2\pi \sqrt{a^3/\mu }$ is the orbital period.

For an optically thin population of blackbodies in radiative equilibrium, a fragment absorbs $\sigma L_{\star }/(4\pi r^2)$ and reradiates the same luminosity. Thus the fractional luminosity of a phase-mixed eccentric wire isThe specific angular momentum gives $dt/r^2=d\theta /h$, so the integral is $2\pi /h$. Since $Ph=2\pi a^2\sqrt{1-e^2}$,This derives the intended result after explicitly correcting the density to $1/v$.

The error is consequential. If one instead imposes the literal density $d\sigma /ds\propto v$, its time weighting is $v^{2}dt$. Using $⟨v^2⟩=\mu /a$, $⟨r^{-2}⟩=[a^2\sqrt{1-e^2}{]}^{-1}$ and $⟨r^{-3}⟩=[a^3(1-e^2{)}^{3/2}{]}^{-1}$ givesAt $e=1/2$ this is $5/3$ times the printed result. The inverse-radius averages can also be obtained directly with eccentric anomaly $E$, $r=a(1-e\cos E)$ and $dt=(1-e\cos E)dE/n$.

Measure angles from periapsis in the direction of motion. The true anomaly $f$ is the angle at the stellar focus between the pericentre ray and the star-particle radius vector. Draw the auxiliary circle of radius $a$ around the ellipse's center; the point on that circle with the same projected major-axis coordinate as the particle defines the eccentric anomaly $E$ at the center.

The mean anomaly $M=n(t-\tau )$ is a uniformly advancing time phase, with $M=0$ at pericentre passage $t=\tau$. It is the angle swept by a fictitious uniformly moving circular phase marker, not generally the actual position angle. The mean motion is the angular frequencyThe second formula is Kepler's equation, which converts the uniform time phase into geometric position.

The sudden potential change leaves the star's position and velocity unchanged. Its new energy is thereforeThus $E_1>0$ precisely when $r<a$. Parametrize the orbit by eccentric anomaly $\eta$:A phase-mixed orbit is uniform in the mean anomaly ${\theta }_r$. For $e>0$, the condition $r<a$ is $-\pi /2<\eta <\pi /2$ modulo one period. The corresponding mean-anomaly interval has lengthHenceThe energy increase $G(M_0-M_1)/r$ is largest near periapsis, so stars are more readily unbound there than near apoapsis. A perfectly circular orbit is the measure-zero marginal case $E_1=0$ at every phase.

For each point $(\Delta {\varphi }_1,e)$, sample $\delta \in [-\Delta {\varphi }_1,\Delta {\varphi }_1]$ and the three conjunction branchesFor every branch, solve Kepler's equation for the eccentric anomaly, evaluate $r=a(1-e\cos E)$, and minimize the planet-planetesimal separation over $\delta$ and $k$. Mark the point as encounter-capable when this minimum is below a chosen $R_{\mathrm{enc}}$, naturally the Hill radius for strong scattering. Repeating this calculation on a grid traces the boundary in the $\Delta {\varphi }_1$--$e$ plane.

Direct integrations of the circular restricted three-body problem can then refine the geometric map by allowing the resonant argument, eccentricity, and conjunction kicks to evolve self-consistently. The integrations distinguish merely orbit-crossing initial data from trajectories that actually enter the encounter region, and reveal chaotic layers near the boundary.

Consider the exact two-body problem with masses $m_1,m_2$. In the inertial centre of mass frame, its internal kinetic energy is $K=m|\dot{\mathbf{r}}{|}^2/2$, where $\mathbf{r}$ is their relative position, $m=m_1{m}_2/(m_1+m_2)$ is the reduced mass, and $U=-\kappa /r$ with $\kappa =G_N{m}_1{m}_2$. Take a noncollision elliptic Kepler orbit with semi-major axis $a$ and orbital eccentricity $0<e<1$. Write $\eta$ for eccentric anomaly, to avoid confusing it with total energy $E$. With mean motion $n=\sqrt{\kappa /(m{a}^3)}$ and orbital period $P=2\pi /n$, Kepler's equation and the orbit geometry giveThe time origin is periapsis. Since $dt=(1-e\cos \eta )d\eta /n$, the time average of $1/r$ over one orbital period is exactly $1/a$. Therefore $⟨K{⟩}_P=\kappa /(2a)$ and $⟨U{⟩}_P=-\kappa /a$, independently of orbital eccentricity. Equal steps in eccentric anomaly are not equal time steps; unweighted averaging of a drawing of the ellipse gives the wrong answer.

The original dimensionless example in both diagrams has $m=a=\kappa =1$ and $e=0.6$. At periapsis, $r=0.4$, $K=2$, $U=-2.5$; at apoapsis, $r=1.6$, $K=0.125$, $U=-0.625$. Both have $E=-0.5$, but $2K+U$ is respectively $1.5$ and $-0.375$. Only the orbital period averages are $⟨K{⟩}_P=0.5$ and $⟨U{⟩}_P=-1$. The instantaneous zero of $2K+U$ occurs at $r=a$, not at every point of the orbit.

The endpoint term can also be evaluated exactly: $G=m\mathbf{r}\cdot \dot{\mathbf{r}}=mn{a}^{2}e\sin \eta$. Starting at periapsis, a window ending at eccentric anomaly $\eta$ hasAt $T=P/4$, Kepler's equation gives $\eta \simeq 2.09133$ and a positive residual $\simeq 0.3314\kappa /a$. At $T=P/2$ or $P$, $G$ returns to zero and the residual vanishes exactly. The lower panel shows the running residual and its endpoint expression coinciding; the shaded quarter-orbital period window explains why a short observation need not satisfy the infinite-time theorem.