关于三维重构WSe2/WS2莫尔超晶格中莫尔平带成像的研究报告
本研究由来自美国加州大学伯克利分校、亚利桑那州立大学以及日本国家材料科学研究所等多个机构的科研团队共同完成，于2021年7月发表在《自然·材料》（*Nature Materials*）期刊上，标题为“在三维重构的WSe2/WS2超晶格中成像莫尔平带”。
一、 学术背景与研究动机
本研究属于凝聚态物理和二维材料科学领域，聚焦于过渡金属二硫化物（Transition Metal Dichalcogenide, TMD）异质结中的莫尔超晶格。近年来，TMD异质结中的莫尔超晶格因其能形成极窄的“平带”电子能带而成为研究强关联量子现象（如关联绝缘体、维格纳晶体等）的新兴平台。平带意味着电子的动能极小，从而电子间的库仑相互作用主导系统行为，为探索新颖的量子态提供了机会。
然而，对TMD莫尔超晶格中平带性质的深入理解和精确调控，严重依赖于对原子尺度下莫尔超晶格结构的微观认知，特别是三维结构重构（如面外翘曲和面内应变重分布）如何影响电子能带。此前的研究大多基于简化的理论模型，这些模型未能充分计入复杂的三维重构效应，导致对平带空间分布和能量特性的预测可能与实际情况存在偏差。因此，该研究的核心目标是：通过结合高精度的实验测量与第一性原理计算，定量揭示WSe2/WS2莫尔异质结在原子尺度下的三维重构细节，并阐明这种重构如何决定其关键的低能量莫尔平带特性。
二、 详细研究流程与方法
本研究整合了先进的实验表征与大规模理论模拟，主要流程包括器件制备与表征、三维结构探测与重构分析、电子态谱学测量以及第一性原理计算验证。
高质量异质结器件的制备与初步表征：
- 研究对象与制备： 研究团队使用机械剥离和干法转移技术，制备了高质量的WSe2/WS2异质双层结构。与化学气相沉积（CVD）生长的样品相比，这种技术能更好地控制堆叠角度（本研究为近零度对齐）和界面质量。器件以六方氮化硼（hBN）为衬底，并创新性地使用石墨烯纳米带阵列作为顶部接触电极，硅衬底作为背栅以调控载流子浓度。这种设计使得样品表面既有暴露的TMD区域，也有被单层石墨烯覆盖的区域，这对于后续区分几何形貌和电子态效应至关重要。
- 表征手段： 利用原子力显微镜（AFM）在室温大气环境下对器件表面形貌进行了初步检查，确认了异质结构和石墨烯电极的分布。随后，在超低温（5.4 K）、超高真空环境下，使用扫描隧道显微镜（STM）对样品进行原子级成像，清晰观测到了周期约为8.16纳米的莫尔超晶格图案，证明了器件的高质量。
莫尔超晶格三维重构的揭示：
- 实验观测与矛盾： 高分辨STM形貌图显示，在每个莫尔晶胞内，存在一个显著的“谷”被六个“峰”环绕的蜂窝状高度调制。为了排除电子态密度（LDOS）变化对表观高度的影响（这是一个常见挑战），研究者巧妙地利用了石墨烯覆盖区域。在TMD带隙内的偏压下扫描石墨烯覆盖区域，由于石墨烯LDOS相对均匀，测得的形貌更真实地反映了底层TMD的几何起伏。结果证实，WSe2/WS2异质结确实呈现“六峰环绕一谷”的形貌。
- 理论模拟与突破： 传统的理解认为，不同堆叠构型（如AA堆叠、AB堆叠）的层间距差异是高度调制的主因，AA堆叠因原子位阻应有最大层间距（即表现为“峰”）。然而，仅考虑垂直松弛的力场模拟结果（预测AA为峰，AB为谷）与实验观察完全相反。
- 第一性原理计算揭示三维重构： 为了解决这一矛盾，团队进行了大规模的第一性原理计算。模拟首次全面揭示了WSe2/WS2莫尔超晶格存在强烈的三维重构：包括显著的面内应变重分布和巨大的面外翘曲。具体而言，由于WSe2（晶格常数较大）和WS2（晶格常数较小）之间的晶格失配，系统倾向于扩大能量更低的AB堆叠区域。这导致WSe2层在AB区域受到压缩，而残余的张应变局域在AA区域。为了部分释放应变，异质双层发生了同相的面外翘曲。计算得到的面外高度分布图完美复现了STM观测到的“六峰环绕一谷”结构，其中AB堆叠位置（BSe/W和BW/S）因翘曲而抬起成为“峰”，AA堆叠位置成为“谷”。理论计算还表明，hBN衬底的存在只会轻微减弱这种翘曲效应。
莫尔平带电子结构的扫描隧道谱学（STS）研究：
- 空间分辨谱学测量： 研究团队在超低温下对异质结不同莫尔位点（AA, BSe/W, BW/S）进行了扫描隧道谱（dI/dV谱）测量，获取了局域电子态密度信息。
- K点与Γ点平带的区分： 区分来源于单层布里渊区K点和Γ点的电子态是关键挑战。研究者利用了两个特征：
- 衰减常数差异： K点态具有更大的面内动量，其波函数在TMD层外衰减更快。通过测量不同探针-样品间距下的dI/dV谱，发现约-1.5 V处的峰比约-1.7 V处的峰对间距更敏感，表明前者源于K点。
- 原子尺度干涉图案： K点态主要由钨（W）的d轨道（角动量m=±2）贡献，会在相邻W原子上产生2π/3的相位 winding，导致原子尺度的建设性和破坏性干涉图案。高分辨dI/dV mapping显示，-1.5 V处的信号在原子尺度上振荡，与WSe2晶格匹配，而-1.7 V处的信号变化平滑，从而确认-1.5 V峰来自K点价带边态，-1.7 V峰来自Γ点态。
- 平带空间分布成像： 通过在大范围区域进行固定能量下的dI/dV mapping，直接可视化了不同能量平带态的空间局域情况。研究发现，位于价带顶（约-1.52 V，标记为K1）的K点莫尔平带高度局域在BSe/W位点。该平带极其狭窄，半高宽经解卷积后仅为10 meV ± 1 meV。在略低的能量（约-1.59 V，标记为K2），态密度分布发生变化，在BSe/W位点中心出现节点，形成环状分布。类似的局域行为也在Γ点平带（γ1局域于BSe/W，γ2在BSe/W出现节点）中被观察到。
第一性原理计算验证与物理机制阐释：
- 能带结构与态密度计算： 基于力场重构后的莫尔超晶格原子结构，研究者进行了大规模密度泛函理论（DFT）计算。计算得到的价带迷你能带结构显示，最顶部的价带（源自K点）宽度仅约10 meV，且与下一个能带间隔约30 meV，与STS测量的窄平带特征定量吻合。
- 理论态密度分布与实验对比： 计算了不同能量区间（E1-E4，对应实验中的K1, K2, γ1, γ2）的局域态密度空间分布。结果显示，E1（顶部K点带）的态密度强烈局域在BSe/W区域；E2（次深K点带）在BSe/W中心出现节点；E3（顶部Γ点带）局域在BSe/W；E4（次深Γ点带）在BSe/W出现节点。这些计算结果与实验dI/dV mapping图像高度一致。
- 平带起源分析： 计算进一步揭示了K点平带主要源于单层WSe2在三维重构下产生的形变（应变和翘曲），而非层间杂化势的调制。相反，Γ点平带则主要来源于莫尔超晶格内不均匀的层间杂化。这一发现突破了传统上认为层间相互作用变化是莫尔能带调制主导机制的认知。
三、 主要研究结果
- 实验直接观测到原子级三维重构： 结合STM形貌测量和第一性原理计算，首次定量揭示了WSe2/WS2莫尔异质结存在强烈的面内应变重分布和面外翘曲，形成了“六峰（AB位）环绕一谷（AA位）”的独特三维重构结构。
- 发现并表征了极其狭窄的K点莫尔平带： STS测量直接观测到位于价带顶、宽度仅约10 meV的K点莫尔平带。该平带在能量上与更深能带隔离约50 meV，且空间上高度局域于BSe/W位点。
- 绘制了多组平带的空间分布图： 实验观测到不止一个平带，包括来自K点和Γ点的多个平带。它们的空间局域特性各不相同，并随能量变化，例如从在BSe/W位点的强局域峰变为围绕该位点的环状分布或节点。
- 理论与实验高度定量一致： 基于三维重构结构的第一性原理计算，不仅复现了实验观测的形貌，其预测的平带能量分离（K1与γ1约180 meV vs. 实验205 meV）、带宽（~10 meV）以及各平带态密度的空间分布模式，均与STS实验结果高度吻合。
- 阐明了三维重构的主导作用： 研究明确指出，传统的、忽略三维重构的简化模型（预测AA位点局域）与实验矛盾。而全面的第一性原理计算表明，应变重分布和三维翘曲是决定WSe2/WS2异质结最低能量莫尔平带（尤其是K点平带）的关键因素，其影响力超过了单纯的层间相互作用变化。
四、 研究结论与意义
本研究的核心结论是：在WSe2/WS2莫尔超晶格中，三维结构重构（应变和翘曲）主导了低能电子结构，导致在价带顶产生了一个宽度仅为10 meV且高度局域的K点莫尔平带。这一发现纠正了先前简化理论模型的认知。
其科学价值在于： * 提供了关键微观认知： 首次在原子尺度和meV能量尺度上，定量揭示了TMD莫尔超晶格的三维原子重构及其对平带电子结构的决定性影响，填补了该领域长期缺失的关键微观知识。 * 为强关联物理研究奠定基础： 对莫尔平带宽、空间局域性和起源的精确测定，为在TMD莫尔体系中模拟哈伯德模型、理解和调控新奇的强关联量子现象（如广义维格纳晶体、关联绝缘体等）提供了至关重要的输入参数和物理图像。 * 展示了先进研究方法论： 本研究成功结合了高质量器件制备、高分辨扫描探针技术以及大规模第一性原理计算，为研究其他二维材料异质结的莫尔物理树立了方法论范例。
五、 研究亮点
- 突破性发现： 直接实验观测并理论证实了TMD莫尔超晶格中强烈的三维原子重构，并发现其对最关键的平带特性起主导作用。
- 技术创新： 利用石墨烯覆盖层区分几何与电子效应对STM形貌的影响，巧妙解决了STM测量中的固有难题。结合原子尺度干涉图案和衰减常数差异，可靠地区分了K点和Γ点电子态。
- 理论与实验深度交融： 研究并非简单对比，而是通过实验发现与简化理论的矛盾，驱动了包含全三维重构的第一性原理计算，最终实现了从原子结构到电子性质的定量、自洽解释，形成了完整的证据链。
- 研究对象的精确性： 使用高质量、近零扭角的剥离型异质结器件，并采用绝缘hBN衬底和栅极控制，最大程度减少了基底干扰和掺杂不确定性，获得了本征的物理信息。
六、 其他有价值的内容
研究还指出，K点平带主要源于单层的形变，而Γ点平带源于层间杂化，这深化了对不同对称性轨道电子态在莫尔势中响应机制的理解。此外，观测到的平带空间分布模式（如基态局域、激发态出现节点或环状分布）类似于量子势阱中的波函数特征，为理解莫尔势阱的形状和强度提供了直观图像。这些细节对于未来设计通过应变或电场调控莫尔平带和关联物态具有重要意义。