 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3229 | 回復: 5 | ||||||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||||
LinaInverse木蟲 (小有名氣)
|
[求助]
計算時候總是不能收斂如何解決？
|
|||||
|
分子構型似乎沒發(fā)生明顯變化，但是就是不能收斂，急求各位大俠解答 給個截圖 |
回復此樓» 收錄本帖的淘帖專輯推薦
 過渡態(tài)計算 | 量化計算 | 第一性原理 | 計算化學 |
» 猜你喜歡
求調劑，本科吉大，一志愿華工339
已經有2人回復
-
意大利博洛尼亞大學招聘研究助理2名（2D材料/鋰離子與鈉離子電池方向）
已經有1人回復
-
物理化學論文潤色/翻譯怎么收費?
已經有80人回復
-
關于碳量子點抗菌的研究
已經有15人回復
-
南大化學調劑
已經有11人回復
-
0703調劑
已經有4人回復
-
1
已經有21人回復
-
0703 調劑
已經有5人回復
-
西南某211招27屆博士
已經有30人回復
-
大連工業(yè)杰青/長江團隊-生物質材料-儲能電池方向招收2026級博士生
已經有0人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- 關于VASP計算收斂的判斷
已經有25人回復
- 為什么會出現l9999錯誤及過渡態(tài)初始結構怎么設定？
已經有14人回復
- 怎么優(yōu)化這個函數
已經有5人回復
- 算例第一個時間步就不收斂？
已經有14人回復
- 四個收斂值為什么不收斂？
已經有6人回復
- 被積函數很快收斂到0，為什么積分不收斂。
已經有3人回復
- castep scf不收斂怎么辦，怎么辦��！ 急急急��！
已經有3人回復
- ansys接觸計算不收斂問題
已經有16人回復
- g09中優(yōu)化構型的能量如何表示
已經有7人回復
- 靜態(tài)計算能收斂，同樣的參數弛豫不能收斂，這是為什么？
已經有4人回復
- 自己使用QHA時計算退出
已經有23人回復
- 計算表面時：Vacuum真空層和slab層厚度收斂性測試，保持超晶胞體積不變嗎？？
已經有10人回復
- 如何畫化學實驗用具？燒杯、錐形瓶......
已經有15人回復
- siesta 表面不收斂
已經有3人回復
- 【求助】計算界面，總不收斂，請教如何處理？
已經有5人回復
- 【求助】用Orca計算一個片簇的單點能，不收斂。
已經有7人回復
- 【求助】多相流計算中只有連續(xù)性方程不收斂，求收斂策略
已經有38人回復
- 【求助】贗勢基組問題,急!
已經有5人回復
- 【求助】scan計算出現收斂失敗怎么辦/
已經有3人回復
- 【討論】dlpoly能量不收斂
已經有3人回復
LinaInverse
木蟲 (小有名氣)
- 應助: 3 (幼兒園)
- 金幣: 3419.6
- 紅花: 1
- 帖子: 98
- 在線: 58.5小時
- 蟲號: 1024554
- 注冊: 2010-05-21
- 性別: GG
- 專業(yè): 配位化學
|
居然忘了…… 不過不是你說的問題。dft，ub3lyp基組 >>>>>>>>>> Convergence criterion not met. SCF Done: E(UB3LYP) = -2200.00081496 A.U. after 351 cycles Convg = 0.2434D-03 -V/T = 2.0293 KE= 2.137326153433D+03 PE=-8.831064536717D+03 EE= 2.691309091686D+03 Annihilation of the first spin contaminant: S**2 before annihilation 6.0241, after 6.0002 Convergence failure -- run terminated. Error termination via Lnk1e in /gpfssan1/home/hyan//g09/l502.exe at Wed Dec 28 17:12:48 2011. |
3樓2012-01-02 15:56:30
nyzhaoyin
新蟲 (正式寫手)
- 應助: 27 (小學生)
- 金幣: 1044.9
- 散金: 112
- 紅花: 9
- 帖子: 681
- 在線: 241.5小時
- 蟲號: 563016
- 注冊: 2008-05-24
- 性別: GG
2樓2012-01-02 12:43:10
nyzhaoyin
新蟲 (正式寫手)
- 應助: 27 (小學生)
- 金幣: 1044.9
- 散金: 112
- 紅花: 9
- 帖子: 681
- 在線: 241.5小時
- 蟲號: 563016
- 注冊: 2008-05-24
- 性別: GG
【答案】應助回帖
★ ★
感謝參與，應助指數 +1
zhou2009(金幣+2): 2012-01-02 21:00:56
感謝參與，應助指數 +1
zhou2009(金幣+2): 2012-01-02 21:00:56
|
（2.3） scf不收斂怎么辦？幾何優(yōu)化不收斂怎么辦？引用回帖: 雖然二者有聯系，但實際上是兩個相對獨立的問題，要分開回答。 這里提供的只是解決方案，所有的方法都用上了還是不收斂得情況也是有的。 首先要分清scf不收斂和幾何構型優(yōu)化不收斂：scf不收斂指的是自洽場疊代不收斂（什么？沒聽說過什么叫自洽場？那還是回去學習些量化基礎知識再開展計算吧�。�，可以認為是對指定結構的波函數不斷優(yōu)化的過程，是為了找到這個某個指定結構下能量最低的波函數，而幾何構型優(yōu)化是對結構的優(yōu)化的過程，是為了找到某個指定的組分下能量極小結構（注意，不一定是能量最小結構）。在量子化學計算的幾何構型優(yōu)化中，每一步的幾何構型優(yōu)化都包含的很多次的scf計算。 1、scf不收斂的解決方案。 (1) 可以加大scf的循環(huán)次數，默認的循環(huán)次數是128次，通過scf=(maxcycle=n)來設置最大循環(huán)次數n。建議不要超過512，更多的循換沒有必要。 (2) 如果加大循環(huán)次數不管用，在分子有對稱性的情況下，使用scf=dsymm關鍵詞來強制密度對稱，有時可以收斂。另外，此關鍵詞很多時候對"scf is confused”這種錯誤很管用。 (3) 使用scf=symm關鍵詞，使用的前提同上，有時可以收斂。 (4) 如果(2)(3)兩步都不行，可以將對稱的分子中的某幾個原子的位置微調，使分子喪失對稱性。這等效于nosymm關鍵詞，但個人經驗，這種方式比nosymm好用的多。 (5) 如果還不行，只能拿出殺手锏了，就是使用qc，但不建議直接使用，而是使用xqc關鍵詞，比如scf=(maxcycle=80,xqc)，意思是如果scf正常計算(dc)在80個循環(huán)之內不收斂才進行昂貴的qc計算，因為scf不收斂多數在幾個優(yōu)化的過程中出現，無法判斷哪一步優(yōu)化的時候會出現scf不收斂，所以用xqc比純粹使用qc要省時的多。 (6) 中級用戶可以在輸入文件的井號“#”開頭那一行井號后面加上字母"p"來輸出更多的信息，其中就有自洽場疊代的信息，分析原因可能會對采用什么方法提供指導。 (7) 前面有蟲子提到一個有用的方案但沒說清楚，我這里補充一下：如果用用小基組計算，scf可以收斂，那么保存好檢查點文件，換成大基組的時候從檢查點文件中讀取初始猜測(使用guess=read關鍵詞)，有時可以算過去。 (8) 如果你計算中使用了含有彌散函數的基組（比如6-31+G*基組中的“+”，Aug-cc-pVTZ基組中的Aug等），可以先去掉彌散函數（比如例子中的兩個基組變?yōu)?-31G*和cc-pVTZ），這樣很容易收斂，需要的話，換成帶有彌散函數的基組時，從檢查點文件的讀初始猜測。原理與(7)類似，有時候管用。 注意：(7)和(8)不適合很大的體系，此時，使用小基組或者不帶彌散函數的基組與后來使用的大基組在基組數量上差得太多（如果你本來就遇到了scf收斂問題，而使用這兩個方案時，機組數量差距超過512個，多數時候要出問題），反倒容易引起scf不收斂，因為讀檢查點文件時，對那些多出來的基組的猜測是空的，而直接計算時最起碼還有個半經驗的計算提供了相對合理的初始猜測。 (9) 上述方法有時可以組合使用，如果經過各種組合處理都不收斂，那么放棄吧，你的分子的電子結構太差了。 2、幾何構型優(yōu)化不收斂的解決方案。引起這個問題的原因比較多，可能也說不全，盡量做到全面。 (1) 如果是很小的分子（10原子以內），初始結構可能離平衡結構比較遠，又在輸出文件優(yōu)化的最后一步判斷是否收斂的位置（可以通過查找Threshold字段找到，它下面有四個判斷項，都是YES才代表優(yōu)化收斂）看到了“-- Number of steps exceeded, NStep=xxx”，可以通過加大優(yōu)化的循環(huán)次數來解決問題。使用關鍵詞opt=(maxcycle=n). (2) 優(yōu)化到如果四個收斂標準中前兩項早就收斂了，而后兩項尤其是第三項不收斂，這需要判斷原因。如果每一步優(yōu)化過程中的能量始終在下降，那么可以繼續(xù)讓他算，超過了最大步數停掉了的話把結構拿出來，重新提交就行了。如果能量忽大忽小出現跳躍，說明遇到了比較平坦的勢能面，或者優(yōu)化的算法不好。此時建議按如下順序處理： a. 使用opt=maxstep=n來縮小最大步長，即原子移動的距離，n的默認值是30，只能設置整數。 b. 使用關鍵詞opt=gdiis，很多時候可以很快收斂。 c. 使用opt=calcfc關鍵詞，在幾何優(yōu)化的第一步計算力常數，為優(yōu)化定一個大方向，類似在第一個三岔路口選一條路。 d. 使用opt=calcall關鍵詞，在幾何優(yōu)化計算的每一步都計算力常數，類似于在每一個三岔路口都來選一條路。這個關鍵詞不到萬不得已不要使用，非常耗時。 (3) 四個判斷標準的后兩項早就收斂了，但前兩項似乎固定在某個數值上了。遇到這種情況，首先去檢查你優(yōu)化第一步的時候給出的每個軌道得本征值部分，也就是軌道的能量，如果HOMO和LUMO的能量相同或者非常接近，那么你的設置有問題，即自旋多重度不對。如果HOMO-LUMO能量沒問題，那么將目前的計算停掉，提出結構，使用opt=gdiis關鍵詞，很多時候管用。 (4) 如果是在幾何構型優(yōu)化的過程中出現因為scf不收斂導致無法優(yōu)化，可以先用上面scf不收斂的解決方案來處理，如果都不行，可以通過在#后面添加p輸出更多信息，如果看到每次scf能量都能收斂到10的-6次方以下然后開始出現能量的跳躍，就是收斂不到10的-8次方（這種情況還比較常見），可以使用scf=(conver=6)來讓幾何優(yōu)化繼續(xù)下去，等優(yōu)化好了，再將此關鍵詞去掉，我的經驗，優(yōu)化的結構一般不會有太大的差異。如果不到10的-6能量就開始跳越，不適合此方法。 (5) 如果幾何優(yōu)化過程中發(fā)現鍵變得異常短，而且此時幾何優(yōu)化總是不收斂。遇到這種情況，問題出在你的基組使用，肯定是使用了贗勢基組卻沒有讀入贗勢。由于贗勢基組將內層電子和原子核用一個有效核心勢來表示，你沒有讀入贗勢就說明沒有讀入這個有效核心勢，那么相當于把原子的內層電子都給除掉了，巨大的核引力把原子之間的距離變得異常短也不足為奇。 比較通俗的解答請看2.10 scf加xqc吧 或者nodiis |
4樓2012-01-02 16:18:07
nyzhaoyin
新蟲 (正式寫手)
- 應助: 27 (小學生)
- 金幣: 1044.9
- 散金: 112
- 紅花: 9
- 帖子: 681
- 在線: 241.5小時
- 蟲號: 563016
- 注冊: 2008-05-24
- 性別: GG
5樓2012-01-02 16:19:55
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿西安交通大學材料工程專業(yè) 282分求調劑 +10 | 楓橋ZL 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 材料工程專碩 348分求調劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
|
[考研] 工科0856求調劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
|
[考研] 313求調劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 307求調劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
|
[考研] 材料 271求調劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
|
[考研] 330求調劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
|
[考研] 材料學學碩080502 337求調劑-一志愿華中科技大學 +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
|
[考研] 求調劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
|
[考研] 346求調劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
|
[考研] 304求調劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
|
[考研] 304求調劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
|
[考研] 290求調劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
|
[考研] 求調劑一志愿南京航空航天大學289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 材料學碩297已過四六級求調劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
|
[考研] 生物學調劑招人�。�！ +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 277調劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 308求調劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 有沒有道鐵/土木的想調劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|