计算化学公社

标题: 求助gaussian优化，因力小于收敛限制100倍造成opt和freq都有no [打印本页]

作者
Author: shashengfo 时间: 2025-5-23 20:22
标题: 求助gaussian优化，因力小于收敛限制100倍造成opt和freq都有no
本帖最后由 shashengfo 于 2025-5-23 20:32 编辑

求助各位老师，gaussian16c02优化此分子由于力小于收敛限制100倍就直接收敛了，虽然没有虚频，但造成opt和freq都有no。Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法 - 思想家公社的门口：量子化学·分子模拟·二次元就是sob老师这篇博文里3.1说的那种情况，我不知道后面应该怎么处理了，有看到Gaussian做结构优化时遇到受力小于收敛标准100倍而收敛但位移还未达到收敛标准的问题 - 量子化学 (Quantum Chemistry) - 计算化学公社这篇帖子，但是不知道我这个分子能不能按柔性大分子那样不考虑位移收敛。希望能得到各位老师的指点，谢谢！

首先用#P PBE1PBE/def2SVP EmpiricalDispersion=GD3BJ freq opt=recalc=8优化后得到结果是opt四个yes，freq两个no。
Maximum Force          0.000034    0.000450    YES
RMS    Force          0.000004    0.000300    YES
Maximum Displacement    0.000490    0.001800    YES
RMS    Displacement    0.000099    0.001200    YES

Maximum Force          0.000034    0.000450    YES
RMS    Force          0.000004    0.000300    YES
Maximum Displacement    0.017973    0.001800    NO
RMS    Displacement    0.004174    0.001200    NO

再此结构上使用#p opt=(GDIIS,Recalc=5) freq def2svp empiricaldispersion=gd3bj int=superfine pbe1pbe 继续计算得到的结果是opt和freq显示最大位移是no且超出收敛限两倍，数值一样。（也试了calcall结果差不多）

Maximum Force          0.000001    0.000450    YES
RMS    Force          0.000000    0.000300    YES
Maximum Displacement    0.005438    0.001800    NO
RMS    Displacement    0.000542    0.001200    YES

作者
Author: chands 时间: 2025-5-23 21:24
你的分子柔性很大，势能面平缓，做结构优化未必能做到四个YES。现在做的已经很好了，注意点柔性分子的势能面可能有多个极小值，可能有多个稳定的结构。

作者
Author: shashengfo 时间: 2025-5-23 21:54

chands 发表于 2025-5-23 21:24
你的分子柔性很大，势能面平缓，做结构优化未必能做到四个YES。现在做的已经很好了，注意点柔性分子的势能 ...

感谢老师指教！

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