在服务器上的高斯16把优化好的结构进行频率计算时为什么总是出现错误QPErr

xuziwen

求助各位老师，这个频率输入文件我已经改了很多次了，只要freq后面不加上=calcfc就可以正常运行，虽然输出结果无虚频，但是[color=rgba(0, 0, 0, 0.85)][backcolor=rgba(0, 0, 0, 0.04)]最大力、RMS 力、最大位移、RMS 位移这些不收敛
[color=rgba(0, 0, 0, 0.85)][backcolor=rgba(0, 0, 0, 0.04)]所以我想试试加上=calcfc再进行一下频率计算看看情况，但是却出现了图片中的报错
这是我的输入文件
%mem=10GB
%nprocshared=8
%chk=01_freq.chk
#p freq=calcfc b3lyp/genecp em=gd3


iop(1/8=50) scrf=(smd, solvent=methylcyclohexane)


temperature=407.15, pressure=29.61


01_freq


0 1
Ti                 0.59588700   -0.37130000    0.88763000
C                  0.76669000   -0.65484100   -1.46075900
C                  2.11056300   -0.66434600   -1.00221100
C                  2.39327000    0.62032400   -0.43980500
C                  1.22689600    1.42036500   -0.56874400
C                  0.18910800    0.63504400   -1.17984500
C                  3.71721900    1.06619700    0.11173300
C                  3.09238900   -1.79401400   -1.12938900
C                  0.06219400   -1.82039500   -2.10482100
C                  1.11159800    2.85017300   -0.10394800
Si                -1.59930200    0.82111400   -0.60980800
N                 -1.27671200    0.06009900    0.93414500
C                  0.78298300   -2.43229000    1.28592000
C                  1.51628400    0.46334500    2.59505300
H                  4.28603400    0.22299000    0.51261800
H                  3.59025400    1.78635100    0.92335800
H                  4.33116300    1.53650000   -0.66650400
H                  3.81980100   -1.78137300   -0.31302600
H                  3.65275800   -1.72542800   -2.06999700
H                  2.59413400   -2.76561800   -1.10332100
H                 -0.75091300   -2.20630700   -1.48202800
H                  0.75250800   -2.64660100   -2.28586400
H                 -0.36533900   -1.52580200   -3.06845200
H                  0.58837300    2.93645800    0.85428800
H                  0.56536200    3.44962400   -0.83652800
H                  2.10016600    3.29722600    0.02411000
H                  0.22301500   -2.62897600    2.21464700
H                  0.38691300   -3.09728700    0.50900600
H                  1.83047900   -2.71180900    1.46715000
H                  0.94398200    0.13952100    3.47825700
H                  1.55599200    1.55961700    2.60966600
H                  2.54430600    0.08960000    2.70728500
C                 -2.95001100   -0.23211900   -1.41183300
C                 -3.79772400   -0.89000700   -0.33576700
C                 -3.54311100   -0.99805400    0.97634400
C                 -2.34450100   -0.49918300    1.75351500
C                 -2.23990700    2.59065400   -0.51114100
H                 -2.55244000   -0.99673700   -2.08810300
H                 -3.59491600    0.40095900   -2.03706000
H                 -4.73603200   -1.32139300   -0.68718500
H                 -4.28832900   -1.49762800    1.59593300
H                 -1.95888200   -1.33214000    2.36127900
H                 -2.69639000    0.25034700    2.48364900
H                 -3.28960700    2.57158900   -0.19616400
H                 -2.19417800    3.08515300   -1.48854800
H                 -1.68416100    3.19833400    0.20613100


Ti     0
S    1   1.00
      4.3720000              1.0000000
S    1   1.00
      1.0980000              1.0000000
S    1   1.00
      0.4178000              1.0000000
S    1   1.00
      0.0872000              1.0000000
S    1   1.00
      0.0314000              1.0000000
P    1   1.00
     12.5200000              1.0000000
P    1   1.00
      1.4910000              1.0000000
P    1   1.00
      0.4859000              1.0000000
P    1   1.00
      0.0530000              1.0000000
P    1   1.00
      0.0160000              1.0000000
D    3   1.00
     20.2100000              0.0341682
      5.4950000              0.1710005
      1.6990000              0.4405849
D    1   1.00
      0.4840000              1.0000000
D    1   1.00
      0.1157000              1.0000000
****
C H Si N 0
6-31G(d,p)
****


TI     0
TI-ECP     2     10
d potential
  3
1    265.3263909            -10.0000000
2     47.7687815            -51.8427816
2     11.8903334             -9.1429145
s-d potential
  5
0     81.4730696              3.0000000
1     72.6496724             19.4825579
2     31.8128213            207.3349279
2      6.1664468            235.6744501
2      5.8268347           -166.8784387
p-d potential
  5
0     50.2966943              5.0000000
1     63.5089754              5.5348822
2     26.0996084            177.8419384
2      5.6022573            107.4207153
2      5.2171069            -71.9065902

paramecium86

你这是不是哪里看的教程看错了。opt=calcfc 是在优化的时候第一步精确地计算fc 。freq要算频率肯定是要算fc的。所以freq不会有这个选项。解决你这个问题通常是比如写 opt=calcall 算一次频率一般能解决不是4个yes的情况。

xuziwen

paramecium86 发表于 2025-9-16 10:04
你这是不是哪里看的教程看错了。opt=calcfc 是在优化的时候第一步精确地计算fc 。freq要算频率肯定是要算fc ...

你好，感谢回复，我的结果优化确实是已经写了opt=calcfc的，内容如下#p opt=(calcfc,cartesian,notrustupdate,gediis) b3lyp/genecp  em=gd3 iop(1/8=50)
尽管结构优化收敛都是四个YES但是进行频率计算是却是四个NO，所以也不知道该怎么办了，但是刚刚看到了一篇sob老师的帖子有专门讲这种情况

paramecium86

xuziwen 发表于 2025-9-16 10:46
你好，感谢回复，我的结果优化确实是已经写了opt=calcfc的，内容如下#p opt=(calcfc,cartesian,notrustup ...

这种可以直接写opt=calcall把当前已经优化过得结构重新算一下。大多可以解决。

zjxitcc

你这不会是通过AI修改的gjf文件吧。。。cartesian,notrustupdate,gediis和iop(1/8=50)这些往往都用不到。另外，粘贴输入文件时最好使用代码框，或者直接上传gjf文件作为附件。你现在直接贴在帖子内容里，空行数目都不对了。

xuziwen

paramecium86 发表于 2025-9-16 10:55
这种可以直接写opt=calcall把当前已经优化过得结构重新算一下。大多可以解决。

是的，按照您的方法解决了！感谢！

xuziwen

zjxitcc 发表于 2025-9-16 11:07
你这不会是通过AI修改的gjf文件吧。。。cartesian,notrustupdate,gediis和iop(1/8=50)这些往往都用不到。另 ...

是朋友教的。。。可能他也是Ai搜的吧，谢谢您的建议