diff options
| -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 |