| 查看: 4583 | 回复: 20 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
关于弹性常数矩阵结果求助 已有3人参与
|
|||
|
VASP自从出了5.2版本以后，根据手册和网上以往的帖子，能够直接计算出弹性常数了，弹性常数以矩阵的形式表示出来。 第一次计算， 我设置 ISIF=3 IBRION=6 POTIOM=0.01 截断能取了400eV. 计算完成后在OUTCAR里的结果如下：  坛子里面有些朋友说计算得到的弹性常数中不能出现负数，于是在别人的建议下。 我做了第二次计算， 提高了截断能至600eV. POTIOM=0.1 计算完成后在OUTCAR里得到了第二次的结果如下：  对比第一次和第二的结果，可以看到，负数值没有消除，另外两次计算得到的弹性常数数值偏差很大。 我想请教坛子里的高手：第一，矩阵中的弹性常数能否为负数？ 第二，为什么提高了截断能和 POTIOM后，计算的结果偏差如此之大？ 如果稍微更改下参数，结果就不同了，后续计算的准确性也就无法保证了。 希望有“大虾“能够答疑解惑！ 不甚感谢！  |
回复此楼» 猜你喜欢
313求调剂
已经有4人回复
-
材料调剂
已经有13人回复
-
材料专硕283求调剂
已经有13人回复
-
295求调剂
已经有7人回复
-
266分，一志愿电气工程，本科材料，求材料专业调剂
已经有7人回复
-
307分材料专业求调剂
已经有5人回复
-
一志愿北交大材料工程总分358求调剂
已经有5人回复
-
求调剂
已经有4人回复
-
材料与化工371求调剂
已经有11人回复
-
材料求调剂
已经有4人回复
emilyoyang
木虫 (正式写手)
- 应助: 232 (大学生)
- 金币: 4340.4
- 红花: 42
- 帖子: 667
- 在线: 189.9小时
- 虫号: 1630291
- 注册: 2012-02-20
- 性别: MM
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★
zyqh: 金币+5, ★有帮助, 谢谢！ 2015-12-22 16:05:03
liliangfang: 金币+1, 谢谢交流 2015-12-24 20:40:16
zyqh: 金币+5, ★有帮助, 谢谢！ 2015-12-22 16:05:03
liliangfang: 金币+1, 谢谢交流 2015-12-24 20:40:16
|
Read the IBRION tag carefully. It is of great help to understand how VASP derive elastic tensor and then you will figure out. I agree with zyl雨田's statement that you need to do elastic stability analysis before you make a judgement on your results. Please refers to this paper. http://arxiv.org/abs/1410.0065 |
15楼2015-12-22 01:41:09
2楼2015-12-20 22:35:29
jimsmart
捐助贵宾 (职业作家)
- 应助: 109 (高中生)
- 贵宾: 0.2
- 金币: 3583.7
- 散金: 100
- 红花: 5
- 帖子: 3110
- 在线: 568.1小时
- 虫号: 57301
- 注册: 2005-01-29
- 性别: GG
- 专业: 光谱分析
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zyqh: 金币+5, ★有帮助 2015-12-21 12:55:53
liliangfang: 金币+1, 谢谢交流 2015-12-24 20:39:57
感谢参与，应助指数 +1
zyqh: 金币+5, ★有帮助 2015-12-21 12:55:53
liliangfang: 金币+1, 谢谢交流 2015-12-24 20:39:57
|
你算出的结果,是一种状态的显示,既然算出来了结果,就有一定的合理性; 关键是你怎么去解释它和解读它所带来的信息; 你改变了参数,结果发生变化也是自然而然的事情,没有变化那才叫奇怪; 弹性常数能否为负数, 要用发展的眼光看. 如果别人说可以或不可以,只是一家之辞; 若错了,你会怪他,可那又何必呢? 若对了,也不过是适合了一个问题,对提高你的计算水平又没有太大的帮助. |
3楼2015-12-21 08:32:02
4楼2015-12-21 12:48:33
