计算化学公社

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

[Gaussian/gview] 求助-如何计算pt催化的油酸加成反应的过渡态及能垒

[复制链接 Copy URL]

10

帖子

0

威望

21

eV
积分
31

Level 2 能力者

我想通过gaussian 09w和gaussview计算铂催化的油酸(含有一个双键的18碳脂肪酸)的氢化反应。目前已经完成结构的优化和频率分析,想对过渡态和反应能垒进行分析,但是不知道接下来该如何做,求大神指点。1、我该如何画gjf文件,氢 铂和油酸是一个什么样的状态?怎么才能体现催化剂的作用?反应机理如图。

2、我自己尝试着用QST2做了一次过渡态分析,结果出来有唯一虚频,但是out文件的分子结构不太正确。

3.、我在哪里可以设置反应的温度?


不是学化学的,有很多地方不明白,请大神指点,感谢。


out文件图片.png (59.01 KB, 下载次数 Times of downloads: 43)

out文件,末端H原子没有键,pt也链接着H

out文件,末端H原子没有键,pt也链接着H

氢化.png (64.94 KB, 下载次数 Times of downloads: 34)

氢化.png

qinghua2.gjf

8.01 KB, 下载次数 Times of downloads: 2

qinghua2.chk

5.35 MB, 下载次数 Times of downloads: 0

1633

帖子

4

威望

4093

eV
积分
5806

Level 6 (一方通行)

喵星人

2#
发表于 Post on 2020-8-12 21:01:18 | 只看该作者 Only view this author
你觉得催化剂是单原子Pt?

5万

帖子

99

威望

5万

eV
积分
112440

管理员

公社社长

3#
发表于 Post on 2020-8-13 08:21:44 | 只看该作者 Only view this author
1 至少得用一个Pt团簇才能描述Pt金属
什么叫团簇参考
使用量子化学程序基于簇模型计算金属表面吸附问题
http://sobereva.com/540http://bbs.keinsci.com/thread-16532-1-1.html

2 毫无有效信息,别人也不知道具体是什么结构、虚频什么模式。没事甭用QST2,这里说了
简谈Gaussian里找过渡态的关键词opt=TS和QST2、3
http://sobereva.com/460http://bbs.keinsci.com/thread-12093-1-1.html

3 温度根本不影响找到的过渡态,仔细看
谈谈温度、压力、同位素设定对量子化学计算结果产生的影响
http://sobereva.com/423http://bbs.keinsci.com/thread-10159-1-1.html

算热力学量,比如算与反应速率关系密切的自由能垒,才体现温度,看
使用Shermo结合量子化学程序方便地计算分子的各种热力学数据
http://sobereva.com/552http://bbs.keinsci.com/thread-17494-1-1.html
谈谈如何通过势垒判断反应是否容易发生
http://sobereva.com/506http://bbs.keinsci.com/thread-14313-1-1.html

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

10

帖子

0

威望

21

eV
积分
31

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2020-8-13 14:38:31 | 只看该作者 Only view this author
喵星大佬 发表于 2020-8-12 21:01
你觉得催化剂是单原子Pt?

见笑了,外行人实在是不明白,需要怎么改进还请指点,感谢

10

帖子

0

威望

21

eV
积分
31

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2020-8-13 14:39:48 | 只看该作者 Only view this author
sobereva 发表于 2020-8-13 08:21
1 至少得用一个Pt团簇才能描述Pt金属
什么叫团簇参考
使用量子化学程序基于簇模型计算金属表面吸附问题

多谢sob老师!我以为直接用原子pt就能代表,关于团簇的知识我会认真学习的

10

帖子

0

威望

21

eV
积分
31

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2020-8-13 20:18:16 | 只看该作者 Only view this author
sobereva 发表于 2020-8-13 08:21
1 至少得用一个Pt团簇才能描述Pt金属
什么叫团簇参考
使用量子化学程序基于簇模型计算金属表面吸附问题

sob老师,我用gaussview该怎么做团簇呢,有没有推荐的资料,我看你发的那个是其他软件的。谢谢老师

1633

帖子

4

威望

4093

eV
积分
5806

Level 6 (一方通行)

喵星人

7#
发表于 Post on 2020-8-13 22:33:01 | 只看该作者 Only view this author
Fiee 发表于 2020-8-13 14:38
见笑了,外行人实在是不明白,需要怎么改进还请指点,感谢

