“…As the spin and spatial variables cannot be separated in the stochastic Liouville approach, one constructs a basis being an outer product of the spin and rotational variables: |Oα)=|sans-serifΣ,σ)⊗|L,K,M), where Σ ranges from 1 to (2S+1), i.e., 3 for 14 N, σ=−sans-serifΣ,…,sans-serifΣ while the quantum numbers L , K , M correspond to the indices of Wigner rotation matrices; for practical calculations it is sufficient to set L = 1, …, 8, K = − L , …, L , M = − L , …, L . Then the relaxation rate, R1HN, is given as [22,49,50,51,52]: R1HN=Re{[T11]+[M]−1[T11]}where the matrix [ M ] is defined in the {|Oα)(Oβ|}…”