From 5957b1b38a3bc309afaca0612e65935f98f412ff Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Tue, 10 May 2016 11:49:47 -0400 Subject: get rid of unused parameters, added 2photon detuning scan, doppler shift readiness code added --- .../Nlevels_with_MOR.xmds | 31 ++++++++++++---------- 1 file 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 ]]> @@ -98,8 +99,8 @@ - - + + @@ -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 --> - + @@ -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; -- cgit v1.2.3