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 | |
| parent | 8d51e505a60afeaac7b450ddcd92aa9f452e1747 (diff) | |
| download | multi_mode_eit-d35714d49d3512aeba4e0ee009aa23ad3930fcfb.tar.gz multi_mode_eit-d35714d49d3512aeba4e0ee009aa23ad3930fcfb.zip | |
3 level system for debugging
| -rw-r--r-- | compass.m | 12 | ||||
| -rw-r--r-- | liouville.m | 11 |
2 files changed, 13 insertions, 10 deletions
@@ -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,12 +26,12 @@ field_description; %tune probe frequency detuning_p=0; N_detun_steps=100; -detuning_p_min=-B_field*gmg*4; % span +/-4 Zeeman splitting -detuning_p_max=-detuning_p_min; +%detuning_p_min=-B_field*gmg*4; % span +/-4 Zeeman splitting +%detuning_p_max=-detuning_p_min; detuning_freq=zeros(1,N_detun_steps+1); kappa_p =zeros(1,N_detun_steps+1); kappa_m =zeros(1,N_detun_steps+1); -detun_step=(detuning_p_max-detuning_p_min)/N_detun_steps; +%detun_step=(detuning_p_max-detuning_p_min)/N_detun_steps; fprintf (stderr, "calculating atom properties\n"); fflush (stderr); @@ -65,6 +65,7 @@ E_field_lab.left = E_field_zero + (0.00000+0.00000i)*E_field_probe + (0.00000 phi=pi*2/8; % theta is angle between lab z axis (light propagation direction) and magnetic field axis (z') theta=pi/2; +theta=0; % we define light as linearly polarized E_field_lab.x=cos(phi)*E_field_lab.linear; @@ -102,7 +103,6 @@ for phi=phis; wp0=w_pf1; wd=w_pf1-w_hpf_ground; - detuning_p=detuning_p_min+detun_step*(detuning_p_cntr-1); detuning_p=detuning_freq(detuning_p_cntr); wp=wp0+detuning_p; wm=wd-(wp-wd); 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; |