From 13f7d144ea868547b7f21000f1fddd836f8b38e9 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 15 Nov 2012 16:31:17 -0500 Subject: fields definitions are placed --- .../realistic_Rb_and_fields.xmds | 641 +++++++++++++++++++++ 1 file changed, 641 insertions(+) create mode 100644 xmds2/realistic_Rb_and_fields/realistic_Rb_and_fields.xmds (limited to 'xmds2') diff --git a/xmds2/realistic_Rb_and_fields/realistic_Rb_and_fields.xmds b/xmds2/realistic_Rb_and_fields/realistic_Rb_and_fields.xmds new file mode 100644 index 0000000..bfed139 --- /dev/null +++ b/xmds2/realistic_Rb_and_fields/realistic_Rb_and_fields.xmds @@ -0,0 +1,641 @@ + + + + --Ndens=1e15 --Lcell=10.0e-2 --Temperature=1e-3 --Pwidth=0.4e-6 --delta1=0 --delta2=0 --delta3=0 --E1o=0 --E2o=1e2 --E3o=0 --E4o=0 + + + + + + realistic_Rb_and_fields + + Eugeniy Mikhailov + + License GPL. + + Solving simplified Rb atom model + with fields propagation along spatial axis Z + with Doppler broadening. + + + We assume four-wave mixing condition when w3-w4=w2-w1 i.e. fields E3 and E4 drive the same + resonance as fields E2 and E1. + + + * --------------- | F=1, 2P_3/2 > + * \ \ + * \ E3_r \ -------- | F=2, 2P_+1/2 > + * \ E4_r \ / \ + * \ \ / E2_l \ + * \ / \ E1_l + * | F=2, 2S_1/2 > -------------- \ + * \ \ + * \ \ + * ------------- | F=1, 2S_1/2 > + * + + + We are solving + dE/dz+(1/c)*dE/dt=i*eta*rho_ij, where j level is higher than i. + Note that E is actually a Rabi frequency of electromagnetic field not the EM field + in xmds terms it looks like + dE_dz = i*eta*rhoij - 1/c*L[E], here we moved t dependence to Fourier space + + VERY IMPORTANT: all Rabi frequency should be given in [1/s], if you want to + normalize it to something else look drho/dt equation. + No need to renormalizes eta as long as its express through + the upper level decay rate in the same units as Rabi frequency. + + + + + section + double v_thermal_averaged=0; + // Maxwell distribution velocities range to take in account in [m/s] + double V_maxwell_min = 0, V_maxwell_max = 0; + + // repopulation rate (atoms flying in/out the laser beam) in [1/s] + const double gt=0.01 *(2*M_PI*1e6); + + // Larmor frequency + double WL=0; + + + // inner use variables + double probability_v; // will be used as p(v) in Maxwell distribution + + ]]> + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + 0 to provide range for Maxwell velocity distribution\n"); + v_thermal_averaged=sqrt(k_boltzmann*Temperature/mass); + // Maxwell distribution velocities range to take in account in [m/s] + // there is almost zero probability for higher velocity p(4*v_av) = 3.3e-04 * p(0) + V_maxwell_min = -4*v_thermal_averaged; V_maxwell_max = -V_maxwell_min; + ]]> + + + + + + + + + + + + + + z + + + + + + + + + + + Ep1 Em1 Ep2 Em2 Ep3 Em3 Ep4 Em4 + + + + + + + + probability_v + + + + + + + + probability_v_norm + + Maxwell_distribution_probabilities + + + + + + + + Ep1a Em1a Ep2a Em2a Ep3a Em3a Ep4a Em4a + + E_field Maxwell_distribution_probabilities Maxwell_distribution_probabilities_norm + + + + + + + + +r0101 r0113 r0202 r0214 r0303 r0309 r0315 r0404 r0410 r0416 r0505 r0511 r0602 r0606 r0614 r0703 r0707 r0709 r0715 r0804 r0808 r0810 r0816 r0909 r0915 r1010 r1016 r1111 r1313 r1414 r1515 r1616 + + + + + + // this is because if r11=1, everything else is 0 and then every small increment + // seems to be huge and adaptive solver makes smaller and smaller steps. + // As quick and dirty fix I reshuffle initial population + // so some of the population sits at the second ground level |2> + // TODO: Fix above. Make the equation of motion for r11 + // and express other level, let's say r44 + // through population normalization + + //read from Mathematica generated RbInits.txt +//---------------- RbInits.txt starts ------------------ +r0101 = 0.125; +r0109 = 0; +r0110 = 0; +r0111 = 0; +r0112 = 0; +r0113 = 0; +r0114 = 0; +r0115 = 0; +r0116 = 0; +r0201 = 0; +r0202 = 0.125; +r0209 = 0; +r0210 = 0; +r0211 = 0; +r0212 = 0; +r0213 = 0; +r0214 = 0; +r0215 = 0; +r0216 = 0; +r0301 = 0; +r0302 = 0; +r0303 = 0.125; +r0309 = 0; +r0310 = 0; +r0311 = 0; +r0312 = 0; +r0313 = 0; +r0314 = 0; +r0315 = 0; +r0316 = 0; +r0401 = 0; +r0402 = 0; +r0403 = 0; +r0404 = 0.125; +r0409 = 0; +r0410 = 0; +r0411 = 0; +r0412 = 0; +r0413 = 0; +r0414 = 0; +r0415 = 0; +r0416 = 0; +r0501 = 0; +r0502 = 0; +r0503 = 0; +r0504 = 0; +r0505 = 0.125; +r0509 = 0; +r0510 = 0; +r0511 = 0; +r0512 = 0; +r0513 = 0; +r0514 = 0; +r0515 = 0; +r0516 = 0; +r0601 = 0; +r0602 = 0; +r0603 = 0; +r0604 = 0; +r0605 = 0; +r0606 = 0.125; +r0609 = 0; +r0610 = 0; +r0611 = 0; +r0612 = 0; +r0613 = 0; +r0614 = 0; +r0615 = 0; +r0616 = 0; +r0701 = 0; +r0702 = 0; +r0703 = 0; +r0704 = 0; +r0705 = 0; +r0706 = 0; +r0707 = 0.125; +r0709 = 0; +r0710 = 0; +r0711 = 0; +r0712 = 0; +r0713 = 0; +r0714 = 0; +r0715 = 0; +r0716 = 0; +r0801 = 0; +r0802 = 0; +r0803 = 0; +r0804 = 0; +r0805 = 0; +r0806 = 0; +r0807 = 0; +r0808 = 0.125; +r0809 = 0; +r0810 = 0; +r0811 = 0; +r0812 = 0; +r0813 = 0; +r0814 = 0; +r0815 = 0; +r0816 = 0; +r0909 = 0; +r0912 = 0; +r0913 = 0; +r0914 = 0; +r0915 = 0; +r0916 = 0; +r1009 = 0; +r1010 = 0; +r1012 = 0; +r1013 = 0; +r1014 = 0; +r1015 = 0; +r1016 = 0; +r1109 = 0; +r1110 = 0; +r1111 = 0; +r1112 = 0; +r1113 = 0; +r1114 = 0; +r1115 = 0; +r1116 = 0; +r1212 = 0; +r1312 = 0; +r1313 = 0; +r1412 = 0; +r1413 = 0; +r1414 = 0; +r1512 = 0; +r1513 = 0; +r1514 = 0; +r1515 = 0; +r1612 = 0; +r1613 = 0; +r1614 = 0; +r1615 = 0; +r1616 = 0; +//---------------- RbInits.txt ends ------------------ + ]]> + + + + + + + + + + + + + + + + + 100 + + + + + + density_matrix + E_field_avgd + + + + + + + + + Lt + + + E_field + density_matrix + + + + + + + + + + + + E_field_avgd + Ip1_out Im1_out Ip2_out Im2_out Ip3_out Im3_out Ip4_out Im4_out + + + + + + + + -- cgit v1.2.3