目录

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

其中列出了: