二维材料形成能计算

一、项目背景与科学问题

以单层过渡金属硫属化物（TMD）为代表的二维材料因原子级厚度、可调带隙与优异的光电性能而备受关注。实验合成二维材料时，其热力学稳定性与缺陷行为直接决定样品质量与器件性能。形成能是衡量材料稳定性的核心热力学量：形成能为负表示化合物相对于单质稳定，绝对值越大越容易合成；空位形成能则反映晶格中原子空位的产生难易，直接影响载流子散射与催化活性位密度。本项目的科学问题包括：MoS2、MoSe2、WS2、WSe2等单层TMD的形成能各为多少？硫空位与金属空位的形成能随化学势如何变化？利用密度泛函理论系统计算可给出定量答案。

二、计算方法与理论背景

二维材料形成能定义为E_f=E_tot-E_bulk_ref的标准化形式，反映自标准态单质生成化合物所释放的能量。以MoS2为例，其每化学式单元形成能可写作E_f=E(MoS2)-E(Mo_bulk)-2×E(S_bulk)并除以单元数，其中单质参考态需注意取最稳定的体相结构。空位形成能的计算采用超胞法：构造含一个空位的大超胞，E_v=E(defect)-E(perfect)+E(host_atom)，其中移走的原子放回其化学势参考库。由于二维材料存在两个元素端（贫硫与富硫），空位形成能强烈依赖生长环境的化学势，需在允许的热力学窗口内扫描化学势并绘制形成能随化学势的变化曲线。

三、计算设置与模型构建

以单层TMD为对象构建原胞，随后扩为4×4×1超胞并沿垂直方向加入15埃真空层，以消除层间周期性镜像相互作用。VASP计算采用PAW方法与PBE泛函，平面波截断能450 eV。k点采样使用Gamma中心网格，超胞取标准精度，结构优化至原子受力小于0.02 eV每埃。单质参考态分别计算体相Mo、W的BCC结构与斜方S、六方Se的总能。空位模型在4×4×1超胞中移除中心S或Mo原子后重新弛豫，计算形成能，并让硫化学势在贫硫端与富硫端之间变化以覆盖不同生长条件。

四、结果分析与案例图解读

图3的左侧面板（a）给出四种单层TMD的形成能柱状图。MoS2的形成能约为-1.86 eV每化学式单元，MoSe2约-1.62 eV，WS2约-1.81 eV，WSe2约-1.58 eV。全部为负值说明这四种二维材料在热力学上均可稳定存在；硫化物（MoS2、WS2）的形成能绝对值大于相应的硒化物，表明S体系成键更强、更易合成，与实验中MoS2、WS2易于剥离制备的现象一致。同族对比显示W基与Mo基形成能接近，反映二者离子半径与电负性相似，化学成键性质相当。

图3的右侧面板（b）给出硫空位与金属空位形成能随硫化学势的变化。横轴硫化学势从-1.4 eV（贫硫端，硫活性低）变化到0 eV（富硫端，硫活性高），红色曲线为硫空位形成能，蓝色虚线为金属空位形成能。可见硫空位形成能随硫化学势升高而增大：在贫硫条件下硫原子容易脱出晶格，形成能低至约1.9 eV；富硫条件下硫化学势高，空位形成能升高至约2.7 eV。金属空位形成能则随硫化学势升高而下降，因为富硫环境有利于金属原子脱离并进入富硫气氛。两条曲线交叉说明通过调节生长气氛可以控制主导缺陷类型，贫硫条件倾向产生硫空位，富硫条件倾向抑制硫空位、促进金属空位，这对理解二维材料的n型掺杂来源与催化活性调控具有重要意义。

两面板结合，左面板给出了材料本征的合成稳定性，右面板揭示了缺陷形成的热力学窗口，共同完成了二维材料形成能的系统计算，为实验生长参数优化提供了理论依据。

五、结论与展望

本项目用VASP计算了四种单层TMD的形成能与空位形成能，确认硫化物比硒化物更稳定，空位形成能显著依赖化学势并可据此调控缺陷类型。案例图以形成能柱状图与空位形成能-化学势曲线双面板呈现。后续可拓展：(1)计算晶界与边缘形成能，理解二维材料形貌演化；(2)引入杂化泛函考察缺陷态的电子结构；(3)将形成能计算与机器学习势结合，实现大体系缺陷的高通量筛选。