关于二维非公度异质结弛豫与畴形成的学术研究报告

本文旨在介绍一篇发表于《物理评论B》期刊上的原创性研究论文。该研究由来自哈佛大学、明尼苏达大学和堪萨斯大学的研究团队共同完成，并于2018年12月5日正式发表。论文题为《非公度二维异质结中的弛豫与畴形成》，由Stephen Carr、Daniel Massatt、Steven B. Torrisi、Paul Cazeaux、Mitchell Luskin和Efthimios Kaxiras合作撰写。这项研究属于凝聚态物理和材料科学领域，特别是二维材料堆叠结构与力学性质的研究前沿。

一、 研究背景与目标

二维层状材料，如石墨烯、二硫化钼等，因其独特的电学、光学和机械性质而备受关注。当两层这样的材料以微小的扭转角或晶格常数失配堆叠在一起时，会形成周期远大于原子间距的莫尔图案。在这种“非公度”或“近公度”的堆叠中，层间相互作用会诱导原子位置发生弛豫，以最小化系统的总能量。这种弛豫过程常常导致均匀堆叠区域（畴）的形成，畴与畴之间由狭窄的应变孤子（即畴壁）分隔。这种现象在实验上已被观测到，例如在小角扭转双层石墨烯中观察到的规则三角畴网络。理解并预测这种弛豫行为对于调控二维异质结的电子性质至关重要，因为畴壁可能承载拓扑保护的边缘态，在量子信息等领域具有潜在应用价值。

然而，现有的理论模型，特别是基于实空间的连续介质模型，在处理严格非公度的双层或多层体系时面临挑战。对于多层体系，各层两两之间产生的莫尔周期可能互不匹配，导致整个体系在实空间中缺乏一个公共的周期性超胞，这使得传统的实空间连续介质方法难以直接推广。本研究的主要目标，正是为了克服这一困难，提出一个全新的、基于“构型空间”的普适性理论框架，用于高效、准确地计算非公度二维异质结的力学弛豫模式，并深入阐释畴形成的物理机制。

二、 研究方法与工作流程

本研究主要包含三个核心部分：理论框架的构建、关键能量参数的从头算获取、以及针对具体材料体系的弛豫计算与结果分析。研究未涉及传统意义上的实验样本，其“研究对象”是特定的二维材料体系及其理论模型。

1. 理论框架：构型空间方法的引入 这是本研究的核心创新点。研究团队摒弃了直接在实空间描述原子位移的传统思路，转而引入“构型空间”这一概念。构型空间的核心思想是：将双层结构中每一个原子位置所对应的局部堆垛环境（即两层原子间的相对位移矢量）作为基本变量。对于非公度堆叠，尽管原子在实空间中的排布是非周期的，但所有可能的局部堆垛构型却密集地分布在一个紧凑的区域内——即底层晶格的单位晶胞所定义的环面（二维环面可理解为将一个矩形对边连接形成的面包圈形状）。通过一个由晶格失配矩阵定义的线性映射，可以将实空间中的位置映射到构型空间中的一个点。

基于此，研究将系统的总能量（包括层内弹性能和层间堆垛能）在构型空间中进行重新表述。层内弹性能采用线性各向同性连续介质弹性理论描述，依赖于材料的体弹性模量和剪切模量。层间相互作用则通过“广义层错能”曲面来描述，该曲面给出了单位晶胞内任意相对堆垛构型所对应的能量。最终，弛豫问题转化为在构型空间的环面上寻找一个平滑的位移场u(b)，使得总能量泛函取最小值。通过傅里叶级数展开和标准优化算法，可以高效地求解这个位移场。

该方法的优势在于：首先，构型空间本身是周期性的，无论实空间中的堆叠是否具有周期性，这为处理严格非公度体系提供了自然的数学框架。其次，它将弛豫问题与材料的本征性质（弹性模量和GSFE）直接关联，无需依赖经验性的原子间势函数，这对于缺乏可靠经验势的材料（如二硫化钼）尤为重要。最后，该方法计算效率高，能够处理非常小的扭转角情况。

2. 材料参数的第一性原理计算 为了应用上述理论模型，需要获取具体材料的弹性模量和广义层错能（GSFE）。研究选取了三个代表性体系进行深入分析：小角扭转双层石墨烯，以及二硫化钼（MoS2）的两种高对称堆叠构型（0°和180°扭转）。

• 弹性模量计算：对单层材料施加各向同性和各向异性的面内应变（在x和y方向施加±1.5%的应变，步长0.3%），基于密度泛函理论（DFT）计算其基态能量随应变的变化。通过对能量-应变关系进行线性拟合，得到材料的体弹性模量k和剪切模量g。
• 广义层错能计算：对于双层体系，固定所有原子的面内位置，允许其垂直方向弛豫，使用包含范德华力修正（vdW-DFT方法，具体采用SCAN+rVV10泛函）的DFT计算。将顶层相对于底层在单位晶胞内进行平移（在一个9x9的网格上采样），计算每个堆垛构型下的能量，从而获得GSFE能量曲面。随后，利用一个基于对称性考虑的傅里叶级数形式对离散的能量点进行拟合，得到解析的V_GSFE(b)函数表达式。论文中提供了石墨烯和两种MoS2堆叠的弹性模量及GSFE傅里叶系数表格。

