菲頂頂

銅蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 32.6
散金: 45
紅花: 10
帖子: 631
在線: 250.2小時
蟲號: 1462090
注冊: 2011-10-26
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] vasp表面計算詳細解說已有1人參與

想和大家交流下vasp計算表面的一些問題：
1. 計算表面的步驟：
                           a。先建立一個表面超胞模型(可以加入對稱性嗎？因為如果引入對稱性，就變成超胞的兩邊是真空層了)；
                           b。然后測試真空層和原子層數(shù)，這里的收斂標準是能量小于10eV（測試原子層數(shù)的時候可以不對稱嗎？）；
                           c。然后進行表面結(jié)構(gòu)的弛豫，固定中間幾層原子，觀察層間距的變化，若固定的層數(shù)之間有較大的位移（原子既然固定了，固定原子層之間還會有有位移嘛？），就增加弛豫原子層數(shù)（是固定原子的層數(shù)還是沒固定原子的層數(shù)），在觀察層間距的變化，知道變化很小為止（這里是做靜態(tài)計算嗎？然后從OUTcar里讀取原子位置，看層之間的變化？這里有一個什么標準嗎？）；
                        d。此外，還要測試表面能，這時要保證slab足夠厚，當slab層數(shù)變大時，表面能變化很小就行了（這里表面能使需要用公式計算的吧~~~，這個和原子層數(shù)測試是同時進行的嗎？邊看總能收不收斂，還要看表面能收不收連？而且說要保證晶胞大小一致，那這樣真空層不就是變了嗎，那前面測試真空層不就沒有意義了）。
                        e。計算表面的性質(zhì)
括號里面的就是我的疑問,還有一個問題就是k點的測試在什么時候測，k點在垂直晶面的方向上為1，那么如果是110面，k點該怎么設(shè)呢？
我也是在學習的起步階段，最近看這些有些糊涂，思路不太清晰，好像小木蟲上也沒有詳細的帖子，希望我們大家交流一下，

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

xiewei娜

» 猜你喜歡

288求調(diào)劑，一志愿華南理工大學071005 已經(jīng)有5人回復(fù)
一志愿中國海洋大學，生物學，301分，求調(diào)劑已經(jīng)有5人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有12人回復(fù)
0856調(diào)劑，是學校就去已經(jīng)有5人回復(fù)
0817 化學工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有13人回復(fù)
梁成偉老師課題組歡迎你的加入已經(jīng)有11人回復(fù)
復(fù)試調(diào)劑已經(jīng)有6人回復(fù)
0703化學調(diào)劑，六級已過，有科研經(jīng)歷已經(jīng)有11人回復(fù)
321求調(diào)劑已經(jīng)有10人回復(fù)
東華理工大學化材專業(yè)26屆碩士博士申請已經(jīng)有8人回復(fù)

1樓 2014-08-25 19:53:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jpchou

專家顧問 (著名寫手)

專家經(jīng)驗: +20
1ST強帖: 5
應(yīng)助: 231 (大學生)
金幣: 2342
散金: 4367
紅花: 103
帖子: 1061
在線: 757.3小時
蟲號: 585915
注冊: 2008-08-08
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
管轄: 第一性原理

???????:

24?: Originally posted by ????? at 2017-02-19 19:03:33
????????????
3??????????????? ??μ?P32??
?????????????????K??mesh?????????
???????????????????????λ????????????????????????????????????Selective Dynamic??
??????P ...

???????????????????
> ??????]??  "???" ????
?????y????Щ?S???

??????????????10meV???????????????????С?????????????????ο????
> ??N??  "????????" ?  ??????????v "???" ?
??????? 10 meV??????????v???? total energy ???  energy per atom
?????v total energy convergence is 10 meV??????????????y?e???????????

?????????f?????????? energy per atom
energy per atom ???? 10 meV
??????????]???????????????????????S???????????
??????f?????@Щ?l?????????????r???]?k??????κ???}

贊一下

回復(fù)此樓

27樓2017-02-20 00:39:03

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 31 個回答

菲頂頂

銅蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 32.6
散金: 45
紅花: 10
帖子: 631
在線: 250.2小時
蟲號: 1462090
注冊: 2011-10-26
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

我先來自己回答一部分解決了的問題吧


                           c。這里是弛豫，所以就是優(yōu)化結(jié)構(gòu)的計算了，固定原子層數(shù)，看固定的原子層數(shù)的間距，如果變化較大的話，說明固定的原子層數(shù)太少，需要增加固定的原子層數(shù)，再弛豫，看層間距變化，知道層間距變化很小為止。
                        d。這里是需要固定不同層數(shù)，看表面能的變化，這個是在和弛豫表面的時候一起做的，所以超晶胞的大小是不變的。

