版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

Yan_Jordan

木蟲 (小有名氣)

QC強帖: 1
應助: 55 (初中生)
金幣: 1944
散金: 221
紅花: 11
帖子: 169
在線: 108.2小時
蟲號: 6034100
注冊: 2017-03-17
性別: GG
專業(yè): 理論和計算化學

【答案】應助回帖

感謝參與，應助指數(shù) +1

你這應該是優(yōu)化計算不收斂，不知道我下不了你的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

這樣的語句吧，
你的問題應該是優(yōu)化不收斂，原因有好幾種可能：（給你貼一些建議）
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).

我的建議是，你先不做優(yōu)化，直接坐energy計算，如果不能收斂，那說明你的初始結構是有問題的。其他則可以參考上面貼出來的信息

贊一下

回復此樓

» 本帖已獲得的紅花（最新10朵）

wittymadder

11樓2017-04-04 15:51:52

已閱關注TA 給TA發(fā)消息送TA紅花 TA的回帖

wuy069

鐵桿木蟲 (正式寫手)

應助: 20 (小學生)
金幣: 7921.9
散金: 15
紅花: 10
帖子: 506
在線: 549.1小時
蟲號: 760723
注冊: 2009-04-30
性別: GG
專業(yè): 理論和計算化學

還有一個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發(fā)消息送TA紅花 TA的回帖

wuy069

鐵桿木蟲 (正式寫手)

應助: 20 (小學生)
金幣: 7921.9
散金: 15
紅花: 10
帖子: 506
在線: 549.1小時
蟲號: 760723
注冊: 2009-04-30
性別: GG
專業(yè): 理論和計算化學

【答案】應助回帖

感謝參與，應助指數(shù) +1

根據你的輸入文件，已經給你計算完了，傳到百度網盤了
http://pan.baidu.com/s/1migCtx2

贊一下

回復此樓

» 本帖已獲得的紅花（最新10朵）

wittymadder

13樓2017-04-04 16:51:22

已閱關注TA 給TA發(fā)消息送TA紅花 TA的回帖

wittymadder

新蟲 (正式寫手)

應助: 1 (幼兒園)
金幣: 2266.8
散金: 150
紅花: 2
帖子: 400
在線: 54.9小時
蟲號: 4440708
注冊: 2016-02-27
性別: MM
專業(yè): 高分子材料的加工與成型

求頂

發(fā)自小木蟲Android客戶端

回復此樓

14樓2017-04-05 19:42:51

已閱關注TA 給TA發(fā)消息送TA紅花 TA的回帖

wittymadder

新蟲 (正式寫手)

應助: 1 (幼兒園)
金幣: 2266.8
散金: 150
紅花: 2
帖子: 400
在線: 54.9小時
蟲號: 4440708
注冊: 2016-02-27
性別: MM
專業(yè): 高分子材料的加工與成型

送紅花一朵

引用回帖:

13樓: Originally posted by wuy069 at 2017-04-04 16:51:22
根據你的輸入文件，已經給你計算完了，傳到百度網盤了
http://pan.baidu.com/s/1migCtx2

感謝天使～請問輸入文件是改成這樣嗎？運行了一下如圖二不知道有沒有開始
（另外還好想知道你怎么算可以在5分鐘算完）

發(fā)自小木蟲Android客戶端

贊一下

回復此樓

15樓2017-04-05 20:23:31

已閱關注TA 給TA發(fā)消息送TA紅花 TA的回帖

wittymadder

新蟲 (正式寫手)

應助: 1 (幼兒園)
金幣: 2266.8
散金: 150
紅花: 2
帖子: 400
在線: 54.9小時
蟲號: 4440708
注冊: 2016-02-27
性別: MM
專業(yè): 高分子材料的加工與成型

送紅花一朵

引用回帖:

11樓: Originally posted by Yan_Jordan at 2017-04-04 15:51:52
你這應該是優(yōu)化計算不收斂，不知道我下不了你的out文件，被瀏覽器禁了，不過我猜測你的文件這個錯誤信息前面應該有類似：
Maximum Force 0.020301 0.000450 NO
RMS Force 0. ...

感謝～我先算能量看看（原諒英語渣看不懂英文）

發(fā)自小木蟲Android客戶端

贊一下

回復此樓

16樓2017-04-05 20:28:41

已閱關注TA 給TA發(fā)消息送TA紅花 TA的回帖

paramecium86

版主 (著名寫手)

QC強帖: 5
應助: 588 (博士)
貴賓: 0.01
金幣: 26384.6
紅花: 349
帖子: 2026
在線: 1442.3小時
蟲號: 2466443
注冊: 2013-05-16
性別: GG
專業(yè): 理論和計算化學
管轄: 量子化學

【答案】應助回帖

看了輸出文件提示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 再試試

如果加多步數(shù)不行就可以試試 opt里面不用modredandunt 而使用calcall 這個會很慢但是有時候能幫助找到很難找的TS態(tài)

贊一下

回復此樓

17樓2017-04-05 22:56:13

已閱關注TA 給TA發(fā)消息送TA紅花 TA的回帖

wuy069

鐵桿木蟲 (正式寫手)

應助: 20 (小學生)
金幣: 7921.9
散金: 15
紅花: 10
帖子: 506
在線: 549.1小時
蟲號: 760723
注冊: 2009-04-30
性別: GG
專業(yè): 理論和計算化學

引用回帖:

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文件，你就不用寫幾何了，電荷和多重度寫上就行。當然你也可以把我給你優(yōu)化好的幾何直接取出來，你自己再優(yōu)化下，應該很快就收斂了

贊一下

回復此樓

18樓2017-04-06 16:16:59

已閱關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉我要訂閱樓主 wittymadder 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 268求調劑 +3	簡單點0 2026-03-17	3/150	2026-03-17 23:15 by wangjihu
[考研] 能源材料化學課題組招收碩士研究生8-10名 +3	脫穎而出 2026-03-16	6/300	2026-03-17 21:19 by z1z2z3879
[考研] 326求調劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[碩博家園] 湖北工業(yè)大學生命科學與健康學院-課題組招收2026級食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 332求調劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 梁成偉老師課題組歡迎你的加入 +8	一鴨鴨喲 2026-03-14	10/500	2026-03-17 15:07 by 一鴨鴨喲
[考研] 085600材料與化工 +4	安全上岸！ 2026-03-16	4/200	2026-03-17 14:02 by 勇敢太監(jiān)王公公
[考研] 302求調劑 +4	小賈同學123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學123
[碩博家園] 深圳大學碩士招生（2026秋，傳感器方向，僅錄取第一志愿） +4	xujiaoszu 2026-03-11	9/450	2026-03-17 10:29 by xujiaoszu
[考研] 材料與化工304求B區(qū)調劑 +7	邱gl 2026-03-11	8/400	2026-03-17 09:36 by 努力學習賺彩禮
[考研] 274求調劑 +5	時間點 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 318求調劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 0703一志愿211 285分求調劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 求老師收留調劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 289求調劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 308 085701 四六級已過求調劑 +7	溫喬喬喬喬 2026-03-12	14/700	2026-03-14 10:49 by JourneyLucky
[考研] 一志愿中科院，化學方向，295求調劑 +4	一氧二氮 2026-03-11	4/200	2026-03-13 22:35 by JourneyLucky
[考研] 材料工程調劑 +9	咪咪空空 2026-03-12	9/450	2026-03-13 22:05 by 星空星月
[考研] 290求調劑 +9	ADT 2026-03-11	9/450	2026-03-13 21:55 by JourneyLucky