diff options
Diffstat (limited to 'liouville.m')
| -rw-r--r-- | liouville.m | 14 |
1 files changed, 10 insertions, 4 deletions
diff --git a/liouville.m b/liouville.m index fd1c840..fe49313 100644 --- a/liouville.m +++ b/liouville.m @@ -41,6 +41,10 @@ rhoLiouville=zeros(N,1); Nlevels, ... H0, g_decay, g_dephasing, dipole_elements ... ); + +atom_properties_fname='atom_properties.mat'; +save( atom_properties_fname, 'L0m', 'polarizability_m', 'dipole_elements' ) ; + for detuning_p_cntr=1:N_detun_steps+1; wp0=w12; detuning_p=detuning_p_min+detun_step*(detuning_p_cntr-1); @@ -52,9 +56,10 @@ for detuning_p_cntr=1:N_detun_steps+1; E_field =[0, Ep, Ed, Epc, Edc, 0, 0 ]; freq_index=freq2index(wp,modulation_freq); - atom_field_problem.L0m = L0m; - atom_field_problem.polarizability_m = polarizability_m; - atom_field_problem.dipole_elements = dipole_elements; + atom_field_problem.atom_properties_fname = atom_properties_fname; + %atom_field_problem.L0m = L0m; + %atom_field_problem.polarizability_m = polarizability_m; + %atom_field_problem.dipole_elements = dipole_elements; atom_field_problem.E_field = E_field; atom_field_problem.modulation_freq = modulation_freq; atom_field_problem.freq_index = freq_index; @@ -67,7 +72,8 @@ for detuning_p_cntr=1:N_detun_steps+1; endfor % once we define all problems the main job is done here -kappa_p=cellfun( @susceptibility_steady_state_at_freq, problems_cell_array); +%kappa_p=cellfun( @susceptibility_steady_state_at_freq, problems_cell_array); +kappa_p=parcellfun(2, @susceptibility_steady_state_at_freq, problems_cell_array); figure(1); plot(detuning_freq, real(kappa_p)); title("probe dispersion"); |