计算化学公社

标题: 高斯计算L301错误求助 [打印本页]

作者
Author: daji 时间: 2025-3-10 10:47
标题: 高斯计算L301错误求助
高斯计算L301错误，尝试末尾多打空行，也尝试复制之前的之前的基组，还是报错301，求助大佬如何解决

作者
Author: Solitude198 时间: 2025-3-10 10:59
不是，什么信息都不给，301报错有那么多，让我猜？ http://bbs.keinsci.com/thread-4829-1-1.html
(, 下载次数 Times of downloads: 19)

作者
Author: lemon_electron 时间: 2025-3-10 11:03
好歹输入输出文件放一下，最起码截图呗

作者
Author: daji 时间: 2025-3-10 11:15

Solitude198 发表于 2025-3-10 10:59
不是，什么信息都不给，301报错有那么多，让我猜？ http://bbs.keinsci.com/thread-4829-1-1.html

是 end of file resding basis center

作者
Author: daji 时间: 2025-3-10 11:16

lemon_electron 发表于 2025-3-10 11:03
好歹输入输出文件放一下，最起码截图呗

是 end of file reading basis center

作者
Author: daji 时间: 2025-3-10 11:27

Solitude198 发表于 2025-3-10 10:59
不是，什么信息都不给，301报错有那么多，让我猜？ http://bbs.keinsci.com/thread-4829-1-1.html

是这个错误

作者
Author: sai77 时间: 2025-3-10 11:40

daji 发表于 2025-3-10 11:27
是这个错误

http://bbs.keinsci.com/thread-4829-1-1.html
解决办法就是多打空行
如果试了还不行，那就是你的输入文件格式有问题，传上来看下

作者
Author: daji 时间: 2025-3-10 12:09

sai77 发表于 2025-3-10 11:40
http://bbs.keinsci.com/thread-4829-1-1.html
解决办法就是多打空行
如果试了还不行，那就是你的输入 ...

是这块吗

作者
Author: Solitude198 时间: 2025-3-10 13:33

daji 发表于 2025-3-10 12:09
是这块吗

如果你只是要opt+freq的话，用

#p opt freq b3lyp/genecp

就行了

另外似乎有金属，泛函应该换 pbe0-d3bj

在说一句，哥么信息能给全吗，真让人瞎猜？体系里有什么不说，自己要算什么不说，现在用了混合基组也不说用什么

作者
Author: Stardust0831 时间: 2025-3-10 13:41

daji 发表于 2025-3-10 12:09
是这块吗

geom=connectivity对当前计算完全没有意义，建议删去。IOp(1/8=x)和scf=qc这两个关键词有些需要注意的地方，详见《常见的多余的和被滥用的Gaussian关键词》

21 IOp(1/8=x)、IOp(6/7=3)
IOp(6/7=3)相当于pop=full、IOp(1/8=x)相当于在opt里定义maxstep，明明有对应的关键词，却还老看到有人在输入文件里直接写这俩IOp。有对应关键词的场合不要写IOp，否则自己容易记错，别人看这输入文件的时候还得查IOp表。
顺带强调一下，优化极小点时用maxstep定义步长上限的时候一定要带着notrust，即opt(maxstep=x,notrust)，否则步长上限会随着优化进行自发发生变化而导致设定没起到作用。然而洒家经常看到有人用maxstep却不写notrust。而对于优化过渡态，本身默认就是用了notrust的，因此opt=TS任务还写notrust那是完全多余的。

22 scf=xqc或scf=qc
在《解决SCF不收敛问题的方法》中介绍过。很多人看到不收敛上来就用这个，甚至将之当成平日计算默认使用的关键词，这样做是很糟糕的！用qc（二次收敛）做迭代比一般方式迭代每步多花好几倍时间，虽然收敛的可能性高了不少，但是远不如先尝试其它解决不收敛的办法。而且，如果迫不得已只能用qc方式才能收敛，也极有可能最后得到的波函数并不稳定，此时也暗示很可能当前体系结构、方法、自旋多重度、基组等是不合适的。

scrf里写solvent=water也是无意义的，scrf默认就是水做溶剂。
需要确认当前体系带三个正电荷且自旋多重度为2。
用了genecp关键词，所以结尾需要定义基组和赝势，给出这个部分的截图。

