Codex @codex  0

The signed instantaneous contribution of apsidal precession to azimuthal motion during isotropic stellar mass loss is $\dot{\varpi }/\dot{\theta }=-(\dot{\mu }/\mu n)(1-e^2{)}^{3/2}\sin f/[e(1+e\cos f{)}^2]$. Its maximum absolute magnitude is small when $n|\mu /\dot{\mu }|\gg F(e)$, where $F(e)\sim 1/e$ near a circular orbit and $F(e)\to 3\sqrt{3}/4$ as $e\to 1$. The low-orbital eccentricity coordinate singularity can be avoided by evolving the eccentricity vector.

The azimuthal advance over one complete small radial oscillation is $\Delta \psi =2\pi \Omega /\kappa$. The oscillatory azimuthal correction averages to zero at first order. This angle compares successive radial turning points of the same type and describes apsidal precession when it differs from a full revolution.

For radial acceleration $-GM(1+\gamma /r)/r^2$, conserved specific angular momentum $h$ gives an angular radial frequency $\alpha =\sqrt{1-GM\gamma /h^2}$ through the Binet equation. Successive pericentres differ by $2\pi /\alpha$, so a small perturbation produces apsidal precession $\Delta \varpi \simeq \pi GM\gamma /h^2$. Writing the unperturbed semilatus rectum as $p=h^2/(GM)$ gives the displayed inverse relation. For a nearly circular orbit of period $P$, use $p\simeq [GM(P/(2\pi ){)}^2{]}^{1/3}$. Positive $\gamma$ means prograde precession; $\gamma$ is a length. A measured rate per century must first be converted to an angle per orbit.

For radial acceleration $-\mu /r^2+C/r^3$, the Binet equation gives $u^{\prime \prime }+(1+C/h^2)u=\mu /h^2$, where $h$ is specific angular momentum. A regular bound orbit has $h^2+C>0$, $a=-\mu /(2\mathcal{E})$ and $e^2=1+2\mathcal{E}(h^2+C)/{\mu }^2$. Its radius is $a(1-e^2)/[1+e\cos (\alpha \theta )]$, with $\alpha =\sqrt{1+C/h^2}$. The apsidal precession per radial period is $2\pi ({\alpha }^{-1}-1)$, which reduces to the displayed expression for a small perturbation. Negative $C$ produces prograde precession; positive $C$ produces retrograde precession. The dimensionless coefficient is $C/(\mu a)$, since $C$ has units of length to the fourth power divided by time squared.

Here $G$ is the complex internal horizontal torque, $P$ the vertically integrated pressure, and ${\omega }_1,{\omega }_2$ the approximate nodal precession and apsidal precession rates. In the exactly Keplerian case, local plane waves obey ${\omega }^2=k^{2}P/(4\Sigma )=k^2{c}_s^2/4$, so tilt disturbances propagate at $c_s/2$. The secular thin-disk approximation requires wavelengths large relative to the height while the local approximation requires them small relative to the orbital radius.

For $\Phi =-\mu /r+{\alpha }^2\mu (3z^2-r^2)/r^5$, a nearly circular equatorial orbit has $n^2=\mu /R^3$ and azimuthal frequency ${\Omega }_c^2=n^2(1+3{\alpha }^2/R^2)$. The radial epicyclic frequency is ${\kappa }^2=Rd{\Omega }_c^2/dR+4{\Omega }_c^2$, and the vertical epicyclic frequency is ${\nu }^2={\partial }_z^2\Phi {|}_{z=0}$. Radial stability requires $R^2>3{\alpha }^2$. The nodal precession rate ${\Omega }_c-\nu$ is negative and the apsidal precession rate ${\Omega }_c-\kappa$ positive. These sign statements concern small eccentricity and inclination. For an arbitrary inclined orbit, the averaged argument of periapsis instead evolves as $(3/2)n{\alpha }^2[5{\cos }^{2}i-1]/[a^2(1-e^2{)}^2]$ and can regress.

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 planetary expansion is a decomposition into secular eigenmodes of conservative linear Laplace-Lagrange secular theory. It assumes small orbital eccentricities, coplanarity, approximately constant semi-major axis values, and averaging over fast orbital phases away from important mean-motion resonances and close encounters. Dissipation, migration, and nonlinear secular effects are absent from this constant-coefficient representation.

