24小时热门版块排行榜    

查看: 2057  |  回复: 2

随便注册

银虫 (初入文坛)

[求助] OUTCAR中stress怎么读取 已有1人参与

导师让重复一篇论文,需要读取OUTCAR中的stress数据。论文中为“stress sigma(希腊符号),单位为N/m” ,现在我在OUTCAR中找到相关的部分,但是不知道读取那个数据。

  FORCE on cell =-STRESS in cart. coord.  units (eV/reduce length):
  Direction    X        Y        Z        XY       YZ       ZX
  --------------------------------------------------------------------------------------
  Alpha Z      8.63      8.63      8.63
  Ewald     2147.14   2147.14  -2819.40      0.00      0.00      0.00
  Hartree   2151.32   2151.32  -2102.86      0.00      0.00      0.00
  E(xc)      -79.02    -79.02    -82.16      0.00      0.00      0.00
  Local    -4568.58  -4568.58   4629.03      0.00      0.00      0.00
  n-local     -3.65     -3.14     -7.45      0.14      0.00      0.00
  augment     -0.92     -0.92     -0.90      0.00      0.00      0.00
  Kinetic    301.03    297.10    374.72      0.08      0.00      0.00
  -------------------------------------------------------------------------------------
  Total      -45.77    -45.77     -0.40      0.00      0.00      0.00
  in kB     -218.32   -218.32     -1.90      0.00      0.00      0.00
  external pressure =     -146.18 kB  Pullay stress =        0.00 kB
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huisheng

铜虫 (小有名气)

热爱生活,热爱学习
2楼2016-05-12 16:52:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awdawdawdd

金虫 (小有名气)

【答案】应助回帖

in KB 那一行就是cell的stress了啊 1 kB= 0.1GPa,如果你要cell本身的stress记得所有的数字要乘以-1 因为这里是force on cell

stress的单位应该是Pa 或者 GPa吧,其中1Pa=N/(m^2) , lz应该漏了一个m,不是N/m

另外 lz结构的stress好大 是不是没有relax?或者是固定了cell模拟高压态么?
3楼2016-05-16 03:07:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 随便注册 的主题更新
信息提示
请填处理意见