24小时热门版块排行榜    

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

zyqh

铁杆木虫 (著名写手)

[求助] 关于弹性常数矩阵结果求助 已有3人参与

VASP自从出了5.2版本以后,根据手册和网上以往的帖子,能够直接计算出弹性常数了,弹性常数以矩阵的形式表示出来。
第一次计算, 我设置 ISIF=3 IBRION=6 POTIOM=0.01 截断能取了400eV.
计算完成后在OUTCAR里的结果如下:
关于弹性常数矩阵结果求助

坛子里面有些朋友说计算得到的弹性常数中不能出现负数,于是在别人的建议下。
我做了第二次计算, 提高了截断能至600eV. POTIOM=0.1
计算完成后在OUTCAR里得到了第二次的结果如下:
关于弹性常数矩阵结果求助-1


对比第一次和第二的结果,可以看到,负数值没有消除,另外两次计算得到的弹性常数数值偏差很大。
我想请教坛子里的高手:第一,矩阵中的弹性常数能否为负数?
                                       第二,为什么提高了截断能和 POTIOM后,计算的结果偏差如此之大?
如果稍微更改下参数,结果就不同了,后续计算的准确性也就无法保证了。
希望有“大虾“能够答疑解惑!
                                     不甚感谢!
回复此楼
雨后方能见彩虹
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zyqh: 金币+5, 有帮助 2015-12-21 15:47:36
liliangfang: 金币+1, 谢谢交流 2015-12-24 20:40:03
Try to reduce the POTIM. e.g POTIM = 0.005

POTIM refers to the ion displacement step. If you choose too large step (POTIM), a large stressing force will be applied to the crystal structure, and accordingly, the crystal structure will be broken. That is why you get negative values in tensor matrix. Note that, a negative value in tensor matrix does not always indicate you get the wrong results. In some cases (low symmetry), the negative value in tensor matrix will show up. Keep in mind that, the main tensor should never be negative.(Cii, i=1...6).
6楼2015-12-21 15:07:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 21 个回答

tianlan9057

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-12-21 08:01:54
zyqh: 金币+5, 有帮助 2015-12-21 12:48:47
1、你算的是什么材料的弹性常数矩阵啊?不同对称性的材料,矩阵不同
2、弹性常数矩阵的一些分量可以为负值,但XX分量必为正!
2楼2015-12-20 22:35:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jimsmart

捐助贵宾 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zyqh: 金币+5, 有帮助 2015-12-21 12:55:53
liliangfang: 金币+1, 谢谢交流 2015-12-24 20:39:57
你算出的结果,是一种状态的显示,既然算出来了结果,就有一定的合理性;
关键是你怎么去解释它和解读它所带来的信息;
你改变了参数,结果发生变化也是自然而然的事情,没有变化那才叫奇怪;
弹性常数能否为负数, 要用发展的眼光看.
如果别人说可以或不可以,只是一家之辞;
若错了,你会怪他,可那又何必呢?
若对了,也不过是适合了一个问题,对提高你的计算水平又没有太大的帮助.
3楼2015-12-21 08:32:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyqh

铁杆木虫 (著名写手)

引用回帖:
2楼: Originally posted by tianlan9057 at 2015-12-20 22:35:29
1、你算的是什么材料的弹性常数矩阵啊?不同对称性的材料,矩阵不同
2、弹性常数矩阵的一些分量可以为负值,但XX分量必为正!

首先,很感谢您的回答,我算的是钙钛矿 三斜(Tric)晶系的。
我重文献上了解了,不同相,矩阵是不一样的。
但是我不明白我稍微改动一点参数,为什么弹性常数矩阵的某些数值变化会很大?
雨后方能见彩虹
4楼2015-12-21 12:48:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见