计算化学公社

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

[新手求助] 计算自由能差的问题

[复制链接 Copy URL]

230

帖子

0

威望

2653

eV
积分
2883

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 一颗赛艇 于 2020-8-17 13:28 编辑

最近干活的时候遇到了个问题,反应方程如下:

A + B -> AB -> TS

A和B先形成复合物AB,然后是TS。因为要算KIE,用到速率常数所以要算自由能差。问题是复合物AB有很多构像,想准确点的话我打算通过跑IRC把 AB -> TS 路径上的所有反应相关的全找出来然后用Boltzmann分布平均一下。但是这样太麻烦,我们老板说:你就用TS的自由能直接减A与B的自由能之和。

我心想这不对吧,不能这么减吧?毕竟AB和TS的自由能包括3N-6个振动,A+B的则相当于3N-12

请问各位这么做对不对?错的话该怎么做?谢谢!

-------------------------------------------

另:
以前在付献彩的物理化学书上的热力学数据表找了个例子:NO2 和 N2O4 的 标准自由能。 我用M06-2X,无溶剂和使用频率校正因子的情形下按照上述的“直接”减法算了一下自由能差,发现和实验吻合得很好。也许这个不具有代表性?现在真是一头雾水

1万

帖子

0

威望

7400

eV
积分
18164

Level 6 (一方通行)

2#
发表于 Post on 2020-8-17 18:08:22 | 只看该作者 Only view this author
取决于实验是怎么做的。
如果实验是通过测量A和B的消耗来测的反应速率(且AB的稳态浓度远低于A和B的浓度),那么应该用TS-A-B
如果实验是通过测量复合物AB的消耗来测的反应速率,那么应该用TS-AB
如果实验是通过测量产物的生成来测的反应速率,那么需要根据AB的稳态浓度确定,如果AB稳态浓度很低,那么用TS-A-B,如果AB稳态浓度很高(就是A和B先几乎全部复合成AB再往下反应),那么用TS-AB
从数值上讲,这两种算法的唯一区别就是差了A+B->AB的平衡同位素效应EIE,如果AB是非共价复合物,那么EIE会很接近1,所以怎么算问题都不大。但是如果是算活化Gibbs自由能本身,就需要特别注意是从AB为起点还是从A+B为起点算了。

评分 Rate

参与人数
Participants 3
eV +11 收起 理由
Reason
一颗赛艇 + 5 谢谢
喵星大佬 + 3
thanhtam + 3

查看全部评分 View all ratings

BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

164

帖子

0

威望

1939

eV
积分
2103

Level 5 (御坂)

3#
发表于 Post on 2020-8-17 18:19:06 | 只看该作者 Only view this author
不需要管振动模式个数的问题

230

帖子

0

威望

2653

eV
积分
2883

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2020-8-18 09:18:24 | 只看该作者 Only view this author
thanhtam 发表于 2020-8-17 18:19
不需要管振动模式个数的问题

在宏观下N很大,无论是6还是12都无所谓。那么在微观下也不用管是为什么呢?谢谢

164

帖子

0

威望

1939

eV
积分
2103

Level 5 (御坂)

5#
发表于 Post on 2020-8-18 10:26:52 | 只看该作者 Only view this author
一颗赛艇 发表于 2020-8-18 09:18
在宏观下N很大,无论是6还是12都无所谓。那么在微观下也不用管是为什么呢?谢谢

抱歉,我没有想过这个问题的合理性,一直认为是常识。有什么物理原理规定要振动模式数量相同才能把自由能相减吗?

1万

帖子

0

威望

7400

eV
积分
18164

Level 6 (一方通行)

6#
发表于 Post on 2020-8-18 11:05:15 | 只看该作者 Only view this author
thanhtam 发表于 2020-8-18 10:26
抱歉,我没有想过这个问题的合理性,一直认为是常识。有什么物理原理规定要振动模式数量相同才能把自由能 ...

确实没有这种限制,不然如果关心的就是A+B->AB这个反应本身,岂不是没法算了
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 12:01 , Processed in 0.170878 second(s), 21 queries , Gzip On.

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