| 5 | 1/1 | 返回列表 |
| 查看: 2124 | 回復(fù): 7 | |||
| 【懸賞金幣】回答本帖問(wèn)題,作者北科大小學(xué)渣將贈(zèng)送您 500 個(gè)金幣 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
北科大小學(xué)渣新蟲 (正式寫手)
|
[求助]
VASP計(jì)算石墨烯-Ni(111)界面,結(jié)構(gòu)跑散 已有2人參與
|
||
|
各位蟲友好,最近我在用VASP計(jì)算石墨烯-Ni(111)界面的結(jié)合能和相互作用。我用了4層Ni金屬層和單層石墨烯,根據(jù)多個(gè)文獻(xiàn)中的報(bào)道,Ni原子排列在C原子上方和六元環(huán)中間的上方最穩(wěn)定,我按這個(gè)構(gòu)建了超晶胞,并設(shè)置了20埃的真空層,這種排列Ni(111)晶格失配度在1.2%。結(jié)構(gòu)優(yōu)化計(jì)算中,我固定了石墨烯層中所有的碳原子,并且固定了晶胞參數(shù),只允許所有金屬原子自由弛豫。但是弛豫計(jì)算結(jié)束后,我期待的是所有金屬原子應(yīng)該在x-y面上的位置幾乎不變,應(yīng)該在z軸上弛豫,并且每層金屬層上的Ni原子應(yīng)該還會(huì)保持在同一高度(z坐標(biāo)一樣),這樣才會(huì)和文獻(xiàn)中報(bào)道的一樣,才能得到穩(wěn)定的界面距離和界面結(jié)構(gòu)。 但是,當(dāng)結(jié)束結(jié)構(gòu)弛豫計(jì)算后,結(jié)構(gòu)似乎“跑散”了,即越是遠(yuǎn)離石墨烯的金屬層上的Ni原子不僅會(huì)在x-y面上弛豫,并且也會(huì)在z方向上弛豫,導(dǎo)致弛豫后Ni原子不在C原子的正上方或六元環(huán)的中心的上方,并且z方向上,同一金屬層的Ni原子也不再同一高度了,越遠(yuǎn)離石墨烯的金屬層,這種現(xiàn)象更明顯。這是為什么?怎么才能得到文獻(xiàn)中一樣的效果?我的計(jì)算問(wèn)題出在哪兒? 感謝各位提供寶貴的意見。 |
木蟲 (正式寫手)
|
你的Ni(111)層在遠(yuǎn)離石墨烯時(shí)“跑散”可能是因?yàn)椋?br />
Ni 原子自由度過(guò)高,應(yīng)固定底部2層Ni,或在x-y方向約束部分Ni原子。VASP 的 Selective Dynamics 選項(xiàng)可以用于固定某些原子的 x-y 方向,保證其只沿 z 方向弛豫。例如,在 POSCAR 里設(shè)置: Selective dynamics Direct x1 y1 z1 F F T x2 y2 z2 F F T ... 弛豫參數(shù)(IBRION, POTIM, EDIFFG)需要優(yōu)化,嘗試 IBRION = 2、POTIM = 0.1–0.15、EDIFFG = -0.01。 檢查K點(diǎn)密度和贗勢(shì)選擇,以確保精度足夠。試試 EDIFFG = -0.01 或更嚴(yán)格的值,看看是否能得到更合理的結(jié)構(gòu)。 |
新蟲 (正式寫手)
新蟲 (正式寫手)
|
感謝你的回答。我是只固定了石墨烯層,允許所有金屬層在三個(gè)方向上自由弛豫,結(jié)果金屬層上的Ni就會(huì)“跑散”,因?yàn)槲铱吹轿墨I(xiàn)中有的文獻(xiàn)報(bào)道的是固定單層石墨烯和晶胞參數(shù),允許所有金屬原子弛豫;還有的文獻(xiàn)是使用3+3金屬層slab,固定遠(yuǎn)離石墨烯的三層金屬層,讓靠近石墨烯的三層金屬和石墨烯自由弛豫。如果文獻(xiàn)中報(bào)道的正確的話,這兩種方法都不應(yīng)該使結(jié)構(gòu)“跑散”。你說(shuō)的固定x, y方向,只允許部分Ni原子在z方向弛豫,會(huì)不會(huì)導(dǎo)致得不到最穩(wěn)定的構(gòu)型?我理解的是即便不固定金屬Ni原子,它們也應(yīng)該大概率只會(huì)沿著z方向弛豫,因?yàn)槲掖罱ǖ臉?gòu)型,Ni原子的排列方式應(yīng)該是最穩(wěn)定的構(gòu)型,所有我期待著不需要多少離子步就該結(jié)束計(jì)算,因?yàn)樵冀Y(jié)構(gòu)已經(jīng)很接近最穩(wěn)定的構(gòu)型,只是z方向上的距離需要優(yōu)化。 |
金蟲 (正式寫手)
|
First, check if you have reproduced all the lattice constants reported in the former papers Second, if it is "yes" for the first question, try using the exactly same methods described in the former papers. In addition, you can consider constraining the x- and y-coordinates of Ni in the coarse calculations of relaxation, and in the more accurate relaxations, remove the constraints to let the Ni atoms be fully relaxed. |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 26超級(jí)電容器申博 +4 | dhdjdjend 2026-02-25 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 中國(guó)科學(xué)技術(shù)大學(xué)材料與化工281求調(diào)劑,有科研和獲獎(jiǎng)經(jīng)歷 +9 | wsxw 2026-03-02 | 10/500 |
|
|
[考研]
材料325求調(diào)劑
30+5
|
mariusuki 2026-03-02 | 10/500 |
|
|
[考研] 289求調(diào)劑 +8 | yang婷 2026-03-02 | 10/500 |
|
|
[考研]
|
15779376950 2026-03-01 | 6/300 |
|
|
[考研] 一志愿314求調(diào)劑 +4 | 202111120625 2026-03-03 | 4/200 |
|
|
[考研] 085600材料與化工 298 調(diào)劑 +3 | 小西笑嘻嘻 2026-03-03 | 3/150 |
|
|
[考研] 一志愿中科大能動(dòng)297求調(diào)劑,本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
[考研] 一志愿天津大學(xué)085600 319分 材料與化工 金屬方向 求調(diào)劑 +4 | 青科11 2026-03-02 | 4/200 |
|
|
[考研] 求調(diào)劑 +4 | Guo_yuxuan 2026-03-02 | 5/250 |
|
|
[考研] 環(huán)境工程專碩307 一志愿211 四六級(jí)已過(guò) 求調(diào)劑 +3 | ccc! 2026-03-03 | 3/150 |
|
|
[考研]
材料工程專碩283求調(diào)劑
5+8
|
,!? 2026-03-02 | 10/500 |
|
|
[考研] 化工335求調(diào)劑 +5 | 摸摸貓貓頭 2026-03-02 | 5/250 |
|
|
[考研] 298求調(diào)劑 +7 | axyz3 2026-02-28 | 8/400 |
|
|
[考研] 261求調(diào)劑 +3 | 陸lh 2026-03-01 | 3/150 |
|
|
[考博] 誠(chéng)招農(nóng)業(yè)博士 +3 | 心欣向榮 2026-02-28 | 3/150 |
|
|
[考研] 284求調(diào)劑 +10 | 天下熯 2026-02-28 | 11/550 |
|
|
[考研] 279求調(diào)劑 +3 | dua1 2026-03-01 | 4/200 |
|
|
[考研] 0856材料求調(diào)劑 +4 | 麻辣魷魚 2026-02-28 | 4/200 |
|
|
[考研] 313求調(diào)劑 +3 | 水流年lc 2026-02-28 | 3/150 |
|