where Y 0 2 is the normalized spherical harmonic with ℓ = 2 and m = 0, r ij (t) is the nuclear spin separation, and θ ij (t) is the polar angle of the direction r ij with respect to laboratory axes, assuming that the applied static magnetic field B 0 is parallel to e z . Then, the NMR relaxation rate R 1 is calculated using