计算化学公社

标题: 含过渡金属体系，能量一直震荡，不收敛 [打印本页]

作者
Author: huo777 时间: 6 hour ago
标题: 含过渡金属体系，能量一直震荡，不收敛
本帖最后由 huo777 于 2026-7-27 11:13 编辑

我的体系含Re，已经根据社长的推文改过一次参数，但是还是很震荡，不收敛
参数如下，原子顺序是Si Al O Re C
SYSTEM = Re/ZSM5
ISTART = 0
ICHARG = 2
SYSTEM = Re/ZSM5
ISTART = 0
ICHARG = 2
PREC = Normal
ENCUT = 520
ISMEAR = 0
SIGMA = 0.05
ISPIN = 2
ISYM = 0
LREAL = .FALSE.
ALGO = Fast
NELM = 200
NELMIN = 8
LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = -1 -1 -1 -2 -1
LDAUU = 0.0 0.0 0.0 4.0 0.0
LDAUJ = 0.0 0.0 0.0 0.0 0.0
AMIX = 0.2
BMIX = 0.0001
AMIX_MAG = 0.8
BMIX_MAG = 0.0001
LMAXMIX = 6
LASPH = .TRUE.
IBRION = 2
NSW = 400
ISIF = 2
POTIM = 0.1
EDIFF = 1E-6
EDIFFG = -0.02
IVDW = 11
LWAVE = .FALSE.
LCHARG = .FALSE.

作者
Author: xiehua 时间: 3 hour ago
用的是Re_sv_GW赝势吗，Re的另外两个赝势做Δ测试，结果很差，如果不是，可以换Re_sv_GW试试

作者
Author: huo777 时间: 2 hour ago

xiehua 发表于 2026-7-27 13:42
用的是Re_sv_GW赝势吗，Re的另外两个赝势做Δ测试，结果很差，如果不是，可以换Re_sv_GW试试

不是，用的就是Re，我用这个试试，谢谢！

作者
Author: xiehua 时间: 2 hour ago

huo777 发表于 2026-7-27 14:51
不是，用的就是Re，我用这个试试，谢谢！

参数：LMAXMIX = 6
GGA_COMPAT = .FALSE.记得设置，我做PAW赝势的时候，一些元素没设置这两个参数，也遇到出现不收敛和能量发散的情况，设置了就好了，希望对你有帮助

作者
Author: huo777 时间: 2 hour ago

xiehua 发表于 2026-7-27 15:01
参数：LMAXMIX = 6
GGA_COMPAT = .FALSE.记得设置，我做PAW赝势的时候，一些元素没设置这两个参数，也遇 ...

嗯嗯谢谢！

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