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

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

建议加入额外参数：

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说的也很清楚哦。

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

用户类似问题：

问题1：按斯莱脱规则计算K.Cu.I的最外层电子感受到的有效...

太难了 >>详细

问题2：交流电有效值是不是平均值 计算通过电荷时为什么用...

交流电有效值是指方均根值（有效值），不是平均值。 计算通过电荷时，也是采用方均根值，不会是平均值，交流电在一个周期内，通过电荷的平均值是零。 >>详细

问题3：计算极化电荷形成的偶极矩

偶极矩是通过专门的实验测定的结果，没有具体的算法。我们可以通过分子的极性大小来比较不同分子的偶极矩。通过互溶程度也可以，相似相溶。 >>详细

问题4：有效核电荷数 计算问题

自己算算就好了 >>详细

问题5：计算电容的公式中的电荷C是自由电荷还是自由电荷与...

自由电荷 >>详细

【本站声明】:文章内容由网络用户投稿，内容仅供参考，不构成投资建议，也不代表本网赞同其观点和对其真实性负责。本网站不承担任何形式的责任。特此声明。如有侵权请联系QQ：100966298 QQ: 77100680 我们会在24小时内处理。