计算化学公社

标题: 求助：跑带锕系元素IRC时，中途SCF不收敛 [打印本页]

作者
Author: Xx_DerEnTh 时间: 2025-2-17 13:14
标题: 求助：跑带锕系元素IRC时，中途SCF不收敛
本帖最后由 Xx_DerEnTh 于 2025-2-17 19:33 编辑

各位老师好！我跑带有锕系元素的IRC时（用guess=read geom=check读取的过渡态结构和chk文件）中途某个点出现了SCF不收敛，根据sob老师的解决SCF不收敛问题的方法 - 思想家公社的门口：量子化学·分子模拟·二次元博文，尝试了novaracc,noincfock,fermi（以往类似体系加这几个关键词解决过问题），没有解决问题。其中IRC步长第一次设置的8，IRC曲线可以走到能量为-703.9左右，后面一步就会出现SCF不收敛，于是根据sob老师在Gaussian中计算IRC的方法和常见问题 - 思想家公社的门口：量子化学·分子模拟·二次元的博文，尝试减小了步长，减小到了6，结果还是在-703.9左右出现SCF不收敛（如图2，图中大多数点都是20-40圈左右就收敛了，只有中间某一个点用了160多圈才收敛。图3为RMS Gradient Norm图）。最后一次将步长设置成了1，仍然是在中途某个点SCF不收敛（如图1,在这个点之前SCF基本上迭代20几圈就能收敛），奇怪的是这次减小步长后能量反而走不到步长为8或6的能量（-703.9左右），只走到了能量大约为-703.857的位置。目前我已想不到其他的解决办法，请教各位老师，有没有其他的建议。谢谢各位老师！

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