本研究由宾夕法尼亚大学物理与天文系的 Christophe De Beule、E. J. Mele，新加坡国立大学先进二维材料中心的 Gayani N. Pallewela、Shaffique Adam，哈佛大学工程与应用科学学院的 Mohammed M. Al Ezzi，以及华盛顿大学圣路易斯分校物理系的 Liangtao Peng 等人合作完成，作为 Letter 发表于 Physical Review B 113, L241402 (2026)，并被编辑推荐（Editors’ Suggestion）。该研究属于凝聚态物理与二维材料领域，聚焦于小扭转角下莫尔（moiré）超晶格中的晶格弛豫（lattice relaxation）问题。
莫尔范德瓦尔斯材料通过晶格失配或层间扭转形成远大于原子晶格常数的超胞，从而实现能带平坦化并诱发超导、关联绝缘态、轨道陈绝缘体等丰富量子物态。在刚性晶格近似下，莫尔结构仅由几何扭转角决定；但实验和理论均表明，二维材料层更像弹性膜，层间范德瓦尔斯堆垛能与层内弹性形变能之间的竞争会导致显著的原子重构，即晶格弛豫。在较大扭转角下，弹性能占主导，晶格弛豫为弱耦合微扰；而随着扭转角减小，莫尔超胞面积增大，粘附能增益超过弹性惩罚，体系进入所谓的“边际扭转”（marginal twist）强耦合区域。此时，能量有利堆垛区域扩大并覆盖大部分莫尔原胞，不利堆垛区域收缩为拓扑缺陷，形成由畴壁（domain wall）构成的周期性网络。此前对该强耦合区域的理论描述主要依赖密度泛函理论、分子动力学模拟或连续弹性理论数值求解，缺乏解析理论。
该论文的核心目标是建立边际扭转双层体系中晶格弛豫的解析理论，涵盖三角孤子网络（triangular soliton network）和蜂窝孤子网络（honeycomb soliton network）两种拓扑结构。理论框架基于连续弹性模型，以层间位移场 φ® = mr + u® 为基本自由度，其中 m 描述刚性莫尔结构，u 为声学位移。体系能量包含弹性项和粘附项，弹性项由单层 Lamé 常数 μ 和 λ 表征，粘附项由堆垛能 v(φ) 给出。对扭转莫尔体系，耦合常数 η = (l²/a²)v₁/μ ≈ c₁/θ² 控制弛豫强度，其中 l 为莫尔周期，a 为晶格常数，v₁ 为粘附势的主导傅里叶分量。|η| ≪ 1 对应弱耦合大角度区域，|η| ≫ 1 对应强耦合小角度区域。对于具有 D₆ 堆垛对称性的双层体系，当 v₁ > 0 时，AB/BA 堆垛为两个简并能量极小值，AA 堆垛为能量极大值，体系形成三角孤子网络；当 v₁ < 0 时，AA 堆垛为唯一能量极小值，体系形成蜂窝孤子网络。
在研究流程上，作者首先将位移场分解为纵向分量 u‖_g 和横向分量 u⊥g，两者在弹性能量中解耦。由于扭转莫尔主要依靠“局部扭转”（twirling）弛豫，横向分量占主导，因此理论分析设定 λ = 0 以简化而不改变核心物理。对于三角网络，孤子为 AB/SP/BA 部分位错，位移场的解析解为一系列横向傅里叶分量的叠加，形式为 u® = (√3a/2π)∑ᵢ ẑ × ĝᵢ ∑ₙ (u⊥{n,0}/n) sin(ngᵢ·r)。在 η → ∞ 极限下，u⊥{n,0} → 2/3，位移场收敛为具有尖锐畴壁的解析表达式。对于有限大 η，畴壁具有有限宽度，作者通过将畴壁斜率与孤立 AB/SP/BA 孤子的斜率匹配，得到 u⊥{n,0}(θ) = (2⁄3)sech²(2nθ/3√c₁)。该函数随 n 指数衰减，求和收敛于 n_max ≈ η。若粘附势包含高阶星形傅里叶分量，则畴壁斜率修正为 2√(c₁ − 8c₂ + 9c₃)，理论可通过等效参数 c₁,eff = c₁ − 8c₂ + 9c₃ 纳入这些贡献。
对于蜂窝网络，孤子为 AA/SP/AA 全位错，由于网络非本原，解更为复杂。数值分析表明，主导贡献除 u⊥{n,0} 外，还包括三组横向分量：u⊥{2n,n}、u⊥{2n+1,n} 和 u⊥{2n+2,n}。作者通过拟合连续弹性数值解，给出了这些分量的经验解析形式 u⊥{2n+i,n} = (aᵢ + bᵢ/√|η|)sech²[n(rᵢ + sᵢ/√|η|)]，并列表给出拟合参数。同时，u⊥{n,0} 表现为振荡衰减形式 −0.26 sin(2.1n)csch(0.27n)。
在结果验证方面，作者将解析理论与 LAMMPS 分子动力学模拟及连续弹性数值解进行了系统对比。图 1 展示了 AA 堆垛 WSe₂ 在 θ = 0.25° 时的局部扭转角分布，理论预测与两种数值方法高度一致。图 2 展示了不同 η 值下各壳层横向傅里叶分量的演化，理论曲线（虚线）与连续弹性结果（实线）在强耦合区域吻合良好。图 3 显示从 LAMMPS 模拟提取的对齐扭转 WSe₂ 主导分量与 (2⁄3)sech²[2n/(3√η)] 拟合完美，拟合所得 η 随角度呈 1/θ² 依赖，畴壁斜率约为 0.044。该研究进一步将理论应用于对齐扭转 MoTe₂，该材料因承载分数量子反常霍尔态而备受关注。作者从机器学习力场数据中提取声学位移场，发现小角度区域由 u⊥_{n,0} 主导，与理论预测一致。通过大角度区域数据拟合出 c₁、c₂、c₃ 参数后，无需额外拟合参数即能准确预测小角度下的弛豫行为，理论曲线与数据点在小角度区域吻合显著。
该研究的主要结论是：在边际扭转强耦合区域，晶格弛豫由一个依赖于扭转角的单一参数 η 表征，三角和蜂窝孤子网络的位移场均可通过解析表达式描述。该理论为理解小角度扭转双层材料的电子性质提供了系统框架，可广泛应用于扭曲双层石墨烯、对齐与反对齐堆垛的扭曲 WSe₂ 和 MoTe₂，以及其他具有相同堆垛对称性的二维同质双层体系。研究还揭示了一个重要现象：在强耦合极限下存在“涌现普适性”（emergent universality），即大角度下具有不同性质和对称性的模型在小角度下表现出相似的物理行为，其性质仅由孤子网络的对称性和莫尔周期与畴壁宽度之比决定。这一普适性已在实验中得到暗示：扭曲 TMD 材料在 θ ≲ 3° 时强耦合物理即已重要，密度泛函预测与从大角度外推的连续模型在该区域出现显著偏离；同时，在扭曲 WSe₂（3°）和 MoTe₂（2.1°）中观察到的量子自旋霍尔态之间的相似性也可归因于两者均处于强耦合区域，电子性质由涌现的三角孤子网络主导。
该研究的亮点在于首次给出了边际扭转区域晶格弛豫的解析解，将弱耦合微扰理论与强耦合解析理论结合起来，形成了覆盖全部扭转角的统一理论框架。其方法上的创新在于将畴壁斜率匹配条件与傅里叶分量解析形式相结合，以最少的可调参数（仅畴壁斜率或等效 c₁）实现了对复杂位移场的高精度描述。此外，该理论框架还可推广至具有正方或矩形对称性的体系（如 PdSe₂ 和 GeSe），并可直接用于改进基于刚性晶格假设的电子连续模型，为计算完全弛豫的莫尔电子结构提供了可行路径。