本研究的主要作者为Hyobin Yoo、Rebecca Engelke、Stephen Carr等，通讯作者为Philip Kim（邮箱：pkim@physics.harvard.edu）。参与机构包括哈佛大学物理系、明尼苏达大学航空航天工程与力学系、堪萨斯大学数学系、密歇根大学材料科学与工程系、滑铁卢大学量子计算研究所及化学系、日本国立材料科学研究所、首尔国立大学物理与天文系及材料科学与工程系、明尼苏达大学数学系以及哈佛大学约翰·保尔森工程与应用科学学院。这项研究于2019年5月发表在《自然-材料》（*Nature Materials*）期刊上，卷18，页码448-453。

一、 研究背景 本研究属于凝聚态物理和材料科学领域，聚焦于二维范德华（van der Waals, vdW）异质结的界面物理。二维材料层间通过较弱的范德华力结合，允许任意堆叠和旋转，形成摩尔（moiré）超晶格。扭转双层石墨烯（twisted bilayer graphene, TBG）是这一领域的明星体系，通过控制层间扭转角（θ），可以连续调节摩尔超晶格的周期，从而改变其电子能带结构。特别是当扭转角接近理论预言的“魔角”（magic angle，约1.1°）时，电子能带变得极平，态密度发散，导致电子关联效应显著增强，已观测到莫特绝缘体和非常规超导等关联电子态。

然而，之前的魔角物理研究大多基于一个关键假设：两层石墨烯是刚性扭转的，即原子结构不发生重构。但理论上，层间的范德华相互作用会倾向于形成能量更低的AB或BA型贝尔纳堆垛，并与层内的晶格弹性形变相互竞争，在特定条件下可能导致原子结构发生重构，形成由公度域（AB/BA堆垛）和畴界组成的结构。这种原子尺度的重构预计会显著改变体系的对称性和电子结构。本研究旨在通过实验直接探测并阐明小扭转角（特别是小于约1°）下TBG中原子重构与电子重构的关联，及其对电子能带结构和输运性质的决定性影响，从而深化对vdW异质结界面工程的理解。

二、 研究流程 本研究采用了一种高度协同的多学科方法，结合了原子尺度结构表征、第一性原理计算和电学输运测量，对一系列具有精确可控扭转角（0° < θ < 4°）的TBG样品进行了系统研究。核心流程包括样品制备、结构表征、理论模拟和电学测量四大环节。

1. 样品制备与表征平台搭建：研究人员首先利用一种基于机械剥离和堆叠的微纳操作方法，精确制备了扭转角可控的TBG样品。随后，为了进行透射电子显微镜（Transmission Electron Microscopy, TEM）研究，他们将TBG样品夹在多层六方氮化硼（hBN）之间，并转移到薄的非晶氮化硅（SiN）支撑膜上。这种封装方式既保护了石墨烯，又使其适用于TEM观察。通过光学显微镜和显微探针技术，他们能够制备并研究尺度超过200纳米、包含超过10^3个摩尔单胞、尺寸变化小于10%的大面积均匀TBG区域。

2. 原子尺度结构表征（TEM与SAED）：这是研究原子重构的核心实验环节。

   • 暗场TEM成像：研究者利用TEM的暗场成像模式，选择性地过滤特定的石墨烯布拉格衍射斑点（g = 1010），来映射TBG的堆垛顺序。对于极小的扭转角（如0.1°），暗场图像显示出一系列周期性排列的三角形区域，它们具有交替的明暗衬度，对应于AB和BA型贝尔纳堆垛畴。这些尖锐的畴界结构直接证明了界面发生了原子重构，形成了公度域阵列。随着扭转角增大（如0.8°、1.2°），三角形畴的衬度减弱，逐渐转变为单向的条纹，这表明重构的强度减弱，结构趋近于未发生重构的简单摩尔图案。
   • 选区电子衍射分析：研究者系统地测量了不同扭转角TBG样品的选区电子衍射（Selected-Area Electron Diffraction, SAED）图。在重构的样品中（如θ=0.4°），他们观察到主布拉格峰周围出现了卫星峰。这些卫星峰的位置对应于由扭转角决定的摩尔超晶格倒易矢量。研究的关键定量分析在于测量卫星峰强度（I_sat）与相邻主布拉格峰强度（I_bra）的比值（I_sat/I_bra）。他们发现，对于11个不同θ的样品，I_sat/I_bra随θ增大呈指数衰减关系（I_sat/I_bra ~ e^(-αθ)，α ≈ 2.75 deg⁻¹）。当θ > 4°时，卫星峰强度低于检测限，表明在大角度下未发生显著重构。这一强度比值是衡量重构强度的关键实验指标。

