24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1773  |  回复: 1

但丁

木虫 (小有名气)

[求助] vasp计算弹性常数晶胞选择已有1人参与

我算单斜晶系弹性常数,单胞(88原子)和原胞(44原子)算出来的弹性常数Cij差别很大,但是计算出的B,G,E都比较接近,所以想请教下:
1 两个结果都可以用么,还是单胞原胞之间Cij可以用公式换算
2 我单胞88原子,选用4 8 4 K点,K点数量30个左右,32核服务器大概要跑40天,这样的精确度是否合适
3 我用弹性常数结果进行buckingham两体势进行势函数拟合,GULP拟合的时候都是用原胞来进行拟合的,但是我用的是单胞弹性常数的结果,这样是否正确

单胞(GPa)
Bh=139     Gh=73     Eh=187
原胞(GPa)
Bh=144     Gh=76     Eh=194
C11kBar                C22kBar            C33kBar     C44kBar        C55kBar         C66kBar          C12kBar         C13kBar                C23kBar                C15kBar        C25kBar        C35kBar        C46kBar
1749.2206        2255.9654  2927.5723  906.8144        927.8352  693.7768        868.6353        1130.4024        1049.7312        255.4382        -1.1155        14.5271        -46.9599   (单胞)
2044.4453        2295.7081  2837.4896  952.3175        928.1164  620.9653        821.9845        1111.5495        1095.0755        81.6558        52.9417        23.5467        16.9443     (原胞)
回复此楼
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

Fly330

银虫 (初入文坛)

【答案】应助回帖


感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-10-12 09:33:39
原胞和单胞对称性不一样,所以可能Cij会有不同,但是弹性常数是统计性质能表现一致性。这样去计算下结构的状态方程EOS,看看算出来的体弹模量B是否与前面一致。一直的话结果还是能用的。

发自小木虫IOS客户端
2t2l
2楼2015-10-11 12:27:29
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 但丁 的主题更新
信息提示
请填处理意见