| 查看: 6299 | 回复: 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学习资料 |
» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
国社科又开始会评了，不知道这次命运如何
已经有9人回复
-
能否申诉？
已经有5人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
-
面上意见出来了
已经有18人回复
-
找导师
已经有7人回复
-
29号明天会评吗
已经有3人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有11人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 用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