diff options
| -rw-r--r-- | field_description.m | 13 | ||||
| -rw-r--r-- | liouville.m | 4 | ||||
| -rw-r--r-- | three_levels.m | 14 |
3 files changed, 17 insertions, 14 deletions
diff --git a/field_description.m b/field_description.m new file mode 100644 index 0000000..901bfc2 --- /dev/null +++ b/field_description.m @@ -0,0 +1,13 @@ +1; + +%EM field definition +Ep=1.0; +Epc=conj(Ep); +Ed=10; +Edc=conj(Ed); +wd=w13; +wp=w12; +modulation_freq=[0, wp, wd, -wp, -wd, wp-wd, wd-wp]; +E_field =[0, Ep, Ed, Epc, Edc, 0, 0 ]; +%modulation_freq=[0, wp, -wp]; +%E_field =[0, Ep, Epc]; diff --git a/liouville.m b/liouville.m index 5987925..bc9d524 100644 --- a/liouville.m +++ b/liouville.m @@ -6,9 +6,13 @@ useful_functions; hbar=1; im_one=0+1i; + +% load atom energy levels and decay description three_levels; %two_levels; +% load EM field description +field_description; Nfreq=length(modulation_freq); diff --git a/three_levels.m b/three_levels.m index eb9c24e..82a216c 100644 --- a/three_levels.m +++ b/three_levels.m @@ -15,20 +15,6 @@ w13=w12-w_hpf; % \ % ___________ |2> -%EM field definition -Ep=1.0; -Epc=conj(Ep); -Ed=10; -Edc=conj(Ed); -wd=w13; -wp=w12; -modulation_freq=[0, wp, wd, -wp, -wd, wp-wd, wd-wp]; -E_field =[0, Ep, Ed, Epc, Edc, 0, 0 ]; -%modulation_freq=[0, wp, -wp]; -%E_field =[0, Ep, Epc]; - - - % unperturbed Hamiltonian energy levels levels_energy=[ w12, 0, w_hpf]; |