diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2010-03-23 15:14:51 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2010-03-23 15:14:51 +0000 |
| commit | 2d670145e1aa87c929658f184728e2727eefec5d (patch) | |
| tree | 8c3f61c42d1956bc4cd9421ff07a27702247a45d | |
| parent | 6c8564266dba940a290a468e8595bc71ef54bf64 (diff) | |
| download | multi_mode_eit-2d670145e1aa87c929658f184728e2727eefec5d.tar.gz multi_mode_eit-2d670145e1aa87c929658f184728e2727eefec5d.zip | |
q
| -rw-r--r-- | compass_lin.m (renamed from compass.m) | 8 | ||||
| -rw-r--r-- | field_description.m | 2 | ||||
| -rw-r--r-- | rb87_D1_line.m | 2 |
3 files changed, 8 insertions, 4 deletions
diff --git a/compass.m b/compass_lin.m index 5015bbf..2113f91 100644 --- a/compass.m +++ b/compass_lin.m @@ -68,7 +68,11 @@ E_field_lab_pos_freq.linear = E_field_zero + (1.00000+0.00000i)*E_field_probe + phi=pi*2/8; % theta is angle between lab z axis (light propagation direction) and magnetic field axis (z') theta=0/4; -theta=pi/2; +theta=65/180*pi; +%small ellipticity angle psi (0 to pi/2) +% 0 will give linear polarization +% pi/2 circular polarization +psi=pi/2; @@ -80,7 +84,7 @@ detuning_freq=[-.075, -.05, -.025, 0 , .025, .05 , .075, .1]; problem_cntr=1; N_angle_steps=15; -min_angle=0; max_angle=pi/2; +min_angle=0; max_angle=pi; phis=min_angle:((max_angle-min_angle)/N_angle_steps):max_angle; for phi=phis; for detuning_p_cntr=1:length(detuning_freq); diff --git a/field_description.m b/field_description.m index a9c41ea..3a58492 100644 --- a/field_description.m +++ b/field_description.m @@ -3,7 +3,7 @@ %EM field definition Ed=0.1; %drive Edc=conj(Ed); -Ep=0.08*Ed; %probe +Ep=0.8*Ed; %probe Epc=conj(Ep); Em=-Ep; % opposite sideband (resulting from EOM modulation of drive) Emc=conj(Em); diff --git a/rb87_D1_line.m b/rb87_D1_line.m index 6edbcdd..96e9774 100644 --- a/rb87_D1_line.m +++ b/rb87_D1_line.m @@ -105,7 +105,7 @@ dipole_elements.right+=dipole_elements.right'; maximum_dipole_elements=abs(dipole_elements.linear) + abs(dipole_elements.left) + abs(dipole_elements.right); %defasing matrix -g_deph=0; +g_deph=100; g_dephasing=zeros(Nlevels); % dephasing only for non zero dipole elements (am I right?) g_dephasing=g_deph*(abs(maximum_dipole_elements) != 0); |