diff options
Diffstat (limited to 'liouville.m')
| -rw-r--r-- | liouville.m | 14 |
1 files changed, 8 insertions, 6 deletions
diff --git a/liouville.m b/liouville.m index 4c78557..4ce714c 100644 --- a/liouville.m +++ b/liouville.m @@ -1,4 +1,5 @@ 1; +clear all; t0 = clock (); % we will use this latter to calculate elapsed time @@ -30,14 +31,15 @@ kappa_p =zeros(1,N_detun_steps+1); kappa_m =zeros(1,N_detun_steps+1); detun_step=(detuning_p_max-detuning_p_min)/N_detun_steps; +% now we create Liouville indexes list +[N, rhoLiouville_w, rhoLiouville_r, rhoLiouville_c]=unfold_density_matrix(Nlevels,Nfreq); +rhoLiouville=zeros(N,1); % calculate E_field independent properties of athe atom % to be used as sub matrix templates for Liouville operator matrix [L0m, polarizability_m]=L0_and_polarization_submatrices( ... - Nlevels*Nlevels, ... - H0, g_decay, g_dephasing, dipole_elements, ... - E_field, ... - modulation_freq, rhoLiouville_w, rhoLiouville_r, rhoLiouville_c ... + Nlevels, ... + H0, g_decay, g_dephasing, dipole_elements ... ); for detuning_p_cntr=1:N_detun_steps+1; wp0=w12; @@ -77,8 +79,8 @@ for detuning_p_cntr=1:N_detun_steps+1; %rho_d=rhoOfFreq(rhoLiouville, 3, Nlevels, Nfreq); % drive frequency %rho_m=rhoOfFreq(rhoLiouville, 4, Nlevels, Nfreq); % opposite sideband frequency - kappa_p(detuning_p_cntr)=sucseptibility(2, rhoLiouville, dipole_elements, Nlevels, Nfreq); - %kappa_m(detuning_p_cntr)=sucseptibility(4, rhoLiouville, dipole_elements, Nlevels, Nfreq); + kappa_p(detuning_p_cntr)=susceptibility(2, rhoLiouville, dipole_elements, Nlevels, Nfreq); + %kappa_m(detuning_p_cntr)=susceptibility(4, rhoLiouville, dipole_elements, Nlevels, Nfreq); detuning_freq(detuning_p_cntr)=detuning_p; %kappa_p_re=real(kappa_p); |