24小时热门版块排行榜    

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

我爱天津

金虫 (正式写手)


引用回帖:
13楼: Originally posted by guangrongbao at 2015-11-02 10:39:14
写的不错...

可是你只用castep不用vasp呀
15楼2015-11-02 12:45:16
已阅   回复此楼   关注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   回复  
五星好评  顶一下,感谢分享!
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 统计一下:硕士毕业答辩后的谢师宴是学生出钱,还是老师出钱? +18 苏东坡二世 2024-06-02 20/1000 2024-06-02 21:44 by 苏东坡二世
[论文投稿] Scientific reports 投稿 +4 chencome12 2024-06-02 4/200 2024-06-02 21:39 by lizhengke06
[微米和纳米] 银浆叠压烧结粘片的问题有没有好的解决方法 +3 微纳材料 2024-05-30 5/250 2024-06-02 21:17 by manmansboy
[教师之家] 职能部门工作人员态度不好是普遍的吗?怎么让他们态度好一些? +5 河西夜郎 2024-06-01 5/250 2024-06-02 19:19 by Quakerbird
[基金申请] +5 河马の史诗 2024-06-02 5/250 2024-06-02 19:05 by 328838485
[教师之家] 选择 +8 459582015 2024-05-31 11/550 2024-06-02 18:30 by auvauv
[基金申请] 化学B02口青基 代表作都是什么水平的?向大佬求助 +12 arthas_007 2024-06-01 14/700 2024-06-02 16:18 by 夏木荫浓
[基金申请] 讨论贴:好本子一定会被打招呼的淘汰吗 +40 ioe3937 2024-05-27 51/2550 2024-06-02 13:52 by Howard28
[基金申请] 九部门发文:不得将专利授权数量作为人才评价、项目评审、职称评定、高校评价等的条件 +17 sjtu2012 2024-05-28 21/1050 2024-06-02 13:43 by 欢乐颂叶蓁
[论文投稿] 编辑对文章提的意见 40+3 qa8386136 2024-06-01 3/150 2024-06-02 13:28 by bobvan
[考博] 导师不让硕转博,让我去国外读博,能理解吗? +12 萧山幽谷 2024-05-29 20/1000 2024-06-02 12:01 by yuan0806
[考博] 24or25材料专业申博 +4 农夫三拳有点痛 2024-05-30 11/550 2024-06-01 14:45 by Napoleonsky
[基金申请] 博后特别资助状态变化 +24 随梦而飞2017 2024-05-30 35/1750 2024-06-01 10:10 by 青岛阳仔
[考博] 申博求助?本硕双非一篇三区和四区去985工科非天坑专业有没有希望? +4 LYS1200 2024-05-29 6/300 2024-06-01 08:50 by 学术渣渣神
[教师之家] 中年 (金币+3) +18 459582015 2024-05-28 19/950 2024-06-01 00:41 by 沈婉婷.Girl
[论文投稿] 《控制理论与应用》期刊收版面费吗? +4 ygj2015 2024-05-29 4/200 2024-05-30 06:59 by nono2009
[论文投稿] 有没有老师需要发表论文 +4 金老师论文助理- 2024-05-29 4/200 2024-05-29 16:51 by liuyupu132
[论文投稿] 核心初审被拒,理由是“选题的意义不明确,文章写得不像是科技论文”,怎么改 5+3 工藤雷花樱 2024-05-27 8/400 2024-05-29 10:09 by topedit
[基金申请] E10开始送了,希望有好运 +5 sail 2024-05-27 5/250 2024-05-28 18:36 by 芝小芝
[基金申请] 面上基金会评专家,有回避机制吗? +4 huang1991js 2024-05-27 4/200 2024-05-27 19:08 by 星火12
信息提示
请填处理意见