diff options
Diffstat (limited to 'xmds2')
| -rw-r--r-- | xmds2/realistic_Rb/realistic_Rb.xmds | 26 |
1 files changed, 8 insertions, 18 deletions
diff --git a/xmds2/realistic_Rb/realistic_Rb.xmds b/xmds2/realistic_Rb/realistic_Rb.xmds index 1b11fd4..1111b33 100644 --- a/xmds2/realistic_Rb/realistic_Rb.xmds +++ b/xmds2/realistic_Rb/realistic_Rb.xmds @@ -61,8 +61,10 @@ const double ha2 = 5.323020344462938e8; const double hb2 = 7.85178251911697e7; const double g2 = 3.8117309832741246e7; - const double lambda1 = 0.00007949788511562656; - const double lambda2 = 0.00007802412096860508; + const double lambda1 = 7.949788511562656e-7; + const double lambda2 = 7.802412096860509e-7; + const double eta1 = 5.450949336831401e-6; + const double eta2 = 5.5397657647874e-6; const double rt6 = 2.449489742783178; const double rt3 = 1.7320508075688772; const double rt2 = 1.4142135623730951; @@ -75,11 +77,6 @@ // Simplified k-vectors const double Kvec1 = Kvec, Kvec2=Kvec, Kvec3=Kvec; - const double Gamma_super=6*(2*M_PI*1e6); // characteristic decay rate of upper level used for eta calculations expressed in [1/s] - // eta will be calculated in the <arguments> section - double eta = 0; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] - double eta1=0, eta2=0, eta3=0; - // --------- Atom and cell properties ------------------------- // range of Maxwell distribution atomic velocities const double mass = (86.909180527 * 1.660538921e-27); // atom mass in [kg] @@ -137,13 +134,6 @@ // Maxwell distribution velocities range to take in account in [m/s] // there is almost zero probability for higher velocity p(4*v_av) = 3.3e-04 * p(0) V_maxwell_min = -4*v_thermal_averaged; V_maxwell_max = -V_maxwell_min; - - // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] - eta = 3*lambda*lambda*Ndens*Gamma_super/8.0/M_PI; - // !FIXME over simplification: we should use relevant levels linewidths - eta1 = eta; - eta2 = eta; - eta3 = eta; ]]> </arguments> <bing /> @@ -404,10 +394,10 @@ <dependencies>density_matrix</dependencies> <![CDATA[ //read from Mathematica generated RbPropEquations.txt - dE1_dz = 0.16666666666666666*i*(2.449489742783178*conj(r0309) + 4.242640687119286*conj(r0410) + 6.*conj(r0511))*eta1 - Lt[E1]; - dE2_dz = -0.4082482904638631*i*(conj(r0709) + conj(r0810))*eta1 - Lt[E2]; - dE3_dz = -0.3333333333333333*i*(1.7320508075688772*conj(r0614) + 3.*conj(r0715) + 4.242640687119286*conj(r0816))*eta2 - Lt[E3]; - dE4_dz = (i*(2.449489742783178*conj(r0113) + 3*conj(r0214) + 3*conj(r0315) + 2.449489742783178*conj(r0416))*eta2)/3. - Lt[E4]; + dE1_dz = 0.16666666666666666*i*Ndens*(2.449489742783178*conj(r0309) + 4.242640687119286*conj(r0410) + 6.*conj(r0511))*eta1 - Lt[E1]; + dE2_dz = -0.4082482904638631*i*Ndens*(conj(r0709) + conj(r0810))*eta1 - Lt[E2]; + dE3_dz = -0.3333333333333333*i*Ndens*(1.7320508075688772*conj(r0614) + 3.*conj(r0715) + 4.242640687119286*conj(r0816))*eta2 - Lt[E3]; + dE4_dz = (i*Ndens*(2.449489742783178*conj(r0113) + 3*conj(r0214) + 3*conj(r0315) + 2.449489742783178*conj(r0416))*eta2)/3. - Lt[E4]; ]]> </operators> </integrate> |