来源论文: https://arxiv.org/abs/2606.16536v1 生成时间: Jun 20, 2026 16:01

平带超导体中的涡旋核尺寸与量子几何：基于Mielke棋盘格子的BdG自洽理论深度解析

0. 执行摘要

在传统的Bardeen-Cooper-Schrieffer (BCS) 超导理论中，超导相干长度 $\xi_{BCS}$（其在实空间中直接对应于超导涡旋核的物理尺寸）是由费米速度 $v_F$ 决定的：$\xi_{BCS} \sim \hbar v_F / \Delta_0$。然而，这一经典图景在**平带系统（Flat-Band Systems）**中彻底失效。平带具有无限大的单粒子有效质量，其能带色散完全消失，这意味着费米速度 $v_F = 0$，常规的 BCS 相干长度发散至无穷大。这引发了一个根本性的物理问题：在完全缺乏动能和能带色散的平带超导体中，超导序参量的空间非均匀性尺度（如涡旋核尺寸）究竟由什么物理机制决定？

针对这一前沿课题，M. Özkurt 与 M. Iskin 在其最新工作中，以吸引力 Hubbard 模型的 Mielke 棋盘格子（Mielke checkerboard lattice） 为原型体系，开展了系统性的理论与数值计算研究。该工作不仅在动量空间中推导出了相干长度的闭合解析表达式，证明其完全由布洛赫态的**量子度规（Quantum Metric）**和粒子-孔穴对称的配对密度决定；同时，作者在实空间中利用自洽的 Bogoliubov-de Gennes (BdG) 方程对有限大小圆盘上的单量子涡旋进行了精确求解。

本博客将对该论文进行深度的技术解析，重点阐明量子度规如何通过虚态带间跃迁过程产生有限的“几何配对质量”（Geometric Pair Mass），从而在没有能带动能的前提下支撑起局域的超导涡旋结构。我们将详细推导其理论框架，展示关键的自洽计算 benchmark 数据，并提供基于贝塞尔函数展开法的实空间自洽 BdG 算法复现指南，最后对该工作的物理局限性与应用前景进行批判性评述。

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

1.1 核心科学问题：色散缺失下的尺度涌现

当超导电性发生在色散完全平坦的能带中时，传统的输运和配对图景被颠覆。由于单粒子动能为零，Cooper 对无法通过常规的带内动力学进行传输。然而，量子力学波动性赋予了布洛赫波函数独特的几何结构，即定义在布洛赫流形上的复投影空间度规——量子度规（Quantum Metric）。正是布洛赫态的量子几何性质，使得 Cooper 对能够通过虚带间跃迁（virtual interband processes）获得非零的有效质量，进而产生有限的超流刚度（Superfluid Weight）和实空间相干尺度。本工作的核心任务在于证明：实空间自洽涡旋核尺寸 $\xi_v$ 是否真正等于由体相布洛赫态量子几何预测的相干长度 $\xi_B$，并厘清其在不同填充物和相互作用强度下的标度律。

1.2 理论基础：Mielke 棋盘格子与低能连续介质模型

Mielke 棋盘格子是一种典型的二子格（A, B）线图（Line-graph）晶格，其几何结构如图 1(a) 所示。它包含一个完全平坦的下能带以及一个与其在特定点相切的色散上能带。非相互作用的紧束缚哈密顿量在晶格动量空间中表示为：

其中，基矢为二组分亚晶格波函数 $\phi_{\sigma \mathbf{k}} = (c_{A\sigma\mathbf{k}}, c_{B\sigma\mathbf{k}})^T$，不含自旋的布洛赫哈密顿量可写为：

其中 $\tau_0$ 是单位矩阵，$\tau_{x,z}$ 是泡利矩阵。其系数定义为：

