| 查看: 2328 | 回复: 6 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
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  |
回复此楼» 猜你喜欢
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有20人回复
-
今天务委会开完了，明天出结果吗
已经有22人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】pwscf计算得到的能带图很乱
已经有16人回复
- 【求助】pwscf计算中的warning要紧吗？如何消除呢？
已经有9人回复
- 【求助】pwscf错误提示diagonalization(ZHEGV*) falied
已经有8人回复
- 【求助】求助pwscf一些参数的物理意义和设定经验
已经有13人回复
- 【求助】利用abinit或pwscf软件在单机上可以计算半导体材料的哪些性质?
已经有8人回复
- 【求助】画能带图
已经有9人回复
- 【求助】关于PWSCF介电常数和声子谱的计算
已经有7人回复
- 【求助】pwscf的若干问题
已经有4人回复
- 【求助】pwscf中的问题
已经有4人回复
- 【求助】如何使用VASP画费米面。
已经有13人回复
wujiandongde
木虫 (小有名气)
- 应助: 11 (小学生)
- 金币: 4265.4
- 散金: 30
- 红花: 5
- 帖子: 182
- 在线: 704.8小时
- 虫号: 941604
- 注册: 2010-01-13
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
7楼2018-08-27 14:37:30
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
2楼2012-08-19 11:12:03
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
3楼2012-08-19 11:14:36
4楼2012-08-19 11:41:04
20