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

1、第一步优化，正常优化结构，注意精度要高一些输入文件如下图建议加入额外参数：ISTART=0ICHARG=2IBRION =2NSW =40NELM=40EDIFF=1E-6EDIFFG=-0.005ISMEAR=0SIGMA=0.05ISIF=3ENCUT=500PREC=AccurateLREAL=AutoLCHARG=.TURE.NPAR=4

2、第二步，静态计算，也就是自洽计算，INCAR如下图：ISTART=0ICHARG=2IBRION =-1NSW =0NELM=40EDIFF=1E-6EDIFFG=-0.005ISMEAR=0SIGMA=0.05ISIF=2ENCUT=500PREC=AccurateLREAL=AutoLCHARG=.TURE.NPAR=4

3、建立三个文件夹x, y, z把静态的INCAR，POSCAR，POTCAR，KPOINTS，POTCAR，CHGCAR,WAVECAR分别复制到文件夹x, y, z中。

4、第三步，做电极化计算，要分三个方向做，也就是3个方向分别做一次静态计算，INCAR参数：system=RhIBRION =-1NSW =0!POTIM=0.1NELM=40EDIFF=1E-6EDIFFG=-0.005ISMEAR=0SIGMA=0.05ISIF=2ENCUT=500PREC=AccurateLMAXMIX=4LREAL=AutoLCHARG=.TURE.NPAR=4LBERRY = TRUEIGPAR = 1NPPSTR = 2DIPOL = 0.2 0.2 0.2

5、注意上面INCAR在最后加入了LBERRY = TRUEIGPAR = 1NPPSTR = 2DIPOL = 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、先到这里吧，下一节小编写个脚本放上去，小伙伴就能直接处理数据了