diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-21 23:07:15 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-21 23:07:15 +0000 |
| commit | a85028715e3a1460f9f6697592e62b07341b4f9f (patch) | |
| tree | 3b61bfb70a10752191bfa5bdfb3cfc423145553a | |
| parent | c680b440e6963ab1f507ca0fcefdc3ce982db8af (diff) | |
| download | multi_mode_eit-a85028715e3a1460f9f6697592e62b07341b4f9f.tar.gz multi_mode_eit-a85028715e3a1460f9f6697592e62b07341b4f9f.zip | |
ground levels decay improved
| -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 |