24小时热门版块排行榜    

查看: 1513  |  回复: 14
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 宏hong 的 28 个金币 ,回帖就立即获得 1 个金币,每人有 1 次机会
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

宏hong

木虫 (正式写手)

[交流] 【求助】表面能计算值为负?

我正在用Dmol3做一个氧化物的表面能,采用表面2层驰豫,
     5层结构优化后得出的total energy为-18795.187382Ha,
     3层结构优化后得出的total energy为-12429.5106102Ha,
我采用的公式如下:
                 S=[E(n)-n*Ub]/2A
                 Ub=[U(n)-U(n-2)]/2
算出的Ub=-3182.838385Ha
带入S公式,则得出E(n)-n*Ub=-18795.187382Ha-(-3182.838385Ha)*5=      -2880.995455
     很显然,算出来的结果是负的,有没有哪位大侠有相关经验的,或者知道哪错的?十分感谢了
回复此楼

» 猜你喜欢

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

liumiao05

至尊木虫 (知名作家)

秦武汉风


宏hong(金币+1):谢谢参与
宏hong(金币+8):十分感谢 2010-01-30 08:36
如8楼和12楼所说,那个公式对单质成立,对化合物的话不是每种表面都可以,因为选slab的时候不能同时兼顾原子比例和表面组成(就是说两个上下两个面都是你计算的表面)
举个例子来说,对于CsCl结构的B2晶体来说,如果你计算(100)面的表面能的话,你选的slab就不可能保证原子比1:1,外层原子肯定要多
所以需要计算化学势,具体的话就涉及到hgye说的formation energy 了,你可以看看关于point defect(点缺陷)的文献
比如这篇
Concentrations of Atomic Defects in B2FexAl1-x. An Ab-Initio Study
J. MAYER, ELSASSER, M. FAHNLE
phys. stat. sol. (b) 191, 283 (1995)
天下英雄出我辈,一入江湖岁月催。鸿图霸业谈笑间,不胜人生一场醉。提剑跨骑挥鬼雨,白骨如山鸟惊飞。尘世如潮人如水,只叹江湖几人回。
13楼2010-01-29 20:41:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

lfhuang

木虫 (著名写手)

★ ★ ★
宏hong(金币+2):谢谢 2010-01-28 10:08
ice_rain(金币+1):谢谢讨论交流 2010-01-29 20:09
S,E(n),U(n),Ub,A都是什么?考虑全表面的自能了吗?
先进使役材料计算(https://aidme.nimte.ac.cn)
2楼2010-01-28 10:04:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宏hong

木虫 (正式写手)

引用回帖:
Originally posted by lfhuang at 2010-01-28 10:04:28:
S,E(n),U(n),Ub,A都是什么?考虑全表面的自能了吗?

S是表面能,E(n)是结构优化后从out文件中找到的total energy,Ub如题中的公式,A是表面的面积。(这是我的理解)
我不知道你指的全表面的自能是指什么?因为我看文献中计算表面能时,采用的也是这个公式。
十分感谢了
3楼2010-01-28 10:08:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hgye

金虫 (小有名气)

宏hong(金币+5): 2010-01-28 18:09
宏hong(金币+1):我试一下,谢谢 2010-01-28 18:11
aylayl08(金币+1):谢谢讨论 2010-01-28 19:04
原因可能是你的层数太少,计算出的Ub不准。特别是三层的模型不够,导致E(3)偏大(绝对值小),Ub的绝对值太大。
另外表面能的计算可以不这样求Ub, 直接用具有相同原子数的体材能量代替,更多的人应该是这样计算的,比你的方法计算量小,不用计算两个表面模型。你可以先这样检验一下你的结果。
4楼2010-01-28 17:45:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
提示: 如果您在30分钟内回复过其他散金贴,则可能无法领取此贴金币
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见