Codex @codex  0

The argument of periapsis is the angle in the orbital plane from the ascending node to periapsis, measured in the direction of motion. It differs from the longitude of periapsis $\varpi =\Omega +\omega$. For an inclined orbit this sum is a conventional longitude, not generally the azimuth of the projected periapsis. The argument is undefined for a circular orbit and the node-based description is singular when the orbit and reference plane coincide.

The particular, forced part of a complex inclination solution. For a fixed, circular, inclined inner planet and an exterior test ring, linear Laplace-Lagrange secular theory gives $\dot{w}=-i\nu (w-w_p)$, so $w_f=w_p$ and $w=w_p+(w(0)-w_p)e^{-i\nu t}$. The ring describes a circle about the planetary tilt in an Argand diagram; relative to the planet, its mutual inclination stays constant while its ascending node precesses.

In a right-handed reference frame, the rotation matrix from an orbital plane frame whose first axis points along the ascending node is $R_z(\Omega )R_x(I)$. Its columns are the orbital frame's orthonormal basis expressed in reference coordinates. The unit vector normal is $(\sin \Omega \sin I,-\cos \Omega \sin I,\cos I)$; its dot product with another such normal gives the cosine of their mutual inclination.

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.

For a circular orbit the origin of mean anomaly is arbitrary because there is no distinguished pericentre. Choose $M=0$ at the ascending node, so $M$ is also the argument of latitude. If a different periapsis convention is used, replace $M$ below by $M+\omega$. The orbital-frame rotation from inclination and node givesIn particular the crossing at $M=0$ has increasing $z$ for prograde $n>0$ and $0<I<\pi$, confirming the ascending node convention.

Put $ϵ=a/a_{\mathrm{pl}}$ and $c=\cos \psi =\hat{\mathbf{r}}\cdot {\hat{\mathbf{r}}}_{\star }$. The binomial series expansion of the direct disturbing function isThe indirect disturbing function subtracts $G{M}_{\star }ac/a_{\mathrm{pl}}^2$, cancelling the dipole term. The constant $G{M}_{\star }/a_{\mathrm{pl}}$ has no spatial derivative and does not perturb the particle. After omitting that constant, the leading stellar quadrupole disturbing function isThus the PDF's equality is understood modulo a dynamically irrelevant additive constant and higher multipoles.

Define $\ell =M+\Omega$ and let ${\ell }_{\star }$ be the star's orbital longitude as seen from the planet. Its orbital plane has $I_{\mathrm{pl}},{\Omega }_{\mathrm{pl}}$. The vector from planet to star is the negative of the planet's heliocentric vector: if ${\ell }_{\mathrm{pl}}$ describes that latter vector, use ${\ell }_{\star }={\ell }_{\mathrm{pl}}+\pi$. This sign matters in $\cos \psi$ but drops out of ${\cos }^2\psi$. Expanding the two rotated unit vectors givesFor example, the particle's horizontal correction is $-(I^2/2)\sin (\ell -\Omega )(-\sin \Omega ,\cos \Omega )$, while its vertical component is $I\sin (\ell -\Omega )$; their dot products with the corresponding stellar components produce these terms. No cubic terms occur because simultaneous reversal of both small orbital inclinations leaves the dot product unchanged.

The two mean longitudes are uniform for the circular Kepler orbits, so orbit averaging is uniform in both $\ell$ and ${\ell }_{\star }$. Set $c_0=\cos (\ell -{\ell }_{\star })$ and $c_2=\cos \psi -c_0$ through quadratic order. The needed averages follow by integrating sine and cosine products over $[0,2\pi {]}^2$:For the last identity, expand $c_0=\cos \ell \cos {\ell }_{\star }+\sin \ell \sin {\ell }_{\star }$ and factor the two integrations, each using $⟨{\sin }^2\ell ⟩=⟨{\cos }^2\ell ⟩=1/2$. Squaring $c_0+c_2$ now yieldswhere the error includes all fourth-order combinations of the two orbital inclinations. ThereforeMultiplying this by $G{M}_{\star }{a}^2/a_{\mathrm{pl}}^3$ gives the requested secular disturbing function. Equivalently, the exact angular factor for two circular rings is $(3{\cos }^2{I}_m-1)/8$, where $I_m$ is their mutual inclination; its small-angle expansion gives the same result. This is the double-averaged circular tidal potential.