这里 $t$ 是跃迁强度，$a$ 是晶格常数。该能带的能量谱为 $\varepsilon_{s,\mathbf{k}} = d_{\mathbf{k}}^0 + s d_{\mathbf{k}}$（其中 $s=\pm$，$d_{\mathbf{k}} = |\mathbf{d}_{\mathbf{k}}|$）。由于在该结构中 $d_{\mathbf{k}}^0 = d_{\mathbf{k}}$，因此下能带（$s=-$）的能量恒为零（$\varepsilon_{-,\mathbf{k}} = 0$），即实现了完美的平带；而上能带（$s=+$）则具有二次色散关系，两带在动量空间的四个对称点 $\mathbf{k} = (\pm \pi/a, 0)$ 和 $(0, \pm \pi/a)$ 处发生二次接触（Quadratic Band Touching）。

为了精确提取量子几何的物理效应并建立解析图景，作者在二次接触点附近对哈密顿量进行二阶展开，构建了低能有效连续介质模型：

其中 $k^2 = k_x^2 + k_y^2$。对该哈密顿量进行对角化，得到：

这完美再现了平坦下能带与各向同性二次色散上能带的接触行为。

1.3 技术难点：奇异量子度规与紫外/红外截断的自洽匹配

在低能极限下，归一化矢量 $\hat{\mathbf{d}}_{\mathbf{k}} = \mathbf{d}_{\mathbf{k}}/|\mathbf{d}_{\mathbf{k}}| = \frac{1}{k^2}(k_x^2 - k_y^2, 0, 2k_x k_y)$。量子度规张量定义为布洛赫波函数参数导数的实部：

代入可得解析形式：

显然，当 $k \to 0$ 时，量子度规 $g_{ij}^{\mathbf{k}} \sim 1/k^2$ 呈现幂律发散。这种在带接触点处的数学奇异性是平带量子几何计算的核心难点：

1. 紫外截断（UV Cutoff）的物理标定：连续介质模型在高动量下发散，必须通过物理上不含拟合参数的方法确定动量截断 $k_c$。作者通过将连续谱中的状态数与实空间有限圆盘（半径 $R$）内包含的晶胞数 $N_c$ 进行严格匹配（State-Counting），即：

   $$N_c = \sum_{\mathbf{k}} 1 = \frac{A k_c^2}{4\pi} = \frac{A}{2a^2}$$

   从而唯一确定了紫外动量截断 $k_c = \frac{\sqrt{2\pi}}{a}$，对应的能量截断为 $E_c = 2ta^2 k_c^2 = 4\pi t$。这提供了一个无拟合参数的宏观-微观桥梁。

2. 红外截断（IR Cutoff）与有限尺寸效应：有限圆盘边界条件天然地对系统施加了离散的动量分辨，这相当于引入了一个正比于 $1/R$ 的红外截断。在弱耦合极限下，该红外截断会强烈地压制低能处的对数发散，导致数值解与体相解析解之间出现偏差。如何在计算中识别并定量分析这一红外截断效应是研究的关键技术瓶颈。

1.4 方法细节：体相解析推导与实空间 BdG 自洽循环

A. 动量空间 BEC 映射与相干长度公式

在极稀薄极限下（填充度 $F \to 0$），超导 Cooper 对可被近似映射为弱相互作用的玻色-爱因斯坦凝聚体（BEC），其玻色子密度为 $n_B = F_B / a^2$，有效相互作用为 $U_B \approx 2U$。对应的 BEC 治愈长度（Healing Length）具有如下形式：

为了求解玻色子质量 $M_B$，需要对双体束缚态（Two-body bound state）的逆有效质量张量 $M_{2b}^{-1}$ 进行计算。该张量由带内（Intraband）和带间（Interband）贡献共同组成：$M_{2b}^{-1} = (M_{2b}^{-1})^{intra} + (M_{2b}^{-1})^{inter}$。其中，带间贡献直接由量子度规 $g_{ij}^{\mathbf{k}}$ 决定：

在弱耦合极限下（$|E_b| \ll E_c$），通过解析积分布洛赫空间，作者推导出了各向同性的逆双体质量：

这里，$\ln(E_c/U)$ 项体现了由于带接触点的二次奇异性导致的带间贡献对数增强。将该质量代入 $\xi_B$ 的表达式，并引入普适的粒子-孔穴对称对密度插值公式 $F_B = F(1-F)/2$，最终得到了平带超导相干长度的闭合解析公式：

