沿高对称线计算磁化率画图Python代码

import numpy as np
import matplotlib.pyplot as plt

chi = np.loadtxt("chi0.txt")
Nk = chi.size
KPoints = [0, Nk/4, Nk/2, 3/5 * Nk, 2*Nk/3, 4/5*Nk, Nk]
Q= [Nk/4]
plt.plot(chi, color = 'black')
plt.xticks(KPoints, (r'${\Gamma}$', 'q1', 'M', 'q2', 'K', 'q3', r'${\Gamma}$'))
for k in KPoints:
    plt.vlines(k, 0, 3, colors = "c", linestyles = "dashed")
#plt.ylim(0, 3)
plt.show()

峰值依赖于费米面。