计算化学公社

标题: 求助：MedeA软件算氢嵌入β-MnO2的过渡态扩散曲线力不收敛怎么解决 [打印本页]

作者
Author: dzyang0426 时间: 2026-3-6 17:20
标题: 求助：MedeA软件算氢嵌入β-MnO2的过渡态扩散曲线力不收敛怎么解决
本帖最后由 dzyang0426 于 2026-3-6 17:25 编辑

各位老师好！我在算氢嵌入β二氧化锰的扩散曲线。初态和末态我使用的是相差一个周期的结构优化过的等价位点，我用的软件是MedeA，通过计算后的job.out看到第四步的力突然飙升，另外中间的image我设了9个，是系统线性插值生成的。
1.我导出来力飙升前的那一步的所有image图像，发现有的image图像会有部分锰锰键生成，不知道是软件显示原因还是真的八面体结构被破坏了，锰锰键键长大概在2.5左右；
2.是不是氢比较轻，容易和氧原子成键然后迁移过程再断键，整个过程更倾向于反应路径而不是扩散路径；
3.力飙升的问题是不是出现在氢氧距离不合理，比如氢距离氧的距离特别近导致受力过大；
4.MedeA软件进行过渡态搜索我找不到怎么设置POTIM这些功能或者不支持，包括这个软件弹簧弹力只能设置为正数，不能设置为-5。
应该怎么做才能使力收敛得到一条正确的扩散曲线，恳请各位老师帮帮忙指点，万分感谢！
(, 下载次数 Times of downloads: 17)
(, 下载次数 Times of downloads: 18)
(, 下载次数 Times of downloads: 17)
(, 下载次数 Times of downloads: 18)

作者
Author: 卡开发发 时间: 2026-4-18 12:51
我不用MedeA，可能有几个点可以尝试排查：
1、针对能量和力不正常的情况最好先检查下前几步不同的image是否存在SCF（电子步）迭代不收敛或者自洽场参数不合理的情况，以及NEB初猜路径和优化算法相关参数是否合理，例如步长太大结构崩溃导致后续自洽场不正常也有可能；
2、反应和扩散我看法是没有明确界限，但过渡态计算确实限定在近独立粒子等统计近似上；
3、4同1。

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