summaryrefslogtreecommitdiff
path: root/liouville.m
diff options
context:
space:
mode:
authorEugeniy Mikhailov <evgmik@gmail.com>2009-12-09 06:01:39 +0000
committerEugeniy Mikhailov <evgmik@gmail.com>2009-12-09 06:01:39 +0000
commit8aa7f5a5ec78704fa9a86c0cddbda53f36457e03 (patch)
treea698b5858d629ca1c8d929e21093a064f32256c7 /liouville.m
parent35af744d95df0d687bccdaaa4dabda6fe830bed7 (diff)
downloadmulti_mode_eit-8aa7f5a5ec78704fa9a86c0cddbda53f36457e03.tar.gz
multi_mode_eit-8aa7f5a5ec78704fa9a86c0cddbda53f36457e03.zip
decay part moved to the end result
Diffstat (limited to 'liouville.m')
-rw-r--r--liouville.m4
1 files changed, 1 insertions, 3 deletions
diff --git a/liouville.m b/liouville.m
index 67dde6b..1dc12b0 100644
--- a/liouville.m
+++ b/liouville.m
@@ -75,8 +75,6 @@ for p=1:N
decay_part=\
( decay_total(g_decay,k)/2 + decay_total(g_decay,j)/2 + g_dephasing(j,k) )* kron_delta(j,m)*kron_delta(k,n) \
- kron_delta(m,n)*kron_delta(j,k)*g_decay(m,j) ;
- decay_part=decay_part*hbar/im_one;
- L0+=decay_part;
Lt=L0;
else
Li= ( dipole_elements(j,m)*kron_delta(k,n)-dipole_elements(n,k)*kron_delta(j,m) )*E_field(w3i);
@@ -84,7 +82,7 @@ for p=1:N
endif
%Lt=-im_one/hbar*Lt*kron_delta(w_jk-w_iner,w_mn); % above if should be done only if kron_delta is not zero
% no need for above kron_delta since the same conditon checked in the outer if statement
- Lt=-im_one/hbar*Lt;
+ Lt=-im_one/hbar*Lt - decay_part;
endif
endfor
if ((p == s))