diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-10 17:21:37 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-10 17:21:37 +0000 |
| commit | b126435d36a59b740c3b4dd0ad6afd3d2417dda7 (patch) | |
| tree | e4d2b49e33055d22fdb27fd0541409f915d87fe3 /liouville.m | |
| parent | dbbb852d00038ac7c7873783c4442724dee1031b (diff) | |
| download | multi_mode_eit-b126435d36a59b740c3b4dd0ad6afd3d2417dda7.tar.gz multi_mode_eit-b126435d36a59b740c3b4dd0ad6afd3d2417dda7.zip | |
fixes for 4 level calculation model
Diffstat (limited to 'liouville.m')
| -rw-r--r-- | liouville.m | 15 |
1 files changed, 10 insertions, 5 deletions
diff --git a/liouville.m b/liouville.m index 6fe83be..204c1a3 100644 --- a/liouville.m +++ b/liouville.m @@ -9,7 +9,8 @@ useful_functions; useful_constants; % load atom energy levels and decay description -three_levels; +four_levels; +%three_levels; %two_levels; % load EM field description @@ -34,7 +35,11 @@ wp0=w12; detuning_p=detuning_p_min+detun_step*(detuning_p_cntr-1); wp=wp0+detuning_p; wm=wd-(wp-wd); -modulation_freq=[0, wp, wd, wm, -wp, -wd, -wm, wp-wd, wd-wp]; +%modulation_freq=[0, wp, wd, wm, -wp, -wd, -wm, wp-wd, wd-wp]; +%E_field =[0, Ep, Ed, Em, Epc, Edc, Emc, 0, 0 ]; +modulation_freq=[0, wp, wd, -wp, -wd, wp-wd, wd-wp]; +E_field =[0, Ep, Ed, Epc, Edc, 0, 0 ]; +Nfreq=length(modulation_freq); % now we create Liouville indexes list [N, rhoLiouville_w, rhoLiouville_r, rhoLiouville_c]=unfold_density_matrix(Nlevels,Nfreq); @@ -65,7 +70,7 @@ rhoLiouville=L\rhoLiouville_dot; %rho_m=rhoOfFreq(rhoLiouville, 4, Nlevels, Nfreq); % opposite sideband frequency kappa_p(detuning_p_cntr)=sucseptibility(2, rhoLiouville, dipole_elements, Nlevels, Nfreq); -kappa_m(detuning_p_cntr)=sucseptibility(4, rhoLiouville, dipole_elements, Nlevels, Nfreq); +%kappa_m(detuning_p_cntr)=sucseptibility(4, rhoLiouville, dipole_elements, Nlevels, Nfreq); detuning_freq(detuning_p_cntr)=detuning_p; %kappa_p_re=real(kappa_p); @@ -74,7 +79,7 @@ detuning_freq(detuning_p_cntr)=detuning_p; endfor figure(1); plot(detuning_freq, real(kappa_p)); title("probe dispersion"); figure(2); plot(detuning_freq, imag(kappa_p)); title("probe absorption"); -figure(3); plot(detuning_freq, real(kappa_m)); title("off resonant sideband dispersion"); -figure(4); plot(detuning_freq, imag(kappa_m)); title("off resonant absorption"); +%figure(3); plot(detuning_freq, real(kappa_m)); title("off resonant sideband dispersion"); +%figure(4); plot(detuning_freq, imag(kappa_m)); title("off resonant absorption"); elapsed_time = etime (clock (), t0) |