| 查看: 5152 | 回复: 20 | |||
zyqh铁杆木虫 (著名写手)
|
[求助]
关于弹性常数矩阵结果求助 已有3人参与
|
||
|
VASP自从出了5.2版本以后，根据手册和网上以往的帖子，能够直接计算出弹性常数了，弹性常数以矩阵的形式表示出来。 第一次计算， 我设置 ISIF=3 IBRION=6 POTIOM=0.01 截断能取了400eV. 计算完成后在OUTCAR里的结果如下：  坛子里面有些朋友说计算得到的弹性常数中不能出现负数，于是在别人的建议下。 我做了第二次计算， 提高了截断能至600eV. POTIOM=0.1 计算完成后在OUTCAR里得到了第二次的结果如下：  对比第一次和第二的结果，可以看到，负数值没有消除，另外两次计算得到的弹性常数数值偏差很大。 我想请教坛子里的高手：第一，矩阵中的弹性常数能否为负数？ 第二，为什么提高了截断能和 POTIOM后，计算的结果偏差如此之大？ 如果稍微更改下参数，结果就不同了，后续计算的准确性也就无法保证了。 希望有“大虾“能够答疑解惑！ 不甚感谢！  |
回复此楼» 收录本帖的淘帖专辑推荐
 八龙尊者 |
» 猜你喜欢
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有13人回复
-
国社科系统bug了，是不是要放榜了？
已经有10人回复
-
现代”学阀”该如何界定
已经有9人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有3人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有11人回复
-
上海工程技术大学激光智能制造课题组-博士研究生
已经有10人回复
-
申博发邮件
已经有11人回复
-
我的奶奶
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 过程装备特殊零部件应力分析
已经有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 (小学生)
- 金币: 6750.4
- 散金: 1454
- 红花: 22
- 帖子: 2388
- 在线: 593.5小时
- 虫号: 2762045
- 注册: 2013-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2015-12-21 12:48:33
zyqh
铁杆木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 6750.4
- 散金: 1454
- 红花: 22
- 帖子: 2388
- 在线: 593.5小时
- 虫号: 2762045
- 注册: 2013-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2015-12-21 12:55:44
zyqh
铁杆木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 6750.4
- 散金: 1454
- 红花: 22
- 帖子: 2388
- 在线: 593.5小时
- 虫号: 2762045
- 注册: 2013-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
7楼2015-12-21 15:47:26
zyqh8楼2015-12-21 21:33:44
zyqh
铁杆木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 6750.4
- 散金: 1454
- 红花: 22
- 帖子: 2388
- 在线: 593.5小时
- 虫号: 2762045
- 注册: 2013-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
9楼2015-12-21 22:19:15
10楼2015-12-21 22:34:43