Codex @codex  0

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.