[论文解读] Theory of Twisted Bilayer Graphene Near Commensuration

Hridis K. Pal, Stephen Carter+1 more|arXiv (Cornell University)|Sep 6, 2014

Graphene research and applications被引用 3

一句话总结

本文为近共格扭转角附近的扭转双层石墨烯构建了一个长波长有效理论，通过利用与共格结构的邻近性，将现有框架扩展至非小角度区域。该理论通过层间耦合的高阶傅里叶分量，捕捉了关键的共格效应（如能带隙和拓扑序），使得在共格角度附近的有限角度窗口内可观测到这些现象，而不仅限于精确共格点。

ABSTRACT

Incommensurately twisted graphene bilayers are described by long-wavelength theories, but to date such theories exist only at small angles of interlayer rotation. We construct a long wavelength theory without such a restriction, instead requiring nearness to commensuration. The theory inherits its energy scale from the exactly commensurate bilayer that it is close to. It is a spatial interpolation between the low-energy theories of commensurate structures with the two possible sublattice exchange (SE) symmetries: SE even and SE odd. In addition to generalizing existing theories, our theory brings into experimental reach so far elusive commensuration physics in graphene such as band gaps and nontrivial band topology.

研究动机与目标

• 将现有适用于小角度区域的扭转双层石墨烯长波长理论推广至非小角度区域。
• 使实验可观测到以往因严格扭转角控制要求而难以实现的共格驱动的奇异现象（如能带隙和拓扑序）。
• 通过在SE偶和SE奇共格结构的低能理论之间进行空间插值，统一描述近共格角度下的非共格与共格双层石墨烯。
• 为任何具有莫尔超晶格的双层异质结提供通用框架，无论莫尔图案的起源或底层对称性如何。
• 证明共格效应不仅存在于精确共格系统中，且在共格角度附近有限角度范围内依然存在，显著降低实验约束。

提出的方法

• 通过插值两种具有相反子晶格交换（SE）对称性的共格双层结构（SE偶与SE奇）的低能理论，构建空间变化的有效哈密顿量。
• 从层间耦合的傅里叶分量推导有效哈密顿量，特别强调高阶分量对能带曲率和能隙的贡献，如先前共格研究所示。
• 应用规范变换以对齐两层的狄拉克哈密顿量，将角度依赖性转移到层间耦合项中，该分量以层间势的傅里叶分量表示。
• 在K和K'点附近采用连续近似，层间耦合由高斯势参数化：$ V(\mathbf{r}_i, \mathbf{r}_j) = t_0 e^{-(|\mathbf{r}_i - \mathbf{r}_j|/l_0)^2} $，以支持数值验证。
• 在 $ \theta = 35.57^\circ $（接近共格角 $ \theta_c = 38.21^\circ $）进行数值紧束缚计算，验证密度振荡和类似能隙特征的预测。
• 通过将密度振荡幅度 $ \delta\rho/\rho_0 $ 与 $ \mathcal{V}^2 $ 匹配来验证理论预测，其中 $ \mathcal{V} $ 由共格情况推导得出，确认无自由参数下的一致性。

实验结果

研究问题

• RQ1 能否将扭转双层石墨烯的长波长理论推广至小角度区域之外，以涵盖大角度附近的近共格构型？
• RQ2 层间耦合的高阶傅里叶分量在近共格角度处电子结构中起何种作用？
• RQ3 共格效应（如能带隙和拓扑序）是否不仅可在精确共格角度处观测到，也可在它们附近的有限角度窗口内观测到？
• RQ4 共格角度附近的角窗宽度如何依赖于共格结构的超原胞尺寸和对称性？
• RQ5 在多大程度上可基于最近精确共格结构的性质来预测近共格区域的电子响应？

主要发现

• 该理论成功预测了在 $ \theta = 35.57^\circ $ 处局域态密度的波矢为 $ \delta\mathbf{K} $ 且具有三角对称性的密度振荡，其形式为 $ \delta\rho(\mathbf{r})/\rho_0 \propto \mathcal{V}^2 $，与预测一致。
• 密度振荡的振幅与层间跃迁强度 $ t_0 $ 呈二次方关系，与理论预测一致，并通过数值计算得到验证。
• 增大层间耦合衰减长度 $ l_0 $ 会抑制共格效应，使得 $ \delta\rho/\rho_0 \propto \mathcal{V}^2 $ 指数衰减，与高阶傅里叶分量的衰减一致。
• 在 $ \theta = 35.57^\circ $ 处 $ \delta\rho/\rho_0 $ 的 $ l_0 $ 依赖性与从 $ \theta_c = 38.21^\circ $ 处共格结构提取的 $ \mathcal{V}^2 $ 的 $ l_0 $ 依赖性完全匹配，验证了理论的预测能力。
• 该理论预测在精确共格结构 $ \theta_c = 38.21^\circ $ 处狄拉克点存在 $ 0.09t_0 $ 的能带隙，而在 $ \theta = 35.57^\circ $ 处该能隙缺失，证实了共格在能隙生成中的关键作用。
• 在 $ \theta = 35.57^\circ $ 处数值计算的态密度与基于共格情况推导的 $ \mathcal{V} $ 的理论预测高度一致，无自由参数，其误差在理论精度 $ \delta\theta/\theta_c $ 范围内。