Gaussian中利用pop=MK方式计算原子电荷能够接受的最低计算级别是什么

让你变成回忆

在Gaussian中采用pop=MK iop(6/33=2,6/41=10,6/42=17)的方式计算有机小分子的ESP原子电荷，由于需要计算的体系实在太多，所以想请教，最低的能够接受的计算级别大概是什么级别？HF计算出来的结果能用吗？

在sob老师的文章里面（Acta Phys. -Chim. Sin. 2012, 28 (1), 1-18）采用HF/6-31G*的计算级别，我这里采用同样的方法可以吗？

snljty

HF并不会比很多杂化泛函明显更快，在用了RI加速后更是显著慢于纯泛函。但是精度却非常低。
老的计算往往是还没有好的隐式溶剂模型，用HF高估的偶极矩近似代替溶剂模型，现在没有这个必要。
如果真的是非常大的体系，可以考虑用ORCA开RI使用纯泛函，然后导出gbw转成miolden，用Multiwfn读入molden的波函数，调用cubegen或者利用http://bbs.keinsci.com/thread-16499-1-1.html的方法计算拟合静电势，可以比Gaussian直接用HF快得多的多而且精度也好得多。

让你变成回忆

snljty 发表于 2020-6-23 20:42
HF并不会比很多杂化泛函明显更快，在用了RI加速后更是显著慢于纯泛函。但是精度却非常低。
老的计算往往是 ...

体系并不大，主要是计算的体系非常多。
而且是自己写的程序去调用高斯计算，所以目前主要还是考虑用高斯来计算。

thanhtam

想快就orca呗

snljty

让你变成回忆 发表于 2020-6-23 20:49
体系并不大，主要是计算的体系非常多。
而且是自己写的程序去调用高斯计算，所以目前主要还是考虑用高斯 ...

自己写程序调用orca和Multiwfn又不比Gaussian麻烦多少。

sobereva

那篇文章用HF/6-31G*只是作为对比电荷计算方法的目的，不是为了电荷有多准确。实际计算显然不能用HF/6-31G*。

对于拟合静电势电荷，基组建议不低于6-311G**或def-TZVP。如果是Gaussian，泛函用B3LYP或PBE0
如果计算大量体系，计算资源又有限，用ORCA优化和计算波函数，结合Multiwfn调用cubegen算拟合静电势电荷绝对是最划算的，只不过是shell脚本里加几行语句的事。此时可以用又便宜又不错的B97-3c，要求更高可以用PBEh-3c。

相关信息见
详谈Multiwfn支持的输入文件类型、产生方法以及相互转换
http://sobereva.com/379（http://bbs.keinsci.com/thread-6020-1-1.html）
Multiwfn现已可以调用cubegen使静电势分析耗时有飞跃式的下降！
http://sobereva.com/435（http://bbs.keinsci.com/thread-10682-1-1.html）


后记：如今已经不推荐调用cubegen了，基于Multiwfn自己的静电势代码计算拟合静电势电荷比调用cubegen还更快，一定要看下文
Multiwfn使用的高效的静电势算法的介绍文章已于PCCP期刊发表！
http://sobereva.com/614（http://bbs.keinsci.com/thread-24940-1-1.html）

让你变成回忆

sobereva 发表于 2020-6-24 16:05
那篇文章用HF/6-31G*只是作为对比电荷计算方法的目的，不是为了电荷有多准确。实际计算显然不能用HF/6-31G* ...

sob老师您好，计算的体系不算大，主要是体系非常多。