diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2016-05-10 11:49:47 -0400 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2016-05-10 11:49:47 -0400 |
| commit | 5957b1b38a3bc309afaca0612e65935f98f412ff (patch) | |
| tree | 5d3a34b2390ccdb1b440aa6f8ffe3460943c059d | |
| parent | 4ec6e1666ea62a6f906e64fc9a297bd0e918355d (diff) | |
| download | noisy_eit_xmds-5957b1b38a3bc309afaca0612e65935f98f412ff.tar.gz noisy_eit_xmds-5957b1b38a3bc309afaca0612e65935f98f412ff.zip | |
get rid of unused parameters, added 2photon detuning scan, doppler shift readiness code added
| -rw-r--r-- | MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds | 31 |
1 files changed, 17 insertions, 14 deletions
diff --git a/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds b/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds index d486176..dfe2e71 100644 --- a/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds +++ b/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds @@ -71,8 +71,7 @@ //const double Ga=3.0 *(2*M_PI*1e6); //const double G4=3.0 *(2*M_PI*1e6); - complex g = Gamma_super; - complex gbc = (2*M_PI)*1e3; // 1/s untits + const double gbc = (2*M_PI)*1e3; // 1/s untits const complex Split = 0; const complex noise = 0; @@ -84,6 +83,8 @@ // inner use variables double probability_v; // will be used as p(v) in Maxwell distribution + // --------- Simulation parameters ------------------------- + const double t_span = 2e-5; // span of the simulation time lattice ]]> </globals> <validation kind="run-time"/> <!--allows to put ranges as variables--> @@ -98,8 +99,8 @@ <!-- Fields detuning in [1/s] --> - <argument name="delta_one_ph" type="real" default_value="0.0" /> - <argument name="delta_two_ph" type="real" default_value="0.0" /> + <argument name="delta_one_ph0" type="real" default_value="0.0" /> + <argument name="delta_two_ph0" type="real" default_value="0.0" /> <!--Pulse duration/width [s] --> <argument name="Pwidth" type="real" default_value="0.1e-6" /> <!-- Atom and cell properties --> @@ -152,7 +153,7 @@ * lattice="10000" domain="(-1e-6, 1e-6)" works for Doppler averaging in up to 400K for Rb when lasers are zero detuned --> - <dimension name="t" lattice="10000" domain="(-1e-6, 1e-6)" /> + <dimension name="t" lattice="10000" domain="(-t_span/2, t_span/2)" /> <dimension name="v" lattice="3" domain="(V_maxwell_min, V_maxwell_max)" /> </transverse_dimensions> </geometry> @@ -295,12 +296,14 @@ double DopplerShiftDr = KvecDr * v; double DopplerShiftPr = KvecDr * v; + double delta_one_ph = delta_one_ph0; + double delta_two_ph = delta_two_ph0 + abs(Gamma_super)*10*((t)/(t_span/2)); - Gab=g+i*(delta_one_ph + delta_two_ph + DopplerShiftPr + 0*noise); // probe right detuning - GAB=g+i*(delta_one_ph + DopplerShiftDr + 0*noise); // drive right detuning + Gab=Gamma_super+i*(delta_one_ph + delta_two_ph + DopplerShiftPr + 0*noise); // probe right detuning + GAB=Gamma_super+i*(delta_one_ph + DopplerShiftDr + 0*noise); // drive right detuning - Gca=g-i*(delta_one_ph + DopplerShiftDr); // drive left detuning - GCA=g-i*(delta_one_ph + delta_two_ph + DopplerShiftDr); // probe left detuning + Gca=Gamma_super-i*(delta_one_ph + DopplerShiftDr); // drive left detuning + GCA=Gamma_super-i*(delta_one_ph + delta_two_ph + DopplerShiftDr); // probe left detuning Gcb=gbc+i*( Split + delta_two_ph + (DopplerShiftPr - DopplerShiftDr) + 0*noise); GCB=gbc+i*(-Split - delta_two_ph - (DopplerShiftPr - DopplerShiftDr) - 0*noise); @@ -317,14 +320,14 @@ complex EpLac = conj(EpLa); complex EpRac = conj(EpRa); - draa_dt = -i*EpRac*rab+i*EpRa*rba-i*EdLac*rac+i*EdLa*rca-2*g*raa + draa_dt = -i*EpRac*rab+i*EpRa*rba-i*EdLac*rac+i*EdLa*rca-2*Gamma_super*raa -i*EdRac*rAB+i*EdRa*rBA-i*EpLac*rAC+i*EpLa*rCA; - drbb_dt = i*EpRac*rab-i*EpRa*rba+g*raa-gbc*rbb+gbc*rcc; - drBB_dt = i*EdRac*rAB-i*EdRa*rBA+g*raa-gbc*rBB+gbc*rCC; + drbb_dt = i*EpRac*rab-i*EpRa*rba+Gamma_super*raa-gbc*rbb+gbc*rcc; + drBB_dt = i*EdRac*rAB-i*EdRa*rBA+Gamma_super*raa-gbc*rBB+gbc*rCC; - drcc_dt = i*EdLac*rac-i*EdLa*rca+g*raa-gbc*rcc+gbc*rbb; - drCC_dt = i*EpLac*rAC-i*EpLa*rCA+g*raa-gbc*rCC+gbc*rBB; + drcc_dt = i*EdLac*rac-i*EdLa*rca+Gamma_super*raa-gbc*rcc+gbc*rbb; + drCC_dt = i*EpLac*rAC-i*EpLa*rCA+Gamma_super*raa-gbc*rCC+gbc*rBB; drab_dt = -Gab*rab+i*EpRa*(rbb-raa)+i*EdLa*rcb; drca_dt = -Gca*rca+i*EdLac*(raa-rcc)-i*EpRac*rcb; |