含铷体系高斯优化结构出现多次l9999、l103报错

xkzb

利用高斯计算含铷体系的结构优化，菜鸟新上手出现了很多问题，也参考了sob老师的几何优化收敛的方法进行了多次调整，包括M062x改b3lyp、6-31G*改def2svp、gdiis改cartesian等等，附件是最新一次调整的输入和输出文件，输出仍报错l9999，依然出现了opt stopped。想问一下各位老师，这种情况下要怎么办呢？

cokie

可以看一下帖子：Gaussian 常见报错及解决方法；新手求助报错时的注意事项（http://bbs.keinsci.com/thread-4829-1-1.html）

事实上你的优化任务是“完成”的，只不过没能在规定步数内收敛（在前面提到的帖子里搜索L9999就可以找到）

建议你考虑以下修改：1. 改用赝势（http://sobereva.com/373），超过第五周期的元素优化的时候建议给相应元素添加赝势，节省优化资源；2. 不清楚为什么要加入两个配体，在你给的输入文件里，另一个配体和Rb根本就没有关系；3. 优化只会优化到你给的初始几何结构附近的势能面极小点，在初始结构并非极其合理的前提下，优化任务无法直接找到全局势能面最小点，建议你如果对两个配体和Rb是怎么配位的很清楚，初始摆位时就应该把配体和Rb摆在合适的位置。

zjxitcc

你这个Rb的摆放位置是认真的吗

xkzb

cokie 发表于 2024-5-11 13:33
可以看一下帖子：Gaussian 常见报错及解决方法；新手求助报错时的注意事项（http://bbs.keinsci.com/thread ...

好的好的，谢谢老师，我去仔细看一下

xkzb

zjxitcc 发表于 2024-5-11 16:48
你这个Rb的摆放位置是认真的吗

老师您好，我是先通过genmer组件随机生成了200种初始复合物构型，然后用xtb进行优化，对结果再用isostat进行分类，选了其中5个较低能量的构型进行的进一步优化，然后在这一步出现了各种问题来着

cokie

xkzb 发表于 2024-5-11 18:02
老师您好，我是先通过genmer组件随机生成了200种初始复合物构型，然后用xtb进行优化，对结果再用isostat ...

虽然看似你做的很严谨，但是既然Rb的配位的位置以及配位数你是清楚的，那你在用genmer生成初始结构的时候，起码要把Rb和与Rb配位的几个原子位置固定住。。。你这最后其中一个配体的三个氢元素和Rb“配位”，另一个配体自由玩耍，你稍微请教一下身边学过化学的人就知道，这个结构完全不可能在你进一步优化的选择范围内。

xkzb

cokie 发表于 2024-5-11 18:31
虽然看似你做的很严谨，但是既然Rb的配位的位置以及配位数你是清楚的，那你在用genmer生成初始结构的时候 ...

其实老师，我只是知道配位数来着，不太清楚配位的具体位置

wzkchem5

xkzb 发表于 2024-5-11 11:02
老师您好，我是先通过genmer组件随机生成了200种初始复合物构型，然后用xtb进行优化，对结果再用isostat ...

有没有看过用xtb算Rb的benchmark文章？Rb是比较少研究的元素，xtb的训练集里很可能很少有含Rb的分子，因此计算结果未必可靠。要结合化学常识，用肉眼确认搜出来的构型是否合理，不能优化出来是什么样就是什么样

xkzb

wzkchem5 发表于 2024-5-11 19:53
有没有看过用xtb算Rb的benchmark文章？Rb是比较少研究的元素，xtb的训练集里很可能很少有含Rb的分子，因 ...

谢谢老师回复，没怎么看过呢，我去搜索看看，菜鸟刚上手，不会的东西特别多，出错的地方肯定不少