24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1042  |  回复: 3

wzjk

金虫 (小有名气)

[求助] QE计算应力应变关系的问题已有1人参与

在用第一性原理计算压力的时候,发现在优化晶胞参数以及原子位置时候,输出结果中有(使用参数为vc-relax)
          total   stress  (Ry/bohr**3)                   (kbar)     P=   -0.18
  -0.00000120  -0.00000000  -0.00000000         -0.18     -0.00     -0.00
   0.00000000  -0.00000120   0.00000000          0.00     -0.18      0.00
   0.00000000   0.00000000  -0.00000120          0.00      0.00     -0.18

。。。。。。。如果我要计算压力与晶体结构变化的关系,是不是应该改变晶胞参数,然后用relax参数计算出新的
total stress,然后与优化的构型算出的stress相减,差值才为作用到晶胞上的压力呢?
=======================
另外,我发现在计算时候可以在&CELL里用 press= 参数来设置外界压力,但是这个压力有方向没有呢?有没有人知道的求解答?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wzjk: 金币+10, ★★★★★最佳答案 2015-07-13 10:07:37
1 0.1kbar 是计算误差内的,可以不用管,当做0就好了,如果想更小可以提高收敛标准
2 这个压力各向同性

» 本帖已获得的红花(最新10朵)

2楼2015-07-07 16:14:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wzjk

金虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by KalaShayminS at 2015-07-07 16:14:48
1 0.1kbar 是计算误差内的,可以不用管,当做0就好了,如果想更小可以提高收敛标准
2 这个压力各向同性

您好 我还有个问题想问下。

在用QE算压力的时候 我得到输出结果
          total   stress  (Ry/bohr**3)                   (kbar)     P=    7.90
   0.00000146   0.00000000  -0.00000000          0.22      0.00     -0.00
   0.00000000   0.00000146  -0.00000000          0.00      0.22     -0.00
  -0.00000000  -0.00000000   0.00015816         -0.00     -0.00     23.27
我的理解为23.27为z轴方向的压力,而0.22和0.22为x和y轴的压力。然后P=7.90又是指的什么呢?
3楼2015-07-13 10:08:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)

引用回帖:
3楼: Originally posted by wzjk at 2015-07-13 10:08:49
您好 我还有个问题想问下。

在用QE算压力的时候 我得到输出结果
          total   stress  (Ry/bohr**3)                   (kbar)     P=    7.90
   0.00000146   0.00000000  -0.00000000          0.22  ...

xyz三个值(对角元)的平均值。
4楼2015-07-13 10:33:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wzjk 的主题更新
信息提示
请填处理意见