求助：用cp2k计算MoS2时使用CI-NEB的收敛不了怎么回事？

查看数: 1714 | 评论数: 9 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

发布时间: 2024-10-22 20:06

用cp2k的计算两层2H-MoS2滑移时能垒，但使用的CI-NEB收敛不了怎么回事？我使用的sobNEB插了七个点，然后用multiwfn输出NEB.inp，其中设置了周期性为XYZ，理论方法PBE，开启了smearing，k点设置为3，3，1，还加入了 ...

回复 Reply

sobereva 发表于 Post on 2024-12-25 18:21:31

jifeng131381 发表于 2024-12-25 17:46
老师问一下楼主计算两层2H-MoS2滑移时能垒这种体系，结构优化和CI-NEB使用XY周期性可以吗,和XYZ周期性 ...

比较薄的slab体系用二维周期性结合MT，比三维周期性结合更划算

细节可以参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里“能量的计算及相关问题”的“低维体系的计算”部分。

jifeng131381 发表于 Post on 2024-12-25 17:46:39

sobereva 发表于 2024-10-26 23:49
显然说明中间有极小点，拿第7个点优化极小点
之前我的回复里说了参考另一个帖子里我的回复，搞清楚原理
...

老师问一下楼主计算两层2H-MoS2滑移时能垒这种体系，结构优化和CI-NEB使用XY周期性可以吗,和XYZ周期性相比是否有差别 (我是金属表面脱氯)

xxxx531 发表于 Post on 2024-10-27 19:18:56

sobereva 发表于 2024-10-26 23:49
显然说明中间有极小点，拿第7个点优化极小点
之前我的回复里说了参考另一个帖子里我的回复，搞清楚原理
...

好的，谢谢老师

sobereva 发表于 Post on 2024-10-26 23:49:13

xxxx531 发表于 2024-10-26 23:29
老师，我插了12个点，还是出现两个峰的状态，而且始末态和过渡态的单点能和我neb算出来的能垒的值不一样 ...

显然说明中间有极小点，拿第7个点优化极小点
之前我的回复里说了参考另一个帖子里我的回复，搞清楚原理

sobereva 发表于 Post on 2024-10-24 12:27:48

xxxx531 发表于 2024-10-24 11:07
老师，我现在跑了164步之后收敛了，但是画出来的能垒图有两个最高点，而且和我之前算的起末态和中间过渡 ...

点数太少时，如果期间还有中间体，显然体现不出来

参考http://bbs.keinsci.com/thread-49211-1-1.html里我的回复

xxxx531 发表于 Post on 2024-10-23 16:15:40

sobereva 发表于 2024-10-23 12:39
可以尝试多跑100步
很多情况100步NEB都不足以达到收敛。从BAND TOTAL ENERGY [au]可以看到能量还在单调变 ...

好的老师，我再尝试一下。

sobereva 发表于 Post on 2024-10-23 12:39:53

可以尝试多跑100步
很多情况100步NEB都不足以达到收敛。从BAND TOTAL ENERGY [au]可以看到能量还在单调变化而没有震荡

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register