diff options
| -rw-r--r-- | xmds2/Nlevels_no_dopler_with_z/Nlevels_no_dopler_with_z.xmds | 17 |
1 files changed, 9 insertions, 8 deletions
diff --git a/xmds2/Nlevels_no_dopler_with_z/Nlevels_no_dopler_with_z.xmds b/xmds2/Nlevels_no_dopler_with_z/Nlevels_no_dopler_with_z.xmds index ea12b3f..ccde248 100644 --- a/xmds2/Nlevels_no_dopler_with_z/Nlevels_no_dopler_with_z.xmds +++ b/xmds2/Nlevels_no_dopler_with_z/Nlevels_no_dopler_with_z.xmds @@ -30,12 +30,12 @@ * We moved to dimensionless units - t -> t*g ,time - z -> z*g/c , distance - rabi_frequency -> rabi_frequency/g - eta -> eta*c/g^2 , coupling constant - gij -> gij/g - Wij -> Wij/g + t -> t*DecayRateNormalization ,time + z -> z*DecayRateNormalization/c , distance + rabi_frequency -> rabi_frequency/DecayRateNormalization + eta -> eta*c/DecayRateNormalization^2 , coupling constant + gij -> gij/DecayRateNormalization + Wij -> Wij/DecayRateNormalization where g is 1MHz rate @@ -63,11 +63,11 @@ <auto_vectorise /> <globals> <![CDATA[ - const double UpperDecayRate=10e6; // measured in s^-1 + const double DecayRateNormalization=10e6; // measured in s^-1 const double c=3e8; const double lambda=794.7e-9; //wavelength in m const double N=1e9; //number of particles per cubic cm - const double eta = 3*lambda*lambda*(N/1e-6)*UpperDecayRate/8.0/3.14*(c/UpperDecayRate/UpperDecayRate); //dimensionless coupling constant + const double eta = 3*lambda*lambda*(N/1e-6)*DecayRateNormalization/8.0/3.14*(c/DecayRateNormalization/DecayRateNormalization); //dimensionless coupling constant // repopulation rate (atoms flying in/out the laser beam) in MHz const double gt=0.01/2; @@ -94,6 +94,7 @@ </globals> </features> + <!-- Let's convert 'z' and 't' to dimentionless units over here --> <geometry> <propagation_dimension> z </propagation_dimension> <transverse_dimensions> |