spin = 3/2;
wq = (2*Pi)*28.1*10^6;
alphaQ = 0; betaQ = 0; gammaQ = 0; eta = 0; cQ = 2*wq;
H0 = HQ[spin];

(* Sequence of three experiments to generate the pseudopure state 11 *)

  ρEq = -H0; (* equilibrium *)

w1 = wq/1000; thetaRF = Pi/2; wRF = wq;
phiRF = 0; phase = 0; hx = HRF[spin];
phiRF = Pi/2; phase = Pi/2; hy = HRF[spin];
H1 = hx + hy;
TP = Pi/(Sqrt[3]*w1);
  ρS = Pulse[-H0]; (* single quantum *)

w1 = wq/1000; thetaRF = Pi/2; wRF = wq/2;
phiRF = 0; phase = 0; hx = HRF[spin];
phiRF = Pi/2; phase = Pi/2; hy = HRF[spin];
H1 = hx + hy;
TP = Pi*wq/(2*Sqrt[3]*w1^2);
  ρD = Pulse[-H0]; (* double quantum *)

ρ3 = (ρEq + ρS + ρD); (* pseudopure state |11> *)