sigma1=880, 1<,81, 0<<; sigma2=880,-I<,8I, 0<<; sigma3=881, 0<,80,-1<<; idt=881, 0<,80, 1<<
HRzpi=MatrixExp@-I*Pi*sigma32DL;HRzp=MatrixExp@-I*Pi2*sigma32DL; HRzm=MatrixExp@I*Pi2*sigma32DL MatrixForm
HRxp=MatrixExp@-I*Pi2*sigma12DL;
HRxm=MatrixExp@I*Pi2*sigma12DL MatrixForm HRyp=MatrixExp@-I*Pi2*sigma22DL MatrixForm HRym=MatrixExp@I*Pi2*sigma22DL MatrixForm HHmd=1Sqrt@2D 881, 1<,81,-1<<L MatrixForm H=KroneckerProduct@sigma1, sigma1D +
KroneckerProduct@sigma2, sigma2D +KroneckerProduct@sigma3, sigma3D HSqrSWAP=MatrixExp@-I*Pi4*2*14*HDL MatrixForm
Exp@I Pi4DKroneckerProduct@idt, HmdD.SqrSWAP.
KroneckerProduct@Rzpi, idtD.SqrSWAP.KroneckerProduct@Rzp, Rxp.RzpD
%SimplifyMatrixForm