Gap structure and phase diagram of twisted bilayer cuprates from a microscopic perspective

该研究通过紧束缚晶格模型计算，揭示了扭转双层铜氧化物超导体的时间反演对称破缺态与范霍夫奇点位置的关联，阐明了扭转角、层间隧穿、掺杂及温度对序参数的影响，并计算了约瑟夫森临界电流以解释现有实验中的矛盾结果。

要点

这篇论文探讨了一个非常前沿且迷人的物理现象：当两层超导体像“千层酥”一样被扭转在一起时，会发生什么？

为了让你轻松理解，我们可以把这篇论文的研究过程想象成**“在微观世界里玩乐高积木和调音”**。

1. 背景：两个“乐高”层被扭在一起

想象你有两块非常薄的、由乐高积木拼成的正方形板（代表铜氧化物超导体）。

• 正常情况：如果你把这两块板完全对齐叠在一起（0 度），它们就像普通的三明治，电流可以顺畅流动。
• 扭转实验：科学家把上面那块板相对于下面那块板旋转了一个角度（比如转了 45 度）。这就好比把两块乐高板错开拼在一起，形成了一种新的、复杂的“莫尔条纹”图案（就像透过两层纱窗看东西时出现的波纹）。

之前的困惑：
几年前，理论物理学家预测，当这两层板扭转到特定角度（特别是 45 度左右）时，它们会进入一种神奇的**“时间反演对称破缺”（TRSB）**状态。

• 通俗解释：想象电流像水流。在普通超导体里，水流是单向的。在这种神奇状态下，水流会自发地产生微小的漩涡，就像水在杯子里自己转起来了，而且这种旋转打破了“时间倒流”的对称性（就像看录像带倒放，水流方向会乱套）。
• 实验的矛盾：有的实验组（如 Zhao 等人）在扭转的晶体中看到了这种“漩涡”的迹象（电流在 45 度时几乎消失了）；但另一组实验（如 Zhu 等人）用更薄的“乐高片”做实验时，却发现电流在 45 度时依然很大，根本没有消失。
• 核心问题：为什么同样的实验，结果却截然不同？

2. 科学家的“显微镜”：微观模型

这篇论文的作者们没有直接去造新的实验设备，而是用超级计算机构建了一个**“微观乐高模型”**。

• 他们把材料看作是一个个原子点阵（像乐高网格）。
• 他们模拟了电子在这些网格上跳跃、相互作用，以及两层板之间“隧穿”（电子穿过中间缝隙）的过程。
• 关键变量：他们像调音师一样，微调了三个旋钮：
  1. 扭转角度（转多少度？）。
  2. 掺杂量（板子里有多少电子？就像调节水的咸淡）。
  3. 层间连接强度（两层板贴得有多紧？隧道有多宽？）。

3. 主要发现：不仅仅是 45 度

通过计算，他们发现了一个非常有趣的“相图”（就像一张天气地图，显示不同条件下会出现什么天气）：

• 不仅仅是 45 度：以前大家只盯着 45 度看。但作者发现，“漩涡”状态（TRSB 态）的出现，其实取决于电子能级中的“范霍夫奇点”（Van Hove Singularity）是否正好落在化学势附近。

  • 比喻：想象一个停车场（电子能级）。当“范霍夫奇点”（一个特别拥挤的停车位）正好停在“化学势”（入口）旁边时，电子们就会开始疯狂地“跳舞”（形成特殊的超导态）。
  • 这个“停车位”的位置，不仅取决于你加了多少电子（掺杂），还取决于两层板贴得有多紧（隧穿强度）。

• 两种不同的“漩涡”：
  他们发现了两种看起来像“漩涡”但本质不同的状态：

  1. $d + id'$ 态：这是大家预期的那种，电流在 45 度时会急剧下降（像 Zhao 的实验）。
  2. $d + is$ 态：这是一种混合了“漩涡”和“普通水流”的状态。在这种状态下，即使扭转到 45 度，电流依然可以很大（像 Zhu 的实验）。

4. 解开谜题：为什么实验结果不同？

这是论文最精彩的部分。作者提出，实验结果的差异，可能源于“隧道”的宽度不同。

