二维费米面的nesting计算

在六角形区域中，随机采样2W个点，根据公式，计算费米面nesting，使用Orgin画图公式对多边形插值处理

主函数

delta = 0.01;
T = 0.0026;
N = 20000;
k=GetKPoints(N);
sus = zeros(N,1);
for i=1:N                    
    for j=1:N
        Hk1=Hamiltonian_TaS2_k(k(j,1),k(j,2));
        [~,Ek1]=eig(Hk1);
        E1 = Ek1(1,1);
        f1 = (exp(E1/T)+1)^-1;
        Hk2=Hamiltonian_TaS2_k(k(j,1)+k(i,1),k(j,2)+k(i,2));
        [~,Ek2]=eig(Hk2);
        E2 = Ek2(1,1);
        f2 = (exp(E2/T)+1)^-1;
        sus(i,1) = sus(i,1) + f1*(1-f2)/(E2 - E1 + 1i*delta);
    end
end
plot3(k(:,1),k(:,2),real(sus),'.');
data = [k,real(sus)];

优化升级方案：

取更多的点N = 4W

计算下面平行四边形区域的布里渊区，均匀取点