应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 关于高斯不能正常终止求助
182/2返回列表上一页12
查看: 2038  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

Yan_Jordan

木虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
你这应该是优化计算不收敛，不知道我下不了你的out文件，被浏览器禁了，不过我猜测你的文件这个错误信息前面应该有类似：
Maximum Force            0.020301     0.000450     NO
RMS     Force            0.007068     0.000300     NO
Maximum Displacement     0.078972     0.001800     NO
RMS     Displacement     0.023716     0.001200     NO

这样的语句吧，
你的问题应该是 优化不收敛，原因有好几种可能：（给你贴一些建议）
Geometry optimizations usually fail to converge for one of a few reasons.
If your initial starting structure is not good, then you should ask if you can provide a better starting structure, for example, one optimized at a lower level of theory. However, if it looks as if the structure is converging to what you want, as seen in your visualizer of choice, then one should restart the optimization from the last step, for example by using geom=allcheck in the route line, and it generally is a good idea to also use opt=CalcFC in these situations if it is not too expensive (say HF or DFT).
If your starting force constants matrix (Hessian) is poor, use a better one. This typically manifests itself when they vary a lot between levels, or if there is a large geometry change during the optimization. One can carry out a series of linked jobs (--Link1--). If you have a previous job, then usually Opt=ReadFC works well, but occasionally Opt=CalcFC, or rarely Opt=CalcAll are needed. In these cases, the forces are often converged, but the steps are not, and the final output will look like
         Item               Value     Threshold  Converged?
Maximum Force            0.000401     0.000450     YES
RMS     Force            0.000178     0.000300     YES
Maximum Displacement     0.010503     0.001800     NO
RMS     Displacement     0.003163     0.001200     NO
Rarely, the coordinate system itself may be at fault. If z-matrix coordinates are being used, then in some cases, poor choices can be made, which result in angles, or three consecutive atoms of the four atoms used to define a torsion angle (dihedral angle) may be collinear (the angle is close to 0 or 180 degrees), which can give problems. In this case one can either formulate a better z-matrix or use the default redundant internal coordinates.
If these methods fail, another option would be to change the optimization method from the default to another type, such as opt=ef (if the number of variables is less than 50) or opt=gdiis (for floppy molecules).

我的建议是，你先不做优化，直接坐energy计算，如果不能收敛，那说明你的初始结构是有问题的。其他则可以参考上面贴出来的信息
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

wittymadder
11楼2017-04-04 15:51:52
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wuy069

铁杆木虫 (正式写手)
还有一个l9999错误吧
Optimization stopped.   -- Number of steps exceeded, NStep= N [..] Error termination request processed by link 9999. Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l9999.exe

Solution: Maximum number of optimization steps is twice the number of variables to be optimized. Try increasing the value by specifying OPT=(MaxCycle=N) in .inp file, where N is the number of optimization steps (see OPT keyword). Alternatively, try to start optimization from different geometry.
赞一下 回复此楼
12楼2017-04-04 16:18:34
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wuy069

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
根据你的输入文件，已经给你计算完了，传到百度网盘了
http://pan.baidu.com/s/1migCtx2
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

wittymadder
13楼2017-04-04 16:51:22
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wittymadder

新虫 (正式写手)
求顶

发自小木虫Android客户端
回复此楼
14楼2017-04-05 19:42:51
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wittymadder

新虫 (正式写手)
送红花一朵
引用回帖:
13楼: Originally posted by wuy069 at 2017-04-04 16:51:22
根据你的输入文件，已经给你计算完了，传到百度网盘了
http://pan.baidu.com/s/1migCtx2

感谢天使～请问输入文件是改成这样吗？运行了一下如图二不知道有没有开始
（另外还好想知道你怎么算可以在5分钟算完）
关于高斯不能正常终止求助


关于高斯不能正常终止求助-1



发自小木虫Android客户端
赞一下 回复此楼
15楼2017-04-05 20:23:31
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wittymadder

新虫 (正式写手)
送红花一朵
引用回帖:
11楼: Originally posted by Yan_Jordan at 2017-04-04 15:51:52
你这应该是优化计算不收敛，不知道我下不了你的out文件，被浏览器禁了，不过我猜测你的文件这个错误信息前面应该有类似：
Maximum Force            0.020301     0.000450     NO
RMS     Force            0. ...

感谢～我先算能量看看（原谅英语渣看不懂英文）

发自小木虫Android客户端
赞一下 回复此楼
16楼2017-04-05 20:28:41
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

看了输出文件 提示link9999，说明是opt没有converge。具体可以查看36800行左右最后一个opt cycle的结果是
  Item                                   Value     Threshold  Converged?
Maximum Force                  0.006337     0.000450     NO
RMS     Force                       0.001761     0.000300     NO
Maximum Displacement     0.036834     0.001800     NO
RMS     Displacement         0.006511     0.001200      NO

Predicted change in Energy=-1.049415D-03
Optimization stopped.
    -- Number of steps exceeded,  NStep= 100
    -- Flag reset to prevent archiving

可能是猜的ts结构不够准确导致的。

检查每一个cycle的能量一直是降低的，没有振荡 所以可以尝试opt里面设置maxcyc=250 再试试

如果加多步数不行就可以试试 opt里面不用modredandunt 而使用calcall 这个会很慢但是有时候能帮助找到很难找的TS态
赞一下 回复此楼
17楼2017-04-05 22:56:13
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wuy069

铁杆木虫 (正式写手)
引用回帖:
15楼: Originally posted by wittymadder at 2017-04-05 20:23:31
感谢天使～请问输入文件是改成这样吗？运行了一下如图二不知道有没有开始
（另外还好想知道你怎么算可以在5分钟算完）


...

第一步l9999后，第二部是直接读的chk文件，
命令行# opt(ts,modredun,maxstep=50,maxcycle=200,noeigen) freq b3lyp/6-31g(d) geom=check guess=read

这样读chk文件，你就不用写几何了，电荷和多重度写上就行。当然你也可以把我给你优化好的几何直接取出来，你自己再优化下，应该很快就收敛了
赞一下 回复此楼
18楼2017-04-06 16:16:59
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wittymadder 的主题更新
182/2返回列表上一页12
信息提示
关闭
请填处理意见
关闭