应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】求助：高斯输出结果中的配分函数
121/2返回列表12下一页
查看: 2121  |  回复: 11
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

Jldbaobei

木虫 (著名写手)

按时毕业找到好工作

[交流] 【求助】求助：高斯输出结果中的配分函数已有2人参与

高斯输出结果：
                     E (Thermal)             CV                S
                      KCal/Mol        Cal/Mol-Kelvin    Cal/Mol-Kelvin
Total                  170.849             73.054            131.471
..................
.................
                      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.151066-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

最下面两个是平动配分函数和转动配分函数
请问在这个输出结果中哪一个是  或  如何得到振动配分函数？
前面的Vib(Bot)和Vib(V=o)分别是什么意思？


===================================================
非常感谢大家的回答，但是对高斯输出结果中
Total Bot       0.116086D-96
Total V=0       0.316193D+21
及
Vib (Bot)       0.151066-112      
Vib (Bot)  1    0.386265D+01         
具体指的是什么仍然不明白
我把 Vib (Bot)  1 ........ Vib (Bot)  31乘积算了下，发现结果和输出结果中的任何数值都不相等.........

请人告知上述符号具体代表什么

===================================================

似乎对高斯输出结果中的配分函数迷糊的不只我一个人
在以往的帖子里搜索亦没有找到答案

[ Last edited by Jldbaobei on 2010-4-14 at 18:48 ]
回复此楼

» 猜你喜欢
1楼2010-03-05 17:27:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

c111999

铜虫 (小有名气)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+10): 2010-03-06 10:22
1。震动佩芬函数数量级的差别来自
http://www.gaussian.com/g_whitepap/thermo.htm
the overall vibrational partition function is[McQuarrie, §4-7, Eq. 4.46]:
                                            e−Θv,K /2T
                                 qv =II----------------------------
                                          1 − e−Θv,K /T
e−Θv,K /2T

http://www.chem.arizona.edu/~salzmanr/480b/statt02/statt02.html
Characteristic vibrational temperatures are usually several thousands of Kelvins except for very "soft" or low frequency vibrational modes.

例如 Θv,K=10T， Exp[-10]=0.0000453999，
vib1*vib2*……vib31会产生极大差异


2。反应速率
假设 A+B=C+D
向前速率 r=k[A][B]， []浓度
则需要算k， k=k0Exp[-ΔEact/RT]
ΔEact=Eprod-Ereac，可以用H，G，A……
k0=qA*qB
如果有系数不为一，则（qA）^x
http://www.gaussian.com/g_whitepap/thermo.htm
Gaussian uses the bottom of the well as the zero of energy (BOT) to determine the other thermodynamic quantities。
如果算速率，建议用vib（bot）

3。qA，qB根据具体情况计算。
如在表面，则qtran将会比气态减少一个自由度。

[ Last edited by c111999 on 2010-3-5 at 23:28 ]
赞一下(3人) 回复此楼
4楼2010-03-05 23:19:22
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jldbaobei

木虫 (著名写手)

按时毕业找到好工作
再顶一顶
依然求助在计算化学反应速率常数时，高斯输出文件中振动配分函数该取哪一个
看到以前也有人问这个问题，但一直没有得到明确的回答
望知道的虫友告知
谢谢啦
赞一下(1人) 回复此楼
9楼2010-04-10 21:33:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

c111999

铜虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+1): 2010-03-06 21:15
看你的参考点（何处能量为零）
http://www.gaussian.com/g_whitepap/thermo.htm
Gaussian uses the bottom of the well as the zero of energy (BOT) to determine the other thermodynamic quantities, but also prints out the V=0 partition function. Ultimately, the only difference between the two references is the additional factor of Θv,K /2, (which is the zero point vibrational energy) in the equation for the internal energy Ev . In the expressions for heat capacity and entropy, this factor disappears since you differentiate with respect to
temperature (T ).
赞一下(2人) 回复此楼
2楼2010-03-05 21:08:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jldbaobei

木虫 (著名写手)

按时毕业找到好工作
请问能说的更清楚一点吗

比如在使用配分函数计算化学反应速率的时候过渡态的全配分函数为转动、平动和振动配分函数的乘积

平动、转动配分函数高斯结果中均直接给出了，那么振动配分函数该使用Vib(Bot)还是Vib(V=0)呢

这两个不同参考点的配分函数数量级怎么差这么多呢？
赞一下 回复此楼
3楼2010-03-05 21:25:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jldbaobei

木虫 (著名写手)

按时毕业找到好工作
多谢
那高斯输出结果中
Total Bot       0.116086D-96
Total V=0       0.316193D+21
指的是两个不同参量下所得总的振动配分函数还是什么能量呢？

Vib (Bot)       0.151066-112       又是什么？
Vib (Bot)  1    0.386265D+01          等是指每一个原子对应的配分函数？

新手入门，请多多指教
赞一下 回复此楼
5楼2010-03-06 15:35:18
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c111999

铜虫 (小有名气)
★ ★
yjcmwgk(金币+2): 2010-04-11 17:18
有兴趣的话试一试把所有的q都乘起来看看等于啥。
可以找个小分子，快。

下面瞎说
total = 总配分函数
vib = vib 总震动
vib1=vib1 频率一
赞一下(1人) 回复此楼
6楼2010-03-07 00:37:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

coolrainbow

木虫 (著名写手)

未来国家冻凉
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
你把本科物理化学的统计力学部分好好看看你就明白了
赞一下(1人) 回复此楼
技术博客：http://hi.baidu.com/coolrainbow/blog
7楼2010-03-07 13:18:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c111999

铜虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
请问coolrainbow：
  
如何计算反应速率常数，如 A2+B=C
k=k0Exp(-Eact/kT)
尤其是系数k0怎么算？
Eact是用 delta U（0k），H，A还是G？
看似k0与Eact的选择很有关，

物理化学上没找到
多谢
赞一下(1人) 回复此楼
8楼2010-03-07 14:07:28
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jldbaobei

木虫 (著名写手)

按时毕业找到好工作
再顶一下
希望明白的虫友明确告知
赞一下 回复此楼
10楼2010-04-12 21:30:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Jldbaobei 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