“…where r ij = r i − r j are the position vectors connecting moments i and j, µ 0 is the vacuum permeability, γ i are the gyromagnetic ratios [12], Ĵi = ( Ĵi,x , Ĵi,y , Ĵi,z ) are the angular momentum vector operators, n ij = r ij /|r ij |, and the prefactor (1/2) avoids double counting. Quadrupolar coupling to nuclear spins with I ≥ 1 from finite electric field gradients has recently been discussed in detail [20,21], however, these interactions have been neglected in our model since they are relatively weak for 7 Li (I = 3/2) so will only affect the µ + polarization at longer times [6,7].…”