diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-21 19:36:02 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-21 19:36:02 +0000 |
| commit | d30c969a69096e0244800782d4636e98bb9bcf55 (patch) | |
| tree | 10501904e588b8c8e0447f2a5066466d47502eb8 /rb87_D1_line.m | |
| parent | 4ba3b5b9c045fda822ceecf227de98ecd73ef28c (diff) | |
| download | multi_mode_eit-d30c969a69096e0244800782d4636e98bb9bcf55.tar.gz multi_mode_eit-d30c969a69096e0244800782d4636e98bb9bcf55.zip | |
proper dipole elements for Rb87D1 line
Diffstat (limited to 'rb87_D1_line.m')
| -rw-r--r-- | rb87_D1_line.m | 25 |
1 files changed, 10 insertions, 15 deletions
diff --git a/rb87_D1_line.m b/rb87_D1_line.m index 1739df2..ea20e11 100644 --- a/rb87_D1_line.m +++ b/rb87_D1_line.m @@ -82,22 +82,15 @@ for j=1:Nlevels for k=1:Nlevels if ( abs(ang_momentum(j) - ang_momentum(k)) == 1) %transition allowed for L =L' +/- 1 - % incorrect Clebsch-Gordan coefficients but should show some ideas % but they correct within a factor, but not sign - if ( ((m(j)- m(k))==1) && ( H0(j,j) < H0(k,k)) ) - dipole_elements.left(j,k)=1; - endif - if ( ((m(j)- m(k))==-1) && ( H0(j,j) < H0(k,k)) ) - dipole_elements.right(j,k)=1; - endif - if ( ((m(j)- m(k))==0) && ( H0(j,j) < H0(k,k)) ) - 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 - endif + if ( ( H0(j,j) < H0(k,k)) ) + de=dipole_elementRb87D1line(... + total_momentum(j),m(j), ... + total_momentum(k),m(k) ... + ); + dipole_elements.left(j,k) = de.left; + dipole_elements.right(j,k) = de.right; + dipole_elements.linear(j,k) = de.linear; endif endif endfor @@ -141,3 +134,5 @@ endfor energy+=B_field * ( gm .* m); % convert frequency to energy units H0=diag(energy)*hbar; + +% vim: ts=2:sw=2:fdm=indent |