diff options
Diffstat (limited to 'liouville.m')
| -rw-r--r-- | liouville.m | 14 |
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"); |