用户工具

站点工具

本页面的其他翻译:
  • zh

adf:sressandforce

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/bohr3
  • 压强(给定固体结构,三个方向压强未必相同),因此也是一个3*3张量,平均值为-1629.20kbar
adf/sressandforce.txt · 最后更改: 2024/02/06 16:00 由 liu.jun

© 2014-2022 费米科技(京ICP备14023855号