如何獲得雙軸應(yīng)力拉伸后的相應(yīng)stress值?
請問對二維結(jié)構(gòu)施加雙軸應(yīng)力后,在OUTCAR中如何獲取相應(yīng)的stress的結(jié)果呢?如果是單軸應(yīng)力取相應(yīng)方向的值就行了,但是對雙軸應(yīng)力該怎么取值呢?
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
返回小木蟲查看更多
京公網(wǎng)安備 11010802022153號
應(yīng)該就是這個吧
-26.53829
別忘記乘以真空層厚度
謝謝回復(fù)!最終結(jié)果是需要乘以真空層厚度的!現(xiàn)在還有一個問題就是,在對單軸a施加應(yīng)力后,優(yōu)化另一個軸b,用的是ISIF=3和OPTCELL結(jié)合,固定住a和c方向,但是得到的CONTCAR中,b的晶格參數(shù)是沒變化的,不知道是什么原因。(vasp是根據(jù)劉錦程的教程重新編譯過的,可以固定真空層優(yōu)化結(jié)構(gòu))
,
已解決