4.4.1.- S=1/2 con defasaje - Schrödinger

Sz={{1,0},{0,-1}};
L= Sz;(*operador de Linblad de defasaje*)
Ldg=ConjugateTranspose[L]; (*Ldagger*)
H2={{(Subscript[w, 1]),0},{0,(Subscript[w, 2])}}; (*Hamiltoniano 2-dimensional*)
Ro={{(Subscript[ro, 11]),(Subscript[ro, 12])},{(Subscript[ro, 21]),(Subscript[ro, 22])}}; 

(*Calcular el Right Hand Side de la ecuación de Lindblad con defasaje para los cuatro componentes de Ro*)
RHSd11=Simplify[(-I (H2.Ro-Ro.H2)-(d/2) ((Ldg .L .Ro)[[1,1]]+(Ro.Ldg .L)[[1,1]] -2 (L. Ro.Ldg)))[[1,1]] ,d∈Reals]; 
RHSd12=Simplify[(-I (H2.Ro-Ro.H2)-(d/2) ((Ldg .L .Ro)[[1,2]]+(Ro.Ldg .L)[[1,2]] -2 (L. Ro.Ldg)))[[1,2]] ,d∈Reals];
RHSd21=Simplify[(-I (H2.Ro-Ro.H2)-(d/2) ((Ldg .L .Ro)[[2,1]]+(Ro.Ldg .L)[[2,1]] -2 (L. Ro.Ldg)))[[2,1]] ,d∈Reals];
RHSd22=Simplify[(-I (H2.Ro-Ro.H2)-(d/2) ((Ldg .L .Ro)[[2,2]]+(Ro.Ldg .L)[[2,2]] -2 (L. Ro.Ldg)))[[2,2]] ,d∈Reals];

(*Left Hand Side = Right Hand Side. Las ecuaciones que surgen para cada componente de Ro de la ecuación de Lindblad*)
Print[ Subscript["Ro'",11],"(t)=" ,RHSd11] 
Print[ Subscript["Ro'",12],"(t)=" ,RHSd12]
Print[ Subscript["Ro'",21],"(t)=" ,RHSd21]
Print[ Subscript["Ro'",22],"(t)=" ,RHSd22]

(*Resolución de las acuaciones diferenciales acopladas obtenidas en el paso anterior*)

DSolve[{ro11'[t]==0, ro12'[t]==ro12[t](-2d-I (Subscript[w, 1])+I (Subscript[w, 2])),ro21'[t]==ro21[t](2d+I (Subscript[w, 1])-I (Subscript[w, 2])),ro22'[t]==0}, {ro11[t],ro12[t],ro21[t],ro22[t]} ,t]