24小时热门版块排行榜    

查看: 647  |  回复: 3
当前主题已经存档。

doudouding

木虫 (小有名气)

[交流] 【求助】castep和gulp计算出来的能量分别是什么

用castep和gulp计算的MgO的能量结果分别是
castep:
Final energy, E             =  -456.4272031579     eV
Final free energy (E-TS)    =  -456.4272031579     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)      =  -456.4272031579     eV

gulp:
Components of energy :

--------------------------------------------------------------------------------
  Interatomic potentials     =           7.17772161 eV
  Monopole - monopole (real) =          -3.84680224 eV
  Monopole - monopole (recip)=         -44.08559419 eV
  Monopole - monopole (total)=         -47.93239642 eV
--------------------------------------------------------------------------------
  Total lattice energy :
    Primitive unit cell      =         -40.75467482 eV
    Non-primitive unit cell  =        -163.01869927 eV
--------------------------------------------------------------------------------


castep计算结果也是一个Primitive unit cell的,这两个有什么不同? manual里都说是单点能的计算,结果怎么会差别怎么大。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjmaxpayne

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
qasd(金币+1,VIP+0):谢谢~ 11-2 19:48
能量的绝对值没有太大意义,取决于不同的方法还有不同的基组选择等。
2楼2009-11-02 15:28:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

doudouding

木虫 (小有名气)

我是想知道这两个能量从物理概念上有什么差别

如果我想用castep来计算Total lattice energy 该怎么做呢?
3楼2009-11-02 15:45:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjmaxpayne

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
freshgirl(金币+8,VIP+0):非常详细,谢谢分享! 11-4 22:22
恩,这样的话,你可能需要计算晶体中各个离子的电荷,比如,可以通过CASTEP进行布局数分布,然后用:

Born-Landé 方程



进行计算

公式中各个符号的含义参见:http://en.wikipedia.org/wiki/Born-Land%C3%A9_equation

n是Born Exponent,一般是和固体的可压缩性有关,具体如何从理论来计算,参考:
Cotton, F. Albert; Wilkinson, Geoffrey; (1966). Advanced Inorganic Chemistry (2d Edn.) New York:Wiley-Interscience.

[ Last edited by yjmaxpayne on 2009-11-3 at 10:17 ]
4楼2009-11-03 10:15:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 doudouding 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 9/450 2026-08-01 18:43 by Kddn1pB6pAzX
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 9zw1gfyi7x 2026-08-01 7/350 2026-08-01 18:42 by Kddn1pB6pAzX
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 9zw1gfyi7x 2026-08-01 6/300 2026-08-01 18:37 by Kddn1pB6pAzX
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-08-01 8/400 2026-08-01 18:31 by Kddn1pB6pAzX
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 6/300 2026-08-01 18:31 by Kddn1pB6pAzX
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:23 by Kddn1pB6pAzX
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 8/400 2026-08-01 18:22 by Kddn1pB6pAzX
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 9zw1gfyi7x 2026-07-31 11/550 2026-08-01 18:22 by Kddn1pB6pAzX
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-07-31 10/500 2026-08-01 18:20 by Kddn1pB6pAzX
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:11 by Kddn1pB6pAzX
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:10 by Kddn1pB6pAzX
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:10 by Kddn1pB6pAzX
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 9zw1gfyi7x 2026-07-31 10/500 2026-08-01 18:03 by Kddn1pB6pAzX
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 8/400 2026-08-01 18:02 by Kddn1pB6pAzX
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 jzg3xk1po8 2026-07-31 7/350 2026-08-01 17:51 by Kddn1pB6pAzX
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 jzg3xk1po8 2026-07-31 12/600 2026-08-01 17:50 by Kddn1pB6pAzX
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 jzg3xk1po8 2026-07-31 9/450 2026-08-01 17:43 by Kddn1pB6pAzX
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +16 低垂的野花 2026-07-31 22/1100 2026-08-01 16:39 by Cathyau
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 8/400 2026-07-31 18:13 by beefly
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
信息提示
请填处理意见