diff options
| -rw-r--r-- | useful_functions.m | 13 |
1 files changed, 7 insertions, 6 deletions
diff --git a/useful_functions.m b/useful_functions.m index 4b239a5..913b505 100644 --- a/useful_functions.m +++ b/useful_functions.m @@ -66,6 +66,7 @@ function [L0m, polarizability_m]=L0_and_polarization_submatrices( ... % now we create Liouville indexes list [Ndummy, rhoLiouville_w_notused, rhoLiouville_r, rhoLiouville_c]=unfold_density_matrix(Nlevels,1); + kron_delta_m=eye(Nlevels); % note that L0 and decay parts depend only on combination of indexes % jk,mn but repeats itself for every frequency L0m=zeros(rho_size); % (Nlevels^2)x(Nlevels^2) matrix @@ -86,17 +87,17 @@ function [L0m, polarizability_m]=L0_and_polarization_submatrices( ... n=rhoLiouville_c(s); % calculate unperturbed part (Hamiltonian without EM field) - L0m(p,s)=H0(j,m)*kron_delta(k,n)-H0(n,k)*kron_delta(j,m); + L0m(p,s)=H0(j,m)*kron_delta_m(k,n)-H0(n,k)*kron_delta_m(j,m); decay_part_m(p,s)= ... ( ... 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) ; - polarizability_m.linear(p,s)= ( dipole_elements.linear(j,m)*kron_delta(k,n)-dipole_elements.linear(n,k)*kron_delta(j,m) ); - polarizability_m.left(p,s)= ( dipole_elements.left(j,m)*kron_delta(k,n)-dipole_elements.left(n,k)*kron_delta(j,m) ); - polarizability_m.right(p,s)= ( dipole_elements.right(j,m)*kron_delta(k,n)-dipole_elements.right(n,k)*kron_delta(j,m) ); + )* kron_delta_m(j,m)*kron_delta_m(k,n) ... + - kron_delta_m(m,n)*kron_delta_m(j,k)*g_decay(m,j) ; + polarizability_m.linear(p,s)= ( dipole_elements.linear(j,m)*kron_delta_m(k,n)-dipole_elements.linear(n,k)*kron_delta_m(j,m) ); + polarizability_m.left(p,s)= ( dipole_elements.left(j,m)*kron_delta_m(k,n)-dipole_elements.left(n,k)*kron_delta_m(j,m) ); + polarizability_m.right(p,s)= ( dipole_elements.right(j,m)*kron_delta_m(k,n)-dipole_elements.right(n,k)*kron_delta_m(j,m) ); endfor endfor L0m=-im_one/hbar*L0m - decay_part_m; |