计算化学公社

标题: 求助gaussian计算丙酸根与钴的结构无法收敛 [打印本页]

作者
Author: 徐徐徐 时间: 2024-6-28 08:52
标题: 求助gaussian计算丙酸根与钴的结构无法收敛
我想计算其结合能但是，一直无法收敛具体输入与输出如图，麻烦哪位老师可以指正一下错误

作者
Author: Dyyyzz 时间: 2024-6-28 09:09
解决SCF不收敛问题的方法
http://sobereva.com/61

不要用m06-2x优化过渡金属

作者
Author: wxyhgk 时间: 2024-6-28 09:48
本帖最后由 wxyhgk 于 2024-6-28 09:49 编辑

1. 首先打开 out 文件

(, 下载次数 Times of downloads: 13)

从提示中可以知道是 L502 模块报错，接着查询《高斯16说明书》第八页，这里我整理了，发现是 SCF 不收敛的报错

(, 下载次数 Times of downloads: 14)

2. 检查 gjf 输入文件

(, 下载次数 Times of downloads: 16)

含有过渡金属 Co，打开基组网站，https://www.basissetexchange.org/

(, 下载次数 Times of downloads: 11)

可以看到你的 6-311++G(d,p) 是没有包含 Co 在里面的，所以基组的选择不是很合理，可以考虑用 SDD 或者 Lanl2DZ 等赝势基组。

作者
Author: sobereva 时间: 2024-6-28 09:56
Gaussian做计算对Co用6-311G系列没有逻辑上的问题，虽然有更好选择，看下文（lanl2DZ就太过时了，不建议用）。Gaussian里的6-311G支持的元素远比BSE上的多，因为Gaussian为了让6-311G系列能支持更多元素，对一些6-311G基组原文没有定义的元素在用6-311G关键词时会使用其它人定义的基组。具体看Gaussian手册里基组部分。

谈谈赝势基组的选用
http://sobereva.com/373（http://bbs.keinsci.com/thread-5625-1-1.html）

作者
Author: 徐徐徐 时间: 2024-6-28 14:40

Dyyyzz 发表于 2024-6-28 09:09
解决SCF不收敛问题的方法
http://sobereva.com/61

好的，谢谢老师，麻烦问一下为什么M062X不可以用于过渡金属吗？

作者
Author: 徐徐徐 时间: 2024-6-28 14:42

wxyhgk 发表于 2024-6-28 09:48
1. 首先打开 out 文件
从提示中可以知道是 L502 模块报错，接着查询《高斯16说明书》第八页，这里我整理了 ...

谢谢老师，我之前也尝试过Lanl2TZ基组进行计算，可出现的依旧是l502。我在尝试一下

作者
Author: Dyyyzz 时间: 2024-6-28 15:12

徐徐徐发表于 2024-6-28 14:40
好的，谢谢老师，麻烦问一下为什么M062X不可以用于过渡金属吗？

因为m06-2x是根据主族元素设计的参数，见
简谈量子化学计算中DFT泛函的选择 http://bbs.keinsci.com/thread-536-1-1.html

作者
Author: wxyhgk 时间: 2024-6-28 23:45

徐徐徐发表于 2024-6-28 14:40
好的，谢谢老师，麻烦问一下为什么M062X不可以用于过渡金属吗？

其实首先要考虑是结构问题，不知道你这个结构是对的吗

作者
Author: xsc6 时间: 2024-6-29 01:38
你好，看起来你应该是刚接触量化计算，实话实说就是你这个输入文件写的不好，很多不合理的地方很多，没有加#p也不知道收敛到什么程度差多少，而且看起来结构也不太好，Co在两个羧基中间卡住应该更好，综上，我的建议是赶紧买http://www.keinsci.com/workshop/KEQC_content.html讲义学习，时间很宝贵，也建议不要购买其他任何机构例如“某算科技”昂贵的视频课，看完都是副作用

作者
Author: 徐徐徐 时间: 2024-6-29 18:29

wxyhgk 发表于 2024-6-28 23:45
其实首先要考虑是结构问题，不知道你这个结构是对的吗

老师，我对C H O做了6-311g基组，Co使用了SDD基组，同样出现了L502错误，这个如果结构错位，请问我该如何更改呢？是调换一下两个丙酸根位置吗？

作者
Author: 徐徐徐 时间: 2024-6-29 18:32

xsc6 发表于 2024-6-29 01:38
你好，看起来你应该是刚接触量化计算，实话实说就是你这个输入文件写的不好，很多不合理的地方很多，没有加 ...

好的老师，我是在默认128步内没有收敛，也尝试SOB老师的收敛办法也失败，请问卡主是指Co放在两个丙酸根的脱氢羧基中间吗？

作者
Author: 徐徐徐 时间: 2024-6-29 18:36

sobereva 发表于 2024-6-28 09:56
Gaussian做计算对Co用6-311G系列没有逻辑上的问题，虽然有更好选择，看下文（lanl2DZ就太过时了，不建议用 ...

>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -980.203544444    A.U. after  129 cycles
         NFock=128  Conv=0.74D-03    -V/T= 2.0642
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 1.0716 S= 0.6496
<L.S>= 0.000000000000E+00
KE= 9.211067476491D+02 PE=-4.536746934631D+03 EE= 1.521490143039D+03
Annihilation of the first spin contaminant:
S**2 before annihilation    1.0716, after    0.7513
Convergence failure -- run terminated.
Error termination via Lnk1e in K:\G16W\l502.exe at Sat Jun 29 17:19:11 2024.
Job cpu time:    0 days  2 hours 16 minutes 10.0 seconds.
Elapsed time:    0 days  2 hours 16 minutes  7.3 seconds.
File lengths (MBytes):  RWF= 170 Int=    0 D2E=    0 Chk=    79 Scr=    1
老师我加了SDD赝势基组之后依旧L502,同时我还尝试了您的解决收敛的方法，可是依旧没有在128内收敛，是否我可以认为是结构问题

作者
Author: wxyhgk 时间: 2024-6-29 20:36

徐徐徐发表于 2024-6-29 18:32
好的老师，我是在默认128步内没有收敛，也尝试SOB老师的收敛办法也失败，请问卡主是指Co放在两个丙酸根的 ...

请告诉我具体的什么情况，你这个我没法判断。

作者
Author: xsc6 时间: 2024-6-29 22:01

徐徐徐发表于 2024-6-29 18:32
好的老师，我是在默认128步内没有收敛，也尝试SOB老师的收敛办法也失败，请问卡主是指Co放在两个丙酸根的 ...

你试试形成一个四配位你先用2-zeta基组试一下记得带上极化，别6-311g 加#p 看看收敛还差多少

作者
Author: sobereva 时间: 2024-6-30 05:27

徐徐徐发表于 2024-6-29 18:36
>>>>>>>>>> Convergence criterion not met.
SCF Done: E(UB3LYP) = -980.203544444 A.U. after ...

没具体输入和输出文件没法准确判断
说不定gjf文件还此存在什么硬伤，别人看不到文件无从得知

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