来源论文: https://arxiv.org/abs/2606.04958v1 生成时间: Jun 05, 2026 06:21

探索4e超导电性与分数涡旋：受挫约瑟夫森结骰子格点阵列的深度理论解析

0. 执行摘要

在凝聚态物理与量子信息科学的交叉前沿，寻找非常规超导相并实现对其拓扑分数激发的控制，是构建硬件级拓扑保护量子计算（如 Kitaev 表面码和 $\mathbb{Z}_2$ 规范场模型）的关键路径。传统的超导电性由电荷为 $2e$ 的库珀对凝聚驱动，而电荷为 $4e$ 的“库珀四重子（Cooper quartets）”凝聚则代表了一种全新的多体关联相。

本技术博客深度解析了一项关于在磁受挫约瑟夫森结骰子（Dice）格点阵列中诱导并确证 $4e$ 超导相和分数涡旋（半涡旋，Half-vortices）的突破性理论与数值模拟工作。该研究针对在每个菱形面元中注入 $f = 1/3$ 磁通量子（$\Phi_0/3$）的受挫 Dice 晶格，综合运用了绝热松弛动力学模拟（Relaxation Dynamics）、经典晶格蒙特卡洛（Monte Carlo）、二维无限张量网络（iPEPS & CTMRG）以及量子涨落二次玻色子微扰展开等前沿算法。

研究的关键结论包括：

1. 在 $f = 1/3$ 受挫点处，系统呈现出显著的、对约瑟夫森能量无序和磁通无序极具鲁棒性的临界电流非对称峰值，这与最新的实验观测高度吻合。
2. 局域超流体的谐波分解确证了在六角配位子格点上，第二谐波 $I_2$ 强度（对应 $4e$ 载流子）在 $f = 1/3$ 至 $f=1/2$ 的宽受挫区间内处于绝对主导地位，强度可达第一谐波 $I_1$ 的数倍。
3. 螺旋模量（Helicity Modulus）的对数尺寸修正拟合确立了相变温度为 $T_{BKT} \approx 0.227 E_J$，对应的相变类型是由半涡旋解禁闭驱动的非常规 Berezinskii-Kosterlitz-Thouless (BKT) 相变，其模量跳变点为 $8T/\pi$，而非传统 $2e$ 超导的 $2T/\pi$。
4. 二维无限张量网络（CTMRG）收缩算法精确提取了在 $T < T_{BKT}$ 的低温区内，常规 $2e$ 关联函数呈指数衰减（相干性被畴壁彻底破坏），而 $4e$ 关联函数呈幂律（代数）衰减，无可辩驳地确证了 $4e$ 幂律相长程秩序的存在。
5. 考虑有限充电能 $E_C$ 引起的量子起伏后，零点能微扰计算揭示了“无序诱导有序（Order-by-Disorder）”效应对零温基态简并度的解除作用，并预测了 $4e$ 相将在有限温度区间 $T^* < T < T_{BKT}$ 内稳定存在。

1. 核心科学问题、理论基础与技术难点

1.1 核心科学问题：多粒子凝聚与分数激发的自然涌现

自 BCS 理论创立以来，凝聚态物理学家长期致力于探索电荷为 $4e$、甚至 $6e$ 等更高阶多粒子复合体的凝聚状态。$4e$ 超流体不仅展现出独特的物理输运特性（例如在磁通偏置下呈现 $\pi$ 周期而非 $2\pi$ 周期的电流-相位关系），更重要的是，其拓扑缺陷是带有半个磁通量子的分数涡旋（Half-vortices, 环流为 $\pm 1/2$）。

从拓扑量子计算的角度来看，$4e$ 超流体可以被视为著名的 $\mathbb{Z}_2$ 规范场理论和 Kitaev 表面码的“经典前驱体”，因为半涡旋是与 Kitaev 表面码中磁性激发具有相同拓扑类别的拓扑激发。半涡旋与普通的 $2e$ 库珀对之间存在着拓扑非平凡的关系：当一个库珀对绕行半涡旋一周时，会获得一个 $e^{i\pi} = -1$ 的阿哈罗诺夫-玻姆（Aharonov-Bohm, AB）相位。这种非平凡的相位关系是构建硬件级抗噪量子存储器的物理基石。然而，如何在可扩展的物理器件中实现这种相，并在不需要人工引入 $\pi$ 周期相互作用（如人为设计的 Lee-Grinstein 模型）的前提下让其自然涌现，是一个长期的理论与实验难题。

