请教如何加快ccsd进行结构优化和频率计算

xdyjiay

各位老师，导师让我用CCSD计算B、N加一个卤素原子（F,Cl,Br）的体系，整个体系一共12个原子，我测试了CCSD/b3lyp*，用gaussian算的，用了12个核，大概两天完成了结构优化，频率计算到现在算了一两个星期，还没有出结果。请问各位老师，有什么方法可以尽快计算完？

wjc404

"CCSD/b3lyp*"是什么奇怪的方法？
化学式是怎样的？有多少个成单电子？

xdyjiay

化学式是BF(N5)2，带一个负电荷，平面结构，计算目的是探究这样的结构能否稳定存在，因为整个体系原子数不多，所以采用ccsd方法，计算条件比较差，所以没加弥散，用个最小可接受级别基组来进行计算。

wjc404

用CCSD的话，基组至少cc-pVTZ才恰当，此时你的体系三百多个基函数，双电子积分可能占100GB以上，物理内存要尽量给大。
高斯CCSD的频率是用梯度差分做的，对于你的体系需要算几十个解析梯度才能出频率。CFOUR程序支持解析频率，但计算速度相对于数值频率没有显著优势。
如果做频率只是要获得热力学量，其实可以有其他近似的方法，节省机时，事半功倍。比如用B3LYP-D3(BJ)/def-TZVP做优化和振动分析，然后用DLPNO-CCSD(T)/aug-cc-pVTZ算单点能校正。

sobereva

绝对不要用CCSD/b3lyp*这种有严重歧义而且非标准化的描述，完全是给别人带来麻烦

也绝对不要用什么“最小可接受级别基组来进行计算”，诸如写文章computational details的时候还这么写绝对要挨批。“最小可接受”是高度主观的

静态相关不很强的情况用双杂化做opt freq的结果通常就已经足够有说服力了（除非极小众的诸如18碳环。但对于18碳环也有可以正确描述的双杂化泛函）。如果真是静态相关强的情况，用CCSD做也没多大意义，还得用CASSCF/多参考

xdyjiay

感谢两位老师的解答和错误纠正，下次会更加注意避免使用非标准化的描述。这个体系静态相关并不强，根据您的回复，反复思考了一下，打算放弃ccsd，用def2TZVP基组，再用b3lyp,revDSDPBEP86，M062X等泛函测试一下我这个体系大概什么泛函精度可以。