3. 理论模拟与结构建模：为了深入理解实验观测到的衍射数据并反推原子构型，研究团队进行了有限元方法（Finite-Element Method, FEM）模拟。他们建立了一个结合层内弹性变形和层间范德华相互作用的连续介质-原子模型，模拟了不同扭转角下TBG的能量最小化结构。利用模拟得到的原子坐标，他们计算了理论SAED图案，其I_sat/I_bra随θ变化的趋势与实验数据高度吻合。更重要的是，FEM模型揭示了一个关键的转变：在特征交叉角θ_c ≈ 1°附近，结构重构行为发生显著变化。当θ < θ_c时，形成近乎公度的AB/BA畴，由宽度恒定（与θ无关）的尖锐畴界（称为孤子，soliton）分隔。当θ > θ_c时，畴内原子堆垛显著偏离理想的贝尔纳堆垛，畴界变宽，从一种堆垛到另一种的过渡变得平缓。这种结构上的“孤子态”与“摩尔态”的交叉转变，是本研究的核心物理发现之一。

4. 电学输运测量：为了探测原子重构对电子结构的实际影响，研究者制备了双栅极（顶栅和底栅）封装的TBG霍尔巴器件，并测量了其电导随载流子密度（通过栅压调节）和垂直磁场的变化。

   • 能带结构与态密度：他们测量了代表性扭转角（如接近第一魔角的θ≈1.1°和接近理论第二魔角的θ≈0.5°）样品的电导G(n)。在θ≈1.1°的样品中，观测到了对应于单粒子能隙和关联莫特能隙的特征。在θ≈0.47°的样品中，电导最小值出现在n/n0 = 0, ±4, ±8处（n0为每个摩尔单胞一个电子），但未出现完全的能隙和热激活行为，暗示了新的狄拉克点（Dirac point）的出现。
   • 电子结构计算对比：基于FEM模型得到的重构原子结构，研究团队进行了紧束缚模型和从头算连续模型计算，得到了重构后的能带结构和态密度（DOS）。计算结果显示，在θ≈1.1°时，重构主要改变了单粒子能隙的大小和带宽，但未破坏魔角物理（如平带导致的DOS尖峰）的定性特征。然而，在θ≈0.5°时，重构彻底改变了电子结构：未重构模型预测的能带结构杂乱无章，DOS具有许多不稳定尖锐特征；而重构模型则给出了清晰的、具有二次狄拉克锥的能带结构，其DOS特征（在n/n0 = 0, ±4, ±8处的极小值）与实验测量的电导特征高度吻合。这强有力地证明了在θ < θ_c时，原子重构导致了显著的“电子重构”，出现了次级狄拉克光谱。
   • 一维拓扑通道输运：通过独立调节顶栅和底栅电压，研究者可以在TBG中施加垂直方向的位移场（D）。他们比较了三种情况：AB堆垛双层石墨烯（θ=0°）、大角度TBG（θ≈2.8° > θ_c）和小角度TBG（θ≈0.47° < θ_c）。在θ≈0.47°的样品中，当施加位移场打开AB和BA畴的能隙时，他们观测到纵向电阻R_xx随|D|增大而增加并趋于饱和的现象。这与理论预言的图像一致：在位移场下，AB和BA畴变成绝缘体，而它们之间的畴界处形成受拓扑保护的一维导电通道。电流被限制在这个通道网络中传输，其饱和电阻值与基于量子电阻（R_q = h/4e^2）的简单网络模型估算值量级相符。
   • 朗道扇形图与维度转变：在垂直磁场下测量θ≈0.47°样品的R_xx随载流子密度和磁场的变化（朗道扇形图），观测到了源于主狄拉克点和次级狄拉克点的多重分形量子霍尔特征（Hofstadter蝴蝶谱）。特别地，当每个摩尔单胞的磁通量等于一个磁通量子（B = 5.4 T）时，出现了水平的电阻特征。随着位移场|D|增大，这些量子霍尔特征被显著抑制。这被解释为电子输运从二维体输运模式向一维拓扑通道网络模式转变的证据，因为磁场下的朗道能级形成需要二维扩展态，而一维通道中的载流子被限制，不利于朗道能隙的形成。

