24小时热门版块排行榜    

查看: 22039  |  回复: 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的回帖

a159951zhao

新虫 (小有名气)


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

bandkit.x运行的语句是什么啊 lunix小白求教
270楼2017-12-27 14:31:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 我爱天津 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[微米和纳米] 银浆叠压烧结粘片的问题有没有好的解决方法 +3 微纳材料 2024-05-30 5/250 2024-06-02 21:17 by manmansboy
[教师之家] 选择 +8 459582015 2024-05-31 11/550 2024-06-02 18:30 by auvauv
[基金申请] 离职后国自然项目不能变更单位,在新单位还能申请新的国自然项目吗 5+3 宜兰 2024-05-31 10/500 2024-06-02 17:50 by 月五年子庚
[基金申请] 化学口B0109(高分子合成),拿青年基金一般需要怎样的文章水平? +14 salmon95 2024-05-30 23/1150 2024-06-02 16:37 by 款款飞
[基金申请] 化学B02口青基 代表作都是什么水平的?向大佬求助 +12 arthas_007 2024-06-01 14/700 2024-06-02 16:18 by 夏木荫浓
[基金申请] 为什么美国要封锁华为而不封锁丘C同呢? +6 zju2000 2024-06-02 7/350 2024-06-02 16:12 by sxjc500
[基金申请] 九部门发文:不得将专利授权数量作为人才评价、项目评审、职称评定、高校评价等的条件 +17 sjtu2012 2024-05-28 21/1050 2024-06-02 13:43 by 欢乐颂叶蓁
[考博] 导师不让硕转博,让我去国外读博,能理解吗? +12 萧山幽谷 2024-05-29 20/1000 2024-06-02 12:01 by yuan0806
[考博] 申请2024或2025年博士研究生 +5 嘟噜嘟1 2024-05-29 12/600 2024-06-01 22:36 by 嘟噜嘟1
[教师之家] 研究生在毕业答辩时挂了,遗憾 +10 spacr2022 2024-05-27 13/650 2024-06-01 17:45 by keyaner23
[论文投稿] 求Sci期刊推荐 10+4 甄小鱼 2024-05-30 7/350 2024-06-01 10:41 by bobvan
[基金申请] 数理的人才答辩通知发了吗 +7 yzy3327 2024-05-30 7/350 2024-06-01 09:15 by 撒哈拉西部
[硕博家园] 求问endnote上Wiley引用格式在哪下载 +3 Elon. 2024-05-30 3/150 2024-05-31 10:13 by crazy peng
[文学芳草园] 物是人非 +4 myrtle 2024-05-30 4/200 2024-05-30 15:05 by mapenggao
[论文投稿] EI学报,一审返修后,为啥不再送审,直接终审中? +4 qweasd12345 2024-05-27 6/300 2024-05-29 00:02 by dut_ameng
[基金申请] E05青基有几个评审 +4 KYXY123 2024-05-28 4/200 2024-05-28 19:25 by popt2t
[论文投稿] EI期刊审稿人邮箱问题 5+3 shier妈妈 2024-05-27 4/200 2024-05-28 14:53 by topedit
[有机交流] 奇怪的物质 100+4 桃桃PXS 2024-05-27 7/350 2024-05-28 10:22 by 091602
[硕博家园] 课题做到一半担心做不出来怎么办 +4 mkpooooo 2024-05-27 7/350 2024-05-27 11:56 by kanmiaolucky
[硕博家园] 周日 +6 1加油哦棒 2024-05-26 9/450 2024-05-27 10:30 by hahamyid
信息提示
请填处理意见