Interplay of interlayer distance and in-plane lattice relaxations in encapsulated twisted bilayers

本文提出了一个理论模型，证明了封装界面的刚度显著影响扭转双层材料中的晶格弛豫，特别是提高了从弱弛豫向强弛豫转变的临界扭转角，并实现了与实验数据更好的吻合。

要点

想象一下，你拥有两张具有图案的、非常粘的壁纸。如果你将一张纸直接叠在另一张之上，但稍微旋转一下，它们的图案就不会完美对齐。相反，它们会创造出一个巨大的、重复的“影子”图案，称为莫列波纹（moiré pattern）。

在量子材料的世界中，科学家们通过旋转这些原子层来创造新的电子特性。然而，原子是“懒惰”的；它们想要安顿在最舒适、最节省能量的位置。因此，当我们旋转这些层时，原子并不仅仅是留在原地——它们会四处挪动、拉伸和挤压，以寻找最佳的契合点。这种挪动被称为晶格弛豫（lattice relaxation）。

问题所在：“悬浮”与“三明治”

长期以来，科学家研究这些扭转层的方式，都是假设它们是悬浮在真空中的（悬浮态）。他们知道，在某些很小的角度下，原子会发生剧烈的挪动（强弛豫），从而创造出由应力墙分隔开的完美对齐岛屿。而在较大的角度下，它们的挪动则较少（弱弛豫）。

但在真实的实验中，这些层并不是悬浮着的。它们通常被夹在其他保护层（如六方氮化硼）之间以保持稳定。这被称为封装（encapsulation）。

这篇论文提出了一个问题：这个“三明治”结构是否改变了原子的挪动方式？

发现：“硬质三明治”效应

作者 V. V. Enaldiev 构建了一个数学模型来回答这个问题。他意识到，这个三明治的保护性“面包”（封装层）起到了一个坚硬约束的作用。

这里有一个类比：

• 扭转层： 想象两块具有蜂窝图案的柔软、有弹性的橡胶垫。当你旋转它们时，蜂窝结构试图通过对齐来达到完美的匹配。
• 封装层： 现在，想象你将这些橡胶垫夹在两块非常坚硬、刚性的木板之间。
• 结果： 在中间（即两块橡胶垫接触的地方），橡胶想要上下挤压以寻找完美的契合点。但上下两块硬板却在说：“不，保持平整！”这些木板阻碍了橡胶上下移动。

论文发现，由于“木板”（封装层）是坚硬的，它们抑制了原子的垂直运动。原子无法像它们希望的那样进行挤压。

主要发现：改变“临界点”

由于原子无法那么容易地挤压，因此需要一个更小的扭转角度才能迫使它们开始水平方向上的挪动，以寻找舒适区。

可以把这想象成一个跷跷板：

1. 悬浮态（自由漂浮）： 原子可以自由地上下移动。它们只有在扭转角度非常小时（大约 1° 到 2.5°）才会开始水平挪动。
2. 封装态（三明治结构）： 原子在垂直方向上被固定住了。因为它们不能利用“上下”移动的小技巧来节省能量，所以它们被迫在更早的时候（在更大的扭转角度下）就开始进行水平挪动。

论文计算出，对于一个完全刚性的三明治结构，这个“临界点”（即原子开始显著挪动的角度）从大约 3.8° 移动到了 4.5°。

这为什么重要

作者展示了，通过调整他模型中的一个数字（代表三明治的硬度），他的预测就能与现实世界的实验完美吻合。

• 现实世界的证明： 实验表明，夹在三明治结构中的扭转层，其行为与悬浮的层不同。
• 模型的成功： 该模型解释了原因：三明治使层在垂直方向上变得更“硬”，这改变了原子决定重新排列自身的角度。

简而言之

这篇论文解释了，当你用一个保护壳包裹住扭转的原子层时，这个壳就像一个坚硬的夹具。这个夹具阻止了原子上下移动，从而迫使它们在不同的角度下进行侧向位置的重新排列，而如果它们是自由悬浮的话，情况则会有所不同。这种简单的“硬度”变化，解释了为什么现实中的实验看起来与那些忽略了保护壳的旧理论不同。

技术摘要

技术摘要：封装扭转双层结构中层间距与面内晶格弛豫的相互作用

