计算化学公社

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

[VASP] VASP算dimer得到的过渡态结构有俩虚频

[复制链接 Copy URL]

34

帖子

0

威望

1000

eV
积分
1034

Level 4 (黑子)

我首先用CINEB进行了粗收敛计算(EDIFF=1E-6, EDIFFG=-0.2)。反应路径以及能垒都比较正常。然后用了neb2dim.pl 生成了dimer的计算文件,把CINEB得到的能量最高的image的CONTCAR复制进去当POSCAR计算。MODECAR就用了脚本自动生成的。然后按标准设置了dimer计算的参数进行计算。很多反应的确得到了合理解决。就是现在有个N-N键解离的过程,我dimer算出来做频率校正发现了俩虚频。其中一个虚频是同断键的方向的(下图虚频379 cm-1)。那我接下来要怎么操作可以继续往一阶鞍点收敛呢?用原先dimer得到的CONTCAR以及MODECAR继续进行dimer计算可以吗?还是说需要根据频率校正得到的俩虚频重新制作一个MODECAR,然后用CONTCAR结构跑dimer呢?

下为虚频情况


   


越学越难顶。

3

帖子

0

威望

13

eV
积分
16

Level 1 能力者

16#
发表于 Post on 2024-9-20 09:12:36 | 只看该作者 Only view this author
小雨呢喃 发表于 2022-6-14 13:30
第二个链接里边的image.py没有呢能发一下吗 谢谢

就是exfreq.py

41

帖子

0

威望

862

eV
积分
903

Level 4 (黑子)

15#
发表于 Post on 2023-12-11 09:10:34 | 只看该作者 Only view this author
Jiongci-21 发表于 2023-12-10 12:59
这个怎么说呢,我觉得这是一种不错的尝试吧。如果你先用0.5粗收敛一个路径,化学上看着也合理,可以把收 ...

收到,非常感谢

22

帖子

0

威望

828

eV
积分
850

Level 4 (黑子)

14#
发表于 Post on 2023-12-11 09:08:55 | 只看该作者 Only view this author
频率做MODECAR,重新计算就行了

34

帖子

0

威望

1000

eV
积分
1034

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2023-12-10 12:59:57 | 只看该作者 Only view this author
本帖最后由 Jiongci-21 于 2023-12-10 13:03 编辑
CHLllooll 发表于 2023-11-25 20:22
您好,我想请教一下在作过渡态搜索的时候粗收敛用-0.5精度够吗,我看刘老师之前发的说是粗收敛可以用-0.5 ...

这个怎么说呢,我觉得这是一种不错的尝试吧。如果你先用0.5粗收敛一个路径,化学上看着也合理,可以把收敛参数调高一点看看能不能收敛得到准确的结构和能量。成功的话事半功倍。
直接用精确的电子收敛和力收敛参数去跑初猜的过渡态路径和结构肯定是不太可能一次性就拿到结果的,程序很容易就把结构跑崩了。
我之前比较喜欢用0.1-0.2的力收敛参数用NI-NEB去优化我的初猜,然后拿到合理的过渡态结构直接跑Dimer,计算成本小,成功率还高,主要是优化几轮后你可以检查下任务,如果结构出问题了可以及时kill了再调整。
越学越难顶。

41

帖子

0

威望

862

eV
积分
903

Level 4 (黑子)

12#
发表于 Post on 2023-11-25 20:22:52 | 只看该作者 Only view this author
Jiongci-21 发表于 2022-6-14 21:27
是的,我也没找到,所以只是学习一下方法,没用他的脚本。我直接用的官方neb2dim.pl 脚本生成的dimer计算 ...

您好,我想请教一下在作过渡态搜索的时候粗收敛用-0.5精度够吗,我看刘老师之前发的说是粗收敛可以用-0.5,但我在-0.5精度下收敛时的eps看着很正常并且结构也正常;同样的过渡态搜索如果缩小精度到-0.05的时候一直不收敛并且eps的弹簧力图很不平滑还有谷那种的,到那时一直是不收敛的,并且还有可能中间的几个结构跑乱导致中间的能量会有一个下降的值。

12

帖子

0

威望

227

eV
积分
239

Level 3 能力者

11#
发表于 Post on 2022-6-15 16:00:37 | 只看该作者 Only view this author
Jiongci-21 发表于 2022-6-14 21:27
是的,我也没找到,所以只是学习一下方法,没用他的脚本。我直接用的官方neb2dim.pl 脚本生成的dimer计算 ...

我也是这么做的

34

帖子

0

威望

1000

eV
积分
1034

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2022-6-14 21:27:13 | 只看该作者 Only view this author
小雨呢喃 发表于 2022-6-14 13:30
第二个链接里边的image.py没有呢能发一下吗 谢谢

是的,我也没找到,所以只是学习一下方法,没用他的脚本。我直接用的官方neb2dim.pl 脚本生成的dimer计算文件
越学越难顶。

12

帖子

0

威望

227

eV
积分
239

Level 3 能力者

9#
发表于 Post on 2022-6-14 13:30:42 | 只看该作者 Only view this author
ahxb 发表于 2022-6-6 09:10
去看VTST官网上对Dimer算法的介绍。

谢谢

12

帖子

0

威望

227

eV
积分
239

Level 3 能力者

8#
发表于 Post on 2022-6-14 13:30:06 | 只看该作者 Only view this author
Jiongci-21 发表于 2022-6-6 10:48
去官网看,然后公社里也有个挺好的帖子
VTST官网:
http://theory.cm.utexas.edu/vtsttools/dimer.html ...

第二个链接里边的image.py没有呢能发一下吗 谢谢

34

帖子

0

威望

1000

eV
积分
1034

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2022-6-6 10:48:11 | 只看该作者 Only view this author
小雨呢喃 发表于 2022-6-5 15:59
请问NEB后如何跑dimer---有教程吗

去官网看,然后公社里也有个挺好的帖子
VTST官网:
http://theory.cm.utexas.edu/vtsttools/dimer.html

计算化学公社帖子:
vasp做过渡态经验<二>:Dimer
http://bbs.keinsci.com/forum.php ... 3&fromuid=19287

越学越难顶。

125

帖子

0

威望

2168

eV
积分
2293

Level 5 (御坂)

6#
发表于 Post on 2022-6-6 09:10:08 | 只看该作者 Only view this author
小雨呢喃 发表于 2022-6-5 15:59
请问NEB后如何跑dimer---有教程吗

去看VTST官网上对Dimer算法的介绍。

12

帖子

0

威望

227

eV
积分
239

Level 3 能力者

5#
发表于 Post on 2022-6-5 15:59:13 | 只看该作者 Only view this author
请问NEB后如何跑dimer---有教程吗

34

帖子

0

威望

1000

eV
积分
1034

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2022-4-19 21:05:53 | 只看该作者 Only view this author
本帖最后由 Jiongci-21 于 2022-4-19 21:07 编辑
易水曦 发表于 2022-4-19 09:37
请问楼主解决了嘛 想知道你怎么处理的

没有找到很好的方案,就是重新调了下结构和收敛标准跑NEB,再同样跑dimer找到了正确的过渡态。
至于怎么在上述dimer跑出来的具有俩虚频的结构实现去虚频找到过渡态希望可以得到高人指点。

越学越难顶。

1

帖子

0

威望

277

eV
积分
278

Level 3 能力者

3#
发表于 Post on 2022-4-19 09:37:57 | 只看该作者 Only view this author
请问楼主解决了嘛 想知道你怎么处理的

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

GMT+8, 2024-11-24 01:52 , Processed in 0.223219 second(s), 31 queries , Gzip On.

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