该公式表明，$\xi_B$ 在半填充（$F=0.5$）时取得极小值，而在能带边缘（$F \to 0, 1$）由于凝聚体配对密度枯竭而呈现 $1/\sqrt{F_B}$ 形式的发散。而对于相互作用强度 $U$，则仅呈现极其微弱的对数依赖性。

B. 实空间自洽 BdG 算法

为了在实空间中检验上述解析公式的正确性，必须求解包含空间非均匀序参量 $\Delta(\mathbf{r}) = \Delta(r)e^{-i\theta}$ 的自洽 Bogoliubov-de Gennes 方程：

自洽循环的步骤如下：

1. 基于当前迭代步的实空间波函数 $u_n^S(\mathbf{r})$ 和 $v_n^S(\mathbf{r})$，自洽更新配对势（Gap Equation）：

   $$\Delta(\mathbf{r}) = -U \sum_n u_n^S(\mathbf{r}) v_n^{S*}(\mathbf{r}) f(E_n)$$

2. 类似地，计算局部自旋向上与自旋向下密度，得到总电荷密度分布：

   $$\rho(\mathbf{r}) = \sum_{n,S} \left[ |u_n^S(\mathbf{r})|^2 f(E_n) + |v_n^S(\mathbf{r})|2 (1 - f(E_n)) \right]$$

3. 通过调整化学势 $\mu$，使得全局平均填充度精确锁定在目标值 $F$：

   $$\frac{1}{A} \int d^2 r \rho(\mathbf{r}) = F/a^2$$

4. 重复上述步骤直至满足收敛标准：$\max_{\mathbf{r}} |\Delta^{(i)}(\mathbf{r}) - \Delta^{(i-1)}(\mathbf{r})| / \max_{\mathbf{r}} |\Delta^{(i-1)}(\mathbf{r})| < 10^{-6}$。

2. 关键 Benchmark 体系、计算数据与性能分析

2.1 自洽体相基态数据（动量空间）

作为数值自洽的基准，图 2 展示了在零温下，利用动量空间自洽方程 (11)-(12) 求解得到的超导配对能隙 $\Delta_0$ 和化学势 $\mu$ 随相互作用强度 $U/t$ 及能带填充度 $F$ 的演化关系。这一热力学基态数据是所有空间非均匀计算的渐近边界条件（Asymptotic boundary conditions）。

• 配对能隙 $\Delta_0$：在弱耦合极限下（$U/t \to 0$），$\Delta_0$ 在半填充 $F=0.5$ 处呈现出明显的极大值，并向两端对称衰减。其渐近行为由 $\Delta_0 \approx \frac{U}{2}\sqrt{F(1-F)}$ 给出，这展现了极强的平带配对效率，完全没有传统色散能带中弱耦合超导的指数压制（BCS 典型的 $\Delta_0 \sim \hbar \omega_D e^{-1/N(0)U}$）。
• 化学势 $\mu$：在弱耦合下表现出完美的粒子-孔穴对称性，即以 $F=0.5$（$\mu=0$）为中心对称。在极稀薄极限 $F \to 0$ 下，无论耦合强弱，化学势均收敛至 $\mu \to E_b/2$，其中 $E_b$ 为双体束缚能。这一特征表明平带超导体在任意弱的吸引相互作用下都会直接进入强配对的 BEC 物理区间。

2.2 实空间单量子涡旋解与电荷密度重分布

本工作最显著的贡献在于，通过在有限圆盘（半径 $R=45a$）上施加 $2\pi$ 位相绕数（Phase winding），数值求解了单量子涡旋的实空间剖面。计算所得的代表性数据分析如下：

A. 序参量空间剖面 $\Delta(r)$ 随相互作用 $U/t$ 的变化（图 4a）

在填充度 $F = 0.49$ 时，绘制了相互作用从极弱耦合（$U/t = 0.25$）变到中等耦合（$U/t = 1.0$）的超导序参量径向分布：