问题陈述
扭转范德华异质结，特别是扭转过渡金属二硫化物（TMD）双层结构，展现出介观尺度的莫尔超晶格，这显著改变了其电学和光电性质。超晶格的形成驱动原子弛豫以最小化总能量，将刚性的莫尔图案转化为由应变吸收畴壁分隔的共格畴（commensurate domains）阵列。虽然悬浮双层的从刚性莫尔图案到弛豫畴结构的转变已得到表征（平行排列的交叉角 $\theta^* \approx 2.5^\circ$，反向排列的 TMD 约为 $\approx 1^\circ$），但近期的实验表明，将 TMD 双层封装在六方氮化硼（hBN）中会增强晶格弛豫，并将这些交叉角推向更高的数值。现有的晶格弛豫理论模型并未考虑封装层所施加的机械约束，特别是界面处层间距的调制作用。

方法论
作者开发了一个显式包含封装效应的扭转 TMD 双层晶格弛豫理论模型。该方法包括：

1. 粘附能公式化：利用依赖于面内偏移 $r_0$ 和层间距 $d$ 的堆叠相关粘附能密度 $W_{ad}(r_0, d)$。
2. 封装建模：将扭转双层视为嵌入在多层平板之间的结构（封装）。扭转双层与封装平板之间的相互作用通过一个有效参数 $k$ 表征的谐振势进行建模，该参数代表了双层/平板界面的刚度。
3. 层间距调制：通过平衡来自扭转层的堆叠相关压力 ($P_{in}$) 与来自封装平板的回复压力 ($P_{out}$)，推导出新的层间距调制表达式 $d_{enc}(r_0)$。这导致了与悬浮双层相比被抑制的层间距变化。
4. 面内弛豫计算：通过最小化包含面内弹性能和局部粘附能（使用修正后的 $d_{enc}$）的总能量泛函，来确定面内位移场。计算采用基于莫尔超晶格倒格矢的傅里叶级数展开。
5. 量化分析：通过定义 $R_\nu(\theta)$（即 XX 堆叠区域相对于莫尔超原胞的面积比），在各种扭转角和界面刚度参数下分析晶格弛豫的强度。

主要贡献与结果

• 模型开发：论文引入了一个单一的唯象参数 $k$，用于表征扭转双层与封装材料之间界面的刚度。该参数控制着层间距调制的抑制程度。
• 封装效应：模型证明，封装抑制了可达到的层间距范围。这种抑制增加了能量上有利（如菱方 3R 堆叠）与不利（XX 堆叠）堆叠之间的能量差。
• 交叉角偏移：主要结果是，增加界面刚度 ($k$) 会提高弱弛豫与强弛豫机制之间的交叉扭转角 ($\theta^*$)。对于悬浮双层 ($k=0$)，$\theta^*$ 计算为 $\approx 3.8^\circ$。对于完全刚性界面 ($k \gg 4\varepsilon$)，$\theta^*$ 增加至 $\approx 4.5^\circ$。
• 定量一致性：通过调整刚度参数以模拟完全刚性界面（可能是由于 TMD 与 hBN 之间的晶格不共格性），该模型预测了交叉角的变化 $\Delta\theta^* \approx 0.75^\circ$。这与实验观察到的平行排列 WSe$_2$ 双层中 $\Delta\theta^* \approx 1^\circ$ 的变化高度吻合。
• 机制行为：研究证实，在强弛豫机制下（小扭转角），XX 堆叠区域的大小与扭转角无关，遵循 $\theta^2$ 的弛豫强度参数标度律，这种行为在不同的封装强度下保持普适性。

意义
该论文声称其意义在于提供了一个理论框架，能够弥合悬浮双层模型与封装系统实验结果之间的差距。通过识别封装界面的刚度会抑制层间距调制，作者解释了为什么封装系统表现出增强的晶格弛豫和更高的交叉扭转角。该模型通过将界面视为有效刚性，成功实现了与现有实验数据的定量一致，这表明 TMD/hBN 界面处的晶格不共格性是决定此类异质结机械响应的关键因素。作者指出，该框架可以扩展到其他扭转电子学异质结构（如 WX$_2$/MoX$_2$ 系统），以预测应变幅度和畴壁网络行为的变化。