liyuanhe211 发表于 2020-2-23 08:02
和BSSE没关系。

正向势垒如何？看看IRC曲线中（最好跑长一点）情况如何？虽然说有电子能垒无自由能垒本 ...

正向能垒（CH2OCO-H2-1 → CH2OCO-H2-TS）8.3 kcal/mol。我跑IRC是把H2分子的两个H定义了phase以后分成正向（文件名字带reverse的那个，去往末态）和逆向（文件名字带foward的那个，回到初态）两段同时跑的(stepsize=20,maxpoint=200)，已经截图加在原帖里了。正向（文件名字带reverse的那个，去往末态）那个是全部跑完的。但为了节省机时，而且逆向（右半边，文件名字带foward的那个forward）那段我看到能量曲线也已经变得相对平坦了，所以没有等全部跑完我就停掉了
用6-311+G(d)的基组的确不是最合理，开始这个课题的时候，仅仅是参照文献的方法，把文献用的B3LYP/6-311+G*换了方法。

和BSSE没关系。

正向势垒如何？看看IRC曲线中（最好跑长一点）情况如何？虽然说有电子能垒无自由能垒本身不算是个异常情况，但放在这个反应里直观上有点奇怪。

另外从基组匹配的角度上，你给烷基末端的碳都加了弥散，却不给直接参与反应的Co-H中的氢原子加弥散，虽然问题不大，而且跟你现在的问题没有直接关系，但不算非常匹配的基组。