VASP计算二维体系应力中的单位问题：KB or N/M（牛/米）? - 第一性原理 - Vasp&MedeA

It is a shame for the soul to be first to give way in this life,when thy body does not give way.

11楼2016-04-12 13:44:16

dongdong33

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 305.4
帖子: 232
在线: 81.8小时
虫号: 472038
注册: 2007-12-04
性别: MM
专业: 半导体材料

in kB的负号表示什么意思？拉应力还是压应力？

12楼2016-04-29 13:51:02

ZHC_Y

金虫 (小有名气)

Dr.

应助: 7 (幼儿园)
金币: 2771.1
红花: 2
帖子: 290
在线: 245.9小时
虫号: 1930069
注册: 2012-08-08
性别: GG
专业: 物理学I

我觉得楼主这里的计算是不是有个地方值得讨论一下. 你的X方向的应力为-150.63161 kB, 而Y方向还有 -19.32335 kB的应力. 从这个计算上来看, 我推测楼主的计算过程是这样的:
1. 把石墨烯的单包由六角晶格变成一个长方形晶格
2. 改变X方向的晶格常数, 然后采用ISIF=2 去优化
3. 最后提取结果
如果是这样计算的话, 那么查看应力的大小的结果就应该像楼主的结果那样, Y方向还有那么大的应力存在.但问题是, Y方向也是应该优化的才对吧? 那样Y方向的应力应该要很小接近0 kB才对.

13楼2017-01-15 21:10:07

往事往亦

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 1868.1
红花: 1
帖子: 680
在线: 213.6小时
虫号: 3159289
注册: 2014-04-24
性别: GG
专业: 凝聚态物性 II ：电子结构

除以3.4这个石墨层间距，我感觉应该是不对的，另外我通过修改了constr_cell_relax.F源代码（只优化x方向-石墨烯），就是对Y方向施加拉伸应变（拉伸x%后固定，并用的ISIF=3），这样做后得到的应力-应变数据为：
in kB -3.23705 -70.14611 -2.60080 0.00000 0.00000 0.00000 1.02
in kB -0.15934 -27.23946 -4.29416 0.00000 0.00000 0.00000 1.04
in kB 0.13776 -57.52228 -4.20366 0.00000 0.00000 0.00000 1.06
in kB -0.00191 -84.64178 -4.00855 0.00000 0.00000 0.00000 1.08
in kB 0.03199 -108.67303 -3.94575 0.00000 0.00000 0.00000 1.10
in kB -0.13334 -128.35889 -3.79733 0.00000 0.00000 0.00000 1.12
in kB 0.00870 -145.29313 -3.87864 0.00000 0.00000 0.00000 1.14
in kB 0.16160 -160.36230 -3.73108 0.00000 0.00000 0.00000 1.16
in kB 0.23975 -172.93999 -3.73327 0.00000 0.00000 0.00000 1.18
in kB -0.13232 -181.16550 -3.70646 0.00000 0.00000 0.00000 1.20
也发现了一些问题：
（1）我想要的效果是，石墨烯的晶格常数b，z值不变，但允许a值变化，并且原子坐标都可以变化，（没达到,,改过代码后，晶格常数b,z不变达到了，a值变化也达到了，但是原子坐标貌似也只优化了x方向的分量，，y,z方向的分量没变化，，）
（2）计算的数据中1.02怎么出现了不遵循从小到大的变化了，，（好奇怪）
（3）我的真空层取得是上下各10A的真空厚度，都乘以20A后，，（我感觉像石墨烯这种高对称，我做的是Y方向拉伸，但仍是应该和你的数据符合哎，但怎么会有这么的差别呢，奇怪，

）

赞一下(1人)

青春不留遗憾。

14楼2017-05-05 17:59:25

曾曾曾伟大

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 17
帖子: 5
在线: 10.3小时
虫号: 7324332
注册: 2017-10-18
专业: 半导体材料

引用回帖:

14楼: Originally posted by 往事往亦 at 2017-05-05 17:59:25
除以3.4这个石墨层间距，我感觉应该是不对的，另外我通过修改了constr_cell_relax.F源代码（只优化x方向-石墨烯），就是对Y方向施加拉伸应变（拉伸x%后固定，并用的ISIF=3），这样做后得到的应力-应变数据为：
in ...

请问你读应力的数据是读第一行还是最后一行的数据呀这个问题困扰我好久了望指教

15楼2018-04-24 09:16:53

嘟嘟格格

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 183.2
帖子: 16
在线: 9.8小时
虫号: 3501629
注册: 2014-10-27
专业: 理论和计算化学

你好，请教一下，你会用VASP计算3D的应力应变曲线吗？

16楼2018-07-20 19:57:09