计算化学公社

标题: 在服务器上的高斯16把优化好的结构进行频率计算时为什么总是出现错误QPErr [打印本页]

作者
Author: xuziwen 时间: 2025-9-16 09:20
标题: 在服务器上的高斯16把优化好的结构进行频率计算时为什么总是出现错误QPErr
求助各位老师，这个频率输入文件我已经改了很多次了，只要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再进行一下频率计算看看情况，但是却出现了图片中的报错 (, 下载次数 Times of downloads: 0)
这是我的输入文件
%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

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

作者
Author: xuziwen 时间: 2025-9-16 10:46

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老师的帖子有专门讲这种情况

作者
Author: paramecium86 时间: 2025-9-16 10:55

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

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

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

作者
Author: xuziwen 时间: 2025-9-17 08:19

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

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

作者
Author: xuziwen 时间: 2025-9-17 08:20

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

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

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