diff options
| -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; |