diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-15 18:59:42 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-15 18:59:42 +0000 |
| commit | d2133f9b0bb3323fe43dd0b6cc94c1f586fe07d0 (patch) | |
| tree | 1f17f939051185e32ed54509a49ed453cc647561 | |
| parent | db34a46f1d99709ee6f9701579fbb6a7ece49ebc (diff) | |
| download | multi_mode_eit-d2133f9b0bb3323fe43dd0b6cc94c1f586fe07d0.tar.gz multi_mode_eit-d2133f9b0bb3323fe43dd0b6cc94c1f586fe07d0.zip | |
selection rule for m=0 -> m'=0 fixed
| -rw-r--r-- | rb87_D1_line.m | 7 |
1 files changed, 4 insertions, 3 deletions
diff --git a/rb87_D1_line.m b/rb87_D1_line.m index ac66de8..78b8ffa 100644 --- a/rb87_D1_line.m +++ b/rb87_D1_line.m @@ -19,7 +19,7 @@ useful_constants; % ---- ---- ---- |S,F=1> % m=-1 m=0 m=1 -B_field=.2; +B_field=.3; w_hpf_ground=6834; w_hpf_exited=817; @@ -90,11 +90,12 @@ for j=1:Nlevels dipole_elements.right(j,k)=1; endif if ( ((m(j)- m(k))==0) && ( H0(j,j) > H0(k,k)) ) - if ( ( total_momentum(j) == 1 ) && (total_momentum(k) == 1 ) ) + dipole_elements.linear(j,k)=1; + if ( (m(j)==0) && ( total_momentum(j) == 1 ) && (total_momentum(k) == 1 ) ) % remember m=0 ->m'=0 is forbidden when F=1 and F'=1 + % so we drop it back to zero dipole_elements.linear(j,k)=0; else - dipole_elements.linear(j,k)=1; endif endif endif |