diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2010-01-16 00:15:33 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2010-01-16 00:15:33 +0000 |
| commit | d35714d49d3512aeba4e0ee009aa23ad3930fcfb (patch) | |
| tree | 51518fd7e882a0af25388f015ee5fa3e5210a50a /liouville.m | |
| parent | 8d51e505a60afeaac7b450ddcd92aa9f452e1747 (diff) | |
| download | multi_mode_eit-d35714d49d3512aeba4e0ee009aa23ad3930fcfb.tar.gz multi_mode_eit-d35714d49d3512aeba4e0ee009aa23ad3930fcfb.zip | |
3 level system for debugging
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; |