and used for numerical simulation of 2H NMR spectra of the paramagnetic samples. Hemisphere was divided by ca. 19,800 points using the ASG tiling scheme [61]. The eight irreducible spin-operators Ti;j with rank i = 1, 2 were used as basis operators [45]. A complex, non-Hermitian exchange matrix relevant for molecular dynamics was diagonalized numerically by the QR algorithm [62]. The eight-step phase cycling ð2/2 2/3 2/4Þ was used in simulation of the PIQE, APIQE and APIQE2 spectra. The cycle was two for the QE sequence.