| 查看: 22059 | 回复: 300 | |||||||||||||||||||||
| 【奖励】 本帖被评价224次，作者我爱天津增加金币 174.8 个 | |||||||||||||||||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||||||||||||||||
[资源]
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
» 收录本帖的淘贴专辑推荐
 纳米技术与能源及模拟 | 精华网帖收集 | 仿真建模与计算 | VASP and MS |
| VASP | 第一性原理相关文档 | VASP软件从入门到高深 | vasp计算资料 |
| 计算-vasp | 第一原理 | 计算 | 资源软件 |
| 分子模拟专题 | 材料的门道 | VASP学习 | 各种各样的软件 |
| 科研 | 第一性原理学习 | 第一性原理计算 |
» 本帖已获得的红花（最新10朵）
guosophia» 本帖@通知
» 猜你喜欢
有没有机械的前辈分享一下评上海优都是什么成果啊
已经有3人回复
-
迷茫毕业季
已经有10人回复
-
交友
已经有3人回复
-
计算机类，请问eswa与kbs哪个相对来说容易一些啊？
已经有4人回复
-
博士找工作真难
已经有22人回复
-
中者为王-观某校级项目为领导瓜分有感
已经有30人回复
-
副教授跳槽
已经有24人回复
-
求助
已经有9人回复
-
看到大家一直在说代表作NC，不明觉厉
已经有24人回复
-
能在两广找到什么工作呢
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
★★★★★ 五星级,优秀推荐
|
bandkit小程序读取带隙还是比较方便的，原来vasp5.3的版本可以，现在5.4.1的读取带隙是负， 明显不对了，应该版本输出的格式有变化，楼主有时间看看，谢谢！ > bandkit.x Fermi level is read from DOSCAR as: -3.829716 max of VBM is 12.000000 eV,locates at the 123th kpoint min of CBM is -20.520606 eV,locates at the 29th kpoint band gap is -32.520606 eV |
233楼2017-03-09 22:43:05
268楼2017-11-25 17:35:24
50