24小时热门版块排行榜    

查看: 1219  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gynan2008

金虫 (正式写手)

[求助] abinit计算弹性常数的几个小问题

关于abinit计算二阶,三阶弹性常数的计算我有几个问题不明白,想问各位大侠一下,   
问题1:原胞中原子的位置一定,我们算A1时,原胞中施加应力方向的原子位置数据改变。[0.9,0,0;0,1,0;,0,0,1]加应变[0.02,0,0;0,0,0;0,0,0].那么计算其他方向加应变时,如A6,我的输入文件是这样的[1,0,0;0,1,0;0,0,1]加应变[0,0.02,0.02;0.02,0,0.02;0.02;0.02;0]不知道是不是正确的?
   问题2: 我们没有设置温度参数(其实我还不会),我们计算的是不是0K下的结构?
   问题3: 我看BTO的相转变图0K下不是立方相也不是四方相,有的文献说是三方相,这样计算有没有可比性?因为,根据原胞的参数,我们默认的是立方相,计算A1~A6的结果就可以了。
   问题4:x轴坐标如何确定,计算A1可以从(0.9~1.1),那么计算A6时呢?是不是用矢量叠加的方法?
问题5:我用多项式拟合的指数函数的形式,然后再用Taylor级数展开得到y=a+A1x+A2*x^2+。。。的形式严密不严密?
请各位大侠指教,坐等帮忙

公示.png



[ Last edited by gynan2008 on 2012-10-25 at 10:58 ]
回复此楼

» 猜你喜欢

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

苔花如米小,也学牡丹开!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tangosnow

木虫 (小有名气)

不明白你所说A1,A6是什么意思
4楼2012-11-01 05:48:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gynan2008 的主题更新
信息提示
请填处理意见