对于接触过磁性材料和强关联电子体系的计算的打工人来说，DFT+U 方法想必一定不会陌生，本方法中我们基于线性响应近似的方法计算自洽 Hubbard U 值。

具体也是将VASP官方的计算方法进行了集成，省去了对POTCAR、INCAR的修改，一键直接进行U值的计算并进行线性拟合。

具体操作如下：

1. 准备必要文件

这里必要的文件是：POSCAR (你也可以自己生成KPOINTS、POTCAR、INCAR文件，也可以通过本脚本自动生成，注意自动生成INCAR后需要手动设置其中的MAGMOM值)

2. 在程序中输入本机运行VASP的代码

这一步是必须的，运行脚本后输入0，0后，输入服务器运行VASP的代码。

随后即生成run.wsy文件，如果不想通过程序生成，也可事先自己写好run.wsy文件。

3. 计算DFT基态

运行程序后后输入1，这里会要求你输入需要计算U值原子在POSCAR中的位置以及U值加入的轨道。

之后程序自动创建文件，调用vaspkit生成新的POTCAR文件，并调VASP进行计算。

(此时会生成input.wsy文件，不要修改)

4. 进行+U的自洽与非自洽计算

5. 运行脚本后输入2，程序自动调VASP进行自洽与非自洽计算。

5. 计算U值

上述计算全部完成之后，运行脚本后输入3，计算U值，并生成output.wsy文件。

打开output.wsy，最后一行即为该原子所对应的U值。

也可如VASP网中通过上述的第一列和最后两列数据进行线性拟合得到。

(由于赝势选取的不同，计算结果与官网有些许差别)

初代版本，还请多多支持，后续还会加入新功能。

最新程序下载地址：https://github.com/Code-WSY/Code-WSY