24小时热门版块排行榜

返回列表

wangsihang

木虫 (正式写手)

应助: 74 (初中生)
金币: 2281.4
红花: 5
帖子: 649
在线: 185.3小时
虫号: 1261538
注册: 2011-04-11
专业: 催化化学

[求助] 利用Gaussian计算速率系数

请问大家按照下面的公式计算速率系数
k(T)=λ(kBT/h)(Q≠/QAQB)exp(-Ea/RT).
λ为量子轨道效应较正系数，其经验表达式为
λ=1+(hν≠/kBT)2/24.
式中，kB为玻尔兹曼常数，h为普朗克常数，Ea为反应活化能，ν≠为过渡态振动虚频；Q≠为过渡态的配分函数；QA与QB为反应物的配分函数；Q为全配分函数；为平动配分函数Qt、转动配分函数Qr和振动配分函数Qv的连乘积：Q=QtQrQv
根据Gaussian频率计算的输出结果，
                  Q          Log10(Q)          Ln(Q)
Total Bot    0.116086D-96       -96.935218    -223.201589
Total V=0    0.316193D+21       20.499952       47.202884
Vib (Bot)    0.151066D-112    -112.820834    -259.779570
Vib (Bot)  1 0.386265D+01       0.586886       1.351355
Vib (Bot)  2 0.342876D+01       0.535138       1.232200
...............
...............

Vib (Bot) 31 0.242768D+00       -0.614809       -1.415650
Vib (V=0)    0.411469D+05       4.614337       10.624903
Vib (V=0)  1 0.439488D+01       0.642947       1.480440
..................
Vib (V=0) 30 0.106140D+01       0.025879       0.059588
Vib (V=0) 31 0.105582D+01       0.023590       0.054318
Electronic    0.300000D+01       0.477121       1.098612
Translational 0.249756D+09       8.397517       19.335997
Rotational    0.102560D+08       7.010978       16.143372
Qv（振动配分函数）应该去哪个值呢？

回复此楼