格林函数数值计算

假设能谱已经算出来，那么算格林函数是非常容易的

可以直接用向量化编程实现，这样的计算效率最高。

nk = 1024
nw = 1024
delta = 0.01
ks = np.linspace(-pi,pi,nk)
omega = np.linspace(-2,2,nw)
Ks, Omega = np.meshgrid(ks,omega)
Ek = -2*np.cos(Ks)
Gk = np.zeros((nw,nk),dtype='double')
Gk = np.power((Omega+1j*delta-Ek),-1)
Ak = -Gk.imag/pi
plt.pcolormesh(Ks, Omega, Ak)
plt.colorbar()

这里积分只对k进行，我们令\omga = 0，这个公式就可以写成两个格林函数的卷积

如何用卷积的方式解决这个计算。

chi0 = signal.fftconvolve(Gk,Gk,'same')/nkx