计算化学公社

标题: Error termination request processed by link 9999. [打印本页]

作者
Author: 小和尚 时间: 2015-5-26 20:36
标题: Error termination request processed by link 9999.
我在做计算时，一直是link9999，把得到坐标优化，还是link 9999，请问要怎么解决？谢谢你啊
输入文件：
%chk=HAc2H2O.chk
%mem=1200MB
%nprocshared=4
#p m062x/6-311++g(d,p) opt=cartesian freq scf=tight

HAc2H2O-

0 1
C                -1.92379300 -1.28262400 -0.20283000
H                -1.78241300 -2.02526200 0.58759800
H                -2.94760400 -0.91490500 -0.13215000
H                -1.75099200 -1.75090000 -1.16701600
C                -0.97828000 -0.13260300 0.03901700
O                -1.14315300 0.71659600 0.87291000
O                0.10156200 -0.17163000 -0.75744000
H                0.63722800 0.63709900 -0.58205000
O                2.47573600 -1.59937300 0.24764000
H                2.04681800 -1.77699300 1.08533500
H                1.74987300 -1.44956600 -0.36346000
O                0.91302000 2.40279900 -0.16352000
H                0.47616100 2.23722500 0.68234300
H                0.20605300 2.74753400 -0.72443900

输出文件：
Sum of Mulliken charges with hydrogens summed into heavy atoms = 0.00000
Electronic spatial extent (au):  <R**2>=          980.2459
Charge=             0.0000 electrons
Dipole moment (field-independent basis, Debye):
X=          -2.5674 Y=             0.9954 Z=             0.1196  Tot=             2.7562
Quadrupole moment (field-independent basis, Debye-Ang):
XX=          -37.8173 YY=          -40.1711 ZZ=          -35.6015
XY=          10.2842 XZ=          -6.6516 YZ=             0.9341
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
XX=             0.0460 YY=          -2.3078 ZZ=             2.2618
XY=          10.2842 XZ=          -6.6516 YZ=             0.9341
Octapole moment (field-independent basis, Debye-Ang**2):
  XXX=          -80.6972  YYY=          -0.9853  ZZZ=          -2.8515  XYY=          -7.8483
  XXY=             1.2969  XXZ=          -0.8869  XZZ=             5.7317  YZZ=          -5.0852
  YYZ=             4.3959  XYZ=             7.7175
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX=          -688.9990 YYYY=          -376.5389 ZZZZ=          -56.9865 XXXY=          105.8595
XXXZ=          -58.7804 YYYX=          34.7233 YYYZ=          -3.5740 ZZZX=          -10.4710
ZZZY=             2.9648 XXYY=          -203.9045 XXZZ=          -112.0123 YYZZ=          -70.9250
XXYZ=          -5.5140 YYXZ=          -16.0232 ZZXY=          -5.9286
N-N= 2.516051771521D+02 E-N=-1.402490360665D+03  KE= 3.807275908183D+02
No NMR shielding tensors so no spin-rotation constants.
Leave Link  601 at Tue May 26 19:59:48 2015, MaxMem=  157286400 cpu:    0.4
(Enter /home/g09root/g09/l9999.exe)

Error termination request processed by link 9999.
Error termination via Lnk1e in /home/g09root/g09/l9999.exe at Tue May 26 19:59:48 2015.
Job cpu time:  0 days  1 hours 35 minutes  7.9 seconds.
File lengths (MBytes):  RWF=    31 Int=    0 D2E=    0 Chk=    3 Scr=    1

作者
Author: sobereva 时间: 2015-5-26 22:09
scf=tight是多余的。用笛卡尔坐标下优化通常没好处。

9成是几何优化不收敛，处理方法这里有详细讨论
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

具体试哪些关键词要看实际收敛趋势，结构怎么变的、能量怎么变的

作者
Author: 小和尚 时间: 2015-5-27 15:16

sobereva 发表于 2015-5-26 22:09
scf=tight是多余的。用笛卡尔坐标下优化通常没好处。

9成是几何优化不收敛，处理方法这里有详细讨论

非常感谢啊！

作者
Author: 王若曦 时间: 2023-3-25 16:24
您好，sob老师，我在做您给出的IRC实例计算分析。我在对SiH2+H2反应的过渡态搜索时先用柔性扫描初猜过渡态，也出现了上述这种的错误输出。我把输入、输出文件附上，请您看下。
输入文件

%nprocshared=4
%mem=600MW
%chk=SiH2+H2_relax.chk
#p opt=modredundant nosymm pm7

SiH2+H2_relax

0 1
Si
H                1          B1
H                1          B2 2          A1
H                1          B3 3          A2 2          D1 0
H                1          B4 3          A3 2          D2 0

B1          1.47000000
B2          2.83328061
B3          1.46999999
B4          2.98742337
A1          80.08090348
A2          87.30876442
A3          11.43867245
D1       -121.08892459
D2          -70.39968655

1 5 S  50 -0.04

输出文件

A COMPUTER IS LIKE A HORSE, IT WILL SENSE WEAKNESS.
   -- GREG WETTSTEIN
Error termination request processed by link 9999.
Error termination via Lnk1e in /home/yin/program/g16/g16/l9999.exe at Sat Mar 25 12:22:40 2023.
Job cpu time:    0 days  0 hours 12 minutes 18.8 seconds.
Elapsed time:    0 days  0 hours  4 minutes 21.6 seconds.
File lengths (MBytes):  RWF=    6 Int=    0 D2E=    0 Chk=    3 Scr=    1

作者
Author: wzkchem5 时间: 2023-3-25 16:29

王若曦发表于 2023-3-25 09:24
您好，sob老师，我在做您给出的IRC实例计算分析。我在对SiH2+H2反应的过渡态搜索时先用柔性扫描初猜过渡态 ...

上传完整输出文件
你给的这几行几乎没有有效信息，不能唯一确定报错原因

作者
Author: 王若曦 时间: 2023-3-25 16:42
老师，这是完整输出文件。

作者
Author: wzkchem5 时间: 2023-3-25 16:47

王若曦发表于 2023-3-25 09:42
老师，这是完整输出文件。

报错原因写得很清楚了
Optimization stopped.
-- Number of steps exceeded, NStep= 23

作者
Author: 王若曦 时间: 2023-3-25 16:58

wzkchem5 发表于 2023-3-25 16:47
报错原因写得很清楚了
Optimization stopped.
-- Number of steps exceeded, NStep= 23

谢谢老师，我下次认真看下输出文件！

作者
Author: rorssharsh 时间: 2023-5-27 10:52

小和尚发表于 2015-5-27 15:16
非常感谢啊！

您好，请问一下您是如何解决的，我也遇到了同样的问题，不过我是计算TD-DFT。

作者
Author: Yx0936 时间: 2023-7-28 18:48

rorssharsh 发表于 2023-5-27 10:52
您好，请问一下您是如何解决的，我也遇到了同样的问题，不过我是计算TD-DFT。

请问解决了吗

作者
Author: sobereva 时间: 2023-7-28 19:07

Yx0936 发表于 2023-7-28 18:48
请问解决了吗

仔细看
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

其他人有类似问题一律看此文

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