| 5 | 1/1 | 返回列表 |
| 查看: 2534 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
lulutiantian新蟲 (小有名氣)
|
[求助]
Au 自由原子計(jì)算的結(jié)果,大家?guī)兔聪聦?duì)不對(duì)?
|
||
|
因?yàn)橐?jì)算一個(gè)體系的結(jié)合能,so計(jì)算一個(gè)金原子在自由態(tài)的能量,建立一個(gè)15*15*15A的空格子,將一個(gè)Au原子放在中心。未做截?cái)嗄軠y(cè)試,選取截?cái)嗄躤cut為360ev,進(jìn)行結(jié)構(gòu)優(yōu)化,然后靜態(tài)計(jì)算,k點(diǎn)設(shè)為3*3*3,以Gamma點(diǎn)為中心產(chǎn)生,靜態(tài)計(jì)算后結(jié)果 E=-0.10956204 而候博士的入門指南中以Al為例說單原子能量在1mev-10mev之間。也就是vasp計(jì)算的單個(gè)原子能量與原子參考組態(tài)的能量差。在vasp所計(jì)算的總能是扣除了原子參考組態(tài)的能量,即POTCAR中EATOM的值。 問題 1.做自由原子能量計(jì)算不測(cè)試截?cái)嗄芸梢詥幔?br /> 2.我的這個(gè)E=-0.10956204很明顯不在1mev-10mev之間。這個(gè)值可以作為Au自由原子的能量嗎? 3.我在計(jì)算中的發(fā)現(xiàn) Au 自由原子優(yōu)化結(jié)算結(jié)果 E= -0.116817 Au自由原子靜態(tài)計(jì)算結(jié)果E=-0.10956204 發(fā)現(xiàn)這兩個(gè)值是有差別的。那么應(yīng)該是以靜態(tài)計(jì)算的結(jié)果為結(jié)合能。 希望做過這方面計(jì)算的蟲友賜教! |
銅蟲 (小有名氣)
|
本人用的不是vasp,所以針對(duì)一般情況隨便說說,不對(duì)的話樓主指正: 1)求結(jié)合能應(yīng)該分別得出Au原子的能量和Au晶體的能量,然后做差。 2)選定贗勢(shì)以后,對(duì)于Au原子和Au晶體,至少有一個(gè)必須測(cè)定截?cái)嗄,截(cái)嗄鼙厝徊荒茈S便亂選的,最好自己測(cè)測(cè),或者有前人的經(jīng)驗(yàn)也行。Au原子和Au晶體可以使用同樣大小的截?cái)嗄堋?br /> 3)K點(diǎn)跟正格子的大小有關(guān)系,對(duì)于Au原子來說,需要使用一個(gè)大格子模擬,所以一個(gè)K點(diǎn)基本是夠的,Au晶體格子小,需要密集一些的K點(diǎn)。 4)Au原子直接放在格子什么地方都行,一個(gè)格子里一個(gè)Au原子,不需要做結(jié)構(gòu)優(yōu)化,但格子應(yīng)該取足夠大,以減少在周期性條件下各個(gè)格子之間的相互作用。Au晶體的話,對(duì)于一般的fcc的Au來說,需要求出最低能量的晶格常數(shù),可以做varible-cell relaxation,或者直接看各個(gè)晶格常數(shù)下Au晶體的能量,然后選取對(duì)應(yīng)最低能量的晶格常數(shù)。 |
金蟲 (小有名氣)

新蟲 (小有名氣)
金蟲 (正式寫手)
納米寶貝
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +7 | zjl的號(hào) 2026-03-16 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 070300化學(xué)319求調(diào)劑 +4 | 錦鯉0909 2026-03-17 | 4/200 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
[考研] 材料與化工求調(diào)劑 +5 | 為學(xué)666 2026-03-16 | 5/250 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +3 | 困于星晨 2026-03-17 | 3/150 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 285化工學(xué)碩求調(diào)劑(081700) +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +7 | Ncdx123456 2026-03-13 | 8/400 |
|
|
[考研] 341求調(diào)劑 +4 | 番茄頭--- 2026-03-10 | 4/200 |
|
|
[考研] 材料專碩288分求調(diào)劑 一志愿211 +4 | 在家想你 2026-03-11 | 4/200 |
|
|
[考研] 一志愿中科院,化學(xué)方向,295求調(diào)劑 +4 | 一氧二氮 2026-03-11 | 4/200 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|