24小时热门版块排行榜    

查看: 1901  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lzb8628491

铜虫 (小有名气)

[交流] 【求助】castep计算TIC的形成能

为了求TIC的形成能,做了计算,但是计算结果和文献相差很大,请大家指点那里算错了
。计算过程如下:
1、计算单个碳原子能量
建立立个10*10*10的简单立方晶胞,在(0.5,0.5,0.5)处添加一个C原子,做能量计算。计算设置为 GGA-PBE, the plane-wave cut off 350 eV, K point [8*8*8],其他保持默认。计算的到 C原子能量为-145.84eV。
2、计算单个钛原子能量
   建立立个10*10*10的简单立方晶胞,在(0.5,0.5,0.5)处添加一个C原子,做能量计算。计算设置为 GGA-PBE, the plane-wave cut off 350 eV/400eV, K point [8*8*8],其他保持默认。
      计算得到Ti原子能量为-1375.68eV。
3、计算TiC 晶胞的能量
建立晶胞 (FM-3m, a=4.3, Ti(0,0,0), C(0.5,0.5,0.5)),进行几何优化,勾选优化晶胞项目GGA-PBE, the plane-wave cut off 350 eV, K point [8*8*8],其他保持默认。得到的能量为-1759.78eV。
4、计算形成能
1759.78-1375.68-145.84=238.26。
而文献报道为7点几,



计算结果单位就是ev,这个没有问题。TIC计算的时候使用粗胞,也就是一个C原子,一个Ti原子。另外我用类似的方法计算金属AL的形成能结果和文献报道的很接近。


我比较怀疑是我计算TiC时不对,有没有人优化过TiC的?优化这个结构是很快的。

[ Last edited by lzb8628491 on 2010-9-2 at 16:49 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

鲤鱼飞

铁虫 (正式写手)


lzb8628491(金币+1):谢谢参与
支持。形成能是用得到的东东。
另外,请问楼主。TiC为FCC结构,类似NaCl的构型,你这晶包(FM-3m, a=4.3, Ti(0,0,0), C(0.5,0.5,0.5)),咋感觉是C在体心的位置,Ti占据顶角位置,是简单立方结构哦?
CAEP
7楼2010-09-04 13:44:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

水座

铁杆木虫 (小有名气)

★ ★
lzb8628491(金币+1):谢谢参与
zzy870720z(金币+1):谢谢提示 2010-09-02 17:25:37
会不会是单位不一样,计算得到的单位是Ha,而文献报道的单位是eV,而1Ha=27.2eV,将238.26Ha/27.2=8.76eV,这就与文献报道的数据较接近。
2楼2010-09-02 12:58:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lzl8181

铁杆木虫 (著名写手)

单身浏览汉

★ ★
lzb8628491(金币+1):谢谢参与
zzy870720z(金币+1):谢谢交流 2010-09-02 20:39:34
lzb8628491(金币+1):谢谢有价值的讨论 2010-09-03 15:01:37
C的能量应该是对的,我算过石墨的每个能量为-621.7/4;金刚石的-310.8/2都比较接近楼主的,因为形成了稳定结构所以比单个C的更低一些。和楼主的想法一样估计是TiC有问题。
    计算时10*10*10的单原子大小和晶胞a=4.3不一致是否有影响?

[ Last edited by lzl8181 on 2010-9-2 at 18:53 ]
版权乜有,盗用不究!
4楼2010-09-02 18:42:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

s20080718

金虫 (著名写手)

发现问题、分析问题、解决问题


lzb8628491(金币+1):谢谢参与
学习  顶贴  需要继续讨论   楼主问题你怎么解决的?
解决问题的关键在于大家一起用心交流!
6楼2010-09-04 09:56:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见