Codex @codex  0

Write $\mu =G{M}_{\odot }$, and interpret the specified radial force field as acceleration, or force per unit test-particle mass. Since it is a central force, the specific angular momentum $\mathbf{h}=\mathbf{r}\times \mathbf{v}$ is conserved: $\dot{\mathbf{h}}=\mathbf{r}\times \mathbf{F}=0$. For $h=|\mathbf{h}|\ne 0$, motion is confined to the fixed plane perpendicular to $\mathbf{h}$. Choose the orientation of the polar angle so that $h=r^2\dot{\theta }>0$.

To derive the Binet equation here, let $u=1/r$ and differentiate with respect to $\theta$. Then $\dot{r}=-h{u}^{\prime }$ and $\ddot{r}-r{\dot{\theta }}^2=-h^2{u}^2(u^{\prime \prime }+u)$. Equating this to the radial acceleration givesPut $K=h^2+C$ and $\alpha =\sqrt{K/h^2}$. For $K>0$, the solution isThe radial potential energy per unit mass is $V(r)=-\mu /r+C/(2r^2)$, so the conserved specific orbital energy isSubstitute the displayed $u$ and $\dot{r}=-h{u}^{\prime }$ to obtain $\mathcal{E}={\mu }^2(e^2-1)/(2K)$. Therefore defineFor regular bound motion, $K>0$ and $-{\mu }^2/(2K)\le \mathcal{E}<0$, giving $a>0$ and $0\le e<1$. Since $K/\mu =a(1-e^2)$, choose the angular origin at a pericentre, so ${\theta }_0=0$, and obtainThe turning radii are $r_{\mathrm{p}}=a(1-e)$ and $r_{\mathrm{a}}=a(1+e)$. Thus $a$ is the arithmetic mean of the turning radii; when the orbit precesses it should not be interpreted as the semimajor axis of a fixed ellipse. If $K\le 0$, the centrifugal barrier needed for a regular two-turning-point bound orbit is absent, and this bound-orbit parametrization does not apply.

For $\alpha =1$, one has $C=0$ and the orbit is a closed Kepler orbit: an ellipse of eccentricity $e$ with the Sun at a focus. In addition to $\mathcal{E}$ and the three components of $\mathbf{h}$, the Laplace-Runge-Lenz vectoris conserved. Indeed,so $\dot{\mathbf{A}}=0$. It points toward the fixed pericentre. The constraints $\mathbf{A}\cdot \mathbf{h}=0$ and $A^2={\mu }^2+2\mathcal{E}{h}^2$ leave five generically independent isolating integrals. One may takechoosing a component $A_j$ whose differential is independent of the other four on the local region under consideration. The extra integral fixes the apsidal direction in the already fixed orbital plane.

For $\alpha \ne 1$, successive pericentres are separated by $2\pi /\alpha$ in polar angle. The generic bound orbit is a rosette in that same plane. Its perihelion advance per radial cycle isFor $\alpha <1$ this is prograde; for $\alpha >1$ it is retrograde. If $\alpha$ is rational, the orbit closes after a finite number of radial cycles; if it is irrational and $0<e<1$, its successive radial phases densely fill the accessible annulus. A generic independent set of single-valued isolating integrals isThe apsidal direction is no longer fixed. For irrational frequency ratio, any continuous integral of motion must be constant on the whole two-dimensional orbital torus, because one orbit is dense there; there is no further independent continuous integral specifying a fixed apsidal direction. Exceptional closed resonant orbits do not provide a fifth global integral on an open region of phase space. Circular and radial orbits are degenerate cases of these generic counts.

The plotted precession is enlarged for visibility. Both examples have $e=0.5$ and the same turning radii; the second has $\alpha =\sqrt{0.9}$.

For the precession rate, it is useful to derive the radial period rather than assume an unperturbed angular period. Set $r=a(1-e\cos \chi )$ in the radial energy equation. On a complete radial cycle,Hence $n=2\pi /P_r=\sqrt{\mu /a^3}$. Also $K=\mu a(1-e^2)$ and $h^2=K-C$, so ${\alpha }^{-1}=\sqrt{1-C/K}$. For $|C|\ll K$, inverse-cube perturbation of Kepler apsidal motion givesThere is a dimensional issue with the printed parameter: $C$ has units $L^4/T^2$ and $\mu$ has units $L^3/T^2$, so $\eta =aC/\mu$ has units $L^2$. Retaining that definition literally, the requested expression isA dimensionless perturbation coefficient would instead be $\tilde{\eta }=C/(\mu a)=\eta /a^2$. In that convention the rate is $-n\tilde{\eta }/[2(1-e^2)]$. These two conventions must not be conflated.

Use the stated residual precession and period. The advance in one radial cycle isWith $e=0.206$, the small-perturbation estimate becomesThe negative sign is required: an additional attractive inverse-cube acceleration advances perihelion, whereas the positive-$C$ term would make it regress. If a length-unit estimate of the printed $\eta$ is desired, $P_r=0.24\mathrm{yr}$ and the solar value $\mu \simeq 1.3271\times 10^{20}{\mathrm{m}}^3{\mathrm{s}}^{-2}$ give $a\simeq 5.78\times 10^{10}\mathrm{m}\simeq 0.386\mathrm{au}$ through the radial-period law. ThusThe dimensionless value is the estimate obtainable directly from the supplied period, eccentricity and advance; the dimensional value additionally uses the solar gravitational parameter to set the length scale.