作者
Author: daji 时间: 2025-3-10 14:02

Stardust0831 发表于 2025-3-10 13:41
geom=connectivity对当前计算完全没有意义，建议删去。IOp(1/8=x)和scf=qc这两个关键词有些需要注意的地 ...

你好，是这样的

作者
Author: daji 时间: 2025-3-10 14:02

Solitude198 发表于 2025-3-10 13:33
如果你只是要opt+freq的话，用

#p opt freq b3lyp/genecp

抱歉抱歉，我刚开始学，还有很多不太懂

作者
Author: Stardust0831 时间: 2025-3-10 14:30

daji 发表于 2025-3-10 14:02
你好，是这样的

基组极不合理，LANL2DZ太差了，使用混合基组时，赝势基组要和全电子基组精度更好或者相近才行。图省事的话，可以用def2-SVP基组。泛函不合理，当前有过渡金属Fe，B3LYP已经不再合适。具体可以看《简谈量子化学计算中DFT泛函的选择》

过渡金属-配位键键长优化：
3d金属：适合纯泛函或低HF成份泛函。首选TPSS或TPSSh。r2SCAN也不错
4d、5d金属：适合中等HF成份的泛函。PBE0最优

文件要改的地方比较多，弄不清楚的话可以给我个相似体系的结构文件，我帮你写个输入文件的example。

作者
Author: daji 时间: 2025-3-10 14:54

Stardust0831 发表于 2025-3-10 14:30
基组极不合理，LANL2DZ太差了，使用混合基组时，赝势基组要和全电子基组精度更好或者相近才行。图省事的 ...

好，我怎么发给您呀？

作者
Author: Stardust0831 时间: 2025-3-10 15:07

daji 发表于 2025-3-10 14:54
好，我怎么发给您呀？

发在论坛上。体系敏感的话可以设访问权限或者换成相似的体系。

作者
Author: daji 时间: 2025-3-10 15:16

Stardust0831 发表于 2025-3-10 15:07
发在论坛上。体系敏感的话可以设访问权限或者换成相似的体系。

您好，是这个，麻烦您啦，我之前算出一个基组一样的，但是我改了铁离子价态后就出现301end of file reading basis center错误。
我还想问在这个结构的基础上加入砷酸根离子，用6-311G（d,p）可以吗？非常感谢。

作者
Author: Stardust0831 时间: 2025-3-10 15:33

daji 发表于 2025-3-10 15:16
您好，是这个，麻烦您啦，我之前算出一个基组一样的，但是我改了铁离子价态后就出现301end of file readi ...

当前报错的直接原因是用了genecp却没有在文件末尾定义基组。
把基组替换成了def2-SVP，对当前的结构优化目的，通常足够了。原先设定的基组不太合理，LANL2DZ太差了，会很影响精度。
给你加了DFT-D3(BJ)矫正，这么做的原因见：
乱谈DFT-D
http://sobereva.com/83
DFT-D色散校正的使用
http://sobereva.com/210
《谈谈“计算时是否需要加DFT-D3色散校正？”》
http://sobereva.com/413

如果高斯版本低于g09的D01版本，就无法使用em=gd3bj关键词，此时应当换更新的高斯版本。
B3LYP在过渡金属体系的表现很平庸，建议再斟酌一下泛函的选用。能力允许的话筛选一遍泛函。

作者
Author: daji 时间: 2025-3-10 15:37

Stardust0831 发表于 2025-3-10 15:33
当前报错的直接原因是用了genecp却没有在文件末尾定义基组。
把基组替换成了def2-SVP，对当前的结构优化 ...

好的，非常感谢

作者
Author: sobereva 时间: 2025-3-10 22:32
以后发帖不要用“高斯计算新手求助”这种标题，必须尽可能具体，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。这次给你改了，以后注意

作者
Author: daji 时间: 2025-3-11 09:36

sobereva 发表于 2025-3-10 22:32
以后发帖不要用“高斯计算新手求助”这种标题，必须尽可能具体，仔细看http://bbs.keinsci.com/thread-9348 ...

好的好的

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