seems too start outputing reasonable values at least for 2 levels system

1 files changed, 67 insertions, 0 deletions

diff --git a/three_levels.m b/three_levels.m
new file mode 100644
index 0000000..beb70b3
--- /dev/null
+++ b/three_levels.m
@@ -0,0 +1,67 @@
+levels=1:3;
+
+%      -----------  |1>
+%         /     \
+%    E_d /       \ 
+%       /         \ E_p
+%      /           \
+%    -------- |3>   \
+%                    \
+%               ___________ |2>
+
+
+
+Nlevels=3;
+w12=1e9;
+w_hpf=6800;
+w13=w12-w_hpf;
+
+
+% unperturbed Hamiltonian energy levels
+levels_energy=[ w12, 0,  w_hpf];
+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
+g_decay(1,3)=gamma; %upper level decay
+g_decay(3,2)=gamma_23; % lower levels mixing
+g_decay(2,3)=gamma_23; % lower levels mixing
+
+%defasing matris
+g_deph=0;
+g_dephasing=zeros(Nlevels);
+g_dephasing(1,2)=g_deph;
+g_dephasing(2,1)=g_dephasing(1,2);
+g_dephasing(1,3)=g_deph;
+g_dephasing(3,1)=g_dephasing(1,3);
+
+
+
+% dipole matrix
+dipole_elements=zeros(Nlevels);
+dipole_elements(1,2)=1;
+dipole_elements(2,1)=dipole_elements(1,2);
+dipole_elements(3,1)=1;
+dipole_elements(1,3)=dipole_elements(3,1);
+
+
+%EM field definition
+Ep=10.0;
+Epc=conj(Ep);
+Ed=0;
+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];
+