From 35af744d95df0d687bccdaaa4dabda6fe830bed7 Mon Sep 17 00:00:00 2001
From: Eugeniy Mikhailov <evgmik@gmail.com>
Date: Wed, 9 Dec 2009 05:34:20 +0000
Subject: seems too start outputing reasonable values at least for 2 levels
 system

---
 three_levels.m | 67 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 67 insertions(+)
 create mode 100644 three_levels.m

(limited to 'three_levels.m')

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];
+
-- 
cgit v1.2.3