计算化学公社

标题: 求助-金属钠切面离子步收敛困难问题 [打印本页]

作者
Author: Penson 时间: 2025-8-22 12:14
标题: 求助-金属钠切面离子步收敛困难问题
本帖最后由 Penson 于 2025-8-22 12:17 编辑

咨询下各位大佬，计算7*7四层的金属Na（110）切面的结构优化时，电子步收敛非常正常，但离子步力和位移的收敛困难，每个离子步的能量也是一直波动状态（如下图）。
(, 下载次数 Times of downloads: 7) (, 下载次数 Times of downloads: 6)
优化了500个离子步后结构看起来也比较正常，请教一下问题可能出在什么地方。以下为初始结构及优化500个离子步后的结构图，附件附上了计算输入文件和结构。
(, 下载次数 Times of downloads: 6) (, 下载次数 Times of downloads: 7)

作者
Author: 1138711019 时间: 2025-8-22 13:28
试试先不加D3呢

作者
Author: Penson 时间: 2025-8-22 13:36

1138711019 发表于 2025-8-22 13:28
试试先不加D3呢

D3也会影响这个吗，我试试看哈

作者
Author: sobereva 时间: 2025-8-23 12:32
Na结合GTH赝势基组做优化时应当把CUTOFF设得很高，比如800，否则可能受力计算精度明显不理想。也可以改用GAPW结合诸如pob-DZVP-rev2全电子基组，这样CUTOFF 400就够。
D3对这个的优化的收敛性不会有什么影响
其它可以考虑的是优化算法用CG，或者减小步长上限。也可以考虑原胞优化+足够多的k点，优化完了再扩胞成实际需要的slab，这样快得多

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