| 24小時(shí)熱門版塊排行榜 |
| 8 | 1/1 | 返回列表 |
| 查看: 1493 | 回復(fù): 7 | ||
lm214019銀蟲(chóng) (小有名氣)
|
[求助]
過(guò)渡態(tài)
|
| 計(jì)算過(guò)渡態(tài)時(shí),一直出現(xiàn)9999錯(cuò)誤(超過(guò)5個(gè)),應(yīng)該加個(gè)什么關(guān)鍵詞可以避免? |


銀蟲(chóng) (小有名氣)
|
我的輸入文件是這樣的: %chk=5Tcuzsm-5.chk %mem=6MW %nproc=1 #N opt=(ts,z-matrix,noeigentest,nrscale) freq b3lyp/3-21g optcyc=300 Title Card Required 2 1 Si 0.00000000 0.00000000 0.00000000 Si 0.00000000 0.00000000 6.41770817 Al 1.29801685 0.00000000 3.21107229 O 0.33237084 0.00010657 4.73012718 O 0.34620934 -0.04882568 1.68408882 O 2.76575769 -1.21674610 3.16045999 Si 2.92082130 -2.93769984 3.07576355 Si 2.33863405 3.04541412 3.53816768 O 2.62399645 1.35919585 3.27713846 H 0.75407383 -1.08709958 -0.66331527 H 0.42704873 1.30726459 -0.54689002 H -1.45070396 -0.19388262 -0.21969879 H 1.00301304 0.85066069 7.09648737 H 0.09712464 -1.38610032 6.92730870 H -1.35773233 0.53666734 6.66050660 H 4.34792661 -3.31939068 2.98586063 H 2.41839893 -3.48644048 1.79636703 H 2.54807112 -3.59559693 4.34801696 H 3.30496835 3.73939673 2.65783471 H 2.57380310 3.28988594 4.97876793 H 0.92877746 3.22178395 3.12392369 Al 4.09876604 0.15116763 3.15956939 Cl 6.01828139 0.19570053 4.68140303 Cl 4.46132208 0.12056725 1.01014697 Mg -1.43180142 -0.02717497 3.22700087 Si 8.26631441 -0.28422264 4.12834776 C 6.81483360 -0.25970718 2.59567709 H 6.13745546 -1.00011244 2.96697067 H 6.35381359 0.70439162 2.64928705 H 7.06218171 -0.48070527 1.57838717 Cl 9.47291120 1.07354540 3.07450236 Cl 8.49773348 -0.03973570 6.20118913 Cl 9.01568130 -2.18821704 3.65584951 |
銀蟲(chóng) (小有名氣)
|
出錯(cuò)部分: Sum of Mulliken charges= 2.00000 Electronic spatial extent (au): Charge= 2.0000 electrons Dipole moment (field-independent basis, Debye): X= 8.7970 Y= 5.8364 Z= -1.2550 Tot= 10.6314 Quadrupole moment (field-independent basis, Debye-Ang): XX= -101.3779 YY= -152.1076 ZZ= -129.0925 XY= -2.8804 XZ= -40.2458 YZ= 1.6128 Traceless Quadrupole moment (field-independent basis, Debye-Ang): XX= 26.1481 YY= -24.5816 ZZ= -1.5665 XY= -2.8804 XZ= -40.2458 YZ= 1.6128 Octapole moment (field-independent basis, Debye-Ang**2): XXX= 11.3080 YYY= 60.5999 ZZZ= -58.2902 XYY= 72.0253 XXY= 39.9780 XXZ= 110.9472 XZZ= 74.3849 YZZ= 21.3764 YYZ= -39.7445 XYZ= -15.4819 Hexadecapole moment (field-independent basis, Debye-Ang**3): XXXX= -9084.0365 YYYY= -2720.6517 ZZZZ= -2253.0993 XXXY= -192.6681 XXXZ= -447.6416 YYYX= 80.3685 YYYZ= 0.6904 ZZZX= -245.0311 ZZZY= -26.2051 XXYY= -1943.6951 XXZZ= -1936.7095 YYZZ= -889.3074 XXYZ= -51.2733 YYXZ= -120.6486 ZZXY= 14.7472 N-N= 3.193312440948D+03 E-N=-1.564641506277D+04 KE= 3.902425821285D+03 THERE IS NO SUBJECT, HOWEVER COMPLEX, WHICH, IF STUDIED WITH PATIENCE AND INTELLIGIENCE WILL NOT BECOME MORE COMPLEX QUOTED BY D. GORDON ROHMAN Error termination request processed by link 9999. Error termination via Lnk1e in /home/g03/l9999.exe at Tue May 15 07:16:07 2012. Job cpu time: 0 days 10 hours 40 minutes 22.0 seconds. File lengths (MBytes): RWF= 55 Int= 0 D2E= 0 Chk= 21 Scr= 1 |
木蟲(chóng) (初入文壇)

