求助长链有机金属化合物SCF计算的基组选择

lujun

各位专家老师：
  大家好！向大家请教个问题：我目前处理一个有机金属化合物（见附图），有机配体与过渡金属Ir的配比为3：1，有机配体是一个长链的有机化合物，含C，H，O，P。整个有机金属化合物的原子数为180个，所以我对C，H，O，P选择cc-pVDZ级别的基组，Ir选择带有赝势的cc-pVDZ级别基组做结构优化，但是在Gaussian里面，SCF的计算速度非常慢。我也尝试用PySCF来做，计算速度仍然很慢。是不是我的基组选择有问题呀，cc-pVDZ不适合这样的长链有机化合物？
  谢谢各位专家！

GoldenBaby

直接上orca解决一切烦恼。
对于你这在情况，我更倾向于是gaussian没有判断出来对称性（如果有的话）。没有判断出来对称性的gaussian计算确实非常慢，这很正常。

Boer

基组的选择可以参考社长的帖子（http://sobereva.com/336）

这么大的柔性应该考虑做一下分子动力学或者构象搜索吧，不然就要简化一下结构。(另外，P507好像是以二聚体的形式配位)

dnlx

ORCA建议用PBEh-3c。
过渡金属化合物的优化我试过B97-3c，键长和频率不准确。r2SCAN-3c也可能不准确。PBEh-3c和PBE0-D4/TZ结果比较接近。这些组合方法的基组都考虑相对论效应了。

sobereva

先说清楚用的什么理论方法。对于DFT本来就不该用cc-pVDZ，下文明确说了
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）
cc-pVDZ换成6-31G*或def2-SV(P)。

“带有赝势的cc-pVDZ级别基组”是不当描述。赝势基组那叫cc-pVDZ-PP。换成SDD或lanl2TZ赝势基组，看下文
谈谈赝势基组的选用
http://sobereva.com/373（http://bbs.keinsci.com/thread-5625-1-1.html）

不说计算资源的前提下说快慢没意义。像样的双路服务器上对这个体系用Gaussian做opt freq完全跑得动，绝对算不上“非常慢”。

ORCA虽然凭借RIJ/RIJCOSX使几何优化每一步比Gaussian更快，但几何优化和SCF的收敛性弱于Gaussian，而且振动分析速度慢，故从实际性价比来讲未必比用Gaussian在PBE0-D3(BJ)/6-31G*&SDD级别下搞这玩意有优势。所以我不赞同说成 “上orca解决一切烦恼”