二维过渡金属硫族化合物莫尔超晶格中的旋转与膨胀重构研究
一、 研究团队与发表信息
本研究的主要作者包括Madeline Van Winkle（第一作者）、Isaac M. Craig（共同第一作者）、Stephen Carr、Medha Dandu、Karen C. Bustillo、Jim Ciston、Colin Ophus、Takashi Taniguchi、Kenji Watanabe、Archana Raja、Sinéad M. Griffin 以及通讯作者 D. Kwabena Bediako。研究团队主要来自美国加州大学伯克利分校、劳伦斯伯克利国家实验室分子铸造厂、布朗大学以及日本国立材料科学研究所。
该研究成果以论文形式发表于 Nature Communications 期刊，于2023年5月在线发表（论文编号：10.1038/s41467-023-38504-7）。
二、 学术背景与研究目标
本研究属于凝聚态物理与材料科学交叉领域，聚焦于二维范德华材料的莫尔超晶格。由过渡金属硫族化合物（TMDs，如MoS₂、WSe₂）构成的双层莫尔结构，因其可通过层间转角（θm）和/或晶格失配（δ）灵活调控电子能带结构，从而产生新颖的光电特性（如莫尔激子、关联电子相）而备受关注。在这些系统中，晶格重构（lattice reconstruction）及其伴随的应变累积对于定义其电子结构起着关键作用。然而，此前对TMD莫尔重构的理解大多基于扫描探针显微术或电子显微成像的定性观察，以及理论模拟的预测。重构发生的具体物理机制——即原子层如何通过局部机械形变（如旋转、拉伸/压缩）来弛豫到更低能量状态——缺乏直接的、定量的实验测量。
因此，本研究旨在直接、定量地揭示TMD莫尔双层（包括同质双层和异质双层）中驱动晶格重构的层内机械形变机制。具体目标包括：1）区分并量化局部旋转（rotation）和局部膨胀/收缩（dilation）在重构过程中的相对重要性；2）研究转角依赖性和晶格失配的影响；3）探究封装层（如六方氮化硼，hBN）对重构模式（面内形变 vs. 面外起伏）的调控作用；4）探索外部单轴异质应变（heterostrain）对重构应变场的操控能力。
三、 详细研究流程与方法
本研究采用了一套结合先进样品制备、高精度四维扫描透射电子显微镜（4D-STEM）成像以及定制化数据分析算法的完整工作流程。
1. 样品制备： 研究制备了两种主要类型的样品：小转角扭曲双层MoS₂（同质双层） 和 WSe₂/MoS₂异质双层。样品均采用机械剥离和干法转移技术制备。对于扭曲同质双层，采用“撕裂-堆叠”（tear-and-stack）方法精确引入目标莫尔转角（θm）。异质双层则通过堆叠两个独立的TMD单层并使其边缘对齐获得。所有莫尔结构均被封装在薄层（5-10 nm）六方氮化硼（hBN）中，以保护样品并研究封装效应。此外，还制备了部分封装（仅单面有hBN）和完全悬浮（无hBN）的WSe₂/MoS₂异质双层样品作为对照。利用二次谐波产生光谱（SHG）确认了堆叠取向（平行，P，~0°；反平行，AP，~60°）。
2. 数据采集：布拉格干涉测量法（Bragg Interferometry） 核心实验技术是基于4D-STEM的布拉格干涉测量法。这是一种新颖的成像方法，其原理如下： * 设备：使用像差校正的Thermo Fisher Scientific Titan 80-300显微镜（TEAM I），在80 kV电压下操作，配备Gatan K3直接电子探测相机。 * 过程：电子束在样品上逐点扫描。在每个扫描位置，探测器不仅记录透射束强度（如传统STEM），而是记录完整的二维衍射图案（即一个4D数据集：二维实空间扫描位置 × 二维倒易空间衍射强度）。 * 干涉原理：当电子束穿过双层TMD时，会被上下两层晶体衍射。来自两层的衍射电子波在重叠的布拉格衍射盘区域发生干涉。干涉图案的强度调制取决于局部层间原子位移矢量u，即一个原子层相对于另一个的局部偏移。 * 数据分析基础：通过拟合每个扫描点处十二个布拉格盘重叠区域的强度（Ij），利用推导出的三角函数关系式（论文中的公式(1)），可以反推出每个像素点对应的二维位移矢量u(x, y)。该公式是对此前用于中心对称材料（如石墨烯）方法的扩展，以适应TMD非中心对称及异质双层的情况。 * 优势：此方法不依赖于精确追踪衍射盘位置（在莫尔系统中因信号重叠和信噪比低而困难），且能“穿透”hBN封装层进行测量，克服了扫描探针技术通常需要暴露样品表面的限制。
3. 数据处理与应变映射： 获得位移场u(x, y)后，通过一系列定制化算法进行深入分析： * 位移场解包裹（Unwrapping）：由于位移矢量u的周期性（加上整数倍晶格矢量后物理等效）和符号模糊性，原始位移场是不连续、非平滑的。研究团队开发了混合整数规划（mixed-integer program）算法来精修区域边界，将位移场“解包裹”成一个平滑、连续的矢量场，这是后续应变计算的关键前提。 * 应变张量计算：对平滑后的位移场进行数值微分（采用中心三点有限差分法），得到位移梯度张量。假设形变为无穷小应变，将该张量分解为一个对称的应变张量（ε）和一个反对称的旋转矩阵（ω）。由此计算出关键的力学形变量： * 重构旋转（θr）：局部固定体旋转（θt）与预设的全局莫尔转角（θm）的一半之差（θr = θt - θm/2），反映了为降低层间堆叠能而发生的局部晶格旋转。 * 膨胀/收缩（Dilation）：应变张量的迹（Tr(ε)），代表局部面内体积的相对变化，即局部晶格常数的相对改变。 * 最大工程剪切应变（γmax）：应变张量最大与最小特征值之差，表征局部剪切形变的强度。 * 为了分离出重构本身的贡献，从计算出的总形变中减去了一个具有相同转角、晶格失配或异质应变的“刚性”莫尔结构（即无重构）所对应的理论形变背景。 * 堆叠区域分类与统计分析：根据位移矢量u将每个实空间像素归类到特定的高对称堆叠序列（如MMXX, XM/MX, SP用于P堆叠；XMMX, MM, XX, SP用于AP堆叠）。然后对属于同一堆叠类型的像素区域进行平均，得到该堆叠类型对应的平均应变值（θr, 膨胀，γmax）和相对面积占比。 * 计算模拟验证：研究还进行了密度泛函理论（DFT）计算和基于弹性理论的模拟，以验证实验观察到的应变场，并理论估算面外起伏对测量膨胀的影响。
四、 主要研究结果
1. 同质双层中的旋转主导重构： 对于扭曲双层MoS₂（同质双层），应变映射结果明确显示，局部旋转（θr）是驱动晶格重构的主导机制。 * 旋转场分布：在P堆叠中，高堆叠能区域（MMXX）显示θr > 0°，导致该区域收缩；而低堆叠能区域（XM/MX）显示θr < 0°，导致其扩展为三角形畴。在AP堆叠中，高能XX区域θr > 0°（收缩），低能XMMX区域θr < 0°（扩展为六角形畴）。这种旋转场的分布与基于层间堆叠能差异的弛豫预期完全一致。 * 膨胀作用微弱：与旋转场形成鲜明对比的是，在同质双层中，未发现膨胀与堆叠序之间的系统性关联。这表明面内的体积应变（晶格常数变化）对同质双层重构贡献甚微。 * 剪切应变积累：旋转主导的重构自然产生了层内剪切应变场（γmax）。实验测量发现，剪切应变高度集中在具有相同θr符号区域之间的边界（即鞍点/畴壁，soliton/domain walls）处。这与旋转重构模型的模拟预测高度吻合，证实了剪切应变是旋转重构的内在产物，并可能对电子能带的平坦化和载流子局域化起关键作用。 * 转角依赖性：研究发现AP堆叠存在两个重构区间。当θm大于一个临界角（θc ≈ 1.25–1.5°）时，MM区域具有轻微负的θr（略微扩张）；当θm小于θc时，MM区域的θr变为正值（收缩），以容纳XMMX区域的大幅扩张。相比之下，P堆叠在所测转角范围内只观察到一个重构区间。总体而言，P堆叠比AP堆叠重构更强烈，在θm ≈ 2°时，其高能MMXX区域的面积已缩减至刚性模型的36%，而AP堆叠的XX区域仍保持刚性模型的79%。
2. 异质双层中的膨胀主导重构： 对于WSe₂/MoS₂异质双层（具有~3.96%的固有晶格失配），重构机制发生了根本性变化。 * 膨胀场分布：实验观察到了显著的、依赖于堆叠序的膨胀/收缩模式。在高堆叠能区域（AP中的XX，P中的MMXX），测量到正的膨胀（晶格被拉伸）；在低堆叠能区域（AP中的XMMX，P中的XM/MX），测量到负的膨胀（晶格被压缩）。这意味着，系统通过局部调整两层材料的有效晶格失配（即局部晶格常数差）来实现能量弛豫：低能区域通过收缩晶格失配来扩张其面积，高能区域则相反。 * 旋转作用的演变：在近零转角（θm = 0.13°）的AP异质双层中，局部旋转很弱且与堆叠序无关。然而，当引入较大的转角（θm = 1.07°）时，旋转场显著增强，其分布变得与AP同质双层相似。这表明，异质双层可以同时容纳膨胀和旋转两种重构机制，其中旋转的贡献随转角的增大而增强。有趣的是，在非零转角的P异质双层中未观察到明显的堆叠依赖旋转，这可能由于转角未达阈值，或面外起伏削弱了旋转重构（见下文）。
3. 封装层对重构路径的调控： 通过比较完全封装、部分封装和完全悬浮的WSe₂/MoS₂异质双层，研究揭示了hBN封装层的关键作用。 * 增强面内重构：实验发现，完全封装的区域表现出最强、最局域化的膨胀和旋转场；而完全悬浮的区域，这些形变最弱、最弥散。这表明，hBN封装层抑制了面外起伏（corrugation），从而迫使系统更多地通过面内的旋转和膨胀/收缩来进行弛豫。 * 机制验证：研究者构建了一个理论模型来解释面外起伏如何影响测量的面内位移投影。起伏会减小原子投影的视在间距，从而在测量上产生类似负膨胀的效应。通过对比实验测量的悬浮与封装样品间的膨胀差异（δdil），发现其变化趋势与DFT计算的面外起伏模型预测相符。这直接证明，悬浮样品经历了面内和面外混合的弛豫，而封装样品则几乎完全通过面内形变进行重构。此外，封装显著提高了重构超晶格的均匀性。
4. 异质应变对重构应变场的重分布： 研究还在扭曲同质双层中引入了单轴异质应变（ε，一层相对于另一层被拉伸），以模拟可控的晶格失配。 * 旋转机制依然主导：即使异质应变高达~1.4%（对应的有效晶格失配δ约0.5%），也未能诱发显著的膨胀重构。局部旋转仍然是主导的弛豫机制。这说明，只有当晶格失配足够大（如WSe₂/MoS₂的~4%）时，膨胀机制才会变得突出。 * 剪切应变的重构：异质应变的主要效应是重组由旋转重构产生的本征层内剪切应变场。虽然异质应变是均匀施加的，但其影响是各向异性的：应变能集中在畴壁（SP）区域，并放大和扭曲了已有的剪切应变分布。通过改变异质应变的方向（拉伸或压缩莫尔原胞），可以产生条纹状、方形或三角形等不同的剪切应变图案。这为实现通过外部机械力操控莫尔势能景观和载流子局域化提供了新途径。
五、 研究结论与意义
本研究通过发展并应用基于4D-STEM的布拉格干涉测量法，首次在原子尺度上直接、定量地绘制并区分了TMD莫尔双层中两种基本的晶格重构机制——局部旋转和局部膨胀/收缩。研究明确得出结论： 1. 在扭曲TMD同质双层中，重构由局部旋转主导，并伴随产生局域在畴壁的剪切应变场。 2. 在具有足够大固有晶格失配的TMD异质双层中，局部膨胀/收缩成为主要的重构途径，旋转贡献随转角增大而出现。 3. hBN封装通过抑制面外起伏，显著增强并局域化了面内重构，突显了样品设计与基底效应在调控莫尔系统性质中的关键作用。 4. 外部异质应变不能在同质双层中诱发膨胀重构，但能有效操控和重分布由旋转产生的剪切应变场，为应变工程调控莫尔电子态提供了实验依据。
本研究的科学价值在于，将莫尔重构的研究从定性观察和理论模拟推进到了定量测量和机制区分的阶段，为理解莫尔超晶格中晶格-电子耦合提供了坚实的实验基础。所揭示的不同重构机制及其竞争关系，对于解释和预测不同TMD莫尔体系（如同质 vs. 异质、不同转角、不同封装条件）的独特光电行为至关重要。
应用价值 体现在为“莫尔工程”提供了新的调控维度：通过选择材料组合（控制晶格失配）、设计转角、施加封装或外部应变，可以主动设计所需的晶格形变模式（旋转或膨胀主导的应变场），进而定制莫尔势的深度、形状以及电子能带的平坦度，这对于实现拓扑物态、关联电子相和新型激子器件具有重要意义。
六、 研究亮点
七、 其他有价值内容
本研究建立的实验与分析方法具有很好的普适性和扩展性。文中指出，该方法不仅适用于本文研究的MoS₂和WSe₂，也应广泛适用于其他h相TMDs（如MoTe₂, WS₂等），甚至可能推广到磁性二维莫尔超晶格（如CrI₃双层）。作为一种基于衍射的成像技术，布拉格干涉法对悬浮和封装结构均兼容，未来有望应用于更复杂的垂直异质结构（如包含栅电极的器件），实现晶格重构与（光）电学性能的直接关联测量，从而在原子尺度上建立结构-物性关联，推动二维莫尔材料的基础研究与器件应用。