Clear[Axs]
Clear[Ixs]
Clear[AIx]
Clear[Ro1]
Clear[Ro]
Clear[R1]
Clear[R2]
Clear[R3]
Clear[FF]
Clear[y]
Clear[H]
Clear[ge]
Clear[mb]
Clear[Bz]
Clear[Apar]
Clear[Aperp]
Clear[DNV]
Clear[H]
Clear[Detuning]
Clear[t]
(*NV papameters*)


Ro:=(({{1, 0}, {0, 0}}))

sx:=PauliMatrix[1];
sy:=PauliMatrix[2];
sz:=PauliMatrix[3];
ix:=PauliMatrix[1];
iy:=PauliMatrix[2];
iz:=PauliMatrix[3];
plus:=(({{0, 1}, {0, 0}}))
minus:=(({{0, 0}, {1, 0}}))
un:=(({{1, 0}, {0, 1}}))

R1:=MatrixExp[-0.5*I*(ix)].Ro.MatrixExp[0.5*I*(ix)]
R2:=MatrixExp[-I*3*2*(iz)].R1.MatrixExp[I*3*2*(iz)]
R3:=MatrixExp[-1*I*(ix)].R2.MatrixExp[1*I*(ix)]
R4:=MatrixExp[-I*t*3*(iz)].R3.MatrixExp[I*t*3*(iz)]
Res1:=R4
Pop=Re[Part[Res1,2,1]];

R1b:=MatrixExp[-0.5*I*(ix)].Ro.MatrixExp[0.5*I*(ix)]
R2b:=MatrixExp[-I*20*2*(iz)].R1b.MatrixExp[I*20*2*(iz)]
R3b:=MatrixExp[-1*I*(ix)].R2b.MatrixExp[1*I*(ix)]
R4b:=MatrixExp[-I*t*20*(iz)].R3b.MatrixExp[I*t*20*(iz)]
Res1:=R4+R4b
Pop=Re[Part[Res1,2,1]];

R1c:=MatrixExp[-0.5*I*(ix)].Ro.MatrixExp[0.5*I*(ix)]
R2c:=MatrixExp[-I*10*2*(iz)].R1b.MatrixExp[I*10*2*(iz)]
R3c:=MatrixExp[-1*I*(ix)].R2b.MatrixExp[1*I*(ix)]
R4c:=MatrixExp[-I*t*10*(iz)].R3b.MatrixExp[I*t*10*(iz)]
Res1:=R4+R4b+R4c
Pop=Re[Part[Res1,2,1]];

R1d:=MatrixExp[-0.5*I*(ix)].Ro.MatrixExp[0.5*I*(ix)]
R2d:=MatrixExp[-I*2*2*(iz)].R1b.MatrixExp[I*2*2*(iz)]
R3d:=MatrixExp[-1*I*(ix)].R2b.MatrixExp[1*I*(ix)]
R4d:=MatrixExp[-I*t*2*(iz)].R3b.MatrixExp[I*t*2*(iz)]
Res1:=R4+R4b+R4c+R4d
Pop=Re[Part[Res1,2,1]];


Plot[Pop, {t,0,3}, PlotPoints->100]