计算化学公社

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

[新手求助] gaussian能否模拟超临界状态下反应机理

[复制链接 Copy URL]

29

帖子

0

威望

133

eV
积分
162

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 xuexi123. 于 2021-12-24 16:46 编辑

大家好,我专业是材料化工的,最近课题牵扯到一个超临界状态下的一个反应,老师的意思是通过分子模拟判断是否可以实现如下目标,之前从来没有接触过这些方面的内容和软件应用,课题组也从来没有老师做过这一块,没人有可以给我指导,所以有点迷茫,希望大家可以给出意见(第一次发帖,不知道这样合格不):超临界状态下进行羟醛缩合反应,想要优化反应物,产物的结构,寻找过渡态,找到一条合适的反应路径,并计算活化能,反应速率常数。因为我这是超临界状态下,物料是处于气液不分的状态,物性,热力学参数肯定会发生变化,所以通过高斯模拟能得到以下数据吗?
1.反应物,过渡态,产物优化结构(这点gaussian可以完成);
2.计算得到超临界状态下的反应路径,计算反应能垒,反应速率常数。
3.计算过渡态的时候,甲醛质子化这一步存在过渡态吗?
最主要是问题是:因为此反应是超临界状态下,在此状态下,物质的物性,热力学参数应该会发生变化,我不知道计算的时候,计算什么东西的时候,该怎么去考虑超临界这个因素。
谢谢大家



82

帖子

0

威望

1197

eV
积分
1279

Level 4 (黑子)

2#
发表于 Post on 2021-12-24 14:57:25 | 只看该作者 Only view this author
这么模糊不清的标题肯定是要被骂的

29

帖子

0

威望

133

eV
积分
162

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2021-12-24 16:43:22 | 只看该作者 Only view this author
Rosefinch 发表于 2021-12-24 14:57
这么模糊不清的标题肯定是要被骂的

我改了  我确实不好描述我的问题,你可以帮忙看一下我的问题么,就是超临界状态下进行羟醛缩合反应,想要优化反应物,产物的结构,寻找过渡态,找到一条合适的反应路径,并计算活化能,反应速率常数。因为我这是超临界状态下,物料是处于气液不分的状态...

1

帖子

0

威望

51

eV
积分
52

Level 2 能力者

4#
发表于 Post on 2021-12-24 17:00:04 | 只看该作者 Only view this author
设置一下温度,压力,溶液也许可以?

29

帖子

0

威望

133

eV
积分
162

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2021-12-25 16:36:02 | 只看该作者 Only view this author
山有木兮 发表于 2021-12-24 17:00
设置一下温度,压力,溶液也许可以?

温度 压力对于结构优化,频率分析等没有多大影响,请问gaussian可以做动力学研究么?我需要先判断能不能,再去考虑如何做

6万

帖子

99

威望

6万

eV
积分
125136

管理员

公社社长

6#
发表于 Post on 2021-12-25 19:25:16 | 只看该作者 Only view this author
Gaussian的计算没法考虑超临界状态
Gaussian做几何优化的对象是微观结构,超临界状态根本不从这个计算层面来体现。而且本身势能面也不依赖于外环境。
注意看此文
谈谈温度、压力、同位素设定对量子化学计算结果产生的影响
http://sobereva.com/423http://bbs.keinsci.com/thread-10159-1-1.html

至于热力学性质,那是受压力和温度影响的,但Gaussian算热力学量(也包括Shermo等计算热力学量的程序)都是把体系当理想气体看待(顶多再同时考虑隐式溶剂模型),没法表现超临界状态。只有有没有一些额外的经验方式去考虑(诸如用特殊参数的隐式溶剂模型之类),搜相关研究文献。

质子化有没有过渡态问题,看谁把甲醛质子化。

还要搞清楚,活化能不是算出来的,而是拟合出来的,仔细看
谈谈如何通过势垒判断反应是否容易发生
http://sobereva.com/506http://bbs.keinsci.com/thread-14313-1-1.html

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

29

帖子

0

威望

133

eV
积分
162

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2021-12-29 21:04:52 | 只看该作者 Only view this author
sobereva 发表于 2021-12-25 19:25
Gaussian的计算没法考虑超临界状态
Gaussian做几何优化的对象是微观结构,超临界状态根本不从这个计算层面 ...

