关于表面吸附分子的偶极校正问题

从头学VASP

如题，本人最近遇到了slab吸附小分子的偶极校正问题，查了网上很多的内容也没有一个比较确切的说法。所以在这里发个贴记录一下，同时希望有会的同学指导一下。
对于单纯的slab，因为固定一半放开一半的缘故，加上了偶极校正进行结构优化，但是只加LDIPOL=T IDIPOL=3得话就是死活不收敛，最后又加上了IDIPOL=0.5 0.5 0.5，最后功函数画出来是一个两边水平的阶梯。
对于吸附分子的slab，我一开始是把分子吸附到slab上面，然后INCAR加上LDIPOL=T IDIPOL=3，也还是死活不收敛。看了网上说，需要先进行结构优化和自洽，优化好的结构再加上偶极校正进行结构优化和自洽，不知道这种情况只加LDIPOL=T IDIPOL=3能不能收敛，如果不能收敛的话，是不是还得要加上IDIPOL=0.5 0.5 0.5。
另外，关于IDIPOL= xyz的设置，网上有说设置几何中心的，那么对于slab吸附分子的情况来说只需要算一下z方向的平均值就可以了吗。这个平均值应该是相对于优化好的结构来说的吗？

从头学VASP

对了 补充一下 用的是Vasp

sobereva

从头学VASP 发表于 2023-10-11 21:25
对了 补充一下 用的是Vasp

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。