diff options
| -rw-r--r-- | MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds | 17 |
1 files changed, 10 insertions, 7 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 c58cb11..8c11b5a 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 @@ -73,7 +73,6 @@ //const double G4=3.0 *(2*M_PI*1e6); const double gbc = (2*M_PI)*1e3; // 1/s untits - const complex Split = 0; const complex noise = 0; complex Gab, GAB, Gca, GCA, Gcb, GCB; @@ -102,6 +101,10 @@ <!-- Fields detuning in [1/s] --> <argument name="delta_one_ph0" type="real" default_value="0.0" /> <argument name="delta_two_ph0" type="real" default_value="0.0" /> + + <!-- Ground levels differential split (Mimics magnetic field) [1/s] --> + <argument name="Split0" type="real" default_value="0.0" /> + <!--Pulse duration/width [s] --> <argument name="Pwidth" type="real" default_value="0.1e-6" /> <!-- Atom and cell properties --> @@ -304,14 +307,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); // probe right detuning - GAB=Gamma_super+i*(delta_one_ph + delta_two_ph + DopplerShiftPr); // probe left detuning + Gab=Gamma_super+i*( Split0 + delta_one_ph + delta_two_ph + DopplerShiftPr); // probe right detuning + GAB=Gamma_super+i*(-Split0 + 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 + DopplerShiftDr); // drive right detuning + Gca=Gamma_super-i*(-Split0 + delta_one_ph + DopplerShiftDr); // drive left detuning + GCA=Gamma_super-i*( Split0 + delta_one_ph + DopplerShiftDr); // drive right detuning - Gcb=gbc+i*( Split + delta_two_ph + (DopplerShiftPr - DopplerShiftDr) ); - GCB=gbc+i*(-Split + delta_two_ph + (DopplerShiftPr - DopplerShiftDr) ); + Gcb=gbc+i*( 2*Split0 + delta_two_ph + (DopplerShiftPr - DopplerShiftDr) ); + GCB=gbc+i*(-2*Split0 + delta_two_ph + (DopplerShiftPr - DopplerShiftDr) ); complex rba=conj(rab); complex rac=conj(rca); |