4.3.1.- S=1/2 con pérdidas y ganancias - Schrödinger

Sminus={{0,0},{1,0}};(*operador de Linblad de pérdidas. SplusDagger*)
Splus={{0,1},{0,0}};(*operador de Linblad de ganancias. SminusDagger*)
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 pérdidas y ganancias para los cuatro componentes de Ro*)
RHSpg11=Simplify[(-I (H2.Ro-Ro.H2)-(p/2) ((Splus .Sminus .Ro)+(Ro.Splus.Sminus) -2 (Sminus.Ro .Splus))-(g/2) ((Sminus .Splus .Ro)+(Ro.Sminus .Splus) -2 (Splus. Ro.Sminus)))[[1,1]] ,{p∈Reals,g∈Reals}]; 
RHSpg12=Simplify[(-I (H2.Ro-Ro.H2)-(p/2) ((Splus .Sminus .Ro)+(Ro.Splus.Sminus) -2 (Sminus.Ro .Splus))-(g/2) ((Sminus .Splus .Ro)+(Ro.Sminus .Splus) -2 (Splus. Ro.Sminus)))[[1,2]] ,{p∈Reals,g∈Reals}]; 
RHSpg21=Simplify[(-I (H2.Ro-Ro.H2)-(p/2) ((Splus .Sminus .Ro)+(Ro.Splus.Sminus) -2 (Sminus.Ro .Splus))-(g/2) ((Sminus .Splus .Ro)+(Ro.Sminus .Splus) -2 (Splus. Ro.Sminus)))[[2,1]] ,{p∈Reals,g∈Reals}]; 
RHSpg22=Simplify[(-I (H2.Ro-Ro.H2)-(p/2) ((Splus .Sminus .Ro)+(Ro.Splus.Sminus) -2 (Sminus.Ro .Splus))-(g/2) ((Sminus .Splus .Ro)+(Ro.Sminus .Splus) -2 (Splus. Ro.Sminus)))[[2,2]] ,{p∈Reals,g∈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)=" ,RHSpg11] 
Print[ Subscript["Ro'",12],"(t)=" ,RHSpg12]
Print[ Subscript["Ro'",21],"(t)=" ,RHSpg21]
Print[ Subscript["Ro'",22],"(t)=" ,RHSpg22]

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

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