24小时热门版块排行榜    

查看: 2505  |  回复: 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的回帖

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的回帖
查看全部 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的回帖

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 00:52 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
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 范进中举一文的中心思想 +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
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 时间戳又变了 +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 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 朋友圈看到的 +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 -大大大大大-
信息提示
请填处理意见