计算化学公社

标题: 求助：小团簇结构优化不收敛 [打印本页]

作者
Author: jimulation 时间: 2021-5-20 11:53
标题: 求助：小团簇结构优化不收敛
从MD模拟轨迹中，提取出一个由Li+和3个分子（离子）配位形成的团簇结构，在opt B3LYP-D3(BJ)/TZVP下优化结构始终不收敛，改为opt=loose B3LYP-D3(BJ)/SVP，依旧不收敛，到了132步的时候报l9999错误，Number of steps exceeded。输出文件显示Force不是很大，但是Displacement一直降不下来，并且有震荡。
请大家帮忙看看，有什么方法能帮助这个团簇几何优化收敛吗？感谢

作者
Author: jiangning198511 时间: 2021-5-20 13:46
你可要考虑先让部分结构优化好后，再添加一些结构

作者
Author: hdhxx123 时间: 2021-5-20 13:55
这样相对小的体系可以考虑calcall和recalc
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

作者
Author: wzkchem5 时间: 2021-5-20 14:54
梯度接近0但是步长降不下来的情况，可以减小步长试试

作者
Author: jimulation 时间: 2021-5-20 18:14

hdhxx123 发表于 2021-5-20 13:55
这样相对小的体系可以考虑calcall和recalc
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.c ...

谢谢，加了recalc=3还是没能收敛...

作者
Author: hdhxx123 时间: 2021-5-20 18:37

jimulation 发表于 2021-5-20 18:14
谢谢，加了recalc=3还是没能收敛...

还有一些方法在博文里，都试试吧，这种情况只能耐心一个个试了

作者
Author: jimulation 时间: 2021-5-20 19:19

hdhxx123 发表于 2021-5-20 18:37
还有一些方法在博文里，都试试吧，这种情况只能耐心一个个试了

是的，方法比较多，正在尝试中。

作者
Author: jimulation 时间: 2021-5-20 21:18

wzkchem5 发表于 2021-5-20 14:54
梯度接近0但是步长降不下来的情况，可以减小步长试试

谢谢

作者
Author: jimulation 时间: 2021-5-20 21:20

jiangning198511 发表于 2021-5-20 13:46
你可要考虑先让部分结构优化好后，再添加一些结构

谢谢，我看博文中没有这个方法，这个在操作时有什么技巧不？

作者
Author: sobereva 时间: 2021-5-20 22:00
优化时候甭用SMD。用默认的IEFPCM
不要用def-SVP，一律用def2-SVP
四个硫在一起的结构非常诡异，怀疑这个模型的化学意义

作者
Author: 喵星大佬 时间: 2021-5-21 01:59
S4是个什么鬼。。。。你怕不是把两个S8各砍了一半。。。。

作者
Author: wzkchem5 时间: 2021-5-21 03:56

喵星大佬发表于 2021-5-20 18:59
S4是个什么鬼。。。。你怕不是把两个S8各砍了一半。。。。

我看了一下楼主的输入输出文件，这个不是问题，因为总电荷是-3，说明S4是-2价离子。
比如在锂硫电池里，S8被还原成两个[S4]2-，很合理。

作者
Author: jimulation 时间: 2021-5-21 10:26

sobereva 发表于 2021-5-20 22:00
优化时候甭用SMD。用默认的IEFPCM
不要用def-SVP，一律用def2-SVP
四个硫在一起的结构非常诡异，怀疑这个 ...

好的，谢谢社长。S4是模拟的带2个负电的多硫离子

作者
Author: jimulation 时间: 2021-5-21 10:48

wzkchem5 发表于 2021-5-21 03:56
我看了一下楼主的输入输出文件，这个不是问题，因为总电荷是-3，说明S4是-2价离子。
比如在锂硫电池里， ...

是的

这个构型几何优化真是不容易收敛

作者
Author: wzkchem5 时间: 2021-5-21 15:43

jimulation 发表于 2021-5-21 03:48
是的这个构型几何优化真是不容易收敛

本来溶液里面每个[S4]2-至少得有2个Li+，你抠出一个Li+和两个[S4]2-，是不是也不太符合真实情况

作者
Author: jimulation 时间: 2021-5-21 16:23

wzkchem5 发表于 2021-5-21 15:43
本来溶液里面每个[S4]2-至少得有2个Li+，你抠出一个Li+和两个[S4]2-，是不是也不太符合真实情况

是从MD中提取的，可能确实难以稳定存在，估计也是优化不收敛的原因

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