24小时热门版块排行榜    

查看: 22061  |  回复: 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朵)

» 本帖@通知

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘瑞蕊

金虫 (正式写手)


★★★★★ 五星级,优秀推荐

顶!d=====( ̄▽ ̄*)b
8楼2015-11-02 08:37:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 301 个回答

我爱天津

金虫 (正式写手)


怎么只有我自己顶 来人顶顶啊 来人啊
4楼2015-11-02 04:13:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我爱天津

金虫 (正式写手)


5楼2015-11-02 08:24:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gsq920306278

至尊木虫 (著名写手)


★ 一星级,一般

请问这脚本的具体细节操作,是在VASP软件里安装,还是不用安装,就像P4VASP软件,直接读取数据呢?
7楼2015-11-02 08:34:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
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   回复  
五星好评  顶一下,感谢分享!
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 应届大龄土博教职选择,石河子大学和广西大学 10+3 茅山老道士 2024-06-10 17/850 2024-06-11 19:07 by rose380
[论文投稿] 给大家分享一下扎心的投稿经历 +10 lizhengke06 2024-06-08 12/600 2024-06-11 16:01 by licao0415
[考博] 还能读上博士吗 +9 小伙走起 2024-06-09 16/800 2024-06-11 15:36 by 小伙走起
[基金申请] 有没有机械的前辈分享一下评上海优都是什么成果啊 +3 wulala800 2024-06-10 3/150 2024-06-11 14:34 by andywei1028
[硕博家园] 博士找工作真难 +17 sdsk47 2024-06-05 22/1100 2024-06-11 13:36 by 心宁
[硕博家园] 求助 +6 LYWwrz 2024-06-09 9/450 2024-06-11 13:12 by powerhours
[考博] 博导推荐 +3 落雨吃鸡蛋 2024-06-07 5/250 2024-06-11 11:12 by 安塔瓦拉多
[食品] 食品博士导师 +3 小李醒yy 2024-06-11 4/200 2024-06-11 10:50 by 小李醒yy
[硕博家园] 迷茫的毕业季 +6 独孤老狗 2024-06-08 6/300 2024-06-11 10:43 by 浅唱,唱歌好听
[论文投稿] water research状态咨询 5+3 Flyyawa 2024-06-10 6/300 2024-06-11 09:45 by bobvan
[教师之家] 郑雅婷工作后被辞退无业中,于是报复硕导,证据不足失败后可能报复父母 +17 苏东坡二世 2024-06-05 22/1100 2024-06-10 23:40 by chong2024
[基金申请] 请问评审专家医学口今年函审是哪天结束的? +3 hyzs6688 2024-06-08 16/800 2024-06-10 10:14 by hyzs6688
[基金申请] B01化学口的,近几年要什么样的评价才能上会呀? +4 salmon95 2024-06-06 5/250 2024-06-09 09:45 by chem-cai
[论文投稿] Analytica Chimica Acta 投稿系统中 abstract 字符数限制600什么意思? 3+3 3921 2024-06-05 5/250 2024-06-09 09:45 by bobvan
[有机交流] 席夫碱反应ph敏感 控制ph在多少呢? 10+4 段sy 2024-06-05 5/250 2024-06-08 23:15 by 段sy
[基金申请] 5B能上会么? +12 Tenfu 2024-06-05 14/700 2024-06-06 22:34 by mch20120521
[论文投稿] under consideration是被拒了了吗? +3 alive.... 2024-06-06 3/150 2024-06-06 12:21 by 知识产权服务
[基金申请] E口YQ答辩 +3 firepick 2024-06-05 3/150 2024-06-05 17:12 by kaka198888
[基金申请] 各位,F01口的YQ有收到通知的吗? +3 zmm891258 2024-06-05 5/250 2024-06-05 17:06 by godlovewho
[论文投稿] Industrial lurication and tribology初次投稿被退回求助 5+3 lmjz66 2024-06-05 3/150 2024-06-05 16:32 by wjykycg
信息提示
请填处理意见