谢谢老师,我还有问题想请教你,我主要是模拟超临界态下羟醛缩合反应机理,比如丙酮和甲醛反应,此状态下无需添加催化剂,甲醛会发生自身歧化反应生成一份子甲醇和甲酸,然后由于甲酸的存在,所以整个溶液显酸性,也就是说甲醛质子化是因为甲酸的存在,如果我要研究丙酮和甲醛的反应(甲醛质子化是其中一步反应)机理,寻找过渡态,验证反应物和产物之间的关系。是不是需要从甲醛自身歧化生成甲酸这一步开始计算啊。

6万

帖子

99

威望

6万

eV
积分
125136

管理员

公社社长

8#
发表于 Post on 2021-12-30 00:13:54 | 只看该作者 Only view this author
xuexi123. 发表于 2021-12-29 21:04
谢谢老师,我还有问题想请教你,我主要是模拟超临界态下羟醛缩合反应机理,比如丙酮和甲醛反应,此状态下 ...

这个不是重点,可以忽略。而且如果能搜到有文章算过歧化的反应机理,那就更没必要自己算了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

29

帖子

0

威望

133

eV
积分
162

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2022-3-25 23:12:16 | 只看该作者 Only view this author
sobereva 发表于 2021-12-30 00:13
这个不是重点,可以忽略。而且如果能搜到有文章算过歧化的反应机理,那就更没必要自己算了

老师你好,首先甲醛质子化是存在酸,如果我不去模拟酸是如何产生的,直接模拟后面甲醛质子化这样也是可以的嘛?

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

10#
发表于 Post on 2022-3-26 00:16:04 | 只看该作者 Only view this author
xuexi123. 发表于 2022-3-25 16:12
老师你好,首先甲醛质子化是存在酸,如果我不去模拟酸是如何产生的,直接模拟后面甲醛质子化这样也是可以 ...

不建议这样做,因为这样的话你不知道质子化甲醛和非质子化甲醛的浓度比,搞不好虽然算出来反应似乎可以发生,但是因为质子化甲醛的平衡浓度太低,实际几乎不反应。
可以不算质子化的过渡态,但是质子化的热力学是必须要算的,或者用可靠的实验结果代替。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

29

帖子

0

威望

133

eV
积分
162

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2022-3-26 00:27:28 | 只看该作者 Only view this author
wzkchem5 发表于 2022-3-26 00:16
不建议这样做,因为这样的话你不知道质子化甲醛和非质子化甲醛的浓度比,搞不好虽然算出来反应似乎可以发 ...

多谢回复,这个实验我们已经做过了,是可以发生了,所以甲醛质子化就是反应机理中其中一步。

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

12#
发表于 Post on 2022-3-26 01:57:40 | 只看该作者 Only view this author
xuexi123. 发表于 2022-3-25 17:27
多谢回复,这个实验我们已经做过了,是可以发生了,所以甲醛质子化就是反应机理中其中一步。

那可以不算,但是也要留意一下质子化甲醛的浓度是否会影响后续反应的分支比,比如有的物种和质子化甲醛可以反应,和非质子化甲醛也可以反应,那么甲醛质子化的比例必然会影响这两个途径的比重。
而且即便不需要,算一下甲醛质子化这一步也有利于读者对反应整体有一个更好的把握。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

29

帖子

0

威望

133

eV
积分
162

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2022-3-26 19:19:09 | 只看该作者 Only view this author
wzkchem5 发表于 2022-3-26 01:57
那可以不算,但是也要留意一下质子化甲醛的浓度是否会影响后续反应的分支比,比如有的物种和质子化甲醛可 ...

emm你说的不算是指的是甲醛自身歧化生成甲酸这一步不用算吗?那我模拟甲醛质子化这一步直接添加一个氢离子吗

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

14#
发表于 Post on 2022-3-26 19:44:13 | 只看该作者 Only view this author
xuexi123. 发表于 2022-3-26 12:19
emm你说的不算是指的是甲醛自身歧化生成甲酸这一步不用算吗?那我模拟甲醛质子化这一步直接添加一个氢离 ...

如果甲酸是甲醛歧化产生的话,那还是需要多考虑一些事情,比如酸的浓度是随时间变化的,这个对你的结论会不会有影响需要考虑。我之前是假设酸的浓度是固定的
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

29

帖子

0

威望

133

eV
积分
162

Level 3 能力者

15#
 楼主 Author| 发表于 Post on 2022-3-27 00:42:29 | 只看该作者 Only view this author
wzkchem5 发表于 2022-3-26 19:44
如果甲酸是甲醛歧化产生的话,那还是需要多考虑一些事情,比如酸的浓度是随时间变化的,这个对你的结论会 ...

酸的浓度对于用gaussian模拟来说也会有影响嘛?

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 12:51 , Processed in 0.186608 second(s), 26 queries , Gzip On.

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