summaryrefslogtreecommitdiff
path: root/rb87_D1_line.m
diff options
context:
space:
mode:
authorEugeniy Mikhailov <evgmik@gmail.com>2009-12-21 23:07:15 +0000
committerEugeniy Mikhailov <evgmik@gmail.com>2009-12-21 23:07:15 +0000
commita85028715e3a1460f9f6697592e62b07341b4f9f (patch)
tree3b61bfb70a10752191bfa5bdfb3cfc423145553a /rb87_D1_line.m
parentc680b440e6963ab1f507ca0fcefdc3ce982db8af (diff)
downloadmulti_mode_eit-a85028715e3a1460f9f6697592e62b07341b4f9f.tar.gz
multi_mode_eit-a85028715e3a1460f9f6697592e62b07341b4f9f.zip
ground levels decay improved
Diffstat (limited to 'rb87_D1_line.m')
-rw-r--r--rb87_D1_line.m15
1 files changed, 12 insertions, 3 deletions
diff --git a/rb87_D1_line.m b/rb87_D1_line.m
index ea20e11..3165df1 100644
--- a/rb87_D1_line.m
+++ b/rb87_D1_line.m
@@ -123,9 +123,18 @@ endfor
gamma_hpf=.001;
for i=1:Nlevels
for j=1:Nlevels
- if ( abs( H0(i,i) - H0(j,j)) == w_hpf_ground*hbar )
- % ground hyperfine are mixed together
- g_decay(i,j)=gamma_hpf;
+ % it would be better to introduce a level corresponding to the bath
+ % but this slows calculations
+ % So
+ % ground hyperfine are equally mixed together
+ if ( (abs( H0(i,i) - H0(j,j)) == w_hpf_ground*hbar ) || ...
+ (abs( H0(i,i) - H0(j,j)) == 0 ) )
+ % i and j correspond to 2 ground levels
+ % we cannot decay to itself
+ if ( i != j )
+ % total eight ground levels F=2 and F=1 so we decay in 7 channels
+ g_decay(i,j)=gamma_hpf/7;
+ endif
endif
endfor
endfor