| 24小時(shí)熱門(mén)版塊排行榜 |
| 查看: 5820 | 回復(fù): 15 | ||||
Illusionist銀蟲(chóng) (正式寫(xiě)手)
|
[求助]
如何做vasp摻雜
|
|||
|
本人想用vasp做摻雜計(jì)算,剛開(kāi)始學(xué) 1.直接建立3x3x3的超晶格模型 1)關(guān)于K的選擇 3x3x3 -1019.397709 4x4x4 -1019.397873 所以選擇了3x3x3 2)ENCUT 280 4093.144030 300 4093.144030 320 4093.144030 340 4093.144030 這里我最后選擇了500ENCUT去進(jìn)行結(jié)構(gòu)優(yōu)化。 2.結(jié)構(gòu)優(yōu)化INCAR: SYSTEM=SnO2 ENCUT=500 ISMEAR=-5 GGA=91 VOSKOWN=1 ISYM=0 ISIF=3 NSW=100 PTIOM=0.1 IBRION=2 這里我選擇打破對(duì)稱性不知道是否合理 ISYM=0設(shè)置前后影響挺大,起碼晶格常數(shù)會(huì)變化。LREAL=auto設(shè)置后,晶格常數(shù)也會(huì)有明顯變化,請(qǐng)問(wèn)一下一般如何處理這樣的問(wèn)題 另外我想在這一步結(jié)束后摻入替位的Sb,但是Sb有+3和+5價(jià)之分,該如何設(shè)定價(jià)態(tài)呢 |
VASP | 科研 |
鐵桿木蟲(chóng) (著名寫(xiě)手)

鐵桿木蟲(chóng) (著名寫(xiě)手)
|
首先對(duì)你的數(shù)據(jù): 2)ENCUT 280 4093.144030 300 4093.144030 320 4093.144030 340 4093.144030 表示懷疑 其次對(duì)稱性沒(méi)必要打破,設(shè)置isif=3就可以了 |

銀蟲(chóng) (小有名氣)
|
先在1x1x1體系下用ISIF=7優(yōu)化得到粗略的晶格常數(shù) 然后做ENCUT和KPOINTS的收斂性測(cè)試 按滿足收斂條件的ENCUT和KPOINTS值,計(jì)算不同晶格常數(shù)下體系的能量 用murn.f擬合得到平衡態(tài)晶格常數(shù)。 把1x1x1體系按比例放大,得到2x2x2/3x3x3/4x4x4等體系。 對(duì)3x3x3體系: KPOINTS設(shè)成 然后進(jìn)行ENCUT收斂性測(cè)試,記錄不同ENCUT情況下能量和cpu time 然后在4x4x4x體系下做同樣的事情,記錄不同ENCUT情況下能量和cpu time 再重復(fù)至5x5x5,6x6x6,7x7x7 計(jì)算時(shí)間超過(guò)48小時(shí)的任務(wù)就可以忽略了。 綜合考慮性能和精度的影響,確定最后選用的摻雜體系,例如我就選用4x4x4體系。 在此體系下,加入摻雜原子,按照之前收斂性測(cè)試的結(jié)果設(shè)置ENCUT 接下來(lái)就開(kāi)始幾何優(yōu)化了。這里一般要求保持原胞的形狀和體積,不過(guò)你也可以先優(yōu)化原胞的體積,一般不要優(yōu)化原胞的形狀。不過(guò)這也說(shuō)不定,例如極高壓情況中體系的晶格會(huì)改變(例如從體心立方變成六角密堆,所以原胞形狀是否改變跟你所做的具體工作有關(guān)。) 優(yōu)化好原胞以后再做離子遲豫。你可以自己寫(xiě)腳本計(jì)算,也可以讓vasp自己算。 |
木蟲(chóng) (正式寫(xiě)手)
|
首先打破對(duì)稱性是可以的,只不過(guò)計(jì)算可能會(huì)多費(fèi)點(diǎn)時(shí)間; 至于替位摻雜,要考慮一下?lián)诫s原子與被摻雜原子是否是等價(jià)的,若果是非等價(jià)摻雜要考慮電荷補(bǔ)償?shù)膯?wèn)題,相關(guān)參數(shù)是NELECT,具體關(guān)于這方面之前本版有很多的相關(guān)帖子可以參考,自己搜一下。 |
銀蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
銀蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)

銀蟲(chóng) (小有名氣)
|
摻雜的問(wèn)題,建議把K格點(diǎn)設(shè)成只有Gamma點(diǎn),再用更大的晶格(例如4x4x4,5x5x5,6x6x6等)。 至于收斂性測(cè)試,要注意收斂條件是 ,而不是 在優(yōu)化幾何結(jié)構(gòu)的時(shí)候,建議不要使用ISIF=3選項(xiàng)!速度會(huì)非常慢。在同樣的計(jì)算條件下,建議用loop.sh方案。 一般來(lái)說(shuō),先在1x1x1體系下在實(shí)驗(yàn)測(cè)得的晶格常數(shù)基礎(chǔ)上用ISIF=7優(yōu)化得到體系的晶格常數(shù)。 然后作等比例放大,得到nxnxn體系的POSCAR,然后用循環(huán)靜態(tài)計(jì)算這個(gè)晶格常數(shù)附近幾個(gè)值對(duì)應(yīng)的體系的能量,用murn方程擬合得到體系的平衡態(tài)晶格常數(shù)。 |
銀蟲(chóng) (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300一志愿211,312分求調(diào)劑院校 +6 | 小黃鴨寶 2026-03-30 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 085601 329分調(diào)劑 +4 | yzsa12 2026-03-31 | 4/200 |
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +10 | 深V宿舍吧 2026-03-29 | 10/500 |
|
|
[考研] 抱歉 +4 | 田洪有 2026-03-30 | 4/200 |
|
|
[考研] 一志愿中海洋材料357 +3 | 麥恩莉. 2026-03-30 | 3/150 |
|
|
[考研] 085600 295分求調(diào)劑 +7 | W55j 2026-03-30 | 9/450 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +9 | hanamiko 2026-03-27 | 9/450 |
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
[考研] 297求調(diào)劑 +17 | 田洪有 2026-03-26 | 18/900 |
|
|
[考研] 328求調(diào)劑 +8 | 嗯滴的基本都 2026-03-27 | 8/400 |
|
|
[考研] 342求調(diào)劑 +4 | 加油a李zs 2026-03-26 | 4/200 |
|
|
[考研] 298求調(diào)劑 +3 | 什么是胖頭魚(yú) 2026-03-30 | 5/250 |
|
|
[考研] 本科211生物醫(yī)學(xué)工程085409求調(diào)劑339分 +3 | 里子木yy 2026-03-29 | 3/150 |
|
|
[考研] 282求調(diào)劑 +4 | wcq131415 2026-03-24 | 4/200 |
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|