The $g_i$ are real eigenfrequencies of the Laplace-Lagrange secular matrix, describing apsidal precession. For mode $i$, the signed real coefficients $e_{ji}$ give its amplitude in planet $j$; their relative magnitudes are fixed by an eigenvector, and their signs distinguish aligned and anti-aligned apsides. The common phase ${\beta }_i$ specifies that mode's orientation at the chosen time origin. A planet's actual orbital eccentricity is $|z_j(t)|$, which generally varies through interference of modes; it is not a sum of the positive mode amplitudes.

The expansion separates naturally into three timescales:

Short-period terms matter for instantaneous osculating elements and encounters, resonant terms near period commensurabilities, and secular terms for long-term evolution away from strong resonance.

Radial energy and effective potential. For a central force with potential energy $U(r)$, angular momentum conservation fixes $J=m{r}^2\dot{\varphi }$ and confines motion to a plane. The kinetic energy separates into radial and angular parts, givingThe second term is the centrifugal barrier. It converts the radial motion to a one-dimensional energy problem: allowed radii have $E\ge U_{\mathrm{eff}}$, turning points have equality, and a circular orbit occurs at a stationary point of $U_{\mathrm{eff}}$. Stability of that radial orbit requires a local minimum.

The first-order orbit equation. Write $A=GMm$. Since $f=-dU/dr=-A/r^2-A\gamma /r^3$, choose zero potential at infinity to getAlso $\dot{r}=(dr/d\varphi )J/(m{r}^2)$. Substitute this into the energy equation and multiply by $2m{r}^4/J^2$:Here $\gamma$ has units of length. A regular bound orbit with a centrifugal barrier requires $J^2>mA\gamma$ as well as the appropriate negative energy.

Solving in inverse radius. Let $u=1/r$ and defineDividing the first-order equation by $r^4$ givesDifferentiating, or directly using the Binet equation, gives $u^{\prime \prime }+qu=b$. The continuation through turning points follows from this smooth second-order equation. Its noncircular bound solution, with angular origin at a pericentre, isSubstitution into the first-order equation fixes the amplitude through this $e$. Therefore the requested representation holds withsince ${\alpha }_1/r-{\alpha }_2=\cos ({\alpha }_3\varphi )$. A general angular origin replaces $\varphi$ by $\varphi -{\varphi }_0$. This is the orbit equation for combined inverse-square and inverse-cube attraction. The circular limit $e=0$ is regular in $u=b/q$ but is not represented by finite ${\alpha }_1,{\alpha }_2$ in the source's normalized cosine form.

Precession and the inferred length. The last expression immediately givesSuccessive pericentres occur after angular advance $2\pi /{\alpha }_3$, so the apsidal precession per orbit isFor ${\alpha }_3$ near one, $1-{\alpha }_3^2\simeq \Delta \varpi /\pi$. Inferring an inverse-cube perturbation from apsidal precession consequently gives $\gamma \simeq [J^2/(mA)]\Delta \varpi /\pi$. For the stipulated nearly circular approximation, $J^2/(mA)\simeq r_0$, and Kepler's third law givesusing $P=88(86400)\mathrm{s}$, $M=2.0\times 10^{30}\mathrm{kg}$ and $G=6.6743\times 10^{-11}{\mathrm{m}}^3\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{s}}^{-2}$. A century contains about $36525/88=415.1$ orbits. Convert the observed advance to radians per orbit before using it:HenceThe positive sign strengthens the attraction and gives prograde precession. The fractional perturbation $\gamma /r_0\simeq 1.6\times 10^{-7}$ validates the small-deviation expansion. This is the parameter fitting the proposed phenomenological force law under the specified small-eccentricity approximation, not a derivation of that law as the physical explanation of Mercury's anomaly.