计算化学公社

标题: 求助：SCF在加dipole correction后不收敛 [打印本页]

作者
Author: diffcwmol 时间: 2022-2-23 21:33
标题: 求助：SCF在加dipole correction后不收敛
请教各位同仁：

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

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

谢谢。

作者
Author: 卡开发发 时间: 2022-2-24 18:18
本帖最后由卡开发发于 2022-2-24 18:24 编辑

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

作者
Author: diffcwmol 时间: 2022-2-26 11:57

卡开发发发表于 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。

作者
Author: 卡开发发 时间: 2022-2-26 15:13

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

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

作者
Author: diffcwmol 时间: 2022-3-5 04:47

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

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

作者
Author: 卡开发发 时间: 2022-3-5 06:16
本帖最后由卡开发发于 2022-3-5 06:29 编辑

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

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

作者
Author: onemoretry 时间: 2022-3-9 18:59
发个POSCAR看一下，有可能是模型在cell中的位置导致的。

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