求助 计算自由能时过渡态比产物自由能低

fanmingren

单点能计算级别： # b2plyp/def2tzvpp
几何优化频率计算级别：#p opt freq 6-31+g(d,p) m062x
IRC:#p irc(calcall,maxpoint=20,noeigen) 6-31+g(d,p) m062x

恳请大佬求助

fanmingren

sobereva 发表于 2024-7-12 23:39
认真看此文
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/ ...

抱歉社长，我这个是推算反应机理的最后一步，两端的结构我都优化好了，单点能计算完毕后用shermo计算的自由能。excel表是根据初始态计算的自由能差值，单点能，结构优化以及IRC我都放在压缩文件里了。就是依据IRC和虚频判断都是存在过渡态的，但是自由能计算的结果却是产物的自由能比过渡态的自由能还要高

fanmingren

kimariyb 发表于 2024-7-12 23:01
这么抽象，直接用中文把计算过程描述清楚不就行了，贴这么多没用的图，一个不就够了，还有 Excel 也不知道 ...

抱歉，我这个是推算反应机理的最后一步，两端的结构我都优化好了，单点能计算完毕后用shermo计算的自由能。excel表是根据初始态计算的自由能差值，单点能，结构优化以及IRC我都放在压缩文件里了

sobereva

认真看此文
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620（http://bbs.keinsci.com/thread-25787-1-1.html）

kimariyb

这么抽象，直接用中文把计算过程描述清楚不就行了，贴这么多没用的图，一个不就够了，还有 Excel 也不知道是要表达什么。IRC 跑出来的结果不就是反应的路径。你直接拿两个独立的分子的能量加起来，能量比稳定复合物、过渡态的能量高也很正常。你要算此反应的能垒，应该要优化 IRC 两端，并且算单点。而不是用两个独立的分子。