版块导航: 正在加载中...

客户端APP下载

论文辅导

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (6853)
>虫友互识 (989)
>休闲灌水 (575)
>文献求助 (470)
>基金申请 (392)
>论文投稿 (344)
>招聘信息布告栏 (310)
>博后之家 (153)
>硕博家园 (122)
>导师招生 (105)
>电化学 (93)
>材料综合 (89)
>教师之家 (75)
>找工作 (64)
>考博 (56)
>海外博后 (43)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

HSE06

新虫 (小有名气)

应助: 3 (幼儿园)
金币: 7.3
红花: 2
帖子: 71
在线: 67小时
虫号: 1660094
注册: 2012-03-02
性别: GG
专业: 凝聚态物性 II ：电子结构

[交流] 学习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，如下所示：
可以看到我做出来的能带图有些毛刺，这是怎么产生的？求教。多谢！

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	第一性原理计算	科研经验

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

pwscf是如何计算debye温度的已经有7人回复
pwscf计算integrated charge出错及不收敛已经有13人回复
pwscf并行计算时总是隔一段时间重新开始已经有5人回复
pwscf能不能计算电荷转移即计算每个原子的mulliken charge 已经有7人回复
请问Pwscf可以在单机上运行吗已经有4人回复
求PWSCF使用说明书已经有3人回复
PWscf计算homo lumo或者valence conduction band，问题已解决已经有14人回复
【求助】pwscf结构优化结果已经有4人回复
【求助】计算得到能量后，参考点是什么呢？已经有6人回复
【求助】画能带图已经有9人回复
【原创】PWSCF中的nbnd和degauss 已经有11人回复
【求助】关于PWSCF介电常数和声子谱的计算已经有7人回复
【求助】pwscf中的问题已经有4人回复
【求助】用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算已经有23人回复
【求助】请教PWSCF晶体建模和自洽计算，不收敛已经有21人回复
【讨论】交流讨论下用VASP计算不同温度下材料的热力学性能已经有19人回复

1楼2013-06-06 17:22:09

已阅回复此楼关注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的回帖

chuanghua304

禁虫 (职业作家)

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

本帖内容被屏蔽

10楼2013-06-26 17:04:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

chuanghua304

禁虫 (职业作家)

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

本帖内容被屏蔽

11楼2013-06-26 17:14:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

chuanghua304

禁虫 (职业作家)

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

本帖内容被屏蔽

12楼2013-06-27 11:41:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

chuanghua304

禁虫 (职业作家)

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

本帖内容被屏蔽

14楼2013-06-28 10:12:44

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

chuanghua304

禁虫 (职业作家)

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

本帖内容被屏蔽

18楼2013-07-06 08:17:09

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 HSE06 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[公派出国] 去英国的小伙伴儿都在哪儿租的房子呀？ +5	65syn 2024-06-14	7/350	2024-06-17 19:03 by 鱼翔浅底1
[高分子] 寻找聚酯反应釜 +3	茕茕恭煮 2024-06-15	4/200	2024-06-17 18:18 by sunliyan123
[找工作] 药学硕士找不到工作，打算去做科研助理了 +5	pom戴墨镜 2024-06-14	11/550	2024-06-17 17:16 by pom戴墨镜
[基金申请] 博后面上今天有bug可以看到是否资助？ +22	lyfbangong 2024-06-12	33/1650	2024-06-17 16:17 by carolloo
[基金申请] 化学会年会改成一年开一次了？ +6	babu2015 2024-06-12	14/700	2024-06-17 14:41 by 我想更懂你29
[硕博家园] 博士毕业高校和就业的相关问题 +6	SCITOPPP 2024-06-14	10/500	2024-06-17 10:08 by SCITOPPP
[找工作] 应届大龄土博教职选择，石河子大学和广西大学 10+4	茅山老道士 2024-06-10	34/1700	2024-06-16 22:05 by yugege2009
[基金申请] 博士后创新人才支持计划公示 +9	aishida144 2024-06-14	15/750	2024-06-16 09:52 by msjy
[考博] 上海交大招收材料化学方向科研助理/“申请考核”博士（已招满） +3	灵梦and紫 2024-06-12	4/200	2024-06-15 20:58 by 1822836277
[论文投稿] 求机械类四区sci推荐 5+3	迷茫小旷 2024-06-14	4/200	2024-06-15 11:25 by bobvan
[基金申请] 博后基金，以往的结果点不开，怎么回事呢？最后一次机会了，两次都没中前面。 +7	kyukitu 2024-06-14	13/650	2024-06-15 06:46 by 我是王小帅
[论文投稿] 投了一篇4区的SCI，审稿人一个拒稿，一个小修，编辑给了大修。 +9	安稳22123 2024-06-13	10/500	2024-06-14 23:45 by jurkat.1640
[食品] 食品博士导师 +6	小李醒yy 2024-06-11	9/450	2024-06-14 23:37 by 小李醒yy
[基金申请] E12面上申请 +4	汉风之遗 2024-06-13	4/200	2024-06-14 15:28 by 天外飞去来
[有机交流] ππ堆积会发生在有机溶剂中吗 5+3	zibuyu0420 2024-06-13	4/200	2024-06-14 14:17 by 小肉干
[基金申请] 国自然基金公布的时候基金号有吗 +8	潇洒怡惜 2024-06-13	11/550	2024-06-14 11:24 by JRfei
[基金申请] 工材E10口函评结束了吗 10+3	我1的飞翔 2024-06-13	5/250	2024-06-14 06:35 by nono2009
[有机交流] 原料反应完了，怎么知道是产物还是中间体 +6	小胡在努力 2024-06-11	8/400	2024-06-13 13:33 by 091602
[论文投稿] with editor日期变更 +3	慎独的小花卷 2024-06-12	8/400	2024-06-13 11:00 by 慎独的小花卷
[基金申请] 博后特助这周出结果吗？往年都是啥时候啊？ +13	jsqy 2024-06-12	17/850	2024-06-12 19:55 by Lynn212