def external_field_log_plot():
file2 = FileOperations('Time_5000_117_from0,04to1,2.txt')
file = FileOperations('Field_0,1_2000_16_from1,0to16,0.txt')
_, x1, y1 = file.reading()
y_0 = file2.searching(0.1)
x1 = [i**2 for i in x1]
y1 = [np.log(i / y_0) for i in y1]
plt.plot(x1, y1, 'y')
file = FileOperations('Field_0,2_2000_16_from1,0to16,0.txt')
_, x2, y2 = file.reading()
y_0 = file2.searching(0.2)
x2 = [i**2 for i in x2]
y2 = [np.log(i / y_0) for i in y2]
plt.plot(x2, y2, 'r')
file = FileOperations('Field_0,3_2000_16_from1,0to16,0.txt')
_, x3, y3 = file.reading()
y_0 = file2.searching(0.3)
x3 = [i**2 for i in x3]
y3 = [np.log(i / y_0) for i in y3]
plt.plot(x3, y3, 'b')
file = FileOperations('Field_0,4_2000_16_from1,0to16,0.txt')
_, x4, y4 = file.reading()
y_0 = file2.searching(0.4)
x4 = [i**2 for i in x4]
y4 = [np.log(i / y_0) for i in y4]
plt.plot(x4, y4, 'g')
plt.xlim(0, 260)
plt.ylim(0, 7)
plt.title('Spin relaxation time dependence on the external magnetic field',
fontsize=26)
plt.xlabel('(Ωτ)^2', fontsize=20)
plt.ylabel('Logarithm of ratio of relaxation times', fontsize=20)
plt.legend(['0.1', '0.2', '0.3', '0.4'],
title='ωτ =',
fontsize=20,
title_fontsize=20,
loc='upper left')
plt.tick_params(axis='both', which='major', labelsize=18)
plt.plot(x1, y1, 'yo')
plt.plot(x2, y2, 'ro')
plt.plot(x3, y3, 'bo')
plt.plot(x4, y4, 'go')
plt.savefig(r'image_files\external_field_log_plot.png')
plt.show()
def relaxation_plot():
file = FileOperations('Time_5000_117_from0,04to1,2.txt')
_, x, y = file.reading()
x1, x2, y1, y2 = [], [], [], []
for i in range(len(x)):
if x[i] < 1:
x1.append(x[i])
y1.append(y[i])
if x[i] > 1:
x2.append(x[i])
y2.append(y[i])
plt.plot(x1, y1, 'r', linewidth=2)
plt.plot(x2, y2, 'r', linewidth=2)
plt.plot(x, y, 'ro', markersize=8)
plt.xlim(0, 1.23)
plt.ylim(0, 42)
plt.title('Spin relaxation time dependence on the external magnetic field',
fontsize=26)
plt.xlabel('Time between collisions ( 1/Ω )', fontsize=20)
plt.ylabel('Spin relaxation time ( 1/Ω )', fontsize=20)
plt.tick_params(axis='both', which='major', labelsize=18)
plt.savefig(r'image_files\relaxation_plot.png')
plt.show()
def relaxation_log_log_plot():
file = FileOperations('Time_5000_117_from0,04to1,2.txt')
_, x, y = file.reading()
x = [np.log(i) for i in x]
y = [np.log(i) for i in y]
for i in range(len(x)):
print(x[i], y[i])
x1, x2, y1, y2 = [], [], [], []
for i in range(len(x)):
if x[i] < 0:
x1.append(x[i])
y1.append(y[i])
if x[i] > 0:
x2.append(x[i])
y2.append(y[i])
plt.plot(x1, y1, 'r', linewidth=2)
plt.plot(x2, y2, 'r', linewidth=2)
plt.plot(x, y, 'ro', markersize=8)
plt.xlim(-3.3, 0.3)
plt.ylim(0, 4)
plt.title('Spin relaxation time dependence on the external magnetic field',
fontsize=26)
plt.xlabel('Logarithm of time between collisions', fontsize=20)
plt.ylabel('Logarithm of spin relaxation time', fontsize=20)
plt.tick_params(axis='both', which='major', labelsize=18)
plt.savefig(r'image_files\relaxation_log_log_plot.png')
plt.show()
def external_field_plot():
file2 = FileOperations('Time_5000_117_from0,04to1,2.txt')
file = FileOperations('Field_0,1_2000_16_from1,0to16,0.txt')
_, x1, y1 = file.reading()
plt.plot(x1, y1, 'y')
file = FileOperations('Field_0,2_2000_16_from1,0to16,0.txt')
_, x2, y2 = file.reading()
plt.plot(x2, y2, 'r')
file = FileOperations('Field_0,3_2000_16_from1,0to16,0.txt')
_, x3, y3 = file.reading()
plt.plot(x3, y3, 'b')
file = FileOperations('Field_0,4_2000_16_from1,0to16,0.txt')
_, x4, y4 = file.reading()
plt.plot(x4, y4, 'g')
plt.xlim(0, 17)
plt.ylim(0, 80)
plt.title('Spin relaxation time dependence on the external magnetic field',
fontsize=26)
plt.xlabel('Time between collisions ( 1/ω )', fontsize=20)
plt.ylabel('Spin relaxation time ( 1/ω )', fontsize=20)
plt.legend(['0.1', '0.2', '0.3', '0.4'],
title='ωτ =',
fontsize=20,
title_fontsize=20,
loc='upper left')
plt.tick_params(axis='both', which='major', labelsize=18)
plt.plot(x1, y1, 'yo')
plt.plot(x2, y2, 'ro')
plt.plot(x3, y3, 'bo')
plt.plot(x4, y4, 'go')
plt.savefig(r'image_files\external_fields_plot.png')
plt.show()
'''
'''
file = FileOperations('Field_0,3_2000_16_from1,0to16,0.txt')
_, x1, y1 = file.reading()
file = FileOperations('Field_spec200_0,3_2000_16_from1,0to16,0.txt')
_, x2, y2 = file.reading()
file = FileOperations('Field_spec50_0,3_2000_16_from1,0to16,0.txt')
_, x3, y3 = file.reading()
x = [np.log(i) for i in x1]
y1 = [np.log(i) for i in y1]
y2 = [np.log(i) for i in y2]
y3 = [np.log(i) for i in y3]
plt.plot(x1, y1, 'ro')
plt.plot(x2, y2, 'bo')
plt.plot(x3, y3, 'go')
plt.show()
'''
file = FileOperations('Field_0,3_2000_16_from1,0to16,0.txt')
_, x1, y1 = file.reading()
file = FileOperations('2_Field_0,3_2000_17_from0,0to16,0.txt')
_, x2, y2 = file.reading()
plt.plot(x1, y1, 'ro')
plt.plot(x2, y2, 'bo')
plt.show()
print("Print file name:", end=' ')
file_name = input()
else:
# most useful file for few next plots; write anything instead 0 and 1 to call it
file_name = 'Field_0,2_2000_16_from1,0to16,0.txt'
print("Print 1 for theoretical plot:", end=' ')
theor_existence = input()
if theor_existence == '1':
theor_existence = '_theor_'
else:
theor_existence = '_'
file = FileOperations(file_name)
inf, x, y = file.reading()
# creating plot
plt.figure(figsize=(19, 9.5))
default_label_size = 20
default_title_size = 26
y_0 = 0
y_theor = []
if plot_type == 'T':
y_theor = time_theor()
plt.title('Spin relaxation time dependence on time between collisions',
fontsize=default_title_size)
plt.xlabel('Time between collisions ( 1/Ω )', fontsize=default_label_size)
plt.ylabel('Spin relaxation time ( 1/Ω )', fontsize=default_label_size)