求助：SCF在加dipole correction后不收敛

diffcwmol

请教各位同仁：

我计算的是二维材料slab模型，原子数较多（~45个）。
由于上下表面静电势能有差异需要使用dipole correction，但是加了dipole correction之后SCF就不收敛了，dE一直居高不下。
但是我取消dipole correction之后，SCF很快就收敛了。

问题：我必须要使用dipole correction计算静电势能，有啥方法可以让SCF收敛吗？

谢谢。

卡开发发

偶极修正的情况建议可以试一下读取电荷密度或波函数，如ICHARG=11。

diffcwmol

卡开发发 发表于 2022-2-24 18:18
偶极修正的情况建议可以试一下读取电荷密度或波函数，如ICHARG=11。

谢谢。
不用dipole correction跑完，得到WAVECAR和CHGCAR；然后加上dipole correction，设置ICHARG=11，显示CHGCAR无法读取，ERROR: charge density could not be read from file CHGCAR for ICHARG>10。

卡开发发

diffcwmol 发表于 2022-2-26 11:57
谢谢。
不用dipole correction跑完，得到WAVECAR和CHGCAR；然后加上dipole correction，设置ICHARG=11， ...

设置ICHARG=11前重新自洽一下看看呢？

diffcwmol

卡开发发 发表于 2022-2-26 15:13
设置ICHARG=11前重新自洽一下看看呢？

谢谢卡老师的建议，这个方法之前试过了，自洽还是不能收敛。我猜估计是上下表面真空静电势差别太大了，大概几个eV，您觉得有没有这种可能？

卡开发发

diffcwmol 发表于 2022-3-5 04:47
谢谢卡老师的建议，这个方法之前试过了，自洽还是不能收敛。我猜估计是上下表面真空静电势差别太大了，大 ...

你可能会错意了，是在不考虑偶极矩修正的情况进行重新自洽得到电子密度和波函数，你可以试着读这个波函数来做偶极矩修正的计算，如果还不行就把电荷密度也读进来。

onemoretry

发个POSCAR看一下，有可能是模型在cell中的位置导致的。