| 查看: 3499 | 回复: 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 | 第一性原理 | 表界面计算 |
| 模拟，理论 |
» 猜你喜欢
268求调剂
已经有8人回复
-
26申博
已经有4人回复
-
211本，11408一志愿中科院277分，曾在中科院自动化所实习
已经有7人回复
-
求调剂，总分315，考的生物医药，一志愿湖南师范大学。调剂到任何专业都可以
已经有3人回复
-
材料专硕306英一数二
已经有7人回复
-
材料工程专硕274一志愿211求调剂
已经有6人回复
-
一志愿南京大学，080500材料科学与工程，调剂
已经有4人回复
-
332求调剂
已经有5人回复
-
271求调剂
已经有14人回复
-
275求调剂
已经有4人回复
★
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
marklau
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 247.8
- 散金: 10
- 帖子: 122
- 在线: 37.4小时
- 虫号: 619105
- 注册: 2008-10-07
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2010-03-01 08:48:59
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
50