| 5 | 1/1 | 返回列表 |
| 查看: 2123 | 回復: 7 | |||
| 【懸賞金幣】回答本帖問題,作者北科大小學渣將贈送您 500 個金幣 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[求助]
VASP計算石墨烯-Ni(111)界面,結(jié)構(gòu)跑散 已有2人參與
|
|||
|
各位蟲友好,最近我在用VASP計算石墨烯-Ni(111)界面的結(jié)合能和相互作用。我用了4層Ni金屬層和單層石墨烯,根據(jù)多個文獻中的報道,Ni原子排列在C原子上方和六元環(huán)中間的上方最穩(wěn)定,我按這個構(gòu)建了超晶胞,并設(shè)置了20埃的真空層,這種排列Ni(111)晶格失配度在1.2%。結(jié)構(gòu)優(yōu)化計算中,我固定了石墨烯層中所有的碳原子,并且固定了晶胞參數(shù),只允許所有金屬原子自由弛豫。但是弛豫計算結(jié)束后,我期待的是所有金屬原子應該在x-y面上的位置幾乎不變,應該在z軸上弛豫,并且每層金屬層上的Ni原子應該還會保持在同一高度(z坐標一樣),這樣才會和文獻中報道的一樣,才能得到穩(wěn)定的界面距離和界面結(jié)構(gòu)。 但是,當結(jié)束結(jié)構(gòu)弛豫計算后,結(jié)構(gòu)似乎“跑散”了,即越是遠離石墨烯的金屬層上的Ni原子不僅會在x-y面上弛豫,并且也會在z方向上弛豫,導致弛豫后Ni原子不在C原子的正上方或六元環(huán)的中心的上方,并且z方向上,同一金屬層的Ni原子也不再同一高度了,越遠離石墨烯的金屬層,這種現(xiàn)象更明顯。這是為什么?怎么才能得到文獻中一樣的效果?我的計算問題出在哪兒? 感謝各位提供寶貴的意見。 |
|
感謝你的回答。我是只固定了石墨烯層,允許所有金屬層在三個方向上自由弛豫,結(jié)果金屬層上的Ni就會“跑散”,因為我看到文獻中有的文獻報道的是固定單層石墨烯和晶胞參數(shù),允許所有金屬原子弛豫;還有的文獻是使用3+3金屬層slab,固定遠離石墨烯的三層金屬層,讓靠近石墨烯的三層金屬和石墨烯自由弛豫。如果文獻中報道的正確的話,這兩種方法都不應該使結(jié)構(gòu)“跑散”。你說的固定x, y方向,只允許部分Ni原子在z方向弛豫,會不會導致得不到最穩(wěn)定的構(gòu)型?我理解的是即便不固定金屬Ni原子,它們也應該大概率只會沿著z方向弛豫,因為我搭建的構(gòu)型,Ni原子的排列方式應該是最穩(wěn)定的構(gòu)型,所有我期待著不需要多少離子步就該結(jié)束計算,因為原始結(jié)構(gòu)已經(jīng)很接近最穩(wěn)定的構(gòu)型,只是z方向上的距離需要優(yōu)化。 |
木蟲 (正式寫手)
|
你的Ni(111)層在遠離石墨烯時“跑散”可能是因為: Ni 原子自由度過高,應固定底部2層Ni,或在x-y方向約束部分Ni原子。VASP 的 Selective Dynamics 選項可以用于固定某些原子的 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點密度和贗勢選擇,以確保精度足夠。試試 EDIFFG = -0.01 或更嚴格的值,看看是否能得到更合理的結(jié)構(gò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ā)表 | |
|---|---|---|---|---|
|
[考研] 0856材料與化工,270求調(diào)劑 +12 | YXCT 2026-03-01 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|
|
[考博] 26超級電容器申博 +4 | dhdjdjend 2026-02-25 | 4/200 |
|
|
[考研] 環(huán)境工程學碩288求助調(diào)劑 +7 | 多吃億口芝士 2026-03-02 | 7/350 |
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
[考博] 26申博 求博導 +3 | 愛讀書的小帥 2026-02-28 | 5/250 |
|
|
[考研] 0703化學306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|
|
[基金申請] 請問大家,研究風險與應對措施那里, 大家都怎么寫呢 ? +3 | cauasen 2026-03-02 | 3/150 |
|
|
[考研] 理學,工學,農(nóng)學調(diào)劑,少走彎路,這里歡迎您! +8 | likeihood 2026-03-02 | 11/550 |
|
|
[考研] 環(huán)境工程專碩307 一志愿211 四六級已過 求調(diào)劑 +3 | ccc! 2026-03-03 | 3/150 |
|
|
[考研] 化學0703求調(diào)劑 學碩 理/工科均可 總分279 +3 | 1一11 2026-03-03 | 5/250 |
|
|
[考研] 考研復試調(diào)劑,過國家線的同學都可報名 +7 | 黑!在干嘛 2026-02-28 | 8/400 |
|
|
[考研] 306分材料調(diào)劑 +5 | chuanzhu川燭 2026-03-01 | 6/300 |
|
|
[考研] 材料調(diào)劑 +3 | 恒順自然 2026-03-02 | 3/150 |
|
|
[考研] 272求調(diào)劑 +6 | 田智友 2026-02-28 | 6/300 |
|
|
[考研] 化工299分求調(diào)劑 一志愿985落榜 +5 | 嘻嘻(*^ω^*) 2026-03-01 | 5/250 |
|
|
[考博] 26申博 +4 | 想申博! 2026-02-26 | 6/300 |
|
|
[考研] 328求調(diào)劑 +3 | aaadim 2026-03-01 | 5/250 |
|
|
[考研] 290求調(diào)劑 +9 | 材料專碩調(diào)劑; 2026-02-28 | 11/550 |
|
|
[考研] 317一志愿華南理工電氣工程求調(diào)劑 +6 | Soliloquy_Q 2026-02-28 | 11/550 |
|