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

xxxx531 · 发表于 Post on 2024-10-22 20:06:26

用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这一项，其他三项基本还差得多。请问是步数设置不够大吗?还是其他参数没有设置好的原因？

插点文件：

插点文件.zip (631.71 KB, 下载次数 Times of downloads: 4)
计算文件：

neb.zip (26.22 KB, 下载次数 Times of downloads: 7)

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

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

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

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

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

xxxx531 · 发表于 Post on 2024-10-24 11:07:38

本帖最后由 xxxx531 于 2024-10-24 11:11 编辑

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

老师，我现在跑了164步之后收敛了，但是画出来的能垒图有两个最高点，而且和我之前算的起末态和中间过渡态的单点能不一致，这是怎么回事？起末态和中间过渡态对于的插点是第一个、第七个、第四个点。

单点能和neb文件如下：

单点能.zip (30.15 KB, 下载次数 Times of downloads: 1)

neb3.zip (34.1 KB, 下载次数 Times of downloads: 2)

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-26 23:29:20

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

单点能.zip (30.15 KB, 下载次数 Times of downloads: 0)

2.xyz (1.3 KB, 下载次数 Times of downloads: 0)

3.xyz (1.3 KB, 下载次数 Times of downloads: 0)

4.xyz (1.3 KB, 下载次数 Times of downloads: 0)

5.xyz (1.3 KB, 下载次数 Times of downloads: 0)

6.xyz (1.3 KB, 下载次数 Times of downloads: 0)

8.xyz (1.3 KB, 下载次数 Times of downloads: 0)

9.xyz (1.3 KB, 下载次数 Times of downloads: 0)

10.xyz (1.3 KB, 下载次数 Times of downloads: 0)

11.xyz (1.3 KB, 下载次数 Times of downloads: 0)

conf1.xyz (1.34 KB, 下载次数 Times of downloads: 0)

conf2.xyz (1.34 KB, 下载次数 Times of downloads: 0)

TS_guess.xyz (1.34 KB, 下载次数 Times of downloads: 0)

traj.xyz (15.55 KB, 下载次数 Times of downloads: 0)

neb.inp (6.8 KB, 下载次数 Times of downloads: 0)

neb.out (336.3 KB, 下载次数 Times of downloads: 0)

neb-1.ener (149.07 KB, 下载次数 Times of downloads: 0)

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

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

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

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

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

好的，谢谢老师

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

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