计算化学公社

标题: 求助：用cp2k计算MoS2时使用CI-NEB的收敛不了怎么回事？ [打印本页]

作者
Author: xxxx531 时间: 2024-10-22 20:06
标题: 求助：用cp2k计算MoS2时使用CI-NEB的收敛不了怎么回事？
用cp2k的计算两层2H-MoS2滑移时能垒，但使用的CI-NEB收敛不了怎么回事？
我使用的sobNEB插了七个点，然后用multiwfn输出NEB.inp，其中设置了周期性为XYZ，理论方法PBE，开启了smearing，k点设置为3，3，1，还加入了范德华修正：
&VDW_POTENTIAL
      POTENTIAL_TYPE  PAIR_POTENTIAL
      &PAIR_POTENTIAL
      R_CUTOFF 9    #Range of potential. The cutoff will be 2 times of this value
      TYPE       DFTD3
      PARAMETER_FILE_NAME  dftd3.dat       REFERENCE_FUNCTIONAL PBE
      LONG_RANGE_CORRECTION T
      &END PAIR_POTENTIAL
跑了15个小时100步自动停止了，但是输出结果只有一个能够达到官方收敛标准的始终只有RMS FORCE这一项，其他三项基本还差得多。请问是步数设置不够大吗?还是其他参数没有设置好的原因？
(, 下载次数 Times of downloads: 0)
插点文件： (, 下载次数 Times of downloads: 4)
计算文件： (, 下载次数 Times of downloads: 7)

作者
Author: sobereva 时间: 2024-10-23 12:39
可以尝试多跑100步
很多情况100步NEB都不足以达到收敛。从BAND TOTAL ENERGY [au]可以看到能量还在单调变化而没有震荡

作者
Author: xxxx531 时间: 2024-10-23 16:15

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

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

作者
Author: xxxx531 时间: 2024-10-24 11:07
本帖最后由 xxxx531 于 2024-10-24 11:11 编辑

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

老师，我现在跑了164步之后收敛了，但是画出来的能垒图有两个最高点，而且和我之前算的起末态和中间过渡态的单点能不一致，这是怎么回事？起末态和中间过渡态对于的插点是第一个、第七个、第四个点。
(, 下载次数 Times of downloads: 0)
(, 下载次数 Times of downloads: 0)
(, 下载次数 Times of downloads: 0)
单点能和neb文件如下： (, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 2)

作者
Author: sobereva 时间: 2024-10-24 12:27

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

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

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

作者
Author: xxxx531 时间: 2024-10-26 23:29
本帖最后由 xxxx531 于 2024-10-26 23:38 编辑

sobereva 发表于 2024-10-24 12:27
点数太少时，如果期间还有中间体，显然体现不出来

参考http://bbs.keinsci.com/thread-49211-1-1.html ...

老师，我插了12个点，还是出现两个峰的状态，而且始末态和过渡态的单点能和我neb算出来的能垒的值不一样呢？起末态和中间过渡态对于的插点是第1个、第12个、第7个点。
(, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0)
(, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0)

作者
Author: sobereva 时间: 2024-10-26 23:49

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

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

作者
Author: xxxx531 时间: 2024-10-27 19:18

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

好的，谢谢老师

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