| 查看: 6490 | 回复: 12 | |||||||||
| 【奖励】 本帖被评价7次，作者xinqing98增加金币 5.2 个 | |||||||||
[资源]
以石墨烯为例Pwscf计算能带的详细步骤和要点
|
|||||||||
|
使用PWSCF计算能带需要以下几个步骤： 1. 构建好结构以后进行scf计算，以附件中的Graphene.scf.in为例，其中重要的是在calculation='scf' 这是告诉系统进行自洽计算。构建好以后，执行： pw.x < Graphene.scf.in > Graphene.scf.out & 2. Graphene.band.in 在这个文件中，calculation='bands'，这是告诉系统做能带计算。另外还要在 &system中加入nbnd = 这是绘制的能带数目。 这是最关键的一个步骤：其中K_POINTS一定要写成“K_POINTS crystal”，这样就可以直接采用xcrystalden中给 出的倒空间坐标（和vasp做能带时的K点坐标一样），但如果在K_POINTS后面什么也不写， 则程序默认的是“K_POINTS tpiba”，这是以2pi/a位单位的cartesian坐标（如果还输入倒空间坐标，产生的K点就不是我们需要的K点，这样会得到错误的能带。） 这里我选的高对称点为： G(0.0,0,0,0.0) K (0.3333,0.3333,0) M (0,0.5,0) ,K点是处在BZ顶角的点， M是处在BZ边上中点的点。 执行：pw.x < Graphene.band.in > Graphene.band.out & 3. 编辑bands.in文件，注意其中的prefix和outdir的数值一定要和前两步完全一样，否则将无法读取前两步 的结果。 执行：bands.x < bands.in > bands.out & 4. 最后一步就可以绘制能带图了： 执行：plotband.x 输入bands.in中指定的*.dat文件，然后根据它自动生成的能带范围写入绘制能带的范围，可大可小。 然后自己明名*.xmgr文件和*.ps文件，再输入费米能级，第一次可以输入0.0即可，然后输入能带 纵坐标间隔和费米能级，全部结束！ 注意，绘制完以后，可以找出自己需要的费米能级，gnuplot 绘制*.xmgr文件，可以找出自己需要的费米能级，然后再执行plotband.x重新绘制一次能带即可。 [ Last edited by xinqing98 on 2013-8-17 at 22:09 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 仿真建模与计算 | 第一性原理 | 第一性原理相关文档 | 第一性原理计算经验 |
| 第一性原理计算 | 方法不错 | QE学习资料 |
» 猜你喜欢
极限拉扯，如何抉择
已经有15人回复
-
申请科研项目纯靠实力能成功吗？
已经有8人回复
-
中科院海洋所招收海洋材料相关领域博士研究生及研究助理
已经有4人回复
-
40+中年人失业后，还能找到好工作吗？
已经有3人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 用wien2k计算反铁磁详细步骤，举例说明！
已经有12人回复
2楼2013-08-18 10:08:05
3楼2013-08-18 16:43:49
4楼2013-08-18 18:38:45
正要用这个呢 5楼2013-08-21 08:03:55
6楼2013-08-21 16:10:12
7楼2013-08-23 15:30:26
8楼2016-09-30 20:39:00
很赞！话说附件，没看到
9楼2018-10-15 10:21:22
10楼2019-01-02 11:46:25
11楼2020-05-24 08:08:35
12楼2020-05-24 17:58:25
13楼2020-05-24 19:44:14
10