| 查看: 1512 | 回复: 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人回复
hgye
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1406
- 帖子: 180
- 在线: 46.1小时
- 虫号: 461024
- 注册: 2007-11-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2010-01-28 17:45:52
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
宏hong
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2754.9
- 散金: 305
- 红花: 1
- 沙发: 1
- 帖子: 649
- 在线: 237.6小时
- 虫号: 842750
- 注册: 2009-09-07
- 性别: GG
- 专业: 燃烧学
|
恩，你说的应该是这个公式吧： 构造一个slab，假设n个原子，计算它的能量E（n），同时计算该材料块体时单个原子对应的能量E'，表面能为 S＝（E（n）－n×Ebulk）/2s, s为slab一面的面积。 按照这个公式，Ebulk为体相能量，即从自带的structure文件中导入优化单胞后得到的数值，为Ebulk=-6365.716872Ha，其中单胞含有8个原子，切面五层超晶胞结构优化后含20个原子，5层结构优化后得出的total energy为E=-18795.187382Ha，带入公式如下： S=-18795.187382+（20/8）*6265.716872.......... 后面的我就不写了（不好写...），这样计算公式括号内的减法部分还是负值（-3130.8952），请问我这样计算错在哪儿了？是不是我对这个公式的理解错了？ 十分感谢！ |
5楼2010-01-28 20:30:45