Rabi oscillations

F.P. Laussy

Manipulate[Plot[Evaluate[{
((Abs[(α Exp[I(Sqrt[(((Δ))^2)+(g^2)])t]+β Exp[-I(Sqrt[(((Δ))^2)+(g^2)])t])])^2),
((Abs[(1/g)(α Exp[I(Sqrt[(((Δ))^2)+(g^2)])t]((Sqrt[(((Δ))^2)+(g^2)])-Δ)-β Exp[-I(Sqrt[(((Δ))^2)+(g^2)])t]((Sqrt[(((Δ))^2)+(g^2)])+Δ))])^2),
((Abs[(α Exp[I(Sqrt[(((Δ))^2)+(g^2)])t]+β Exp[-I(Sqrt[(((Δ))^2)+(g^2)])t])])^2)+((Abs[(1/g)(α Exp[I(Sqrt[(((Δ))^2)+(g^2)])t]((Sqrt[(((Δ))^2)+(g^2)])-Δ)-β Exp[-I(Sqrt[(((Δ))^2)+(g^2)])t]((Sqrt[(((Δ))^2)+(g^2)])+Δ))])^2)
}/.{α->-((-I (Sqrt[1-(A^2)])Exp[I θ] g+A (Δ+(Sqrt[(g^2)+(Δ^2)])))/(2 (Sqrt[(g^2)+(Δ^2)]))),β->-((I (Sqrt[1-(A^2)])Exp[I θ] g+A (-Δ+(Sqrt[(g^2)+(Δ^2)])))/(2 (Sqrt[(g^2)+(Δ^2)])))}],{t,0,10},PlotRange->{0,All},
Frame->True,FrameLabel->{"Time t","|(Subscript[c, a])(t)((|)^2), |(Subscript[c, a])(t)((|)^2)"},ImageSize->800,BaseStyle->Large,PlotStyle->Thick
]
,{{Δ,0, "Detuning"},-5,5},{{A,(1/(Sqrt[2])), "Amplitude A"},0,1},{{θ,0,"Angle θ"},0,2π},{{g,1, "Coupling strength g" },0,10}]