Twisted bilayer graphene from first-principles: structural and electronic properties

本文利用密度泛函理论，对宽范围扭转角下的扭转双层石墨烯进行了全面的从头算研究，提供了完全弛豫的原子结构和详细的电子性质，为未来的多体研究奠定了基础的从头算参考。

要点

想象两张石墨烯（一种由单层碳原子按蜂窝状排列构成的材料）相互堆叠。现在，想象将其中一张相对于另一张轻微扭转。这就形成了“扭转双层石墨烯”（tBLG）。

当你以恰到好处的角度（一个特定的“魔角”）扭转它们时，神奇的事情发生了：内部的电子不再四处飞驰，而是被固定在原地，形成一片平坦、宁静的能量之海。这种状态允许出现诸如超导性（零电阻导电）等奇异行为。

本文就像是一次高分辨率的微观制图探险。作者旨在利用名为“第一性原理”的强大计算机模拟，精确理解这种扭转结构的外观及其内部电子的行为。

以下是他们旅程与发现的分解说明，辅以简单的类比：

1. 挑战：“像素”问题

通常，模拟这些扭转的片层就像试图用一个仅擅长处理小方块的计算机程序，去绘制一幅巨大而复杂的挂毯。“扭转”会产生一个巨大的重复图案（称为莫尔条纹），且随着角度变小，该图案会变得极其巨大。标准的计算机方法（如“平面波”密度泛函理论 DFT）就像是用一把粗刷子去画壁画；它们虽然准确，但过于笨重缓慢，难以处理大片扭转片层上的微小细节。

解决方案： 作者使用了一种特殊的、经过优化的“局域基组”方法（使用 SIESTA 代码）。这就像使用一把笔尖精细、灵活的画笔，可以聚焦于特定的原子，而无需一次性描绘整个宇宙。这使得他们能够模拟包含数万个原子的片层，达到非常小的扭转角度（低至约 1 度），而这些角度此前因难以精确建模而未被触及。

2. 校验地图：“两支画笔是否一致？”

在信任他们那支新的精细画笔之前，他们将其与旧的、笨重的画笔（使用 VASP 代码）在中等扭转角度（2.45 度）下进行了对比。

• 结果： 两种方法几乎完美吻合。原子处于相同的位置，作用于它们的力也完全相同。这证明他们的新方法足够准确，可以信赖用于处理更大、更困难的任务。

3. 扭转的形状：“起皱的毯子”

当你扭转两张片层时，它们不会保持完全平坦。它们会起皱并移动，以找到最舒适的位置，就像毯子落在床上一样。

• 发现： 作者精确计算了原子的移动方式。他们发现，原子主要在特定的位置（称为"AA 位点”，即蜂窝图案完美对齐的点）周围发生位移。
• 类比： 他们将详细的原子地图与“连续弹性模型”进行了比较，后者就像是一张平滑的数学橡胶片近似模型。他们发现，即使在他们模拟的最小角度下，详细的原子地图也与平滑的橡胶片模型完美匹配。这意味着科学家可以使用更简单的橡胶片模型来预测原子的排列方式，从而节省时间。

4. 电子的速度：“交通堵塞”

在这些扭转片层中，电子通常具有“费米速度”（即它们的移动速度）。在“魔角”处，这种速度应降至接近零，从而形成电子被卡住的平坦能带。

• 发现： 作者将他们的结果与一个高度精确的数学模型（“精确 k·p 模型”）进行了比较。他们发现趋势是一致的：随着角度接近魔角，电子的速度减慢。
• 转折： 然而，存在一个微小的“偏移”。他们模拟中的电子在略微不同于数学模型预测的角度上减速。这就像两名赛跑者瞄准同一个终点线，但从略微不同的起跑线出发。作者认为，这种差异源于他们如何处理层间的“胶水”（范德华力）以及描述电子相互作用的特定数学方法。

5. 电子的“纹理”：“波图案”

他们最酷的工作之一是观察电子的“波函数”。想象电子不是一个小球，而是池塘中的一圈涟漪。

• 发现： 他们在三维空间中绘制了这些涟漪的分布。他们发现，涟漪的形状会随着扭转角度的变化而改变。
  • 在较大角度下，涟漪看起来像是紧抱着不同区域之间的“墙壁”。
  • 随着角度变小（接近魔角），涟漪转变为紧抱着图案对齐的“中心”。
• 手性检查： 他们还检查了材料中两个不同点处这些涟漪的“手性”（chirality）。在普通石墨烯中，这些点具有相反的手性（就像左手和右手）。而在扭转双层石墨烯中，他们发现这两个点具有相同的手性。这是该材料独特的指纹，解释了其为何具有如此特殊的拓扑性质。

总结

