In these equations, m(r i ) and m(r j ) are the moments at sites i and j, and rij = (r i -r j )/|r i -r j | is the unit vector joining the positions i and j. Another essential ingredient of the theoretical setup is the presence of the radiation feedback field from the coil and acting on the nuclear spins.