24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 6236  |  回复: 30
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

菲顶顶

铜虫 (正式写手)

[求助] vasp表面计算详细解说 已有1人参与

想和大家交流下vasp计算表面的一些问题:
1. 计算表面的步骤:
                             a。先建立一个表面超胞模型(可以加入对称性吗?因为如果引入对称性,就变成超胞的两边是真空层了);
                             b。然后测试真空层和原子层数,这里的收敛标准是能量小于10eV(测试原子层数的时候可以不对称吗?);
                             c。然后进行表面结构的弛豫,固定中间几层原子,观察层间距的变化,若固定的层数之间有较大的位移(原子既然固定了,固定原子层之间还会有有位移嘛?),就增加弛豫原子层数(是固定原子的层数还是没固定原子的层数),在观察层间距的变化,知道变化很小为止(这里是做静态计算吗?然后从OUTcar里读取原子位置,看层之间的变化?这里有一个什么标准吗?);
                            d。此外,还要测试表面能,这时要保证slab足够厚,当slab层数变大时,表面能变化很小就行了(这里表面能使需要用公式计算的吧~~~,这个和原子层数测试是同时进行的吗?边看总能收不收敛,还要看表面能收不收连?而且说要保证晶胞大小一致,那这样真空层不就是变了吗,那前面测试真空层不就没有意义了)。
                            e。计算表面的性质
括号里面的就是我的疑问,还有一个问题就是k点的测试在什么时候测,k点在垂直晶面的方向上为1,那么如果是110面,k点该怎么设呢?
我也是在学习的起步阶段,最近看这些有些糊涂,思路不太清晰,好像小木虫上也没有详细的帖子,希望我们大家交流一下,
回复此楼

» 收录本帖的淘帖专辑推荐

好书!好读! 仿真建模与计算 第一性原理相关文档 第一性原理计算-tg
第一性原理和电化学 表界面计算

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
菲顶顶: 金币+10, ★★★很有帮助 2014-08-26 16:10:39
fzx2008: 金币+5, 谢谢指导 2014-08-27 13:35:33
a。先建立一个表面超胞模型(可以加入对称性吗?因为如果引入对称性,就变成超胞的两边是真空层了);
> 可加入,也可不加入,看計算題目需求
   "週期" 性系統的計算,不管你加不加對稱,超胞的兩邊永遠都是真空層

b。然后测试真空层和原子层数,这里的收敛标准是能量小于10eV(测试原子层数的时候可以不对称吗?);
>  收斂標準能量小於 10 eV??? 你確定嬤???
    能量小於 1E-04 ?  1E-05?  
    同 a,要不要對稱要看你現在做的題目,目標是做什麼
    大多數的計算當然是要考慮對稱,但某些情況是會特意去 broken symmetry

k点的测试在什么时候测?
>  在做 bulk 計算的時後測的
    建構表面結構,其中一個重要參數就是 lattice constant
    要求 lattice constant ,就要做 bulk 計算
    做 bulk 計算之前就要測過 ENCUT and k-points

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

» 本帖已获得的红花(最新10朵)

3楼2014-08-26 15:49:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 31 个回答

菲顶顶

铜虫 (正式写手)

我先来自己回答一部分解决了的问题吧

                           
                             c。这里是弛豫,所以就是优化结构的计算了,固定原子层数,看固定的原子层数的间距,如果变化较大的话,说明固定的原子层数太少,需要增加固定的原子层数,再弛豫,看层间距变化,知道层间距变化很小为止。
                            d。这里是需要固定不同层数,看表面能的变化,这个是在和弛豫表面的时候一起做的,所以超晶胞的大小是不变的。

a,b和k点的问题还没有解决,希望有人来解答。
2楼2014-08-26 09:57:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)

送红花一朵
引用回帖:
3楼: Originally posted by jpchou at 2014-08-26 15:49:21
a。先建立一个表面超胞模型(可以加入对称性吗?因为如果引入对称性,就变成超胞的两边是真空层了);
> 可加入,也可不加入,看計算題目需求
   "週期" 性系統的計算,不管你加不加對稱,超胞的兩邊永 ...

非常感谢你的回答
b,我写错了是10mev
那在后面测试表面能的时候不是要求slab的原子数比和块体材料的一样吗?那如果不对称的话原子比不就不一样了吗,不一样的话用表面能公式来计算的话算迟来的结果不就没有意义了吗?
关于表面能的测试我说的是对的吗?这几天查看了一些资料,感觉有的人说是在原子层测试的时候,增加原子层数,测试表面能;有的说是弛豫的时候,改变固定原子层数,看表面能是否收敛。到底是怎样的呢?
可能是我才疏学浅,有时候看完了以后理解不到位,仍然希望大家能够帮我解答,我们一起讨论。
4楼2014-08-26 16:09:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
菲顶顶: 金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂求调剂 +7 121. 2026-04-02 7/350 2026-04-02 19:15 by dick_runner
[论文投稿] chinese chemical letters英文版投稿求助 120+4 Yishengeryi 2026-03-30 6/300 2026-04-02 17:19 by Yishengeryi
[考研] 302求调剂一志愿华中师范大学 +5 小江小江江江 2026-04-02 5/250 2026-04-02 15:47 by zikuzi1853
[考研] 材料调剂 +6 一样YWY 2026-04-02 6/300 2026-04-02 14:56 by liu823948201
[考研] 326求调剂 +3 崽崽仔 2026-04-02 3/150 2026-04-02 14:54 by 1753564080
[考研] 学硕机械工程303求调剂 +4 无名所以叫吴明 2026-03-30 4/200 2026-04-02 14:06 by asdfzly
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +4 哇呼哼呼哼 2026-04-02 4/200 2026-04-02 13:10 by yulian1987
[考研] 材料化工340求调剂 +5 jhx777 2026-03-30 5/250 2026-04-02 12:45 by smileboy2006
[考研] 289求调剂 +23 新时代材料 2026-03-27 26/1300 2026-04-02 10:29 by oooqiao
[考研] 【求调剂】新能源材料本科,一志愿211,初试321 +6 求调剂学校, 2026-04-02 6/300 2026-04-02 09:41 by 晴空210210
[考研] 考研调剂 +12 Amber00 2026-03-31 12/600 2026-04-02 09:04 by sanrepian
[考研] 一志愿 南京航空航天大学 ,080500材料科学与工程学硕 +7 @taotao 2026-03-30 7/350 2026-04-01 14:30 by chenqifeng666
[考研] 安徽大学专硕生物与医药专业(086000)324分,英语已过四六级,六级521,求调剂 +10 美味可乐鸡翅 2026-03-26 12/600 2026-04-01 11:17 by syh9288
[考研] 318一志愿吉林大学生物与医药 求调剂 +6 笃行致远. 2026-03-28 6/300 2026-04-01 09:28 by oooqiao
[考研] 299求调剂 +8 嗯嗯嗯嗯2 2026-03-27 8/400 2026-03-31 18:23 by lizhi8172
[考研] 318求调剂 +10 陈晨79 2026-03-30 10/500 2026-03-31 17:37 by 544594351
[考研] 0856 335分 +9 cccchenso 2026-03-29 9/450 2026-03-31 16:37 by lishahe
[考研] 调剂求院校招收 +7 鹤鲸鸽 2026-03-28 7/350 2026-03-31 11:21 by oooqiao
[考研] 348求调剂 +6 小懒虫不懒了 2026-03-28 6/300 2026-03-30 10:29 by Evan_Liu
[考研] 求调剂 +10 张zz111 2026-03-27 11/550 2026-03-30 09:17 by 无际的草原
信息提示
请填处理意见