a，b和k點的問題還沒有解決，希望有人來解答。

贊一下

回復(fù)此樓

2樓2014-08-26 09:57:52

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jpchou

專家顧問 (著名寫手)

專家經(jīng)驗: +20
1ST強帖: 5
應(yīng)助: 231 (大學生)
金幣: 2342
散金: 4367
紅花: 103
帖子: 1061
在線: 757.3小時
蟲號: 585915
注冊: 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

a。先建立一個表面超胞模型(可以加入對稱性嗎？因為如果引入對稱性，就變成超胞的兩邊是真空層了)；
> 可加入，也可不加入，看計算題目需求
"週期" 性系統(tǒng)的計算，不管你加不加對稱，超胞的兩邊永遠都是真空層

b。然后測試真空層和原子層數(shù)，這里的收斂標準是能量小于10eV（測試原子層數(shù)的時候可以不對稱嗎？）；
>  收斂標準能量小於 10 eV??? 你確定嬤???
能量小於 1E-04 ?  1E-05?
同 a，要不要對稱要看你現(xiàn)在做的題目，目標是做什麼
大多數(shù)的計算當然是要考慮對稱，但某些情況是會特意去 broken symmetry

k點的測試在什么時候測？
>  在做 bulk 計算的時後測的
建構(gòu)表面結(jié)構(gòu)，其中一個重要參數(shù)就是 lattice constant
要求 lattice constant ，就要做 bulk 計算
做 bulk 計算之前就要測過 ENCUT and k-points

k點在垂直晶面的方向上為1，那么如果是110面，k點該怎么設(shè)呢？
> k 點選取跟你所計算的面是沒有關(guān)係的
不管是 110, 100, 111, 223, 233..... 任何表面
真空層方向的 k 點就是 1
k 點選取是與 cell 大小成反比
cell 三個邊 (x, y, z) 邊長越長，k 點需要越少

贊一下(2人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

菲頂頂

3樓2014-08-26 15:49:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

菲頂頂

銅蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 32.6
散金: 45
紅花: 10
帖子: 631
在線: 250.2小時
蟲號: 1462090
注冊: 2011-10-26
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

送紅花一朵

引用回帖:

3樓: Originally posted by jpchou at 2014-08-26 15:49:21
a。先建立一個表面超胞模型(可以加入對稱性嗎？因為如果引入對稱性，就變成超胞的兩邊是真空層了)；
> 可加入，也可不加入，看計算題目需求
"週期" 性系統(tǒng)的計算，不管你加不加對稱，超胞的兩邊永 ...

非常感謝你的回答
b，我寫錯了是10mev
那在后面測試表面能的時候不是要求slab的原子數(shù)比和塊體材料的一樣嗎？那如果不對稱的話原子比不就不一樣了嗎，不一樣的話用表面能公式來計算的話算遲來的結(jié)果不就沒有意義了嗎？
關(guān)于表面能的測試我說的是對的嗎？這幾天查看了一些資料，感覺有的人說是在原子層測試的時候，增加原子層數(shù)，測試表面能；有的說是弛豫的時候，改變固定原子層數(shù)，看表面能是否收斂。到底是怎樣的呢？
可能是我才疏學淺，有時候看完了以后理解不到位，仍然希望大家能夠幫我解答，我們一起討論。

贊一下(1人)

回復(fù)此樓

4樓2014-08-26 16:09:52

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 0817 化學工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +13	rare12345 2026-03-18	13/650	2026-03-19 17:30 by jean5056
[考博] 東華理工大學化材專業(yè)26屆碩士博士申請 +8	zlingli 2026-03-13	8/400	2026-03-19 16:32 by 輕松不少隨
[考研] 材料專碩274一志愿陜西師范大學求調(diào)劑 +8	薛云鵬 2026-03-13	8/400	2026-03-19 15:36 by haoshis
[考研] 267一志愿南京工業(yè)大學0817化工求調(diào)劑 +10	SUICHILD 2026-03-12	10/500	2026-03-19 09:51 by Delta2012
[考研] 0703化學調(diào)劑，求各位老師收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 311求調(diào)劑 +4	冬十三 2026-03-18	4/200	2026-03-18 21:47 by 盡舜堯1
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7	步川酷紫123 2026-03-13	7/350	2026-03-18 17:12 by 盡舜堯1
[考研] 0703化學調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[考研] 268求調(diào)劑 +6	簡單點0 2026-03-17	6/300	2026-03-18 09:04 by 無際的草原
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 296求調(diào)劑 +5	大口吃飯身體健 2026-03-13	5/250	2026-03-17 21:05 by 不惑可樂
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 機械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 26考研一志愿中國石油大學(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 266求調(diào)劑 +4	學員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz