| 查看: 4690 | 回复: 20 | |||
zyqh铁杆木虫 (著名写手)
|
[求助]
关于弹性常数矩阵结果求助 已有3人参与
|
||
|
VASP自从出了5.2版本以后，根据手册和网上以往的帖子，能够直接计算出弹性常数了，弹性常数以矩阵的形式表示出来。 第一次计算， 我设置 ISIF=3 IBRION=6 POTIOM=0.01 截断能取了400eV. 计算完成后在OUTCAR里的结果如下：  坛子里面有些朋友说计算得到的弹性常数中不能出现负数，于是在别人的建议下。 我做了第二次计算， 提高了截断能至600eV. POTIOM=0.1 计算完成后在OUTCAR里得到了第二次的结果如下：  对比第一次和第二的结果，可以看到，负数值没有消除，另外两次计算得到的弹性常数数值偏差很大。 我想请教坛子里的高手：第一，矩阵中的弹性常数能否为负数？ 第二，为什么提高了截断能和 POTIOM后，计算的结果偏差如此之大？ 如果稍微更改下参数，结果就不同了，后续计算的准确性也就无法保证了。 希望有“大虾“能够答疑解惑！ 不甚感谢！  |
回复此楼» 收录本帖的淘帖专辑推荐
 八龙尊者 |
» 猜你喜欢
评审感受-评审感受-评审感受
已经有6人回复
-
现在不知道怎么办，感觉很痛苦
已经有5人回复
-
导师各种操作恶心咋办
已经有9人回复
-
生物学博士
已经有6人回复
-
窗边初夏的小雨
已经有8人回复
-
护理论文 晋升
已经有3人回复
-
2026年申博-电池方向
已经有11人回复
-
26年申博自荐-计算机视觉
已经有4人回复
-
2026博士申请求助
已经有5人回复
-
研究生做的很差，你们会让毕业吗？
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 过程装备特殊零部件应力分析
已经有144人回复
- 非线性分析中的ansys跟踪显示
已经有83人回复
- 站里关于MS Analysis计算材料的应力应变与Cubic以及Trigonal cell弹性常数的讨论
已经有34人回复
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
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
zyqh
铁杆木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 6714.6
- 散金: 1454
- 红花: 22
- 帖子: 2380
- 在线: 593.5小时
- 虫号: 2762045
- 注册: 2013-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2015-12-21 12:48:33
zyqh
铁杆木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 6714.6
- 散金: 1454
- 红花: 22
- 帖子: 2380
- 在线: 593.5小时
- 虫号: 2762045
- 注册: 2013-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2015-12-21 12:55:44
zyqh
铁杆木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 6714.6
- 散金: 1454
- 红花: 22
- 帖子: 2380
- 在线: 593.5小时
- 虫号: 2762045
- 注册: 2013-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
7楼2015-12-21 15:47:26
zyqh8楼2015-12-21 21:33:44
zyqh
铁杆木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 6714.6
- 散金: 1454
- 红花: 22
- 帖子: 2380
- 在线: 593.5小时
- 虫号: 2762045
- 注册: 2013-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
9楼2015-12-21 22:19:15
10楼2015-12-21 22:34:43