 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術(shù)科研互動平臺 » 計算模擬區(qū) » 量子化學 » Gaussian » 過渡態(tài)
| 8 | 1/1 | 返回列表 |
| 查看: 1492 | 回復: 7 | ||
lm214019銀蟲 (小有名氣)
|
[求助]
過渡態(tài)
|
| 計算過渡態(tài)時，一直出現(xiàn)9999錯誤（超過5個），應(yīng)該加個什么關(guān)鍵詞可以避免？ |
回復此樓» 猜你喜歡
基金正文30頁指的是報告正文還是整個申請書
已經(jīng)有9人回復
-
誠邀有志之士加入江蘇大學環(huán)境學院合成生物學課題組！(長期接收學碩專碩、碩博連讀)
已經(jīng)有16人回復
-
物理化學論文潤色/翻譯怎么收費?
已經(jīng)有53人回復
-
研究生招生
已經(jīng)有0人回復
-
一志愿南大理論與計算化學方向（070300）342分，求調(diào)劑
已經(jīng)有1人回復
-
歡迎化學、化工、物理、材料、或計算機背景的考研咨詢
已經(jīng)有0人回復
-
哈爾濱工業(yè)大學（深圳）綠色化學團隊招收2026年秋季入學博士生
已經(jīng)有11人回復
-
澳大利亞昆士蘭科技大學PhD博士獎學金
已經(jīng)有0人回復
-
山西大學資源與環(huán)境工程研究所化工分離/環(huán)境功能材料小組招收2026年博士生
已經(jīng)有2人回復
2樓2012-05-16 14:35:01
3樓2012-05-16 14:47:44
lm214019
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 389
- 散金: 20
- 帖子: 143
- 在線: 74.5小時
- 蟲號: 1155115
- 注冊: 2010-11-24
- 性別: MM
- 專業(yè): 理論和計算化學
|
我的輸入文件是這樣的： %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 |
4樓2012-05-16 15:00:09
lm214019
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 389
- 散金: 20
- 帖子: 143
- 在線: 74.5小時
- 蟲號: 1155115
- 注冊: 2010-11-24
- 性別: MM
- 專業(yè): 理論和計算化學
|
出錯部分： 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 |
5樓2012-05-16 15:05:01
baoa9992
木蟲 (初入文壇)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 1914.9
- 帖子: 34
- 在線: 141.2小時
- 蟲號: 611071
- 注冊: 2008-09-24
- 性別: GG
- 專業(yè): 理論和計算化學
6樓2012-05-17 08:42:16
meteoric30
木蟲 (正式寫手)
- 應(yīng)助: 43 (小學生)
- 金幣: 1719.5
- 散金: 12
- 紅花: 2
- 帖子: 436
- 在線: 320.7小時
- 蟲號: 1046846
- 注冊: 2010-06-24
- 性別: GG
- 專業(yè): 催化化學
7樓2012-05-18 09:33:45
fishdewow
銅蟲 (初入文壇)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 136.1
- 帖子: 9
- 在線: 34.7小時
- 蟲號: 1184708
- 注冊: 2011-01-06
- 專業(yè): 理論和計算化學
【答案】應(yīng)助回帖
★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
lm214019: 金幣+5, ★★★很有幫助 2012-05-18 10:41:37
感謝參與，應(yīng)助指數(shù) +1
lm214019: 金幣+5, ★★★很有幫助 2012-05-18 10:41:37
|
本班有帖子說過，optcyc設(shè)置為200，優(yōu)化不出來的話，盲目加大循環(huán)次數(shù)意義不大，幾何構(gòu)型優(yōu)化不收斂的解決方案。引起這個問題的原因比較多，可能也說不全，盡量做到全面。 (1) 如果是很小的分子（10原子以內(nèi)），初始結(jié)構(gòu)可能離平衡結(jié)構(gòu)比較遠，又在輸出文件優(yōu)化的最后一步判斷是否收斂的位置（可以通過查找Threshold字段找到，它下面有四個判斷項，都是YES才代表優(yōu)化收斂）看到了“-- Number of steps exceeded, NStep=xxx”，可以通過加大優(yōu)化的循環(huán)次數(shù)來解決問題。使用關(guān)鍵詞opt=(maxcycle=n). (2) 優(yōu)化到如果四個收斂標準中前兩項早就收斂了，而后兩項尤其是第三項不收斂，這需要判斷原因。如果每一步優(yōu)化過程中的能量始終在下降，那么可以繼續(xù)讓他算，超過了最大步數(shù)停掉了的話把結(jié)構(gòu)拿出來，重新提交就行了。如果能量忽大忽小出現(xiàn)跳躍，說明遇到了比較平坦的勢能面，或者優(yōu)化的算法不好。此時建議按如下順序處理： a. 使用opt=maxstep=n來縮小最大步長，即原子移動的距離，n的默認值是30，只能設(shè)置整數(shù)。 b. 使用關(guān)鍵詞opt=gdiis，很多時候可以很快收斂。 c. 使用opt=calcfc關(guān)鍵詞，在幾何優(yōu)化的第一步計算力常數(shù)，為優(yōu)化定一個大方向，類似在第一個三岔路口選一條路。 d. 使用opt=calcall關(guān)鍵詞，在幾何優(yōu)化計算的每一步都計算力常數(shù)，類似于在每一個三岔路口都來選一條路。這個關(guān)鍵詞不到萬不得已不要使用，非常耗時。 (3) 四個判斷標準的后兩項早就收斂了，但前兩項似乎固定在某個數(shù)值上了。遇到這種情況，首先去檢查你優(yōu)化第一步的時候給出的每個軌道得本征值部分，也就是軌道的能量，如果HOMO和LUMO的能量相同或者非常接近，那么你的設(shè)置有問題，即自旋多重度不對。如果HOMO-LUMO能量沒問題，那么將目前的計算停掉，提出結(jié)構(gòu)，使用opt=gdiis關(guān)鍵詞，很多時候管用。 (4) 如果是在幾何構(gòu)型優(yōu)化的過程中出現(xiàn)因為scf不收斂導致無法優(yōu)化，可以先用上面scf不收斂的解決方案來處理，如果都不行，可以通過在#后面添加p輸出更多信息，如果看到每次scf能量都能收斂到10的-6次方以下然后開始出現(xiàn)能量的跳躍，就是收斂不到10的-8次方（這種情況還比較常見），可以使用scf=(conver=6)來讓幾何優(yōu)化繼續(xù)下去，等優(yōu)化好了，再將此關(guān)鍵詞去掉，我的經(jīng)驗，優(yōu)化的結(jié)構(gòu)一般不會有太大的差異。如果不到10的-6能量就開始跳越，不適合此方法。 (5) 如果幾何優(yōu)化過程中發(fā)現(xiàn)鍵變得異常短，而且此時幾何優(yōu)化總是不收斂。遇到這種情況，問題出在你的基組使用，肯定是使用了贗勢基組卻沒有讀入贗勢。由于贗勢基組將內(nèi)層電子和原子核用一個有效核心勢來表示，你沒有讀入贗勢就說明沒有讀入這個有效核心勢，那么相當于把原子的內(nèi)層電子都給除掉了，巨大的核引力把原子之間的距離變得異常短也不足為奇。 上述解決方案轉(zhuǎn)自本版 http://www.gaoyang168.com/bbs/viewthread.php?tid=1566146 |
8樓2012-05-18 10:10:38
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 lm214019 的主題更新
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 269專碩求調(diào)劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
|
[考研] 考研調(diào)劑 +3 | 呼呼？~+123456 2026-03-21 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
|
[考研] 材料學學碩080502 337求調(diào)劑-一志愿華中科技大學 +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
|
[考研] 化學求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計算機技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
|
[考研] 南京大學化學376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿 南京航空航天大學大學 ，080500材料科學與工程學碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
|
[考研] 一志愿西安交通大學 學碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 0703化學調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿南京大學，080500材料科學與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|