| 查看: 1509 | 回复: 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 很显然，算出来的结果是负的，有没有哪位大侠有相关经验的，或者知道哪错的？十分感谢了 |
回复此楼» 猜你喜欢
博士延得我，科研能力直往上蹿
已经有7人回复
-
退学或坚持读
已经有27人回复
-
面上基金申报没有其他的参与者成吗
已经有5人回复
-
有70后还继续奋斗在职场上的吗？
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复
-
多组分精馏求助
已经有6人回复
acridine
木虫 (正式写手)
- 1ST强帖: 11
- 应助: 13 (小学生)
- 金币: 4901.6
- 红花: 13
- 帖子: 706
- 在线: 389.5小时
- 虫号: 305478
- 注册: 2006-12-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★
aylayl08(金币+3):谢谢兄台指点 2010-01-29 09:07
宏hong(金币+5): 2010-01-29 09:09
aylayl08(金币+3):谢谢兄台指点 2010-01-29 09:07
宏hong(金币+5): 2010-01-29 09:09
|
得到的数值这么大肯定是原子数目的问题了，我怀疑是不是你选择的表面模型不对，导致模型内的原子数与体材对应不上？ 简单点说就是：你计算Eslab-n*Ebulk的时候，这个Eslab体系所含有的原子种类和数目是否等于n个bulk模型内的原子种类和数目？ 注意原子种类也要一致，否则出现E("12个A原子+8个B原子" 的slab)-5*E("2个A原子+2个B原子" 的bulk) 这种减法，此时得到的结果是没有意义的。 原子种类和数目对应正确，仅仅原子之间结合有不同，是不会导致Eslab和n*Ebulk的总能产生如此大的差别的  [ Last edited by acridine on 2010-1-29 at 01:03 ] |
8楼2010-01-29 00:57:31
lfhuang
木虫 (著名写手)
- 应助: 39 (小学生)
- 金币: 4344.2
- 散金: 1546
- 红花: 36
- 帖子: 1403
- 在线: 676.7小时
- 虫号: 505824
- 注册: 2008-02-17
- 性别: GG
- 专业: 金属材料表面科学与工程
2楼2010-01-28 10:04:28
宏hong
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2754.9
- 散金: 305
- 红花: 1
- 沙发: 1
- 帖子: 649
- 在线: 237.6小时
- 虫号: 842750
- 注册: 2009-09-07
- 性别: GG
- 专业: 燃烧学
3楼2010-01-28 10:08:24
hgye
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1406
- 帖子: 180
- 在线: 46.1小时
- 虫号: 461024
- 注册: 2007-11-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2010-01-28 17:45:52