Dmol3优化MgO共吸附CO2和C，结构优化不收敛

安好 · 发表于 Post on 2023-4-17 16:47:57

本帖最后由安好于 2023-4-17 16:58 编辑

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

下面是.outmol文件的结果，没有收敛。

含光君 · 发表于 Post on 2023-4-17 19:13:14

参数基本没太大问题，可以自己查看outmol中结构优化是否有收敛趋势、并确认没有跑散再决定是否继续优化。

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

安好 · 发表于 Post on 2023-4-17 22:44:24

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

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

下面是我优化前结构的侧视图

下面是我优化后但未收敛的.xsd文件

含光君 · 发表于 Post on 2023-4-18 00:07:33

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

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

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

aaq2800 · 发表于 Post on 2023-4-18 07:30:51

本帖最后由 aaq2800 于 2023-4-18 07:31 编辑

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

安好 · 发表于 Post on 2023-4-18 16:13:32

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

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

安好 · 发表于 Post on 2023-4-18 16:13:53

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

谢谢

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CASTEP/Dmol3/MS] Dmol3优化MgO共吸附CO2和C，结构优化不收敛