计算化学公社

标题: vsapsol溶剂化下电子步不收敛 [打印本页]

作者
Author: xmuqyj 时间: 2023-12-20 10:51
标题: vsapsol溶剂化下电子步不收敛
我在优一个单层的C3N4结构，结构优化是很快收敛的了，但是加上vaspsol参数就不能收敛了。奇怪的是C3N4负载上金属原子vaspsol又能收敛了。。。请问是怎么回事呢？

INCAR
System = C3N4

ISMEAR = 0
SIGMA = 0.05
ALGO = FAST
ENCUT = 450
EDIFF = 1E-5

#ISPIN = 2
#MAGMOM = 0
#LORBIT = 11
LDIPOL = .TRUE.
IDIPOL = 3

NWRITE = 0
#LWAVE = .FALSE.
LCHARG = .FALSE.
LVTOT =  .FALSE.
LVHAR =  .FALSE.
LELF =  .FALSE.

IBRION = 2
POTIM = 0.15
EDIFFG = -0.01
NSW = 500
NCORE = 8
#ISIF = 3
IVDW = 12

# Solvation
LSOL = .TRUE.
EB_K = 78.4

作者
Author: 白皛白 时间: 2023-12-20 10:59
因为看起来还是有收敛的趋势的，所以可以通过调整NELM参数来多跑几步看看能否收敛。
或者算一个不加溶剂化的单点写一个WAVECAR，然后加溶剂化计算的时候读取这个WAVECAR，一般也会加快SCF收敛。

作者
Author: 123wo 时间: 2023-12-20 14:37
可以给出完整的OSZICAR吗，C3N4和负载了单原子的都要

作者
Author: 卡开发发 时间: 2023-12-20 14:44
如下经验上的建议：
1、ALGO换成N。
2、LDIPOL=T和IDIPOL=3。
这两个会导致难收敛，额外适当调整AMIX和BMIX（通常是降低减小震荡）。如果不行真空下面先做个波函数然后读取进去。额外可能需要考虑PREC=A或适当提升ENAUG，PB方程对这个会有点敏感。

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