先说清楚用的什么理论方法。对于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解决一切烦恼”

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

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

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

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