| 查看: 5042 | 回复: 20 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
zyqh铁杆木虫 (著名写手)
|
[求助]
关于弹性常数矩阵结果求助 已有3人参与
|
|
|
VASP自从出了5.2版本以后，根据手册和网上以往的帖子，能够直接计算出弹性常数了，弹性常数以矩阵的形式表示出来。 第一次计算， 我设置 ISIF=3 IBRION=6 POTIOM=0.01 截断能取了400eV. 计算完成后在OUTCAR里的结果如下：  坛子里面有些朋友说计算得到的弹性常数中不能出现负数，于是在别人的建议下。 我做了第二次计算， 提高了截断能至600eV. POTIOM=0.1 计算完成后在OUTCAR里得到了第二次的结果如下：  对比第一次和第二的结果，可以看到，负数值没有消除，另外两次计算得到的弹性常数数值偏差很大。 我想请教坛子里的高手：第一，矩阵中的弹性常数能否为负数？ 第二，为什么提高了截断能和 POTIOM后，计算的结果偏差如此之大？ 如果稍微更改下参数，结果就不同了，后续计算的准确性也就无法保证了。 希望有“大虾“能够答疑解惑！ 不甚感谢！  |
回复此楼» 猜你喜欢
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有17人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
今天务委会开完了，明天出结果吗
已经有21人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
zyqh
铁杆木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 6735.4
- 散金: 1454
- 红花: 22
- 帖子: 2385
- 在线: 593.5小时
- 虫号: 2762045
- 注册: 2013-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
18楼2015-12-23 17:14:41
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
zyqh
铁杆木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 6735.4
- 散金: 1454
- 红花: 22
- 帖子: 2385
- 在线: 593.5小时
- 虫号: 2762045
- 注册: 2013-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2015-12-21 12:48:33
20