计算化学公社

标题: NEB跑大概结构计算虚频，使用虚频跑DIMMER，得到的结构受力收敛但无虚频是什么原因？ [打印本页]

作者
Author: 大耳朵涂涂 时间: 2024-11-27 11:16
标题: NEB跑大概结构计算虚频，使用虚频跑DIMMER，得到的结构受力收敛但无虚频是什么原因？
本帖最后由大耳朵涂涂于 2024-11-27 11:16 编辑

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

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

有两个问题想请教一下大家：
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]

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