diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-21 20:48:03 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-21 20:48:03 +0000 |
| commit | 9e39b7b6b7a74414551a70ea8ccdb11a98023ef5 (patch) | |
| tree | 83ff6e83d990431ebb61a15c9cb20c99a8ddf4d1 /liouville.m | |
| parent | 23dface6fc5ff3b0cf280846b8cacae26b1608cd (diff) | |
| download | multi_mode_eit-9e39b7b6b7a74414551a70ea8ccdb11a98023ef5.tar.gz multi_mode_eit-9e39b7b6b7a74414551a70ea8ccdb11a98023ef5.zip | |
move field strength from the for loop
Diffstat (limited to 'liouville.m')
| -rw-r--r-- | liouville.m | 17 |
1 files changed, 9 insertions, 8 deletions
diff --git a/liouville.m b/liouville.m index 5a8cef2..29b07f9 100644 --- a/liouville.m +++ b/liouville.m @@ -50,6 +50,14 @@ atom_properties.L0m=L0m; atom_properties.polarizability_m=polarizability_m; atom_properties.dipole_elements=dipole_elements; +%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.linear = E_field_zero + (0.00001+0.00001i)*E_field_probe + (0.00001+0.00001i)*E_field_drive; +E_field.right = E_field_zero + (1.00001+0.00001i)*E_field_probe + (1.00001+0.00001i)*E_field_drive; +E_field.left = E_field_zero + (0.00001+0.00001i)*E_field_probe + (0.00001+0.00001i)*E_field_drive; + fprintf (stderr, "tuning laser in forloop to set conditions vs detuning\n"); fflush (stderr); for detuning_p_cntr=1:N_detun_steps+1; @@ -59,14 +67,7 @@ for detuning_p_cntr=1:N_detun_steps+1; wp=wp0+detuning_p; wm=wd-(wp-wd); %modulation_freq=[0, wp, wd, wm, -wp, -wd, -wm, wp-wd, wd-wp]; - %E_field =[0, Ep, Ed, Em, Epc, Edc, Emc, 0, 0 ]; - 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.linear = E_field_zero + (0.00001+0.00001i)*E_field_probe + (0.00001+0.00001i)*E_field_drive; - E_field.right = E_field_zero + (1.00001+0.00001i)*E_field_probe + (1.00001+0.00001i)*E_field_drive; - E_field.left = E_field_zero + (0.00001+1.00001i)*E_field_probe + (0.00001+1.00001i)*E_field_drive; + modulation_freq=[0, wp, wd, -wp, -wd, wp-wd, wd-wp]; freq_index=freq2index(wp,modulation_freq); atom_field_problem.E_field = E_field; |