计算化学公社

标题: 求助-关于orca的溶剂化模型 [打印本页]

作者
Author: Somnus荣小荣 时间: 2019-5-6 09:38
标题: 求助-关于orca的溶剂化模型
求助各位老师，大神~最近在算亚硫酸根小分子的液相反应，使用M062X方法，SMD溶剂化模型进行优化，频率计算，后选择CCSD(T)-F12方法，orca3.03程序进行更高级别的单点能校正，但是发现经M062X计算的反应自由能变和CCSD(T)-F12方法校正后的自由能变相差十几kcal，我怀疑是orca中的SMD溶剂化模型有问题，使用smd模型的关键词是，！cosmos"water" %cosmos smd true 有些不明白为何使用smd模型要通过cosmos来实现，而不是直接写入！smd“water” ，后经查手册发现是因为smd模型中的静电贡献来自cosmos。在orca4.0以后的版本中，使用smd模型的关键词是！cpcm"water" %cpcm smd true ，此时静电贡献来自于cpcm，而后我也进行了气相下M062X方法优化，频率计算和CCSD(T)-F12方法校正，发现校正前后自由能变相差4kcal左右，这个值应该是正常范围的，所以我觉得这可能是因为orca中的smd模型中的静电贡献不是来自smd而是来自cosmos或cpcm，所以造成高斯的M062X方法和orca的CCSD(T)-F12方法校正后，自由能变相差很大的原因，我这样理解是否正确，？

作者
Author: sobereva 时间: 2019-5-6 09:42
原本的SMD的静电部分是基于IEFPCM计算的，非静电部分是SMD方法自己定义的。然而ORCA不支持IEFPCM这么好的计算静电部分的溶剂模型，只能支持更糙的COSMO（从ORCA4开始已经废弃）和CPCM，所以ORCA里用SMD时只能基于COSMO或CPCM算SMD溶剂模型中的静电部分。显然原理上没有Gaussian的SMD严格，但也还算可以用。

作者
Author: Somnus荣小荣 时间: 2019-5-6 16:59
嗯，明白了，谢谢sob老师！还想请问一下，有哪个程序可以支持CCSD(T)-F12方法使用SMD（IEFPCM计算静电部分）溶剂化模型吗？或者如果我使用高斯进行CCSD(T)-aug-cc-pvqz，进行更高级别的单点校正合理吗，优化和频率计算使用的是M062X/6-311++G(3df，3pd)，体系大小就是亚硫酸根这样的离子。

作者
Author: sobereva 时间: 2019-5-7 08:00

Somnus荣小荣发表于 2019-5-6 16:59
嗯，明白了，谢谢sob老师！还想请问一下，有哪个程序可以支持CCSD(T)-F12方法使用SMD（IEFPCM计算静电部分 ...

“CCSD(T)-F12结合溶剂模型”并没有什么意义
如下文所说，标准的计算溶剂下自由能的方法是单点在气相下计算，溶解自由能单独用M05-2X/6-31G*结合SMD来算
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html）

作者
Author: Somnus荣小荣 时间: 2019-5-7 09:35
好嘞，明白了，多谢sob老师！

作者
Author: yaoshi1206 时间: 2020-10-27 20:32
请教一问题，亚硫酸H2SO3的基态是单线态还是三线态？
可以从什么书上获得相应知识。

作者
Author: wzkchem5 时间: 2020-10-27 21:14

yaoshi1206 发表于 2020-10-27 20:32
请教一问题，亚硫酸H2SO3的基态是单线态还是三线态？
可以从什么书上获得相应知识。

建议学习一下Lewis理论，以及共振式的画法。
常见主族化合物的基态自旋态，一般画出共振式就知道了，个别少数情况下还需要用分子轨道理论才能判断。
还是自己学习、掌握判断方法比较好，不要问单独的某个分子的自旋态，不然要么就是遇到啥都要问人，要么就是做不恰当的类推，判断错。

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