The routine EnsembleAverage may be used to generate arbitrary distributions of nuclear spin interactions. In an example posted on the forum, I used EnsembleAverage to generate the spectrum of 133Cs (I=7/2) in an inhomogeneous gel. The same notebook includes a simulation of quadrupolar echoes in spin-7/2 with curious fractional timings.
Sep 19