• Zhao 的实验（晶体）：使用的是大块晶体，层与层之间的接触可能比较“粗糙”或者距离较远，导致隧穿强度较弱。在这种条件下，系统倾向于进入那种电流会消失的"$d + id'$"状态。
• Zhu 的实验（薄片）：使用的是极薄的 flakes（薄片），放在基底上。作者推测，由于基底的不平整或制造过程中的微小变形，导致局部区域的隧穿强度变强了（就像把两层板压得更紧，或者隧道修得更宽）。
  • 比喻：当隧道变宽（隧穿增强）时，电子们不再执着于转圈圈（$d+id'$），而是开始混合着普通水流一起跑（$d+is$ 甚至纯 s 波），导致电流没有消失，看起来就像没有发生“时间反演对称破缺”。

5. 结论与意义

这篇论文告诉我们：

1. 没有绝对的“标准答案”：扭转双层超导体的行为非常敏感，稍微改变一下“贴得紧不紧”（隧穿强度）或者“电子多少”（掺杂），整个系统的状态就会发生剧变。
2. 解释了矛盾：之前的实验之所以打架，很可能是因为他们的样品在微观层面上，电子“隧道”的宽度不一样。
3. 未来的方向：要真正搞清楚这个现象，我们需要更精确地测量两层材料之间的“接触细节”，而不仅仅是看宏观的电流。

一句话总结：
这篇论文就像给两个打架的科学家做了一次“法医鉴定”，指出他们看到的“不同现象”其实是因为他们手里的“乐高积木”虽然看起来一样，但拼合的紧密程度（隧穿强度）不同，导致电子们选择了不同的“舞蹈步法”。这为未来制造更稳定的量子器件提供了重要的理论地图。

技术摘要

这是一篇关于**扭曲双层铜氧化物超导体（Twisted Bilayer Cuprates）**微观物理性质的研究论文。作者通过紧束缚晶格模型（Tight-binding lattice model），系统研究了扭转角、层间隧穿强度、掺杂浓度和温度对超导序参量及相图的影响，旨在解释近期实验中关于时间反演对称性破缺（TRSB）态观测结果的矛盾。

以下是该论文的详细技术总结：

1. 研究背景与问题 (Problem)

• 背景： 自 Can 等人预测扭曲双层铜氧化物中可能存在破坏时间反演对称性（TRSB）的 $d+id'$ 拓扑超导态以来，实验界试图探测该状态，但结果存在显著矛盾。
  • Zhao 等人（2023）在扭曲的双层 BSCCO 晶体结中观察到临界电流在 45°扭转角附近显著抑制，并观测到 $\pi$ 相移的约瑟夫森衍射图样，支持 TRSB 态的存在。
  • Zhu 等人（2021）在超薄 BSCCO 薄片构成的扭曲结中，发现临界电流对扭转角依赖很弱，且在 45°附近没有显著抑制，暗示不存在拓扑态。
• 核心问题： 造成这些实验差异的原因是什么？是样品差异（如晶体 vs 薄片）、实验条件（如隧穿势垒厚度、表面粗糙度）还是微观物理机制（如层间耦合强度）的不同？目前的连续介质模型缺乏微观细节，无法完全解释这些差异。

2. 方法论 (Methodology)

• 模型构建： 作者采用了一个微观晶格模型，包含两个相互扭曲的正方晶格层。
  • 哈密顿量： 包含层内最近邻和次近邻跃迁（$t, t'$）、层内最近邻吸引相互作用（$V$，用于模拟 $d$ 波配对）以及层间隧穿项（$g_{ij}$）。
  • 层间隧穿： 采用 phenomenological 形式 $g_{ij} = g_0 e^{-(r_{ij}-c)/\rho}$，其中 $g_0$ 为隧穿强度，$c$ 为层间距。
  • 对称性处理： 对于共格扭转角（Commensurate angles），定义莫尔超胞（Moiré unit cell）。研究发现扭曲双层系统的对称群为 $D_4$（而非单层正方晶格的 $D_{4h}$）。
• 计算方法：
  • 在平均场近似下，自洽求解超导能隙方程。
  • 序参量投影： 将收敛后的序参量投影到 $D_4$ 群的不可约表示（Irreps, 如 $A_1, B_1, B_2, E$ 等）上，以确定序参量的对称性。
  • TRSB 判定： 通过规范变换（Gauge transformation）检查实空间序参量的虚部是否可消除。若不可消除，则确认为 TRSB 态。
  • 物理量计算： 计算态密度（DOS）、子晶格谱函数、费米面以及约瑟夫森临界电流（$J_c$）。

