summaryrefslogtreecommitdiff
path: root/rb87_D1_line.m
diff options
context:
space:
mode:
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