 Solution

Let $K_{ij}={\partial }_j{u}_i$, $C_{ij}=⟨R_i{R}_j⟩$, $\alpha =2H/\zeta$ and $D_R=2k_{B}T/\zeta$. The configuration equation of the Hookean dumbbell model becomesHere the paper's $k$ denotes the Boltzmann constant $k_B$. Integrating the second moment equation against $R_i{R}_j$ and discarding boundary terms at infinite separation givesFor example the spring drift differentiates each factor $R_i,R_j$, while the Laplacian gives ${\nabla }_{\mathbf{R}}^2(R_i{R}_j)=2{\delta }_{ij}$. The Hookean dumbbell conformation equation isIn the row-vector convention of $\mathbf{R}\cdot \nabla \mathbf{u}$, the displayed gradient tensor is $K^T$, giving exactly the transposes in the requested equation. The dimensionless conformation tensor $A=HC/(k_{B}T)$ consequently has upper-convected derivative $-(A-I)/\lambda$, with polymer relaxation time $\lambda =\zeta /(4H)$. This normalization matters: the unscaled $C$ has equilibrium value $(k_{B}T/H)I$, rather than $I$.