|
本帖最后由 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的位置。目前我已想不到其他的解决办法,请教各位老师,有没有其他的建议。谢谢各位老师!
|
-
1.png
(22.78 KB, 下载次数 Times of downloads: 4)
步长为1
-
2.png
(103.51 KB, 下载次数 Times of downloads: 3)
步长为6
-
3.png
(97.06 KB, 下载次数 Times of downloads: 4)
步长为6的RMS Gradient Norm图
|