| 查看: 3802 | 回复: 7 | ||||||||
| 当前主题已经存档。 | ||||||||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 monowolf 的 1 个金币 | ||||||||
monowolf铁杆木虫 (著名写手)
|
[交流]
【求助】表面能计算
|
|||||||
|
我发表论文时表面能的计算公式是这样定义的: The surface energy is defined as: r=(Eslab(Me-Al)-n(Al)Ebulk(Al)-n (Me )Ebulk(Me) /2A (1) whereEslab(Me-Al) is the total energy of Me-Al slab, while nAl and nMe are the number of Al and doped Me atoms per computational cell, and the and are the reference bulk total energy per atom, respectively. The factor 1/2 takes into account that the top and bottom of the slab are equivalent, and the A is the area of the Me-Al surface. ( )表示下标, 计算模型为5层原子，最底下一层固定，其余四层可以自由弛豫。 审稿专家提出意见说表面能的计算公式不对： p.5 - the authors define their surface energy (equation (1)) by simply dividing the total energy difference by 2. This implies that both surfaces of the slab are equal, and that is presumably not true, sinceone side is relaxed and the other is not. A better way to calculate the surface energy would be to do acalculation on a static 5-layer slab and deduce the "static" surface energy in that manner. The energy difference in equation (1) could then be adjusted for this static surface energy, and the division by 2 could be eliminated. I wouldn't expect this modification to significantly change the calculated surface energy trends, but it might change the quantitative values to some extent. 我没有看明白，他希望我怎么来计算这个表面能，希望看明白的给我建议，看怎么计算表面，怎么修改合适。谢谢。 附件为这问题的word，主要是论坛不能写公式。 请大家帮忙 |
回复此楼» 收录本帖的淘帖专辑推荐
 零点能 | Photochemistry | 第一性原理 | 表界面计算 |
| 模拟，理论 |
» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
marklau
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 247.8
- 散金: 10
- 帖子: 122
- 在线: 37.4小时
- 虫号: 619105
- 注册: 2008-10-07
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2010-03-01 08:48:59
★
monowolf(金币+1):谢谢 2010-03-01 09:53
ice_rain(金币+1):谢谢交流讨论 2010-03-03 19:25
monowolf(金币+1):谢谢 2010-03-01 09:53
ice_rain(金币+1):谢谢交流讨论 2010-03-03 19:25
|
本帖内容被屏蔽 |
3楼2010-03-01 09:32:47
hgye
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1406
- 帖子: 180
- 在线: 46.1小时
- 虫号: 461024
- 注册: 2007-11-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★
monowolf(金币+3):谢谢 2010-03-01 09:53
monowolf(金币+3):采用你的建议 2010-03-02 11:05
ice_rain(金币+2):谢谢指导 2010-03-03 19:25
monowolf(金币+3):谢谢 2010-03-01 09:53
monowolf(金币+3):采用你的建议 2010-03-02 11:05
ice_rain(金币+2):谢谢指导 2010-03-03 19:25
|
你原来的计算肯定是不对的，因为你上下两个面一个被固定，一个允许弛豫，两个面是不等价的。专家的建议是你在做一个5层全部固定的静态计算，用你的公式（1）处理，得到没有弛豫的表面能。在你原来的计算中不要除以2，相当于得到了一个弛豫面加上一个不弛豫面的表面能。两次计算相减就可以得到一个弛豫后的表面能。 个人建议，你还可以这样做，把你建的模型中间部分固定，上下两个面都让弛豫，这样两个面就是等价的了，用你原来的方法计算就没有问题。 |
4楼2010-03-01 09:45:45
monowolf
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 贵宾: 3.05
- 金币: 5268.7
- 散金: 4010
- 红花: 2
- 帖子: 1865
- 在线: 145.1小时
- 虫号: 75148
- 注册: 2005-06-16
- 性别: GG
- 专业: 金属材料的制备科学与跨学
5楼2010-03-01 09:53:13
hgye
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1406
- 帖子: 180
- 在线: 46.1小时
- 虫号: 461024
- 注册: 2007-11-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2010-03-01 10:05:04
7楼2010-03-03 18:14:18
宏hong
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2754.9
- 散金: 305
- 红花: 1
- 沙发: 1
- 帖子: 649
- 在线: 237.6小时
- 虫号: 842750
- 注册: 2009-09-07
- 性别: GG
- 专业: 燃烧学
8楼2010-03-09 10:18:15
