应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 学习PWSCF(一)：能带计算步骤与讨论
291/3返回列表123下一页
查看: 9354  |  回复: 28
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

HSE06

新虫 (小有名气)

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

从昨天下午开始学习pwscf，倒腾了这么长时间，终于把能带图做出来了，没想到这么难搞。把做法贴上来和大家分享。作为一个初学者，在这个过程中也有许多问题，请大家指正。
第一步： scf计算

输入文件si.scf.in内容如下：
CODE:
&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

执行命令：pw.x < si.scf.in > si.scf.out

第二步：nscf计算（这一步和上一步scf有什么关系？需要读入上面的什么文件信息吗？）
输入文件si.nscf.in内容如下：
CODE:
&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

运行命令:pw.x < si.nscf.in > si.nscf.out

第三步：band计算（这一步和上一步nscf有什么关系？为什么需要nscf这一步？）
输入文件内容si.bands.in
CODE:
&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

运行命令:pw.x < si.bands.in > si.bands.out

第四步：（我不知道这一步应该叫什么名字）
编辑文件bands.x内容如下
CODE:
&bands
    prefix  = 'silicon'
    outdir='./'
    filband = 'bands.dat'
/

运行命令:bands.x < bands.in > bands.out

第五步:貌似这一步的目的是画出能带图？
运行plotband.x
CODE:
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

苍天啊大地啊，终于做出来了能带图，即si.bands.ps，如下所示：
可以看到我做出来的能带图有些毛刺，这是怎么产生的？求教。多谢！学习PWSCF(一)：能带计算步骤与讨论
1.jpg

