diff options
Diffstat (limited to 'xmds2')
| -rw-r--r-- | xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds | 21 |
1 files changed, 13 insertions, 8 deletions
diff --git a/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds b/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds index 1e5b25c..d63704f 100644 --- a/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds +++ b/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds @@ -49,14 +49,11 @@ const double k_boltzmann= 1.3806505e-23; // Boltzmann knostant in [J/K] const double lambda=794.7e-9; //wavelength in m const double Kvec = 2*M_PI/lambda; // k-vector - const double N=1e09*(1e6); //number of particles per cubic m i.e. density 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 Temperature=5; // cell temperature in [K] ! make sure it is not equal to zero! - // TODO: looks like Temperature > 10 K knocks solver, I am guessing detunings are too large and thus it became a stiff equation 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) @@ -85,7 +82,7 @@ ]]> </globals> - <validation kind="run-time"/> <!--allows to put ranges as wariables--> + <validation kind="run-time"/> <!--allows to put ranges as variables--> <benchmark /> <arguments> <!-- Rabi frequency divided by 2 in [1/s] --> @@ -97,10 +94,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" /> |