 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1976 | 回復(fù): 5 | ||
vigaryang木蟲 (正式寫手)
|
[求助]
溶劑化構(gòu)型優(yōu)化問題
|
|
在溶劑化環(huán)境下優(yōu)化總是很難得到優(yōu)化構(gòu)型（SCF可收斂），看來一下log文件，一般都是四條標(biāo)準(zhǔn)中的第二條YES，而其它三條是NO。請問對于這種情況，最有效的解決辦法是什么？謝謝！ 試過iop(1/8)=5或者opt(maxstep=5)等方法，不太管用。 |
回復(fù)此樓» 猜你喜歡
廣西師范大學(xué)化學(xué)與藥學(xué)學(xué)院楊秀林課題組2026年擬招收博士研究生1名
已經(jīng)有25人回復(fù)
-
26年儲能、電池方向博士申請
已經(jīng)有1人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有66人回復(fù)
-
求調(diào)劑，本科吉大，一志愿華工339
已經(jīng)有2人回復(fù)
-
意大利博洛尼亞大學(xué)招聘研究助理2名（2D材料/鋰離子與鈉離子電池方向）
已經(jīng)有1人回復(fù)
-
關(guān)于碳量子點抗菌的研究
已經(jīng)有15人回復(fù)
-
南大化學(xué)調(diào)劑
已經(jīng)有11人回復(fù)
-
0703調(diào)劑
已經(jīng)有2人回復(fù)
-
1
已經(jīng)有14人回復(fù)
-
0703 調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- ds中的分子優(yōu)化
已經(jīng)有5人回復(fù)
- 怎么能用量子化學(xué)計算含氮雜環(huán)芳香化合物與水分子的相互作用
已經(jīng)有17人回復(fù)
- Gaussian對分子構(gòu)型優(yōu)化加溶劑化與否對最后結(jié)果有影響嗎？
已經(jīng)有8人回復(fù)
- 請教初始結(jié)構(gòu)優(yōu)化問題
已經(jīng)有10人回復(fù)
- 【求助】優(yōu)化：怎樣找到一個兩個分子之間的相對位置
已經(jīng)有26人回復(fù)
- 【求助】過渡態(tài)能計算溶劑化效應(yīng)嗎
已經(jīng)有4人回復(fù)
- 【求助】用ADF優(yōu)化γ【Mo8O26】4-構(gòu)型不理想
已經(jīng)有7人回復(fù)
- 【活動】征集：高斯09能做但高斯03做不了的研究實例
已經(jīng)有8人回復(fù)
書萬里
木蟲 (著名寫手)
歷史控+原著控
- 應(yīng)助: 70 (初中生)
- 金幣: 3295.6
- 散金: 742
- 紅花: 8
- 沙發(fā): 4
- 帖子: 1519
- 在線: 466.8小時
- 蟲號: 1244505
- 注冊: 2011-03-25
- 性別: GG
- 專業(yè): 中國通史
2樓2012-08-30 18:53:55
lishijunzong
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 214 (大學(xué)生)
- 金幣: 5034
- 紅花: 10
- 帖子: 574
- 在線: 195.1小時
- 蟲號: 1121020
- 注冊: 2010-10-13
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
3樓2012-08-31 12:49:31
yangwenpeng
新蟲 (正式寫手)
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 670.5
- 紅花: 5
- 帖子: 583
- 在線: 91.5小時
- 蟲號: 1468602
- 注冊: 2011-10-31
- 專業(yè): 電化學(xué)
【答案】應(yīng)助回帖
★ ★
gmy1990: 金幣+2 2012-09-05 15:23:06
gmy1990: 金幣+2 2012-09-05 15:23:06
|
(1) 可以加大scf的循環(huán)次數(shù)，默認(rèn)的循環(huán)次數(shù)是128次，通過scf(maxcycle=n)來設(shè)置最大循環(huán)次數(shù)n。建議不要超過512，更多的循換沒有必要。 (2) 如果加大循環(huán)次數(shù)不管用，在分子有對稱性的情況下，使用scf(dsymm)關(guān)鍵詞來強(qiáng)制密度對稱，有時可以收斂。另外，此關(guān)鍵詞很多時候?qū)?quot;scf is confused”這種錯誤很管用。 (3) 使用scf(symm)關(guān)鍵詞，使用的前提同上，有時可以收斂。 (4) 如果(2)(3)兩步都不行，可以將對稱的分子中的某幾個原子的位置微調(diào)，使分子喪失對稱性。這等效于nosymm關(guān)鍵詞，但個人經(jīng)驗，這種方式比nosymm好用的多。 |
4樓2012-09-05 11:21:56
 |
