diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-09 05:34:20 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-09 05:34:20 +0000 |
| commit | 35af744d95df0d687bccdaaa4dabda6fe830bed7 (patch) | |
| tree | 79fe4fbcc33b0e2ee8166a1f6f7d00b4c7c19826 /two_levels.m | |
| parent | e48d9efc84dd2e701fd67e154cac2569d7e2eced (diff) | |
| download | multi_mode_eit-35af744d95df0d687bccdaaa4dabda6fe830bed7.tar.gz multi_mode_eit-35af744d95df0d687bccdaaa4dabda6fe830bed7.zip | |
seems too start outputing reasonable values at least for 2 levels system
Diffstat (limited to 'two_levels.m')
| -rw-r--r-- | two_levels.m | 54 |
1 files changed, 54 insertions, 0 deletions
diff --git a/two_levels.m b/two_levels.m new file mode 100644 index 0000000..18860cb --- /dev/null +++ b/two_levels.m @@ -0,0 +1,54 @@ +1; + +% ----------- |1> +% \ +% \ +% \ E_p +% \ +% \ +% \ +% ___________ |2> + + + +Nlevels=2; +w12=1e4; +w_hpf=68; +w13=w12-w_hpf; + + +% unperturbed Hamiltonian energy levels +levels_energy=[ w12, 0]; +levels_energy=levels_energy*hbar; +H0=zeros(Nlevels); +H0=diag(levels_energy); +%for i=1:Nlevels + %H0(i,i)=levels_energy(i); +%endfor + +% decay matrix g(i,j) correspnds to decay from i-->j +gamma=6; +gamma_23=.01; +g_decay=zeros(Nlevels); +g_decay(1,2)=gamma; %upper level decay + +%defasing matris +g_deph=0; +g_dephasing=zeros(Nlevels); +g_dephasing(1,2)=g_deph; +g_dephasing(2,1)=g_dephasing(1,2); + + +% dipole matrix +dipole_elements=zeros(Nlevels); +dipole_elements(1,2)=1; +dipole_elements(2,1)=dipole_elements(1,2); + + +%EM field definition +Ep=1.0; +Epc=conj(Ep); +wp=w12; +modulation_freq=[0, wp, -wp]; +E_field =[0, Ep, Epc]; + |