From 9e39b7b6b7a74414551a70ea8ccdb11a98023ef5 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Mon, 21 Dec 2009 20:48:03 +0000 Subject: move field strength from the for loop --- liouville.m | 17 +++++++++-------- 1 file changed, 9 insertions(+), 8 deletions(-) (limited to 'liouville.m') 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; -- cgit v1.2.3