解离焓中的RT如何算。

xieyi

算解离焓的时候，那个RT是什么怎么算的，我看到一篇中文文献有说R 为常数8 . 3 1 4 ；T 为温度，但是这个是怎么算的呢？虽然虽然自由基和母体的RT相互抵消，但是H的RT还在，不知道要怎么算？The thermal corrections contain the scaled ZPE + 3/2RT(translation) + 3/2RT (rotation) plus a small amount (vibrational excitation) for phenol and phenoxyl, and add RT to convert fromenergy to enthalpy (H =E + PV = E + RT). For the H atom,adding 3/2RT (translation) + RT (enthalpy conversion)。这篇英文文献说二分之五RT=0.00094hartree.

sobereva

R=8.3144648 J/mol/K，乘上T就是能量了。
不清楚你的意思，你就用比如高斯做freq任务计算每个片段和母体的焓，然后求差就行了，都用不着自己手算这项。

xieyi

哦哦，看到这个R的单位，我就明白了~~~加上RT似乎是能量转变成焓的意思。

取freq任务之后输出的焓就完事了（这是比较粗略的做法，更好的做法是用尽可能高精度的方法算单点能，再加上中或低精度算的热力学校正量，这样最划算）。
您上次在另一个帖子中给我的回复。
所以我采取您的建议，用单点能和零点能算，然后零点能再用校正因子校正一下。得到的就是文献上的羟基BDE算法

xieyi

sobereva 发表于 2015-8-3 12:03
R=8.3144648 J/mol/K，乘上T就是能量了。
不清楚你的意思，你就用比如高斯做freq任务计算每个片段和母体的 ...

老师，很奇怪，这篇文章里面设定的是298K。但是他却说二分之五RT是0.00094hartree（附复制过来的文章H-atom total energy =-0.50027hartree in all four basis sets and the correction 5/2RT = 0.00094 hartree）这是什么情况？

liyuanhe211

xieyi 发表于 2015-8-3 15:19
老师，很奇怪，这篇文章里面设定的是298K。但是他却说二分之五RT是0.00094hartree（附复制过来的文章H-at ...

RT = 8.314 J/mol`k * 298 k  = 2.5 kJ/mol (这个数最好记住) = 0.6 kCal/mol
1 Hartree = 2625.5 kJ/mol (这个数也最好记住，实际上是2625.499)
所以RT= 9.5E-4，和你说的一样，有什么问题吗？

xieyi

文献中写的是二分之五RT，不是RT。