diff options
Diffstat (limited to 'test.m')
| -rw-r--r-- | test.m | 113 |
1 files changed, 104 insertions, 9 deletions
@@ -11,14 +11,91 @@ useful_functions; useful_constants; % load atom energy levels and decay description -four_levels; -%three_levels; -%two_levels; +%four_levels; +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +Nlevels=4; +w1=1e9; +w2=0; +w_hpf=6800; +w3=w_hpf; +w4=w_hpf+.1; % separation of levels |3> and |4> somewhat like Zeeman splitting +w12=w1-w2; +w13=w1-w3; + +% +% ----------- |1> +% / \ +% E_d / \ +% / \ E_p +% / \ +% -------- |3> \ +% -------- |4> \ +% \ +% ___________ |2> + + +% unperturbed Hamiltonian energy levels +levels_energy=[ w1, 0, w3, w4]; +levels_energy=levels_energy*hbar; +H0=zeros(Nlevels); +H0=diag(levels_energy); +%for i=1:Nlevels + %H0(i,i)=levels_energy(i); +%endfor + +% decay matrix g(i,j) correspnds to decay from i-->j +gamma=6; +gamma_23=.001; +g_decay=zeros(Nlevels); +g_decay(1,2)=gamma; %upper level decay +g_decay(1,3)=gamma; %upper level decay +g_decay(1,4)=gamma; %upper level decay +g_decay(3,2)=gamma_23; % lower levels mixing +g_decay(2,3)=gamma_23; % lower levels mixing +g_decay(4,2)=gamma_23; % lower levels mixing +g_decay(2,4)=gamma_23; % lower levels mixing + +%defasing matris +g_deph=0; +g_dephasing=zeros(Nlevels); +g_dephasing(1,2)=g_deph; +g_dephasing(2,1)=g_dephasing(1,2); +g_dephasing(1,3)=g_deph; +g_dephasing(3,1)=g_dephasing(1,3); +g_dephasing(1,4)=g_deph; +g_dephasing(4,1)=g_dephasing(4,1); + + + +% dipole matrix +dipole_elements=zeros(Nlevels); +dipole_elements(1,2)=1; +dipole_elements(2,1)=dipole_elements(1,2); +dipole_elements(1,3)=1; +dipole_elements(3,1)=dipole_elements(1,3); +dipole_elements(1,4)=1; +dipole_elements(4,1)=dipole_elements(1,4); +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + + +%EM field definition +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + +Ep=0.01; %probe +Epc=conj(Ep); +Ed=.1; %drive +Edc=conj(Ed); +Em=-Ep; % opposite sideband (resulting from EOM modulation of drive) +Emc=conj(Em); +%wd=w13; +%wp=w12; +%wm=wd-(wp-wd); +%modulation_freq=[0, wp, wd, wm, -wp, -wd, -wm, wp-wd, wd-wp]; +E_field =[0, Ep, Ed, Em, Epc, Edc, Emc, 0, 0 ]; +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -% load EM field description -field_description; -Nfreq=length(modulation_freq); +%Nfreq=length(modulation_freq); @@ -31,7 +108,10 @@ detuning_freq=zeros(1,N_detun_steps+1); kappa_p =zeros(1,N_detun_steps+1); kappa_m =zeros(1,N_detun_steps+1); -wp0=w12; +w_pf1=1e9; +w_hpf_ground=6800; +wd=w_pf1-w_hpf_ground; +wp0=w_pf1; wp=wp0+detuning_p_min; wm=wd-(wp-wd); %modulation_freq=[0, wp, wd, wm, -wp, -wd, -wm, wp-wd, wd-wp]; @@ -44,6 +124,16 @@ Nfreq=length(modulation_freq); [N, rhoLiouville_w, rhoLiouville_r, rhoLiouville_c]=unfold_density_matrix(Nlevels,Nfreq); rhoLiouville=zeros(N,1); +dipoles=dipole_elements; +Efld=E_field; +clear dipole_elements; +dipole_elements.right=dipoles; +dipole_elements.left=dipoles; +dipole_elements.linear=dipoles; +clear E_field; +E_field.right=Efld; +E_field.left=Efld; +E_field.linear=Efld; % calculate E_field independent properties of athe atom % to be used as sub matrix templates for Liouville operator matrix [L0m, polarizability_m]=L0_and_polarization_submatrices( ... @@ -73,8 +163,13 @@ L_new=L; rhoLiouville_new=rhoLiouville; % uncomment to update reference file -% save 'L_and_rhoL_referenced.mat' L +%save 'L_and_rhoL_referenced.mat' L rhoLiouville + +clear L rhoLiouville; + +load 'L_and_rhoL_referenced.mat' +diff_with_reference_L=sum(sum(abs(L_new-L))) -diff_with_reference=sum(sum(abs(L_new-L))) +diff_with_reference_rhoL=(sum(abs(rhoLiouville_new -rhoLiouville))) elapsed_time = etime (clock (), t0) |