Ressonância entre um pente de frequências e um átomo de dois níveis - programa no Maple

Voltar

Programa no Maple para o estudo da ressonância da interação entre um pente de frequências e um átomo de dois níveis.

restart;

q[22] := 2*Pi*5e6;        {taxa de relaxação da população}
q[12] := 0.5*q[22];        {taxa de relaxação da coerência}
fr := 100e6;            {taxa de repetição do laser em Hz}
Tp := 100e-15;            {largura temporal do pulso em s}
Pulsos := 50;            {quantidade de pulsos}
Omega := 0.01e*q[22]/(fr*Tp);    {frequencia de rabi por modo}
delta:=0;            {dessintonia}

rho22n:=0;            {população inicial}
sigma12n:=0;            {coerência inicial}
m:=0;


{equações de Bloch}
a:=diff(rho[22](t),t)=I*A*sigma[12](t)-I*A*conjugate(sigma[12](t))-q[22]*rho[22](t);
b:=diff(sigma[12](t),t)=(I*delta-q[12])*sigma[12](t)-I*A*(1-2*rho[22](t));


{interação com os pulsos do trem}
for i from 3 to 2*Pulsos + 3 do;

if (i mod 2)=0 then
A:=Omega;
eq:=dsolve({a,b,rho[22](0)=rho22n,sigma[12](0)=sigma12n},{rho[22](t),sigma[12](t)},type=numeric,method=classical[rk4],stepsize=20e-15);
m:=m+100e-15;
rho22n:=solve(eq(100e-15)[2],rho[22](t));
sigma12n:=solve(eq(100e-15)[3],sigma[12](t));
end if;

if (i mod 2)=1 then
A:=0;
eq:=dsolve({a,b,rho[22](0)=rho22n,sigma[12](0)=sigma12n},{rho[22](t),sigma[12](t)},type=numeric,method=classical[rk4],stepsize=10e-12);
m:=m+10e-9;
rho22n:=solve(eq(10e-9)[2],rho[22](t));
sigma12n:=solve(eq(10e-9)[3],sigma[12](t));
end if;

if (i mod 1) = 0 then
B[i,1]:=m*1e9;
B[i,2]:=rho22n;
B[i,3]:=Re(sigma12n);
B[i,4]:=Im(sigma12n);
end if;

printf("%12.10f %12.10f %12.10f %12.10f\n",B[i,1],B[i,2],B[i,3],B[i,4]);
end do;

b:=Matrix(B,2*Pulsos+3,4);
ExportMatrix("C:/users/dell/Desktop/tempo.dat",b);