从昨天下午开始学习pwscf,倒腾了这么长时间,终于把能带图做出来了,没想到这么难搞。把做法贴上来和大家分享。作为一个初学者,在这个过程中也有许多问题,请大家指正。
第一步: scf计算
输入文件si.scf.in内容如下: 执行命令:pw.x < si.scf.in > si.scf.out
第二步:nscf计算(这一步和上一步scf有什么关系?需要读入上面的什么文件信息吗?)
输入文件si.nscf.in内容如下: 运行命令:pw.x < si.nscf.in > si.nscf.out
第三步:band计算(这一步和上一步nscf有什么关系?为什么需要nscf这一步?)
输入文件内容si.bands.in 运行命令:pw.x < si.bands.in > si.bands.out
第四步:(我不知道这一步应该叫什么名字)
编辑文件bands.x内容如下 运行命令:bands.x < bands.in > bands.out
第五步:貌似这一步的目的是画出能带图?
运行plotband.x 苍天啊大地啊,终于做出来了能带图,即si.bands.ps,如下所示:
可以看到我做出来的能带图有些毛刺,这是怎么产生的?求教。多谢!
1.jpg
[ Last edited by HSE06 on 2013-6-6 at 17:23 ] |