计算化学公社

标题: 几何优化设定maxcyc=1，力已经算出来了，为何还是9999错误 [打印本页]

作者
Author: agent99 时间: 2020-2-11 07:28
标题: 几何优化设定maxcyc=1，力已经算出来了，为何还是9999错误
跑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.

作者
Author: zjxitcc 时间: 2020-2-11 10:33
恐怕你对“没有任何错误信息”有所误解，输出文件很清楚地写着

Number of steps exceeded

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

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

作者
Author: agent99 时间: 2020-2-12 09:49

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

Number of steps exceeded

谢谢回答。

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

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

作者
Author: sobereva 时间: 2020-2-13 17:55
这种目的就应该用force关键词，我怀疑是程序开发者不知道Gaussian有这个关键词

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