24小时热门版块排行榜    

CyRhmU.jpeg
查看: 5962  |  回复: 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 的主题更新
信息提示
请填处理意见