diff options
Diffstat (limited to 'liouville.m')
| -rw-r--r-- | liouville.m | 11 |
1 files changed, 7 insertions, 4 deletions
diff --git a/liouville.m b/liouville.m index 233ba53..8280e9c 100644 --- a/liouville.m +++ b/liouville.m @@ -10,10 +10,10 @@ useful_functions; useful_constants; % load atom energy levels and decay description -rb87_D1_line; +%rb87_D1_line; %four_levels_with_polarization; %four_levels; -%three_levels; +three_levels; %two_levels; % load EM field description @@ -26,7 +26,8 @@ field_description; %tune probe frequency detuning_p=0; N_detun_steps=200; -detuning_p_min=-B_field*gmg*4; % span +/-4 Zeeman splitting +%detuning_p_min=-B_field*gmg*4; % span +/-4 Zeeman splitting +detuning_p_min=-1.1; detuning_p_max=-detuning_p_min; detuning_freq=zeros(1,N_detun_steps+1); kappa_p =zeros(1,N_detun_steps+1); @@ -61,7 +62,8 @@ E_field_lab.left = E_field_zero + (0.00000+0.00000i)*E_field_probe + (0.00000 % phi is angle between linear polarization and axis x phi=pi*2/8; % theta is angle between lab z axis (light propagation direction) and magnetic field axis (z') -theta=pi/2; +theta=0/2; +phi=0/4; % we define light as linearly polarized % where phi is angle between light polarization and axis x @@ -107,6 +109,7 @@ fflush (stderr); %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); total_absorption_vs_detuning=parcellfun(2, @total_field_absorption, problems_cell_array); +%total_absorption_vs_detuning=cellfun(@total_field_absorption, problems_cell_array); %save 'xi_vs_detuning.mat' detuning_freq xi_linear xi_left xi_right ; save '/tmp/total_absorption_vs_detuning.mat' detuning_freq total_absorption_vs_detuning; |