计算化学公社

标题: 求助：团簇原子数过多用高斯好像算不动怎么办 [打印本页]

作者
Author: Celia_Chen 时间: 2024-7-12 10:14
标题: 求助：团簇原子数过多用高斯好像算不动怎么办
如题，正在计算一个具有18个锆原子的金属团簇分子，一共有570个原子，使用b3lyp/def2svp基组，但是发现算不动（使用的是学院64核的超算）；
此分子结构如下：
(, 下载次数 Times of downloads: 17)

在听取同组师兄的建议之后，将正丁氧基改成正甲氧基后，有376个原子，仍然使用b3lyp/def2svp基组，但是发现还是算不动；

分子结构如下：
(, 下载次数 Times of downloads: 12)
再然后，尝试使用更小的基组，C、H、O使用6-31G**基组，Zr使用def2svp基组，然后还是算不动...
在这里求助各位大佬，请问遇到这种情况应该怎么办，打扰了。

作者
Author: dzdhp 时间: 2024-7-12 12:46
你要不要试试先用半经验方法初步优化，然后在用orca，这么大的体系orca要快很多

作者
Author: sobereva 时间: 2024-7-12 13:34
300多个原子，B3LYP结合2-zeta基组用Gaussian在像样的双路服务器上对此体系优化能算得动，只不过时间很长而已。不重要的地方可以用更小的基组节约耗时（混合基组。另外，当前情况氢根本没必要带极化）。虽然改用ORCA靠纯泛函做优化每一步会更快，但是之后振动分析耗时巨高，而且优化达到同样收敛限所需要的步数很可能明显更多。

如果仅仅需要优化结构而之后不需要做振动分析，用scf=conver=7结合opt=loose能节约优化过程可能一半的时间。

对于C、O等，6-31G**和def2-SVP耗时基本没差别，甚至后者的原高斯函数还更少一丁点，明显 “尝试使用更小的基组，C、H、O使用6-31G**基组” 是错误的描述。

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