关于《Physical Review Research》期刊2021年论文《测量转角双层石墨烯局部莫尔晶格异质性》的学术研究报告

一、 研究团队与发表信息

本项研究的主要作者包括来自荷兰莱顿大学莱顿物理研究所Huygens-Kamerlingh Onnes实验室的Tjerk Benschop、Tobias A. de Jong、Vincent Stalman、Sense Jan van der Molen和Milan P. Allan，以及来自西班牙巴塞罗那科学技术学院光子科学研究所的Petr Stepanov、Xiaobo Lu和Dmitri K. Efetov。其中，Milan P. Allan为通讯作者。该研究以《Measuring local moiré lattice heterogeneity of twisted bilayer graphene》为题，于2021年2月16日发表在《Physical Review Research》第3卷第013153页。

二、 学术背景与研究目的

本研究隶属于凝聚态物理和纳米材料科学领域，具体聚焦于魔角石墨烯及相关莫尔超晶格材料的结构表征。近年来，转角双层石墨烯因其在特定“魔角”（约1.1°）下展现出的超导、关联绝缘态等一系列新奇物理现象而成为研究热点。这些电子性质对莫尔超晶格的周期和均匀性极为敏感，微小的扭转角变化或应变都会显著改变其电子结构。然而，在样品制备过程中，由于“撕裂-堆叠”技术引入的应力、污染以及不可避免的异质性，不同器件之间甚至同一器件内部都可能存在扭转角的不均匀性和应变，这被认为是导致器件性能差异、限制其质量的关键因素。尤其是对于没有顶部六方氮化硼封装的“开放式”器件，尽管理论上可能具有更干净的界面，但从未观测到超导现象，其结构均匀性一直缺乏系统的纳米尺度表征。

因此，本研究旨在开发一种新型的、高空间分辨率的表征方法，以连续映射的方式可视化并定量分析转角双层石墨烯器件中莫尔晶格的异质性，包括局部有效扭转角和局部各向异性（与异质应变相关）。研究目标包括：1）量化开放式转角双层石墨烯器件中扭转角的变化范围；2）区分扭转角效应与应变效应；3）将开放式器件的均匀性与已报道的封装器件进行比较，探讨开放式器件缺乏超导性的可能结构原因。

三、 详细研究流程与方法

本研究主要包含三个核心流程：样品制备、扫描隧道显微镜成像以及新开发的空间锁相算法数据分析。

流程一：样品制备与表征 研究团队采用改进的“撕裂-堆叠”技术制备了目标扭转角在1.5°至2.0°之间的开放式转角双层石墨烯器件。具体步骤包括：将同一石墨烯薄片预切割成两半以确保初始晶体学对准；利用六方氮化硼和聚二甲基硅氧烷/聚碳酸酯印章依次拾取和旋转其中一半石墨烯，以实现目标转角；最终将转角双层石墨烯/六方氮化硼异质结构转移到带有预图案化电极和石墨烯栅极的SiO2/Si基底上。制备完成后，器件在插入扫描隧道显微镜前，于超高真空环境中进行350°C、12小时的退火处理以清洁表面。所有扫描隧道显微镜测量均在4.2 K的低温下进行。

流程二：扫描隧道显微镜形貌成像 研究使用扫描隧道显微镜获取了大面积、高质量的拓扑形貌图。成像的关键在于能够同时分辨顶层石墨烯的原子晶格和莫尔超晶格。如图2(a)所示，在小尺度图像中可清晰看到原子和莫尔周期。为了进行异质性研究，需要包含大量莫尔晶胞的大视野图像，如图2©和图4(a)所示，其尺度达数百纳米。这些图像是后续算法分析的原始数据输入。

流程三：空间锁相算法开发与应用（核心创新方法） 本研究最大的创新在于开发并应用了一种名为“空间锁相算法”的新型图像分析方法，用于从扫描隧道显微镜形貌图中提取纳米级的局部晶格变形信息。该算法的工作流程如下：

1. 参考波矢确定：首先，对获取的大尺度拓扑图像进行快速傅里叶变换，识别出莫尔超晶格的布拉格峰（如图4(b)中的q1, q2, q3）。通过拟合高斯峰，确定这些倒空间波矢的平均位置。为了分析各向异性，算法强制设定三个参考波矢具有相等的模长和彼此60°的夹角，从而定义一个理想的、各向同性的三角参考晶格。

