算β-H消除过渡态结构不收敛

lzhu

老师，您好，我在计算β-H消除过渡态一直没办法收敛，尝试了博文中的
尝试不同的优化方法（GDIIS，RFO，GEDIIS）使用更好的Hessian矩阵（opt=calcall）调整步长上限（opt=maxcyc=90）
改变收敛标准（opt=loose）
还进行了势能面扫描，没有发现明显的最高点。这些都没有成功，想问问老师有没有其他的方法帮助收敛或者使用其他的点表征一下这个过渡态的能量。

sai77

opt=maxcyc=90是调整的优化步数，调整步长上限是opt=maxstep=N。

我之前也碰到过一模一样的问题，花了一个多月才找到这个过渡态。
就是因为势能面非常非常平缓，柔性扫描从反应物扫产物或者从产物扫反应物都尝试一遍，步长要设低一点才能看到拐点。
找过渡态也结合calcall和maxstep=1，最后才找到。

sobereva

有时calcall（嫌昂贵用recalc）同时结合较小步长上限效果才好。分清楚maxcyc和maxstep
提到博文时必须明确给出文章链接。假定你说的是
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164
建议改用PBE0-D3(BJ)