费米面计算,这里我们以Cu为例,见vaspkit/examples/Cu_fermi_surface
。
第一步,准备好FCC Cu
的POSCAR,注意一定是原胞(primitive cell )和用于静态计算的INCAR文件;
第二步,运行VASPKIT
,输入261
命令,产生用于计算计算费米面的KPOINTS
和POTCAR
文件,执行界面如下图所示;
------------>>
261
+-------------------------- Warm Tips --------------------------+
* Accuracy Levels: (1) Medium: 0.02~0.01;
(2) Fine: < 0.01;
* 0.01 is Generally Precise Enough!
+---------------------------------------------------------------+
Input KPT-Resolved Value (unit: 2*PI/Angstrom):
------------>>
0.008
-->> (1) Reading Structural Parameters from POSCAR File...
-->> (2) Written KPOINTS File!
第三步,提交VASP作业;
第四步,待VASP计算完成后,再次运行VASPKIT
,输入262
命令,生成包含费米面数据的FERMISURFACE.bxsf
文件,执行界面如下图所示;
------------>>
262
-->> (1) Reading Input Parameters From INCAR File...
-->> (2) Reading Structural Parameters from POSCAR File...
-->> (3) Reading Fermi-Level From DOSCAR File...
ooooooooo The Fermi Energy will be set to zero eV oooooooooooooooo
-->> (4) Reading Energy-Levels From EIGENVAL File...
-->> (5) Written FERMISURFACE.bxsf File!
第五步,运行XcrySDen
软件,按照下图步骤。我们发现第五条能带穿过费米面,因此我们选择Band Number 5
,然后点击Selected
,然后得到Cu的费米面。推荐使用FermiSurfer( https://fermisurfer.osdn.jp/ )软件可视化费米面。
对于半导体不存在费米面,如果想可视化导带或价带中某个能级处的等能面,新建FERMI_ENERGY.in文件然后在该文件第二行(第一行为注释行)中人为指定能级大小,注意该值一定要保证穿过价带或导带。接下来调用vaspkit-252或263命令提取可视化。
如果您使用VASPKIT,请记得引用哦!
V. Wang, N. Xu, J.-C. Liu, G. Tang, W.-T. Geng, VASPKIT: A User-Friendly Interface Facilitating High-Throughput Computing and Analysis Using VASP Code, Computer Physics Communications 267, 108033, (2021), https://doi.org/10.1016/j.cpc.2021.108033
欢迎关注VASPKIT公众号。
「感觉有帮助?一键投喂 牛奶/咖啡/冰阔乐!」
(๑>ڡ<)☆哇~太棒了!
使用微信扫描二维码完成支付