计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1172|回复 Reply: 2
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] 求助:VASP结构优化Fe7S8/C3N4异质结构不收敛

[复制链接 Copy URL]

4

帖子

0

威望

107

eV
积分
111

Level 2 能力者

本帖最后由 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

初始结构:



优化三百多离子步后的结构:



最近一百多步能量变化情况:

最近一百多步最大力变化情况:

最近几步最大力变化情况:




88

帖子

0

威望

1346

eV
积分
1434

Level 4 (黑子)

2#
发表于 Post on 2023-12-15 17:44:28 | 只看该作者 Only view this author
快收敛了EDIFF可以调小一点
另外你这个失配率多少,可能应力太大导致变形了

4

帖子

0

威望

107

eV
积分
111

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2023-12-15 17:55:32 | 只看该作者 Only view this author
Deepast 发表于 2023-12-15 17:44
快收敛了EDIFF可以调小一点
另外你这个失配率多少,可能应力太大导致变形了

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-25 15:50 , Processed in 0.289563 second(s), 29 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list