5樓2012-09-06 08:09:43
hlq3175
木蟲 (正式寫手)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 2378.8
- 散金: 107
- 紅花: 4
- 帖子: 444
- 在線: 301.2小時
- 蟲號: 870062
- 注冊: 2009-10-12
- 性別: MM
- 專業(yè): 高分子合成化學(xué)
【答案】應(yīng)助回帖
|
scf不收斂怎么辦？幾何優(yōu)化不收斂怎么辦？引用回帖: 雖然二者有聯(lián)系，但實際上是兩個相對獨(dú)立的問題，要分開回答。 這里提供的只是解決方案，所有的方法都用上了還是不收斂得情況也是有的。 首先要分清scf不收斂和幾何構(gòu)型優(yōu)化不收斂：scf不收斂指的是自洽場疊代不收斂（什么？沒聽說過什么叫自洽場？那還是回去學(xué)習(xí)些量化基礎(chǔ)知識再開展計算吧�。�，可以認(rèn)為是對指定結(jié)構(gòu)的波函數(shù)不斷優(yōu)化的過程，是為了找到這個某個指定結(jié)構(gòu)下能量最低的波函數(shù)，而幾何構(gòu)型優(yōu)化是對結(jié)構(gòu)的優(yōu)化的過程，是為了找到某個指定的組分下能量極小結(jié)構(gòu)（注意，不一定是能量最小結(jié)構(gòu)）。在量子化學(xué)計算的幾何構(gòu)型優(yōu)化中，每一步的幾何構(gòu)型優(yōu)化都包含的很多次的scf計算。 1、scf不收斂的解決方案。 (1) 可以加大scf的循環(huán)次數(shù)，默認(rèn)的循環(huán)次數(shù)是128次，通過scf=(maxcycle=n)來設(shè)置最大循環(huán)次數(shù)n。建議不要超過512，更多的循換沒有必要。 (2) 如果加大循環(huán)次數(shù)不管用，在分子有對稱性的情況下，使用scf=dsymm關(guān)鍵詞來強(qiáng)制密度對稱，有時可以收斂。另外，此關(guān)鍵詞很多時候?qū)?quot;scf is confused”這種錯誤很管用。 (3) 使用scf=symm關(guān)鍵詞，使用的前提同上，有時可以收斂。 (4) 如果(2)(3)兩步都不行，可以將對稱的分子中的某幾個原子的位置微調(diào)，使分子喪失對稱性。這等效于nosymm關(guān)鍵詞，但個人經(jīng)驗，這種方式比nosymm好用的多。 (5) 如果還不行，只能拿出殺手锏了，就是使用qc，但不建議直接使用，而是使用xqc關(guān)鍵詞，比如scf=(maxcycle=80,xqc)，意思是如果scf正常計算(dc)在80個循環(huán)之內(nèi)不收斂才進(jìn)行昂貴的qc計算，因為scf不收斂多數(shù)在幾個優(yōu)化的過程中出現(xiàn)，無法判斷哪一步優(yōu)化的時候會出現(xiàn)scf不收斂，所以用xqc比純粹使用qc要省時的多。 (6) 中級用戶可以在輸入文件的井號“#”開頭那一行井號后面加上字母"p"來輸出更多的信息，其中就有自洽場疊代的信息，分析原因可能會對采用什么方法提供指導(dǎo)。 (7) 前面有蟲子提到一個有用的方案但沒說清楚，我這里補(bǔ)充一下：如果用用小基組計算，scf可以收斂，那么保存好檢查點文件，換成大基組的時候從檢查點文件中讀取初始猜測(使用guess=read關(guān)鍵詞)，有時可以算過去。 (8) 如果你計算中使用了含有彌散函數(shù)的基組（比如6-31+G*基組中的“+”，Aug-cc-pVTZ基組中的Aug等），可以先去掉彌散函數(shù)（比如例子中的兩個基組變?yōu)?-31G*和cc-pVTZ），這樣很容易收斂，需要的話，換成帶有彌散函數(shù)的基組時，從檢查點文件的讀初始猜測。原理與(7)類似，有時候管用。 注意：(7)和(8)不適合很大的體系，此時，使用小基組或者不帶彌散函數(shù)的基組與后來使用的大基組在基組數(shù)量上差得太多（如果你本來就遇到了scf收斂問題，而使用這兩個方案時，機(jī)組數(shù)量差距超過512個，多數(shù)時候要出問題），反倒容易引起scf不收斂，因為讀檢查點文件時，對那些多出來的基組的猜測是空的，而直接計算時最起碼還有個半經(jīng)驗的計算提供了相對合理的初始猜測。 (9) 上述方法有時可以組合使用，如果經(jīng)過各種組合處理都不收斂，那么放棄吧，你的分子的電子結(jié)構(gòu)太差了。 2、幾何構(gòu)型優(yōu)化不收斂的解決方案。引起這個問題的原因比較多，可能也說不全，盡量做到全面。 (1) 如果是很小的分子（10原子以內(nèi)），初始結(jié)構(gòu)可能離平衡結(jié)構(gòu)比較遠(yuǎn)，又在輸出文件優(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)化到如果四個收斂標(biāo)準(zhǔn)中前兩項早就收斂了，而后兩項尤其是第三項不收斂，這需要判斷原因。如果每一步優(yōu)化過程中的能量始終在下降，那么可以繼續(xù)讓他算，超過了最大步數(shù)停掉了的話把結(jié)構(gòu)拿出來，重新提交就行了。如果能量忽大忽小出現(xiàn)跳躍，說明遇到了比較平坦的勢能面，或者優(yōu)化的算法不好。此時建議按如下順序處理： a. 使用opt=maxstep=n來縮小最大步長，即原子移動的距離，n的默認(rè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) 四個判斷標(biāo)準(zhǔn)的后兩項早就收斂了，但前兩項似乎固定在某個數(shù)值上了。遇到這種情況，首先去檢查你優(yōu)化第一步的時候給出的每個軌道得本征值部分，也就是軌道的能量，如果HOMO和LUMO的能量相同或者非常接近，那么你的設(shè)置有問題，即自旋多重度不對。如果HOMO-LUMO能量沒問題，那么將目前的計算停掉，提出結(jié)構(gòu)，使用opt=gdiis關(guān)鍵詞，很多時候管用。 (4) 如果是在幾何構(gòu)型優(yōu)化的過程中出現(xiàn)因為scf不收斂導(dǎo)致無法優(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)層電子和原子核用一個有效核心勢來表示，你沒有讀入贗勢就說明沒有讀入這個有效核心勢，那么相當(dāng)于把原子的內(nèi)層電子都給除掉了，巨大的核引力把原子之間的距離變得異常短也不足為奇。 比較通俗的解答請看2.10 scf加xqc吧 或者nodiis 轉(zhuǎn)自http://www.gaoyang168.com/bbs/viewthread.php?tid=4007785 |
6樓2012-09-18 17:05:23
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 vigaryang 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機(jī)用戶 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4 | ioodiiij 2026-03-17 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
|
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研]
 11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|