 Hookean dumbbell conformation equation with anisotropic mobility

Let $C=⟨\mathbf{R}{\mathbf{R}}^T⟩$ be the unnormalized conformation tensor, $G_{ij}={\partial }_j{u}_i$ the velocity gradient, and $M/\zeta$ a symmetric positive hydrodynamic mobility matrix independent of the individual connector. The connector Fokker-Planck equation has drift $G\mathbf{R}-(2H/\zeta )M\mathbf{R}$ and diffusion matrix $(2k_{B}T/\zeta )M$. Integrating its second moment twice by parts gives the displayed upper-convected derivative, $C^{\nabla }=D_{t}C-GC-C{G}^T$. If $M$ is a polynomial in $C$, the two matrices commute and the relaxation term simplifies to $-(4H/\zeta )M(C-k_{B}TI/H)$. No Gaussian distribution assumption is required for this moment closure.