| 10 | 1/1 | 返回列表 |
| 查看: 2391 | 回復(fù): 9 | ||||
[交流]
【求助】對于難收斂的體系 已有6人參與
|
| 各位,對于vasp優(yōu)化難收斂的體系,可不可以這樣做,先用低的收斂標(biāo)準(zhǔn)使其收斂,然后讀入波函數(shù)和電荷再采用高的收斂標(biāo)準(zhǔn)進(jìn)行計算,可行否?多謝指教! |
VASP優(yōu)化不收斂 |
捐助貴賓 (正式寫手)
木蟲 (正式寫手)
Dr.

榮譽(yù)版主 (知名作家)
老和山猥瑣派九段
木蟲 (正式寫手)
Dr.
|
如果已經(jīng)開始弛豫了好幾步了,就copy過來,如果第一步都沒過去,那就不必了。 收斂參數(shù)的調(diào)整量子化學(xué)論壇以前有個老帖 你自己看看有沒有你能用的 對于一些磁性體系、鑭系和錒系元素及相關(guān)化合物的靜態(tài)計算(電子迭代),經(jīng)常會遇到“難收斂”的問題。 下面給出幾個相關(guān)Flag及設(shè)置方法:(有些話從說明書上摘錄的,不大好翻譯,就copy了,請見諒) 1、LMAXMIX Default: LMAXMIX = 2 An additional flag controls up to which l quantum number the onsite PAW charge densities are passed through the charge density mixer. Higher l-quantum numbers are usually not handled by the mixer. In order to obtain fast convergence to the groundstate, you can try the following setting: LMAXMIX = 4 for d elements LMAXMIX = 6 for f elements 這個FLAG對于含d電子和f電子的體系是非常重要的,很大一部分體系的收斂問題可以通過設(shè)置合適的LMAXMIX值來解決。 2、ALGO, IALGO, LDIAG If the self-consistency loop does not converge within 40 steps, it will probably not converge at all. In this case you should reconsider the tags IALGO, LDIAG, and the mixing-parameters. 這是說明書上的建議。 一般情況下,或使用IALGO=48時遇到收斂問題的話,可以考慮設(shè)IALGO為38(4.5以前的版本可設(shè)為8),或設(shè)置ALGO=Normal or Fast (in VASP.4.5 and later versions)。 3、NELMDL NELMDL gives the number of non-selfconsistent steps at the beginning; if one initializes the wave functions randomly the initial wave functions are far from anything reasonable. The resulting charge density is also 'nonsense'. Therefore it makes sense to keep the initial Hamiltonian, which corresponds to the superposition of atomic charge densities, fixed during the first few steps. Choosing a 'delay' for starting the charge density update becomes essential in all cases where the SC-convergence is very bad (e.g. surfaces or molecules/clusters chains). Without setting a delay VASP will probably not converge or at least the convergence speed is slowed down. NELMDL might be positive or negative. A positive number means that a delay is applied after each ionic movement -- in general not a convenient option. A negative value results in a delay only for the start-configuration. 4、mixing-parameters 對于一些難收斂的體系,可以使用“l(fā)inear mixing”,具體詳見VASP說明書中的“Mixing-tags”。 For an initial linear mixing (BMIX ~ 0) an optimal setting for A(AMIX) can be found easily by setting Aopt=Acurrent*Γmean. For the Kerker scheme either A or q0(i.e. AMIX or BMIX) can be optimized, but we recommend to change only BMIX and keep AMIX fixed (you must decrease BMIX if the mean eigenvalue is larger than one, and increase BMIX if the mean eigenvalue is smaller than one). 盡管VASP說明書中給出了調(diào)節(jié)AMIX和BMIX的一些較為明確的建議,但是實(shí)際去調(diào)節(jié)的時候,還是挺難的,但原則上說,是可以通過調(diào)節(jié)這兩個Flag來使得收斂問題得以解決的,只是得有耐心。 5、kmesh, SIGMA 收斂問題還跟kmesh及SIGMA(當(dāng)使用ISMEAR不等于-5 和-4時)的設(shè)置有關(guān)。要達(dá)到同樣的精度,較小的SIGMA則需要較大的kmesh;而且,當(dāng)SIGMA較小時,若kpoints不夠多,也會出現(xiàn)難收斂的情況。 P.S. 對于不同系統(tǒng)的計算,問題的原因不一定一樣,因而可能解決之一問題的方法也不一定會一樣。其實(shí)以上大部分都來自VASP的說明書,只是我在遇到收斂問題時試過這些方法. |

榮譽(yù)版主 (知名作家)
…………

金蟲 (著名寫手)
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 材料專碩306英一數(shù)二 +8 | z1z2z3879 2026-03-16 | 10/500 |
|
|
[考研] 302求調(diào)劑 +6 | 呼呼呼。。。。 2026-03-17 | 6/300 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +4 | 錦鯉0909 2026-03-17 | 4/200 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考研] 0854可跨調(diào)劑,一作一項核心論文五項專利,省、國級證書40+數(shù)一英一287 +3 | 小李0854 2026-03-16 | 3/150 |
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
[碩博家園] 深圳大學(xué)碩士招生(2026秋,傳感器方向,僅錄取第一志愿) +4 | xujiaoszu 2026-03-11 | 9/450 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[基金申請]
有必要更換申報口嗎
20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 快樂的香蕉 2026-03-11 | 4/200 |
|
|
[考研] 308求調(diào)劑 +5 | 是Lupa啊 2026-03-11 | 5/250 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[考研] 290求調(diào)劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
[考研] 工科調(diào)劑 +4 | Jiang191123! 2026-03-11 | 4/200 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|