3. 关键贡献与主要结果 (Key Contributions & Results)

A. 相图与序参量对称性

• 丰富的相图： 作者绘制了扭转角（$\theta$）、掺杂（$\nu$）和隧穿强度（$g_0$）的相图。发现系统存在多种超导态，包括：
  • TRS 态（时间反演对称）： 纯 $B_1$ ($d_{x^2-y^2}$)、纯 $A_1$ ($s$ 波)。
  • TRSB 态（时间反演对称破缺）： $B_1 + iA_1$、$B_1 + iB_2$、$B_2 + iA_1$ 等混合态。
• 对称性转变： 不同于早期连续模型预测的 $d+is$ 或 $d+id'$，微观模型揭示了更复杂的对称性。例如，在 45°附近，随着 $g_0$ 增加，系统可能从 $B_1 + iB_2$ 态转变为 $B_1 + iA_1$ 态，甚至在大 $g_0$ 下转变为纯 $A_1$ ($s$ 波) 态。
• Van Hove 奇点的作用： TRSB 态的稳定性与正常态中 Van Hove 奇点（VHS）相对于化学势的位置密切相关。当 VHS 靠近化学势时，倾向于形成 TRSB 态（如 $B_1 + iA_1$）。掺杂和隧穿强度通过移动 VHS 位置来调控基态。

B. 费米面与谱函数分析

• 通过计算子晶格谱函数，发现特定对称性的序参量倾向于在费米面的反节点（antinodal）区域具有高谱权重。
• 随着 $g_0$ 增加，费米面变得复杂，不同不可约表示（Irreps）的稳定性发生竞争和转变。例如，在 36.87°扭转角下，低 $g_0$ 时 $B_1+iB_2$ 占优，高 $g_0$ 时 $B_2+iA_1$ 或纯 $B_2$ 占优。

C. 约瑟夫森临界电流 ($J_c$) 与实验矛盾的解释

• 小隧穿强度 ($g_0$)： 在 $g_0$ 较小且掺杂接近 Van Hove 奇点时，计算出的 $J_c$ 随扭转角增加而急剧下降，在 45°附近几乎消失。这与 Zhao 等人 在晶体结中的实验结果一致，支持 TRSB 态（如 $B_1+iB_2$）的存在。
• 大隧穿强度 ($g_0$)： 当 $g_0$ 较大时，系统倾向于形成 $B_1+iA_1$ 态或纯 $A_1$ ($s$ 波) 态。在这些状态下，$J_c$ 对扭转角的依赖性很弱，且在 45°附近保持显著数值。这与 Zhu 等人 在超薄薄片实验中的结果一致。
• 结论： 实验差异可能源于有效层间耦合强度的不同。薄片实验中的表面粗糙度、隧穿势垒厚度变化或杂质可能导致局部隧穿增强，从而将系统推入 $s$ 波或 $B_1+iA_1$ 主导的区域，掩盖了拓扑 TRSB 态的特征。

D. 温度依赖性

• 不同不可约表示的序参量具有不同的临界温度 ($T_c$)。随着温度升高，TRSB 态（多分量）通常会先失去其中一个分量，转变为 TRS 态（单分量）。

4. 意义与结论 (Significance & Conclusions)

• 统一解释： 该研究提供了一个统一的微观框架，解释了为何不同实验组在扭曲铜氧化物中观测到了截然不同的临界电流行为。关键在于层间隧穿强度和掺杂水平决定了基态的对称性。
• 微观机制： 强调了 Van Hove 奇点位置和费米面拓扑结构在稳定 TRSB 态中的核心作用。
• 实验指导： 指出要确认拓扑超导态的存在，必须严格控制隧穿势垒的均匀性和厚度。如果隧穿过强，系统可能进入平庸的 $s$ 波或混合态，导致拓扑特征消失。
• 方法论贡献： 展示了利用 $D_4$ 对称性群投影技术处理莫尔超胞中复杂序参量的有效性，超越了简单的连续介质近似。

总结： 本文通过高精度的微观晶格模型计算，揭示了扭曲双层铜氧化物中丰富的超导相图，证明了层间隧穿强度是调控 TRSB 态稳定性的关键参数，并为近期关于该体系实验结果的矛盾提供了合理的物理机制解释（即不同实验条件下的有效隧穿强度差异导致了不同的基态对称性）。