24小时热门版块排行榜    

查看: 7324  |  回复: 13

nwpu2010

铁虫 (初入文坛)

[求助] CASTEP计算弹性常数C11,C12,C44严重偏大

求助高手,最近计算金属间化合物Co基合金,计算出来的弹性常数普遍偏大,特别是C11,请教一下,计算弹性常数是与那些因数有关,如何才能减小这些常数,特别是C11值,下面是我计算的结果以及别人计算的结果与实验值。Co是强磁性金属,考了其自旋极化。       
         This work                Experiment        Other works
                  LDA                              
ɑ0(Å        3.5315                    ɑ3.599          3.571;3.582
C11        394                    271              363.39;264
C12        212                    172               189.92;162
C44        209                    162               211.64;153
另外,我计算的2X2X2的超胞结构,计算的结构与实验相差太远,希望高手能帮忙解释一下可能的原因!
回复此楼

» 收录本帖的淘帖专辑推荐

仿真建模与计算 动力学强帖

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

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

不要太累
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
zzy870720z(金币+2): 谢谢指教 2011-12-18 09:17:06
计算的精度足够了么?CASTEP中弹性常数的计算是strain-stress方法计算的,stress对于计算精度及收敛条件都很敏感,如果精度或者前期收敛做的不好,结果偏差很大是正常的。
至于对付金属,LDA一般够用,并且,不管是LDA还是GGA,偏差到50%,实在有点说不过去~
思想重于技巧,内涵重于表象
5楼2011-12-16 19:32:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nwpu2010

铁虫 (初入文坛)

引用回帖:
楼: Originally posted by 06022045 at 2011-12-16 15:50:09:
如果你算得没错的话,在LDA中换换不同的赝势。 最好用GGa算一下,这样己验证了你LDA是否正确,有可以回答审稿人,算得结果就是这样。

非常谢谢你的回答,我再试试!
不要太累
8楼2011-12-17 14:16:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

06022045

金虫 (小有名气)

【答案】应助回帖

★
感谢参与,应助指数 +1
uuv2010(金币+1): 多谢指导! 2011-12-16 09:46:45
用GGA在算一下,LDA一般都要大一些。而且有时候理论核试验有差距也是正常的。理论毕竟是理想条件的
2楼2011-12-16 08:26:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nwpu2010

铁虫 (初入文坛)

引用回帖:
楼: Originally posted by 06022045 at 2011-12-16 08:26:03:
用GGA在算一下,LDA一般都要大一些。而且有时候理论核试验有差距也是正常的。理论毕竟是理想条件的

谢谢你的回复,我用GGA也算了,加上自旋极化他就不收敛了,体系是超胞,32个原子,这个算出来的结果是在是很纠结,有实验对比,审稿人直接就是与实验相差太大!
不要太累
3楼2011-12-16 11:04:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

06022045

金虫 (小有名气)

【答案】应助回帖

★
zzy870720z(金币+1): 鼓励交流 2011-12-18 09:16:51
如果你算得没错的话,在LDA中换换不同的赝势。 最好用GGa算一下,这样己验证了你LDA是否正确,有可以回答审稿人,算得结果就是这样。
4楼2011-12-16 15:50:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sars518

木虫 (著名写手)

小木虫灌水大队队长

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
zzy870720z(金币+1): 鼓励交流 2011-12-18 09:17:33
fzx2008(金币+1): 谢谢指教 2011-12-18 13:17:37
你计算的晶格常数 多大? 与试验比较 差距如何?
你的Ecut Kpoints 取多少? 收敛精度?
你用的是 UltraSoft 还是 Norm-Conserving ?
这些都是很关键的因素。Castep cij 对 Ecut和KP 特别敏感
人生是一场旅行,在欣赏风景的同时,别忘了你也是别人的风景。
6楼2011-12-16 20:17:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nwpu2010

铁虫 (初入文坛)

内容已删除
不要太累
7楼2011-12-17 14:13:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nwpu2010

铁虫 (初入文坛)

希望知道的进来讨论讨论,别让帖子沉了!
不要太累
9楼2011-12-17 14:19:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夭夭小竹

新虫 (初入文坛)

同问同问……我现在还遇见了个更郁闷的问题,E-cut设置改变k-point自动会变……我都要崩溃了……眼看过几天就中期了,我还什么都没做出来……挂了我这四年就要变五年了啊……
10楼2012-05-04 22:39:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nwpu2010 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 SX5mJmPWDfOO 2026-11-13 4/200 2026-11-15 07:34 by xSrQNSSjcgoX
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 XJXMuRjD7cK7 2026-11-13 4/200 2026-11-15 06:50 by xSrQNSSjcgoX
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-15 05:33 by xSrQNSSjcgoX
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 uKOSCocCb6Fm 2026-11-13 3/150 2026-11-15 04:38 by xSrQNSSjcgoX
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 uKOSCocCb6Fm 2026-11-12 7/350 2026-11-15 03:10 by xSrQNSSjcgoX
[文学芳草园] 厉娜 +4 张丹枫1991 2026-11-09 6/300 2026-11-15 00:48 by Equinoxhua
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 FChKX8meu1ds 2026-11-09 8/400 2026-11-15 00:35 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 YHiw1gou1Eda 2026-11-10 8/400 2026-11-14 22:51 by Equinoxhua
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 YHiw1gou1Eda 2026-11-10 8/400 2026-11-14 22:20 by Equinoxhua
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 FChKX8meu1ds 2026-11-09 8/400 2026-11-14 22:04 by Equinoxhua
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 G2rhSvqvhNzl 2026-11-10 7/350 2026-11-14 21:58 by Equinoxhua
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-14 16:09 by KXypaGgBBlw7
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-13 11:37 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 Z7qfCsdIiD1T 2026-11-11 5/250 2026-11-12 11:44 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 czAiolzR0IJC 2026-11-10 9/450 2026-11-12 11:00 by HajghKOMynmz
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 10:49 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 09:46 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 FChKX8meu1ds 2026-11-09 12/600 2026-11-12 09:43 by HajghKOMynmz
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 05:19 by 9fOTqj074tIg
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 0TEMgxdT1eIw 2026-11-09 8/400 2026-11-12 05:00 by 9fOTqj074tIg
信息提示
请填处理意见