计算化学公社

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

[CP2K] 搜索过渡态时通过什么方法可以预测反应物的第一步反应是什么?

[复制链接 Copy URL]

128

帖子

0

威望

613

eV
积分
741

Level 4 (黑子)

各位好,我的反应物和催化剂的吸附构型如下,由于我的反应路径不明确,我根据一些文章做的实验产物猜测中间产物进行第一步过渡态的搜索计算,可试了一些情况,始终得不到过渡态先增后减的趋势,请问有什么好的办法可以预测第一步反应是发生在哪吗?本人非化学专业出身,基础比较弱,望各位老师不吝赐教,非常感谢!

128

帖子

0

威望

613

eV
积分
741

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2024-2-28 19:10:41 | 只看该作者 Only view this author
Exodus 发表于 2024-2-27 16:29
那把氢原子移动到中间位置做一次单点能看一下,这么短的距离做几次单点能基本就能确认有没有过渡态了

老师,我将上面的反应物和产物结构作为始末态进行过渡态搜索,收敛后的结果显示从反应物到产物的能量是一直递减的,这是不是说明该步反应没有过渡态?接下来我该怎么调整尝试,是不是要尝试换其他的中间产物进行过渡态的搜索呢!

36

帖子

0

威望

2430

eV
积分
2466

Level 5 (御坂)

10#
发表于 Post on 2024-2-27 16:29:08 | 只看该作者 Only view this author
那把氢原子移动到中间位置做一次单点能看一下,这么短的距离做几次单点能基本就能确认有没有过渡态了

128

帖子

0

威望

613

eV
积分
741

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2024-2-27 15:45:08 | 只看该作者 Only view this author
Exodus 发表于 2024-2-27 13:45
从这个结构来看似乎氢原子只移动了不到1埃的样子,量一下中间体O-H的距离,确定O-H键断开了,顺便先看一下这 ...

老师,氢原子移动距离大概才0.4埃左右,优化后的反应物O-H距离为1.0725,H-催化剂的距离为1.4492,能量为-11643.70577;优化后的产物O-H距离为1.4678,H-催化剂距离为1.0554,能量是-11643.7347。计算过中间产物的连接性指数和键级分别为0.29和0.43,说明中间产物的O-H键应该是断开的。请问老师我还可以进行哪些调整尝试呢?

36

帖子

0

威望

2430

eV
积分
2466

Level 5 (御坂)

8#
发表于 Post on 2024-2-27 13:45:47 | 只看该作者 Only view this author
从这个结构来看似乎氢原子只移动了不到1埃的样子,量一下中间体O-H的距离,确定O-H键断开了,顺便先看一下这步的热力学变化。单从这个结构来看变化很小,很可能这个所谓的中间体只是因为吸附导致的O-H键发生略微的变化,实际并没有解离,可以把中间体的分子往上移动一点再优化,看它还会不会掉下来

128

帖子

0

威望

613

eV
积分
741

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2024-2-27 12:12:24 | 只看该作者 Only view this author
Exodus 发表于 2024-2-25 20:44
如果是热催化反应,这样一个完美的表面看起来不像有活性的样子,光催化则不一定;至少应该先优化得到基元反 ...

老师,我的是热催化,初末结构优化完成了的,如上回复中的图片,可以帮我看看嘛,非常感谢

128

帖子

0

威望

613

eV
积分
741

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2024-2-26 10:34:53 | 只看该作者 Only view this author
@Exodus 老师,我的是热催化反应,我依据猜测的机理在反应物优化结构(左图)基础上把酚羟基的O-H键断开,并把H原子靠近催化剂层进行了产物的结构优化(右图),在得到反应物和产物的优化结构后才进行的CI-NEB过渡态搜索,可始终得不到理想的结果。

202402261031563090..png (166.63 KB, 下载次数 Times of downloads: 4)

202402261031563090..png

202402261034089913..png (143.72 KB, 下载次数 Times of downloads: 4)

202402261034089913..png

36

帖子

0

威望

2430

eV
积分
2466

Level 5 (御坂)

5#
发表于 Post on 2024-2-25 20:44:05 | 只看该作者 Only view this author
如果是热催化反应,这样一个完美的表面看起来不像有活性的样子,光催化则不一定;至少应该先优化得到基元反应初末态的结构,再构建过渡态的结构。

128

帖子

0

威望

613

eV
积分
741

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2024-2-25 19:32:06 | 只看该作者 Only view this author
第一步产物是根据这个去除了酚羟基上的氢

202402251931016088..png (19.33 KB, 下载次数 Times of downloads: 6)

202402251931016088..png

128

帖子

0

威望

613

eV
积分
741

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2024-2-25 19:30:43 | 只看该作者 Only view this author
Exodus 发表于 2024-2-24 20:36
实际的催化剂真的是一个perfect的表面吗,你的产物是什么呢,至少应该依据产物来推测一个合理的反应机理再 ...

谢谢老师的解答,请问一下如何确定催化剂是否是perfect表面呢?产物的话我参考一些实验文章猜测的结论,第一步反应是脱除笨环酚羟基上的氢,仅仅是第一步脱除一个H的情况我在真空和引入隐式溶剂水下计算了很久都没有得出想要的结果来。

36

帖子

0

威望

2430

eV
积分
2466

Level 5 (御坂)

2#
发表于 Post on 2024-2-24 20:36:28 | 只看该作者 Only view this author
实际的催化剂真的是一个perfect的表面吗,你的产物是什么呢,至少应该依据产物来推测一个合理的反应机理再开始计算

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

GMT+8, 2024-11-25 05:25 , Processed in 0.228875 second(s), 24 queries , Gzip On.

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