计算化学公社

标题: Dmol3优化MgO共吸附CO2和C，结构优化不收敛 [打印本页]

作者
Author: 安好 时间: 2023-4-17 16:47
标题: Dmol3优化MgO共吸附CO2和C，结构优化不收敛
本帖最后由安好于 2023-4-17 16:58 编辑

大佬们求助，我在用dmol3算CO+CO=C+CO2的过渡态时找不到过渡态，就想用dmol3优化一下我的反应物和产物的结构，但是反应物结构优化收敛了，产物结构在优化的时候不收敛，请问是我的参数设置有问题吗？基底是MgO，附上我的.input和.outmol文件。已经算过MgO单独吸附C、MgO单独吸附CO2的吸附能，并筛选了最佳吸附位点，但是优化MgO共吸附C和CO2时总是失败，不收敛，请问是我的参数设置有什么问题吗？
求大佬指点啊！
(, 下载次数 Times of downloads: 15) (, 下载次数 Times of downloads: 15) (, 下载次数 Times of downloads: 18) (, 下载次数 Times of downloads: 21) (, 下载次数 Times of downloads: 20) 下面是.outmol文件的结果，没有收敛。
(, 下载次数 Times of downloads: 18)

作者
Author: 含光君 时间: 2023-4-17 19:13
参数基本没太大问题，可以自己查看outmol中结构优化是否有收敛趋势、并确认没有跑散再决定是否继续优化。

建模仅给top视角看不清楚，建议补侧视图。

作者
Author: 安好 时间: 2023-4-17 22:44

含光君发表于 2023-4-17 19:13
参数基本没太大问题，可以自己查看outmol中结构优化是否有收敛趋势、并确认没有跑散再决定是否继续优化。
...

您好，下面是补充的内容，结构并没有跑散，我认为可以继续优化，但是，是继续用原来的参数设置继续优化，还是可以改一下相关参数呢？
麻烦您看看有什么建议吗？
(, 下载次数 Times of downloads: 18)

下面是我优化前结构的侧视图
(, 下载次数 Times of downloads: 15)

下面是我优化后但未收敛的.xsd文件
(, 下载次数 Times of downloads: 18) (, 下载次数 Times of downloads: 15)

作者
Author: 含光君 时间: 2023-4-18 00:07

安好发表于 2023-4-17 22:44
您好，下面是补充的内容，结构并没有跑散，我认为可以继续优化，但是，是继续用原来的参数设置继续优化， ...

我觉得对于当前体系，k点精度或许太粗糙了，可改用221或者331
以及如果是我的话，把全电子基组改成DSPP赝势，可以节省计算量，基组改成4.4版本，并且勾选偶极矫正

优化时可固定住下面两层MgO原子，有助于结构优化收敛

作者
Author: aaq2800 时间: 2023-4-18 07:30
本帖最后由 aaq2800 于 2023-4-18 07:31 编辑

density mixing，可以都改成0.1，smearing可以再改小点0.005，displacement step改成0.01，K点设置有点少咯

作者
Author: 安好 时间: 2023-4-18 16:13

含光君发表于 2023-4-18 00:07
我觉得对于当前体系，k点精度或许太粗糙了，可改用221或者331
以及如果是我的话，把全电子基组改成DSPP ...

谢谢您，收敛了，开始用Dmol3搜索过渡态了

作者
Author: 安好 时间: 2023-4-18 16:13

aaq2800 发表于 2023-4-18 07:30
density mixing，可以都改成0.1，smearing可以再改小点0.005，displacement step改成0.01，K点设置有点少咯

谢谢

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