计算化学公社

标题: 求助，想计算表面能但是完全无法收敛 vasp [打印本页]

作者
Author: ZizhenZhou 时间: 2020-11-5 11:58
标题: 求助，想计算表面能但是完全无法收敛 vasp
本帖最后由 ZizhenZhou 于 2020-11-5 12:03 编辑

如题，我对 LiMnO2 conventional cell 切 001 表面。想去进行优化然后计算能量，vasp第一个电子步完全不收敛，能量乱飞。
INCAR 中加入了 IDIPOL， LDIPOL。我有考虑表面可能极化十分严重，可是不知道有什么办法可以解决。求各位指点一二。

作者
Author: wzkchem5 时间: 2020-11-5 13:41
表面极化是一方面，另外锰酸盐是强相关体系，SCF本来就不好收敛。
你加U了吗？

作者
Author: ZizhenZhou 时间: 2020-11-6 11:28

wzkchem5 发表于 2020-11-5 13:41
表面极化是一方面，另外锰酸盐是强相关体系，SCF本来就不好收敛。
你加U了吗？

您好，我有加U和磁矩，下面是我的INCAR，麻烦您看一下。
EDIFFG = -0.03
NSW = 200
NELM = 300

LDIPOL = .True.
IDIPOL = 3
LWAVE = .FALSE.
LCHARG = .FALSE.

ISYM = 0
ISPIN = 2
MAGMOM = 12*0.0 12*2.0 24*0.0
LORBIT = 11

LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = -1 2 -1
LDAUU = 0.0 3.9 0.0
LMAXMIX = 4
AMIX = 0.1
BMIX = 0.0001
AMIX_MAG = 0.8
BMIX_MAG = 0.0001

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