From c680b440e6963ab1f507ca0fcefdc3ce982db8af Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Mon, 21 Dec 2009 22:48:07 +0000 Subject: calculation of total absorption is added --- useful_functions.m | 40 ++++++++++++++++++++++++++++++++++++++++ 1 file changed, 40 insertions(+) (limited to 'useful_functions.m') diff --git a/useful_functions.m b/useful_functions.m index 22c54be..f799f90 100644 --- a/useful_functions.m +++ b/useful_functions.m @@ -258,4 +258,44 @@ function [xi_linear, xi_left, xi_right]=susceptibility_steady_state_at_freq( ato xi_left = xi.left; endfunction +function [dEdz_coef_linear, dEdz_coef_left, dEdz_coef_right] =dEdz_at_freq(wi, rhoLiouville, dipole_elements, E_field) +% complex absorption coefficient +% at given E_field frequency component + Nlevels=( size(dipole_elements.linear)(1) ); + rho=rhoOfFreq(rhoLiouville, wi, Nlevels); + + dEdz_coef_linear = sum( sum( transpose(dipole_elements.linear) .* rho ) ); + dEdz_coef_left = sum( sum( transpose(dipole_elements.left) .* rho ) ); + dEdz_coef_right = sum( sum( transpose(dipole_elements.right) .* rho ) ); + +endfunction + +function [dEdz_linear, dEdz_left, dEdz_right]=dEdz( atom_field_problem) +% complex absorption coefficient for each field + global atom_properties; + L0m = atom_properties.L0m ; + polarizability_m = atom_properties.polarizability_m ; + dipole_elements = atom_properties.dipole_elements ; + + E_field = atom_field_problem.E_field ; + modulation_freq = atom_field_problem.modulation_freq ; + Nfreq=length(modulation_freq); + + rhoLiouville=rhoLiouville_steady_state(L0m, polarizability_m, E_field, modulation_freq); + for wi=1:Nfreq + [dEdz_linear(wi), dEdz_left(wi), dEdz_right(wi)] = dEdz_at_freq(wi, rhoLiouville, dipole_elements, E_field); + endfor +endfunction + +function total_absorption=total_field_absorption(atom_field_problem) +% summed across all frequencies field absorption coefficient + global atom_properties; + [dEdz_linear, dEdz_left, dEdz_right]=dEdz( atom_field_problem); + dEdz_total= dEdz_linear .* conj(dEdz_linear) ... + +dEdz_left .* conj(dEdz_left) ... + +dEdz_right .* conj(dEdz_right); + total_absorption = sum(dEdz_total); +endfunction + + % vim: ts=2:sw=2:fdm=indent -- cgit v1.2.3