diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-21 18:24:42 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-21 18:24:42 +0000 |
| commit | 4ba3b5b9c045fda822ceecf227de98ecd73ef28c (patch) | |
| tree | 64f41e480f2786946f68be99eab1d1e9692e53f6 /rb87_D1_line.m | |
| parent | 23940f2c21ce176c032944db4c330c0f487c7f19 (diff) | |
| download | multi_mode_eit-4ba3b5b9c045fda822ceecf227de98ecd73ef28c.tar.gz multi_mode_eit-4ba3b5b9c045fda822ceecf227de98ecd73ef28c.zip | |
dipole elements for j->k transition where j is lower energy level
Diffstat (limited to 'rb87_D1_line.m')
| -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 78b8ffa..1739df2 100644 --- a/rb87_D1_line.m +++ b/rb87_D1_line.m @@ -77,19 +77,20 @@ H0=diag(energy)*hbar; dipole_elements.left = zeros(Nlevels); dipole_elements.right = zeros(Nlevels); dipole_elements.linear = zeros(Nlevels); +% define dipole elements for transition from level j->k 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)) ) + 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)) ) + 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)) ) + 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 |