| 查看: 5558 | 回复: 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学习资料 |
» 猜你喜欢
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
孩子确诊有中度注意力缺陷
已经有6人回复
-
2026博士申请-功能高分子，水凝胶方向
已经有6人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
-
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
用wien2k计算反铁磁详细步骤，举例说明！
已经有12人回复
3楼2013-08-18 16:43:49
2楼2013-08-18 10:08:05
4楼2013-08-18 18:38:45
正要用这个呢5楼2013-08-21 08:03:55