diff options
| -rw-r--r-- | MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds | 14 |
1 files changed, 10 insertions, 4 deletions
diff --git a/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds b/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds index 8c11b5a..001b424 100644 --- a/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds +++ b/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds @@ -328,14 +328,20 @@ complex EpLac = conj(EpLa); complex EpRac = conj(EpRa); + // ground level repopulation rate, i.e. arrival from other levels decay + // note that a level repopulate itself, thus there is 4 in + // a term like drbb_dt = ... + Repop - 4*gbc*rbb; + // which otherwise decays only in 3 other ground levels + complex Repop = Gamma_super*raa + gbc*(rbb+rcc+rBB+rCC); + draa_dt = -i*EpRac*rab+i*EpRa*rba-i*EdLac*rac+i*EdLa*rca-4*Gamma_super*raa -i*EpLac*rAB+i*EpLa*rBA-i*EdRac*rAC+i*EdRa*rCA; - drbb_dt = i*EpRac*rab-i*EpRa*rba+Gamma_super*raa-gbc*rbb+gbc*rcc; - drBB_dt = i*EpLac*rAB-i*EpLa*rBA+Gamma_super*raa-gbc*rBB+gbc*rCC; + drbb_dt = i*EpRac*rab-i*EpRa*rba + Repop - 4*gbc*rbb; + drBB_dt = i*EpLac*rAB-i*EpLa*rBA + Repop - 4*gbc*rBB; - drcc_dt = i*EdLac*rac-i*EdLa*rca+Gamma_super*raa-gbc*rcc+gbc*rbb; - drCC_dt = i*EdRac*rAC-i*EdRa*rCA+Gamma_super*raa-gbc*rCC+gbc*rBB; + drcc_dt = i*EdLac*rac-i*EdLa*rca + Repop - 4*gbc*rcc; + drCC_dt = i*EdRac*rAC-i*EdRa*rCA + Repop - 4*gbc*rCC; drab_dt = -Gab*rab+i*EpRa*(rbb-raa)+i*EdLa*rcb; drca_dt = -Gca*rca+i*EdLac*(raa-rcc)-i*EpRac*rcb; |