计算化学公社

标题: 关于过渡态搜索中收敛判据的问题 [打印本页]

作者
Author: LYDong 时间: 2023-10-17 10:38
标题: 关于过渡态搜索中收敛判据的问题
本帖最后由 LYDong 于 2023-10-17 10:40 编辑

各位大佬好，新手纯菜鸡求助。我在做一个化学组成为C7H24O3Si4的过渡态搜索时，反应主体部分为乙烯吸收Si上的H原子并且与Si原子成键，理论方法和基组时b3lyp/6-31g(d,p)，输出文件末尾显示Normal termination，做震动分析也显示只有一个虚频，但是四个收敛判据其中的最大位移却显示没有收敛，请问对结果有影响吗？

作者
Author: RAL 时间: 2023-10-17 11:45
http://sobereva.com/278

作者
Author: wzkchem5 时间: 2023-10-17 16:33
本来结构优化收敛的判据就不是四个YES
如果比如说有一个NO，但离收敛很近，且其他三个YES收敛得特别好，那么程序也认为结构优化收敛了

作者
Author: sobereva 时间: 2023-10-18 03:56
2L提到的博文里说了

值得一提的是优化过程中即便位移没收敛，但只要受力小于收敛限的100倍也自动被Gaussian算作收敛。我认为这主要是考虑到势能面非常平缓的柔性大分子或分子团簇，相对于其体系的尺度，把结构收敛到那么精确意义也不大。这种情况下，连opt最后都有NO，显然freq最后也肯定有NO。若对结构精度要求高的话可以把优化的收敛限设严。

精度要求不很高的话可以直接用。如果想让最大位移能小点，可以结合calcall继续优化过渡态，反正体系也不大花不了多少时间

作者
Author: LYDong 时间: 2023-10-19 09:43

sobereva 发表于 2023-10-18 03:56
2L提到的博文里说了

精度要求不很高的话可以直接用。如果想让最大位移能小点，可以结合calcall继续优化 ...

谢谢老师解答

作者
Author: LYDong 时间: 2023-10-19 09:43

wzkchem5 发表于 2023-10-17 16:33
本来结构优化收敛的判据就不是四个YES
如果比如说有一个NO，但离收敛很近，且其他三个YES收敛得特别好，那 ...

谢谢解答

作者
Author: LYDong 时间: 2023-10-19 10:03

RAL 发表于 2023-10-17 11:45
http://sobereva.com/278

谢谢解答

欢迎光临计算化学公社 (http://bbs.keinsci.com/)