木蟲(chóng) (正式寫手)
銅蟲(chóng) (初入文壇)
|
本班有帖子說(shuō)過(guò),optcyc設(shè)置為200,優(yōu)化不出來(lái)的話,盲目加大循環(huán)次數(shù)意義不大,幾何構(gòu)型優(yōu)化不收斂的解決方案。引起這個(gè)問(wèn)題的原因比較多,可能也說(shuō)不全,盡量做到全面。 (1) 如果是很小的分子(10原子以內(nèi)),初始結(jié)構(gòu)可能離平衡結(jié)構(gòu)比較遠(yuǎn),又在輸出文件優(yōu)化的最后一步判斷是否收斂的位置(可以通過(guò)查找Threshold字段找到,它下面有四個(gè)判斷項(xiàng),都是YES才代表優(yōu)化收斂)看到了“-- Number of steps exceeded, NStep=xxx”,可以通過(guò)加大優(yōu)化的循環(huán)次數(shù)來(lái)解決問(wèn)題。使用關(guān)鍵詞opt=(maxcycle=n). (2) 優(yōu)化到如果四個(gè)收斂標(biāo)準(zhǔn)中前兩項(xiàng)早就收斂了,而后兩項(xiàng)尤其是第三項(xiàng)不收斂,這需要判斷原因。如果每一步優(yōu)化過(guò)程中的能量始終在下降,那么可以繼續(xù)讓他算,超過(guò)了最大步數(shù)停掉了的話把結(jié)構(gòu)拿出來(lái),重新提交就行了。如果能量忽大忽小出現(xiàn)跳躍,說(shuō)明遇到了比較平坦的勢(shì)能面,或者優(yōu)化的算法不好。此時(shí)建議按如下順序處理: a. 使用opt=maxstep=n來(lái)縮小最大步長(zhǎng),即原子移動(dòng)的距離,n的默認(rèn)值是30,只能設(shè)置整數(shù)。 b. 使用關(guān)鍵詞opt=gdiis,很多時(shí)候可以很快收斂。 c. 使用opt=calcfc關(guān)鍵詞,在幾何優(yōu)化的第一步計(jì)算力常數(shù),為優(yōu)化定一個(gè)大方向,類似在第一個(gè)三岔路口選一條路。 d. 使用opt=calcall關(guān)鍵詞,在幾何優(yōu)化計(jì)算的每一步都計(jì)算力常數(shù),類似于在每一個(gè)三岔路口都來(lái)選一條路。這個(gè)關(guān)鍵詞不到萬(wàn)不得已不要使用,非常耗時(shí)。 (3) 四個(gè)判斷標(biāo)準(zhǔn)的后兩項(xiàng)早就收斂了,但前兩項(xiàng)似乎固定在某個(gè)數(shù)值上了。遇到這種情況,首先去檢查你優(yōu)化第一步的時(shí)候給出的每個(gè)軌道得本征值部分,也就是軌道的能量,如果HOMO和LUMO的能量相同或者非常接近,那么你的設(shè)置有問(wèn)題,即自旋多重度不對(duì)。如果HOMO-LUMO能量沒(méi)問(wèn)題,那么將目前的計(jì)算停掉,提出結(jié)構(gòu),使用opt=gdiis關(guān)鍵詞,很多時(shí)候管用。 (4) 如果是在幾何構(gòu)型優(yōu)化的過(guò)程中出現(xiàn)因?yàn)閟cf不收斂導(dǎo)致無(wú)法優(yōu)化,可以先用上面scf不收斂的解決方案來(lái)處理,如果都不行,可以通過(guò)在#后面添加p輸出更多信息,如果看到每次scf能量都能收斂到10的-6次方以下然后開(kāi)始出現(xiàn)能量的跳躍,就是收斂不到10的-8次方(這種情況還比較常見(jiàn)),可以使用scf=(conver=6)來(lái)讓幾何優(yōu)化繼續(xù)下去,等優(yōu)化好了,再將此關(guān)鍵詞去掉,我的經(jīng)驗(yàn),優(yōu)化的結(jié)構(gòu)一般不會(huì)有太大的差異。如果不到10的-6能量就開(kāi)始跳越,不適合此方法。 (5) 如果幾何優(yōu)化過(guò)程中發(fā)現(xiàn)鍵變得異常短,而且此時(shí)幾何優(yōu)化總是不收斂。遇到這種情況,問(wèn)題出在你的基組使用,肯定是使用了贗勢(shì)基組卻沒(méi)有讀入贗勢(shì)。由于贗勢(shì)基組將內(nèi)層電子和原子核用一個(gè)有效核心勢(shì)來(lái)表示,你沒(méi)有讀入贗勢(shì)就說(shuō)明沒(méi)有讀入這個(gè)有效核心勢(shì),那么相當(dāng)于把原子的內(nèi)層電子都給除掉了,巨大的核引力把原子之間的距離變得異常短也不足為奇。 上述解決方案轉(zhuǎn)自本版 http://www.gaoyang168.com/bbs/viewthread.php?tid=1566146 |
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
pk3725069 2026-03-19 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 291 求調(diào)劑 +4 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 4/200 |
|
|
[考研] 306求調(diào)劑 +7 | chuanzhu川燭 2026-03-18 | 7/350 |
|
|
[考研]
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
[考研] 315分,誠(chéng)求調(diào)劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
[考研] 求調(diào)劑 +6 | 十三加油 2026-03-21 | 6/300 |
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無(wú)聊 2026-03-21 | 4/200 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開(kāi)心 2026-03-16 | 4/200 |
|