1.2 理论基础：Dice 格点上的 XY 模型与零能畴壁

本研究将目光投向了具有特殊几何结构的骰子（Dice）格点阵列（如图 1 所示）。Dice 格点由两种截然不同的子格点构成：配位数为 6 的核心格点（Hub）和配位数为 3 的边缘格点。整个晶格由菱形面元（Rhombic plaquettes）拼接而成，三个菱形面元共同围绕一个配位数为 3 的顶点，拼成一个大六角形。

描述该系统在经典极限下的哈密顿量为受挫二维 XY 模型：

其中 $E_J$ 为约瑟夫森耦合能量，$\varphi_r$ 是格点 $r$ 处的超导相位，而 $A_{rr'}$ 是通过 Peierls 替换引入的矢势相位，其沿任意菱形面元 $p$ 的定向环流满足：

受挫参数 $f = \Phi / \Phi_0$ 代表穿过每个菱形面元的磁通量。

当 $f = 1/3$ 时，穿过每个大六角形（由三个菱形面元组成）的总磁通恰好为整数量子 $\Phi_0$。这一特定的磁通配置导致了极高的物理简并度。在零温极限下，系统的涡旋配置（每个六角形包含一个整磁通涡旋）可以完美地映射到三角晶格上的经典反铁磁伊辛模型（Ising Model）。在这一简并流形（degenerate manifold）中，不同的基态之间通过所谓的**零能畴壁（zero-energy domain walls）**相互隔离。这些畴壁穿过格点时不消耗额外的约瑟夫森能量，意味着其“线张力”（String Tension）在 $T^* = 0$ 时降为零。在有限温度下，畴壁发生无序增殖。在这些零能畴壁的交汇点，自然会暴露出分数磁通激发——即半涡旋，从而在有限温度下为 $4e$ 相的涌现提供了完美的土壤。

1.3 技术难点

1. 多亚稳态与相滑移阻碍物理收敛：由于 $f = 1/3$ 处基态具有极高简并度，普通的能量极小化方法极易陷入杂乱无章的局部亚稳态中。此外，在绝热改变外加相位偏置时，系统内会发生极其不规则的随机相位滑移（Phase Slips），这导致在大晶格上无法直接提取到干净且具有周期性的超流电流-相位关系。
2. 无限系统收缩与长程关联函数提取：要证明 $4e$ 超流相存在，必须证明在低温下，相算符 $e^{i\varphi_r}$ 的关联 $C_{2e}$（代表 $2e$ 库珀对）呈指数衰减，而 $e^{i2\varphi_r}$ 的关联 $C_{4e}$（代表 $4e$ 库珀四重子）呈幂律（代数）衰减。然而，普通的有限尺寸蒙特卡洛模拟受制于强烈的边界效应，极难在足够长程的尺度上将这两者清晰地分离。
3. 连续变量的张量网络表示：经典 XY 模型的相位变量 $\varphi$ 是连续的，将其表示为局域维度有限的二维张量网络需要非平凡的特征展开（Character Expansion）截断，并需要对庞大的二维网络进行高精度收缩（CTMRG）。

1.4 方法细节

为了攻克上述难点，本项工作设计了一套极具互补性的算法工具链：

• 绝热松弛动力学模拟（Relaxation Dynamics with Simulated Annealing）： 首先对系统进行高精度的模拟退火（Simulated Annealing），温度从 $\beta = E_J / T = 0.1$ 线性降温到 $\beta = 100$，经历 $10^6$ 个蒙特卡洛步。退火至目标温度后，执行局部能量最小化的松弛步（Relaxation Step, RS）：对三个子格点依次进行循环更新，通过直接令偏导数为零，将每个格点 $r$ 的相位 $\varphi_r$ 设为其邻居相位的非线性平均： $$\tan \varphi_r = \frac{\sum_{r'} \sin(\varphi_{r'} + A_{rr'})}{\sum_{r'} \cos(\varphi_{r'} + A_{rr'})}$$ 该更新自然满足局部基尔霍夫电流守恒（Eq. 4），能够极其逼真地模拟实验中的绝热扫描过程，捕获局域和全局的相滑移。
• 无限角转移矩阵重整化群（CTMRG）： 将配分函数 $Z = \int \mathcal{D}\varphi e^{-eta H}$ 映射为二维张量网络。利用修正贝塞尔函数对波尔兹曼权重进行傅里叶级数展开，将连续的连续相位积分化为离散的特征指数求和。通过 CTMRG 算法对该张量网络进行渐进收缩，引入环境键维度（Environment Bond Dimension） $\chi_E$ 作为变分截断参数（取 $\chi_E \in [200, 800]$），在热力学极限下精确逼近该二维系统的各类热力学量及长程关联函数。
• 螺旋模量的 Weber-Minnhagen 尺度修正拟合： 在蒙特卡洛框架下计算超流刚度（螺旋模量 $\rho_s$）。针对半涡旋激发引起的非常规 BKT 相变，使用重整化群流方程推导出的二阶有限尺寸对数修正方程对不同尺寸（自 $6 \times 6$ 至 $72 \times 72$）的模拟数据进行最小二乘拟合，其有限尺寸修正关系写为： $$\rho_s(L) - \rho_s^c = \rho_s^c \frac{1}{2\ln L + C}$$ 以此极其精确地外推出真正的热力学极限临界温度 $T_{BKT}$。

2. 关键 Benchmark 体系与数据深度解析

2.1 临界电流受挫谱：非对称峰与 AB 局域化（Fig. 2 & Fig. 3）

在 $50 \times 50$ 的 Dice 晶格阵列上，计算得到的全局临界电流 $I_c$ 随磁通受挫参数 $f$ 的依赖关系展现出了极具特征的物理图景：

• $f = 1/2$ 处的 AB 局域化山谷：在半受挫点 $f = 1/2$ 处，临界电流降至全局极小值。这是因为在该点处，围绕单个菱形面元传播的常规 $2e$ 库珀对会发生完全的破坏性干涉，即所谓的**阿哈罗诺夫-玻姆笼锁（Aharonov-Bohm Caging）**效应。这使得传统的 $2e$ 载流子被局域化在格点内部，无法参与大范围的电荷输运。这一现象是骰子格点区别于正方形和三角形格点最显著的特征之一。
• $f = 1/3$ 处的非对称峰值（如图 2 所示）：相反，在 $f = 1/3$ 处，$I_c$ 出现了一个明显的局域峰值，且该峰值呈现出高度的不对称性。这一特征与实验中观测到的超导切换电流谱完全吻合，暗示了非平凡电荷输运机制的开启。
• 能量密度的 Cuspy 极小值（如图 3 所示）：通过自由边界条件下的退火松弛，发现 $f = 1/3$ 处的能量密度 $\varepsilon$ 表现出尖锐的、导数不连续的极小值（Cusp）。这对应着零温基态下高度有序的涡旋排布。而 $f = 1/2$ 处对应的则是能量密度的全局最大值。

2.2 谐波分解：4e 载流子超流的绝对主导性（Fig. 5 & Fig. 6）

为了明确超流通道中的电荷载体，研究人员对局域超流电流进行了高阶谐波傅里叶分解。通过在两个相邻的 6 配位 Hub 格点（构成一个菱形的两端）之间引入局域相位差 $\delta\theta$，将输运电流展开为：

其中 $n=1$ 对应常规 $2e$ 库珀对输运， $n=2$ 对应 $4e$ 库珀四重子输运， $n=3$ 对应 $6e$ 复合体输运。

• 在 $f = 1/3$ 点处的 $I_2$ 暴涨（如图 5 所示）：谐波分析谱清晰地显示，随着 $f$ 接近 $1/3$，第一谐波 $I_1$ 的强度急剧下降，而第二谐波 $I_2$ 经历了一次极其强烈的抬升，其幅度的平方 $|I_2|^2$ 达到了 $|I_1|^2$ 的 5 倍以上。这表明，在 $f = 1/3$ 附近，穿过六角配位子格点的超流电流不再由库珀对驱动，而是几乎完全由库珀四重子（4e）介导。这一行为在其他子格点上并不存在，而是六角配位子格点的独有物理属性。
• 局域涡旋动力学的周期的性复原（如图 11 所示）： 这种由第二谐波主导的输运，在实空间中有着生动的分数涡旋动力学对应。在外加磁通偏置 $\Phi_{ext}$ 绝热递增的过程中，系统在大部分时间内处于无激发的简并基态。当偏置达到临界阈值时（如 $\Phi_{ext} \sim 0.28 \times 2\pi$），一对半涡旋和半反涡旋（图中表现为没有普通涡旋的蓝色和绿色六角形激发）会从 Hub 格点外侧核化产生。随着偏置继续增加，这对半涡旋绕行 Hub 岛屿进行旋转，直至偏置达到另一阈值时（$\sim 0.30 \times 2\pi$）发生重新湮灭。这一旋转和湮灭过程导致岛屿周围的涡旋配置发生了 $2\pi/3$ 的旋转，由于等效相位变化了 $\pi$，相当于在 Hub 格点周围引入了一个 Lee-Grinstein 畴壁圈。当外加偏置累积到 $2\pi$ 时，这一过程完整地发生了两次。这表明，在完整的 $2\pi$ 周期内，系统经历了两轮完美的核化-湮灭循环，从而在输运上表现出高度稳健的 $\pi$ 周期性，确证了 $I_2$ 谐波的主导地位。

2.3 无序强韧性（Fig. 7）

实际制备的超导器件中不可避免地存在几何面积涨落（导致磁通无序 $\Delta f$）以及约瑟夫森结屏障涨落（导致耦合能无序 $\Delta E_J$）。

• 约瑟夫森无序的极高耐受度：引入极其严重（$\Delta E_J = 0.3 E_J$）的随机无序后，临界电流在 $f = 1/3$ 处的峰值依然完好无损，仅发生了大约 $12\%$ 的轻微压低（图 7 桃红色曲线）。这表明该 4e 相在无序扰动下具有极高的局域鲁棒性。
• 磁通无序修正虚假极小值：引入标准差为 $\Delta f = 0.04$ 的高斯磁通无序后，计算得到的临界电流曲线在 $f \sim 0.31$ 处的局域极小值被成功抹平（图 7 蓝色实线），使得理论预测的极小值移动到了 $f = 1/2$。这一无序修正效果与实验阵列的实际输运特征更加吻合，从侧面证明了实验器件确实处于中等强度的磁通无序扰动下，而 4e 相在真实实验环境下能够稳定存在。

2.4 残余熵的精准定量（Fig. 8）

利用张量网络（iPEPS）方法，研究者精确计算了 $f = 1/3$ 处的零温熵密度。为了剥离连续超导相位涨落带来的发散背景，研究人员将 $f = 1/3$ 下的熵与无磁通受挫的 $f = 0$ 状态进行相减，提取出纯粹由涡旋物理简并引起的伊辛式离散熵贡献 $\Delta S = s_{f=1/3} - s_{f=0}$。

• 结果显示，在温度降至算法收敛的极限温度 $T/E_J \sim 0.1$ 时，熵差曲线（图 8）平滑且精准地收敛到了 Wannier 针对三角晶格反铁磁伊辛模型给出的解析残余熵常数： $$S_I = \frac{3}{\pi} \int_0^{\pi/6} \ln (2\cos\omega) d\omega \approx 0.323$$ 这一完美符合，在热力学极限下证实了 $f = 1/3$ Dice 晶格在极低温下的基态流形具有与三角晶格反铁磁伊辛模型相同的宏观简并结构。

2.5 螺旋模量相变温度拟合（Fig. 12）

在常规的 $2e$ 超导体中，BKT 相变是由单位强度的整数磁通涡旋（Vorticity $v = \pm 1$）解禁闭引起的，其刚度跃变满足 Nelson-Kosterlitz 公式：

然而，在 $4e$ 超导相中，最基本的解禁闭激发是半涡旋（Vorticity $v = \pm 1/2$）。由于涡旋之间的相互作用强度正比于其荷的平方（$v^2$），半涡旋的势能仅为整数涡旋的 $1/4$。因此，由半涡旋解禁闭驱动的非常规 BKT 相变，对应的螺旋模量刚度跃变应当陡增至：

通过蒙特卡洛计算不同晶格大小下的 $\rho_s$，并用 Weber-Minnhagen 对数修正拟合方程（Eq. C6）进行拟合（图 12），误差极小值（Inset）指向了极其精确的临界相变温度：

这为基于半涡旋去禁闭机制驱动的非常规 BKT 转变提供了强力的数值判据。

2.6 关联函数的长程行为：$2e$ 指数衰减与 $4e$ 代数衰减（Fig. 13 & Fig. 14）

张量网络算法（CTMRG）在热力学极限下的模拟彻底揭示了该系统的双重长程相干行为：

• 常规 2e 相干性的彻底丧失（Fig. 13）： 在 $T = 0.17 E_J < T_{BKT}$ 的超导温度下，常规 Cooper 对相干算符关联函数 $\tilde{C}_{2e}(\mathbf{r}) = \langle \cos(\varphi_r - \varphi_0) \rangle$ 表现出标准的指数衰减（Exponential decay）。这一物理行为的本质是，即便在超导温度以下，系统中已经充斥着大量的零能畴壁。任何一条畴壁的跨越，都会导致超导相位发生 $\pi$ 的突变，从而瞬间将 $2e$ 的远距离相位关联抹杀。
• 常规 4e 长程关联的保存（Fig. 14）： 相反，四重子相干算符关联函数 $\tilde{C}_{4e}(\mathbf{r}) = \langle \cos(2(\varphi_r - \varphi_0)) \rangle$ 在同样的温度和空间范围内，表现出极其稳健的幂律衰减（Power-law decay / Algebraic decay）。因为对于 $4e$ 算符而言，其对应的相位因子是 $2\varphi_r$。即使超导相位跨越畴壁时改变了 $\pi$，乘以 2 后便变为 $2\pi$。这意味着 $4e$ 相位关联函数对零能畴壁的增殖完全免疫！在张量网络收缩中，研究人员展示了随着环境维度 $\chi_E$ 从 200 增大到 800，幂律代数衰减的行为可以在越来越宽的空间尺度上保持。外推证实（Eq. 8），在热力学极限下关联长度 $\xi$ 发散，确证了 $4e$ 长程相干是热力学极限下的真实物理相。

3. 代码实现细节、复现指南与开源工具链

为了帮助科研工作者能够快速复现本项工作中的核心数值结果，以下梳理了详细的代码实现思路、算法流程和适用的开源计算资源。

3.1 松弛动力学算法实现指南

松弛动力学是用于获取临界电流谱（Fig. 2 和 Fig. 7）的核心。可以使用 C++、Julia 或带有 Numba 加速的 Python 实现以下核心逻辑：

1. 格点建立与边界约束： 定义 Dice 晶格：建立三套相互嵌套的子格点。在 y 方向的底部和顶部放置 Source 和 Drain 岛屿，其相位分别锁定为 $0$ 和 $\Phi_{ext}$。根据磁受挫参数 $f$，在各菱形边上分配 Peierls 替换相位 $A_{rr'}$，确保沿每个菱形面元的回路积分满足 Eq.(2)。
2. 退火预平衡： 使用经典极小化退火：初始温度 $\beta = 0.1$，最终 $\beta = 100$。采用经典 Metropolis 算法更新每个内部格点的相位，确保接受率稳定在合理区间。
3. 绝热相位扫描与松弛更新： 锁定外接岛屿相位：$\varphi_{source} = 0, \varphi_{drain} = \Phi_{ext}$。执行局部能量最小化的松弛步（Relaxation Step, RS）：通过直接令能量偏导数为零（Eq. A1），迭代更新内部格点的相位 $\varphi_r$： $$\tan \varphi_r = \frac{\sum_{r'} \sin(\varphi_{r'} + A_{rr'})}{\sum_{r'} \cos(\varphi_{r'} + A_{rr'})}$$ 重复该迭代直到全系统最大的相位改变量 $\delta\varphi < 10^{-7}$。
4. 电流积分与临界电流提取： 绝热扫描外接偏置 $\Phi_{ext}$（自 $0$ 至 $110 \times 2\pi$，步长 $\delta\Phi = 2\pi / 500$）。在每步相位稳定后，计算流向漏极的总电流 $I_{tot}$： $$I_{tot} = \frac{2e}{\hbar} E_J \sum_{r \in Drain, r' \in Bulk} \sin(\varphi_r - \varphi_{r'} - A_{rr'})$$ 在该绝热扫描序列中提取的最大电流 $I_{tot}$，即为对应受挫度 $f$ 下的全局临界电流 $I_c(f)$。

3.2 张量网络（iPEPS/CTMRG）复现技术栈

由于经典 XY 模型的相位 $\varphi$ 连续，无法直接投入传统的离散物理张量网络。其关键解决方案是引入特征展开（Character Expansion）：

1. 特征分解物理张量节点： 将两两格点之间的指数项利用修正第一类贝塞尔函数进行傅里叶分解： $$e^{\beta E_J \cos(\varphi_r - \varphi_{r'} - A_{rr'})} = \sum_{m=-\infty}^{\infty} I_m(\beta E_J) e^{im(\varphi_r - \varphi_{r'} - A_{rr'})}$$ 其中 $I_m(x)$ 为修正第一类贝塞尔函数。在这个级数中，求和指数 $m$ 自然充当了张量网络收缩的辅助键（Bond）。在数值实现中，可以根据温度低截断在合理区间内（例如 $m \in [-M, M]$，通常在低温下 $M=10$ 至 $20$ 已足够精准）。
2. 构建顶点张量： 在 Dice 格点的 Hub 节点和三配位节点上，对连续相位变量 $\varphi$ 进行局域高斯积分。由于： $$\int_0^{2\pi} \frac{d\varphi}{2\pi} e^{i(m_1 + m_2 + \dots) \varphi} = \delta_{m_1 + m_2 + \dots, 0}$$ 该 Kronecker delta 限制将连续变量完全转化为离散的张量。
3. CTMRG 收缩与能隙拟合： 利用角转移矩阵重整化群（CTMRG）算法对该张量网络进行渐进收缩，截断环境键维度（Environment Bond Dimension） $\chi_E \in [200, 800]$，在外推中提取第一能隙 $1/\xi$ 和次导能隙 $\delta$，验证 $\xi$ 的发散。

3.3 开源工具推荐

• ITensor (Julia/C++)：功能极其强大且对物理学家极度友好的张量网络库，内置了高效的 MPS/MPO，且能够方便地推广至二维 PEPS 的变分收缩。其 Julia 语言版本具有优异的高性能计算和多核并行支持。
• peps-torch：基于 PyTorch 实现的二维 PEPS 算法库，非常适合利用自动微分进行 CTMRG 收缩，可作为离散 XY 张量网络的基础框架。

4. 关键引用文献与局限性评述

4.1 关键引用文献

1. Korshunov (2005): S. E. Korshunov, Phys. Rev. Lett. 94, 087001 (2005). 这是本项工作最直接的理论基石。Korshunov 首次预言了在 $1/3$ 受挫的 Dice 格点上，零能畴壁的增殖会导致分数涡旋激发，并提出该点处可能存在非常规 $4e$ 超导相。
2. Lee & Grinstein (1985): D. H. Lee and G. Grinstein, Phys. Rev. Lett. 55, 541 (1985). 经典的文章，引入了具有额外的 $\pi$-周期相互作用项的经典超流模型，作为多粒子超流体的标准理论范式。
3. Bondar et al. (2025): J. D. Bondar et al., arXiv:2510.07412 (2025). 提供了本项工作最直接的实验原型——基于半导体-超导体异质结阵列的可调 Dice 晶格约瑟夫森结器件输运测量。

4.2 本工作理论局限性批判

尽管本项工作通过极其精湛的数值算法，在经典极限下确立了 $4e$ 相的长程幂律秩序，但从前沿科研和实际器件应用的角度审视，仍存在以下重要理论局限：

1. 忽略长程库仑相互作用与局部充电能： 本项工作中的经典 XY 模型完全忽略了各个超导岛屿的自充电能（Self-Charging Energy $E_C$）。在超导-半导体混合结阵列中，随着栅极电压的降低和岛屿尺度的微纳化，局部充电能 $E_C = 4e^2 / C^{self}$ 会显著增加。虽然论文在第 IV 节尝试通过二次谐波 LC 微扰展开来评估其对基态零点能 $E_{zp}$ 的修正（Eq. 12），但这仅仅是一种简化的局域谐振近似。真实的量子起伏会导致电荷-相位共轭不确定性关系（Eq. 11）的完全开启，激发长程库仑相互作用，并驱动系统发生向 Mott 绝缘体的量子相变，从而彻底摧毁 $4e$ 超流相。
2. “无序诱导有序”在极低温下的双刃剑效应： 论文指出，有限的充电能 $E_C$ 会在零温下解除经典基态的简并度，导致系统在 $T=0$ 时选择具有特定空间对称性的有序态（如蜂窝状或条纹状涡旋晶格），而 $4e$ 相仅在有限温度窗口 $T^* < T < T_{BKT}$ 内存在（图 15）。这意味着在实验上，如果为了抑制充电能而将温度降得过低，系统反而会退化为传统的具有 $2e$ 超流相干性的有序晶格态，这极大地压缩了实验上直接观测纯净 $4e$ 相的有效温度窗口（预计在 10mK 到 500mK 之间）。
3. 局域空间物理探针的缺失： 当前的数值模拟（如螺旋模量和关联函数）都是基于全局热力学和输运变量。在实验上，全局输运测量（如临界电流）由于受到热激活相位滑移（TAPS）和外接导线噪声的干扰，往往极难将纯净的 $4e$ 输运与多通道 $2e$ 协同输运（即两个不相关的库珀对同时随时间滑移）清晰地区分开。要彻底证实半涡旋的存在，亟需在实空间中利用高分辨率**扫描 SQUID 显微镜（Scanning SQUID Microscopy）**直接对这些被限制在畴壁交汇处的分数磁通进行局域直接成像。

5. 其他必要补充：拓扑量子计算与硬件级保护的展望

5.1 从半涡旋到非阿贝尔统计：$\mathbb{Z}_2$ 规范场与 Kitaev 表面码

在传统的超导量子比特中，主要的去相干源来源于环境电磁噪声对约瑟夫森相位的扰动，以及单个非平衡准粒子的隧穿。为了对抗这些去相干源，目前的行业主流是在软件算法层面上构建庞大且高冗余度的表面码（Surface Codes），这需要成千上万个物理量子比特才能融合成一个逻辑量子比特。

而本项工作揭示的 Dice 格点 $4e$ 超流相和半涡旋（分数涡旋），则在硬件物理层面上提供了一种自然的拓扑防护机制。在 $4e$ 超流相中：

1. 电荷量子化倍增：系统的物理算符具有 $\pi$ 周期性而非 $2\pi$ 周期性，这意味着任何单电荷（如单电子准粒子 $1e$）的非弹性隧穿过程在能谱上都被极大地抑制。
2. 拓扑阿哈罗诺夫-玻姆保护：一个常规的 $2e$ 库珀对在绕行半涡旋一周时，会获得一个精确的 $e^{i\pi} = -1$ 相位。在量子层面，这允许我们将量子逻辑态编码在这些半涡旋激发的拓扑排布中。由于半涡旋是非局域的拓扑缺陷，任何局域的微扰（如局域的温度抖动、弱磁场起伏）都无法在不改动全局拓扑同伦类的情况下改变该系统的相位态。这使我们能够在微观器件中构建出等效的 Kitaev 表面码 或 $\mathbb{Z}_2$ 拓扑有序态。在这种物相下，量子态的写入、读取和保持都是通过操纵分数半涡旋的核化与运动来实现的，具有极高的硬件级容错能力。

5.2 混合超导-半导体平台上的可编程超导集成电路

随着近年来微纳制备工艺（特别是基于原位铝外延生长的二维 InAs 异质结）的飞跃，利用电极栅极电压对约瑟夫森耦合强度 $E_J$ 进行实时、局域、可控的电学调节已经成为现实。

利用这一高度可扩展的平台，我们不仅可以静态地观察 4e 超导相，更可以构想如下的可编程拓扑量子芯片方案：

• 栅极控制畴壁和分数涡旋阵列编织：通过在 Dice 晶格的不同区域施加局部栅极电压，可以动态地压低特定约瑟夫森结的 $E_J$，从而实现对特定零能畴壁的创建、引导和锁定，甚至通过精细的设计，实现半涡旋在纳米晶格中的位置交换（Braiding）操作。这一技术的发展将大大简化量子比特容错门的设计难度，为未来规模化拓扑超导量子芯片的实现提供关键的物理路径。