计算化学公社

标题: 求助：VASP结构优化Fe7S8/C3N4异质结构不收敛 [打印本页]

作者
Author: jzj688 时间: 2023-12-15 12:32
标题: 求助：VASP结构优化Fe7S8/C3N4异质结构不收敛
本帖最后由 jzj688 于 2023-12-15 12:33 编辑

本人用MS进行建模，异质结构两部分（上Fe7S8（2 0 0），下C3N4（0 0 2）），后在VASP中进行结构优化，优化了三百多离子步（中间停过，后续算），能量一直下降，最大力在远离我设定的收敛精度（0.05）震荡，一直不收敛。目前优化的结构中下面的C3N4（0 0 2）和初始结构相差较大，一部分片段已经向上弯和Fe7S8（2 0 0）原子吸引成键。
目前，有几个问题想请教各位老师：
1、目前优化的结构中下面的C3N4（0 0 2）和初始结构相差较大，我所构建的初始和目前结构是否合理？
2、优化过程中，能量一直下降，最大力在远离我设定的收敛精度震荡，是什么原因造成的？如何解决？（我之前也尝试调整过一些参数，效果不理想）

INCAR:
ISTART =  1
ISPIN  =  1
ICHARG =  1
LREAL  = Auto
LWAVE  = .TRUE.
LCHARG = .TRUE.
ENCUT  =  450
ALGO = Fast
IVDW = 11
NCORE = 4
KPAR = 2

ISMEAR =  0
SIGMA  =  0.2
NELM =  200
EDIFF  =  1E-05

NSW =  300
IBRION =  2
ISIF =  2
EDIFFG = -0.05

KPOINTS:K-Spacing Value to Generate K-Mesh: 0.040
0
Monkhorst-Pack
2 2 1
0.0  0.0  0.0

初始结构：
(, 下载次数 Times of downloads: 17)
(, 下载次数 Times of downloads: 17)

优化三百多离子步后的结构：
(, 下载次数 Times of downloads: 18)
(, 下载次数 Times of downloads: 17)

最近一百多步能量变化情况：
(, 下载次数 Times of downloads: 16)
最近一百多步最大力变化情况：
(, 下载次数 Times of downloads: 15)
最近几步最大力变化情况：
(, 下载次数 Times of downloads: 15)

作者
Author: Deepast 时间: 2023-12-15 17:44
快收敛了EDIFF可以调小一点
另外你这个失配率多少，可能应力太大导致变形了

作者
Author: jzj688 时间: 2023-12-15 17:55

Deepast 发表于 2023-12-15 17:44
快收敛了EDIFF可以调小一点
另外你这个失配率多少，可能应力太大导致变形了

从图上看最大力还差0.3左右收敛，我试一试调小一点
失配率5%多一点，MS建模build layer时没有提示

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