diff options
| -rw-r--r-- | rb87_D1_line.m | 15 |
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 |