| 查看: 15444 | 回复: 15 | |||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||
bjwang铁杆木虫 (正式写手)
|
[求助]
VASP计算二维体系应力中的单位问题：KB or N/M（牛/米）? 已有1人参与
|
||||||
|
利用改变晶格常数的方法计算应力应变的关系（stress-strain curve）：epsilon = (a-a_{0})/a_{0}； 具体计算中：设置ISIF = 2，真空层取1nm。 验证2D的石墨烯在5%的形变下的应力（文献中的计算值近似为150N/M或40GPa） 但我在OUTCAR中找到 “in kB -150.63161”，单位是KB。 1，换算成GPa 按照1KB=0.1GPa，我的结果150.63161KB=15.6GPa，与40GPa不符； 2，换算成N/M 150.63161KB=15.6GPa=15.6*10^9牛/米平方， 那么按照有关文献的说法，算二维体系的应力，考虑单位长度上的力，而不是单位面积上的，即单位要表示为牛/米， 牛/米平方要乘上真空层的厚度（这个不确定）？ 但结果也不对：15.6*10^9牛/米平方*1nm=15.6N/M，和文献中的150N/M刚好错了一个数量级。 但有个帖子（http://muchong.com/html/201406/7592923.html）说，在vasp计算石墨烯的杨氏模量中，单位就是N/M（为什么？？？） 那上面的应力的单位也应是N/M? 如果是这样，结果就和文献符合的很好，但是我上面的换算过程哪里有问题呢？ 还是计算有问题？ 谢谢各位！ 附计算结果 FORCE on cell =-STRESS in cart. coord. units (eV): Direction XX YY ZZ XY YZ ZX -------------------------------------------------------------------------------------- Alpha Z 11.84580 11.84580 11.84580 Ewald 1223.83420 1231.72516 -1753.82933 0.00000 0.00000 0.00000 Hartree 1199.50298 1201.90716 -1177.06490 0.00000 0.00000 0.00000 E(xc) -57.55665 -57.59624 -59.87916 0.00000 0.00000 0.00000 Local -2603.83155 -2612.99666 2720.52797 0.00000 0.00000 0.00000 n-local -12.23578 -13.54492 -12.15478 0.08818 0.00000 0.00000 augment -1.40512 -1.40962 -1.31894 0.00000 0.00000 0.00000 Kinetic 229.49916 238.74199 271.32154 -4.37353 0.00000 0.00000 Fock 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ------------------------------------------------------------------------------------- Total -10.34696 -1.32733 -0.55181 0.00000 0.00000 0.00000 in kB -150.63161 -19.32335 -8.03326 0.00000 0.00000 0.00000 external pressure = -59.33 kB Pullay stress = 0.00 kB 附文献stress-strain curve图：  N每M.png  GPa.png |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理计算-tg | 表界面计算 | VASP计算资料收集byZF | vasp计算方法收集 |
| 第一性 |
» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp加电场弛豫原子位置很怪异
已经有5人回复
- vasp计算能带用origin处理遇到的问题请教
已经有0人回复
- 计算输运的免费软件
已经有6人回复
- 计算含有过渡金属的氧化物二维体系如何设置自旋情况
已经有6人回复
- 求助关于用vasp计算二维石墨烯纳米带所得到的能带结果与报道对不上
已经有7人回复
bjwang
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 5860.3
- 散金: 94
- 红花: 6
- 帖子: 792
- 在线: 556.3小时
- 虫号: 250077
- 注册: 2006-05-08
- 专业: 凝聚态物性 II ：电子结构
送红花一朵 5楼2016-03-29 12:36:15
2楼2016-03-27 20:22:02
bjwang
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 5860.3
- 散金: 94
- 红花: 6
- 帖子: 792
- 在线: 556.3小时
- 虫号: 250077
- 注册: 2006-05-08
- 专业: 凝聚态物性 II ：电子结构
3楼2016-03-27 20:59:15
4楼2016-03-28 14:49:47
