Double-averaged circular tidal potential (${\overline{\mathcal{R}}}_{\star }=\frac{G{M}_{\star }{a}^2}{8a_p^3}(3{\cos }^2{I}_m-1)$)

For two circular Kepler orbits with mutual inclination $I_m$, independently average the stellar quadrupole disturbing function over both mean anomalies. A circular ring with unit normal $\mathbf{k}$ has $⟨{\hat{r}}_i{\hat{r}}_j⟩=({\delta }_{ij}-k_i{k}_j)/2$. Contracting the two ring tensors gives $⟨{\cos }^2\psi ⟩=[1+(\mathbf{k}\cdot {\mathbf{k}}_p{)}^2]/4$, proving the displayed exact angular factor. For small orbital inclinations relative to a common plane, $I_m^2=I^2+I_p^2-2I{I}_p\cos ({\Omega }_p-\Omega )+O(I^4)$. Thus the secular angular factor is $1/4-3(I^2+I_p^2)/8+3I{I}_p\cos ({\Omega }_p-\Omega )/4$. This potential generates the slow stellar contribution to satellite nodal precession.