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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助】表面能计算值为负？
52/1返回列表
查看: 1906  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 宏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
     很显然，算出来的结果是负的，有没有哪位大侠有相关经验的，或者知道哪错的？十分感谢了
回复此楼

» 猜你喜欢
1楼 2010-01-28 09:45:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hgye

金虫 (小有名气)
★ ★
宏hong(金币+5):好，谢谢提示 2010-01-29 14:21
ice_rain(金币+2):谢谢交流 2010-01-29 20:11
8楼说的很对，你的问题就是原子数目不一致。在建立slab的时候要考虑到这个问题，比如是AB型化合物，那么slab中A和B原子数目之比要等于bulk原包中A、B数目之比。如果确实无法做到这一点，就要比较formation energy，你看看关于表面缺陷结构计算的文章，那种情况下原子数目比大多都不一致。
赞一下(1人) 回复此楼
12楼2010-01-29 12:30:49
已阅   回复此楼   关注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都是什么？考虑全表面的自能了吗？
赞一下(7人) 回复此楼
先进使役材料计算（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的回帖
查看全部 15 个回答
提示: 如果您在30分钟内回复过其他散金贴，则可能无法领取此贴金币
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 11:03 by xSrQNSSjcgoX
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 09:46 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-15 07:56 by xSrQNSSjcgoX
[硕博家园] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 SX5mJmPWDfOO 2026-11-13 4/200 2026-11-15 07:34 by xSrQNSSjcgoX
[硕博家园] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 XJXMuRjD7cK7 2026-11-13 4/200 2026-11-15 06:50 by xSrQNSSjcgoX
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-15 05:33 by xSrQNSSjcgoX
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 uKOSCocCb6Fm 2026-11-12 5/250 2026-11-15 03:21 by xSrQNSSjcgoX
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +6 FChKX8meu1ds 2026-11-09 10/500 2026-11-13 11:08 by Equinoxhua
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:51 by ECLTXnJnRSXf
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 TgvkbZZl9tqB 2026-11-11 6/300 2026-11-12 12:42 by ECLTXnJnRSXf
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:11 by HajghKOMynmz
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:00 by HajghKOMynmz
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 7/350 2026-11-12 10:49 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 11/550 2026-11-12 09:57 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +7 FChKX8meu1ds 2026-11-09 12/600 2026-11-12 09:43 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +8 cqv1zKnzd7hr 2026-11-09 12/600 2026-11-12 09:32 by HajghKOMynmz
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 3/150 2026-11-12 08:26 by HajghKOMynmz
[考研] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 cqv1zKnzd7hr 2026-11-09 8/400 2026-11-12 05:08 by 9fOTqj074tIg
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 0TEMgxdT1eIw 2026-11-09 8/400 2026-11-12 05:00 by 9fOTqj074tIg
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-12 04:57 by 9fOTqj074tIg
信息提示
关闭
请填处理意见
关闭