计算化学公社

标题: 表面能的计算公式中的面积项求助 [打印本页]

作者
Author:
张存金    时间: 2019-12-18 19:49
标题: 表面能的计算公式中的面积项求助
计算表面能的公式不太了解,γ=[Eslab-n*Ebulk] /2A 这里的A是指的那一个面的面积,是我切得那个面的面积吗,(我疑惑在不是有很多的面吗,表面能是所有表面的能量吧,这个A是那个面啊)下面是我自己做的001面的结构,一个是2X2的单层一个是2X2的三层结构,然后我取的001面的面积,算出来的表面能一个是6.9eV/nm2,一个5.3 有点迷惑
EslabnEbulkab sinθ2表面能
2x2x3 001
-2.67E+02
12
-2.26E+01
5.8143
5.8143
0.866025404
2
6.952976767
2X2    001
-8.73E+01
4
-2.26E+01
5.8143
5.8143
0.866025404
2
5.377519464

作者
Author:
ghifi37    时间: 2019-12-19 09:59
由于周期性的缘故,只有上下才是表面,所以是2A。两个真空层厚度取得一样厚吗?切完后有没做relax?
作者
Author:
sobereva    时间: 2019-12-19 15:58
如置顶的新社员必读贴和论坛首页的公告栏所示,求助帖必须在帖子标题明确体现出“提问”、“求助”要素及具体内容(仔细看http://bbs.keinsci.com/thread-9348-1-1.html),我已把你的不恰当标题“表面能的计算”改了,以后务必注意
作者
Author:
张存金    时间: 2019-12-19 17:36
sobereva 发表于 2019-12-19 15:58
如置顶的新社员必读贴和论坛首页的公告栏所示,求助帖必须在帖子标题明确体现出“提问”、“求助”要素及具 ...

好的,谢谢

作者
Author:
张存金    时间: 2019-12-19 17:39
ghifi37 发表于 2019-12-19 09:59
由于周期性的缘故,只有上下才是表面,所以是2A。两个真空层厚度取得一样厚吗?切完后有没做relax?

真空层厚度一样,做了驰豫,但是是整体驰豫,是不是应该驰豫表面?
作者
Author:
hakuna    时间: 2019-12-19 19:29
算表面能取这么薄的表面模型,没有意义
作者
Author:
张存金    时间: 2019-12-19 19:52
hakuna 发表于 2019-12-19 19:29
算表面能取这么薄的表面模型,没有意义

不是很理解,大概多厚,根据什么来,希望您能给比较详细的解答,先谢谢您了
作者
Author:
hakuna    时间: 2019-12-19 21:16
张存金 发表于 2019-12-19 19:52
不是很理解,大概多厚,根据什么来,希望您能给比较详细的解答,先谢谢您了

逐渐增加层厚度,看看收敛情况
体系太大算不动就取1x1
作者
Author:
Hmumu    时间: 2023-7-19 10:04
请问你问题解决了吗
作者
Author:
茻苍    时间: 2023-10-18 16:27
应该就是你切之后的表面的面积




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3