计算化学公社

标题: 如何利用dmol频率计算给出的热力学参数得到分子配分函数 [打印本页]

作者
Author: Beautygrape 时间: 2015-1-6 14:50
标题: 如何利用dmol频率计算给出的热力学参数得到分子配分函数
我算的是金属富勒烯体系，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
具体该怎么算呢？

作者
Author: sobereva 时间: 2015-1-6 15:07
你可以搜搜有没有现成的程序，没有的话就只能自己写，公式在此文都有
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=123

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

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

作者
Author: Beautygrape 时间: 2015-1-6 16:28
C:\Documents and Settings\Administrator\桌面\1.bmp

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)