now we treat atom as polarization aware one

1 files changed, 9 insertions, 5 deletions

diff --git a/liouville.m b/liouville.m
index 9bf5753..c6603b4 100644
--- a/liouville.m
+++ b/liouville.m
@@ -10,7 +10,8 @@ useful_functions;
 useful_constants;
 
 % load atom energy levels and decay description
-four_levels;
+four_levels_with_polarization;
+%four_levels;
 %three_levels;
 %two_levels;
 
@@ -55,7 +56,9 @@ for detuning_p_cntr=1:N_detun_steps+1;
 	%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    ];
+	E_field.linear =[0,  0 , 0 ,  0  , 0  , 0,     0    ];
+	E_field.right  =[0,  0 , Ed,  0  , Edc, 0,     0    ];
+	E_field.left   =[0,  Ep, 0 ,  Epc, 0  , 0,     0    ];
 	freq_index=freq2index(wp,modulation_freq);
 
 	atom_field_problem.E_field          = E_field;
@@ -71,11 +74,12 @@ endfor
 
 % once we define all problems the main job is done here
 %kappa_p=cellfun( @susceptibility_steady_state_at_freq, problems_cell_array);
-kappa_p=parcellfun(2, @susceptibility_steady_state_at_freq, problems_cell_array);
+%kappa_p=parcellfun(2, @susceptibility_steady_state_at_freq, problems_cell_array);
+[xi_linear, xi_left, xi_right]=parcellfun(2, @susceptibility_steady_state_at_freq, problems_cell_array);
 
 
-figure(1); plot(detuning_freq, real(kappa_p)); title("probe dispersion");
-figure(2); plot(detuning_freq, imag(kappa_p)); title("probe absorption");
+figure(1); plot(detuning_freq, real(xi_left)); title("probe dispersion");
+figure(2); plot(detuning_freq, imag(xi_left)); 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");