计算化学公社

标题: CP2K冻结框架后优化吸附小分子位置，优化收敛时能量Decrease in energy为NO [打印本页]

作者
Author: 沙耶NEO 时间: 2026-1-29 12:31
标题: CP2K冻结框架后优化吸附小分子位置，优化收敛时能量Decrease in energy为NO
本帖最后由沙耶NEO 于 2026-1-29 21:12 编辑

如图所示，在BFGS收敛时该步Decrease in energy为NO，请问是否需要特别在意这个情况，比如用最后一步结构重新进行一次优化。附上输入文件，删除了原子坐标部分。输出文件已压缩
(, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 0)

作者
Author: UW_0728. 时间: 2026-1-29 15:55
按理来讲是不太正常的，有可能收敛到的并不是真正意义上的极小点；如果晶胞结构看着符合预期的话应该勉强能用

作者
Author: 沙耶NEO 时间: 2026-1-29 16:31
本帖最后由沙耶NEO 于 2026-1-29 17:26 编辑

UW_0728. 发表于 2026-1-29 15:55
按理来讲是不太正常的，有可能收敛到的并不是真正意义上的极小点；如果晶胞结构看着符合预期的话应该勉强能 ...

这个不是最终需要的结构，所以直接在它基础上解冻所有原子优化最终构象应该是没问题吧

作者
Author: Stardust0831 时间: 2026-1-29 17:04
本帖最后由 Stardust0831 于 2026-1-29 17:06 编辑

没关系的，这个只是说明这步的能量变化不是负值，并不是收敛的判据。这一步的能量的变化本身也已经非常小了，直接无视这个no就行。no的原因你可以认为是数值噪音。
最大力（最关键的判据），最大步长、RMS步长、RMS力都已经足够收敛了。除非你要算声子谱（会对精度有明显更高要求，需要改小最大力收敛限），否则这个结构就是够用的。
不可以解冻基底的原子，那些是用来描述晶体内部环境的，如果解冻了就会产生一些不应该有的表面弛豫。
你没法直接传那么大的计算文件是因为要求计算化学论坛要求原始文件先压缩后再上传。

作者
Author: 沙耶NEO 时间: 2026-1-29 17:26

Stardust0831 发表于 2026-1-29 17:04
没关系的，这个只是说明这步的能量变化不是负值，并不是收敛的判据。这一步的能量的变化本身也已经非常小了 ...

感谢

作者
Author: sobereva 时间: 2026-1-29 17:43
仔细看置顶的社员必读贴对上传文件尺寸的要求

优化过程中途偶尔能量上升不用管，这是算法不理想导致或数值噪音的问题。只要最终能达到收敛限，并且最后一步的能量比优化过程中每一步都更低就行。

作者
Author: 沙耶NEO 时间: 2026-1-29 21:10

sobereva 发表于 2026-1-29 17:43
仔细看置顶的社员必读贴对上传文件尺寸的要求

优化过程中途偶尔能量上升不用管，这是算法不理想导致或数 ...

感谢卢老师

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