vasp入门到精通[21]计算电极化，Born有效电荷

1、第一步优化，正常优化结构，注意精度要高一些输入文件如下图

建议加入额外参数：

ISTART=0

ICHARG=2

IBRION =2

NSW =40

NELM=40

EDIFF=1E-6

EDIFFG=-0.005

ISMEAR=0

SIGMA=0.05

ISIF=3

ENCUT=500

PREC=Accurate

LREAL=Auto

LCHARG=.TURE.

NPAR=4

2、第二步，静态计算，也就是自洽计算，INCAR如下图：

ISTART=0

ICHARG=2

IBRION =-1

NSW =0

NELM=40

EDIFF=1E-6

EDIFFG=-0.005

ISMEAR=0

SIGMA=0.05

ISIF=2

ENCUT=500

PREC=Accurate

LREAL=Auto

LCHARG=.TURE.

NPAR=4

3、建立三个文件夹x,  y,   z

把静态的INCAR，POSCAR，POTCAR，KPOINTS，POTCAR，CHGCAR,

WAVECAR分别复制到文件夹x, y, z中。

4、第三步，做电极化计算，要分三个方向做，也就是3个方向分别做一次静态计算，INCAR参数：

system=Rh

IBRION =-1

NSW =0

!POTIM=0.1

NELM=40

EDIFF=1E-6

EDIFFG=-0.005

ISMEAR=0

SIGMA=0.05

ISIF=2

ENCUT=500

PREC=Accurate

LMAXMIX=4

LREAL=Auto

LCHARG=.TURE.

NPAR=4

LBERRY = TRUE

IGPAR = 1

NPPSTR = 2

DIPOL = 0.2 0.2 0.2

5、注意上面INCAR在最后加入了

LBERRY = TRUE

IGPAR = 1

NPPSTR = 2

DIPOL = 0.2 0.2 0.2

其中第一次IGPAR = 1，第二次IGPAR = 2，第三次IGPAR = 3分别在x, y, z 文件夹中去做计算，将得到的OUTCAR里面的结果进行处理

6、因为时间关系小编就不解释怎么处理了，manual里写的挺清楚的

7、如上图，计算后OUTCAR会出现三个项：  ev， bp， ion。小编不做多解释了，希望小伙伴自行百度，了解一下基础知识。

8、把第二步的静态计算的POSCAR坐标文件改动一下原坐标，让他进行一点偏移

9、然后把第二步到第九步重新做一遍，也就是说做一遍偏移坐标后的静态计算，极化计算。

10、处理数据得到born有效电荷：

manual说的也很清楚哦。

11、先到这里吧，下一节小编写个脚本放上去，小伙伴就能直接处理数据了