Laplace-Lagrange inclination matrix ($B$)

For circular planetary orbits with angular momenta $L_j=M_j\sqrt{G{M}_{\star }{a}_j}$, the inclination matrix has positive off-diagonal entries, zero row sums, and $L_j{B}_{jk}=L_k{B}_{kj}$. Positive diagonal symmetrization of a matrix with $W=\mathrm{diag}(L_j)$ makes $W^{1/2}B{W}^{-1/2}$ a symmetric matrix. Moreover,All eigenvalues are real and nonpositive. The zero eigenvalue is a common tilt of every orbital plane; nonzero modes describe relative nodal precession.

If two nearby inner planets couple much more strongly to each other than to a distant third planet, their slow secular eigenmode has a common complex inclination. With $L_{12}=L_1+L_2$, the effective inner-to-outer coupling isAngular momentum balance makes the reverse coupling $c{L}_{12}/L_3$. The reduced matrix has eigenvalues $0$ and $-c(1+L_{12}/L_3)$. If $L_3\gg L_{12}$, the nonzero frequency is approximately $-c$, while the fast inner differential mode has frequency $-(B_{12}+B_{21})$.