-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathlind_quan.py
113 lines (106 loc) · 3.49 KB
/
lind_quan.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
import numpy as np
import matplotlib.pyplot as plt
q_re_rho = [[],[]]
q_im_rho = [[],[]]
l_re_rho = [[],[]]
l_im_rho = [[],[]]
q_re_rho[0].append( np.empty(0, dtype=float) )
q_re_rho[0].append( np.empty(0, dtype=float) )
q_re_rho[1].append( np.empty(0, dtype=float) )
q_re_rho[1].append( np.empty(0, dtype=float) )
q_im_rho[0].append( np.empty(0, dtype=float) )
q_im_rho[0].append( np.empty(0, dtype=float) )
q_im_rho[1].append( np.empty(0, dtype=float) )
q_im_rho[1].append( np.empty(0, dtype=float) )
l_re_rho[0].append( np.empty(0, dtype=float) )
l_re_rho[0].append( np.empty(0, dtype=float) )
l_re_rho[1].append( np.empty(0, dtype=float) )
l_re_rho[1].append( np.empty(0, dtype=float) )
l_im_rho[0].append( np.empty(0, dtype=float) )
l_im_rho[0].append( np.empty(0, dtype=float) )
l_im_rho[1].append( np.empty(0, dtype=float) )
l_im_rho[1].append( np.empty(0, dtype=float) )
tmp = []
with open('Quan/re_rho_00.txt', 'r') as f:
for line in f:
tmp.append(float(line.split()[0]))
q_re_rho[0][0] = np.array( tmp )
tmp = []
with open('Quan/re_rho_01.txt', 'r') as f:
for line in f:
tmp.append(float(line.split()[0]))
q_re_rho[0][1] = np.array( tmp )
tmp = []
with open('Quan/re_rho_10.txt', 'r') as f:
for line in f:
tmp.append(float(line.split()[0]))
q_re_rho[1][0] = np.array( tmp )
tmp = []
with open('Quan/re_rho_11.txt', 'r') as f:
for line in f:
tmp.append(float(line.split()[0]))
q_re_rho[1][1] = np.array( tmp )
tmp = []
with open('Quan/im_rho_01.txt', 'r') as f:
for line in f:
tmp.append(float(line.split()[0]))
q_im_rho[0][1] = np.array( tmp )
tmp = []
with open('Quan/im_rho_10.txt', 'r') as f:
for line in f:
tmp.append(float(line.split()[0]))
q_im_rho[1][0] = np.array( tmp )
tmp = []
with open('Lind/re_rho_00.txt', 'r') as f:
for line in f:
tmp.append(float(line.split()[0]))
l_re_rho[0][0] = np.array( tmp )
tmp = []
with open('Lind/re_rho_01.txt', 'r') as f:
for line in f:
tmp.append(float(line.split()[0]))
l_re_rho[0][1] = np.array( tmp )
tmp = []
with open('Lind/re_rho_10.txt', 'r') as f:
for line in f:
tmp.append(float(line.split()[0]))
l_re_rho[1][0] = np.array( tmp )
tmp = []
with open('Lind/re_rho_11.txt', 'r') as f:
for line in f:
tmp.append(float(line.split()[0]))
l_re_rho[1][1] = np.array( tmp )
tmp = []
with open('Lind/im_rho_01.txt', 'r') as f:
for line in f:
tmp.append(float(line.split()[0]))
l_im_rho[0][1] = np.array( tmp )
tmp = []
with open('Lind/im_rho_10.txt', 'r') as f:
for line in f:
tmp.append(float(line.split()[0]))
l_im_rho[1][0] = np.array( tmp )
dt = 0.002
t_fin = 0.5
K = int(t_fin/dt)+1
t = np.linspace(0,t_fin,K)
plt.plot(t, q_re_rho[0][0] ,'r-', label=r'$\rho_{ee}$')
plt.plot(t, q_re_rho[1][1] ,'b-', label=r'$\rho_{gg}$')
plt.plot(t, q_re_rho[0][1] ,'k-', label=r'Re[$\rho_{eg}$]')
plt.plot(t, q_re_rho[1][0] ,'g-', label=r'Re[$\rho_{ge}$]')
plt.plot(t, q_im_rho[0][1] ,'c-', label=r'Im[$\rho_{eg}$]')
plt.plot(t, q_im_rho[1][0] ,'m-', label=r'Im[$\rho_{ge}$]')
plt.plot(t, l_re_rho[0][0] ,'r--')
plt.plot(t, l_re_rho[1][1] ,'b--')
plt.plot(t, l_re_rho[0][1] ,'k--')
plt.plot(t, l_re_rho[1][0] ,'g--')
plt.plot(t, l_im_rho[0][1] ,'c--')
plt.plot(t, l_im_rho[1][0] ,'m--')
plt.legend(fontsize=20.0)
plt.xlim([0,t_fin])
plt.xlabel(r'$t$', fontsize=20.0)
plt.ylabel(r'$\rho(t)$', fontsize=20.0)
plt.xticks(fontsize=20.0)
plt.yticks(fontsize=20.0)
plt.title(r'$\omega$=3, $\Gamma$=1, $n$=10', fontsize=20.0)
plt.show()