 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 6183 | 回復(fù): 30 | ||||||||
[求助]
vasp表面計(jì)算詳細(xì)解說 已有1人參與
|
||||||||
|
想和大家交流下vasp計(jì)算表面的一些問題： 1. 計(jì)算表面的步驟： a。先建立一個(gè)表面超胞模型(可以加入對(duì)稱性嗎？因?yàn)槿绻雽?duì)稱性，就變成超胞的兩邊是真空層了)； b。然后測(cè)試真空層和原子層數(shù)，這里的收斂標(biāo)準(zhǔn)是能量小于10eV（測(cè)試原子層數(shù)的時(shí)候可以不對(duì)稱嗎？）； c。然后進(jìn)行表面結(jié)構(gòu)的弛豫，固定中間幾層原子，觀察層間距的變化，若固定的層數(shù)之間有較大的位移（原子既然固定了，固定原子層之間還會(huì)有有位移嘛？），就增加弛豫原子層數(shù)（是固定原子的層數(shù)還是沒固定原子的層數(shù)），在觀察層間距的變化，知道變化很小為止（這里是做靜態(tài)計(jì)算嗎？然后從OUTcar里讀取原子位置，看層之間的變化？這里有一個(gè)什么標(biāo)準(zhǔn)嗎？）； d。此外，還要測(cè)試表面能，這時(shí)要保證slab足夠厚，當(dāng)slab層數(shù)變大時(shí)，表面能變化很小就行了（這里表面能使需要用公式計(jì)算的吧~~~，這個(gè)和原子層數(shù)測(cè)試是同時(shí)進(jìn)行的嗎？邊看總能收不收斂，還要看表面能收不收連？而且說要保證晶胞大小一致，那這樣真空層不就是變了嗎，那前面測(cè)試真空層不就沒有意義了）。 e。計(jì)算表面的性質(zhì) 括號(hào)里面的就是我的疑問,還有一個(gè)問題就是k點(diǎn)的測(cè)試在什么時(shí)候測(cè)，k點(diǎn)在垂直晶面的方向上為1，那么如果是110面，k點(diǎn)該怎么設(shè)呢？ 我也是在學(xué)習(xí)的起步階段，最近看這些有些糊涂，思路不太清晰，好像小木蟲上也沒有詳細(xì)的帖子，希望我們大家交流一下， |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 好書！好讀！ | 仿真建模與計(jì)算 | 第一性原理相關(guān)文檔 | 第一性原理計(jì)算-tg |
| 第一性原理和電化學(xué) | 表界面計(jì)算 |
» 本帖已獲得的紅花（最新10朵）
» 猜你喜歡
278求調(diào)劑
已經(jīng)有7人回復(fù)
-
26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑
已經(jīng)有5人回復(fù)
-
280求調(diào)劑
已經(jīng)有9人回復(fù)
-
279求調(diào)劑
已經(jīng)有5人回復(fù)
-
286求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工（0856）304求 B區(qū) 調(diào)劑
已經(jīng)有3人回復(fù)
-
268求調(diào)劑
已經(jīng)有8人回復(fù)
-
265求調(diào)劑
已經(jīng)有14人回復(fù)
-
22 350 本科985求調(diào)劑，求老登收留
已經(jīng)有3人回復(fù)
-
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 關(guān)于（VASP）Fe100表面的SCF收斂
已經(jīng)有19人回復(fù)
- 在用vasp進(jìn)行表面吸附時(shí)，遇到的問題，求高手解惑！
已經(jīng)有12人回復(fù)
- VASP計(jì)算H原子 吸附在Pt(111)表面時(shí)，最后CONTCAR 里少了H 原子
已經(jīng)有3人回復(fù)
- VASP計(jì)算TiO2表面能問題 求助
已經(jīng)有8人回復(fù)
- vasp計(jì)算金屬表面的吸附能問題
已經(jīng)有8人回復(fù)
- vasp表面和吸附分子的結(jié)構(gòu)分別優(yōu)化后，如何將它們放到一起
已經(jīng)有11人回復(fù)
- VASP計(jì)算表面吸附 但是能量總是不收斂/。。。
已經(jīng)有9人回復(fù)
- vasp表面計(jì)算的結(jié)構(gòu)優(yōu)化問題
已經(jīng)有7人回復(fù)
- 關(guān)于vasp計(jì)算表面能中的Ebulk
已經(jīng)有7人回復(fù)
- 用VASP計(jì)算表面吸附問題時(shí)，考慮范德華力時(shí)的問題
已經(jīng)有6人回復(fù)
- VASP表面馳豫真空層設(shè)置
已經(jīng)有9人回復(fù)
- 如何用VASP計(jì)算表面缺陷
已經(jīng)有4人回復(fù)
- 表面帶電缺陷體系VASP計(jì)算參數(shù)設(shè)置和模型搭建
已經(jīng)有8人回復(fù)
- 【求助】vasp為什么采用周期性邊界條件（或超原胞模型）來(lái)處理原子分子表面體系等，有
已經(jīng)有8人回復(fù)
- 【求助】vasp做表面結(jié)構(gòu)
已經(jīng)有11人回復(fù)
- 【求助】用vasp計(jì)算表面的DOS有什么需要注意的？
已經(jīng)有5人回復(fù)
- 【求助】能否發(fā)一份用vasp計(jì)算超晶胞某個(gè)表面的輸入文件給我
已經(jīng)有7人回復(fù)
- 【求助】如何用VASP 計(jì)算體系的表面能
已經(jīng)有11人回復(fù)
- 【求助】vasp做表面計(jì)算時(shí)真空層的厚度
已經(jīng)有13人回復(fù)
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +20 - 1ST強(qiáng)帖: 5
- 應(yīng)助: 231 (大學(xué)生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時(shí)
- 蟲號(hào): 585915
- 注冊(cè): 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 第一性原理
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
菲頂頂: 金幣+10, ★★★很有幫助 2014-08-26 16:10:39
fzx2008: 金幣+5, 謝謝指導(dǎo) 2014-08-27 13:35:33
感謝參與，應(yīng)助指數(shù) +1
菲頂頂: 金幣+10, ★★★很有幫助 2014-08-26 16:10:39
fzx2008: 金幣+5, 謝謝指導(dǎo) 2014-08-27 13:35:33
|
a。先建立一個(gè)表面超胞模型(可以加入對(duì)稱性嗎？因?yàn)槿绻雽?duì)稱性，就變成超胞的兩邊是真空層了)； > 可加入，也可不加入，看計(jì)算題目需求 "週期" 性系統(tǒng)的計(jì)算，不管你加不加對(duì)稱，超胞的兩邊永遠(yuǎn)都是真空層 b。然后測(cè)試真空層和原子層數(shù)，這里的收斂標(biāo)準(zhǔn)是能量小于10eV（測(cè)試原子層數(shù)的時(shí)候可以不對(duì)稱嗎？）； > 收斂標(biāo)準(zhǔn)能量小於 10 eV??? 你確定嬤??? 能量小於 1E-04 ? 1E-05? 同 a，要不要對(duì)稱要看你現(xiàn)在做的題目，目標(biāo)是做什麼 大多數(shù)的計(jì)算當(dāng)然是要考慮對(duì)稱，但某些情況是會(huì)特意去 broken symmetry k點(diǎn)的測(cè)試在什么時(shí)候測(cè)？ > 在做 bulk 計(jì)算的時(shí)後測(cè)的 建構(gòu)表面結(jié)構(gòu)，其中一個(gè)重要參數(shù)就是 lattice constant 要求 lattice constant ，就要做 bulk 計(jì)算 做 bulk 計(jì)算之前就要測(cè)過 ENCUT and k-points k點(diǎn)在垂直晶面的方向上為1，那么如果是110面，k點(diǎn)該怎么設(shè)呢？ > k 點(diǎn)選取跟你所計(jì)算的面是沒有關(guān)係的 不管是 110, 100, 111, 223, 233..... 任何表面 真空層方向的 k 點(diǎn)就是 1 k 點(diǎn)選取是與 cell 大小成反比 cell 三個(gè)邊 (x, y, z) 邊長(zhǎng)越長(zhǎng)，k 點(diǎn)需要越少 |
» 本帖已獲得的紅花（最新10朵）
3樓2014-08-26 15:49:21
zxzdai
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 143
- 帖子: 42
- 在線: 9.9小時(shí)
- 蟲號(hào): 3451422
- 注冊(cè): 2014-10-01
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
30樓2018-12-12 14:00:01
送紅花一朵 |
非常感謝你的回答 b，我寫錯(cuò)了是10mev 那在后面測(cè)試表面能的時(shí)候不是要求slab的原子數(shù)比和塊體材料的一樣嗎？那如果不對(duì)稱的話原子比不就不一樣了嗎，不一樣的話用表面能公式來(lái)計(jì)算的話算遲來(lái)的結(jié)果不就沒有意義了嗎？ 關(guān)于表面能的測(cè)試我說的是對(duì)的嗎？這幾天查看了一些資料，感覺有的人說是在原子層測(cè)試的時(shí)候，增加原子層數(shù)，測(cè)試表面能；有的說是弛豫的時(shí)候，改變固定原子層數(shù)，看表面能是否收斂。到底是怎樣的呢？ 可能是我才疏學(xué)淺，有時(shí)候看完了以后理解不到位，仍然希望大家能夠幫我解答，我們一起討論。 |
4樓2014-08-26 16:09:52
2樓2014-08-26 09:57:52
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +20 - 1ST強(qiáng)帖: 5
- 應(yīng)助: 231 (大學(xué)生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時(shí)
- 蟲號(hào): 585915
- 注冊(cè): 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 第一性原理
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
菲頂頂: 金幣+20, ★有幫助 2014-08-26 19:05:23
菲頂頂: 金幣+20, ★有幫助 2014-08-26 19:05:23
|
????????????????????????slab?????????????????????? > ?? ?]???@???°?... ?y?????????N????K?????????????? ??????????????????????????????????????????????????????????????????????????? > ???????? E(surface) = [ E(surface system) - N*E(bulk energy per atom) ] / 2A E(surface system) - ???????Q?]???P?S???????????? N - ????Y??????????? E(bulk energy per atom) - ?@??????????????????}???? ?? bulk ???????w??? (??? Ag FCC ?Y??) ??]????N???_???? ?? bulk ??????w??? (??? Si diamond ?Y??) ??~ ?????? 2 ??????@???? bulk energy per atom A - ??????e ????~??????????????]?б???c?K????? ?е????????????????????????????????????????? ?е??????????????????????????????????????????????????????? > ?@????????????????N????? ????y??r???????????S?????@?????????^??"???"?? ??S??Щ???J??????y??r????????????? ?@????????????????? ??N????????????Y??????????? ???c?????????? ??N????"???" ??????????N????????? |
5樓2014-08-26 16:54:44
|
E（bulk）是體相的能量，你說的是指體相中的單個(gè)能量。你說的應(yīng)該是單質(zhì)表面能的計(jì)算吧，計(jì)算單質(zhì)沒有問題，如果有多種原子，想計(jì)算他的表面，那就不一樣了啊。那你說的這種計(jì)算應(yīng)該就是他們說的在有多種原子的時(shí)候就需要考慮原子的化學(xué)式了，也就是說用原子能量代替體相能量。那這種情況下就應(yīng)該是不用考慮表面對(duì)不對(duì)稱了。 原子層測(cè)試是判斷大概去多少層原子，這里應(yīng)該是靜態(tài)計(jì)算；而弛豫則是對(duì)表面結(jié)構(gòu)做優(yōu)化判斷中間固定多少層原子才能保證體相結(jié)構(gòu)的性質(zhì)，我認(rèn)為他們是不同的， |
6樓2014-08-26 17:39:55
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +20 - 1ST強(qiáng)帖: 5
- 應(yīng)助: 231 (大學(xué)生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時(shí)
- 蟲號(hào): 585915
- 注冊(cè): 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 第一性原理
★ ★ ★
fzx2008: 金幣+3, 謝謝指導(dǎo) 2014-08-27 13:39:21
fzx2008: 金幣+3, 謝謝指導(dǎo) 2014-08-27 13:39:21
|
E（bulk）是體相的能量，你說的是指體相中的單個(gè)能量。你說的應(yīng)該是單質(zhì)表面能的計(jì)算吧，計(jì)算單質(zhì)沒有問題，如果有多種原子，想計(jì)算他的表面，那就不一樣了啊。那你說的這種計(jì)算應(yīng)該就是他們說的在有多種原子的時(shí)候就需要考慮原子的化學(xué)式了，也就是說用原子能量代替體相能量。那這種情況下就應(yīng)該是不用考慮表面對(duì)不對(duì)稱了。 > 嗯嗯~沒錯(cuò)~你點(diǎn)到問題了 若是合金表面能計(jì)算 表面原子數(shù)與塊材原子數(shù)若是有成比例 那的確會(huì)"比較容易"去計(jì)算合金表面能 但~仍然不是必要的 原子數(shù)不成比例?K不影響表面能計(jì)算 原子層測(cè)試是判斷大概去多少層原子，這里應(yīng)該是靜態(tài)計(jì)算；而弛豫則是對(duì)表面結(jié)構(gòu)做優(yōu)化判斷中間固定多少層原子才能保證體相結(jié)構(gòu)的性質(zhì)，我認(rèn)為他們是不同的 > 嗯嗯~這我沒有意見 原子層測(cè)試用靜態(tài)式計(jì)算或是用弛豫式計(jì)算，我覺得都對(duì)，只是思考角度不同 (我這裡指的是原子是否允許移動(dòng)，靜態(tài)是不允許，弛豫是允許 弛豫指的是 relaxation or optimization 對(duì)吧?) 我想~我們兩討論的的方向似乎不太一樣，呵呵 |
7樓2014-08-26 18:01:05
8樓2014-08-26 19:05:18
9樓2014-08-26 19:20:08
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +20 - 1ST強(qiáng)帖: 5
- 應(yīng)助: 231 (大學(xué)生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時(shí)
- 蟲號(hào): 585915
- 注冊(cè): 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 第一性原理
10樓2014-08-26 19:45:11
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 菲頂頂 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑 +9 | 吃吃吃才有意義 2026-03-19 | 9/450 |
|
|
|
[考研] 299求調(diào)劑 +4 | shxchem 2026-03-20 | 6/300 |
|
|
|
[考研] 材料 271求調(diào)劑 +3 | 展信悅_ 2026-03-21 | 3/150 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +7 | 楓橋ZL 2026-03-18 | 9/450 |
|
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
|
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|