计算化学公社

标题: 使用Gaussian提供的波函数做DLPNO-CCSD(T)计算时如何考虑溶剂模型 [打印本页]

作者
Author: zhouoh 时间: 2020-8-1 13:01
标题: 使用Gaussian提供的波函数做DLPNO-CCSD(T)计算时如何考虑溶剂模型
我计算的反应体系是在乙醚下进行的，但是我查手册发现好像Orca的溶剂模型并不支持乙醚，我想用Gaussian计算考虑SMD溶剂模型的HF波函数，之后传给Orca做电子相关部分的计算。请问这种方法可以吗？我在sob老师的博文中看到，精确的考虑溶剂的效应需要做外迭代，请问不考虑的话误差会很大吗？谢谢各位老师！

作者
Author: hebrewsnabla 时间: 2020-8-1 13:10
orca支持乙醚的，见手册9.35.3节

作者
Author: zhouoh 时间: 2020-8-1 13:29

hebrewsnabla 发表于 2020-8-1 13:10
orca支持乙醚的，见手册9.35.3节

谢谢，我看见支持乙醚了，请问一下这两个程序的溶剂模型会导致收敛的波函数不一样吗

作者
Author: hebrewsnabla 时间: 2020-8-1 13:39

zhouoh 发表于 2020-8-1 13:29
谢谢，我看见支持乙醚了，请问一下这两个程序的溶剂模型会导致收敛的波函数不一样吗

会。

作者
Author: Accelerator 时间: 2020-8-1 17:31
隐式溶剂模型就不是给post-HF直接用的把带着隐式溶剂化模型的HF波函数当做参考态波函数更是没有道理
正确的做法见http://sobereva.com/327

作者
Author: biogon 时间: 2020-8-1 17:33
不进行外迭代就不用带上溶剂模型了

作者
Author: zhouoh 时间: 2020-8-1 18:21

Accelerator 发表于 2020-8-1 17:31
隐式溶剂模型就不是给post-HF直接用的把带着隐式溶剂化模型的HF波函数当做参考态波函数更是没有道理
正确 ...

明白了，谢谢老师

作者
Author: sobereva 时间: 2020-8-2 03:50
ORCA是支持乙醚的。

你当前把Gaussian波函数传给ORCA的做法并没什么意义，直接在ORCA里结合SMD表现的乙醚做DLPNO-CCSD(T)就完了。此时是无法考虑外迭代，溶剂效应只是在HF级别下表现。

如果你是想要得到溶剂下的精确自由能，实际上最好的做法是用Gaussian做opt freq得到自由能热校正量，再在M05-2X/6-31G*下得到溶解自由能，再在ORCA下用DLPNO-CCSD(T)计算气相下精确电子能量，即下文说的
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html）
这种做法不仅结果理想而且也原理上也十分得当，没有任何别扭的地方。

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