| 查看: 21989 | 回复: 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» 本帖@通知
» 猜你喜欢
民办高校入职后稳定吗？
已经有16人回复
-
又见奇葩学生，可咋办
已经有23人回复
-
Nano-Micro Letters投稿可以用PDF上传吗
已经有3人回复
-
求推荐期刊
已经有4人回复
-
Angew拒稿重投经历大修返回又送审了
已经有5人回复
-
请问什么溶剂能溶解二氧化锰
已经有4人回复
-
蹲一个男朋友
已经有28人回复
-
评审规则突发奇想
已经有24人回复
-
换导师
已经有13人回复
-
来用亲身体会，一起说说年轻老师的辛苦
已经有24人回复
» 本主题相关商家推荐: (我也要在这里推广)
4楼2015-11-02 04:13:53
7楼2015-11-02 08:34:31
24楼2015-11-03 08:33:55
31楼2015-11-04 07:54:13
★★★★★ 五星级,优秀推荐
|
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
5楼2015-11-02 08:24:35
8楼2015-11-02 08:37:04
10楼2015-11-02 09:11:16
11楼2015-11-02 10:23:46
13楼2015-11-02 10:39:14
15楼2015-11-02 12:45:16
22楼2015-11-03 02:35:02
23楼2015-11-03 02:42:08
25楼2015-11-03 10:08:03
26楼2015-11-03 10:23:32
28楼2015-11-03 10:40:30
29楼2015-11-03 10:40:44
32楼2015-11-04 08:31:29
35楼2015-11-04 23:19:22
38楼2015-11-05 18:41:11
39楼2015-11-05 19:44:10
42楼2015-11-06 14:27:29
45楼2015-11-08 14:19:16
50楼2015-11-11 22:51:48
简单回复
我爱天津2楼
2015-11-02 04:02
回复
额
我爱天津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
回复
五星好评 顶一下，感谢分享！
liangpei16楼
2015-11-02 14:17
回复
五星好评 顶一下，感谢分享！
小虫虫LG17楼
2015-11-02 14:42
回复
五星好评 顶一下，感谢分享！
get-it18楼
2015-11-02 14:58
回复
五星好评 顶一下，感谢分享！
Nobelbo19楼
2015-11-02 16:34
回复
一般 顶 发自小木虫Android客户端
ifmc123420楼
2015-11-02 18:34
回复
五星好评 顶一下，感谢分享！
漫天飘雪21楼
2015-11-02 21:19
回复
五星好评 
lwk20527楼
2015-11-03 10:24
回复
五星好评 顶一下，感谢分享！
luchhi30楼
2015-11-03 16:37
回复
五星好评 顶一下，感谢分享！
rdwjf33楼
2015-11-04 11:44
回复
五星好评 顶一下，感谢分享！
muzi198434楼
2015-11-04 21:19
回复
五星好评 顶一下，感谢分享！
lion_tiger36楼
2015-11-05 10:19
回复
五星好评 顶一下，感谢分享！
1xiaohua137楼
2015-11-05 14:35
回复
顶一下，感谢分享！
漫天飘雪40楼
2015-11-05 21:27
回复
顶一下，感谢分享！
qnwbwng41楼
2015-11-06 10:16
回复
三星好评 感谢分享 发自小木虫Android客户端
junfeiyan43楼
2015-11-08 11:32
回复
五星好评 顶一下，感谢分享！
ridge888844楼
2015-11-08 13:33
回复
五星好评 顶一下，感谢分享！
cy1302120246楼
2015-11-08 22:06
回复
五星好评 顶一下，感谢分享！
03一号位47楼
2015-11-10 17:59
回复
五星好评 顶一下，感谢分享！
ridge888848楼
2015-11-11 09:29
回复
顶一下，感谢分享！
jymm49楼
2015-11-11 21:44
回复
五星好评 顶一下！ 发自小木虫IOS客户端
48