混合表象下单粒子格林函数计算

G = 1/（z-H(k))

其中H(k)一个方向作傅里叶变换，另一个方向不做， 直接对矩阵进行求逆。

Matlab代码

nw = 1000; delta = 0.002; omega = linspace(-1,1,nw); A = zeros(nw,nk); for i=1:nwfor j=1:nkH0 = Hamiltonian_H0(k(j),t,N);H = H0+U*Nupavg-0.5*U*eye(m,m);G = ((omega(i)+1j*delta)*eye(m,m)-H)^-1;A(i,j) = -imag(trace(G));end end [X,Y] = meshgrid(k,omega); mesh(X,Y,A);

测试通过

总结

以上是生活随笔为你收集整理的混合表象下单粒子格林函数计算的全部内容，希望文章能够帮你解决所遇到的问题。

如果觉得生活随笔网站内容还不错，欢迎将生活随笔推荐给好友。