版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

菲顶顶

铜虫 (正式写手)

应助: 2 (幼儿园)
金币: 32.6
散金: 45
红花: 10
帖子: 631
在线: 250.2小时
虫号: 1462090
注册: 2011-10-26
专业: 凝聚态物性 II ：电子结构

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

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

回复此楼

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

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

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

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

菲顶顶

铜虫 (正式写手)

应助: 2 (幼儿园)
金币: 32.6
散金: 45
红花: 10
帖子: 631
在线: 250.2小时
虫号: 1462090
注册: 2011-10-26
专业: 凝聚态物性 II ：电子结构

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


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

a，b和k点的问题还没有解决，希望有人来解答。

赞一下

回复此楼

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

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

菲顶顶

铜虫 (正式写手)

应助: 2 (幼儿园)
金币: 32.6
散金: 45
红花: 10
帖子: 631
在线: 250.2小时
虫号: 1462090
注册: 2011-10-26
专业: 凝聚态物性 II ：电子结构

送红花一朵

引用回帖:

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

非常感谢你的回答
b，我写错了是10mev
那在后面测试表面能的时候不是要求slab的原子数比和块体材料的一样吗？那如果不对称的话原子比不就不一样了吗，不一样的话用表面能公式来计算的话算迟来的结果不就没有意义了吗？
关于表面能的测试我说的是对的吗？这几天查看了一些资料，感觉有的人说是在原子层测试的时候，增加原子层数，测试表面能；有的说是弛豫的时候，改变固定原子层数，看表面能是否收敛。到底是怎样的呢？
可能是我才疏学浅，有时候看完了以后理解不到位，仍然希望大家能够帮我解答，我们一起讨论。

赞一下(1人)

回复此楼

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

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

菲顶顶

铜虫 (正式写手)

应助: 2 (幼儿园)
金币: 32.6
散金: 45
红花: 10
帖子: 631
在线: 250.2小时
虫号: 1462090
注册: 2011-10-26
专业: 凝聚态物性 II ：电子结构

引用回帖:

5楼: Originally posted by jpchou at 2014-08-26 16:54:44
那在后面测试表面能的时候不是要求slab的原子数比和块体材料的一样吗？
> ?? 沒有這回事吧...
測表面能為什麼要跟塊材原子數比一樣???

那如果不对称的话原子比不就不一样了吗，不一样的话用表面能公式来 ...

E（bulk）是体相的能量，你说的是指体相中的单个能量。你说的应该是单质表面能的计算吧，计算单质没有问题，如果有多种原子，想计算他的表面，那就不一样了啊。那你说的这种计算应该就是他们说的在有多种原子的时候就需要考虑原子的化学式了，也就是说用原子能量代替体相能量。那这种情况下就应该是不用考虑表面对不对称了。
原子层测试是判断大概去多少层原子，这里应该是静态计算；而弛豫则是对表面结构做优化判断中间固定多少层原子才能保证体相结构的性质，我认为他们是不同的，

赞一下

回复此楼

6楼2014-08-26 17:39:55

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

菲顶顶

铜虫 (正式写手)

应助: 2 (幼儿园)
金币: 32.6
散金: 45
红花: 10
帖子: 631
在线: 250.2小时
虫号: 1462090
注册: 2011-10-26
专业: 凝聚态物性 II ：电子结构

引用回帖:

7楼: Originally posted by jpchou at 2014-08-26 18:01:05
E（bulk）是体相的能量，你说的是指体相中的单个能量。你说的应该是单质表面能的计算吧，计算单质没有问题，如果有多种原子，想计算他的表面，那就不一样了啊。那你说的这种计算应该就是他们说的在有多种原子的时候 ...

我们讨论的应该是一个问题吧嘿嘿我们互相交换知识才能学习么
我差不多也明白了，谢谢你哦~~~~~~·
你说的：弛豫指的是 relaxation or optimization 對吧?
是的，弛豫是这个

赞一下

回复此楼

8楼2014-08-26 19:05:18

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