[ Last edited by HSE06 on 2013-6-6 at 17:23 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com
  • 附件 1 : 1.jpg
    • 2013-06-06 17:15:53, 27.37 K

» 收录本帖的淘贴专辑推荐

精华网帖收集 第一性原理 第一性原理相关文档 initio
pwscf 第一性原理计算 科研经验

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2013-06-06 17:22:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

KalaShayminS

铁杆木虫 (著名写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
sunyang1988: 金币+2, 谢谢交流 2013-06-06 18:08:43
nbnd：计算用的能带数。决定能画出的能带数量。
3是权重。一般计算不需要手工设。
nscf计算需要读取上一步产生的outdir中的东西。这一步nscf是用来读取Fermi能量的。
k点指定是作为示例，实际计算k点可以取别的。
毛刺是bands.x本身的缺陷导致的。
deltaE是画图坐标轴的刻度。
赞一下(3人) 回复此楼
2楼2013-06-06 17:42:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

HSE06

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

多谢!
nbnd决定计算能带的数目，是不是scf步就不需要指定了？在band.x指定能带数就可以了吧？
赞一下 回复此楼
3楼2013-06-06 18:32:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
3楼: Originally posted by HSE06 at 2013-06-06 18:32:31
多谢!
nbnd决定计算能带的数目，是不是scf步就不需要指定了？在band.x指定能带数就可以了吧？...

需要在最后一步pw.x的时候指定。
赞一下 回复此楼
4楼2013-06-07 10:32:44
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

HSE06

新虫 (小有名气)
引用回帖:
4楼: Originally posted by KalaShayminS at 2013-06-07 10:32:44
需要在最后一步pw.x的时候指定。...

好的，多谢
回复此楼
5楼2013-06-07 11:08:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuanghua304

禁虫 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖内容被屏蔽
6楼2013-06-26 10:48:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuanghua304

禁虫 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖内容被屏蔽
7楼2013-06-26 10:59:28
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

HSE06

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

通过指定目录就可以了
赞一下 回复此楼
8楼2013-06-26 15:18:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

HSE06

新虫 (小有名气)
引用回帖:
7楼: Originally posted by chuanghua304 at 2013-06-26 10:59:28
我看到的教程不是算能带用calculation='bands' ，你怎么用了两次calculation = 'nscf' ,啊？

第一个nsc是读取scf的文件
第二次nscf是生成si.bands.out
赞一下 回复此楼
9楼2013-06-26 15:19:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuanghua304

禁虫 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖内容被屏蔽
10楼2013-06-26 17:04:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 HSE06 的主题更新
291/3返回列表123下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件) 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 博后面上今天有bug可以看到是否资助？ +21 lyfbangong 2024-06-12 32/1600 2024-06-17 15:11 by 伐木打太极
[教师之家] 饶议：什么制度能保障大学普通教师不用为领导拎包，不用看领导脸色 +10 zju2000 2024-06-12 16/800 2024-06-17 14:54 by 大天尊
[文学芳草园] 累并快乐着 +14 MYHLD521 2024-06-14 14/700 2024-06-17 14:51 by shl2112501
[论文投稿] 审稿问题：为什么荧光激发波长和紫外吸收波长差的大？ 10+5 sdawege 2024-06-14 9/450 2024-06-17 13:15 by 人族首领
[公派出国] 去英国的小伙伴儿都在哪儿租的房子呀？ +3 65syn 2024-06-14 4/200 2024-06-17 11:24 by 爱吃菜的虫
[催化] 镍负载氧化铝的保存问题 8+3 lwn0130 2024-06-15 4/200 2024-06-17 10:48 by adaihao
[基金申请] 面青地会评时间 +7 tanjydd 2024-06-15 7/350 2024-06-17 10:40 by 你想干啥
[博后之家] 山东大学(青岛)“天然药物生物智造”课题组 招聘“博士后”(年薪20.4-55.6万元) +5 第二种态度 2024-06-11 11/550 2024-06-17 10:08 by 小懂事k
[有机交流] 车间生产，真空度很高，温度很高，但减压蒸馏速度很慢。 10+12 召唤鬼泣lL 2024-06-13 37/1850 2024-06-17 02:17 by xmclyn
[找工作] 应届大龄土博教职选择，石河子大学和广西大学 10+4 茅山老道士 2024-06-10 34/1700 2024-06-16 22:05 by yugege2009
[考博] 上海交大招收材料化学方向科研助理/“申请考核”博士（请勿回复帖子或站内投条） +3 灵梦and紫 2024-06-12 4/200 2024-06-15 20:58 by 1822836277
[基金申请] BO4的YQ答辩通知发布了吗？ +6 博学笃行 2024-06-11 6/300 2024-06-15 16:04 by 悲催科研狗
[基金申请] 为什么我的博后基金还在流动站审核中？不会是学院给我卡了吧？ +14 王凯12 2024-06-13 26/1300 2024-06-15 15:22 by 好人与坏人
[论文投稿] 投了一篇4区的SCI，审稿人一个拒稿，一个小修，编辑给了大修。 +9 安稳22123 2024-06-13 10/500 2024-06-14 23:45 by jurkat.1640
[基金申请] 面上基金有一个（两个）C是不是就没戏了.... 5+3 zzzm116 2024-06-13 17/850 2024-06-14 22:26 by lzt8076
[基金申请] 博士后面上项目状态还是专家评审吗 10+9 Thatcheremu 2024-06-13 55/2750 2024-06-13 21:23 by 乌合麒麟
[硕博家园] 科研求助 +5 杲www 2024-06-12 6/300 2024-06-13 16:16 by 姓李名明
[论文投稿] 摩擦磨损论文投稿 +3 jmysan 2024-06-12 3/150 2024-06-13 08:36 by 莱茵润色
[基金申请] 博后特助这周出结果吗？往年都是啥时候啊？ +13 jsqy 2024-06-12 17/850 2024-06-12 19:55 by Lynn212
[硕博家园] 申博 +3 悦悦小小鱼 2024-06-12 3/150 2024-06-12 15:11 by chen5805
信息提示
关闭
请填处理意见
关闭