diff options
| -rw-r--r-- | xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds | 10 | ||||
| -rw-r--r-- | xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds | 10 |
2 files changed, 16 insertions, 4 deletions
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds b/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds index 5102c46..045adac 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds @@ -36,7 +36,7 @@ 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 i + No need to renormalizes eta as long as its express through the upper level decay rate in the same units as Rabi frequency. </description> @@ -49,7 +49,8 @@ const double lambda=794.7e-9; //wavelength in m const double Kvec = 2*M_PI/lambda; // k-vector const double Gamma_super=6*(2*M_PI*1e6); // characteristic decay rate of upper level used for eta calculations expressed in [1/s] - const double eta = 3*lambda*lambda*Ndens*Gamma_super/8.0/M_PI; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] + // eta will be calculated in the sequence -> filter element to use --Ndens parameter + double eta = 0; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] // --------- Atom and cell properties ------------------------- // range of Maxwell distribution atomic velocities @@ -218,6 +219,11 @@ </vector> <sequence> + <filter> + <![CDATA[ + eta = 3*lambda*lambda*Ndens*Gamma_super/8.0/M_PI; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] + ]]> + </filter> <!--For this set of conditions ARK45 is faster than ARK89--> <!--ARK45 is good for small detuning when all frequency like term are close to zero--> <integrate algorithm="ARK45" tolerance="1e-5" interval="Lcell"> diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds b/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds index e0440ef..d662779 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds @@ -41,7 +41,7 @@ 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 i + No need to renormalizes eta as long as its express through the upper level decay rate in the same units as Rabi frequency. </description> @@ -54,7 +54,8 @@ const double lambda=794.7e-9; //wavelength in m const double Kvec = 2*M_PI/lambda; // k-vector const double Gamma_super=6*(2*M_PI*1e6); // characteristic decay rate of upper level used for eta calculations expressed in [1/s] - const double eta = 3*lambda*lambda*Ndens*Gamma_super/8.0/M_PI; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] + // eta will be calculated in the sequence -> filter element to use --Ndens parameter + double eta = 0; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] // --------- Atom and cell properties ------------------------- // range of Maxwell distribution atomic velocities @@ -225,6 +226,11 @@ </vector> <sequence> + <filter> + <![CDATA[ + eta = 3*lambda*lambda*Ndens*Gamma_super/8.0/M_PI; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] + ]]> + </filter> <!--For this set of conditions ARK45 is faster than ARK89--> <!--ARK45 is good for small detuning when all frequency like term are close to zero--> <integrate algorithm="ARK45" tolerance="1e-5" interval="Lcell"> |