24小时热门版块排行榜    

查看: 6510  |  回复: 25

herarysara

木虫 (小有名气)

[求助] 求gaussian计算反应速率常数,指前因子方法。

各位大虾请赐教:
    在站内找了计算指前因子的帖子,可是说得都不是太清楚,没怎么看明白,求各位前辈帮忙解释下,如何在得到过渡态振动频率后,计算指前因子和反应速率常数,主要是针对加成反应之类的多反应物体系,有什么较为有用的参考资料请大家推荐下,不胜感激!
回复此楼

» 收录本帖的淘帖专辑推荐

量化

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

wangth0921

铁虫 (职业作家)

★ ★ ★ ★ ★
小红豆(金币+5): 谢谢 2011-11-15 19:27:01
引用回帖:
10楼: Originally posted by herarysara at 2011-09-21 10:58:02:
你好,能详细告诉我一下具体在输出文件的那块查找,或者是不是计算频率还需要加入什么关键词吗?我在结果里面找不到相关的配分函数。

Sum of electronic and zero-point Energies=           -192.984385
Sum of electronic and thermal Energies=              -192.979453
Sum of electronic and thermal Enthalpies=            -192.978509
Sum of electronic and thermal Free Energies=         -193.011464

                     E (Thermal)             CV                S
                      KCal/Mol        Cal/Mol-Kelvin    Cal/Mol-Kelvin
Total                   58.222             15.948             69.361
Electronic               0.000              0.000              0.000
Translational            0.889              2.981             38.096
Rotational               0.889              2.981             23.983
Vibrational             56.445              9.986              7.281
Vibration  1             0.620              1.896              2.597
                       Q            Log10(Q)             Ln(Q)
Total Bot       0.111379D-27        -27.953195        -64.364611
Total V=0       0.285572D+13         12.455716         28.680345

这里Q就是配分函数,对应你所给的温度。“Total V=0”和"Total Bot"的差别是exp(ZPE/RT),即采用的零点不同。如果你要直接用Gaussian给出的数值,建议在输入文件中加入Scale Factor。配分函数的计算很简单,建议去看看物化书上统计热力学章节。

指前因子是个经验数值,可简单计算为kT/h*Q(TS)/(QA*QB)。更准确的方法是计算几个不同温度下的速率常数后拟合得到,可以对Arrenhenius公式或其他带T^n的公式来拟合。

» 本帖已获得的红花(最新10朵)

21楼2011-11-15 10:08:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

zzti314

银虫 (小有名气)

【答案】应助回帖

★
gmy1990(金币+1): 2011-09-19 18:31:05
好像用gaussian计算不了反应速率常数 指前因子
2楼2011-09-19 16:43:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

herarysara

木虫 (小有名气)

引用回帖:
2楼: Originally posted by zzti314 at 2011-09-19 16:43:23:
好像用gaussian计算不了反应速率常数 指前因子

我看很多文章中都有提到gaussian算的,可是都没说清楚。
3楼2011-09-20 08:10:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

herarysara

木虫 (小有名气)

没有人知道吗?最近各位版主前辈怎么都不在啊,郁闷~
4楼2011-09-20 08:11:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
5楼2011-09-20 08:35:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzti314

银虫 (小有名气)

【答案】应助回帖

引用回帖:
3楼: Originally posted by herarysara at 2011-09-20 08:10:48:
我看很多文章中都有提到gaussian算的,可是都没说清楚。

这个 我感觉是用gaussian算出构型 配分函数之类的  然后需要自己往后算 不能从gaussian中直接得到结果
6楼2011-09-20 10:22:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

herarysara

木虫 (小有名气)

引用回帖:
6楼: Originally posted by zzti314 at 2011-09-20 10:22:37:
这个 我感觉是用gaussian算出构型 配分函数之类的  然后需要自己往后算 不能从gaussian中直接得到结果

是这样的,gaussian怎么算配分函数呢?能详细解释下吗?谢了!
7楼2011-09-20 11:28:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

引用回帖:
7楼: Originally posted by herarysara at 2011-09-20 11:28:07:
是这样的,gaussian怎么算配分函数呢?能详细解释下吗?谢了!

做频率分析的时候,Gaussian自动给出所有的配分函数。
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
8楼2011-09-20 16:21:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

慢三儿

木虫 (小有名气)

【答案】应助回帖

我觉得按照活化能的定义,应该是要计算出反应的过渡态结构,其能量与反应物原料能量的差就是活化能,然后再根据阿伦尼乌斯方程计算反应速率常数吧?
9楼2011-09-20 21:52:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

herarysara

木虫 (小有名气)

引用回帖:
8楼: Originally posted by zhangmt at 2011-09-20 16:21:03:
做频率分析的时候,Gaussian自动给出所有的配分函数。

你好,能详细告诉我一下具体在输出文件的那块查找,或者是不是计算频率还需要加入什么关键词吗?我在结果里面找不到相关的配分函数。
10楼2011-09-21 10:58:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 herarysara 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +3 wangyl_123 2026-10-10 6/300 2026-10-11 15:04 by Abin0406
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 5/250 2026-10-11 09:00 by 好运来~9527
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见