2. 局部位移场提取：算法的核心思想是将测量的形貌图视为一个理想参考晶格经过局部变形后的结果。这个变形由一个二维位移场u®描述。为了提取u®，算法对每个参考波矢q_j执行空间锁相操作：

   • 混频：将测量的形貌图tm®与一个形式为exp(i q_j · r)的复数参考平面波相乘。
   • 低通滤波：对乘积结果进行高斯低通滤波，滤除高频分量（主要是原子晶格信号和倍频信号）。滤波窗口的大小（如图4(a)中蓝色圆圈所示）是关键参数，需要在空间分辨率（窗口小）和信噪比（窗口大）之间取得平衡，通常选取几个莫尔周期的尺度。
   • 相位提取：滤波后得到的复信号的相位图φ_j®即包含了位移场在q_j方向上的投影信息，满足φ_j® = q_j · u®。对相位图进行相位解包裹处理，以消除2π跳变。

3. 变形参数分解：通过至少两个方向的相位图（例如φ1和φ2），可以求解出完整的二维位移场u®。随后，计算位移场的雅可比矩阵J = I + ∇u，其中∇u是位移梯度张量。对雅可比矩阵进行极分解：J = W A = W V D V^T。

   • 旋转矩阵W：描述了整个晶格的刚性旋转（相对于基底）。
   • 缩放与各向异性矩阵A：进一步分解为旋转矩阵V（指示各向异性的主轴方向ψ®）和对角缩放矩阵D = diag(d1, d2)，其中d1 ≥ d2。
   • 有效扭转角θ*®：通过缩放因子的几何平均值√(d1 d2)计算局部莫尔波长λ®，再利用关系式θ® = 2 arcsin(a / (2λ®))（其中a为石墨烯晶格常数）转换为局部有效扭转角。
   • 局部各向异性κ®：定义为κ = d1/d2。κ=1表示各向同性，κ>1表示存在各向异性，其方向由ψ®给出。
   • 异质应变ε®：基于Kerelsky等人（2019）的模型，假设异质性主要来源于其中一层石墨烯的单轴应变，而另一层仅发生旋转。通过分析相对位移场，可以将测得的各向异性κ®和方向ψ®转化为单层所承受的应变大小ε®。

4. 算法验证与误差评估：为验证算法的准确性，研究团队进行了自洽性检查。他们使用Lawler-Fujita算法（另一种基于位移场校正晶格畸变的方法）对原始图像进行“去畸变”校正，生成一个理论上完美的周期晶格图像。然后，将空间锁相算法应用于这个校正后的图像。结果显示，提取出的残余扭转角变化和各向异性比原始结果小一个数量级以上，证明了算法的自洽性和高精度。

四、 主要研究结果

应用上述空间锁相算法，研究团队对开放式转角双层石墨烯器件进行了详细分析，得到了以下关键结果：

1. 局部有效扭转角分布：图4©展示了从图4(a)形貌图提取的局部有效扭转角θ*®的空间分布图。结果显示，扭转角在空间上平滑变化，平均值为2.02°。在排除边界效应后，整个视野内的扭转角标准偏差为0.033°，峰峰值变化为0.23°。在数百纳米尺度的区域内（如图4©中红色方框所示），标准偏差可低至0.02°，峰峰值变化仅为0.08°。这表明，即使在开放式器件中，扭转角的均匀性也相当高。

2. 局部各向异性与异质应变：图4(d)和(e)分别展示了局部各向异性κ®及其方向ψ®的分布。除了一个来源不明的垂直条纹特征外，各向异性分布也较为平滑。基于应变模型计算得到的异质应变ε®如图4(f)所示，平均应变约为0.14%，标准偏差为0.09%。

3. 与封装器件的比较：为了与已报道的封装器件（如使用扫描超导量子干涉器件显微镜的研究）进行公平比较，研究团队考虑了不同技术的点扩散函数差异。他们将提取的扭转角分布图用30纳米宽的高斯滤波器进行模糊化处理，以模拟扫描超导量子干涉器件显微镜的分辨率。处理后，全视野的扭转角峰峰值变化为0.20°，标准偏差为0.036°。这一结果与封装器件中测得的均匀性水平（变化小于4%）相当。这表明，就扭转角均匀性而言，高质量的开放式器件可以与封装器件相媲美。

4. 应变与扭转角异质性的影响对比：研究引用连续模型的计算结果指出，对于接近魔角的样品，约0.1%的异质应变会导致范霍夫奇点分裂约5 meV；而约0.03°的扭转角变化也会产生类似的能带结构扰动。本研究中观测到的应变不均匀性（~0.09% std）和扭转角不均匀性（~0.02° - 0.033° std）对电子结构的影响程度大致相当。因此，在制备高质量样品时，必须同时考虑和控制应变与扭转角两种不均匀性。

5. 方法稳健性验证：为了排除扫描隧道显微镜实验中压电陶瓷漂移对测量结果的干扰，研究团队对同一区域以不同扫描速度和不同时间点采集的多幅形貌图进行了重复分析。结果显示，提取出的主要变形特征（如垂直条纹、各向异性最小值位置）在不同图像中高度一致，而仅有一些与扫描线相关的假象不同。这证实了所观测到的异质性是样品固有的，而非仪器假象。

五、 研究结论与意义

本研究得出结论：通过新开发的空间锁相算法，发现在所研究的高质量开放式转角双层石墨烯器件中，局部扭转角在数百纳米尺度上的变化非常小（峰峰值约0.08°），其均匀性水平与六方氮化硼封装的器件相当。同时，器件中存在可观的局部异质应变（平均0.14%），其对电子结构的影响与观测到的扭转角不均匀性影响程度相似。

这项研究的科学价值在于： 1. 方法学创新：开发了一种通用的、高空间分辨率的空间锁相算法，能够以亚莫尔晶胞的分辨率连续映射莫尔晶格的局部扭转角、各向异性和应变。该方法不限于扫描隧道显微镜，也可应用于原子力显微镜、低能电子显微镜等其他成像技术，以及各种莫尔材料体系。 2. 对转角石墨烯物性的新认识：研究表明，开放式器件在结构均匀性上并不逊色于封装器件。这间接提示，开放式器件中从未观测到超导现象的原因，可能并非源于其莫尔晶格本身更大的不均匀性，而更可能与缺乏顶部六方氮化硼层所导致的对称性破缺、环境屏蔽或界面效应有关。这为理解魔角石墨烯超导的起源和器件工程提供了新的视角。 3. 指导样品制备：定量揭示了应变和扭转角不均匀性对电子结构具有可比性的影响，强调在追求精确魔角的同时，控制和表征应变同样至关重要。

六、 研究亮点

1. 新颖的分析方法：空间锁相算法是本研究最突出的亮点。它实现了对莫尔晶格畸变的高分辨率、定量化、可视化分析，能够分离旋转、缩放和各向异性等不同变形模式，超越了传统傅里叶变换（空间分辨率低）或单点拟合方法（统计性差、对噪声敏感）的局限。
2. 重要的实验发现：首次在纳米尺度上系统量化了开放式转角双层石墨烯器件的扭转角与应变不均匀性，并得出了其均匀性可与封装器件媲美的关键结论，挑战了关于开放式器件质量差的普遍猜想。
3. 明确的物理见解：通过对比理论计算，明确了观测到的应变和扭转角不均匀性对电子结构具有相近量级的影响，为材料制备和物性调控提供了清晰的指导。
4. 技术的普适性与启发性：论文明确指出该算法可广泛应用于其他莫尔材料和探测技术，并建议其潜力可扩展至其他领域，例如通过准粒子干涉测量来研究能带的拓扑性质。

七、 其他有价值内容

论文在附录中提供了大量补充信息，极具参考价值： * 算法数学细节：附录A详细推导了算法的数学框架，包括晶格变形的描述、位移场与相位的关系、雅可比矩阵的分解公式等。 * 莫尔位移与相对位移的关系：附录B建立了测量的莫尔晶格位移场与上下两层石墨烯相对位移场之间的理论联系，并给出了从各向异性推导单轴应变的公式。 * 相位解包裹处理：附录C讨论了相位奇异点的处理策略。 * 器件结构与数据概览：附录D展示了器件示意图，附录G和K汇总了多个器件的数据结果，证明了研究结论的可靠性。 * 算法精度验证：附录E通过Lawler-Fujita校正后的残差分析，定量证明了算法的高精度。 * 压电漂移排除：附录F通过多组数据对比，有力排除了仪器漂移是主要信号来源的可能性。 * 与扫描超导量子干涉器件显微镜的定量比较：附录I详细描述了如何进行点扩散函数修正以实现与扫描超导量子干涉器件显微镜数据的公平对比。 * 与传统应变模型误差对比：附录J通过高斯拟合莫尔斑点的误差分析，暗示了传统方法在低像素密度（本研究约5像素/莫尔周期）下面临的挑战，凸显了新算法的优势。

这项研究不仅为转角二维材料领域提供了一种强大的新型表征工具，还通过精密的实验测量和数据分析，深化了对魔角石墨烯器件结构-物性关系的理解，特别是为探索开放式器件中超导缺失的原因指明了新的方向。