如何利用dmol频率计算给出的热力学参数得到分子配分函数

Beautygrape

我算的是金属富勒烯体系，dmol中频率计算会给出一些热力学参数如下：
   H,Trans:                0.889 kcal/mol
   H,Rot  :                0.889 kcal/mol
   H,pV   :                0.592 kcal/mol
   H,Vib - ZPVE:          21.782 kcal/mol
   S,Trans:               46.987  cal/mol.K
   S,Rot  :               39.440  cal/mol.K   with Symmetry Number:   1
   S,Vib  :              114.768  cal/mol.K
   C,Trans:                4.968  cal/mol.K
   C,Rot  :                2.981  cal/mol.K
   C,Vib  :              183.785  cal/mol.K
   ===================================
   H,Total - ZPVE:        24.152 kcal/mol
   S,Total:              201.194  cal/mol.K
   C,Total (p):          191.734  cal/mol.K
   G,Total:              -35.834 kcal/mol

我想得到分子配分函数，就像高斯中那样
                      Q            Log10(Q)             Ln(Q)
Total Bot       0.171891-206       -206.764748       -476.093426
Total V=0       0.142851D+26         25.154883         57.921258
具体该怎么算呢？

sobereva

你可以搜搜有没有现成的程序，没有的话就只能自己写，公式在此文都有
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=123

卡开发发

上次木虫上面和你说过了，除了振动项是谐振子求和，其余的就是理想气体的统计公式，找本统计（热）力学就能找到这些公式。
即便如此，就算完全无基础不知道如何下手，MS上也已经给出了参考文献Hirano, T., in MOPAC Manual, Seventh Edition, Stewart, J. J. P., Ed. (1993).
就可以查找这个文献，详见附件6.10.2的部分，给出了涉及到的公式，最后再根据配分函数分解的规则换算总配分函数即可。

sobereva

如果实在懒得写，可以用kisthelp，把分子的信息、振动频率数据写成kisthelp的输入文件格式，然后它就直接能给出热力学数据和配分函数。
kisthelp介绍见
使用KiSThelP结合Gaussian基于过渡态理论计算反应速率常数
http://sobereva.com/246

Beautygrape

C:\Documents and Settings\Administrator\桌面\1.bmp

Beautygrape

这是个好方法，只是当遇到开壳层的时候总是提示while reading electronic degeneracy in file .....  2, not an integer,刚才想上传照片，谁知道居然没成功，也不知道怎么传图片