diff options
Diffstat (limited to 'compass.m')
| -rw-r--r-- | compass.m | 41 |
1 files changed, 23 insertions, 18 deletions
@@ -12,10 +12,10 @@ useful_constants; basis_transformation; % load subroutines % 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 @@ -55,19 +55,20 @@ atom_properties.dipole_elements=dipole_elements; Ed=.1; Edc=conj(Ed); Ep=0.8*Ed; Epc=conj(Ep); -%modulation_freq=[0, wp, wd, -wp, -wd, wp-wd, wd-wp]; -E_field_drive =[0, 0 , Ed, 0 , Edc, 0, 0 ]; -E_field_probe =[0, Ep, 0 , Epc, 0 , 0, 0 ]; -E_field_zero =[0, 0 , 0 , 0 , 0 , 0, 0 ]; -E_field_lab.linear = E_field_zero + (1.00000+0.00000i)*E_field_probe + (1.00000+0.00000i)*E_field_drive; -E_field_lab.right = E_field_zero + (0.00000+0.00000i)*E_field_probe + (0.00000+0.00000i)*E_field_drive; -E_field_lab.left = E_field_zero + (0.00000+0.00000i)*E_field_probe + (0.00000+0.00000i)*E_field_drive; +%light_positive_freq = [wp, wd, wp-wd]; +E_field_drive = [0 , Ed, 0 ]; +E_field_probe = [Ep, 0 , 0 ]; +E_field_zero = [0 , 0 , 0 ]; +E_field_lab_pos_freq.linear = E_field_zero + (1.00000+0.00000i)*E_field_probe + (1.00000+0.00000i)*E_field_drive; +%E_field_lab_pos_freq.right = E_field_zero + (0.00000+0.00000i)*E_field_probe + (0.00000+0.00000i)*E_field_drive; +%E_field_lab_pos_freq.left = E_field_zero + (0.00000+0.00000i)*E_field_probe + (0.00000+0.00000i)*E_field_drive; + % 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=0/4; theta=pi/2; -theta=pi/4; @@ -89,17 +90,21 @@ for phi=phis; detuning_p=detuning_freq(detuning_p_cntr); wp=wp0+detuning_p; wm=wd-(wp-wd); - %modulation_freq=[0, wp, wd, wm, -wp, -wd, -wm, wp-wd, wd-wp]; - modulation_freq=[0, wp, wd, -wp, -wd, wp-wd, wd-wp]; + + light_positive_freq=[ wp, wd, wp-wd]; % we define light as linearly polarized % where phi is angle between light polarization and axis x - [E_field_lab.x, E_field_lab.y] = rotXpolarization(phi, E_field_lab.linear, modulation_freq); - E_field_lab.z=E_field_zero; + [E_field_lab_pos_freq.x, E_field_lab_pos_freq.y] = rotXpolarization(phi, E_field_lab_pos_freq.linear); + E_field_lab_pos_freq.z=E_field_zero; + % now we transfor x,y,z, to x',y', and z' with respect to magnetic field az z' axis - coord_transf_m = lin2circ * oldlin2newlin(theta); - E_field.right = coord_transf_m(1,1)*E_field_lab.x + coord_transf_m(1,2)*E_field_lab.y + coord_transf_m(1,3)*E_field_lab.z; - E_field.left = coord_transf_m(2,1)*E_field_lab.x + coord_transf_m(2,2)*E_field_lab.y + coord_transf_m(2,3)*E_field_lab.z; - E_field.linear = coord_transf_m(3,1)*E_field_lab.x + coord_transf_m(3,2)*E_field_lab.y + coord_transf_m(3,3)*E_field_lab.z; + E_field_pos_freq=xyz_lin2atomic_axis_polarization(theta, E_field_lab_pos_freq); + + % we calculate dc and negative frequiencies as well as amplitudes + [modulation_freq, E_field] = ... + light_positive_frequencies_and_amplitudes2full_set_of_modulation_frequencies_and_amlitudes(... + light_positive_freq, E_field_pos_freq); + freq_index=freq2index(wp,modulation_freq); atom_field_problem.E_field = E_field; |