在物理期刊《Physical Review B》第103卷第L121102期，发表于2021年3月1日，来自印度科学研究院物理系凝聚态理论中心的Indrajit Maity、Prabal K. Maiti、H. R. Krishnamurthy以及Manish Jain（通讯作者）共同发表了一项关于扭转过渡金属二硫族化合物（Transition Metal Dichalcogenides, TMDs）双层异质结中莫尔晶格（moiré lattice）重构的开创性研究。这项研究挑战了过去对结构弛豫效应的普遍假设，揭示了在大扭转角度下莫尔超晶格会发生剧烈的重构，并显著改变体系的电子能带结构，为理解和设计基于莫尔材料的量子现象提供了新的物理图像和理论基础。
这项研究的学术背景根植于近年来二维材料莫尔超晶格这一前沿领域。当两层二维材料（如石墨烯或TMDs）以微小角度或特定角度堆叠时，会形成周期远大于原子晶格的莫尔图案。这种莫尔超晶格能产生平带（flat bands），是导致关联电子态、超导、莫特绝缘体等奇异物理现象的关键。为了精确计算这些电子性质，必须考虑原子尺度的结构弛豫效应，即原子位置会从理论上的刚性扭转构型发生偏移以降低系统总能量。然而，以往几乎所有的研究都默认一个基本前提：即使经过弛豫，莫尔晶格的晶格常数（即莫尔超晶胞的尺寸）与未弛豫的刚性扭转结构保持一致。Maity等人的工作正是针对这一关键预设，通过引入全局优化的模拟退火（Simulated Annealing, SA）方法，发现在特定角度区间，这一预设会被打破，发生根本性的“晶格重构”，从而深刻影响材料的物理性质。因此，本研究的目标是系统探究大扭转角度下TMD双层体系的晶格重构现象、其物理起源，以及重构对电子结构的重大影响。
研究工作的流程严谨而系统，主要包含以下几个步骤。首先，研究对象是四种扭转双层TMDs：MoS2、MoSe2、WSe2和WS2。研究聚焦于扭转角（θ）接近0°和60°的区域，因为这两个区域存在不同的高对称堆叠构型。研究构建了包含成千上万个原子的大型莫尔超晶胞进行计算模拟。其次，在研究方法和软件上，作者使用了他们之前开发的Twister代码来构建初始的扭转双层结构。随后，他们采用了两种截然不同的弛豫方法来获得能量最低的结构：第一种是标准的局域能量最小化方法（Standard Relaxation, SR），即从初始构型出发，在势能面上寻找最近的局域极小值点；第二种是模拟退火方法（SA），该方法通过经典的分子动力学模拟，在1 K温度下对系统进行采样，然后冷却至0 K再进行能量最小化。SA方法的优势在于能够帮助系统跨越势能垒，有更高概率找到全局或更优的极小值，这对于具有指数级多个局域极小值的复杂势能面至关重要。分子动力学和能量最小化计算通过LAMMPS软件完成。体系中，层内原子相互作用由Stillinger-Weber势描述，而层间的范德华相互作用则由Kolmogorov-Crespi (KC)势描述，该参数已被证明能准确捕捉TMDs的层间作用。此外，作者还使用修改版的phonopy代码计算了声子谱，并使用基于密度泛函理论（DFT）的SIESTA软件包进行了电子结构计算，以验证重构对电子性质的影响。这些计算在超级计算机上完成，确保了大规模模拟的可行性。
研究取得了一系列重要且相互关联的结果，可以概括为以下几个关键发现。第一，他们成功观测并定义了“晶格重构”。对于扭转角接近0°的情况（例如θ=1°），SR和SA方法得到的结构差异不大，都保持了未弛豫结构的莫尔晶格常数，仅在一些能量不利的AA堆叠处，SA方法观察到了微小的畴壁（Domain Wall, DW）弯曲。然而，当扭转角接近60°时，行为出现了戏剧性的分区。他们将θ→60°的区域分为两个区间：区域I（θ < 58.3°）：SR和SA得到的结构均未改变莫尔晶格常数，畴壁略有弯曲。区域II（θ ≥ 58.4°）：情况发生根本改变。使用SR方法时，能量最有利的AA’堆叠区域面积显著增大，从类似于勒洛三角形（Reuleaux triangle）的形状演变为近似的六边形，连接能量最不利的A’B堆叠的畴壁变得显著弯曲。然而，使用SA方法进行全局优化后，结构发生了剧烈的“重构”：三个A’B堆叠会聚拢形成“三聚体”（trimer）结构单元，这些结构单元本身排列成新的三角形晶格。最重要的是，重构后的莫尔晶格常数（a^SA_m）变成了未弛豫晶格常数（a_m）的√3倍（a^SA_m = √3 a_m）。这明确证明，在大扭转角下，弛豫后的莫尔晶格常数可以与其刚性扭转的原始值不同，即发生了“晶格重构”。通过计算莫尔尺度（以AA/A‘B堆叠点代表莫尔格点）的径向分布函数gm®，他们确认了在区域II，SA得到的结构中A’B堆叠形成了最近邻为2的三角形晶格，而SR结构则保持了六角网络。
第二，他们深入探讨了晶格重构的几何与能量起源。通过分析SR结构中能量最有利的AA’区域几何形状，他们引入了一个关键几何参数：垂直平分线长度（p）与莫尔晶格常数（a_m）的比值（p/a_m）。研究发现，当p/a_m ≥ 1（对应θ约58.5°）时，AA’区域呈现为勒洛三角形，其周长与面积之比是局域最大值。由于畴壁的能量高于AA’区域，这种几何形状在能量上是不利的，从而驱动了重构以降低总能量。从能量角度分析，系统的总能量是层内（应变和弯曲）能量与层间能量的总和。对于接近60°的扭转角，层间能量近似正比于畴壁的总长度L。在区域II，SR方法得到的畴壁是显著弯曲的，其长度L_sr可以通过重构被最小化，使得畴壁在重构结构中变得近乎直线（L_sa < L_sr）。能量计算证实，在区域II，SA方法相比SR方法带来的能量增益远大于接近0°的情况，这主要归功于晶格重构导致的畴壁长度缩短。而在接近0°时，SR得到的畴壁原本就是直线，SA带来的能量增益主要源自畴壁在AA堆叠处的微小弯曲（一种屈曲失稳），与晶格重构无关。声子谱计算进一步支持了重构的机制：在θ=58.5°时，发现一个声子模式显著软化并变得几乎无色散，表现出类似零模的特征，这种软模通常会诱发晶格重构。
第三，研究揭示了晶格重构对电子能带结构的重大影响。作者以包含近2.5万个原子的MoS2扭转双层（θ=58.47°）为对象，对比了SR和SA弛豫后结构的电子能带。结果显示，重构不仅使带隙增加了约20 meV，还显著改变了能带边缘附近的能级间隔。更重要的是，尽管两种结构在价带顶（Valence Band Maximum, VBM）和导带底（Conduction Band Minimum, CBM）附近都产生了超平带（带宽小于1 meV），但对应这些平带的电子波函数局域化位置截然不同。在SR（未重构）结构中，VBM附近的电子态主要局域在单个AA’堆叠区域，其势阱形状近似等边三角形。而在SA（重构）结构中，由于三聚体的形成，被三聚体包围的AA’区域（面积较小且固定）与其他AA’区域（面积较大）的局域环境不同，导致等边三角形势阱的简并被解除。计算发现，VBM附近的前三个能带的波函数强烈局域于被三聚体包围的AA’区域，而第四个能带则局域于更大的AA’区域。CBM附近的电子态则主要局域在AB’堆叠区域。这种电子和空穴的空间分离源于面内应变驱动的莫尔势（moiré potential）在重构后的重新分布。重构后，被三聚体包围的AA’区域莫尔势垒高度最大（约171 meV），而其他AA’区域的势垒则明显降低（约98 meV）。AB’区域的势阱深度也发生了变化。这些发现为实验观察到的在AB’堆叠处存在大的隧穿电流提供了理论解释。
基于上述结果，本研究得出的核心结论是：对于扭转双层TMDs，当扭转角大于约58.5°时，会发生根本性的莫尔晶格重构，导致莫尔晶格常数发生√3倍的扩大，并形成以A‘B堆叠三聚体为基本单元的三角形晶格。这一重构现象源于在特定几何条件下（p/a_m ≥ 1），系统通过畴壁重排来最小化总能量的内在驱动。重构并非简单的局域原子弛豫，而是一种在莫尔尺度上涌现的新的有序现象。更重要的是，这种晶格重构会显著改变体系的电子结构，特别是调控平带电子的空间局域位置和能级，为在莫尔体系中实现和调控关联电子态、激子束缚等提供了全新的自由度。
本研究的科学价值和应用意义重大。在科学价值层面，它首次明确揭示了在特定扭转角区间，莫尔晶格的全局结构可能完全不同于基于刚性模型或标准弛豫的预期，这修正了该领域的一个基本认知，强调在理论研究和高精度计算中必须考虑这种大尺度的重构效应。其次，它建立了几何参数（p/a_m）作为预测重构发生的判据，并阐明了其能量学和软模动力学的微观机制。第三，它建立了从结构重构到电子性质（特别是平带局域化）的清晰物理图像，为理解实验现象（如复杂畴壁网络、空间不均匀性）提供了理论基础。在应用价值上，这项工作为“莫尔工程”提供了新的思路：通过精确控制扭转角度，不仅可以调节莫尔周期，还可以触发晶格重构这一二级相变，从而在更大范围内调节电子态、光学响应和振动性质，为设计新型量子材料和光电器件开辟了新途径。
本研究的亮点在于：第一，重要的新发现：首次报道了扭转双层TMDs中莫尔晶格常数的重构现象，打破了该领域的普遍预设。第二，方法的创新性：通过对比标准局域弛豫（SR）和全局优化的模拟退火（SA）方法，凸显了在复杂势能面中寻找全局最小结构的重要性，该方法具有普适意义。第三，物理图像的完整性：研究从结构表征（层间距图、径向分布函数）、几何分析（p/a_m比）、能量学（畴壁长度与总能量）、动力学（声子软模）到最终的电子性质（DFT能带计算、波函数局域化），构成了一个逻辑严密、证据链完整的物理故事。第四，研究对象的系统性：不仅在MoS2中进行了详细研究，还验证了结论在MoSe2、WSe2、WS2等其他TMDs中的普适性。
此外，文中还提到了一些有价值的延伸内容。例如，在模拟过程中观察到了如扭曲六边形、Kagome晶格等亚稳态过渡结构，并指出在包含更多自由度的超大超胞模拟中，可能存在由三聚体以上单元构成的、更复杂的熵稳定重构结构。这些结构可能因衬底、应变等外部因素而稳定存在，这有助于解释实验中观察到的非均匀六边形畴壁网络和“呼吸”模式。研究也暗示，这种由几何约束驱动的重构机制可能普遍存在于其他具有畴壁网络的莫尔材料中。最后，作者利用DFT对较小体系进行了完全弛豫计算，并证明在施加小的压缩应变后，SA得到的重构结构确实比SR结构更稳定，进一步验证了结论的可靠性。