计算化学公社

标题: 求助VRC-VTST输入文件中VIB选项如何设置 [打印本页]

作者
Author: liurui 时间: 2023-8-1 21:30
标题: 求助VRC-VTST输入文件中VIB选项如何设置
本帖最后由 liurui 于 2023-8-1 21:39 编辑

各位老师，目前关注的反应在考虑零点校正后为负势垒，目前在实验VRC-VTST方法进行计算，但是输入文件中有一项VIB，反应物和产物都有一个这个，猜测可能是平衡结构的频率结果，但是波数和a.u.之间进行换算并未得到E-03这样的数量级，想请问这些值是如何得到的？谢谢各位老师！非常感谢！
# Energy in a.u.
ENERGY -237.692205266
# Frequencies in a.u.; have been scaled
VIB
3.0336311764E-03 5.0499657596E-03 5.5367232386E-03
END

附件中附上了TESTRUN中给的例子，难道是附件中的gaussian输出的这个吗？

作者
Author: tinytiger98 时间: 2023-8-2 08:47
图片中的Q是配分函数。
au到cm-1之间你是如何转换到？

作者
Author: liurui 时间: 2023-8-2 10:29

tinytiger98 发表于 2023-8-2 08:47
图片中的Q是配分函数。
au到cm-1之间你是如何转换到？

老师，我尝试将振动频率cm-1转换单位到a.u.(1 a.u.=2.2055 cm-1),发现和例子中给的对应不上，查找输出文件只有配分函数的量级貌似对应，但是和输入文件要的貌似不是一个，谢谢老师。

作者
Author: dantevinsky 时间: 2023-8-2 11:29

liurui 发表于 2023-8-2 10:29
老师，我尝试将振动频率cm-1转换单位到a.u.(1 a.u.=2.2055 cm-1),发现和例子中给的对应不上，查找输出文 ...

1 a.u.=2.194746E5 cm-1。。。。。

作者
Author: liurui 时间: 2023-8-2 11:55

dantevinsky 发表于 2023-8-2 11:29
1 a.u.=2.194746E5 cm-1。。。。。

老师但是将振动频率换算也无法得到TESTRUN中给的算力中如下的量级E-03

# Frequencies in a.u.; have been scaled
VIB
3.0336311764E-03 5.0499657596E-03 5.5367232386E-03
END
谢谢老师！

作者
Author: dantevinsky 时间: 2023-8-2 12:32
两个六百多波数的，一个一千多波数的，这振动频率不是很正常吗。你自己的体系中几百到几千的波数的振动模式肯定也在这个量级啊

作者
Author: liurui 时间: 2023-8-2 14:33

dantevinsky 发表于 2023-8-2 12:32
两个六百多波数的，一个一千多波数的，这振动频率不是很正常吗。你自己的体系中几百到几千的波数的振动模式 ...

老师，我的意思是例子中给的量级在E-03和我算的几百几十不一样，是否这个不是频率呢？

作者
Author: dantevinsky 时间: 2023-8-2 15:33
这就是频率，毫无疑问，这是a.u.下的频率，换算过去就是几百cm-1。且，其实这里的vibs应该不会影响最终计算结果(按照VRC-VTST公式来看的话)。

作者
Author: liurui 时间: 2023-8-2 16:14

dantevinsky 发表于 2023-8-2 15:33
这就是频率，毫无疑问，这是a.u.下的频率，换算过去就是几百cm-1。且，其实这里的vibs应该不会影响最终计算 ...

好嘞，我去测试一下，谢谢老师！

作者
Author: Melvin 时间: 2023-8-3 00:40

dantevinsky 发表于 2023-8-2 15:33
这就是频率，毫无疑问，这是a.u.下的频率，换算过去就是几百cm-1。且，其实这里的vibs应该不会影响最终计算 ...

按楼上说的换算就行，你还纠结单位可以看Polyrate手册，里面有更改频率或者能量单位的关键词

作者
Author: adminer 时间: 2023-12-29 16:17
这个可以在fu61文件里面找到吧

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