所以你要算的反应,具体是用什么催化,是Pt黑?Pt/C?PtO2?还是均相的配合物?啥都不说别人怎么可能知道你要算的是啥

831

帖子

1

威望

7194

eV
积分
8045

Level 6 (一方通行)

8#
发表于 Post on 2020-8-13 22:42:45 | 只看该作者 Only view this author
Fiee 发表于 2020-8-13 20:18
sob老师,我用gaussview该怎么做团簇呢,有没有推荐的资料,我看你发的那个是其他软件的。谢谢老师

文中用到orca的部分用gaussian代替没有问题,方法基组都给了。

10

帖子

0

威望

21

eV
积分
31

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2020-8-14 11:58:49 | 只看该作者 Only view this author
喵星大佬 发表于 2020-8-13 22:33
所以你要算的反应,具体是用什么催化,是Pt黑?Pt/C?PtO2?还是均相的配合物?啥都不说别人怎么可能知道 ...

是Pt/C,我以为只是计算铂呢,C只是一个载体,也要画出来吗?我不是学化学的,很懵,麻烦你了

10

帖子

0

威望

21

eV
积分
31

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2020-8-14 11:59:17 | 只看该作者 Only view this author
hebrewsnabla 发表于 2020-8-13 22:42
文中用到orca的部分用gaussian代替没有问题,方法基组都给了。

哦哦,我以为是不一样的呢,谢谢你

1633

帖子

4

威望

4093

eV
积分
5806

Level 6 (一方通行)

喵星人

11#
发表于 Post on 2020-8-14 23:08:18 | 只看该作者 Only view this author
Fiee 发表于 2020-8-14 11:58
是Pt/C,我以为只是计算铂呢,C只是一个载体,也要画出来吗?我不是学化学的,很懵,麻烦你了

主要是暴露的表面还有晶形可能不一样,要看怎么处理。
你这个在表面上反应的分子不小,固体的原子也不好算,不管是扣个团簇出来还是用平面波程序处理都很不好整

10

帖子

0

威望

21

eV
积分
31

Level 2 能力者

12#
 楼主 Author| 发表于 Post on 2020-8-17 08:53:54 | 只看该作者 Only view this author
喵星大佬 发表于 2020-8-14 23:08
主要是暴露的表面还有晶形可能不一样,要看怎么处理。
你这个在表面上反应的分子不小,固体的原子也不好 ...

哦哦,我刚开始以为就用一个pt原子就可以表示催化剂呢,我这个是111晶型,我想能不能简化呢,我对计算的精度要求不是很高,只要就是看一下反应能垒和过渡态

10

帖子

0

威望

21

eV
积分
31

Level 2 能力者

13#
 楼主 Author| 发表于 Post on 2020-8-17 09:06:05 | 只看该作者 Only view this author
hebrewsnabla 发表于 2020-8-13 22:42
文中用到orca的部分用gaussian代替没有问题,方法基组都给了。

你好,我想请教一下高斯怎么设置冻结原子呢

831

帖子

1

威望

7194

eV
积分
8045

Level 6 (一方通行)

14#
发表于 Post on 2020-8-17 10:18:40 | 只看该作者 Only view this author
Fiee 发表于 2020-8-17 09:06
你好,我想请教一下高斯怎么设置冻结原子呢

http://sobereva.com/404

1万

帖子

0

威望

7400

eV
积分
18164

Level 6 (一方通行)

15#
发表于 Post on 2020-8-17 11:16:31 | 只看该作者 Only view this author
喵星大佬 发表于 2020-8-14 23:08
主要是暴露的表面还有晶形可能不一样,要看怎么处理。
你这个在表面上反应的分子不小,固体的原子也不好 ...

我觉得油酸截断到顺式-3-己烯问题不大,然后就可以用平面波程序算了。需要注意的是表面氢原子的覆盖度问题,Pt表面光秃秃的只有一两个吸附氢原子可能不太现实,有可能已经吸附了一层或半层氢原子。
另外这么重要且基础的反应,没准已经有文献算过了吧
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

本版积分规则 Credits rule

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

GMT+8, 2024-11-26 07:42 , Processed in 0.177791 second(s), 23 queries , Gzip On.

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