简而言之，本文构建了一个高度详细、原子级精度的扭转石墨烯三维模型。他们证明了这种新的、高效的计算机方法与那些笨重、缓慢的方法同样有效。他们确认了原子以可预测的方式起皱，且这种起皱与简单的橡胶片数学模型相符；他们还精确绘制了电子如何随扭转角度变化而减速并改变其“形状”。这为未来希望研究更复杂效应（例如这些材料如何实现无电阻导电）的科学家提供了坚实可靠的基石。

技术摘要

技术摘要：基于第一性原理的扭曲双层石墨烯：结构与电子性质

问题陈述
扭曲双层石墨烯（tBLG）展现出奇异现象，如超导性、关联绝缘体行为和拓扑态，这些现象由大尺度原子重构与孤立平带之间的相互作用驱动。尽管连续模型（$k \cdot p$）和紧束缚方法提供了有价值的见解，但它们往往缺乏与高分辨率实验进行比较所需的原子尺度结构和电子特征。相反，第一性原理密度泛函理论（DFT）研究受限于计算成本，通常将模拟限制在较大的扭转角，或依赖难以处理小角度所需的大尺寸莫尔超胞的平面波代码。此外，先前的 DFT 研究往往聚焦于可用简单模型描述的性质，而忽略了近期实验测量的详细波函数纹理。

方法论
作者对宽范围扭转角的 tBLG 进行了全面的第一性原理研究，角度低至 $\theta = 0.987^\circ$。该研究采用 SIESTA 代码进行 DFT 计算，使用优化的局域基组（数值原子轨道，NAOs），具体利用了针对石墨烯优化的单-$\zeta$极化（SZP）和双-$\zeta$极化+f（DZPF）基组。这种方法使得模拟包含数万个原子的大尺寸莫尔超胞成为可能，而这对于标准平面波方法而言在计算上是不可行的。

关键方法论步骤包括：

• 结构弛豫：通过将力最小化至 10 meV/Å以下，获得了完全弛豫的公度结构。研究采用了 PBE 交换关联泛函，并结合 DFT-D3 色散校正来处理层间范德华相互作用。
• 验证：将结果与使用 VASP 进行的平面波 DFT 计算（针对 $2.45^\circ$ 扭转角）以及针对晶格弛豫的连续弹性模型进行了基准对比。
• 电子结构：沿 $K-\Gamma-M-K$ 路径计算了能带结构。研究将这些结果与重现从头算紧束缚结果的“精确”$k \cdot p$ 连续模型进行了比较。
• 波函数分析：提取并分析了实空间波函数，以确定其莫尔尺度特征（例如 AA、AB/BA、DW、AAz）和对称性性质，包括狄拉克节点处的手性。

关键结果

1. 结构一致性：对于 $\theta = 2.45^\circ$，局域基组 DFT（SIESTA）获得的原子弛豫与平面波 DFT（VASP）显示出极好的一致性，坐标差异小于弛豫幅度的 1%。此外，晶格弛豫模式与连续弹性模型完全吻合，直至最小的模拟角度 $0.987^\circ$。
2. 电子性质与角度偏移：虽然 DFT 得出的电子性质（费米速度、能带宽度）与精确的 $k \cdot p$ 模型在定性上吻合，但观察到了系统性的扭转角偏移。DFT 能带并未在 $k \cdot p$ 模型中看到的“魔角”处表现出费米速度和能带宽度的极小值；相反，这些特征在 DFT 结果中似乎出现在略小的角度处。为了在魔角附近使能带结构与 DFT 结果对齐，需要在 $k \cdot p$ 计算中引入 $\Delta\theta \approx 0.05^\circ$ 的偏移。
3. 波函数特征：该研究详细描述了波函数纹理随扭转角的演变。它证实，随着角度减小，平带和色散带在特征上经历了连续演变（例如，从类畴壁状到 AA 局域化）。至关重要的是，对称性性质的分析证实，$K$ 和 $K'$ 谷处的狄拉克锥具有相同的手性，这是 tBLG 中脆弱拓扑的一个特征。
4. 基组与泛函敏感性：作者调查了其 DFT 结果与 $k \cdot p$ 模型（基于 SCAN 泛函和 rVV10 范德华校正参数化）之间差异的来源。他们发现，交换关联泛函和范德华校正的选择显著影响层间距，进而影响耦合强度和能带平坦度。SIESTA 中的 SZP 基组产生的层间距略大于 $k \cdot p$ 模型中使用的参数，这表明“魔角”的定义对这些计算细节是敏感的。

意义
本文确立了 tBLG 微观结构和电子性质的稳健 ab initio 参考点。通过证明优化的局域轨道 DFT 能够以原子分辨率准确模拟大尺寸莫尔超胞，该工作弥合了连续模型与高分辨率实验数据之间的差距。作者断言，他们的结果为未来纳入多体效应的研究提供了必要的基础，提供了一条计算上可行的途径，以在以往如此大规模系统中无法达到的精度下模拟具有数万个原子的范德华材料。对扭转角偏移的识别以及对波函数对称性的详细表征，为完善 tBLG 理论模型提供了关键见解。