• 在涡旋中心（$r=0$），由于相位奇异性，$\Delta(r=0) = 0$ 得到严格保证。
• 随着径向距离 $r$ 的增加，$\Delta(r)$ 迅速恢复，并在几个晶格常数内渐近趋近于体相自洽解 $\Delta_0$（由图中的水平黑色虚线标定）。这直接证明了本文实空间连续介质 BdG 框架与动量空间体相理论在边界处的完全自洽。
• 关键结论：随着相互作用 $U/t$ 的增强，$\Delta(r)$ 的空间恢复速度几乎没有发生显著改变。这一不同寻常的现象表明，涡旋核的几何尺寸对相互作用强度不敏感，这与公式 (14) 中对数依赖性 $\xi_B \propto \sqrt{\ln(E_c/U)}$ 的预测高度一致。

B. 实空间电荷局域密度 $\rho(r)$ 的响应（图 4b）

• 尽管超导序参量在涡旋中心被完全压制，但局部电荷密度 $\rho(r=0)$ 并未衰减至零，而是保持在一个可观的有限值。这归因于局域在超导涡旋核内部的 Caroli-de Gennes-Matricon (CdGM) 束缚态，它们在平带背景下部分填充了能隙。
• 随着相互作用 $U/t$ 的增大，涡旋中心的电荷密度耗尽（depletion）程度逐渐加深，这说明强相互作用能够更有效地将非超导准粒子排斥出涡旋核心区域。

C. 填充度 $F$ 调控超导恢复尺度（图 5）

当相互作用固定在 $U/t = 0.5$ 时，调控全局填充度 $F$：

• 当 $F = 0.49$（接近半填充）时，归一化序参量恢复最为陡峭，意味着涡旋核尺寸极其局域（$\xi_v < a$）。

• 当填充度向能带边缘靠拢（如 $F=0.18$ 或 $F=0.82$）时，超导序参量的恢复曲线变得异常平缓，表明涡旋核极度膨胀。这一行为定量地契合了 $1/\sqrt{F_B}$ 的发散标度律。

• 粒子-孔穴对称性：图 5(a) 显示，关于半填充对称的填充对（如 $0.18$ 与 $0.82$，$0.30$ 与 $0.70$）其归一化序参量剖面几乎完美重合。图 5(b) 进一步展示，对于 $F < 0.5$，涡旋中心的电荷密度呈现耗尽特征（$\rho(r) < \rho_0$）；而对于 $F > 0.5$，涡旋中心则表现为电荷积累（$\rho(r) > \rho_0$），两者的偏差曲线围绕 $\rho a^2 = 0.5$ 互为镜像，严格遵循了如下粒子-孔穴变换律：

  $$\rho(\mathbf{r})_{1-F} = a^{-2} - \rho(\mathbf{r})_F$$

2.3 核心性能对比：实空间提取与动量空间解析（图 6）

为了定量回答核心科学问题，作者从自洽得到的 $\Delta(r)$ 中提取涡旋核尺寸 $\xi_v$，提取标准采用无参数拟合定义：$\Delta(r=\xi_v) = \Delta_0 \tanh(1/\sqrt{2}) \approx 0.61 \Delta_0$。图 6 将提取出的 $\xi_v$（红色圆点）与动量空间解析预测 $\xi_B$（蓝色实线）进行了直接对比：

这一 benchmark 对比无可辩驳地证明了：平带超导体的实空间非均匀尺度在本质上就是由布洛赫态的量子度规决定的。

3. 代码实现细节、复现指南与计算流程

本节提供一套完整的基于 Python 与贝塞尔展开法的实空间自洽 BdG 求解器设计指南，以便研究人员高效复现该工作。

3.1 理论矩阵元解析计算：Lommel 积分的威力