For aligned equatorial and planetary orbital planes, set $I_{\mathrm{pl}}=0$. The secular stellar quadrupole disturbing function then has derivative ${\partial }_I{\overline{\mathcal{R}}}_{\star }=-3G{M}_{\star }{a}^{2}I/(4a_{\mathrm{pl}}^3)$. The supplied Lagrange planetary equations for nodal precession givesThe planetary quadrupole coefficient $J_2>0$ contributes ${\partial }_I{\overline{\mathcal{R}}}_{J_2}=-3G{M}_{\mathrm{pl}}{J}_2{R}_{\mathrm{pl}}^{2}I/(2a^3)$, henceBoth are retrograde nodal precession rates. Equating them cancels $n$ and gives the Laplace radiusThe planetary quadrupole dominates inside $a_L$ and the stellar tide outside it. These are small-inclination rates for planet-bound Kepler orbits in the hierarchical regime $a\ll a_{\mathrm{pl}}$. At exactly $I=0$, the node itself is undefined; the rates are the limiting precession frequencies of an infinitesimal tilt.

Take the sky frame and disk frame to be right-handed orthonormal bases, with positive angular displacement measured from $\hat{X}$ toward $\hat{Y}$. In particular, positive $Z$ is toward the observer, so material crossing the sky plane with positive $\dot{Z}$ is at the ascending node. Write $c_d=\cos {\Omega }_d$, $s_d=\sin {\Omega }_d$, $C_d=\cos I_d$, and $S_d=\sin I_d$. The disk frame's first unit vector lies along its ascending node, and its second has positive $Z$ component. The orbital-frame rotation from inclination and node therefore givesEach column is a disk-frame unit vector in sky coordinates, so $\mathbf{X}=T\mathbf{x}$ and the inverse is $\mathbf{x}=T^{\mathsf{T}}\mathbf{X}$. At the disk's ascending node, the $Z$ component of its tangential velocity is positive, confirming the sign of the tilt.

Let $T_p=R_Z({\Omega }_p)R_X(I_p)$ define a provisional planet frame whose first axis is the planet's sky-plane ascending node. Its normal is ${\hat{z}}_p=(\sin {\Omega }_p\sin I_p,-\cos {\Omega }_p\sin I_p,\cos I_p)$. Expressing that normal in the disk orthonormal basis, with $\Delta \Omega ={\Omega }_p-{\Omega }_d$, givesThe last component is the dot product of the two orbital normals, hence the cosine of the mutual inclination. The mutual ascending node points along ${\hat{z}}_d\times {\hat{z}}_p$: motion there has positive disk-normal component. SetThen $\hat{x}\prime =(A\hat{x}+B\hat{y})/\sqrt{A^2+B^2}$, and the quadrant-correct longitude of ascending node is ${\Omega }_m=atan2(B,A)$. ThusThe tangent alone cannot fix the quadrant, and ${\Omega }_m$ is undefined for exactly coincident orbital planes.

For the two rotation routes, define $\psi$ as the oriented angle, within the planet's orbital plane, from its sky-plane ascending node to its disk-plane ascending node. If ${\hat{x}}_p,{\hat{y}}_p$ are the first two columns of $T_p$, take $\psi =atan2(\hat{x}\prime \cdot {\hat{y}}_p,\hat{x}\prime \cdot {\hat{x}}_p)$. Both routes describe the same planet orthonormal basis:The left route goes through the disk frame; the right goes through the sky-node planet frame. Comparing their third columns gives the two boxed relations, while comparison of first columns fixes the extra in-plane angle.

Put $I_d=\pi /2+{\delta }_d$ and $I_p=\pi /2+{\delta }_p$. A second-order Taylor expansion gives $\cos I_m=1-[({\delta }_p-{\delta }_d{)}^2+\Delta {\Omega }^2]/2+O({\delta }^4)$, and therefore the mutual inclination of nearly edge-on orbits satisfiesAll expansion angles are in radians. For plotting in degrees this is $I_m\simeq \sqrt{(I_p-87^{\circ }{)}^2+(4^{\circ }{)}^2}$. The given planet has exact $I_m\simeq 4.119^{\circ }$ and ${\Omega }_m\simeq 76.044^{\circ }$. The exact minimum occurs at $I_p=atan2(\sin I_d\cos \Delta \Omega ,\cos I_d)\simeq 86.993^{\circ }$, with $I_{m,\min }\simeq 3.995^{\circ }$. At $I_p=85^{\circ }$ and $95^{\circ }$, the second-order values are ${\sqrt{20}}^{\circ }\simeq 4.472^{\circ }$ and ${\sqrt{80}}^{\circ }\simeq 8.944^{\circ }$.

