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

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

好久不见了，这一节小编给大家介绍如何计算电极化，born有效电荷。

参照manual，这里使用NaF作为例子

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

建议加入额外参数：

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

第二步，静态计算，也就是自洽计算，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

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

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

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

第三步，做电极化计算，要分三个方向做，也就是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

注意上面INCAR在最后加入了

LBERRY = TRUE

IGPAR = 1

NPPSTR = 2

DIPOL = 0.2 0.2 0.2

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

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

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

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

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

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

manual说的也很清楚哦。

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