计算化学公社

标题: gaussian结构优化和频率计算报错请教 [打印本页]
作者
Author: greedbob    时间: 2020-5-15 14:24
标题: gaussian结构优化和频率计算报错请教
本帖最后由 greedbob 于 2020-5-15 14:37 编辑

各位老师好，我在做对苯二甲酸阴离子的结构优化和振动分析，结构优化可以成功收敛，但继续频率分析时报错OrtVc1 failed #3.

有两个问题想请教一下，非常感谢各位，输入和输出文件已压缩放在附件 (, 下载次数 Times of downloads: 26) 。

1. 我用b3lyp/6-311+g(d,p)优化此阴离子合理吗？
具体输入文件为：
%chk=bpdc-311+gdp-3.chk
%mem=60GB
%nprocshared=12
#p opt=gdiis freq b3lyp/6-311+g(d,p) nosymm

bpdc-311+gdp-3

-2 1
C                 -0.74926500    0.52086000    0.01950800
C                 -1.49160600    1.65058000   -0.34480700
C                 -1.47116100   -0.60897400    0.42256200
C                 -2.87784200    1.64703900   -0.30757000
H                 -0.96705300    2.53361900   -0.69179600
C                  0.72483800    0.52081300   -0.01954000
C                 -2.85740400   -0.60558500    0.45860000
H                 -0.92882500   -1.49187700    0.74137200
C                 -3.59288700    0.52063900    0.09458900
H                 -3.44840600    2.51977500   -0.59846100
C                  1.46720500    1.65103500    0.34316000
C                  1.44670600   -0.60961700   -0.42097500
H                 -3.41193900   -1.47825600    0.77924900
C                 -5.13464500    0.52014100    0.13518800
C                  2.85344200    1.64740700    0.30592900
H                  0.94267700    2.53458200    0.68889400
C                  2.83295000   -0.60630900   -0.45702600
H                  0.90435100   -1.49296800   -0.73851200
O                 -5.67266500   -0.55442500    0.50989400
O                 -5.69254300    1.59422900   -0.21063300
C                  3.56845800    0.52041800   -0.09462600
H                  3.42402900    2.52053900    0.59558100
H                  3.38746800   -1.47944600   -0.77643100
C                  5.11021700    0.51982800   -0.13523200
O                  5.64821000   -0.55528200   -0.50841900
O                  5.66814200    1.59438700    0.20908500

2. 一开始我没有使用opt=gdiis关键词，会在循环一段圈数后报错" Optimization aborted.     -- Gradient out of range.    -- Maximum allowed force =    0.100E+21"，因此用了此参数，opt可以收敛了，但freq会报错"OrtVc1 failed #3." 应该如何处理？

非常感谢各位！

作者
Author: sobereva    时间: 2020-5-15 17:22
关键词本身没什么问题。是bug。拿最后一步结构直接做freq再试
作者
Author: greedbob    时间: 2020-5-16 21:28
sobereva 发表于 2020-5-15 17:22
关键词本身没什么问题。是bug。拿最后一步结构直接做freq再试

非常感谢您的指导，我再次在此opt结果上，继续用g16进行了freq分析，仍然报错"OrtVc1 failed #3."，换用g09 E01版本也是。请问有哪些建议可以尝试吗？我遇到过好几次这个报错了，谷歌也没有找到类似问题的解决方法，有点绝望。

另外opt之后做freq是 只需要方法和基组相同 呢？还是其他所有关键词都要相同，如scf积分精度关键词这些。

我正在用orca尝试计算，正如您在博文中提到过的，Pople系列基组在orca中不划算，花了一个小时左右还没算完。
请问对于有机阴离子体系，我用RI-B3LYP-D3(BJ)/def2-TZVP(-f)基组怎么样，不合适的话您有什么建议吗？

非常感谢您！
作者
Author: sobereva    时间: 2020-5-17 03:36
greedbob 发表于 2020-5-16 21:28
非常感谢您的指导，我再次在此opt结果上，继续用g16进行了freq分析，仍然报错"OrtVc1 failed #3."，换用g ...

上传振动分析任务的带#P的输出文件

影响势能面的设置必须相同。SCF相关设置尽量相同，而且本来也没有理由不相同

不是Pople系列在ORCA里不划算，而是ORCA的杂化泛函下的振动分析本来就远比Gaussian要慢（即便开RI）。当前是阴离子体系，为了避免审稿人吐槽起码给羧基氧带上弥散
作者
Author: greedbob    时间: 2020-5-18 10:33
sobereva 发表于 2020-5-17 03:36
上传振动分析任务的带#P的输出文件

影响势能面的设置必须相同。SCF相关设置尽量相同，而且本来也没有 ...

您好，重复opt计算发现，单独opt有时也会报错，但有时可以算出来。
上传的这次是opt计算收敛，但freq报错。这是我的输入文件和输出文件：

(, 下载次数 Times of downloads: 11)
谢谢您的指导。
作者
Author: sobereva    时间: 2020-5-18 23:59
greedbob 发表于 2020-5-18 10:33
您好，重复opt计算发现，单独opt有时也会报错，但有时可以算出来。
上传的这次是opt计算收敛，但freq报 ...

用你最后结构我做了freq没有任何问题。nosymm和novaracc完全多余，毫无必要写

(, 下载次数 Times of downloads: 37)
(, 下载次数 Times of downloads: 17)

如果用我的输入文件在你那里也不行，8成是你的Gaussian程序版本有bug，或者和机子有兼容性问题
作者
Author: greedbob    时间: 2020-5-19 09:32
sobereva 发表于 2020-5-18 23:59
用你最后结构我做了freq没有任何问题。nosymm和novaracc完全多余，毫无必要写

非常感谢！

还报错，应该如您所说是我机器或软件的问题了。
作者
Author: liyangNK    时间: 2022-10-28 17:29
greedbob 发表于 2020-5-19 09:32
非常感谢！

还报错，应该如您所说是我机器或软件的问题了。

你好！我们现在仪器遇到了相同的问题，请问最后怎么解决的？谢谢
作者
Author: greedbob    时间: 2022-10-30 00:13
liyangNK 发表于 2022-10-28 17:29
你好！我们现在仪器遇到了相同的问题，请问最后怎么解决的？谢谢

没有找到原因，测试了多个linux发行版、虚拟机都有比较冷门的报错，怀疑是3900X+X570平台和我使用的软件版本兼容性的问题，后面我个人改用orca了，特别是orca5.0发布以后，使用体验很不错~
作者
Author: liyangNK    时间: 2022-11-2 09:04
greedbob 发表于 2022-10-30 00:13
没有找到原因，测试了多个linux发行版、虚拟机都有比较冷门的报错，怀疑是3900X+X570平台和我使用的软件 ...

谢谢回复。 我们最近也被这个搞得头大



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3