计算化学公社

标题: 较大体系几何优化速度vs准确率和方法选择的问题 [打印本页]

作者
Author: Shannon 时间: 2014-12-19 07:33
标题: 较大体系几何优化速度vs准确率和方法选择的问题
各位老师们，对于如何优化大体系几何结构，我有一点不明白的地方。
我有一个50原子左右的体系，想优化得到比较准确的结构，然后继续算单点能量和布居分析，用什么方法得到的结构会特别省时间，而且比较准确呢？
结构已经通过PM6预优化了一次，目前在几个选着中纠结。
STO-3G* 和 3-21G*哪一个会比较好？不太敢去删掉极化函数，体系中含有N和O原子。用single zeta的基组对结构影响大吗？
不含 HF交换成分的纯粹的 DFT-GGA 和 HF比较，那一个会比较快呢？我感觉DFT-GGA+外加DFT-D3似乎会比较快，因为HF的exchange部分没办法density fit，耗时是 O(N^4)。但不是很确定。

作者
Author: sobereva 时间: 2014-12-19 08:21
STO-3G*不算标准基组，也不建议用。3-21G*仅对于第三周期的才有极化函数，比如S、Cl，对于N、O是没有极化函数的。
强烈建议用ORCA在GGA级别下结合def2-SVP来计算，由于RI技术，速度远比高斯快得多得多。而且这个基组的质量也比6-31G**好很多（虽然都是2-zeta+极化）。
对于高斯，GGA和HF速度差不多，高斯的密度拟合对GGA的速度没有质的提升。但对于orca，由于库仑拟合，用GGA远比HF或杂化泛函快得多（几倍）。
但HF交换项在ORCA里也可以用COSX来加速计算，因此杂化泛函的计算速度也比高斯快不少。

这方面的问题在此帖有详细讨论
大体系弱相互作用计算的解决之道
http://sobereva.com/214

作者
Author: Shannon 时间: 2014-12-19 08:49
谢谢sob，我这就去试一试。

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