ZZ：PWSCF计算实例分析（3）态密度计算

例子，计算Cu的态密度
一、自洽计算（见前一个例子）
二、非自洽计算，增加k点，并采用四面体方法来确定电子的占有数
&control
calculation='nscf'
prefix='cu',
pseudo_dir = './',
outdir='./'
/
&system
ibrav=2, celldm(1) =6.73, nat=1, ntyp=1,
ecutwfc = 25.0, ecutrho = 300, nbnd=8, occupations='tetrahedra'
/
&electrons
conv_thr = 1.0e-8
mixing_beta = 0.7
/
ATOMIC_SPECIES
Cu 63.55 Cu.pz-d-rrkjus.UPF
ATOMIC_POSITIONS
Cu 0.0 0.0 0.0
K_POINTS {automatic}
12 12 12 0 0 0

三、采用dos.x计算总态密度
&inputpp
outdir='./'
prefix='cu'
fildos='cu.dos',
Emin=-5.0, Emax=25.0, DeltaE=0.1
/

四、采用projwfc.x来计算分波态密度
&inputpp
outdir='./'
prefix='cu'
Emin=-5.0, Emax=25.0, DeltaE=0.1 ngauss=1, degauss=0.02
/

在计算态密度的步骤就是如上面所述：a)，先进行自洽计算，保留输出的势、电荷密度和波函数；b)，然后读入上一步自洽计算得到的势或电荷密度或波函数，进行非自洽计算，其中增加k点网格，并采用四面体方法来确定电子占有数; c)，采用dos.x计算总态密度；d)，采用projwfc.x计算分波态密度。

在自洽计算中occupations关键词的赋值已设置为'tetrahedra'表示采用四面体方法确定电子占有数和费米能级。另外，K_POINTS {automatic} 下面的k点设置已增密，设置为12x12x12，为使得计算的态密度光滑，有可能需设置的更密些。

在采用dos.x总态密度计算中，输入文件中由&inputpp 和'/'来之前的关键词来设置，它的关键词有：
outdir用来设置计算上非自洽计算输出文件的目录，设置为'./'表示是当前目录；
prefix用来标记当前所计算的体系，也确定了上一步非自洽计算输出的势或电荷密度或波函数的文件的名称，此例子中设置为'cu'，注意它们的赋值应该与上一步的非自洽计算中的一致。
fildos用来指明所计算的总态密度将写到哪个文件中，此例子中赋值为'cu.dos',表示总态密度将写到cu.dos文件；
Emin用来设置计算态密度时，能量范围的最小值，赋值为-5.0，表示将从E=-5.0 eV开始输出对应的态密度值；
Emax用来设置计算态密度时，能量范围的最大值，赋值为25.0，表示将到E=25.0eV为止输出对应的态密度；
DeltaE用来设置计算态密度时，按多大的能量间隔输出态密度，这里设为0.1eV输出态密度。

在采用projwfc.x计算态密度时，&inputpp 和'/'来之前的关键词来设置，它的关键词与dos.x的输入文件中的关键词差不多：
ngauss用来设置态密度时展宽的方法，这样是为了使得所计算的态密度看起来光滑，可以赋值：
0，表示采用简单的高斯函数
1，表示采用一阶Methfessel-Paxton函数
-1，表示采用Marzari-Vanderbilt“冷离散“方法，其实就是一种函数形式
-99，表示采用Fermi-Dirac函数

degauss用来设置展宽函数中的展宽系数。