| 7 | 1/1 | 返回列表 |
| 查看: 5333 | 回復(fù): 6 | |||||
貝芷苡銅蟲 (小有名氣)
|
[求助]
摻雜化學(xué)勢(shì)及形成能計(jì)算
|
|
本人在算氧化鋅中摻雜形成能,但有些計(jì)算參數(shù)找不到參照,特發(fā)此貼,希望各位高手給指導(dǎo)一下。一些 氧化鋅形成焓:Hf ZnO=(1/M)EtotZnO-(1/N)EtotZn-1/2EtotO2, where Hf ZnO is the enthalpy of formation。EtotZnO為含M個(gè)鋅原子的超胞總能,EtotZn為N個(gè)鋅原子的超胞總能,超胞均為六方密排結(jié)構(gòu),EtotO2為邊長(zhǎng)為10埃的立方體中放一氧分子的總能,此時(shí)得出的形成焓滿足Uzn+Uo=Hf ZnO,本人計(jì)算各參數(shù)值為Hf ZnO=-3.443eV此值與文獻(xiàn)中的值符合,(1/N)EtotZn=-1.262eV,1/2EtotO2=-4.388eV,那么鋅空位在富氧Uo=0,Uzn=Hf ZnO=-3.443eV及缺氧條件Uzn=0,Uo=Hf ZnO=-3.443eV下的形成能(計(jì)算公式:Ef=Etot(q)-Etot(p)+Uzn,Etot(q) 一下為我的個(gè)輸入文件 INCAR SYSTEM =ZnO ISTART =0 ICHARG =2 ENCUT =380 EDIFF =1E-5 EDIFFG =-0.01 ISMEAR=0 SIGMA=0.1 NSW=100 NGX=62 NGY=62 NGZ=66 ISPIN=2 IBRION = 2 ISIF =3 PREC=Accurte LREAL= Auto KPOINTS SYSTEM =ZnO 0 M 3 3 2 0.0 0.0 0.0 靜態(tài)計(jì)算時(shí)KPOINTS SYSTEM =ZnO 0 M 5 5 3 0.0 0.0 0.0 另EtotZn為N個(gè)鋅原子的超胞總能,超胞均為六方密排結(jié)構(gòu),EtotZn=-1.262(晶胞含2個(gè)鋅原子),EtotZn=-0.944,哪個(gè)更準(zhǔn)確? |
帶電摻雜體系形成能 | vasp計(jì)算資料 | @相關(guān) |

木蟲 (小有名氣)
|
你的q值、超胞大小都和文獻(xiàn)一樣嗎?尤其是超胞大小對(duì)結(jié)果影響較大。PRB, 61, 15019 公式也少了一大項(xiàng),一般情況下Eformation都和費(fèi)米能級(jí)、價(jià)帶頂能量以及矯正位勢(shì)有關(guān)。PRB, 54, 16676 |

木蟲 (小有名氣)

銅蟲 (小有名氣)
|
UGa+UN=UGaN(bulk),中的UGaN(bulk)就是生成焓吧? 在PRB 63,075205這篇文章中有這樣的介紹Under thermal equilibrium UZn+UO=DHZnO , where DHZnO is the formation enthalpy of ZnO,在Journal of Crystal Growth 287 (2006) 58–65中則是這樣描述的The upper bounds defined above also lead to lower bounds by using the thermodynamic stability condition for ZnO,UZn+UO = EtotZnO, (3)where EtotZnO is the total energy of a two-atom unit of bulk ZnO. |

木蟲 (小有名氣)

|
化學(xué)勢(shì)和形成焓在概念上完全不一樣。對(duì)于焓值,當(dāng)不考慮壓力對(duì)體積的影響時(shí),零點(diǎn)內(nèi)能就可以代替焓值,如果你算的不是零點(diǎn)的焓值,這個(gè)時(shí)候就要考慮溫度增加的同時(shí)振動(dòng)能的影響,對(duì)于振動(dòng)能,需要計(jì)算聲子譜及其積分,還需考慮簡(jiǎn)諧近似或者準(zhǔn)簡(jiǎn)諧近似的結(jié)果。對(duì)于各種焓值的定義,也推薦你看看Phys.Rev.B 73.144104(2006)這篇文章。對(duì)于摻雜形成能你可以上網(wǎng)搜一下計(jì)算H分子形成能的腳本文章,里面有一篇PRB的參考文獻(xiàn),好像也是VASP課題組發(fā)的文章,可能對(duì)你計(jì)算能量有所幫助。 對(duì)于化學(xué)勢(shì),是自由能對(duì)原子摩爾數(shù)的微分。在自由能里面已經(jīng)不僅僅是焓的影響,還有熵項(xiàng)的影響,這些熵項(xiàng)來自于組態(tài)、振動(dòng)、電子等方面。在計(jì)算的時(shí)候需要考慮。通常情況下只能通過自由能曲線斜率或者其他方式獲得化學(xué)勢(shì)。 |

| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +9 | 小熊joy 2026-03-14 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +3 | 睿08 2026-03-17 | 5/250 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +6 | 妮妮ninicgb 2026-03-15 | 9/450 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[考研]
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級(jí)475 +5 | Becho359 2026-03-11 | 5/250 |
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] 336求調(diào)劑 +6 | Iuruoh 2026-03-11 | 6/300 |
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|