diff options
Diffstat (limited to 'xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations')
| -rw-r--r-- | xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds | 30 |
1 files changed, 25 insertions, 5 deletions
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 e7e7b4b..e0440ef 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 @@ -50,10 +50,21 @@ <![CDATA[ const double pi = M_PI; const double c=3.e8; + const double k_boltzmann= 1.3806505e-23; // Boltzmann knostant in [J/K] const double lambda=794.7e-9; //wavelength in m - const double N=1e09*(1e6); //number of particles per cubic m i.e. density + 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*N*Gamma_super/8.0/M_PI; // eta constant in the wave equation for Rabi frequency. Units are [1/(m 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)] + + // --------- Atom and cell properties ------------------------- + // range of Maxwell distribution atomic velocities + const double mass = (86.909180527 * 1.660538921e-27); // atom mass in [kg] + // above mass expression is written as (expression is isotopic_mass * atomic_mass_unit) + + // Average sqrt(v^2) in Maxwell distribution for one dimension + const double v_thermal_averaged=sqrt(k_boltzmann*Temperature/mass); + // Maxwell distribution velocities range to take in account in [m/s] + const double V_maxwell_min = -4*v_thermal_averaged, V_maxwell_max = -V_maxwell_min; // there is almost zero probability for higher velocity p(4*v_av) = 3.3e-04 * p(0) // repopulation rate (atoms flying in/out the laser beam) in [1/s] const double gt=0.01 *(2*M_PI*1e6); @@ -72,6 +83,7 @@ complex r21, r31, r41, r32, r42, r43; // density matrix elements ]]> </globals> + <validation kind="run-time"/> <!--allows to put ranges as variables--> <benchmark /> <arguments> <!-- Rabi frequency divided by 2 in [1/s] --> @@ -83,10 +95,18 @@ <argument name="delta1" type="real" default_value="0.0" /> <argument name="delta2" type="real" default_value="0.0" /> <argument name="delta3" type="real" default_value="0.0" /> - <!--Cell length--> - <argument name="Lcell" type="real" default_value="1.5e-2" /> - <!--Pulse duration/width--> + <!--Pulse duration/width [s] --> <argument name="Pwidth" type="real" default_value="0.1e-6" /> + <!-- Atom and cell properties --> + <!--Cell length [m] --> + <argument name="Lcell" type="real" default_value="1.5e-2" /> + <!--Density of atoms [1/m^3] --> + <argument name="Ndens" type="real" default_value="1e15" /> + <!--Atoms temperature [K] --> + <!--TODO: looks like Temperature > 10 K knocks solver, + I am guessing detunings are too large and thus it became a stiff equation--> + <!--! make sure it is not equal to zero!--> + <argument name="Temperature" type="real" default_value="5" /> </arguments> <bing /> <fftw plan="patient" /> |