Quantum ESPRESSO计算材料应力与原子受力

AMS软件不仅提供基于STO基组LCAO方法的周期性DFT计算模块BAND，也集成了开源软件Quantum ESPRESSO，为它提供无缝衔接的输入设置与输出结果分析、显示。

设置适当k点数目（沿某方向尺寸越大，则该方向k点数目越小）、泛函与赝势（会根据作业、体系自动调整可选泛函、赝势）： 根据需要勾选Stress或Force：

保存并提交作业。

SCM → Output打开输出文件，窗口底部搜索“Forces acting on atoms”，可以看到类似如下内容：

     Forces acting on atoms (cartesian axes, Ry/au):

     atom    1 type  1   force =     0.03700647   -0.06115831    0.00218578
     atom    2 type  1   force =    -0.03702985   -0.06112002    0.00181016
     atom    3 type  1   force =     0.04087611    0.06066577    0.00442889
     atom    4 type  1   force =    -0.03592144    0.06254055    0.00092244
     atom    5 type  1   force =    -0.07123133    0.00105399    0.00202996
     atom    6 type  1   force =     0.07157081   -0.00022164    0.00297182
     atom    7 type  2   force =    -0.00239153   -0.00009510   -0.00775077
     atom    8 type  2   force =    -0.00073065    0.00056528   -0.00117456
     atom    9 type  2   force =    -0.00214858   -0.00223053   -0.00542371

     Total force =     0.176396     Total SCF correction =     0.000933


     Computing stress (Cartesian axis) and pressure

          total   stress  (Ry/bohr**3)                   (kbar)     P=    -1629.20
  -0.01104628  -0.00000121   0.00000223        -1624.96       -0.18        0.33
  -0.00000121  -0.01105673  -0.00000762           -0.18    -1626.50       -1.12
   0.00000223  -0.00000762  -0.01112226            0.33       -1.12    -1636.14

• 该体系有9个原子，因此有9行数据，依次列出9个原子的受力情况，三列数字代表在x、y、z三个方向的分量，力的单位为Ry/Bohr，其中Ry是能量单位，Bohr是长度的原子单位。
• stress是二阶张量，因此是一个3*3的矩阵，数值单位为Ry/bohr³
• 压强（给定固体结构，三个方向压强未必相同），因此也是一个3*3张量，平均值为-1629.20kbar