这里会显示出您选择的修订版和当前版本之间的差别。
后一修订版 | 前一修订版 | ||
adf:elastictensor [2023/08/13 00:35] – 创建 liu.jun | adf:elastictensor [2023/08/13 09:38] (当前版本) – [DFT 计算弹性张量、体模量、剪切模量、杨氏模量(三维周期边界)] liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
- | ======DFT 计算弹性张量(三维周期边界)====== | + | ======DFT 计算弹性张量、体模量、剪切模量、杨氏模量====== |
+ | |||
+ | =====参数设置===== | ||
+ | {{ : | ||
+ | 基组方面,我们选择了较小的基组DZP,如果存在重元素,则可以设置较大基组,具体参考:[[adf: | ||
+ | 这里我们选择libXC泛函,计算效率会更高一些,并设置礼包XC泛函为普适性较强,计算量较小的PBE-D3: | ||
+ | {{ : | ||
+ | 用户也可以根据精度需求,设置为其他泛函。 | ||
+ | |||
+ | 勾选弹性张量计算: | ||
+ | {{ : | ||
+ | |||
+ | 一般步长设为0.003,结果可能更好一些,步长对最终结果实际上还是存在较大影响的,用户亦可自行测试该值的影响。 | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | 保存并提交作业。这个作业需要非常长的时间完成,因为需要做大量的形变,计算这些形变结构的能量。 | ||
+ | =====结果查看===== | ||
+ | |||
+ | SCM → Output,翻到最后,可以看到: | ||
+ | < | ||
+ | | ||
+ | |||
+ | |||
+ | | ||
+ | 52.938 | ||
+ | | ||
+ | -0.045 | ||
+ | -0.118 | ||
+ | -0.666 | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | | ||
+ | |||
+ | |||
+ | Bulk modulus | ||
+ | Shear modulus [GPa]: | ||
+ | Young modulus [GPa]: | ||
+ | </ | ||
+ | 其中列出了: | ||
+ | * 弹性张量: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 |