计算化学公社

标题: 一个含碘体系尝试部分方法后仍难以收敛 [打印本页]

作者
Author: 不干正事 时间: 2022-1-5 08:50
标题: 一个含碘体系尝试部分方法后仍难以收敛
大家好，对一个含有I的体系(C,H,N,Al,I,中性CS)用linux版g16做优化和计算频率时，一开始用BP86，Al和I用def2tzvp，CHN用def2SVP，第一轮就不能收敛，随后陆续尝试了博文中以下方法：由于没有用弥散函数，SCF=NoVarAcc没有使用，SCF=noincfock，SCF=vshift=400，guess=huckel也都没有生效，也是第一轮收敛128次报错。因博文中说明尼苏达泛函难收敛，所以就尝试用b3lyp泛函的def2svp时，能完全正常走完，但和得到的晶体结构有键角的误差较大，随后以它作为初猜，泛函改用b3lyp，单独对I使用def2SVp，其余用6-311g*时，迭代好几轮后，再次报了未收敛的错误。计算时也加了#P，但分析不出来什么原因，也上传了输出文件。希望大家帮助，谢谢大家。

作者
Author: sobereva 时间: 2022-1-5 09:03
这种体系没任何理由用BP86。用B3LYP就完了，I用def2-TZVP或lanl08(d)，其它6-311G*。没必要尝试任何其它的计算级别。
def2-svp偏low，对这么小体系没考虑的意义。

作者
Author: 不干正事 时间: 2022-1-5 09:08

sobereva 发表于 2022-1-5 09:03
这种体系没任何理由用BP86。用B3LYP就完了，I用def2-TZVP或lanl08(d)，其它6-311G*。没必要尝试任何其它的 ...

好的，谢谢卢老师，我试试，还有就是除了将用b3lyp，I改为def2tzvp外，还要加其他关键词吗，比如还要加SCF=noincfock，用能级移动法等这些吗？

作者
Author: sobereva 时间: 2022-1-5 10:01

不干正事发表于 2022-1-5 09:08
好的，谢谢卢老师，我试试，还有就是除了将用b3lyp，I改为def2tzvp外，还要加其他关键词吗，比如还要加SC ...

不用

作者
Author: 不干正事 时间: 2022-1-5 10:41

sobereva 发表于 2022-1-5 10:01
不用

好的，卢老师，我试试看

作者
Author: 不干正事 时间: 2022-1-5 18:21

sobereva 发表于 2022-1-5 10:01
不用

卢老师，尝试b3lyp.CHNAl用6311g*，I用def2tzvp后，开始迭代几轮没问题，但后续又重新出现了不收敛的情况。

作者
Author: chemhou 时间: 2022-1-5 18:45

不干正事发表于 2022-1-5 18:21
卢老师，尝试b3lyp.CHNAl用6311g*，I用def2tzvp后，开始迭代几轮没问题，但后续又重新出现了不收敛的情况 ...

这明显是结构不合理引起的，I没有用ecp

作者
Author: 不干正事 时间: 2022-1-5 19:00

chemhou 发表于 2022-1-5 18:45
这明显是结构不合理引起的，I没有用ecp

初始的结构我是用测出来的晶体结构导入的，应该不会有太大问题。

作者
Author: snljty 时间: 2022-1-5 19:21
本帖最后由 snljty 于 2022-1-5 19:48 编辑

不干正事发表于 2022-1-5 19:00
初始的结构我是用测出来的晶体结构导入的，应该不会有太大问题。

你这明显没用对赝势啊。认真看卢老师博文《详解Gaussian中混合基组、自定义基组和赝势基组的输入》http://sobereva.com/60。
不想再回复一次了，【认真】看，卢老师说的是直接写在Route行也就是#开头的那行后面才不用genecp，后面自定义当然要写GenECP。

作者
Author: 不干正事 时间: 2022-1-5 19:45

snljty 发表于 2022-1-5 19:21
你这明显没用对赝势啊。认真看卢老师博文《详解Gaussian中混合基组、自定义基组和赝势基组的输入》http:/ ...

我看博文里的def2基组对第五周期后的元素没有必要特意加ecp，所以就没有加ecp来定义赝势。

作者
Author: chemhou 时间: 2022-1-5 20:10

不干正事发表于 2022-1-5 19:00
初始的结构我是用测出来的晶体结构导入的，应该不会有太大问题。

结构不合理引起的不收敛，而结构不合理是因为I没用ecp
很明白的两句话啊

作者
Author: 不干正事 时间: 2022-1-5 20:40

chemhou 发表于 2022-1-5 20:10
结构不合理引起的不收敛，而结构不合理是因为I没用ecp
很明白的两句话啊

不好意思，还是不太明白。我的理解是这个它晶体结构测出来的，应该是它真实存在的样子，应该是合理的，以及我看到博文上说用def2基组可以不用加ecp定义赝势，就没有加，不知道我的理解哪个地方没想通。谢谢解答，真还比较迷糊

作者
Author: wzkchem5 时间: 2022-1-5 20:50

不干正事发表于 2022-1-5 13:40
不好意思，还是不太明白。我的理解是这个它晶体结构测出来的，应该是它真实存在的样子，应该是合理的，以 ...

“卢老师说的是直接写在Route行也就是#开头的那行后面才不用genecp，后面自定义当然要写GenECP。”这句没理解吗？
“用def2基组可以不用加ecp定义赝势”是有适用条件的，具体来说就是只有当这个基组写在#开头的那行的时候才成立，并不是对于所有情况都成立。要搞清楚博文每一句话的适用条件

作者
Author: 不干正事 时间: 2022-1-5 21:06

wzkchem5 发表于 2022-1-5 20:50
“卢老师说的是直接写在Route行也就是#开头的那行后面才不用genecp，后面自定义当然要写GenECP。”这句没 ...

明白了，就是在最后一行定义的赝势还是不能少，谢谢了

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