Review of the tight-binding method applicable to the properties of moiré superlattices

本综述为原子级紧束缚方法以及用于模拟各种莫尔超晶格的电子、输运和光学性质的数值技术提供了全面的理论与实践指南，同时也阐明了它们与有效低能连续体模型之间的联系。

要点

想象一下，你拥有两片透明的、具有蜂窝状图案的塑料片（就像石墨烯一样）。如果你将它们完美地堆叠在一起，它们看起来就像单层薄片。但如果你稍微旋转其中一层，或者轻微拉伸它，图案就不再对齐了。相反，它们会创造出一种巨大的、旋转的干涉图样，被称为莫尔超晶格（Moiré superlattice）。

这就像你对着光举起两个窗纱，然后扭转其中一个。你会看到一个巨大的、缓慢移动的波纹图案出现，这个图案比窗纱上的单个小孔要大得多。在原子世界中，这些“巨大的波纹”正是产生最神奇、最奇异物理现象的地方，比如无电阻电流流动的超导现象，或是材料表现出的磁性。

然而，研究这些巨大的原子波对于计算机来说是一场噩梦。因为这种图样非常巨大，以至于这个图样的单个“单元”就包含了数千个原子。试图计算这个庞大群体中每个原子的行为，就像试图通过询问体育场里的每一个人来预测每个人的运动轨迹一样——这既耗时又占用内存。

这篇论文是关于一种特定捷径——紧束缚（Tight-Binding, TB）方法的指南。以下是它如何使用简单的类比来解释这一方法的：

1. 问题所在：太多的原子

论文指出，虽然我们有强大的工具来研究小规模原子群（例如密度泛函理论，即 DFT），但对于这些巨大的莫尔超晶格来说，这些工具运行得太慢了。另一方面，简单的数学模型（连续介质模型）虽然很快，但会丢失微小的细节，比如原子是如何物理性地移动和松弛位置的。

2. 解决方案：紧束缚“邻里”地图

紧束缚方法就像是一张邻里地图。它不是同时计算整个体育场的物理情况，而只是观察一个原子如何与其紧邻的邻居（也就是坐在你旁边的人）进行相互作用。

• 工作原理： 它假设一个原子的行为主要取决于它的邻居是谁以及它们之间的距离有多远。
• 为什么它很棒： 它保留了单个原子的细节（因此它可以观察原子是否被挤压或拉伸），同时又足够快，能够处理数千个原子。它是“金发姑娘原则”下的理想选择：既不过于简单，也不过于缓慢。

3. 工具箱：针对不同材料的不同地图

论文回顾了如何为三种主要类型的材料构建这些“邻里地图”：

• 石墨烯（碳蜂窝结构）： 这里的地图相对简单，侧重于电子如何在碳原子之间“跳跃”。论文展示了通过调整地图中原子间的“距离”，科学家可以精确预测材料何时进入“魔角”超导状态。
• TMDs（过渡金属二硫族化合物）： 这些像是包含金属和其他元素的复杂三明治。这里的地图需要更加详细（使用 11 种不同的“轨道”或电子路径）才能准确描述物理特性。
• hBN（六方氮化硼）： 这通常用作其他材料的平滑基底。论文解释了如何绘制石墨烯的碳原子与该基底的硼/氮原子之间的相互作用图谱。

4. 处理数学问题：“随机游走”技巧

当莫尔图样变得巨大（包含数百万个原子）时，即使是邻里地图也太大了，无法直接求解。论文介绍了一种聪明的技巧，称为线性缩放方法（例如核多项式方法，Kernel Polynomial Method）。

• 类比： 想象你想知道体育场里所有人的平均身高。你不需要测量每一个人。相反，你随机挑选几个人，测量他们，然后使用统计公式来推测整个人群的平均值。
• 结果： 这使得科学家能够在标准计算机上模拟拥有数百万个原子的材料，计算诸如光与材料的相互作用或电流如何流动等内容。

5. “魔法”之源：松弛（Relaxation）

论文的一个关键点是，原子并不是静止不动的雕像；它们会扭动并寻找舒适的位置（松弛）。

