函数f4(qz)在qz=1e3*kf时已经为0,将f4在(1,Inf)积分,按说,积分上限取到1e3*kf就可以了,事实上,积分上限在10*kf到1e3*kf之间时,积分结果是一致的,但是如果把上限放的非常大,(1,1e8*kf)这时候,积分结果变化又非常大。
后进末学,求大贤指点一二,感激涕零,顿首再拜!
global EF A D h2 x Epsilon0 Epsilons Epsilonh e m m2 Ns kb T B DP ul ut rho E0 kf K2 KL KT;
kb=1.38e-23;
T=3; %温度
h2=1.055e-34;
x=0.3; %Mg%
Epsilon0=8.8542e-12; %真空介电常数
Epsilons=7.9; % 间接
Epsilonh=3.7; %
a1=(3.250+0.05*x)*1e-10; c1=(5.204-0.17*x)*1e-10; a2=3.250e-10;c2=5.204e-10; %
e=1.6022e-19; %电子电荷
DP=3.8*e;
rho=5760;
ul=5e3; %
ut=2.1e3; %
cl=14.1e10; %
ct=2.47e10; %
h15=-0.37;%压电系数间接引用
h31=-0.62;
h33=0.96;
ela2=(8*(2*h15+h31)^2+12*(2*h15+h31)*h33+15*h33^2)/105; %
eta2=(6*(h33-h15-h31)^2+16*(h33-h15-h31)*h15+48*h15^2)/105;
K2=ela2/(Epsilon0*Epsilons*cl)+eta2/(Epsilon0*Epsilons*ct); %
KL=ela2/(Epsilon0*Epsilons*cl);
KT=eta2/(Epsilon0*Epsilons*ct);
m0=9.11e-31; m=0.26*m0; %自由电子
Vc=0.9*2.145*x*e;%
sigmat=-0.054-0.038*x; sigmai=0.038*x;sigmab=0.054;
deltaz=0.1e-9;%变分间隔
sigmad=1e17/1e-6*e*deltaz;%一
Ns=1.4e15;
B=(12*m*e^2/(Epsilon0*Epsilons*h2^2))^(1/3)*(11*Ns/32)^(1/3);
kf=sqrt(2*pi*Ns);%费米
E0=h2^2*B^2/(8*m)+3*e^2/(Epsilon0*Epsilons*B)*(11*Ns/16);%
EF=E0+pi*h2^2*Ns/m;%
f4=@(qz)B.^6./(B.^2+qz.^2).^3.*DP.^2.*h2.*sqrt(kf.^2+qz.^2)./(2.*rho.*ul).*(1./(exp(h2.*ul.*sqrt(kf.^2+qz.^2)./(kb.*T))-1)).*(1-1./(exp((h2.^2.*kf.^2./(2.*m)+E0+h2.*ul.*sqrt(kf.^2+qz.^2)-EF)./(kb.*T))+1));
>> quadgk(f4,1,10*kf)
ans =
2.8824e-063
>> quadgk(f4,1,100*kf)
ans =
2.8824e-063
>> quadgk(f4,1,1e8*kf)
ans =
4.0394e-094