| 查看: 1840 | 回复: 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 很显然，算出来的结果是负的，有没有哪位大侠有相关经验的，或者知道哪错的？十分感谢了 |
回复此楼» 猜你喜欢
浙江师范大学是怎么坑我的
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
求教各位大神：2026教育部人文社科青年基金项目何时公示呀？
已经有6人回复
-
面上有专家说实验设备不是我们单位的
已经有8人回复
-
科研求助
已经有4人回复
-
博士去二本高校当辅导员，如何调整心态？
已经有11人回复
-
江苏省自然基金 什么时候出结果
已经有6人回复
-
闲聊
已经有3人回复
宏hong
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2754.9
- 散金: 305
- 红花: 1
- 沙发: 1
- 帖子: 649
- 在线: 237.6小时
- 虫号: 842750
- 注册: 2009-09-07
- 性别: GG
- 专业: 燃烧学
3楼2010-01-28 10:08:24
lfhuang
木虫 (著名写手)
- 应助: 39 (小学生)
- 金币: 3649.6
- 散金: 2321
- 红花: 37
- 帖子: 1432
- 在线: 681.5小时
- 虫号: 505824
- 注册: 2008-02-17
- 性别: GG
- 专业: 金属材料表面科学与工程
2楼2010-01-28 10:04:28
hgye
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1406
- 帖子: 180
- 在线: 46.1小时
- 虫号: 461024
- 注册: 2007-11-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2010-01-28 17:45:52
宏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