plt.figure(num=None, figsize=(8, 13), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
#ax.plot(root_array[:,0], np.real(root_array[:,1]), '.', color = 'b')
#ax.plot(root_array[:,0], np.imag(root_array[:,1]), '.', color = 'r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], y_range[-1])
ax.set_ylabel(r'$\bar c_{ph}$', fontsize=30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize=30)
line1 = sol.line_trace_sp(slab.disp_rel,
1,
1.301,
0.01,
0,
4,
args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel,
1,
1.083,
0.01,
0,
4,
args=(slab.K, slab.M_A))
line3 = sol.line_trace_sp(slab.disp_rel,
2.32,
1.89,
0.01,
0.1,
plt.figure(num=None, figsize=(8, 13), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
#ax.plot(root_array[:,0], np.real(root_array[:,1]), '.', color = 'b')
#ax.plot(root_array[:,0], np.imag(root_array[:,1]), '.', color = 'r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], y_range[-1])
ax.set_ylabel(r'$\bar c_{ph}$', fontsize=30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize=30)
line1 = sol.line_trace_sp(slab.disp_rel,
0.64,
0.75,
0.01,
0.01,
4,
args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel,
0.55,
0.565,
0.001,
0,
4,
args=(slab.K, slab.M_A))
line3 = sol.line_trace_sp(slab.disp_rel,
2,
0.914,
0.001,
0.06,
plt.figure(num=None, figsize=(8, 13), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
#ax.plot(root_array[:,0], np.real(root_array[:,1]), '.', color = 'b')
#ax.plot(root_array[:,0], np.imag(root_array[:,1]), '.', color = 'r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], y_range[-1])
ax.set_ylabel(r'$\bar c_{ph}$', fontsize=30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize=30)
line1 = sol.line_trace_sp(slab.disp_rel,
1,
1.262,
0.01,
0.14,
4,
args=(slab.K, slab.M_A))
line2 = sol.line_trace_mp(slab.disp_rel_mp,
2,
0.786 + 0.05 * 1j,
0.01,
0,
4,
args=(slab.K, slab.M_A))
line3 = sol.line_trace_sp(slab.disp_rel,
1,
0.643,
0.01,
0.05,
font = {'size': 15}
matplotlib.rc('font', **font)
plt.figure(num=None, figsize=(8, 11), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
ax.plot(root_array[:,0], np.real(root_array[:,1]), '.', color = 'b')
ax.plot(root_array[:,0], np.imag(root_array[:,1]), '.', color = 'r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], y_range[-1])
ax.set_ylabel(r'$\bar c_{ph}$', fontsize = 30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize = 30)
line1 = sol.line_trace_sp(slab.disp_rel, 0.48, 0.483, 0.01, 0, 4, args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel, 0.32, 0.315, 0.01, 0, 4, args=(slab.K, slab.M_A))
line3 = sol.line_trace_sp(slab.disp_rel, 0.88, 0.66, 0.01, 0, 4, args=(slab.K, slab.M_A))
line4 = sol.line_trace_mp(slab.disp_rel_mp, 1.28, 0.557, 0.01, 0.01, 4, args=(slab.K, slab.M_A))
line5 = sol.line_trace_mp(slab.disp_rel_mp, 3.68, 0.56, 0.01, 0.01, 4, args=(slab.K, slab.M_A))
line6 = sol.line_trace_mp(slab.disp_rel_mp, 3.84, 0.542, 0.01, 0.01, 4, args=(slab.K, slab.M_A))
ax.plot(line1[:,0], line1[:,1], color = 'b', linewidth = 2, linestyle='-', label='Sausage')
ax.plot(line2[:,0], line2[:,1], color = 'g', linewidth = 2, linestyle='-', label='Kink')
ax.plot(line3[:,0], line3[:,1], color = 'b', linewidth = 2, linestyle='-')
ax.plot(line4[:,0], line4[:,1], color = 'g', linewidth = 2, linestyle='-')
ax.plot(line5[:,0], line5[:,1], color = 'b', linewidth = 2, linestyle='-')
ax.plot(line6[:,0], line6[:,1], color = 'g', linewidth = 2, linestyle='-')
plt.legend(loc=4)
plt.figure(num=None, figsize=(8, 11), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
ax.plot(root_array[:, 0], np.real(root_array[:, 1]), '.', color='b')
ax.plot(root_array[:, 0], np.imag(root_array[:, 1]), '.', color='r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], y_range[-1])
ax.set_ylabel(r'$\bar c_{ph}$', fontsize=30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize=30)
line1 = sol.line_trace_sp(slab.disp_rel,
0.72,
0.5,
0.01,
0.01,
4,
args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel,
0.48,
0.4,
0.01,
0,
4,
args=(slab.K, slab.M_A))
line3 = sol.line_trace_mp(slab.disp_rel_mp,
0.96,
0.56,
0.01,
0.01,
plt.figure(num=None, figsize=(8, 13), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
#ax.plot(root_array[:,0], np.real(root_array[:,1]), '.', color = 'b')
#ax.plot(root_array[:,0], np.imag(root_array[:,1]), '.', color = 'r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], 2.13)
ax.set_ylabel(r'$\bar c_{ph}$', fontsize=30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize=30)
line1 = sol.line_trace_sp(slab.disp_rel,
0.8,
0.72,
0.01,
0,
4,
args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel,
0.8,
0.657,
0.01,
0,
4,
args=(slab.K, slab.M_A))
line3 = sol.line_trace_sp(slab.disp_rel,
1.36,
1.58,
0.01,
0,
plt.figure(num=None, figsize=(8, 13), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
#ax.plot(root_array[:,0], np.real(root_array[:,1]), '.', color = 'b')
#ax.plot(root_array[:,0], np.imag(root_array[:,1]), '.', color = 'r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], y_range[-1])
ax.set_ylabel(r'$\bar c_{ph}$', fontsize=30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize=30)
line1 = sol.line_trace_sp(slab.disp_rel,
0.15,
0.215,
0.002,
0,
4,
args=(slab.K, slab.M_A))
line101 = line1[:, 0][np.abs(np.imag(line1[:, 1])) < 1e-5]
line111 = line1[:, 1][np.abs(np.imag(line1[:, 1])) < 1e-5]
line102 = line1[:, 0][np.abs(np.imag(line1[:, 1])) > 1e-5]
line112 = line1[:, 1][np.abs(np.imag(line1[:, 1])) > 1e-5]
break_pt = int(
np.where(
np.array(
[line101[i + 1] - line101[i]
for i in np.arange(len(line101) - 1)]) > 0.1)[0])
line2 = sol.line_trace_sp(slab.disp_rel,
font = {'size': 15}
matplotlib.rc('font', **font)
plt.figure(num=None, figsize=(8, 13), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
#ax.plot(root_array[:,0], np.real(root_array[:,1]), '.', color = 'b')
#ax.plot(root_array[:,0], np.imag(root_array[:,1]), '.', color = 'r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], y_range[-1])
ax.set_ylabel(r'$\bar c_{ph}$', fontsize = 30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize = 30)
line1 = sol.line_trace_sp(slab.disp_rel, 0.18, 0.208, 0.002, 0, 4, args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel, 1, 0.63, 0.01, 0.41, 4, args=(slab.K, slab.M_A))
line3 = sol.line_trace_sp(slab.disp_rel, 1, -0.09, 0.01, 0, 4, args=(slab.K, slab.M_A))
line4 = sol.line_trace_sp(slab.disp_rel, 1, -0.568, 0.01, 0.37, 4, args=(slab.K, slab.M_A))
line5 = sol.line_trace_sp(slab.disp_rel, 1.96, -0.689, 0.001, 1.907, 4, args=(slab.K, slab.M_A))
line6 = sol.line_trace_sp(slab.disp_rel, 1.8, 0.359, 0.01, 0, 4, args=(slab.K, slab.M_A))
line7 = sol.line_trace_sp(slab.disp_rel, 0.8, 0.394, 0.01, 0, 4, args=(slab.K, slab.M_A))
line8 = sol.line_trace_sp(slab.disp_rel, 0.3, 0.381, 0.001, 0.15, 0.403, args=(slab.K, slab.M_A))
ax.plot(line1[:,0][np.abs(np.imag(line1[:,1])) < 1e-5],
line1[:,1][np.abs(np.imag(line1[:,1])) < 1e-5],
color = 'b', linewidth = 2, linestyle='-', label='Stable')
ax.plot(line1[:,0][np.abs(np.imag(line1[:,1])) > 1e-5],
line1[:,1][np.abs(np.imag(line1[:,1])) > 1e-5],
color = 'b', linewidth = 2, linestyle='--', label='Unstable')
ax.plot(line1[:,0][np.abs(np.imag(line1[:,1])) < 1e-5],
plt.figure(num=None, figsize=(8, 13), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
#ax.plot(root_array[:,0], np.real(root_array[:,1]), '.', color = 'b')
#ax.plot(root_array[:,0], np.imag(root_array[:,1]), '.', color = 'r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], y_range[-1])
ax.set_ylabel(r'$\bar c_{ph}$', fontsize=30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize=30)
line1 = sol.line_trace_sp(slab.disp_rel,
1,
0.666,
0.01,
0,
4,
args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel,
1,
0.572,
0.01,
0,
4,
args=(slab.K, slab.M_A))
line3 = sol.line_trace_sp(slab.disp_rel,
1,
-0.41,
0.01,
0,
plt.figure(num=None, figsize=(8, 13), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
#ax.plot(root_array[:,0], np.real(root_array[:,1]), '.', color = 'b')
#ax.plot(root_array[:,0], np.imag(root_array[:,1]), '.', color = 'r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], y_range[-1])
ax.set_ylabel(r'$\bar c_{ph}$', fontsize=30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize=30)
line1 = sol.line_trace_sp(slab.disp_rel,
1,
1.23,
0.01,
0.12,
4,
args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel,
1,
0.863,
0.002,
0,
4,
args=(slab.K, slab.M_A))
line3 = sol.line_trace_sp(slab.disp_rel,
3.2,
0.085,
0.01,
0.1,
plt.figure(num=None, figsize=(8, 11), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
ax.plot(root_array[:, 0], np.real(root_array[:, 1]), '.', color='b')
ax.plot(root_array[:, 0], np.imag(root_array[:, 1]), '.', color='r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], y_range[-1])
ax.set_ylabel(r'$\bar c_{ph}$', fontsize=30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize=30)
line1 = sol.line_trace_sp(slab.disp_rel,
0.4,
0.891,
0.01,
0,
4,
args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel,
0.35,
0.47,
0.001,
0.002,
4,
args=(slab.K, slab.M_A))
line3 = sol.line_trace_mp(slab.disp_rel_mp,
1.36,
0.651,
0.01,
0.01,
plt.figure(num=None, figsize=(8, 13), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
#ax.plot(root_array[:,0], np.real(root_array[:,1]), '.', color = 'b')
#ax.plot(root_array[:,0], np.imag(root_array[:,1]), '.', color = 'r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], y_range[-1])
ax.set_ylabel(r'$\bar c_{ph}$', fontsize=30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize=30)
line1 = sol.line_trace_sp(slab.disp_rel,
1,
1.875,
0.01,
0.87,
4,
args=(slab.K, slab.M_A))
line3 = sol.line_trace_sp(slab.disp_rel,
1,
1.01,
0.01,
0,
4,
args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel,
1,
0.887,
0.01,
0,
plt.figure(num=None, figsize=(8, 13), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
#ax.plot(root_array[:,0], np.real(root_array[:,1]), '.', color = 'b')
#ax.plot(root_array[:,0], np.imag(root_array[:,1]), '.', color = 'r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], y_range[-1])
ax.set_ylabel(r'$\bar c_{ph}$', fontsize=30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize=30)
line1 = sol.line_trace_sp(slab.disp_rel,
1,
1.975,
0.01,
0,
4,
args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel,
1,
1.228,
0.01,
0.1,
4,
args=(slab.K, slab.M_A))
line3 = sol.line_trace_sp(slab.disp_rel,
1,
1.028,
0.01,
0,
font = {'size': 15}
matplotlib.rc('font', **font)
plt.figure(num=None, figsize=(8, 13), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
#ax.plot(root_array[:,0], np.real(root_array[:,1]), '.', color = 'b')
#ax.plot(root_array[:,0], np.imag(root_array[:,1]), '.', color = 'r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], y_range[-1])
ax.set_ylabel(r'$\bar c_{ph}$', fontsize = 30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize = 30)
line1 = sol.line_trace_sp(slab.disp_rel, 1, 0.681, 0.01, 0, 4, args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel, 1, 0.578, 0.01, 0, 4, args=(slab.K, slab.M_A))
line3 = sol.line_trace_sp(slab.disp_rel, 1, -0.308, 0.01, 0, 4, args=(slab.K, slab.M_A))
line4 = sol.line_trace_sp(slab.disp_rel, 1, -0.674, 0.01, 0, 4, args=(slab.K, slab.M_A))
line5 = sol.line_trace_sp(slab.disp_rel, 2, -0.837, 0.01, 1.78, 4, args=(slab.K, slab.M_A))
line6 = sol.line_trace_sp(slab.disp_rel, 3.8, -0.839, 0.001, 3.6, 4, args=(slab.K, slab.M_A))
line7 = sol.line_trace_sp(slab.disp_rel, 2.72, 0.08, 0.002, 0.05, 4, args=(slab.K, slab.M_A))
ax.plot(line1[:,0], line1[:,1], color = 'b', linewidth = 2, linestyle='-', label='Sausage')
ax.plot(line2[:,0], line2[:,1], color = 'g', linewidth = 2, linestyle='-', label='Kink')
ax.plot(line2[:,0], np.imag(line2[:,1]), color = 'r', linewidth = 2, linestyle='-', label='Imaginary')
ax.plot(line3[:,0], line3[:,1], color = 'g', linewidth = 2, linestyle='-')
ax.plot(line4[:,0], line4[:,1], color = 'b', linewidth = 2, linestyle='-')
ax.plot(line5[:,0], line5[:,1], color = 'g', linewidth = 2, linestyle='-')
ax.plot(line6[:,0], line6[:,1], color = 'b', linewidth = 2, linestyle='-')
ax.plot(line7[:,0], line7[:,1], color = 'b', linewidth = 2, linestyle='-')
fontsize=20)
ax.annotate(r'$\bar c_0 + M_A$',
xy=(-0.13, slab.c0 - 0.02),
xycoords='data',
annotation_clip=False,
fontsize=20)
ax.annotate(r'$-\bar c_0 + M_A$',
xy=(-0.16, -slab.c0 - 0.04),
xycoords='data',
annotation_clip=False,
fontsize=20)
line1 = sol.line_trace_sp(slab.disp_rel,
0.045,
0.567,
0.001,
0,
0.216,
args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel,
0.045,
0.539,
0.001,
0,
0.54,
args=(slab.K, slab.M_A))
line3 = sol.line_trace_sp(slab.disp_rel,
1.02,
0.733,
0.001,
0.937,
fontsize=20)
ax.annotate(r'$\bar c_0 + M_A$',
xy=(-0.13, slab.c0 - 0.02),
xycoords='data',
annotation_clip=False,
fontsize=20)
ax.annotate(r'$-\bar c_0 + M_A$',
xy=(-0.16, -slab.c0 - 0.04),
xycoords='data',
annotation_clip=False,
fontsize=20)
line1 = sol.line_trace_sp(slab.disp_rel,
0.0975,
0.8839,
0.005,
0,
0.3232,
args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel,
0.6,
0.8197,
0.01,
-0.01,
1.5,
args=(slab.K, slab.M_A))
line3 = sol.line_trace_sp(slab.disp_rel,
0.6,
0.5205,
0.001,
0,
plt.figure(num=None, figsize=(8, 13), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
#ax.plot(root_array[:,0], np.real(root_array[:,1]), '.', color = 'b')
#ax.plot(root_array[:,0], np.imag(root_array[:,1]), '.', color = 'r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], y_range[-1])
ax.set_ylabel(r'$\bar c_{ph}$', fontsize=30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize=30)
line1 = sol.line_trace_sp(slab.disp_rel,
0.48,
0.748,
0.01,
0.01,
4,
args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel,
0.48,
0.585,
0.01,
0,
4,
args=(slab.K, slab.M_A))
line3 = sol.line_trace_sp(slab.disp_rel,
1.84,
1.5,
0.01,
1,
font = {'size': 15}
matplotlib.rc('font', **font)
plt.figure(num=None, figsize=(8, 13), dpi=80, facecolor='w', edgecolor='k')
ax = plt.subplot()
#ax.plot(root_array[:,0], np.real(root_array[:,1]), '.', color = 'b')
#ax.plot(root_array[:,0], np.imag(root_array[:,1]), '.', color = 'r')
ax.set_xlim(x_range[0], x_range[-1])
ax.set_ylim(y_range[0], y_range[-1])
ax.set_ylabel(r'$\bar c_{ph}$', fontsize = 30, rotation=0)
ax.set_xlabel(r'$k x_0$', fontsize = 30)
line1 = sol.line_trace_sp(slab.disp_rel, 1, 0.635, 0.01, 0.37, 4, args=(slab.K, slab.M_A))
line2 = sol.line_trace_sp(slab.disp_rel, 1, 0.495, 0.01, 0, 4, args=(slab.K, slab.M_A))
line3 = sol.line_trace_sp(slab.disp_rel, 1, -0.383, 0.01, 0, 4, args=(slab.K, slab.M_A))
line4 = sol.line_trace_sp(slab.disp_rel, 1, -0.64, 0.01, 0.55, 4, args=(slab.K, slab.M_A))
line5 = sol.line_trace_sp(slab.disp_rel, 2.32, -0.159, 0.005, 0.06, 4, args=(slab.K, slab.M_A))
line6 = sol.line_trace_sp(slab.disp_rel, 2.24, -0.12, 0.005, 0.17, 4, args=(slab.K, slab.M_A))
ax.plot(line1[:,0], line1[:,1], color = 'b', linewidth = 2, linestyle='-', label='Quasi-Sausage')
ax.plot(line2[:,0], line2[:,1], color = 'g', linewidth = 2, linestyle='-', label='Quasi-Kink')
ax.plot(line2[:,0], np.imag(line2[:,1]), color = 'r', linewidth = 2, linestyle='-', label='Imaginary')
ax.plot(line3[:,0], line3[:,1], color = 'g', linewidth = 2, linestyle='-')
ax.plot(line4[:,0], line4[:,1], color = 'b', linewidth = 2, linestyle='-')
ax.plot(line5[:,0], line5[:,1], color = 'g', linewidth = 2, linestyle='-')
ax.plot(line6[:,0], line6[:,1], color = 'b', linewidth = 2, linestyle='-')
ax.plot([x_range[0], x_range[-1]], [slab.c0 + slab.M_A, slab.c0 + slab.M_A],