24小时热门版块排行榜    

查看: 7145  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 师弟论文见刊大半年才想起来申请专利,还能抢救一下吗? +5 13108017953 2026-09-21 5/250 2026-09-27 16:46 by Leogzhya
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +3 zju2000 2026-09-26 3/150 2026-09-27 16:26 by lala1008
[硕博家园] 上海工程技术大学 激光智能制造课题组 2027级博士研究生招生 +5 oxidpl 2026-09-20 7/350 2026-09-27 10:43 by Equinoxhua
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:42 by 4IW0sJvtEqX8
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:26 by 4IW0sJvtEqX8
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:15 by 4IW0sJvtEqX8
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:03 by 4IW0sJvtEqX8
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:44 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 GxESXJmptOtQ 2026-09-25 3/150 2026-09-27 06:32 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:08 by XgCC8uTcwILl
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:05 by XgCC8uTcwILl
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:56 by XgCC8uTcwILl
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:56 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:53 by XgCC8uTcwILl
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:50 by XgCC8uTcwILl
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:44 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 6/300 2026-09-27 05:42 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 5/250 2026-09-25 11:19 by yhq9807
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 7/350 2026-09-25 11:18 by yhq9807
[教师之家] 课题组招2027级博士 上海工程技术大学 激光智能制造方向 +4 oxidpl 2026-09-20 6/300 2026-09-24 20:37 by oxidpl
信息提示
请填处理意见