关键词Mg3Bi2-xSbx范德华修正弹性常数DFT热电材料一、Mg3Bi2-xSbx热电材料弹性DFT研究背景Mg3Bi2-xSbx作为一种室温高性能热电材料其力学性能对器件可靠性至关重要。虽然已有第一性原理对其弹性特性进行了计算但范德华vdW修正对这一层状体系的弹性参数影响此前未被系统评估。非局域电荷密度涨落引起的vdW相互作用对层状拓扑材料的结构和力学性质有显著影响定量评估不同vdW修正方案的效果具有理论和工程双重价值。中科院力学所与哈工大深圳团队在Materials上发表了基于DFT的系统对比研究建立40原子超胞模型评估了PBE-D3和vdW-DFq两种vdW修正对Mg3Bi2-xSbx0≤x≤2弹性性质的影响。研究发现PBE-D3使晶格常数收缩约1%并更接近实验值除x0.375特殊点外两种vdW修正不改变弹性性质随Sb含量变化的趋势仅影响绝对值通常≤10%。二、计算方法与模型设置所有DFT计算使用VASP软件包完成采用PBE泛函作为交换关联基准。对比了三种模型无vdW修正的PBE原版、PBE-D3色散修正以及vdW-DFq非局域关联泛函。平面波截断能520 eVMonkhorst-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之间。值得注意的是在x0.3756个Sb原子处两种vdW修正对部分弹性模量的趋势产生扰动作者归因为该浓度下特殊的有序排列构型。二维杨氏模量极图揭示了弹性各向异性的浓度演化纯Mg3Bi2x0的各向异性最强随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修正方案筛选。