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

---
 two_levels.m | 54 ++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 54 insertions(+)
 create mode 100644 two_levels.m

(limited to 'two_levels.m')

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