DFT 计算弹性张量、体模量、剪切模量、杨氏模量

基组方面，我们选择了较小的基组DZP，如果存在重元素，则可以设置较大基组，具体参考：BAND：如何设置基组 这里我们选择libXC泛函，计算效率会更高一些，并设置礼包XC泛函为普适性较强，计算量较小的PBE-D3： 用户也可以根据精度需求，设置为其他泛函。

勾选弹性张量计算：

一般步长设为0.003，结果可能更好一些，步长对最终结果实际上还是存在较大影响的，用户亦可自行测试该值的影响。

保存并提交作业。这个作业需要非常长的时间完成，因为需要做大量的形变，计算这些形变结构的能量。

SCM → Output，翻到最后，可以看到：

 Elastic Tensor [GPa] in Voigt notation

  
     176.208      52.938       0.572      -0.045      -0.118      -0.666
      52.938     187.576      -0.746       0.251       0.065       0.680
       0.572      -0.746       3.548       3.295       0.500       0.182
      -0.045       0.251       3.295     -32.958     -36.534      -0.206
      -0.118       0.065       0.500     -36.534      17.300      -0.301
      -0.666       0.680       0.182      -0.206      -0.301      66.113




 Elastic moduli

  
Bulk modulus  [GPa]:    52.540
Shear modulus [GPa]:    31.062
Young modulus [GPa]:    77.846

其中列出了：

• 弹性张量：C$_{11}$ = 176.208 GPa，C$_{22}$ = 187.576 GPa，这里如果材料XY方向各向同性，则一般需要对这 2 个数值进行平均。C$_{12}$ = C$_{21}$ = 52.938 GPa。
• 体模量：52.540 GPa
• 剪切模量：31.062 GPa
• 杨氏模量：77.846 GPa