计算化学公社

标题: la离子9H2O合物优化不收敛，求助 [打印本页]

作者
Author: KYJY 时间: 2023-8-29 22:07
标题: la离子9H2O合物优化不收敛，求助
卢老师，您好，各位老师好，请教老师我在优化La镧离子的九水合物的时候，体系的自旋多重度为2，使用的泛函为PBE0，非金属原子O和H为6-31G*，La离子为SDD赝势基组，优化完的输出文件显示，H原子跑远了，我第一次优化自选多重度为2的化合物，心里没底，我修改关键词# opt=calcall genecp pbe1pbe，取最后一帧结构，再去跑，下载优化的pdb文件，监测发现H原子还是跑远了，我要使用其他泛函进行计算吗？使用其他泛函的话，有推荐吗？老师，输入文件和输出文件上传如下，请各位老师指教，谢谢啦，感激不尽

作者
Author: sobereva 时间: 2023-8-30 00:15
带上溶剂模型
说明稳定结构下就是部分以氢氧根方式配位的

作者
Author: Satoru 时间: 2023-8-30 01:49
La离子一般应该是+3价吧，这个电荷设定是不是有些问题，导致La形式上变成了0价，还原了O-H

作者
Author: KYJY 时间: 2023-8-30 10:44

sobereva 发表于 2023-8-30 00:15
带上溶剂模型
说明稳定结构下就是部分以氢氧根方式配位的

好的，谢谢老师，我加了溶剂模型，关键词# PBE1PBE/genecp scrf=SMD opt freq，请问老师，关键词是这样写的吗？

作者
Author: KYJY 时间: 2023-8-30 10:47

Satoru 发表于 2023-8-30 01:49
La离子一般应该是+3价吧，这个电荷设定是不是有些问题，导致La形式上变成了0价，还原了O-H

这个gaussView自己判定的，双重态是没问题的，体系的电荷我也不太确定，请问老师，怎么确定体系的电荷？

作者
Author: 北大-陶豫 时间: 2023-8-30 11:02
你想算的是3价La对吧？
La离子带3个正电荷，再加一些中性的水分子，总的电荷就是+3啊。

作者
Author: KYJY 时间: 2023-8-30 15:10

北大-陶豫发表于 2023-8-30 11:02
你想算的是3价La对吧？
La离子带3个正电荷，再加一些中性的水分子，总的电荷就是+3啊。

好的，谢谢老师，那我要删除游离的H原子，从刚开始去优化吗？我加了溶剂模型，只有RMS force收敛了额，其余没收敛，这个数据可以用吗？

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