如何获得双轴应力拉伸后的相应stress值？
请问对二维结构施加双轴应力后，在OUTCAR中如何获取相应的stress的结果呢？如果是单轴应力取相应方向的值就行了，但是对双轴应力该怎么取值呢？
Direction XX YY ZZ
-------------------------------------------------------------------------
Alpha Z 9.20340 9.20340 9.66357
Ewald 4061.86411 4061.86410 -4952.10567
Hartree 4140.48553 4140.48553 -4247.34363
E (xc) -96.16427 -96.16427 -97.42057
Local -8608.03234 -8608.03224 8773.729887
n-local 219.87098 219.87098 221.535325
augment -42.24293 -42.24293 -44.897975
Kinetic 308.90367 308.90367 335.978935
Fock 0.00000 0.00000 0.000000
--------------------------------------------------------------------------
Total -6.11176 -6.11176 -0.02577
in kB -26.53829 -26.53829 -0.11751
external pressure = -17.72 kB
返回小木虫查看更多
京公网安备 11010802022153号
应该就是这个吧
-26.53829
别忘记乘以真空层厚度
谢谢回复！最终结果是需要乘以真空层厚度的！现在还有一个问题就是，在对单轴a施加应力后，优化另一个轴b，用的是ISIF=3和OPTCELL结合，固定住a和c方向，但是得到的CONTCAR中，b的晶格参数是没变化的，不知道是什么原因。（vasp是根据刘锦程的教程重新编译过的，可以固定真空层优化结构）
，
已解决