高斯计算铁催化的氢转移过程时SCF能量前后不一致

syx123456

大家好，近期我在用Gaussian 09 D.01版本计算一个铁催化的硅氢转移过程时，遇到了如下问题：
1. 已收敛至不动点的过渡态在计算IRC时，同一结构给出的能量、能量梯度(RMS Gradient)不同，且个人认为超出可接受的误差范围。
2. 对过渡态进行固定优化后，将固定优化所得结构进行opt=TS计算，同一结构前后能量、能量梯度不一致。（使用了固定优化结构38作为过渡态计算初猜）
3. 优化中间体结构时，梯度明显处于不合理范围，但程序报告已收敛到极小点。
使用的泛函/基组：ωB97X-D/TZVP(Fe),6-31G(d)(C,H,N,Si)。目前怀疑此问题与SCF自洽场初猜有关，但是不知道应该如何解决。
输入输出文件在附件中，请问大家有什么比较好的解决办法？

wzkchem5

检查一下收敛的SCF波函数是否符合预期（通过看布居分析、前线轨道等等），如果不符合预期，通过旋转轨道得到符合预期波函数的初猜，重新做SCF

sobereva

那叫驻点，不叫不动点

对优化过渡态最后得到的波函数做波函数稳定性测试，如果稳定，而且IRC第一步的能量高于优化过渡态最后一步的，说明IRC第一步（过渡态结构）没收敛到稳定波函数。用guess=read读优化过渡态最后的chk文件

这种体系并不适合ωB97X-D。ωB97X-D对于含过渡金属体系并没有什么优点

syx123456

sobereva 发表于 2022-5-9 20:26
那叫驻点，不叫不动点

对优化过渡态最后得到的波函数做波函数稳定性测试，如果稳定，而且IRC第一步的能 ...

谢谢老师！