计算化学公社

标题: 计算孤立分子能量时偶极校正是否需要打开 [打印本页]

作者
Author: Sinddo 时间: 2021-2-17 09:31
标题: 计算孤立分子能量时偶极校正是否需要打开
计算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还是否需要打开呢？

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

作者
Author: Sinddo 时间: 2021-2-17 11:22

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

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

作者
Author: 欧拉角 时间: 2021-2-17 11:56

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

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

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

作者
Author: Sinddo 时间: 2021-2-17 11:59

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

前面提到的问 ...

明白了，谢谢~

作者
Author: 小阿龙 时间: 2024-1-5 21:08
请问什么是孤立分子，我在优化甲烷、水、甲醇这些分子算不算

欢迎光临计算化学公社 (http://bbs.keinsci.com/)