几何优化设定maxcyc=1，力已经算出来了，为何还是9999错误

agent99

跑Truhlar组的QMMM软件里的一个test case总是出错，检测到问题来自于gaussian部分，发现gaussian输出文件最后是9999错误，然而再往前翻找不到任何错误信息。这是gaussian部分的输入文件（看起来比较奇怪，其实就是用gaussian计算QM/MM的力（所以maxcyc=1）然后返回到QMMM主程序。那些He原子是QM和MM区域之间的虚原子）：

%chk=Test.chk   
#p  opt=(maxcyc=1,nomicro,cartesian) oniom( hf/midix:amber=softfirst)=(scalecharge=555)
# geom=connectivity nosym      
#       6d                                                                     


Gradient Calculation


   0   1
C-CT-0.0000D+00    0 -0.333633961D+01    0.187514194D+01   -0.134209217D+01   H
H-CT-0.0000D+00    0 -0.372725112D+01    0.138549578D+01   -0.222571788D+01   H
H-CT-0.0000D+00    0 -0.381837794D+01    0.283132479D+01   -0.120396589D+01   H
O-CT-0.0000D+00    0 -0.195296177D+01    0.207934351D+01   -0.137465093D+01   H
H-CT-0.0000D+00    0 -0.155463187D+01    0.131228162D+01   -0.180714799D+01   H
H-HC-0.0000D+00    0 -0.361957080D+01    0.124728296D+01   -0.466939436D+00   H
He-HC--.3831D+00   0 -0.313506052D+01   -0.211709100D+00   -0.183846630D+00   L
He-HC--.3831D+00   0 -0.338286780D+01    0.157179455D+01    0.105156604D+01   L
He-HC--.3831D+00   0 -0.506476134D+01    0.793660580D+00   -0.904621300D-01   L
He-HC-0.3549D+00   0 -0.343431942D+01    0.391422605D+00   -0.149257805D+00   L
He-HC-0.3549D+00   0 -0.355822306D+01    0.128317443D+01    0.468448530D+00   L
He-HC-0.3549D+00   0 -0.439916983D+01    0.894107445D+00   -0.102565555D+00   L


     1
     2
     3
     4
     5
     6     1  1.00
     7
     8
     9
    10
    11
    12


NonBon  0 1

输出文件末尾（完整文件见附件）：
按说力已经算出来了，没有任何错误信息，为何会出现9999错误？

         Item               Value     Threshold  Converged?
Maximum Force            0.120613     0.000450     NO
RMS     Force            0.044326     0.000300     NO
Maximum Displacement     0.114488     0.001800     NO
RMS     Displacement     0.040064     0.001200     NO
Predicted change in Energy=-3.165346D-02
Optimization stopped.
    -- Number of steps exceeded,  NStep=   1
    -- Flag reset to prevent archiving.
Lowest energy point so far.  Saving SCF results.
Largest change from initial coordinates is atom    0       0.000 Angstoms.
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

...
Error termination request processed by link 9999.
Error termination via Lnk1e in /act/apps/Gaussian16a03/g16/l9999.exe at Mon Feb 10 14:58:40 2020.

zjxitcc

恐怕你对“没有任何错误信息”有所误解，输出文件很清楚地写着

Number of steps exceeded

因为你设定1步，而1步有没有达到收敛要求，就会报错说超过步数上限。

关键词force就是算力的，你可以试试这个关键词，无需opt。理论上讲这个更符合你的需求，但是实际操作一下方知可行性。

agent99

zjxitcc 发表于 2020-2-10 20:33
恐怕你对“没有任何错误信息”有所误解，输出文件很清楚地写着

Number of steps exceeded

谢谢回答。

这个输入文件是QMMM软件生成的，我没有修改。我对它的理解是用来算QM区域的力然后返回给QMMM主程序做整体的结构优化，不然跑1步没有意义（当然因为我也是刚开始学这个软件所以我的理解可能不对）。

我给软件作者发了信，目前还未得到答复。

sobereva

这种目的就应该用force关键词，我怀疑是程序开发者不知道Gaussian有这个关键词