对于接触过磁性材料和强关联电子体系的计算的打工人来说,DFT+U 方法想必一定不会陌生,本方法中我们基于线性响应近似的方法计算自洽 Hubbard U 值。
具体也是将VASP官方的计算方法进行了集成,省去了对POTCAR、INCAR的修改,一键直接进行U值的计算并进行线性拟合。
具体操作如下:
- 准备必要文件
这里必要的文件是:POSCAR (你也可以自己生成KPOINTS、POTCAR、INCAR文件,也可以通过本脚本自动生成,注意自动生成INCAR后需要手动设置其中的MAGMOM值)
- 在程序中输入本机运行VASP的代码
这一步是必须的,运行脚本后输入0,0后,输入服务器运行VASP的代码。
随后即生成run.wsy文件,如果不想通过程序生成,也可事先自己写好run.wsy文件。
运行程序后后输入1,这里会要求你输入需要计算U值原子在POSCAR中的位置以及U值加入的轨道。
之后程序自动创建文件,调用vaspkit生成新的POTCAR文件,并调VASP进行计算。
(此时会生成input.wsy文件,不要修改)
-
进行+U的自洽与非自洽计算
-
运行脚本后输入2,程序自动调VASP进行自洽与非自洽计算。
5. 计算U值
上述计算全部完成之后,运行脚本后输入3,计算U值,并生成output.wsy文件。
打开output.wsy,最后一行即为该原子所对应的U值。
也可如VASP网中通过上述的第一列和最后两列数据进行线性拟合得到。
(由于赝势选取的不同,计算结果与官网有些许差别)
初代版本,还请多多支持,后续还会加入新功能。
最新程序下载地址:https://github.com/Code-WSY/Code-WSY