关于扭转双层石墨烯应变场可视化的研究报告
作者与发表信息
本研究的核心作者包括Nathanael P. Kazmierczak、Madeline Van Winkle、Colin Ophus、Stephen Carr以及通讯作者D. Kwabena Bediako。研究团队主要来自加州大学伯克利分校化学系、劳伦斯伯克利国家实验室分子铸造厂的美国国家电子显微中心、布朗大学物理系，以及日本国立材料科学研究所（National Institute for Materials Science）的国际材料纳米结构学研究中心。这项研究成果于2021年4月15日在线发表于国际顶级学术期刊 Nature Materials。
学术背景与研究目标
该研究属于凝聚态物理与二维材料科学的前沿交叉领域。通过堆叠晶格常数不同或存在微小旋转错位的二维范德华（van der Waals, vdW）材料，可以形成莫尔超晶格。这种莫尔图案会在材料上叠加一个纳米级的周期性势场，从而剧烈地改变其电子能带结构，为设计新型电子态提供了强大的平台。一个标志性的例子是“魔角”扭转双层石墨烯（magic-angle twisted bilayer graphene, MA-TBG），其在约1.1°的层间转角下会涌现出超导、关联绝缘体等奇异物相。然而，这些莫尔材料的能带结构极其脆弱，微小的结构形变，尤其是由晶格自发弛豫产生的内部层内应变和重构，就能对其产生决定性影响。如何在高分辨率下直接可视化并量化这些被封装在六方氮化硼（hBN）绝缘层内的纳米级形变和应变场，是该领域长期以来面临的巨大挑战。传统显微技术难以穿透封装层，或只能提取一维应变信息。因此，本研究的核心目标是开发一种全新的显微成像技术，旨在穿透封装层，实现对MA-TBG中二维应变张量场的直接、高分辨率可视化，从而揭示其结构弛豫与电子行为之间的内在联系。
详细研究流程与方法
为实现上述目标，研究者开发并实施了一套名为“布拉格干涉术”（Bragg interferometry）的创新方法，其核心是基于四维扫描透射电子显微镜（4D-STEM）。整个研究流程可分解为以下几个紧密衔接的步骤：
样品制备与电子显微数据采集：研究团队采用标准的“撕裂-堆叠”法精心制备了一系列层间扭转角（θm）在0.1°至1.6°之间的hBN/TBG/hBN异质结样品。这些样品被转移至非晶氮化硅薄膜上，用于后续的电子显微分析。在数据采集阶段，研究者使用位于劳伦斯伯克利国家实验室的像差校正扫描透射电子显微镜，在80 kV电压和能量过滤条件下工作。他们将聚焦的电子束在样品上以0.5 nm的步长进行栅格扫描，并在每个扫描点，利用高速直接电子探测器（K3相机）快速记录下一张汇聚束电子衍射图样。这一过程即为4D-STEM，它收集了一个包含数百到数万个衍射图样的庞大四维数据集，每个图样都蕴含了该束斑位置的局部原子堆叠信息。
创新的布拉格干涉术与位移场拟合：这是本研究最核心的方法学创新。当电子束穿透具有微小转角的双层石墨烯时，其衍射图样中会出现两组相互错开、部分重叠的布拉格衍射盘。研究者发现，这两组衍射盘重叠区域的总衍射强度严格遵循一个与层间局部原子位移矢量相关的余弦平方公式。基于此，他们开发了一套拟合程序：首先从每张衍射图样中去背景噪声，然后手动定义12对相互干涉的布拉格盘重叠区域并计算总强度（Ij）。通过非线性回归算法，再将每个扫描像素点的Ij强度信号拟合为一个局部的原子位移矢量。遍历所有扫描点后，即可构建出一幅完整的位移矢量场图，每个像素的位移矢量由一个独特的色相-亮度图例来表示（例如红色系代表AA堆垛，蓝色系代表AB堆垛）。这项技术在概念上类似于全息干涉，它通过解析电子波的干涉信号，巧妙地规避了直接测量布拉格盘位置偏移的传统应变分析方法的困难，首次实现了在hBN封装条件下对TBG结构的直接可视化。
几何分析与应变张量场计算：得到像素级位移场后，研究者进一步进行了复杂的图像后处理和力学分析。首先，他们通过对位移场中AA堆垛区的三角剖分，在单个AB/BA畴的尺度上计算了局部的扭转角（θm）和单轴异质应变（εh）的分布情况。其次，也是最具突破性的环节，是二维应变张量场的计算。由于从衍射强度直接拟合出的位移矢量存在180°的模糊性和跨畴不连续性（即相位包裹），无法直接求导。为此，研究团队自主开发了一套基于几何的“相位解包裹”算法。该算法首先自动识别AA、AB和鞍点（SP）堆垛区域，然后递归地为每个AB畴质心分配一个参考位移矢量，再将每个像素的原始包裹位移解算至最接近其参考矢量的连续值。经过多次移动窗口平滑处理后，残余的不连续性被消除，最终得到一个连续、可微的位移场。对该场进行全广义变分去噪后，通过微分和应变公式，计算出每个像素点的完整二维应变张量，包括正应变、剪切应变、最大剪切应变（γmax）以及重构旋转场（θr）等关键力学量。
理论计算与模型验证：为了将实验观测的力学弛豫与电子性质联系起来，研究团队使用基于密度泛函理论衍生的第一性原理紧束缚模型，进行了电子能带结构计算。他们基于实验应变图构建了一个参数化的原子重构模型，该模型可以“选择性”地打开或关闭由AA区旋转主导或AB区旋转主导的弛豫模式。通过计算不同模式组合下的电子能带和层间隧穿函数，研究者得以量化两种弛豫模式对不同角度下TBG电子结构，特别是平带形成与破坏的孤立影响，并验证了所提出的弛豫模型的有效性。
主要研究结果
通过上述研究方法，该研究获得了一系列突破性的结果。
首先，通过对局域位移矢量的几何分析，研究者在看似均匀的魔角TBG区域（100 nm × 100 nm）内，精确量化了纳米尺度的短程无序度。他们发现，局部扭转角的标准差约为0.03°，而局部异质应变（εh）的标准差介于0.06%至0.09%之间。这种从单个AB畴到相邻AB畴的微小空间涨落，揭示了MA-TBG本征的短程结构无序，这可能是理解不同实验中观测到电子相变出现差异的关键。
其次，应变场映射揭示了两种截然不同的结构弛豫机制，其分界点位于约θm=0.5°。在大于0.5°的区域内，弛豫由AA堆垛区的正向重构旋转（θaa_r > 0）主导，这收缩了高能AA区面积，并导致鞍点区产生垂直于其边界的剪切应变（sxy > syx），伴随鞍点区的负向旋转（θsp_r < 0）。而在小于0.5°的区域内，弛豫则由AB/BA区为趋向低能伯纳尔堆垛而进行的负向反旋转（θab_r < 0）主导，这促成了真正的剪切孤子壁（soliton walls）的形成，此时鞍点区的剪切应变转为平行于孤子壁（syx > sxy），且θsp_r变为正值。在θm≈0.5°时，两种机制达到平衡，鞍点区呈现出纯剪切状态（θsp_r = 0°）。这一发现完美解释了为何在大于魔角（如1.1°）时平带得以稳定存在，而在更小角度时平带会被破坏。
再次，理论计算完美印证了上述双机制模型。通过选择性弛豫模型计算的电子能带结构清晰显示，AA区旋转主要负责定义平坦的低能能带模式，而AB区反旋转则负责隔离这四个平带并使其更弥散。在大于0.5°时，仅施加AA区旋转的层间隧穿函数与完全重构状态的吻合度更高；而在小于0.5°时，仅施加AB区旋转的吻合度更高。这为建设性调控莫尔电子结构提供了直接的力学-电子关联。
最后，研究揭示了异质应变对重构应变场的显著影响。在施加了约0.7%单轴异质应变的魔角TBG区域，研究者观察到原本六重对称的点状鞍点应变场发生了各向异性积累，融合成了显著的条纹状一维应变结构。他们提出的“位移挤压”模型认为，异质应变通过改变莫尔单胞几何形状，迫使SP区角度变小，从而导致连接性的剪切应变场形成，以使材料维持重构。有限元模拟完美重现了这一现象，并揭示了堆垛能作为驱动力相对于层内应变能的重要性。
研究结论与价值
本研究开发的新型布拉格干涉术为莫尔材料的纳米力学研究树立了新的标杆。其科学价值在于：首次穿透封装层，直接且定量地描绘了MA-TBG中的二维应变张量场；揭示了由AA区旋转和AB区反旋转主导的两个截然不同的结构弛豫机制，并阐明了它们对电子平带形成的孤立影响；在单畴尺度量化了本征的短程结构无序，为理解魔角石墨烯中性质多变的电子相提供了统一框架；同时，阐释了宏观异质应变如何通过“位移挤压”转化为对称性破缺的纳米级一维应变条纹。该方法的应用价值在于，布拉格干涉术不仅适用于其他莫尔体系，也可推广到具有共格倒易格矢的非莫尔异质结，并且有望与原位力学拉伸等装置兼容，为探索应变工程调控量子物相开辟了全新的途径。
研究亮点总结
本研究的亮点在于：第一，方法的首创性，布拉格干涉术是一种前所未有的4D-STEM分析策略，它通过解析布拉格盘的干涉强度而非位置来获取位移信息，从而解决了封装样品和低信噪比下的最大难题。第二，发现的深刻性，双阶段弛豫机制的发现及其对电子结构的分解式理解，将结构力学与电子行为在微观机制层面直接挂钩，超越了以往的唯象描述。第三，视角的独特性，在单莫尔畴尺度量化短程无序，以及对异质应变形成条纹相的机理阐释，为理解物理性质的微观不均一性提供了关键的几何-力学图景。