计算化学公社

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

[CP2K] NEB跑大概结构计算虚频,使用虚频跑DIMMER,得到的结构受力收敛但无虚频是什么原因?

[复制链接 Copy URL]

20

帖子

0

威望

145

eV
积分
165

Level 3 能力者

本帖最后由 大耳朵涂涂 于 2024-11-27 11:16 编辑

各位大佬好,最近在使用CP2K计算固体催化剂催化脱氢反应。计算思路如下:先使用CI-NEB插点跑出大概的构型,然后选取合适的结构计算频率,用虚频跑DIMMER最后拿到过渡态结构。
最近在计算一个表面氧脱氢反应,反应初态和末态如下图,是氧原子脱除亚甲基上的氢原子。


遇到的问题是这样的:使用CI-NEB跑出的大概构型,然后根据构型的虚频跑DIMMER,最后得到的结构却没有虚频。然后我分析了NEB得到的构型及振动方向如下图。最大的虚频是氢原子朝向碳原子方向振动,DIMMER计算的过程也是氢原子慢慢向碳原子靠近,最后稳定在一个距离上(C-H键长1.17A),但是根据最后稳定的这个结构计算频率,却没有氢原子在碳原子和氧原子之间震动的虚频了。


有两个问题想请教一下大家:
1.为什么DIMMER跑到最后没有虚频了,是设置不对吗?
2.这个反应是一个脱氢反应,反应路径是氢原子从碳原子跑到氧原子那里,从NEB的结果看,能量先升高后降低。但是中间构型的虚频方向却是氢原子指向氧原子方向。频率的确切方向,比如向左或是向右,与反应的路径和原子移动的方向是否有关联?PS:我计算了氢原子靠近氧原子的结构的频率,虚频指向是氢原子指向氧原子,也就是所氢原子靠近碳原子的结构虚频方向指向碳原子,靠近氧原子的结构虚频方向指向氧原子,这么看虚频方向与反应方向相反的。

附件为中间态和DIMMER计算结果结构的频率计算设置及结果:
[url=]NEB中间态频率计算.inp[/url]
[url=][url=]NEB中间态频率计算结果.mol[/url][/url]
[url=][url=][url=]DIMMER结果结构频率计算.inp[/url][/url][/url]
[url=][url=][url=][url=]DIMMER结果结构频率计算结果.mol[/url][/url][/url][/url]


DIMMER结果结构频率计算结果.mol

289.51 KB, 下载次数 Times of downloads: 1

DIMMER结果结构频率计算.inp

25.81 KB, 下载次数 Times of downloads: 1

NEB中间态频率计算.inp

25.79 KB, 下载次数 Times of downloads: 1

NEB中间态频率计算结果.mol

289.51 KB, 下载次数 Times of downloads: 1

本版积分规则 Credits rule

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

GMT+8, 2024-11-27 21:37 , Processed in 0.617167 second(s), 24 queries , Gzip On.

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