计算化学公社

标题: 高斯计算出错Unable to choose the S8 parameter [打印本页]

作者
Author: 金斗鱼 时间: 2022-7-7 20:38
标题: 高斯计算出错Unable to choose the S8 parameter
各位大佬好，最近算任务遇到了以下的问题，输入关键字为#P HF/6-31G(d) em=GD3BJ opt freq polar，输出会显示“Unable to choose the S8 parameter”，是色散校正的使用不当嘛，不知道如何修改，请各位指教！gaussian的版本为 Gaussian 09 D.01，输入和输出的截图及文件放在下面了，先谢谢大家！ (, 下载次数 Times of downloads: 21) (, 下载次数 Times of downloads: 15)

(, 下载次数 Times of downloads: 1)

(, 下载次数 Times of downloads: 0)

作者
Author: chands 时间: 2022-7-7 21:00
本帖最后由 chands 于 2022-7-7 21:28 编辑

用HF干啥呢，用DFT不香吗？算极化率要加弥散。简谈量子化学计算中DFT泛函的选择 http://sobereva.com/272

作者
Author: wzkchem5 时间: 2022-7-7 23:46
可能是高斯不支持对HF添加色散校正。但是本来一般也没必要用HF，因为B3LYP、M06-2X等等常见泛函基本上永远比HF好。

作者
Author: sobereva 时间: 2022-7-8 04:26
纯属胡算，以这样的计算发文章碰上稍有理论计算常识的审稿人都会被秒拒。

polar任务甭和其它任何任务一起用

HF优化出的结构是垃圾。好好看此文了解点最基本的方法选用常识
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）

Gaussian里没有自带给HF用的D3参数

6-31G*算极化率是极度荒诞的，好好看下文了解能用的下限是什么
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

作者
Author: 金斗鱼 时间: 2022-7-8 09:01

sobereva 发表于 2022-7-8 04:26
纯属胡算，以这样的计算发文章碰上稍有理论计算常识的审稿人都会被秒拒。

polar任务甭和其它任何任务一 ...

感谢sob哈，我其实是算溶剂化自由能的，是别人的研究中比较了用HF、B3LYP、M062X和M052X算的结果，最好的结果就是用HF优化结构和能量计算得来的。别人的研究和我的相似，只是换了物质，我就进行了尝试。极化率这个错误确实不太懂，因为目前也不用极化率，想着以后可能用到，所有的任务我都带上了polar，我再去好好学习一下吧

作者
Author: 金斗鱼 时间: 2022-7-8 09:02

chands 发表于 2022-7-7 21:00
用HF干啥呢，用DFT不香吗？算极化率要加弥散。简谈量子化学计算中DFT泛函的选择 http://sobereva.com/272

是别人的研究中比较了用HF、B3LYP、M062X和M052X算的结果，最好的结果就是用HF优化结构和能量计算得来的。别人的研究和我的相似，只是换了物质，我就进行了尝试。

作者
Author: 金斗鱼 时间: 2022-7-8 09:04

wzkchem5 发表于 2022-7-7 23:46
可能是高斯不支持对HF添加色散校正。但是本来一般也没必要用HF，因为B3LYP、M06-2X等等常见泛函基本上永远 ...

看来是这样了

作者
Author: wzkchem5 时间: 2022-7-8 14:41

金斗鱼发表于 2022-7-8 02:02
是别人的研究中比较了用HF、B3LYP、M062X和M052X算的结果，最好的结果就是用HF优化结构和能量计算得来的 ...

这是因为误差抵消的原因。因为HF从理论上讲没有任何理由比DFT好，所以文献发现HF更好必然是因为其他地方引入了一个误差（比如基组不够大），HF引入了相反符号的误差，导致误差抵消。对于你的物质，不能保证这两个误差仍然抵消。所以还是要选用合适的泛函，搭配合适的基组。
如果只因为文献用了一个理论级别发现不错，不理解背后原因就擅自拿来用的话，做出来结果没人信的，一定会被审稿人毙掉。

作者
Author: sobereva 时间: 2022-7-9 03:57

金斗鱼发表于 2022-7-8 09:01
感谢sob哈，我其实是算溶剂化自由能的，是别人的研究中比较了用HF、B3LYP、M062X和M052X算的结果，最好的 ...

你拿这么过时、被学术界可以称得上已经被充分唾弃的方法去计算，不被批才怪。SMD原文表13里就有HF/6-31G*和M05-2X/6-31G*的对比，后者整体精度明显强于前者。另一种比较不错的溶剂模型uESE（http://sobereva.com/593）也没有干出用HF优化这种事。
而且，就算某一篇文章里恰好发现用HF优化再算溶解自由能由于误差抵消的巧合而最终结果不错，也不代表这结论就能移植到你的体系。只有当研究者用大批量分子构造的测试集测试发现这现象是普遍的，而且你的体系基本特征和测试集里的相同，才能把文献里的做法搬过来。

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