24小时热门版块排行榜    

查看: 2504  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 00:40 by OEbVnUOu01ol
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 00:28 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 8/400 2026-08-22 23:51 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 10/500 2026-08-22 23:42 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 3/150 2026-08-22 22:03 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-22 21:54 by OEbVnUOu01ol
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[教师之家] 跳槽后在研项目怎么办? +3 简单化xn 2026-08-22 7/350 2026-08-22 19:46 by 简单化xn
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见