计算孤立分子能量时偶极校正是否需要打开

Sinddo

计算Slab吸附一个冠醚分子，（1）在计算孤立分子能量时，是否需要进行偶极校正？论坛里有说需要的，而vasp官方的教程里没有进行校正，偶极校正是否对能量影响不大？ （2）如果需要校正的话，在INCAR文件里IDIPOL = 4即各个方向进行校正，而LDIPOL是否需要打开？我打开LDIPOL电子步没能成功收敛，而只开IDIPOL能够顺利结束，在结果中也有偶极矫正的结果：
DIPCOR: dipole corrections for dipol
direction  1 min pos   228, direction  2 min pos   222, direction  3 min pos   226,
dipolmoment          -0.001618     -0.019035      0.001084 electrons x Angstroem
Tr[quadrupol]        51.752379

energy correction for charged system         0.000000 eV
dipol+quadrupol energy correction            0.000001 eV
added-field ion interaction          0.000000 eV  (added to PSCEN)
可以看出对能量的影响很小，那么LDIPOL还是否需要打开呢？

欧拉角

如果分子极性不强影响会比较小，偶极矩比较大的分子就得考虑修正了。LDIPOL打开的情况下自洽场确实不容易收敛，您可以试着读取不加LDIPOL时候的电荷密度进行校正的计算。

Sinddo

欧拉角 发表于 2021-2-17 10:14
如果分子极性不强影响会比较小，偶极矩比较大的分子就得考虑修正了。LDIPOL打开的情况下自洽场确实不容易收 ...

感谢指导，想问下如果我只需要能量数据的话，是否程序输出E0已经校正了能量？您所说的读取电荷密度校正计算是指计算什么数据，是否与能量有关？

欧拉角

Sinddo 发表于 2021-2-17 11:22
感谢指导，想问下如果我只需要能量数据的话，是否程序输出E0已经校正了能量？您所说的读取电荷密度校正计 ...

已经校正了，但是是仅仅校正了能量，开LDIPOL则是对电荷密度、静电势等还会再进行修正。

前面提到的问题我是针对开LDIPOL情况下自洽场难收敛给的建议，您可以在不开LDIPOL的情况下先进行自洽计算达到收敛，然后读取波函数和电荷密度并开启LDIPOL=T进行后续计算。

Sinddo

欧拉角 发表于 2021-2-17 11:56
已经校正了，但是是仅仅校正了能量，开LDIPOL则是对电荷密度、静电势等还会再进行修正。

前面提到的问 ...

明白了，谢谢~

小阿龙

请问什么是孤立分子，我在优化甲烷、水、甲醇这些分子算不算