| 24小時(shí)熱門版塊排行榜 |
| 查看: 5788 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 306求調(diào)劑 +6 | chuanzhu川燭 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 306求調(diào)劑 +5 | 來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-22 | 5/250 |
|
|
[考研] 305分求調(diào)劑(食品工程) +4 | Sxy112 2026-03-21 | 6/300 |
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
[考研] 269專碩求調(diào)劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 3/150 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 321求調(diào)劑 +9 | 何潤(rùn)采123 2026-03-18 | 11/550 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開(kāi)心 2026-03-16 | 4/200 |
|