24小时热门版块排行榜    

查看: 1891  |  回复: 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):谢谢参与
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的回帖

zzjust

新虫 (小有名气)

★ ★ ★
lzb8628491(金币+1):谢谢参与
zzy870720z(金币+2):谢谢指导 2010-09-02 17:25:46
内容已删除
3楼2010-09-02 14:19:16
已阅   回复此楼   关注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的回帖

lzb8628491(金币+1):谢谢参与
5楼2010-09-02 19:02:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

s20080718

金虫 (著名写手)

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


lzb8628491(金币+1):谢谢参与
学习  顶贴  需要继续讨论   楼主问题你怎么解决的?
解决问题的关键在于大家一起用心交流!
6楼2010-09-04 09:56:34
已阅   回复此楼   关注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的回帖

鲤鱼飞

铁虫 (正式写手)

你按照NaCl的构型计算计算呢?
CAEP
8楼2010-09-04 13:44:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

truemantwo

木虫 (正式写手)

★ ★
lzb8628491(金币+1):谢谢参与
zzy870720z(金币+1):谢谢提示 2010-09-05 15:17:42
会不会是形成能的公式计算有问题吧,貌似不能简单相减,还得看你建立的TiC晶胞里面总共有多少个原子
9楼2010-09-04 16:03:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiao_ma523

铁杆木虫 (著名写手)


lzb8628491(金币+1): 谢谢参与
请教下楼主,单个Ti原子能量计算时是怎么设置的,我计算Ti的时候不收敛。。
10楼2013-10-20 18:23:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lzb8628491 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见