基组方面,我们选择了较小的基组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
其中列出了: