在opt固氮酶活性中心结构时，结构不收敛

大头攒毛

各位老师好，请问一下，我在opt FeMo-co时，总是提示scf不收敛，计算方法和基组也尝试着换过，但发现还是不收敛。请教各位老师，麻烦指导一下，究竟是哪儿的问题？

大头攒毛

wzkchem5 发表于 2023-6-27 14:19
这个体系有反铁磁耦合，必须准备片段初猜才能算，直接算即使能收敛，几乎肯定也会收敛到错误的波函数。建议 ...

谢谢老师答疑解惑

wzkchem5

这个体系有反铁磁耦合，必须准备片段初猜才能算，直接算即使能收敛，几乎肯定也会收敛到错误的波函数。建议查一下Serena DeBeer的文章，了解一下固氮酶活性中心怎么算。
还有一个问题，高斯应该算不了顺磁性体系的核磁吧？据我所知orca可以算但高斯可能不行。顺磁性体系有额外的顺磁屏蔽项的贡献，和温度、EPR g值、零场分裂有关。像orca是单独有一个模块算顺磁性体系的核磁。不过高斯是不是真的不能做这个我不太了解，还请其他人补充。
此外，Mo的配体简化成羟基和氨基恐怕是简化过度了，估计反映不了真实的配体环境。