计算化学公社

标题: Warnings in Gaussian frequency calculation [打印本页]

作者
Author:
ymygc    时间: 2018-2-9 04:53
标题: Warnings in Gaussian frequency calculation
我在计算两条parallel reaction pathways的决速步的吉布斯自由能,想要由此预算实验上不同产物的比例,所以需要精确的吉布斯自由能。command line如下:
# opt=(calcfc,tight,ts,noeigentest) freq m06/6-311G(d,p) int=ultrafine Temperature=423.15

高斯输出文件里有两个警告(加粗处),请问这警告会多大程度影响吉布斯自由能的精确性呢?如此影响较大,该怎样修正?


    Eigenvalues --  ******************************
           X            0.99981  -0.01461   0.01301
           Y            0.01387   0.99837   0.05544
           Z           -0.01380  -0.05525   0.99838
This molecule is an asymmetric top.
Rotational symmetry number  1.
Warning -- assumption of classical behavior for rotation
           may cause significant error
Rotational temperatures (Kelvin)      0.00709     0.00467     0.00405
Rotational constants (GHZ):           0.14773     0.09724     0.08445
    1 imaginary frequencies ignored.
Zero-point vibrational energy    1222745.6 (Joules/Mol)
                                  292.24322 (Kcal/Mol)
Warning -- explicit consideration of  66 degrees of freedom as
           vibrations may cause significant error



作者
Author:
sobereva    时间: 2018-2-9 13:49
这是提醒你,对于那些低频振动模式,当做谐振子方式考虑可能对计算的热力学量产生明显误差,改用诸如受阻转子、自由转子等方式考虑可能会得到明显改善

但一般没必要太在意这个,比如是个大分子反应,反应的地方只涉及局部很小区域,不怎么牵扯柔性区域构象变化,那么那些地方计算前后误差会充分抵消
作者
Author:
ymygc    时间: 2018-2-9 23:31
sobereva 发表于 2018-2-9 13:49
这是提醒你,对于那些低频振动模式,当做谐振子方式考虑可能对计算的热力学量产生明显误差,改用诸如受阻转 ...

谢谢老师!
作者
Author:
devilove    时间: 2021-10-3 15:19
sobereva 发表于 2018-2-9 13:49
这是提醒你,对于那些低频振动模式,当做谐振子方式考虑可能对计算的热力学量产生明显误差,改用诸如受阻转 ...

可是sob老师,这里不是提示“ 1 imaginary frequencies ignored.” 存在虚频嘛?也不用在意影响嘛
作者
Author:
wzkchem5    时间: 2021-10-3 16:13
devilove 发表于 2021-10-3 08:19
可是sob老师,这里不是提示“ 1 imaginary frequencies ignored.” 存在虚频嘛?也不用在意影响嘛

如果这个虚频是平动/转动的频率,就可以不用考虑,否则需要尽量消掉。一般说的虚频必须消掉指的是振动的频率
作者
Author:
QJW    时间: 2025-4-25 01:01
wzkchem5 发表于 2021-10-3 16:13
如果这个虚频是平动/转动的频率,就可以不用考虑,否则需要尽量消掉。一般说的虚频必须消掉指的是振动的 ...

佬,怎么看是什么频率呢?

作者
Author:
wzkchem5    时间: 2025-4-25 09:46
QJW 发表于 2025-4-25 01:01
佬,怎么看是什么频率呢?

搜最近的translation、rotation两个词分别出现在哪,然后对输出文件做阅读理解,不难发现你这个频率是不是平动或振动频率




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3