【求助】MS表面能的计算
求助：计算铝001晶面的表面能
计算中表面能的计算公式为
γ=[Eslab-Ebulk] /2A
其中,A表示表面积的大小,Eslab表示模型能量值,Ebulk表示
晶体内部原子数目相同时系统的能量值.
通过建立铝晶胞（4个铝原子）模型结果优化得到E=-228.79ev
然后建立001表面，其中depth为3.5，设置真空层17.98，优化后得到
E=-399.36ev，共七个铝原子（不知道对不对？？）
用ms中的tools-atom volumes&surfaces计算得到表面积A为18.39 (A-2)
问题是：
1。Ebulk是不是就是铝晶胞的能量，还是乘7/4？？或者是其它的算法？
如何得到晶体内部原子数目相同时的Ebulk呢？
2.这样算的表面积对不对，还需要除2吗？应该如何算表面积？
希望大家可以给出详细的计算过程和解答，参考文献的值是0.92 J/m2，非常感谢！
[ Last edited by freshgirl on 2009-6-26 at 16:44 ]
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
A 应该是 Slab 的 a*b
Ebulk是晶胞能量乘7/4
表面积是slab的a*b
可以试着算一下
试算一下：γ(Al-001)=[Eslab-Ebulk] /2A=[-399.36-(-228.79*7/4)]ev/18.39 A2=0.8896J/m2,与文献值有较好吻合
1eV=1.602x10^(-19)J
另询：建立表面的depth与设置真空层设置有什么依据？是由小到大设置不同值看每原子平均能量在哪个临界尺寸上稳定（i.e.超过该值能量的每原子平均结果基本不变，depth不同时表面所含原子数目也不同）吗，
γ(Al-001)=[Eslab-Ebulk] /2A，是能量差除以表面积，还是除以2倍的表面积？
楼主，4个铝原子是怎么取得？从ms导入后如何改？
感谢大家的回答，我算了一下如果表面积设为a*b,应该是驰豫前的表面积，驰豫后的我觉得应该用MS中计算出来的表面积，这时候算表面能就不用*2，不知道这样理解对不对，请大家指教。谢谢！
附回楼上的，铝的体晶胞实际上就是4个原子的。
晶体内部原子数目相同时系统的能量值.具体是什么意思?
[ Last edited by wangqj1 on 2008-6-30 at 20:41 ]