vaspkit计算三维费米面

Noblemon · 发表于 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。希望各位不吝指教。

vaspkit · 发表于 2021-12-6 16:50:57

必须是金属体系。为了减小输出文件的占用空间，所以只输出了穿过费米能级的的能带，未穿过费米能的能带即时输出也不会显示，因此不需要输出全部能带。

Noblemon · 发表于 2021-12-6 19:16:11

vaspkit 发表于 2021-12-6 16:50
必须是金属体系。为了减小输出文件的占用空间，所以只输出了穿过费米能级的的能带，未穿过费米能的能带即时 ...

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

vaspkit · 发表于 2021-12-6 22:18:08

Noblemon 发表于 2021-12-6 19:16
谢谢您的解答。但是我的材料体系是半导体，想画电子和空穴的pocket。有没有可能通过修改某个脚本实现输出 ...

变通一下，在FERMI_ENERGY.in文件中修改第二行数字从而手动平移费米能到所关心的能量位置，然后再次调用vaspkit提取费米面。

Noblemon · 发表于 2021-12-7 06:59:11

vaspkit 发表于 2021-12-6 22:18
变通一下，在FERMI_ENERGY.in文件中修改第二行数字从而手动平移费米能到所关心的能量位置，然后再次调用va ...

好的，谢谢指点

dlbao · 发表于 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。所以我想确认一下。谢谢！

vaspkit · 发表于 2022-4-9 07:41:15

必须通过修改FERMI_ENERGY.in文件中第二行数字手动平移费米能，因为VASPKIT只输出穿过费米能的能带。

米墨 · 发表于 2022-5-25 00:24:05

1.3.3 版本输出没有输出 FERMI_ENERGY.in，应该如何得到.in 文件呢

vaspkit · 发表于 2022-5-25 16:17:15

运行提取态密度等功能可得到FERMI_ENERGY文件，重新命名为FERMI_ENERGY.in，手动修改第二行数字调整费米能级。

米墨 · 发表于 2022-10-10 22:49:24

大神们我想请教一下，我的体系中有镧系元素，但是体系是半导体，算出的bandgap也在1.5左右，通过调节FERMI_ENERGY.in的fermi能级，绘制等能面，但是发现为什么将fermi能级调到bandgap中间还会有等能面出现呢，理论上这里并没有能带穿过啊，我试了其他体系，穿过bandgap时就没有等能面出现，这是什么原因呢？

		自动登录	找回密码
密码			立即注册