diff options
| -rw-r--r-- | liouville.m | 4 |
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)) |