在完成Hatree-Fork自洽迭代后，計算單粒子格林函數和譜函數很容易，只需要如下計算：

這個計算量并不大，是可以接受的程度

nw = 200 delta = 0.002 nk = 200 ks = np.linspace(0,2*pi,nk) omega = np.linspace(-1,1,nw) Gk_up = np.zeros((nw,nk,m,m),dtype='complex') Gk_down = np.zeros((nw,nk,m,m),dtype='complex') Ak_up = np.zeros((nw,nk),dtype='double') Ak_down = np.zeros((nw,nk),dtype='double') N = 32 for i in range(nw):for j in range(nk):Hk0 = Hamiltonian_H0(ks[j],N)Hk_up = Hk0 + U*(Ndown_avg -0.5*np.eye(m))Hk_down = Hk0 + U*(Nup_avg - 0.5*np.eye(m))Gk_up[i,j,:,:] = inv((omega[i]+1j*delta)*np.eye(m)-Hk_up)Gk_down[i,j,:,:] = inv((omega[i]+1j*delta)*np.eye(m)-Hk_down)Ak_up[i,j] = -np.imag(Gk_up[i,j,:,:].trace())/piAk_down[i,j] = -np.imag(Gk_down[i,j,:,:].trace())/pi X,Y = np.meshgrid(ks,omega) plt.pcolormesh(X, Y, Ak_up)

總結

以上是生活随笔為你收集整理的Hatree-Fork单粒子格林函数计算的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。