三、 主要结果 1. 原子重构的实验证实与定量刻画：通过暗场TEM直接可视化了小扭转角（如0.1°）TBG中AB/BA畴的周期性三角形阵列，为原子重构提供了直接证据。通过SAED定量分析I_sat/I_bra比值，建立了重构强度随扭转角指数衰减的标度关系（I_sat/I_bra ~ e^(-2.75θ)），并确定了发生显著重构的角度上限约为4°。 2. 结构重构的交叉转变：结合FEM模拟，研究发现原子重构行为在特征角θ_c ≈ 1°处发生交叉转变。θ < θ_c时，体系进入“孤子态”，形成清晰的公度畴和尖锐、宽度恒定的孤子边界；θ > θ_c时，体系处于“摩尔态”，重构较弱，畴界模糊。这一转变发生在摩尔周期与畴界宽度可比拟时，是层间相互作用能与层内弹性能竞争的结果。 3. 电子重构与魔角物理的修正：电学测量和理论计算共同证明，原子重构深刻影响电子结构。在θ > θ_c范围（如第一魔角1.1°附近），重构会调整能隙大小和带宽，但不改变平带和相关物理的本质。在θ < θ_c范围（如0.5°附近），重构则导致电子结构发生根本性改变：未重构模型中预测的第二魔角附近的强关联特征被显著抑制，取而代之的是由重构诱导产生的次级狄拉克能带，其DOS峰值比第一魔角处低一个数量级以上，解释了实验中未观测到强关联现象的原因。 4. 一维拓扑通道网络的观测：在小角度TBG（θ < θ_c）中施加垂直位移场，实验上观测到了与一维拓扑边界态网络输运一致的电阻饱和行为。通过分析朗道扇形图中量子霍尔特征随位移场的演化，进一步证实了在高位移场下电子输运从二维体态主导向一维拓扑通道网络主导的维度转变。这实现了对TBG中拓扑导电通道的全局电学探测。

四、 结论与价值 本研究系统揭示了扭转双层石墨烯中原子尺度重构的存在、规律及其对电子结构和输运性质的决定性影响。主要结论是：当扭转角小于约1°的特征交叉角θ_c时，体系会发生从非公度摩尔结构向由孤子边界分隔的公度畴阵列的转变（结构重构）；这一重构进而引发显著的电子重构，使得简单的摩尔能带描述失效，并诱导出次级狄拉克能带；在垂直电场下，被能隙化的三角形畴之间的孤子网络形成受拓扑保护的一维导电通道。

其科学价值在于： 1. 深化了对vdW异质结界面物理的理解：明确了“刚性旋转”近似在极小转角下的局限性，强调了原子弛豫和重构在决定这些异质结构电子性质中的核心作用。 2. 为“魔角”物理提供了更完整的图像：指出只有在θ > θ_c且未发生剧烈重构的角度区域（如第一魔角附近），刚性模型预测的平带和强关联物理才能成立。在更小的角度下，关联效应因重构导致的能带展宽而被抑制。 3. 开辟了新的物态调控途径：展示如何通过控制扭转角和垂直电场，在同一材料体系中实现从二维关联电子态到一维拓扑通道网络的可控切换。孤子网络为研究和操控一维拓扑边界态提供了理想平台。 4. 提供了多尺度研究范本：该工作将原子尺度结构表征（TEM）、微观理论建模（FEM）和宏观电学输运测量有机结合，为研究其他二维摩尔材料体系提供了方法论范例。

五、 研究亮点 1. 发现并系统表征了结构重构的“孤子-摩尔”转变：首次通过实验定量确定了TBG中原子重构的强度标度律和临界转变角θ_c ≈ 1°，明确了两种不同的重构区域。 2. 建立了原子重构与电子重构的直接关联：通过结合实验与计算，清晰展示了原子结构弛豫如何从根本上改变小转角TBG的电子能带，导致次级狄拉克点的出现，并解释了为何在第二魔角附近未观察到强关联现象。 3. 实现了对一维拓扑通道网络的全局电学探测：通过位移场调控，在宏观器件尺度上观测到了源于孤子网络的拓扑保护一维输运特征，并将此与磁场下的维度转变证据相互印证。 4. 创新的多技术联用方法：成功将适用于原子尺度表征的TEM技术与需要高质量封装的低温电输运测量技术结合在同一批可控样品上，确保了结构-物性关联研究的可靠性与精确性。

六、 其他有价值内容 文章还详细描述了用于分析衍射强度的运动学近似、用于能量最小化计算和能带计算的FEM及密度泛函理论方法细节（见“方法”部分）。作者也讨论了之前有研究利用局域栅极形成的腔室效应观测到Aharonov-Bohm振荡来探测一维通道，而本研究则展示了在整个介观器件尺度上实现和探测这种拓扑通道网络的能力。这些细节为同行复现和深入研究提供了重要参考。