在圆盘几何（半径 $R$）下，单粒子哈密顿量的径向-角度微分算符通过基矢展开可转化为代数矩阵。其中，非对角块矩阵元 $B_{mm+2}^{jj'}$ 的核心部分为径向重叠积分：

这是一个经典的Lommel 积分。由于边界条件要求 $\beta_{jm}$ 为 $J_m(x)$ 的零点，利用 Lommel 积分公式，我们可以避开高成本的实空间数值积分，直接得到解析闭合矩阵元（对于 $m \neq -1$）：

这一解析转化使得构建超大维度的 BdG 哈密顿量矩阵变得极其迅速且精度极高。

3.2 自洽计算工作流与复现步骤

+-------------------------------------------------------------+
|                   输入参数: R, M, J, U, F                   |
+-------------------------------------------------------------+
                               |
                               v
+-------------------------------------------------------------+
|             计算贝塞尔函数零点并标定能级 (Eq. A7)            |
+-------------------------------------------------------------+
                               |
                               v
+-------------------------------------------------------------+
|           构造单粒子哈密顿量动能矩阵块 H_mm' (Eq. C17)        |
+-------------------------------------------------------------+
                               |
                               v
+-------------------------------------------------------------+
|            初始化体相超导能隙 Delta(r) = Delta_0            |
+-------------------------------------------------------------+
                               |
                               v
    +----------------------> [循环开始]
    |                          |
    |                          v
    |    +-----------------------------------------------------+
    |    |   使用当前 Delta(r) 计算配对矩阵块 Delta_mm' (Eq. D2)|
    |    +-----------------------------------------------------+
    |                          |
    |                          v
    |    +-----------------------------------------------------+
    |    |          组装完整的超大 Hermite BdG 矩阵             |
    |    +-----------------------------------------------------+
    |                          |
    |                          v
    |    +-----------------------------------------------------+
    |    |  调用 GPU (cuSOLVER) 进行全谱对角化，获取特征值/向量 |
    |    +-----------------------------------------------------+
    |                          |
    |                          v
    |    +-----------------------------------------------------+
    |    |        二分法调整化学势 mu 以精确锁定全局填充 F        |
    |    +-----------------------------------------------------+
    |                          |
    |                          v
    |    +-----------------------------------------------------+
    |    |          依据 Eq. E1-E3 更新 Delta(r) 和 rho(r)      |
    |    +-----------------------------------------------------+
    |                          |
    |                          v
    |    +-----------------------------------------------------+
    |    |               检查 Delta 相对误差是否 < 1e-6         |
    |    +-----------------------------------------------------+
    |                          |
    |            /\            |
    |           /  \-----------+
    |          / 否 \
    +----------\    /
                \  /
                 \/
                 | 是
                 v
+-------------------------------------------------------------+
|             收敛! 输出 Delta(r) 曲线并提取涡旋尺寸           |
+-------------------------------------------------------------+

3.3 核心 Python 伪代码实现

以下伪代码展示了基于上述算法的矩阵构建与自洽循环核心逻辑：

import numpy as np
from scipy.special import jn_zeros, jv
import scipy.linalg as la

def solve_bdg_vortex(R=45.0, M=60, J=50, U=0.5, target_F=0.49, t=1.0):
    # 1. 预计算贝塞尔零点
    beta = {} 
    for m in range(-M, M + 3): # 包含边界通道
        beta[m] = jn_zeros(m, J)
    
    # 2. 标定动量截断 (紫外物理截断)
    Ec = 4.0 * np.pi * t
    # 筛选满足等式 (A7) 的低能模式
    active_modes = []
    for m in range(-M, M + 1):
        for j in range(1, J + 1):
            energy = (beta[m][j-1] / R)**2
            if energy <= 4.0 * np.pi:
                active_modes.append((m, j))
                
    num_modes = len(active_modes)
    print(f"激活的单粒子模数: {num_modes}")
    
    # 3. 构造单粒子矩阵 H
    # 利用对角阵 (C1) 与解析 Lommel 矩阵元 (C17) 拼装
    H_single = np.zeros((2 * num_modes, 2 * num_modes), dtype=complex)
    # [此处省略子空间到全矩阵索引的具体映射代码]
    # 映射逻辑：根据选择定则 hmm' = 0 除非 m' = m 或 m +- 2 进行块稀疏填充
    
    # 4. 初始化 Delta(r)
    Delta_r = np.ones(1000) * 0.1 # 实空间径向网格
    
    # 5. 自洽循环
    tolerance = 1e-6
    error = 1.0
    mu = 0.0
    
    while error > tolerance:
        # A. 计算配对重叠积分 C_jj'
        # 依据公式 (D2) 在径向网格上数值积分
        C_matrix = construct_pairing_matrix(Delta_r, active_modes, beta, R)
        
        # B. 组装 BdG 矩阵 (2阶 Nambu 空间 x 2阶 亚晶格空间)
        # 维度为 4 * num_modes
        BdG_matrix = assemble_full_bdg(H_single, C_matrix, mu)
        
        # C. 求解本征值与本征向量 (推荐使用 cupy / PyTorch 在 GPU 上加速)
        eigenvalues, eigenvectors = la.eigh(BdG_matrix)
        
        # D. 二分法自洽寻找合适的化学势 mu 满足粒子数守恒
        mu = update_chemical_potential(eigenvalues, eigenvectors, target_F)
        
        # E. 更新实空间 Delta(r)
        new_Delta_r = calculate_new_gap(eigenvectors, eigenvalues, R)
        
        # F. 检查收敛性
        error = np.max(np.abs(new_Delta_r - Delta_r)) / np.max(np.abs(Delta_r))
        Delta_r = new_Delta_r
        
    return Delta_r

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

4.1 关键参考文献

1. Törmä, S. Peotta, and B. A. Bernevig, Nat. Rev. Phys. 4, 528 (2022)
   • 物理关联：该综述系统总结了平带系统超导电性与量子几何的理论奠基。本工作是将该综述中提及的“体相超流刚度”向“实空间非均匀局域结构”拓展的关键尝试。
2. S. Peotta and P. Törmä, Nat. Commun. 6, 8943 (2015)
   • 物理关联：开创性地证明了平带超导流刚度下限由布洛赫态的切尔恩数（Chern number）或量子度规决定。本工作推导双体配对质量时所采用的带间/带内分解法正是基于该理论框架。
3. M. Iskin, Phys. Rev. B 110, 144505 (2024)
   • 物理关联：提供了极稀薄极限下平带 BEC 映射和相干长度 $\xi_B$ 的原始推导。本工作直接继承了该文献中的对密度插值思想，并将其应用至实空间 BdG 对抗中。
4. F. Gygi and M. Schlüter, Phys. Rev. B 43, 7609 (1991)
   • 物理关联：奠定了利用自洽实空间 BdG 计算超导体单量子涡旋结构的算法基础。

4.2 局限性批判评述

尽管该工作极其优雅地将动量空间几何理论与实空间微观计算相统一，但站在量子化学与多体物理的前沿视角，该工作仍存在以下不容忽视的局限性：

1. 连续介质近似对微观格点结构的分辨率缺失

作者在实空间求解时，并非在原始的 Mielke 棋盘格子紧束缚格点哈密顿量上求解 BdG，而是采用了展开到二次项的低能有效连续介质算符 $H(\mathbf{r})$（即公式 19）。这一近似虽然允许了极其完美的圆盘圆柱对称性及贝塞尔函数解析展开，但代价是抹平了晶格尺度（$\sim a$）的细节。例如：

• 它无法描述超导涡旋核内可能出现的、由真实格点晶格畸变诱导的 Friedel 振荡以及微观电荷密度波（CDW）竞争。
• 当 $\xi_v < a$（如半填充附近）时，连续介质假设在物理上已经失效，必须要用离散格点 BdG 才能捕捉到真正的格点阻挫效应。

2. 平均场理论（BdG）对强涨落二维系统的过度乐观预测

本工作完全基于零温 BCS-BdG 平均场理论。然而，在二维（2D）极性平带超导体中，由于没有单粒子动能的支撑，相位涨落（Phase Fluctuations）将变得异常剧烈。其真实的超导相变温度并非由平均场能隙决定，而是由 Berezinskii-Kosterlitz-Thouless (BKT) 转变温度控制，即：

其中 $D_s$ 是超流刚度。在极稀薄填充或强耦合区，平均场理论会显著高估超导转变温度，且无法正确处理涡旋-反涡旋对的激发。因此，本文中关于涡旋核在实空间稳定存在的结论在有限温度下必须通过量子蒙特卡洛（QMC）或动力学平均场理论（DMFT）进行多体涨落校正。

3. 二次接触点（Quadratic Band Touching）的普适性限制

Mielke 棋盘格子的量子度规在 $k \to 0$ 呈现奇异发散，这是其逆双体质量出现对数项 $\ln(E_c/U)$ 的根本原因。然而，对于其他平带体系（例如 Lieb 晶格、Kagome 晶格或拓扑平带），能带接触点可能表现为狄拉克接触（线性发散 $g \sim 1/k$）或完全绝缘隔离。在那些体系中，相干长度随相互作用的标度律极有可能退化为非对数的形式（如直接正比于 $1/\sqrt{U}$）。因此，本工作的具体标度规律无法直接无缝推广到所有平带超流体系中。

5. 补充性深层物理剖析：布洛赫态的非局域性与重整化质量

为了帮助量子化学及计算材料学背景的读者更透彻地理解该工作，我们需要进一步阐明：为什么布洛赫态的几何性质能够产生质量？

5.1 从 Wannier 函数重叠看量子度规

在紧束缚近似下，晶格的局域性质由 Wannier 函数 $W_{\mathbf{R}}(\mathbf{r})$ 描述。量子度规 $g_{ij}^{\mathbf{k}}$ 在实空间中的物理对应，实际上是 Wannier 函数的空间局域化程度。一个关键的等式联系了量子度规在第一布里渊区的积分与 Wannier 函数的二阶空间矩：

在平带系统中，尽管单粒子的色散能带是平的，但由于晶格几何阻挫（Geometric Frustration），非平庸的拓扑或几何性质使得我们无法构造出同时具有正交性且在指数级局域化（Exponentially localized）的 Wannier 函数。Wannier 轨道不可避免地在空间中“蔓延”并相互重叠。这种空间重叠正是带间虚跃迁过程的物理起源：Cooper 对通过利用相邻 Wannier 轨道的空间尾部重叠，在空间中实现了量子隧穿，从而获得了有源自量子几何的非零配对质量 $M_B$。

5.2 实验观测展望：如何在真实体系中寻找几何涡旋？

随着实验技术的发展，人工平带体系在以下两个平台中取得了革命性进展，为观测本文预测的几何涡旋核提供了可能：

1. 冷原子光晶格平台：利用激光干涉技术，实验上已经成功构建了人工 Mielke 棋盘光晶格，并装载了具有吸引相互作用的 $^6\text{Li}$ 费米原子。通过调节 Feshbach 共振可以实现从弱耦合到强耦合的连续调控。由于冷原子体系具有极高的洁净度，利用原位成像（In-situ imaging）技术可直接提取物质度规分布，实验检测相干长度随填充度的 $1/\sqrt{F(1-F)}$ 发散趋势。
2. 魔角扭曲双层石墨烯（TBG）及过渡金属硫族化合物（TMDs）：这类莫尔超晶格（Moiré Superlattices）在特定魔角下会涌现出极为平坦的莫尔能带。其实空间超导相干长度极其微小，正好处于晶格常数尺度。通过扫描隧道显微镜（STM）直接对局域态密度（LDOS）进行空间扫描，观测超导涡旋内部 CdGM 束缚态的空间波函数演化，可直接对本文提出的自洽 BdG 密度剖面进行精确的实验印证。

通过将量子度规这一抽象的布洛赫空间几何概念转化为实空间可观测的超导涡旋尺寸，该研究不仅深化了我们对平带超导机制的理解，更为未来设计基于量子几何效应的高温超导器件开辟了全新的道路。