学术研究报告：二硫化钼单层中2H与1T相变机制的理论研究
本报告旨在向国内科研同行详细介绍一篇关于二维过渡金属二硫化物相变机理的重要理论研究论文。该工作深入探讨了二硫化钼单层中半导体性2H相与金属性1T/1T’相之间相变的原子级机制，为理解和调控此类二维材料的物性提供了关键的理论基础。
一、 研究基本信息
本研究由香港理工大学纺织与服装学系的赵文和丁峰*教授共同完成。研究成果以题为 “Energetics and kinetics of phase transition between a 2H and a 1T MoS2 monolayer—a theoretical study” 的论文形式，发表于2017年的 《Nanoscale》 期刊上（2017, 9, 2301）。该期刊由英国皇家化学会出版，论文接收于2016年11月4日，并于2017年1月6日在线发表。
二、 研究背景与目的
科学领域： 本研究属于凝聚态物理、材料科学和计算材料学交叉领域，具体聚焦于二维材料的相变物理与晶体生长动力学。
研究背景与动机： 过渡金属二硫化物作为一类重要的二维材料，因其丰富的相结构和可调的电子特性（如从半导体到金属的转变）而备受关注。其中，二硫化钼（MoS2）单层存在多种晶体结构，最常见的两种是：热力学稳定的半导体性 2H相（金属原子呈三棱柱配位）和亚稳态的金属性 1T相（金属原子呈八面体配位）。1T相还可进一步重构为更稳定的 1T’相（金属原子形成之字形链）。实验上已能通过锂/钠插层、电荷注入或化学功能化等手段诱导2H到1T/1T’的相变。这种可控的相变在电子学、自旋电子学、催化以及能源存储等领域具有巨大的应用潜力。然而，对于相变过程中新相如何在旧相中成核、两相界面的原子结构、以及界面如何迁移扩展等关键物理问题，尚缺乏系统深入的理论理解。传统的第一性原理计算通常只关注最小的周期性原胞，难以捕捉相变这一涉及成核与界面演化的复杂过程。
研究目的： 本研究旨在通过结合经典晶体生长理论与第一性原理计算，系统地阐明MoS2单层中2H到1T/1T’相变的全过程，包括：1）1T相核在2H晶格中的成核行为；2）2H-1T之间可能界面结构的形成能与稳定性；3）相变过程中界面迁移的动力学速率。研究的目标是建立一个完整的理论框架，以解释已有的实验观察（如三角形1T畴的形状、缓慢的界面扩展速度），并预测在不同化学环境下的相变行为。
三、 详细研究方法与流程
本研究是一项完全基于计算模拟的理论工作，其核心流程是理论建模与第一性原理计算分析，不涉及实验样品。研究流程可分为以下几个步骤：
1. 体系建模与计算方法设定： * 研究对象与模型： 研究核心对象是MoS2单层。为研究相变，作者构建了两种主要模型：(1) 畴模型：在大的2H相超胞中嵌入不同尺寸和形状（如三角形）的1T/1T’相区域，用于研究成核和界面形成能；(2) 带状模型：构建包含长直界面的2H-1T’超胞，用于模拟界面（如锯齿型边界）的一维传播动力学。 * 计算方法： 全部计算基于密度泛函理论（Density Functional Theory, DFT）。具体使用维也纳从头算模拟软件包（VASP）进行能量和结构优化计算。交换关联泛函采用广义梯度近似的Perdew-Burke-Ernzerhof（PBE）形式。采用投影级加平面波（PAW）方法处理离子实。由于超胞尺寸较大（>18 Å），布里渊区采样仅使用Γ点。平面波截断能设置为280 eV，几何优化收敛标准为原子间力小于0.01 eV/Å。垂直于单层方向设置了约15 Å的真空层以避免周期性镜像相互作用。 * 关键参数验证： 作者声明对超胞尺寸、k点采样密度和截断能等参数进行了仔细测试，确保了计算精度。
2. 可能界面结构的提出与稳定性分析： * 界面结构设计： 基于2H-MoS2的D3h对称性和实验观察，提出了六种高对称性的2H-1T界面构型。分类依据两个维度：(a) 界面方向：分为armchair（AC，扶手椅）方向和两种zigzag（ZZ，锯齿）方向（分别为1T相以S原子终止的ZZ-S界面和以Mo原子终止的ZZ-Mo界面）；(b) 边界处的硫原子链化学计量：分为富硫型（标记为|+）和缺硫型（标记为|-）。由此组合出六种结构：ZZ-Mo|-, ZZ-Mo|+, ZZ-S|-, ZZ-S|+, AC|-, AC|+。 * 界面形成能计算： 使用畴模型计算不同尺寸三角形1T’畴在2H晶格中的形成能。形成能公式考虑了体系总能量、硫原子化学势（μS）的变化以及相变驱动力（δμ，代表通过掺杂或电荷转移使1T’相相对2H相稳定的能量差）。通过拟合形成能与畴边界长度的关系（公式：G_f = 3ε_v(μ_S) + 3Lγ_b(μ_S) + (√3/4)L^2δ），可以分离出顶点形成能ε_v、单位长度边界能γ_b和体相形成能差δ。边界能γ_b是硫化学势μS的函数，这反映了界面稳定性对生长环境的依赖性。
3. 相稳定性与1T’相选择： * 在开始详细界面研究前，作者首先比较了2H、1T、1T’和1T”四种MoS2相的相对稳定性。计算确认，对于有限尺寸的畴，初始的1T结构在弛豫中会重构为1T”等结构，而1T’结构则保持稳定且能量更低。因此，后续所有关于相变的计算均采用更稳定的1T’相作为金属相的代表。
4. 成核过程的理论分析： * 基于经典成核理论，将1T’相在2H晶格中的形成视为一个热力学过程。通过分析三角形畴形成能G_f与畴大小（n_1t，正比于面积）的关系，可以找到成核能垒G和临界核尺寸n_1t。分析中考虑了相变驱动力δμ的影响，δμ代表了外部刺激（如电荷注入）对1T’相的稳定化作用。
5. 界面传播动力学的模拟： * 模拟方法： 使用带状模型研究界面（尤其是ZZ和AC边界）的逐步扩展机制。该过程被建模为一个可重复的循环，包括两个步骤：(i) 在直边界上形成一个一维（1D）核；(ii) 该核两侧的扭折（kink）向外传播，导致界面整体前进。 * 能量路径计算： 相变过程由界面附近硫原子的迁移驱动。作者系统地移动边界硫原子，计算了每个中间构型的形成能变化（δg_i,j），并使用爬坡式弹性带（Climbing Image Nudged Elastic Band, CI-NEB）方法计算了硫原子迁移的能量势垒。 * 动力学速率推断： 通过比较不同界面构型下形成一维核所需的能量（δg*）和原子迁移势垒，可以定性判断不同方向界面的传播速度快慢。高δg*意味着成核困难，界面迁移速度慢。
6. 形状演化的理论预测： * 结合热力学（边界形成能γ_b）和动力学（界面传播速率v）的分析，分别应用Wulff构造（平衡形状）和动力学Wulff构造（生长形状）理论，预测了1T’畴在不同化学环境（μS）和不同过程（生长 vs. 收缩）中的预期形状。
四、 主要研究结果及其逻辑关联
1. 界面稳定性结果： * 计算得到的六种界面结构的边界形成能γ_b强烈依赖于硫化学势μS。在大多数实验可行的μS范围内（从贫硫到富硫），ZZ-Mo|-、ZZ-S|+ 和 ZZ-Mo|+ 这三种锯齿型边界在能量上最为有利，而AC边界和其他ZZ边界稳定性较差。 * 体相形成能差δ 的计算值（约0.39-0.57 eV/MoS2）与使用最小原胞计算的结果（0.55 eV/单位）吻合良好，验证了模型的可靠性。 * 形状预测（热力学平衡）： 根据Wulff构造，1T’畴的平衡形状取决于μS。在贫硫（Mo-rich）条件下，畴呈以ZZ-Mo|-边界为主的三角形；随着硫含量增加，形状变为由ZZ-Mo|-和ZZ-S|+边界组成的六边形；在极富硫条件下，最终演变为以ZZ-Mo|+边界为主的三角形。这从热力学上解释了实验中常观测到的三角形1T畴结构。
2. 成核过程结果： * 计算表明，即使在1T’相因外部因素（δμ）而变得比2H相更稳定之后，在完美的2H晶格中形成1T’核仍然面临极高的成核能垒（在δμ < 0.84 eV时，能垒>10 eV；即使δμ增大，能垒也仅缓慢降低至约4.0 eV）。这是因为三角形核的形成包含一个恒定的顶点能量惩罚项（3ε_v）。 * 科学意义： 这一结果解释了为什么在实验中，1T相的成核通常从2H相的缺陷位点或边缘开始，因为这些位置可以显著降低成核能垒，与文献报道一致。
3. 界面传播动力学结果： * AC边界传播： AC边界上硫原子的迁移势垒很低（<0.5 eV），且形成一维核的能量（δg*_ac ≈ 0.07 eV）可忽略不计。这表明**AC边界的迁移速度非常快**。 * **ZZ边界传播：** 相比之下，ZZ边界（ZZ-Mo|-, ZZ-S|+, ZZ-Mo|+）上形成一维核需要克服很高的能量势垒（δg* 分别为 4.05 eV, 3.79 eV, 2.68 eV）。这是因为ZZ边界迁移涉及新的Mo之字形链形成和显著的晶格重构。 * **动力学速率排序：** 由此推断界面迁移速度的顺序为：在贫硫环境下（μS < 0.67 eV，稳定界面为ZZ-Mo|-），**v_ZZ-Mo|- < v_ZZ-S|+ < v_AC < v_others**；在富硫环境下（μS > 0.67 eV，稳定界面为ZZ-Mo|+），v_ZZ-S|+ < v_ZZ-Mo|+ < v_AC < v_others。锯齿型边界的传播速度远慢于扶手椅型边界。
4. 畴生长与收缩的形状演化预测： * 生长过程（动力学主导）： 在1T’畴扩张时，传播最慢的边界将主导畴的最终形状。根据动力学Wulff构造，在贫硫环境下，生长畴将呈现以慢速ZZ-Mo|-边界为主的三角形；在富硫环境下，则呈现以慢速ZZ-S|+边界为主的三角形。这与热力学平衡形状预测一致，并得到了实验观察的支持。 * 收缩过程（反向相变）： 在1T’畴收缩（即逆转为2H相）时，情况相反。传播最快的边界（AC及其他高指数边界）会迅速消失，最终导致畴的周长趋于圆形，与具有清晰低指数边的生长形状形成鲜明对比。
结果间的逻辑链条： 首先，界面稳定性计算确定了哪些边界结构在能量上占优，从而预测了热力学平衡形状。接着，成核分析揭示了相变启动的困难性，指向了缺陷辅助成核的必要性。然后，动力学模拟揭示了不同界面迁移速率的巨大差异，这直接决定了非平衡生长过程中的畴形貌。最后，结合热力学和动力学，完整地解释了畴在生长和收缩两个方向上的不同形状演化路径。整个研究从原子尺度构建了从成核到生长/收缩的完整相变物理图像。
五、 研究结论与价值
结论： 本研究通过系统的第一性原理计算，揭示了MoS2单层中2H到1T/1T’相变的原子机制。主要结论包括： 1. 界面稳定性： ZZ-Mo|-、ZZ-S|+和ZZ-Mo|+是能量最稳定的两相界面结构。1T’畴的平衡形状为三角形或六边形，具体取决于硫化学环境。 2. 成核特性： 在完美2H晶格中形成1T’核的能垒极高，这解释了实验观察中相变常始于缺陷或边缘的现象。 3. 动力学行为： AC边界的迁移速度极快，而ZZ边界的迁移则受到高能量势垒的限制，速度缓慢。这直接导致了实验观测到的缓慢界面扩展速率（~1 nm/min）。 4. 形状演化： 生长中的1T’畴由迁移最慢的ZZ边界主导而呈三角形；收缩中的畴则由迁移最快的边界主导而趋于圆形。
科学价值与应用意义： * 科学价值： 本研究提供了一个研究二维材料相变的普适理论框架，将经典晶体生长理论（成核、Wulff构造）与第一性原理计算深度融合，清晰地解析了相变过程的热力学（能量）与动力学（速率）控制因素。它不仅解释了MoS2中的具体现象，其方法也适用于其他MX2材料（如MoSe2, WS2等）乃至更广泛的二维材料体系。 * 应用价值： 理解相变机制对于精准调控二维材料的电子结构至关重要。研究指出的界面类型、成核位点控制以及生长/收缩动力学的各向异性，为设计基于相变工程的新型器件（如相变存储器、催化活性位点调控、异相结电子器件）提供了关键的理论指导。例如，通过控制硫化学势和环境，可以预测和制备具有特定形状和界面的金属性畴，从而优化其电学或催化性能。
六、 研究亮点
七、 其他有价值内容
论文还提供了详细的电子补充信息，包含了公式推导的细节、额外的图表以及计算参数，确保了研究的可重复性和透明度。此外，研究中对不同MX2材料（MoS2, MoSe2, MoTe2, WS2, WSe2, WTe2）四种相（2H, 1T, 1T‘, 1T’‘）的形成能进行了对比计算（见表1），指出WTe2的1T’相能量甚至低于其2H相，这拓宽了研究的视野，暗示了材料本征稳定性对相变行为的根本影响。
这项研究是一项计算材料学领域的范例工作，它通过精湛的理论建模和计算，深刻揭示了二维材料中复杂相变过程的微观物理机制，对相关领域的研究者具有重要的参考价值。