24小时热门版块排行榜    

查看: 1643  |  回复: 4

gengmingIG

至尊木虫 (正式写手)

[求助] 计算表面能疑惑

在计算材料(比如Mg2SiO4)表面能的时候,几乎所有的文献和书籍都给出
E(表面)=[E(slab)-n*E(bulk)]/2A ,

按照我的理解应该是晶体一个unit cell的能量。
但当为了保证表面结构,无法构建成unit cell整数倍厚度的slab模型时, 这个又应该是什么呢?n显然不能取非整数吧?

不知道有没有在这个问题上有经验的,还望指教一下。或者有没有介绍表面能比较详细的文献或书籍给推荐一下。多谢。
回复此楼

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

第一性原理计算-tg 表界面计算 vasp 第一性原理

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

devinhu

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gengmingIG: 金币+15, 有帮助, 多谢,我先看看文章。 2017-03-08 19:29:58
10.1021/acs.jpcc.5b05136
我是按照这篇文章里面来的
2楼2017-03-08 14:53:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gengmingIG: 金币+25, ★★★很有帮助, 看起来意思是变成热力学问题处理,不知道有没有具体的文献或书籍具体介绍这样方法,还请赐教。 2017-03-08 19:32:07
这个要计算 好些相关化合物的结构和总能 :
1. Mg2SiO4,
2. MgO
3. Mg2Si
4. SiO2
5. Si
6. O2
然后做出它们的相图。确定化学势mu_Mg ,mu_Si ,可取的范围。
将Ebulk替换为 n_Mg * mu_Mg + n_Si * mu_Si + n_O * mu_O.
n_x(x=Mg, Si, O)是原子x的个数。 mu_x 是原子x的化学势。
再来计算表面能。
3楼2017-03-08 15:47:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gengmingIG

至尊木虫 (正式写手)

引用回帖:
3楼: Originally posted by valenhou001 at 2017-03-08 15:47:56
这个要计算 好些相关化合物的结构和总能 :
1. Mg2SiO4,
2. MgO
3. Mg2Si
4. SiO2
5. Si
6. O2
然后做出它们的相图。确定化学势mu_Mg ,mu_Si ,可取的范围。
将Ebulk替换为 n_Mg * mu_Mg + n_Si * mu_Si + ...

不知道有没有具体的文献或书籍具体介绍这样方法,还请赐教。
4楼2017-03-09 09:47:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
gengmingIG: 金币+15, 有帮助, 多谢了。 2017-03-09 14:21:36
引用回帖:
4楼: Originally posted by gengmingIG at 2017-03-09 10:47:50
不知道有没有具体的文献或书籍具体介绍这样方法,还请赐教。...

有很多这样的文献(google)。还没有见到有专门的书籍来介绍,有的话,也是一些文章的综合。
上一个帖子(devinhu)提到的就是一篇好的范例,只不过针对的二元化合物。但是可以推广到三维体系。表面的稳定性随压强的影响,一般都是考虑某个元素的气态下的化学势随压强的影响的。

至于针对具体问题(研究案例),需要个别计算操作。
5楼2017-03-09 10:16:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gengmingIG 的主题更新
信息提示
请填处理意见