晶体压力计算

didi_dudu

在用第一性原理计算压力的时候，发现在优化晶胞参数以及原子位置时候，输出结果中有（使用参数为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= 参数来设置外界压力，但是这个压力有方向没有呢？有没有人知道的求解答？

zhou

press是静水压，计算压力与晶体结构变化的关系，press设不同的值结构优化就可以了。
不过这是静态的关系，不能和实验比较，如果考虑温度，是要算自由能的。