计算一个二茂铁在碳环内转一圈的旋转能垒，过渡态收敛但无虚频

ThualishuT

想计算一个二茂铁在碳环内转一圈的旋转能垒，先在PM7下做了柔性扫描和过渡态、极小点结构的预优化以及附带的振动分析、IRC。柔性扫描显示二茂铁转动一圈过程只有一个过渡态，在PM7下做过渡态搜索，虚频方向、IRC均沿正确的路径。取预优化得到的过渡态结构为初猜结构使用关键词Opt(CalcFC,MaxStep=5,TS,NoEigen,Tight) Freq wB97xD/GenECP（C/H用6-31G**，Fe用SDD）进行过渡态搜索，优化过程结构变化不大，结果收敛的较好，但振动分析显示无虚频，可能是什么原因？恳请各位赐教。

ThualishuT

Deepast 发表于 2023-3-24 01:34
不好说，因为你是换了方法再进行计算的
你这个过渡态在PM7仅有13的虚频，虽然你说IRC可以跑通，但换成wB97 ...

感谢指点，我再尝试尝试。

ThualishuT

sobereva 发表于 2023-3-24 01:20
过渡金属配合物别拿半经验方法说事，就算跑出来你期望的结果，也极有可能DFT验证不出来，此时说明半经验方 ...

感谢指点，我再尝试尝试。

Deepast

不好说，因为你是换了方法再进行计算的
你这个过渡态在PM7仅有13的虚频，虽然你说IRC可以跑通，但换成wB97xD不知道过渡态是不是这个结构
根据我的经验我可以提一些建议：
1.过渡态优化前先给初始结构做一下频率分析，一定要有虚频，再看看最大虚频率是否对应反应方向。只有虚频正确，过渡态优化算法才能跟踪到正确的过渡态。频率正确后再readfc也不会增加多少耗时。
2.条件允许还是用wB97xD扫描一遍
3.尝试其它找过渡态方法

我个人也使用过半经验初猜来找过渡态，只能说对于某些体系还是有效果的。不过直接优化的话很多时候是靠运气，毕竟理论方法根本不同。

sobereva

过渡金属配合物别拿半经验方法说事，就算跑出来你期望的结果，也极有可能DFT验证不出来，此时说明半经验方法的结果定性错误。
你的这个体系适当简化，只保留二茂铁附近的部分并冻结边缘，DFT结合这样的简化的模型来计算也花不了多少时间
另外，wB97XD计算此体系没有什么好处，明显不如用PBE0-D3(BJ)或者TPSSh-D3(BJ)来算