| 查看: 22086 | 回复: 300 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[资源]
VASP处理能带小程序，非常方便，一秒种搞定！
|
||
|
送给大家一个home made 小程序，由C语言编写，用来处理由VASP得到的自旋极化(ISPIN=2)或非自旋极化(ISPIN=1)的能带，经过测试，还是非常方便的。以下几点供参考： 1、网上一直盛传的pbnd.x(听说是著名的侯博士写的？)在近些年内方便了广大VASP使用者，对我们的科研事业做出了极大的贡献！这套处理程序首先需要由gk.x读取syml来生成KPOINTS，计算后得到EIGENVAL，再由pbnd.x读取syml和EIGENVAL来生成bnd.dat,以方便使用者可以直接导入origin作能带图【1】。美中不足的是，syml作为这套工具的必备文件，需要使用者自己来写，格式比较复杂，对于初学者很容易望而生畏【2】。 2、VASP.5.X的版本出现后，新增加了line-mode格式的KPOINTS写法【3】。这种写法专为能带计算而设计，不需要像5.0之前的版本一样手动生成高对称点的连线。在实际操作中，用过Material Studio的可视化功能，找到体系的高对称点【4】，并添加到KPOINTS中，VASP可以自动生成各点的连线。所以，需要一个新的工具，来将结果转换成可以直接导入origin作图。 3、这次的能带处理工具bandkit，优点在于方便，不需要手动输入任何参数。它从POSCAR里读取基矢，从KPOINTS里读取用line-mode格式取的高对称点，从EIGENVAL中读取能带数据，从DOSCAR中读取费米能级，这些都是自动完成的，使用者只需要打一个bandkit.x的命令，就可以得到可以直接用在文章里作图的数据了！ 4、输出文件解释:band.dat是直接将EIGENVAL传换成origin可读取的格式，未经任何修改。band_modified.dat在纵坐标上，将费米能级标为0；在横坐标上，修正为倒空间的实际长度，单位为A^-1。程度的标准输出即在屏幕上打出的内容，有带隙的大小，band edge所在的位置是line-mode里第几个K点，还有生成的band_modified.dat文件中，各个高对称点所对应的横坐标(这个在origin里作图的时候要用到的！)。 5、对于半导体，屏幕上输出的带隙值经测试认为正确。对于带隙为0的情况，屏幕上输出的带隙不正确(通常会是很小的接近0的一个数)，这种情况我懒得再改了，因为你带隙都为0了，就不需要来看我这个输出的带隙值了。 6、考虑自旋极化与不考虑两种情况经过测试，都可以正常处理。但考虑自旋极化情况中，屏幕上打出的带隙值是spinup 和spindown各自的band gap。在band.dat和band_modified.dat中，第一列为横坐标(k点)，第二列为spinup ，第三列为spindown。 7、发布一个扣扣群：329143489物理化学材料互助群 329143489，版主高抬贵手，飞蝗芜湖！   8、另外还测试过vtst tools里的bandgap.pl，发现会出错，不知道各位是否试过？ 更多了解，请尝试example。 欢迎同行提出问题和建议，什么都可以提，我不骂人，但是有可能拒绝回答。 程序出现错误请提出，我会积极改正！ 【1】http://baike.baidu.com/link?url= ... hkDBJw4ZHLyhG7njYaa 【2】http://muchong.com/html/201404/7339472.html 【3】http://cms.mpi.univie.ac.at/vasp ... e_calculations.html 【4】在MS中build-symmetry-find symmetry和tools-Brillouin Path两步，分别先找到对称性，然后生成高对称点，参见MS说明书  [ Last edited by 我爱天津 on 2015-11-3 at 18:51 ] 之前上传的程序 总有人说解压不了 其实 是我的失误 我应该把文件命名成bandkit.tar.gz 所以说 你gunzip之后 得到的bandkit文件 其实是一个tar! 你可以直接tar -xf bandkit.gz 也可以先gunzip再tar -xf bandkit 如果不懂操作可以加我扣扣 或 扣扣群 咨询~ [ Last edited by 我爱天津 on 2015-12-1 at 17:40 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com - 附件 1 : bandkit.gz
2015-11-02 03:50:46, 83.27 K
» 本帖已获得的红花（最新10朵）
guosophia» 本帖@通知
» 猜你喜欢
ACS Nano投稿后分配到副编辑手里12天了，能确定送审了吗？
已经有11人回复
-
累并快乐着
已经有15人回复
-
审稿人含糊拒稿，还需要回复吗？如何回复？
已经有20人回复
-
硕博巨婴，也许才刚刚开始
已经有85人回复
-
F口信息学部拿面上，大概需要什么样的成果
已经有10人回复
-
34岁读博士晚吗
已经有31人回复
-
2025考博
已经有12人回复
-
太卷了
已经有23人回复
-
杭电、天津科技、青农和宁波工程学院如何选？
已经有23人回复
-
2024国社科通讯评审
已经有23人回复
» 本主题相关商家推荐: (我也要在这里推广)
103楼2016-01-13 22:29:36
4楼2015-11-02 04:13:53
5楼2015-11-02 08:24:35
7楼2015-11-02 08:34:31
简单回复
我爱天津3楼
2015-11-02 04:13
回复
顶一下，感谢分享！
mink6楼
2015-11-02 08:33
回复
五星好评 顶一下，感谢分享！
刘瑞蕊9楼
2015-11-02 08:45
回复
顶一下，感谢分享！
hakuna12楼
2015-11-02 10:27
回复
五星好评 顶一下，感谢分享！
songyumin14楼
2015-11-02 11:07
回复
五星好评 顶一下，感谢分享！
20