VASPKIT论坛

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
热搜: VASPKIT
查看: 4715|回复: 9

vaspkit计算三维费米面

[复制链接]

1

主题

4

帖子

19

积分

新手上路

Rank: 1

积分
19
发表于 2021-12-6 14:32:12 | 显示全部楼层 |阅读模式
求助各位老师:
我之前想用vasp+vaspkit(261,262)计算我的材料的三维费米面,但是发现得出的bxsf文件是空的,用xcrysden打开也会报错。vaspkit 1.2.5 和 1.3.0均有相同问题。
于是用example里的poscar,kpoints和incar重新跑了一下fcc铜的计算,发现我的计算结果和例子中得到的结果仍然不一样。比如例子中的bxsf用xcrysden打开会显示25条band(附图1),但是我得到的bxsf文件只有1条band(附图2)。个人感觉是我计算的bxsf只输出了穿过费米面的那一条band,但是不知道怎么才能显示全部的band。希望各位不吝指教。


1

1

2

2

vaspkit_26_question.zip

702.78 KB, 下载次数: 0

回复

使用道具 举报

35

主题

593

帖子

3430

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
3430
QQ
发表于 2021-12-6 16:50:57 | 显示全部楼层
必须是金属体系。为了减小输出文件的占用空间,所以只输出了穿过费米能级的的能带,未穿过费米能的能带即时输出也不会显示,因此不需要输出全部能带。
若使用VASPKIT软件请引用文献 Computer Physics Communications 267, 108033, (2021) ,您的支持是对我们付出的最大肯定。
VASPKIT英文网址https://vaspkit.com, 中文网址http://vaspkit.cn, 中文论坛http://bbs.vaspkit.cn, QQ学术交流群:331895604和498683602, 微信公众号:VASPKIT
回复

使用道具 举报

1

主题

4

帖子

19

积分

新手上路

Rank: 1

积分
19
 楼主| 发表于 2021-12-6 19:16:11 | 显示全部楼层
vaspkit 发表于 2021-12-6 16:50
必须是金属体系。为了减小输出文件的占用空间,所以只输出了穿过费米能级的的能带,未穿过费米能的能带即时 ...

谢谢您的解答。但是我的材料体系是半导体,想画电子和空穴的pocket。有没有可能通过修改某个脚本实现输出全部band的功能呢?
回复

使用道具 举报

35

主题

593

帖子

3430

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
3430
QQ
发表于 2021-12-6 22:18:08 | 显示全部楼层
Noblemon 发表于 2021-12-6 19:16
谢谢您的解答。但是我的材料体系是半导体,想画电子和空穴的pocket。有没有可能通过修改某个脚本实现输出 ...

变通一下,在FERMI_ENERGY.in文件中修改第二行数字从而手动平移费米能到所关心的能量位置,然后再次调用vaspkit提取费米面。
若使用VASPKIT软件请引用文献 Computer Physics Communications 267, 108033, (2021) ,您的支持是对我们付出的最大肯定。
VASPKIT英文网址https://vaspkit.com, 中文网址http://vaspkit.cn, 中文论坛http://bbs.vaspkit.cn, QQ学术交流群:331895604和498683602, 微信公众号:VASPKIT
回复

使用道具 举报

1

主题

4

帖子

19

积分

新手上路

Rank: 1

积分
19
 楼主| 发表于 2021-12-7 06:59:11 | 显示全部楼层
vaspkit 发表于 2021-12-6 22:18
变通一下,在FERMI_ENERGY.in文件中修改第二行数字从而手动平移费米能到所关心的能量位置,然后再次调用va ...

好的,谢谢指点
回复

使用道具 举报

0

主题

2

帖子

8

积分

新手上路

Rank: 1

积分
8
发表于 2022-4-9 06:06:00 | 显示全部楼层
本帖最后由 dlbao 于 2022-4-9 06:09 编辑
vaspkit 发表于 2021-12-6 22:18
变通一下,在FERMI_ENERGY.in文件中修改第二行数字从而手动平移费米能到所关心的能量位置,然后再次调用va ...

请问一下,费米面移动的能量单位是什么?我理解vaspkit根据vasp计算结果读取的数值应该是eV为单位的。但是我用FermiSurfer作图,并尝试在其中移动费米能时,发现结果跟我的能带结果不一致。我阅读FermiSurfer的手册时注意到,他们说在菜单里设置Fermi energy shift的单位是 Ry。所以我想确认一下。谢谢!
Ry.png

回复

使用道具 举报

35

主题

593

帖子

3430

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
3430
QQ
发表于 2022-4-9 07:41:15 | 显示全部楼层
必须通过修改FERMI_ENERGY.in文件中第二行数字手动平移费米能,因为VASPKIT只输出穿过费米能的能带。
若使用VASPKIT软件请引用文献 Computer Physics Communications 267, 108033, (2021) ,您的支持是对我们付出的最大肯定。
VASPKIT英文网址https://vaspkit.com, 中文网址http://vaspkit.cn, 中文论坛http://bbs.vaspkit.cn, QQ学术交流群:331895604和498683602, 微信公众号:VASPKIT
回复

使用道具 举报

0

主题

2

帖子

10

积分

新手上路

Rank: 1

积分
10
发表于 2022-5-25 00:24:05 | 显示全部楼层
1.3.3 版本输出 没有 输出 FERMI_ENERGY.in,应该如何得到.in 文件呢
回复

使用道具 举报

35

主题

593

帖子

3430

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
3430
QQ
发表于 2022-5-25 16:17:15 | 显示全部楼层
运行提取态密度等功能可得到FERMI_ENERGY文件,重新命名为FERMI_ENERGY.in,手动修改第二行数字调整费米能级。
若使用VASPKIT软件请引用文献 Computer Physics Communications 267, 108033, (2021) ,您的支持是对我们付出的最大肯定。
VASPKIT英文网址https://vaspkit.com, 中文网址http://vaspkit.cn, 中文论坛http://bbs.vaspkit.cn, QQ学术交流群:331895604和498683602, 微信公众号:VASPKIT
回复

使用道具 举报

0

主题

2

帖子

10

积分

新手上路

Rank: 1

积分
10
发表于 2022-10-10 22:49:24 | 显示全部楼层
大神们我想请教一下,我的体系中有镧系元素,但是体系是半导体,算出的bandgap也在1.5左右,通过调节FERMI_ENERGY.in的fermi能级,绘制等能面,但是发现为什么将fermi能级调到bandgap中间还会有等能面出现呢,理论上这里并没有能带穿过啊 ,我试了其他体系,穿过bandgap时就没有等能面出现,这是什么原因呢?
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|Archiver|手机版|小黑屋|VASPKIT论坛 ( 陕ICP备2020017291号 )

GMT+8, 2024-3-29 07:39 , Processed in 0.109265 second(s), 25 queries .

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表