计算化学公社

标题: 高斯单点算时不同级别计算结果相差大 [打印本页]

作者
Author: yuying348 时间: 2019-8-21 18:04
标题: 高斯单点算时不同级别计算结果相差大
找过渡态时气态结构优化采用用b3lyp/6-31g（d）-D3，找到了两个过渡态，它们的电子能，自由能差值都在1千卡以内。保存此结构采用m062X/6-311+g（d，p）进行单点算，溶剂化用smd，电子能差值到10千卡左右，不知道是什么原因？涉及原子有C H N O S Br。
请教大家帮忙。谢谢。

作者
Author: shalene 时间: 2019-8-21 19:32
用双杂再算算？

作者
Author: sobereva 时间: 2019-8-22 00:23
级别相差这么大，而且一个考虑溶剂效应一个没考虑，差10kcal/mol太正常了

作者
Author: yuying348 时间: 2019-8-22 00:33
本帖最后由 yuying348 于 2019-8-22 00:39 编辑

谢谢老师。现在用小机组气态的自由能值或者电子能计算出来的势垒和实验一致，但是大机组单点算后彻底无语过渡态振动趋势，跑irc都合理。我下一步该怎么做？
小机组计算的两个势垒都是14左右，（大机组单点算能量+小机组自由能校正）计算的势垒一是44，势垒二是34左右

作者
Author: 望月凝香 时间: 2019-8-22 06:23
尝试用b3lyp-d3/def2svp//def2tzvpp (smd)来计算。
使用同一个泛函，同样考虑溶剂效应，只提升基组级别。

作者
Author: yuying348 时间: 2019-8-22 17:25
本帖最后由 yuying348 于 2019-8-22 23:01 编辑

望月凝香发表于 2019-8-22 06:23
尝试用b3lyp-d3/def2svp//def2tzvpp (smd)来计算。
使用同一个泛函，同样考虑溶剂效应，只提升基组级别。

谢谢。采用您推荐的机组单点算，两个过渡态能量差别也是很小，符合想要的结果。但是，计算的体系含有C H N O S Br ,总共101个原子，请教您有采用推荐机组的文章推荐一下吗？

作者
Author: yuying348 时间: 2019-8-22 21:51

sobereva 发表于 2019-8-22 00:23
级别相差这么大，而且一个考虑溶剂效应一个没考虑，差10kcal/mol太正常了

谢谢。通过计算单点发现差别不是溶剂化引起的，是方法引起的。我计算的体系含有C H N O S Br ,总共101个原子，请教您对于我目前研究的体系采用b3lyp-D3大机组溶剂化单点算，可信吗？您推荐采用什么机组合适？

作者
Author: 望月凝香 时间: 2019-8-23 05:02

yuying348 发表于 2019-8-22 17:25
谢谢。采用您推荐的机组单点算，两个过渡态能量差别也是很小，符合想要的结果。但是，计算的体系含有C H ...

看sob关于基组选择的帖子，有引文说pople基组构筑质量的。
其中一个是预测了苯的最低能量构型非平面。

作者
Author: yuying348 时间: 2019-8-23 09:04
好的。谢谢您。

作者
Author: sobereva 时间: 2019-8-23 22:29

望月凝香发表于 2019-8-22 06:23
尝试用b3lyp-d3/def2svp//def2tzvpp (smd)来计算。
使用同一个泛函，同样考虑溶剂效应，只提升基组级别。

你把//符号前后弄反了
量化界的习俗是//后头的是优化，前头的是单点

作者
Author: sobereva 时间: 2019-8-23 22:30

yuying348 发表于 2019-8-22 17:25
谢谢。采用您推荐的机组单点算，两个过渡态能量差别也是很小，符合想要的结果。但是，计算的体系含有C H ...

仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

烂基组就算算出来貌似和实验对得更好的结果也毫无意义，纯属侥幸，发文章时候肯定被审稿人批

用M06-2X/def2-TZVP//B3LYP-D3(BJ)/6-311G*计算，SMD模型考虑溶剂效应（包括opt freq阶段）。不行再考虑双杂化，以及考虑模型/机理合理性问题

作者
Author: yuying348 时间: 2019-8-24 01:44
谢谢您。

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