• 类比： 想象一群人站在网格中。如果你扭转这个网格，中间的人可能会为了节省空间而挤在一起，而边缘的人则会散开。
• 发现： 紧束缚方法之所以特别，是因为它可以考虑这种“挤在一起”的行为。论文表明，如果你忽略这种松弛，你得到的是错误的物理结果。如果包含它，你就能准确预测“平带”（电子能量水平，在此水平下电子会停滞并开始发生强烈相互作用），从而导致超导等奇异现象。

6. 论文中的实际案例

作者用两个具体的故事展示了这种方法：

• 十二边形晶体： 他们研究了一种形成十二边形图案的扭转石墨烯结构。由于这种图案不会以简单的方式重复，标准的数学方法失效了。紧束缚方法利用“随机游走”技巧，成功预测了光和电在该独特形状下的行为。
• 被捕获的激子： 他们研究了一个 WSe2 层覆盖在扭转石墨烯上的系统。他们展示了石墨烯中原子的“挤压”行为如何创造出微小的陷阱，从而捕捉并固定住“里德堡激子”（一种激发态粒子），这解释了实验中观察到的特定信号。

总结

简而言之，这篇论文是一本关于如何构建和使用特定类型计算机模型的说明书，旨在理解巨大的、扭转的原子图案。它认为紧束缚方法是这项工作的最佳工具，因为它达到了完美的平衡：它既有足够的细节来观察单个原子的移动和松弛，又足够快，能够处理这些巨大的莫后超晶格。它架起了简单、快速的理论与缓慢、超高精度模拟之间的桥梁。

技术摘要

适用于莫尔超晶格性质研究的紧束缚方法综述

问题陈述
由两层二维（2D）材料通过相对旋转或晶格失配形成的莫尔超晶格，已成为探索奇异量子现象（包括非常规超导性、关联绝缘态和拓扑效应）的多功能平台。然而，这些系统的理论建模面临着重大挑战。不同于埃级尺度的材料，由于长波长干涉图案的存在，莫尔系统包含大量的原子（通常为数千到数百万个）。此外，其电子结构高度依赖于晶格弛豫和原子重构。虽然密度泛函理论（DFT）提供了极高的精度，但其计算成本使得模拟大规模莫尔超晶格变得效率低下。相反，连续体模型虽然提供了有效的低能描述，但缺乏捕捉局部无序、应变和特定化学效应所需的原子级分辨率。因此，需要一种能够弥合原子级模拟精度与连续体模型效率之间差距的理论框架。

方法论
本综述重点讨论了紧束缚（TB）方法作为建模莫尔材料的核心工具。本文系统地分类了三类主要材料家族的原子级 TB 哈密顿量构建方式：石墨烯基、过渡金属二硫属化物（TMD）基以及六方氮化硼（hBN）基莫尔超晶格。

1. 哈密顿量构建：

   • 石墨烯： 标准方法采用单粒子 $p_z$ 轨道模型。跳跃参数（$t_{ij}$）使用 Slater-Koster (SK) 关系确定，该关系取决于原子间距和衰减因子。模型通过基于弛豫原子位置修改跳跃项来纳入晶格弛豫。电子相互作用通过平均场近似来处理，包括 Hartree 势（用于长程库仑相互作用）、Hartree-Fock 近似（包含交换项）以及 Hubbard-$U$ 项（用于局部关联）。引入了一种缩放策略，以降低大规模超胞自洽计算的计算成本。
   • TMDs： 由于 TMDs 的复杂性，这些模型通常在每个单层中使用 11 个轨道基组（包括来自金属的 5 个 $d$ 轨道和来自硫族元素的 3 个 $p$ 轨道）。自旋轨道耦合（SOC）通过在位项引入。层间相互作用通过在特定轨道之间添加跳跃项进行建模（例如，同质双层中的 $p$-$p$ 相互作用，异质双层中的 $p$-$d_{z^2}$ 相互作用），其参数通常源自对 DFT 数据进行 Wannier 变换。
   • hBN： 用于扭转双层 hBN 和石墨烯/hBN 异质结构的模型利用了硼（B）和氮（N）原子上的 $p_z$ 轨道。模型考虑了 B 和 N 原子之间的在位能量差异，并采用了各种拟合至 DFT 结果的 SK 参数化方法。

