使用Gaussian提供的波函数做DLPNO-CCSD(T)计算时如何考虑溶剂模型

zhouoh · 发表于 Post on 2020-8-1 13:01:06

我计算的反应体系是在乙醚下进行的，但是我查手册发现好像Orca的溶剂模型并不支持乙醚，我想用Gaussian计算考虑SMD溶剂模型的HF波函数，之后传给Orca做电子相关部分的计算。请问这种方法可以吗？我在sob老师的博文中看到，精确的考虑溶剂的效应需要做外迭代，请问不考虑的话误差会很大吗？谢谢各位老师！

hebrewsnabla · 发表于 Post on 2020-8-1 13:10:09

orca支持乙醚的，见手册9.35.3节

zhouoh · 发表于 Post on 2020-8-1 13:29:46

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

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

hebrewsnabla · 发表于 Post on 2020-8-1 13:39:33

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

会。

Accelerator · 发表于 Post on 2020-8-1 17:31:16

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

biogon · 发表于 Post on 2020-8-1 17:33:05

不进行外迭代就不用带上溶剂模型了

zhouoh · 发表于 Post on 2020-8-1 18:21:57

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

明白了，谢谢老师

sobereva · 发表于 Post on 2020-8-2 03:50:56

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）
这种做法不仅结果理想而且也原理上也十分得当，没有任何别扭的地方。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 使用Gaussian提供的波函数做DLPNO-CCSD(T)计算时如何考虑溶剂模型

浏览过的版块