请教：
界面能计算公式是:
σ界面能=（Eint-ΣNiμi+PV-TS）/2A其中Ni为体系中第i个元素的数目，μi为该元素的化学势。
我用的是castep,想问一下，原子的数目是从所建的模型中数出来？还是从Filter中的symmetry system中得到？
还有化学势怎样得到？是查物理化学数据库查的还是castep计算的？如果是计算该如何计算？

请高手指教！多谢了！

[ Last edited by 干冰388 on 2008-5-28 at 09:06 ]

晶体的界面能是指单位面积界面上的能量与无界面时该区域能量之差。在CASTEP中如果用你提供的那个公式肯定是不现实的。我简要介绍一下我们怎么计算的界面能。首先构造超晶胞，超晶胞由A,B两种材料组成，各含m,n个晶胞。超晶胞两侧都有真空层，A,B相连的就是要研究的界面。其次，要优化超晶胞计算得到其总能量E1.。这个超晶胞的能量E1由五项组成：m个A的晶胞能量(体材料计算），n个B的晶胞能量，A材料的一个表面能Ea,B材料的一个表面能Eb,还有就是要计算的界面能了。减去前面四个能量就得到界面能。

形成界面模型的时候是先分别切出两种材料的某一晶面，再形成界面加真空层的这种模型。
1.在这种情况下两种材料中的晶胞个数如何确定？
2.另外，在形成切面时必然会有悬挂键出现，对于界面模型中上下两个表面上的悬挂键你们是如何处理的？我们这的一个老师说需要用H来终结这悬挂键，如果是这样的话，这部分的能量又该如何考虑？

请教一下你们对这方面的问题是如何考虑的！
非常感谢！