2. 数值技术：
   为了处理大规模哈密顿矩阵（通常 $N > 10^5$ 个原子），综述讨论了超越标准对角化方法的手段：

   • 线性定标方法： 重点介绍了核多项式方法（KPM）和紧束缚传播方法（TBPM）。这些随机方法具有 $O(N)$ 的定标特性，允许计算具有数百万个原子的系统的态密度（DOS）、光学电导率和输运性质。
   • 机器学习： 综述了利用深度学习（如 DeepH、HamGNN、DeepTB）从训练数据构建具有 ab initio 质量的 TB 哈密顿量的最新进展，从而加速了参数化过程。

3. 与连续体模型的联系：
   论文详细阐述了如何从原子级 TB 哈密顿量导出低能有效连续体模型。它解释了莫尔势是如何通过层间隧穿的傅里叶分量产生的。综述强调，虽然局部（动量无关）势能够捕捉到平带的出现，但必须引入非局部（动量相关）隧穿项，才能重现在弛豫后的 TB 模型中所观察到的粒子-空穴不对称性及其他特征。

核心贡献与结果

• 全面概述： 本文为用于石墨烯、hBN 和 TMD 莫尔系统的特定 TB 哈密顿量、SK 参数和弛豫势（例如 LAMMPS 势）提供了统一的参考。
• TB 模型验证： 综述证明了经过适当参数化并包含弛豫效应的 TB 模型可以重现关键的实验和 DFT 结果。例如，扭转双层石墨烯（TBG）的 TB 计算成功预测了费米速度消失且平带出现的“魔角”（$\theta \approx 1.05^\circ - 1.1^\circ$）。
• 晶格弛豫的作用： 作者强调，原子弛豫并非微小的修正，而是一个关键因素。在 TBG 中，弛豫导致 AA 和 AB 畴的形成，显著改变了能带结构并诱导了粒子-空穴不对称性。在 TMDs 和 hBN 中，弛豫同样决定了能带的平坦度和带隙的大小。
• 处理相互作用： 综述展示了如何利用平均场 TB 方法（Hartree、Hartree-Fock、Hubbard）来研究平带机制下的关联相，如 Mott 绝缘态和铁磁性。
• 案例研究： 文中提供了两个具体实例来展示 TB 方法的效用：
  1. 十二面体石墨烯准晶体： 通过在拥有 1000 万个原子的系统上使用 TBPM，研究揭示了归因于层间相互作用的 DOS 和光学电导率中的明显峰值，并预测了磁场下的新朗道能级。
  2. WSe$_2$/TBG 中的里德堡莫尔激子： 通过将 TB 与分子动力学相结合，作者解释了 WSe$_2$ 中里德堡激子的空间局域化现象，这是由底层弛豫后的 TBG 中的 AA 区电荷积累引起的，并与实验反射光谱和光致发光光谱相吻合。

意义与范围
本文将紧束缚方法定位为研究莫尔超晶格的核心且实用的工具。其重要性在于其独特的平衡：既保留了模拟晶格弛豫、应变和化学特异性所需的原子级分辨率，又保持了足够的计算效率，足以模拟包含数千到数百万个原子的现实超胞。

作者指出，本综述旨在为研究人员应用 TB 方法研究莫尔系统提供“理论与实践指南”。它弥合了 ab initio DFT 计算与有效连续体模型之间的鸿沟，提供了一种纳入多体相互作用和外部场的系统方法。文章最后指出，尽管 TB 非常强大，但未来仍需在扩展至非六方晶格、改进非正交基组的线性定标算法，以及利用机器学习进行新莫尔材料的自动化发现等方面开展工作。本综述并不声称解决了所有关联问题，而是断言 TB 为理解平带相关现象的起源提供了必要的精确起点。