| 查看: 4324 | 回复: 20 | ||||
[求助]
关于弹性常数矩阵结果求助 已有3人参与
|
||||
|
VASP自从出了5.2版本以后，根据手册和网上以往的帖子，能够直接计算出弹性常数了，弹性常数以矩阵的形式表示出来。 第一次计算， 我设置 ISIF=3 IBRION=6 POTIOM=0.01 截断能取了400eV. 计算完成后在OUTCAR里的结果如下：  坛子里面有些朋友说计算得到的弹性常数中不能出现负数，于是在别人的建议下。 我做了第二次计算， 提高了截断能至600eV. POTIOM=0.1 计算完成后在OUTCAR里得到了第二次的结果如下：  对比第一次和第二的结果，可以看到，负数值没有消除，另外两次计算得到的弹性常数数值偏差很大。 我想请教坛子里的高手：第一，矩阵中的弹性常数能否为负数？ 第二，为什么提高了截断能和 POTIOM后，计算的结果偏差如此之大？ 如果稍微更改下参数，结果就不同了，后续计算的准确性也就无法保证了。 希望有“大虾“能够答疑解惑！ 不甚感谢！  |
回复此楼» 收录本帖的淘帖专辑推荐
 八龙尊者 |
» 猜你喜欢
求个博导看看
已经有17人回复
-
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
带资进组求博导收留
已经有11人回复
-
自荐读博
已经有5人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 过程装备特殊零部件应力分析
已经有144人回复
- 非线性分析中的ansys跟踪显示
已经有83人回复
- 站里关于MS Analysis计算材料的应力应变与Cubic以及Trigonal cell弹性常数的讨论
已经有34人回复
jimsmart
捐助贵宾 (职业作家)
- 应助: 109 (高中生)
- 贵宾: 0.2
- 金币: 3581.7
- 散金: 100
- 红花: 5
- 帖子: 3110
- 在线: 567.6小时
- 虫号: 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
emilyoyang
木虫 (正式写手)
- 应助: 232 (大学生)
- 金币: 4340.4
- 红花: 42
- 帖子: 667
- 在线: 189.9小时
- 虫号: 1630291
- 注册: 2012-02-20
- 性别: MM
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zyqh: 金币+5, ★有帮助 2015-12-21 15:47:36
liliangfang: 金币+1, 谢谢交流 2015-12-24 20:40:03
感谢参与，应助指数 +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
2楼2015-12-20 22:35:29
4楼2015-12-21 12:48:33
5楼2015-12-21 12:55:44
7楼2015-12-21 15:47:26
zyqh8楼2015-12-21 21:33:44
9楼2015-12-21 22:19:15
10楼2015-12-21 22:34:43