3. 弛豫计算与结果分析 在获得材料参数后，研究团队将构型空间方法应用于具体体系。他们重点研究了小扭转角（例如3°）下的弛豫行为。为了更清晰地展示畴形成现象，在部分计算中人为地将层间耦合能放大了100倍以增强视觉效果。通过求解构型空间中的能量最小化问题，得到了最优的位移场u(b)。然后，利用映射关系将构型空间中的位移场转换回实空间，从而获得弛豫后的原子结构。

三、 主要研究结果

1. 构型空间中的弛豫图案 计算结果显示，在构型空间中，弛豫导致低能量堆垛区域（对应于AB/BA堆垛或特定低能堆垛）扩大，而高能量堆垛区域（如AA堆垛）收缩，形成由细线（对应畴壁）和节点（对应高能堆垛点）连接的图案。这一现象直观地揭示了弛豫的驱动力：系统通过原子位移，尽可能增加低能构型的占比，同时将高能构型压缩到最小的区域内。

2. 实空间中的畴结构 将构型空间的解映射回实空间后，可以清晰地看到弛豫形成的畴结构： * 对于石墨烯和0° MoS2：它们的GSFE曲面具有两个简并的最低能量点（AB和BA堆垛）。弛豫后，系统形成由AB和BA两种堆垛构成的三角形畴交替排列的规则图案，畴之间由应变孤子（畴壁）分隔。这与实验上在扭转双层石墨烯中观察到的结果一致。 * 对于180° MoS2：其GSFE曲面只有一个最低能量点。弛豫后，低能堆垛区域扩张，形成六边形的畴，畴壁网络呈蜂窝状。

3. 局域扭转角的变化 研究还发现，弛豫位移场的旋度与构型空间中位移场的散度近似成正比。计算表明，在低能量堆垛畴内，位移场的散度为负，意味着这些区域的局域扭转角减小（“解扭转”效应）；而在高能量堆垛节点处，散度为正，意味着局域扭转角增大。例如，在小角扭转双层石墨烯的AA堆垛节点处，局域扭转角会收敛到一个固定值（如文中指出，当全局扭转角趋近于0°时，AA节点的局域角约为1.7°），这与另一项实空间研究的结果相符。

4. 畴尺寸与扭转角的关系 研究计算了不同扭转角下的弛豫结构。结果表明，随着全局扭转角的减小，构型空间中的畴壁和节点尺寸按比例缩小。然而，由于从构型空间到实空间的映射包含一个与扭转角θ成反比的因子，导致应变孤子在实空间中的宽度在足够小的θ下几乎保持不变，而畴的尺寸（孤子之间的间距）则随θ减小而增大。这解释了实验中观察到的现象：扭转角越小，形成的莫尔超胞越大，畴的尺寸也越大，而畴壁的宽度基本不变。

5. 方法的有效性与推广性 研究证实，基于第一性原理GSFE的构型空间方法能够准确计算MoS2等材料在小扭转角下的弛豫，而无需经验势。该方法自然地解释了畴形成的几何图案（由GSFE的临界点数量和性质决定）和尺寸缩放规律。论文还简要讨论了该方法向多层体系推广的可能性，指出对于p层体系，构型空间将升维为一个2(p-1)维的环面，尽管应变能的表述会变得复杂，但理论框架在原理上是可扩展的。

四、 研究结论与意义

本研究成功发展了一种基于“构型空间”的新理论方法，用于研究非公度二维异质结的力学弛豫。研究得出结论：当两层晶格近乎对齐时，系统会通过原子弛豫形成大面积的均匀堆垛畴，畴间由狭窄的应变孤子壁分隔。畴的几何形状由层间相互作用能（GSFE）的极小值点决定（例如，两个简并极小值导致三角畴，单个极小值导致六边畴），而畴的尺寸则与扭转角成反比缩放。该方法不仅为理解小角扭转双层石墨烯等已知体系提供了新的视角，更重要的是，为预测和计算缺乏经验原子势的复杂二维材料异质结（如过渡金属硫族化合物）的弛豫行为提供了强大且高效的工具。

五、 研究亮点

1. 理论创新：首次将“构型空间”概念系统性地引入二维非公度异质结弛豫问题的研究，建立了全新的理论框架，突破了传统实空间连续介质模型在处理严格非公度及多层体系时的局限性。
2. 方法普适：该方法不依赖于经验原子势，而是直接采用基于第一性原理计算的广义层错能和弹性模量，因此具有很高的准确性和普适性，可应用于广泛的二维材料体系。
3. 物理图像清晰：在构型空间中，弛豫表现为“能量从高向低流动”的直观图像，深刻揭示了畴形成和应变孤子产生的物理本质。同时，建立了局域扭转角变化与构型空间位移场散度的直接联系。
4. 计算高效：通过在紧凑的周期域（构型空间）内求解，能够以极小的计算成本获得高精度的弛豫结构，特别是对于非常小的扭转角情况优势明显。
5. 预测性强：研究明确指出了畴结构对GSFE临界点的依赖关系，以及畴尺寸与扭转角的标度律，这对通过调控堆叠角度来设计特定畴图案和电子性质具有指导意义。

六、 其他有价值的内容

论文还探讨了该方法与经典Frenkel-Kontorova模型（描述一维公度-非公度转变）的潜在联系，并指出在极端小角下可能发生的公度-非公度相变是未来值得深入研究的方向。此外，作者强调了该构型空间方法为后续计算弛豫结构的电子性质（如拓扑边缘态）提供了准确的原子位置输入，架起了结构弛豫与电子性质预测之间的桥梁。文末提及的多层体系推广展望，也为更复杂的范德华异质结研究开辟了新的路径。