For the disk evolution, neglect disk self-gravity, collisions that align the planes, and back-reaction on the planet. In the initial disk plane, use the complex inclination $w=I{e}^{i\Omega }$ for a ring and $w_p=I_m{e}^{i{\Omega }_m}$ for the planet. Linear Laplace-Lagrange secular theory has forced inclination $w_p$ and solution $w(a,t)=w_p[1-e^{-i\nu (a)t}]$ for initially flat rings. For an inner circular perturber, the Laplace coefficient formula is $\nu =n(M_p/M_{\star })\alpha b_{3/2}^{(1)}(\alpha )/4$, with $\alpha =a_p/a$. Its quadrupole approximation uses $b_{3/2}^{(1)}\simeq 3\alpha$, yieldingIn the $(I\cos \Omega ,I\sin \Omega )$ plot referenced to the initial disk, every ring runs clockwise on the same circle centered at $w_p$, starting at the origin. Inner rings advance more quickly. If the axes are instead referenced to the planet's plane, the relative tilt is $w-w_p=-w_p{e}^{-i\nu t}$: these circles are centered at the origin, and the mutual inclination is constant. Both conventions describe the same nodal precession; they should not be mixed.

At an intermediate epoch, rings with $\nu t\gtrsim 1$ have tilted substantially while outer rings with $\nu t\ll 1$ remain near the original disk plane. This radial variation is a planet-induced debris-disk warp. The characteristic affected radius grows like $t^{2/7}$. Conservative differential nodal precession preserves each ring's tilt relative to the planet, so it does not by itself align every orbit; an unresolved inner region can acquire a mean plane near the planetary plane through phase mixing. The secular interpretation is also described in Wyatt's planetary dynamics lectures.

The angleobeys the D'Alembert characteristic. Because reflection in the reference plane makes the disturbing function even in inclination, the nodal coefficient must be even, and thereforeAt astronomical conjunction, ${\lambda }_1={\lambda }_2={\lambda }_c$, soThus $\varphi /q$ is the angular distance from the ascending node of $M_2$ to the conjunction. When $M_2$ is at that node, ${\lambda }_2={\Omega }_2$, andUp to the chosen sign convention, $\varphi /p$ is the angular separation between the inner body and the node when the outer body crosses the reference plane. These two views explain geometrically why libration confines both conjunction and node-crossing phases.

Just after $M_2$ crosses its ascending node, it lies above the nearly fixed orbital plane of $M_1$. At a nearby conjunction, $M_1$ is inward of $M_2$, because $r_1<r_2$. The perturbing acceleration toward $M_1$ therefore has an inward component and a downward normal component. In the view along the node toward the star, the downward component opposes $M_2$'s upward vertical velocity and reduces its orbital inclination; in the face-on view, the inward component perturbs its orbital energy and shifts the conjunction phase. A conjunction before the ascending node reverses the relation between the normal force and vertical velocity. Repeated phase-correlated impulses therefore provide the restoring dynamics of an inclination-type mean-motion resonance.

Resolving $\mathbf{r}=a(\cos f{\mathbf{e}}_1+\sin f{\mathbf{e}}_2)$ into the sky coordinates givesIn particular, $z=0$ at the orbital nodes, and $dz/df=a\sin I>0$ at $f=0$, so that node is indeed the ascending node.

Put the star at the origin, take $\hat{\mathbf{z}}$ along the line of sight toward the observer, and let $\hat{\mathbf{x}}$ point North in the plane of the sky. The $\hat{\mathbf{y}}$ direction is then $90^{\circ }$ anticlockwise from North. Draw the ascending node on the sky plane at position angle $\Omega$, tilt the circular Kepler orbit through the orbital inclination $I$ about that nodal line, and place the planet an angle $f$ along the orbit from the ascending node. The direction of increasing $f$ must cross the sky plane toward $+\hat{\mathbf{z}}$ at $f=0$.

A convenient orthonormal basis in the orbital plane isHere ${\mathbf{e}}_1$ points along the nodal line and ${\mathbf{e}}_2$ fixes the required sense of motion. The planet lies at $\mathbf{r}=a(\cos f{\mathbf{e}}_1+\sin f{\mathbf{e}}_2)$; this relation also specifies all the labels required in the sketch.