diff options
| -rw-r--r-- | MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds | 30 |
1 files changed, 15 insertions, 15 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 2fe5558..65e5922 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 @@ -18,12 +18,12 @@ * / / \ \ * EdL / / \ \ EdR * / / \ \ - * |c> ----------/-- --\----------- |B> + * |c> ----------/-- --\----------- |C> * / \ * / EpL \ EpR * / \ * / \ - * |C> ------------- ------------- |b> + * |B> ------------- ------------- |b> * * * @@ -304,14 +304,14 @@ double delta_two_ph = delta_two_ph0; delta_two_ph += abs(Gamma_super)*1*((t)/(t_span/2)); // two photon sweep to mimic EIT shapes - 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 + Gab=Gamma_super+i*(delta_one_ph + delta_two_ph + DopplerShiftPr); // probe right detuning + GAB=Gamma_super+i*(delta_one_ph + delta_two_ph + DopplerShiftPr); // 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 + GCA=Gamma_super-i*(delta_one_ph + DopplerShiftDr); // drive right detuning - Gcb=gbc+i*( Split + delta_two_ph + (DopplerShiftPr - DopplerShiftDr) + 0*noise); - GCB=gbc+i*(-Split - delta_two_ph - (DopplerShiftPr - DopplerShiftDr) - 0*noise); + Gcb=gbc+i*( Split + delta_two_ph + (DopplerShiftPr - DopplerShiftDr) ); + GCB=gbc+i*(-Split + delta_two_ph + (DopplerShiftPr - DopplerShiftDr) ); complex rba=conj(rab); complex rac=conj(rca); @@ -326,21 +326,21 @@ complex EpRac = conj(EpRa); 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; + -i*EpLac*rAB+i*EpLa*rBA-i*EdRac*rAC+i*EdRa*rCA; 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; + drBB_dt = i*EpLac*rAB-i*EpLa*rBA+Gamma_super*raa-gbc*rBB+gbc*rCC; 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; + drCC_dt = i*EdRac*rAC-i*EdRa*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; drcb_dt = -Gcb*rcb-i*EpRa*rca+i*EdLac*rab; - drAB_dt = -Gab*rAB+i*EdRa*(rBB-raa)+i*EpLa*rCB; - drCA_dt = -Gca*rCA+i*EpLac*(raa-rCC)-i*EdRac*rCB; - drCB_dt = -GCB*rCB-i*EdRa*rCA+i*EpLac*rAB; + drAB_dt = -Gab*rAB+i*EpLa*(rBB-raa)+i*EdRa*rCB; + drCA_dt = -Gca*rCA+i*EdRac*(raa-rCC)-i*EpLac*rCB; + drCB_dt = -GCB*rCB-i*EpLa*rCA+i*EdRac*rAB; ]]> @@ -361,9 +361,9 @@ <dependencies>density_matrix</dependencies> <![CDATA[ dEdL_dz = i*eta*conj(rca) + Lt[EdL] ; - dEdR_dz = i*eta*(rAB) + Lt[EdR] ; - dEpL_dz = i*eta*conj(rCA) + Lt[EpL] ; + dEdR_dz = i*eta*conj(rCA) + Lt[EdR] ; dEpR_dz = i*eta*(rab) + Lt[EpR] ; + dEpL_dz = i*eta*(rAB) + Lt[EpL] ; ]]> </operators> |