请教下结构优化能量一直下降但力的收敛不理想是否应该继续优化

lulu沙

老师您好，我优化了一个石墨烯单胞，使用您提供的脚本检查收敛结果如下图，请问老师这个结果是收敛的嘛？因为Energy of each ionic step出现了两个峰，Max Force of each ionic step又呈上升趋势，所以我拿不准，想请教一下老师。

向往美好生活

白皛白 发表于 2023-1-1 14:29
早晨起来发现已经成功收敛了，这个任务总共耗时一周左右（24核），但最大受力的显著收敛只发生在最后二三十 ...

楼主说的很对，最近再算一个初始结构不知怎么摆放的任务，体系较大（slab有500多个原子，吸附一个有110个原子的体系）。第一分子+slab在结构优化时，受力一直下不去，一直在0.9-0.08之间跳到，那时已经算了210多步离子步了，但与此同时能量一直在单调下降。于是听楼主的建议继续往下算，最先出现的转机是能量下降曲线开始变得平缓了，再往下算果然就收敛了，力的收敛迹象很隐秘，到后面的20-30步离子步才出现.

但是还有一个发现，就是说按照楼主说的，“能量是力的积分，代表总体受力”。所有这个优化过程要从摆放的初始结构优化到最终结构，貌似要走或者说变动一定量的距离，这个变动量没达到时，即使缩小步长（力震荡下不去时缩小步长收敛方法）也没有用，他依然要走完一定的量，过早的缩小步长反而使得走的速度变慢了。这就是第二个分子与slab的优化的过程总结。
第二分子与第三分子结构十分相似，故可比较。第二个分子是0.5步长一直优化到200步，力在震荡，能量在单调下降，但并未出现下降减速拐点，这时等不住就开始减小步长为0.05了，结果步长0.05下算了300步未收敛，能量还在下降；第三个分子是0.5步长优化到240步，力同样在震荡，但能量在单调下降时出现了减速拐点，于是减小步长为0.05，算2个离子步后收敛。我认为这是优化走的路径还不够长就提前减速（减小步长导致的），不懂对不对，请大家批评指正！

最后感谢楼主分享！