| 查看: 2172 | 回复: 11 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
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 ] |
回复此楼» 猜你喜欢
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
-
钴酸锂半电池小倍率容量上不去
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有206人回复
-
吉林大学材料物理本科生求问调剂信息
已经有24人回复
-
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有0人回复
-
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes
已经有2人回复
-
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事
已经有0人回复
-
cif求助下载
已经有1人回复
-
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有2人回复
Jldbaobei
木虫 (著名写手)
按时毕业找到好工作
- 应助: 2 (幼儿园)
- 金币: 3294.4
- 散金: 310
- 帖子: 1019
- 在线: 458.6小时
- 虫号: 507695
- 注册: 2008-02-19
- 专业: 燃烧学
9楼2010-04-10 21:33:50
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+1): 2010-03-06 21:15
小木虫(金币+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楼2010-03-05 21:08:37
Jldbaobei
木虫 (著名写手)
按时毕业找到好工作
- 应助: 2 (幼儿园)
- 金币: 3294.4
- 散金: 310
- 帖子: 1019
- 在线: 458.6小时
- 虫号: 507695
- 注册: 2008-02-19
- 专业: 燃烧学
3楼2010-03-05 21:25:26
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+10): 2010-03-06 10:22
小木虫(金币+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 ] |
4楼2010-03-05 23:19:22