Mg3Bi2-xSbx热电材料弹性DFT仿真复现

关键词：Mg3Bi2-xSbx；范德华修正；弹性常数；DFT；热电材料

一、Mg3Bi2-xSbx热电材料弹性DFT研究背景

Mg3Bi2-xSbx作为一种室温高性能热电材料，其力学性能对器件可靠性至关重要。虽然已有第一性原理对其弹性特性进行了计算，但范德华（vdW）修正对这一层状体系的弹性参数影响此前未被系统评估。非局域电荷密度涨落引起的vdW相互作用对层状拓扑材料的结构和力学性质有显著影响，定量评估不同vdW修正方案的效果具有理论和工程双重价值。

中科院力学所与哈工大（深圳）团队在Materials上发表了基于DFT的系统对比研究，建立40原子超胞模型评估了PBE-D3和vdW-DFq两种vdW修正对Mg3Bi2-xSbx（0≤x≤2）弹性性质的影响。研究发现PBE-D3使晶格常数收缩约1%并更接近实验值，除x=0.375特殊点外，两种vdW修正不改变弹性性质随Sb含量变化的趋势，仅影响绝对值（通常≤10%）。

二、计算方法与模型设置

所有DFT计算使用VASP软件包完成，采用PBE泛函作为交换关联基准。对比了三种模型：无vdW修正的PBE原版、PBE-D3色散修正以及vdW-DFq非局域关联泛函。平面波截断能520 eV，Monkhorst-Pack k点网格分辨率2π×0.03 A^(-1)，力收敛判据0.01 eV/A。弹性常数通过应力-应变法计算，对40原子超胞施加±0.5%的六种线性形变导出Cij张量。

Mg3Bi2与Mg3Sb2均属三方晶系（空间群P-3m1），原胞含5个原子。通过2×2×2超胞扩展为40原子模型，其中16个Bi/Sb位点按不同浓度配比（从0到16个Sb原子，共17种构型）。对每种构型进行完全弛豫后计算弹性常数C11/C12/C13/C33/C44，并导出Voigt-Reuss-Hill平均的多晶弹性模量（体模量B、剪切模量G、杨氏模量E、泊松比ν及Pugh比B/G）。

图1 Mg3Bi2-xSbx晶体结构与17种Sb浓度构型（40原子超胞，复现结果）

图2 不同vdW修正下晶格常数a和c随Sb含量x的变化（复现结果）

三、关键仿真结果与发现

PBE-D3方法使晶格常数a与c均收缩约1%，变化幅度显著大于vdW-DFq（~0.3%）。两种vdW修正均不改变晶格常数随Sb含量增加单调下降的趋势，但PBE-D3给出的a轴结果更接近实验数据，增强了其在层状Bi/Sb体系中的可信度。

弹性常数Cij计算发现PBE-D3使面内C11显著增大（最高达15%），反映了vdW层间结合对层内键合的反作用效应。vdW-DFq的修正幅度介于PBE和PBE-D3之间。值得注意的是，在x=0.375（6个Sb原子）处，两种vdW修正对部分弹性模量的趋势产生扰动，作者归因为该浓度下特殊的有序排列构型。

二维杨氏模量极图揭示了弹性各向异性的浓度演化：纯Mg3Bi2（x=0）的各向异性最强，随Sb含量增加各向异性逐渐减弱。Pugh比B/G均大于1.75，表明Mg3Bi2-xSbx全浓度范围为韧性材料，有利于热电模块的加工与服役可靠性。泊松比在0.25~0.35之间，体现了离子-共价混合键特性。弹性各向异性因子AU随Sb浓度整体呈下降趋势，与2D极图分析结论一致。

图3 弹性常数C11/C12/C13/C33/C44在三种vdW方案下的对比（复现结果）

图4 二维杨氏模量极图与弹性各向异性（4种代表性浓度，复现结果）

图5 Voigt-Reuss-Hill体模量B与剪切模量G随Sb含量变化（复现结果）

图6 泊松比nu与Pugh比B/G（不同vdW修正方案对比，复现结果）

图7 弹性各向异性因子AU随Sb浓度演化趋势（复现结果）

四、复现总结与服务

本复现完整覆盖了VASP中Mg3Bi2-xSbx体系的DFT弹性计算流程：超胞建模与17种浓度构型生成、PBE/PBE-D3/vdW-DFq三种交换关联方案对比、应力-应变弹性常数提取、Voigt-Reuss-Hill多晶模量平均及2D弹性各向异性分析。该评估框架可迁移至其他层状热电材料体系（如Bi2Te3、SnSe等）的vdW修正方案筛选。