[交流] 学习PWSCF(一)：能带计算步骤与讨论 已有14人参与

&CONTROL
                 calculation = 'scf' ,
                      outdir = './' ,
                  pseudo_dir = './' ,
                      prefix = 'silicon' ,
/
&SYSTEM
                       ibrav = 2,
                   celldm(1) = 10.2,
                         nat = 2,
                        ntyp = 1,
                     ecutwfc = 12.0 ,
                        nbnd = 8, （这个参数，能带的个数，是做能带的时候必须指定的吗？）
            exxdiv_treatment = 'gygi-baldereschi' ,
/
&ELECTRONS
/
ATOMIC_SPECIES
   Si   28.08600  Si.vbc.UPF
ATOMIC_POSITIONS
   Si      0.000000000    0.000000000    0.000000000
   Si      0.250000000    0.250000000    0.250000000
K_POINTS
2
   0.250000000    0.250000000    0.750000000      3.000000000    （最后面的这个3是什么意思？感觉很诡异啊）
   0.250000000    0.250000000    0.250000000      1.000000000

&CONTROL
                 calculation = 'nscf' ,
                      outdir = './' ,
                  pseudo_dir = './' ,
                      prefix = 'silicon' ,
/
&SYSTEM
                       ibrav = 2,
                   celldm(1) = 10.2,
                         nat = 2,
                        ntyp = 1,
                     ecutwfc = 12.0 ,
                        nbnd = 8,
            exxdiv_treatment = 'gygi-baldereschi' ,
/
&ELECTRONS
/
ATOMIC_SPECIES
   Si   28.08600  Si.vbc.UPF
ATOMIC_POSITIONS alat
   Si      0.000000000    0.000000000    0.000000000
   Si      0.250000000    0.250000000    0.250000000
K_POINTS tpiba
3
   0.000000000    0.000000000    0.000000000      1.000000000  （这个k点是怎么指定的？为什么这样？）
   1.000000000    0.000000000    0.000000000      2.000000000
   0.500000000    0.500000000    0.500000000      3.000000000

&CONTROL
                 calculation = 'nscf' ,
                      outdir = './' ,
                  pseudo_dir = './' ,
                      prefix = 'silicon' ,
/
&SYSTEM
                       ibrav = 2,
                   celldm(1) = 10.2,
                         nat = 2,
                        ntyp = 1,
                     ecutwfc = 12.0 ,
                        nbnd = 8,
            exxdiv_treatment = 'gygi-baldereschi' ,
/
&ELECTRONS
/
ATOMIC_SPECIES
   Si   28.08600  Si.vbc.UPF
ATOMIC_POSITIONS alat
   Si      0.000000000    0.000000000    0.000000000
   Si      0.250000000    0.250000000    0.250000000
K_POINTS tpiba
36
   0.500000000    0.500000000    0.500000000      1.000000000
   0.400000000    0.400000000    0.400000000      2.000000000
   0.300000000    0.300000000    0.300000000      3.000000000
   0.200000000    0.200000000    0.200000000      4.000000000
   0.100000000    0.100000000    0.100000000      5.000000000
   0.000000000    0.000000000    0.000000000      6.000000000
   0.000000000    0.000000000    0.100000000      7.000000000
   0.000000000    0.000000000    0.200000000      8.000000000
   0.000000000    0.000000000    0.300000000      9.000000000
   0.000000000    0.000000000    0.400000000     10.000000000
   0.000000000    0.000000000    0.500000000     11.000000000
   0.000000000    0.000000000    0.600000000     12.000000000
   0.000000000    0.000000000    0.700000000     13.000000000
   0.000000000    0.000000000    0.800000000     14.000000000
   0.000000000    0.000000000    0.900000000     15.000000000
   0.000000000    0.000000000    1.000000000     16.000000000
   0.000000000    0.100000000    1.000000000     17.000000000
   0.000000000    0.200000000    1.000000000     18.000000000
   0.000000000    0.300000000    1.000000000     19.000000000
   0.000000000    0.400000000    1.000000000     20.000000000
   0.000000000    0.500000000    1.000000000     21.000000000
   0.000000000    0.600000000    1.000000000     22.000000000
   0.000000000    0.700000000    1.000000000     23.000000000
   0.000000000    0.800000000    1.000000000     24.000000000
   0.000000000    0.900000000    1.000000000     25.000000000
   0.000000000    1.000000000    1.000000000     26.000000000
   0.000000000    0.900000000    0.900000000     27.000000000
   0.000000000    0.800000000    0.800000000     28.000000000
   0.000000000    0.700000000    0.700000000     29.000000000
   0.000000000    0.600000000    0.600000000     30.000000000
   0.000000000    0.500000000    0.500000000     31.000000000
   0.000000000    0.400000000    0.400000000     32.000000000
   0.000000000    0.300000000    0.300000000     33.000000000
   0.000000000    0.200000000    0.200000000     34.000000000
   0.000000000    0.100000000    0.100000000     35.000000000
   0.000000000    0.000000000    0.000000000     36.000000000

&bands
    prefix  = 'silicon'
    outdir='./'
    filband = 'bands.dat'
/

input file > bands.dat
Reading 8 bands at 36 k-points
Range: -5.6680 16.4950eV Emin, Emax > -6.0 10.0 （这个范围是指能带的能量窗口吗？）
high-symmetry point: 0.5000 0.5000 0.5000
high-symmetry point: 0.0000 0.0000 0.0000
high-symmetry point: 0.0000 0.0000 1.0000
high-symmetry point: 0.0000 1.0000 1.0000
high-symmetry point: 0.0000 0.0000 0.0000
output file (xmgr) > si.bands.xmgr
bands in xmgr format written to file si.bands.xmgr
output file (ps) > si.bands.ps
Efermi > 6.337          （这个费米能级从哪个文件得到？求关键字……）
deltaE, reference E (for tics) 1.0, 6.337           （deltaE是什么意思啊？后面的E是费米能级吧？）
bands in PostScript format written to file si.bands.ps

2楼: Originally posted by KalaShayminS at 2013-06-06 17:42:47
nbnd：计算用的能带数。决定能画出的能带数量。
3是权重。一般计算不需要手工设。
nscf计算需要读取上一步产生的outdir中的东西。这一步nscf是用来读取Fermi能量的。
k点指定是作为示例，实际计算k点可以取别的。 ...

3楼: Originally posted by HSE06 at 2013-06-06 18:32:31
多谢!
nbnd决定计算能带的数目，是不是scf步就不需要指定了？在band.x指定能带数就可以了吧？...

4楼: Originally posted by KalaShayminS at 2013-06-07 10:32:44
需要在最后一步pw.x的时候指定。...

6楼: Originally posted by chuanghua304 at 2013-06-26 10:48:09
按说计算能带和态密度要读入SCF计算的电荷密度，但是输入文件里面好像没有体现啊？是不是放在同一目录下就行了啊？

7楼: Originally posted by chuanghua304 at 2013-06-26 10:59:28
我看到的教程不是算能带用calculation='bands' ，你怎么用了两次calculation = 'nscf' ,啊？