From 8aa7f5a5ec78704fa9a86c0cddbda53f36457e03 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Wed, 9 Dec 2009 06:01:39 +0000 Subject: decay part moved to the end result --- liouville.m | 4 +--- 1 file changed, 1 insertion(+), 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)) -- cgit v1.2.3