| 查看: 5021 | 回复: 20 | |||
zyqh铁杆木虫 (著名写手)
|
[求助]
关于弹性常数矩阵结果求助 已有3人参与
|
||
|
VASP自从出了5.2版本以后，根据手册和网上以往的帖子，能够直接计算出弹性常数了，弹性常数以矩阵的形式表示出来。 第一次计算， 我设置 ISIF=3 IBRION=6 POTIOM=0.01 截断能取了400eV. 计算完成后在OUTCAR里的结果如下：  坛子里面有些朋友说计算得到的弹性常数中不能出现负数，于是在别人的建议下。 我做了第二次计算， 提高了截断能至600eV. POTIOM=0.1 计算完成后在OUTCAR里得到了第二次的结果如下：  对比第一次和第二的结果，可以看到，负数值没有消除，另外两次计算得到的弹性常数数值偏差很大。 我想请教坛子里的高手：第一，矩阵中的弹性常数能否为负数？ 第二，为什么提高了截断能和 POTIOM后，计算的结果偏差如此之大？ 如果稍微更改下参数，结果就不同了，后续计算的准确性也就无法保证了。 希望有“大虾“能够答疑解惑！ 不甚感谢！  |
回复此楼» 收录本帖的淘帖专辑推荐
 八龙尊者 |
» 猜你喜欢
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有3人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
能否退出参与的面上项目解除限项
已经有17人回复
-
估计是周四
已经有4人回复
-
我面上完蛋了
已经有9人回复
-
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
让我中一个面上吧！
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 过程装备特殊零部件应力分析
已经有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 (小学生)
- 金币: 6730.4
- 散金: 1454
- 红花: 22
- 帖子: 2384
- 在线: 593.5小时
- 虫号: 2762045
- 注册: 2013-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2015-12-21 12:48:33
zyqh
铁杆木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 6730.4
- 散金: 1454
- 红花: 22
- 帖子: 2384
- 在线: 593.5小时
- 虫号: 2762045
- 注册: 2013-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2015-12-21 12:55:44
zyqh
铁杆木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 6730.4
- 散金: 1454
- 红花: 22
- 帖子: 2384
- 在线: 593.5小时
- 虫号: 2762045
- 注册: 2013-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
7楼2015-12-21 15:47:26
zyqh8楼2015-12-21 21:33:44
zyqh
铁杆木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 6730.4
- 散金: 1454
- 红花: 22
- 帖子: 2384
- 在线: 593.5小时
- 虫号: 2762045
- 注册: 2013-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
9楼2015-12-21 22:19:15
10楼2015-12-21 22:34:43
20