24小时热门版块排行榜

返回列表

cj4566

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 17
散金: 24
红花: 3
帖子: 473
在线: 399.8小时
虫号: 356133
注册: 2007-04-26
性别: GG
专业: 无机非金属类光电信息与功

[求助] pwscf画能带图求助

我试着用pwscf画ZnO的能带图，但是出的图很奇怪，没看到能隙，而且纵坐标都是负值。刚学习qe，大家指点一下啊，谢谢了！

下面是输入文件：
1. relax （考虑relax的时候也进行scf，所以后面没有重新scf计算）。
&CONTROL
            calculation  = 'relax',
            restart_mode = 'from_scratch',
            prefix    = "ZnObulk",
            pseudo_dir = "/pkg/suse11/quantum/5.0/pseudo",
            outdir    = "/home/n7388683/QE/tmp",
                  tstress = .true
                  tprnfor = .true
/
&SYSTEM
                     ibrav = 4,
               celldm(1) = 6.211531542958419,
               celldm(3) = 1.61271676300578,
                     nat = 4,
                     ntyp = 2,
               occupations = 'smearing',
                  degauss = 0.02D0,
                  ecutwfc = 30,
                  ecutrho = 300,
/
&ELECTRONS
               conv_thr = 1.D-8,
               mixing_beta = 0.6,
/
&IONS
            ion_dynamics = 'bfgs',
/
ATOMIC_SPECIES
Zn 65.3799972534  Zn.pbe-van.UPF
O 15.9989995956  O.pbe-van_ak.UPF
ATOMIC_POSITIONS crystal
  O 0.3333333333333330 0.6666666666666671 0.3757029511412930
  O  -0.3333333333333330  -0.6666666666666671 0.8757029511412930
Zn 0.3333333333333330 0.6666666666666671 0.0000000000000000
Zn  -0.3333333333333330  -0.6666666666666671 0.5000000000000000
K_POINTS automatic
  10 10 10 0 0 0

2. 能带计算
&CONTROL
            calculation  = 'bands',
            pseudo_dir = "/pkg/suse11/quantum/5.0/pseudo",
            outdir    = "/home/n7388683/QE/tmp",
            prefix    = "ZnObulk",
/
&SYSTEM
                     ibrav = 4,
               celldm(1) = 6.211531542958419,
               celldm(3) = 1.61271676300578,
                     nat = 4,
                     ntyp = 2,
                     nbnd = 18,
                  ecutwfc = 30,
                  ecutrho = 300,
/
&ELECTRONS
         diagonalization = 'david'
/
ATOMIC_SPECIES
Zn 65.3799972534  Zn.pbe-van.UPF
O 15.9989995956  O.pbe-van_ak.UPF
ATOMIC_POSITIONS crystal
O       0.333333333 0.666666667 0.377379183
O    -0.333333333  -0.666666667 0.877379183
Zn    0.333333333 0.666666667  -0.001676232
Zn    -0.333333333  -0.666666667 0.498323768
K_POINTS
31
-0.333333  0.666667  0.500000  1.00
-0.333333  0.666667  0.416667  1.00
-0.333333  0.666667  0.333333  1.00
-0.333333  0.666667  0.250000  1.00
-0.333333  0.666667  0.166667  1.00
-0.333333  0.666667  0.083333  1.00
-0.333333  0.666667  0.000000  1.00
-0.277778  0.638889  0.000000  1.00
-0.222222  0.611111  0.000000  1.00
-0.166667  0.583333  0.000000  1.00
-0.111111  0.555556  0.000000  1.00
-0.055556  0.527778  0.000000  1.00
  0.000000  0.500000  0.000000  1.00
  0.000000  0.416667  0.000000  1.00
  0.000000  0.333333  0.000000  1.00
  0.000000  0.250000  0.000000  1.00
  0.000000  0.166667  0.000000  1.00
  0.000000  0.083333  0.000000  1.00
  0.000000  0.000000  0.000000  1.00
  0.000000  0.000000  0.083333  1.00
  0.000000  0.000000  0.166667  1.00
  0.000000  0.000000  0.250000  1.00
  0.000000  0.000000  0.333333  1.00
  0.000000  0.000000  0.416667  1.00
  0.000000  0.000000  0.500000  1.00
  0.000000  0.083333  0.500000  1.00
  0.000000  0.166667  0.500000  1.00
  0.000000  0.250000  0.500000  1.00
  0.000000  0.333333  0.500000  1.00
  0.000000  0.416667  0.500000  1.00
  0.000000  0.500000  0.500000  1.00

3. 提取数据
&bands
prefix  = 'ZnObulk'
outdir = "/home/n7388683/QE/tmp",
filband = 'ZnObulk.bands.dat'
/

4. 画图
ZnObulk.bands.dat
-10 10
ZnObulk.bands.xmgr
ZnObulk.bands.ps
7.1264
1.0 7.1264

回复此楼