电荷密度波通过计算极化率判断nesting波矢量

计算公式

引入一个小虚部，防止发散

以一维链模型计算，显然在费米面附近发散，这就是费米面嵌套机制

a = 1;
delta = 0.001;
T = 0.026;
sus = [];
for q=0:0.01:2*pi
    x = 0;
    for k=0:0.01:2*pi
        E1 = 2*cos(k*a);
        E2 = 2*cos((k+q)*a);
        f1 = (exp(E1/T)+1)^-1;
        f2 = (exp(E2/T)+1)^-1;
        x = x+ f1*(1-f2)/(E2 - E1 + 1i*delta);
    end
    sus =[sus,x];
end
q=0:0.01:2*pi;
plot(q,real(sus))

波矢量q = pi