求助，用cbs-qb3和b3lyp计算碳碳键断裂能垒，结果弄不明白

xuejunfei1987

文献上用CBS-QB3计算HOCH2CH2O的β碳碳键解离能垒为9.2kcal/mol,计算HOCH2CH2CH2O的β碳碳键解离能垒为11.9kcal/mol， 我用B3LYP/6-311++G(d,p)计算HOCH2CH2O的β碳碳键解离能垒为8.4kcal/mol，而HOCH2CH2CH2O的β碳碳键解离能垒为13.9kcal/mol，为什么HOCH2CH2O的能垒降低了，而HOCH2CH2CH2O[size=14.6667px]的解离能垒升高了呢，按理说应该要降都降啊，不知道还要考虑什么？希望高人给与指点

xuejunfei1987

好的，谢谢老师耐心解答，那我算过渡态和结构优化的时候是不是也要用这个基组啊

sobereva

xuejunfei1987 发表于 2019-4-2 09:13
谢谢sob老师答复，目前我们实验室都用的这种基组，计算条件有限，就算是很烂，我觉得我算个解离能垒还不 ...

正因为计算条件有限，所以才需要用更恰当的级别，避免浪费有限的计算资源
用M06-2X/def-TZVP比你现在这级别更便宜，结果反倒还好得多，何故不用？而且这么小的体系，别说结合def-TZVP了，就算结合def2-TZVP，一般的四核机子都根本没有压力。
B3LYP算有机体系热力学数据比M06-2X差十万八千里，这年头很容易遭到对理论懂行的人的批评，尤其是你算的数据是比较关键的数据的时候。
别再浪费计算资源了，仔细看我上面的博文。是时候刷新观念了，B3LYP算热力学数据是猴年马月的事了

“M06-2x不是算三重激发态吗”完全不懂你的意思，貌似你对M06-2X有天大的误解

xuejunfei1987

sobereva 发表于 2019-4-2 08:21
你用的B3LYP/6-311++G(d,p)计算级别太渣。这种计算级别在现在来看是极度糟糕的选择。这么点体系起码得用M06 ...

谢谢sob老师答复，目前我们实验室都用的这种基组，计算条件有限，就算是很烂，我觉得我算个解离能垒还不至于出错的，算过很多了。和文献里面的对象是对应的，不会出错，文献上就是列表格说的解离能垒。
老师我还要补充一下，那个M06-2x不是算三重激发态吗，我这个是二重态

sobereva

你用的B3LYP/6-311++G(d,p)计算级别太渣。这种计算级别在现在来看是极度糟糕的选择。这么点体系起码得用M06-2X/def2-TZVP。加弥散函数完全没有任何意义，纯粹是白浪费时间，什么时候该加仔细看
谈谈弥散函数和“月份”基组
http://sobereva.com/119

该用什么泛函看
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）

另外，极有可能你的计算流程不对，和文献里取的CBS-QB3计算的对象不对应。CBS-QB3是专门算热力学数据的，你自行计算也得以正确方式算对应的热力学数据