来源论文: https://arxiv.org/abs/2607.06462v1 生成时间: Jul 08, 2026 15:57

解密铜氧化物超导 $T_c$ 控制密码：基于 Emery 模型与 CDMFT 的双机制调控深度解析

0. 执行摘要

自1986年铜氧化物高温超导（High-$T_c$ Cuprates）被发现以来，探索决定并限制其超导转变温度 $T_c$ 的微观物理机制，一直是凝聚态物理和量子化学领域最核心、最困难的课题之一。尽管实验上已经确定了 $T_c$ 与诸多物理量（如超流密度、自旋涨落、电荷能隙、晶格畸变及轨道占据数）之间的关联，但如何从一个统一的、不依赖特定材料微观细节的理论框架出发，定量提取并优化调控 $T_c$ 的普适机制，仍极具挑战性。

本项研究聚焦于铜氧平面的经典微观模型——三带 Emery 模型。研究人员采用非微扰的胞腔动力学平均场理论（CDMFT），结合连续时间量子蒙特卡洛（CT-HYB）杂质求解器，在 Zaanen-Sawatzky-Allen (ZSA) 相图框架下系统扫查了超导临界温度的最大值 $T_c^{\text{max}}$ 对模型微观参数（铜氧能级差、铜位库仑排斥能、空穴掺杂度）的依赖关系。研究不仅证实了 $T_c^{\text{max}}$ 在金属-电荷转移绝缘体边界附近达到最优的传统观点，更意外发现在深入电荷转移绝缘体区（即强关联极限下）同样存在未曾被系统探索的 $T_c^{\text{max}}$ 优化区间。

最核心的科学贡献在于，该研究成功将抽象的模型参数转化为两个可直接观测的物理量：电荷转移能隙大小（Charge-transfer gap size, $\Delta_I$）与氧空穴含量（Oxygen hole content, $2p_p$）。研究表明，$\Delta_I$ 与 $2p_p$ 是独立控制 $T_c^{\text{max}}$ 的双重微观机制：当电荷转移能隙减小且氧空穴含量增加时，$T_c^{\text{max}}$ 呈现单调递增趋势；而梯度场分析进一步定量表明，氧空穴含量是决定 $T_c^{\text{max}}$ 变化的最主导因素。这一成果为理解空穴掺杂铜氧化物超导材料的实验趋势提供了统一的理论基石，并为超冷原子光学晶格模拟 Emery 模型提供了决定性的实验预测。

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

1.1 核心科学问题与理论背景

铜氧化物的高温超导电性发生在典型的二氧化铜（$\text{CuO}_2$）平面内。在微观描述上，最简化的单带 Hubbard 模型虽然能捕捉到莫特（Mott）绝缘体和超导配对的基本物理，但它忽略了铜氧平面的多轨道特征，尤其是氧 $2p$ 轨道与铜 $3d$ 轨道之间的电荷转移行为。因此，能够显式包含这一多轨道杂化特征的三带 Emery 模型，被公认为描述铜氧面电子行为的 canonical 理论框架。

Emery 模型包含极其庞大的多维参数空间，包括：

• 铜 $3d_{x^2-y^2}$ 轨道的在位能 $\epsilon_d$
• 氧 $2p_x, 2p_y$ 轨道的在位能 $\epsilon_p$
• 铜-氧最近邻跃迁（hybridization）强度 $t_{pd}$
• 氧-氧最近邻及次近邻跃迁强度 $t_{pp}, t'_{pp}$
• 铜位上的局部排斥库仑力 $U_d$（氧位排斥一般忽略或设为常数）

寻找这一庞大参数空间中使超导转变温度 $T_c$ 最大化的点（即 $T_c^{\text{max}}$），不仅是一项繁重的计算任务，更因强关联电子相互作用带来的非微扰特性而难上加难。以往的研究多局限于特定的材料参数线，无法提供一幅全景式的、能与实验观测物理量直接挂钩的调控路线图。

1.2 理论基础：Zaanen-Sawatzky-Allen (ZSA) 框架

为了使 Emery 模型的多维参数空间具有可解释性，本研究引入了经典的 ZSA 相图作为制图工具（Map）。在原子极限下（即 $t_{pd} \to 0$），该系统在无掺杂（每个单元胞内 5 个电子，相当于 1 个空穴）状态下的基态由两个关键能量参数决定：

1. 铜位局部排斥能 $U_d$：将一个电子加到已有空穴的铜位上所需的能量成本。
2. 裸电荷转移能 $\Delta = \epsilon_d + U_d - \tilde{\epsilon}_p$：将一个电子从满带的氧 $2p$ 轨道转移到铜的 Hubbard 上能带（UHB）所需的能量成本（其中 $\tilde{\epsilon}_p = \epsilon_p - 2t_{pp}$ 为重整化后的氧能级）。

根据 $U_d$ 与 $\Delta$ 的相对大小，ZSA 框架将关联绝缘体分为两类：

• 当 $U_d < \Delta$ 时，系统能隙由 $U_d$ 决定，属于 Mott-Hubbard 绝缘体（MHI）。
• 当 $U_d > \Delta$ 时，系统能隙由 $\Delta$ 决定，属于 电荷转移绝缘体（CTI）。

由于高 $T_c$ 铜氧化物母体均属于电荷转移绝缘体，本研究将参数探索范围限定在 $U_d > \Delta$ 的 CTI 区域，并靠近金属-绝缘体边界（MET-CTI 边界），重点研究空穴掺杂（Hole doping, $\delta = 5 - n_{\text{tot}}$）后的超导演化。

1.3 技术难点与方法细节：胞腔动力学平均场理论（CDMFT）

由于铜氧面存在极强的局域排斥 $U_d$，传统的单粒子能带理论（如密度泛函理论 DFT-LDA）无法正确描述其绝缘基态及超导态。而一维或小尺寸团簇的计算又难以精确处理二维极限下的超导配对。为此，本研究采用了非微扰的胞腔动力学平均场理论（CDMFT）。

1.3.1 团簇选取与空间几何

CDMFT 将无限维点阵系统映射为一个嵌入在自洽电子浴（bath）中的有限尺寸“团簇杂质模型”。为了处理 $d$ 波配对对称性（$d$-wave superconductivity），必须使用能够包容四个 Cu 位和相应 O 位的最小团簇。本研究使用了包含 12 个格点的团簇：

• $N_d = 4$ 个 Cu 格点
• $N_p = 8$ 个最近邻 O 格点

此 12 胞格点几何完美兼顾了空间关联的分辨率（特别是最近邻铜-铜对关联）与量子蒙特卡洛计算开销的可控性。

1.3.2 哈密顿量与 CDMFT 自洽环路

Emery 哈密顿量为：

其中，$\mathbf{C}^{\dagger}_{\mathbf{k}\sigma} = (d^{\dagger}_{\mathbf{k}\sigma}, p^{\dagger}_{x\mathbf{k}\sigma}, p^{\dagger}_{y\mathbf{k}\sigma})$ 为产生算符向量。非相互作用部分的紧束缚矩阵 $\mathbf{h}_0(\mathbf{k})$ 显式包含：

其中，$V_{dp_\alpha}$ 描述铜-氧杂化跃迁，$W_{p_\alpha p_\beta}$ 描述氧-氧跃迁。CDMFT 的自洽循环步骤如下：

1. 初始猜测：给定一个局域自能矩阵 $\mathbf{\Sigma}_c(i\omega_n)$。
2. 格林函数计算：在全布里渊区内积分，求得局域格林函数： $$\mathbf{G}_l(i\omega_n) = \frac{1}{N_k} \sum_{\mathbf{k}} \left[ (i\omega_n + \mu)\mathbf{I} - \mathbf{h}_0(\mathbf{k}) - \mathbf{\Sigma}_c(i\omega_n) \right]^{-1}$$
3. 有效浴表征：通过 Dyson 方程提取杂质模型的非相互作用格林函数（Cavity Green’s Function）$\mathbf{\mathcal{G}}_0^{-1}(i\omega_n) = \mathbf{G}_l^{-1}(i\omega_n) + \mathbf{\Sigma}_c(i\omega_n)$。
4. 杂质模型求解：使用**连续时间杂化膨胀算法量子蒙特卡洛（CT-HYB）**求解器，计算在非平凡局域排斥 $U_d$ 下的完整杂质格林函数 $\mathbf{G}_c(i\omega_n)$。
5. 自能更新：更新自能 $\mathbf{\Sigma}_{c}^{\text{new}}(i\omega_n) = \mathbf{\mathcal{G}}_0^{-1}(i\omega_n) - \mathbf{G}_c^{-1}(i\omega_n)$。
6. 收敛判断：重复步骤 1-5，直至 $\mathbf{\Sigma}_c$ 满足自洽收敛准则。

1.3.3 超导相的测量与费米子符号问题

为测定超导转变温度 $T_c$，研究引入了超导序参量：

其中 $i, j$ 代表最近邻的 Cu 位。由于 CDMFT 在二维网格上采用有限温度求解，超导转变对应的 $T_c^{\text{CDMFT}}$ 实质上代表了配对序开始显着建立的平均场温度。蒙特卡洛计算中的致命技术障碍是费米子符号问题（Fermion Sign Problem）。在极低温度下，蒙特卡洛平均符号急剧衰减至零。因此，本研究所达到的最低计算温度为 $T = 1/50$（以 $t_{pp} = 1$ 为能量单位）。

2. 关键 Benchmark 体系，计算所得数据，性能数据

2.1 物理体系及基准参数设定

在基准计算中，跃迁强度的基准设定来源于典型空穴掺杂铜氧化物的拟合数据，能量单位统一归一化为氧-氧最近邻跃迁强度 $t_{pp} = 1.0$：

• 氧-氧次近邻跃迁：$t'_{pp} = 1.0$
• 铜-氧跃迁：$t_{pd} = 1.5$
• 铜位在位能：$\epsilon_d = 0.0$
• 铜-氧能级差变化范围：$\tilde{\epsilon}_p - \epsilon_d \in [2.0, 10.0]$
• 铜位排斥力：$U_d \in [7.0, 16.0]$

2.2 ZSA 相图中的物理量提取

在无掺杂（$\delta = 0$，即每个单元胞包含 5 个电子）状态下，研究首先对整个 $(U_d, \Delta)$ 参数空间进行了网格化扫描。计算提取了两个非平凡物理可观测变量：

1. 物理电荷转移能隙 $\Delta_I$： 从系统的化学势 $\mu$ 随总电子数 $n_{\text{tot}}$ 演化的等温线 $\delta(\mu) = 5 - n_{\text{tot}}(\mu)$ 中提取。当发生不可压缩的绝缘体行为时，等温线上会出现明显的阶梯。通过定义平台宽度限制 $d n_{\text{tot}}/d\mu < 0.005$，定量读取平台宽度即为物理能隙 $\Delta_I$。
2. 氧空穴含量 $2p_p$： 定义氧位上的平均空穴数 $2p_p = 2(2-n_p)$，其中 $n_p$ 为氧轨道上的单粒子占据数。由于 $t_{pd}$ 导致的强共价杂化作用，空穴会在铜与氧之间共享（具有混杂的 $d-p$ 轨道特征）。

关键图表数据分析（图2 & 图3）：

• 在 $U_d - \Delta$ 平面内，**等能隙线（$\Delta_I = \text{const}$）与等氧空穴线（$2p_p = \text{const}$）**呈现非平行的几何特征。等氧空穴线近乎垂直，表明氧空穴含量主要受裸电荷转移能 $\Delta$（或能级差 $\tilde{\epsilon}_p - \epsilon_d$）控制；而等能隙线斜率明显。这直接揭示了 $\Delta_I$ 与 $2p_p$ 相互独立，并不存在一维正比关系。
• 局部态密度（LDOS, $N(\omega)$）展示（图3）清楚地表明，对于具有相同物理能隙（如 $\Delta_I = 0.6$）的系统，随着参数向强关联深区（大 $U_d$、小 $\Delta$）演化（点 1 至点 5），电荷转移带与上 Hubbard 带之间的混合程度大不相同，氧轨道的共价混合性在深区显着增强，表现为 $2p_p$ 的系统性上升。

2.3 超导转变温度 $T_c$ 与参数的依赖数据

对选定的 8 个典型母体绝缘状态（点 1 至 5 对应 $\Delta_I = 0.6$；点 6 至 8 对应 $\Delta_I = 0.9$）进行空穴掺杂扫描，获得了对应的超导圆顶（dome）结构（图4 & 图7）：

数据演化物理结论：

• 第一重趋势（固定能隙）：当保持物理电荷转移能隙 $\Delta_I = 0.6$ 不变，逐渐深入电荷转移绝缘区（从 1 到 5）时，随着氧空穴含量 $2p_p$ 从 0.250 提高到 0.330，最高超导转变温度 $T_c^{\text{max}}$ 实现了从 0.0209 到 0.0232 的单调递增。
• 第二重趋势（固定能级差）：对于固定的 $\tilde{\epsilon}_p - \epsilon_d = 6.0$（对比 3 号和 6 号点），当物理能隙 $\Delta_I$ 从 0.6 增加到 0.9（意味着系统更偏离金属-绝缘体边界）时，$T_c^{\text{max}}$ 从 0.0223 显着下降至 0.0207。

2.4 梯度场与双机制的联合定位分析（图12）

为了厘清 $\Delta_I$ 与 $2p_p$ 对 $T_c^{\text{max}}$ 影响的相对强弱，研究人员在 $(\Delta_I, 2p_p)$ 组成的二维物理参数空间中，对数据进行了网格差值和梯度场（Gradient Vector Field）计算：

• 二维面 $T_c^{\text{max}}(\Delta_I, 2p_p)$ 的等高线图显示，在较小 $\Delta_I$ 和较大 $2p_p$ 的区域，超导临界温度最高。
• 梯度矢量 $\nabla T_c^{\text{max}}$ 的方向几乎垂直于 $\Delta_I$ 轴。这表明，相对于能隙减小所带来的微弱贡献，氧空穴含量的增加是导致超导配对强度最主要的、占主导地位的驱动力。

3. 代码实现细节，复现指南，所用的软件包及开源 repo link

3.1 理论复现的技术栈体系

为了完整复现该工作，需要掌握并搭建两阶段的强关联理论计算工作流：

1. 胞腔动力学平均场（CDMFT）计算：用于处理 Emery 模型在 finite temperature 下的自洽环路。
2. 解析延拓（Analytic Continuation）：将蒙特卡洛在松原频率（Matsubara frequency, $i\omega_n$）上的虚时自能及格林函数数据映射到实轴（$\omega$），获得实频态密度（LDOS）。

核心计算软件体系：

• TRIQS (Toolbox for Research on Interacting Quantum Systems): 当前凝聚态物理学界主流的开源多体关联理论平台。它提供了高效率的 C++/Python 接口，极适合搭建多带 DMFT/CDMFT 自洽环路。
  • 开源链接: https://github.com/TRIQS/triqs
• TRIQS/cthyb (CT-HYB Solver): TRIQS 平台下的高性能三维连续时间杂化膨胀量子蒙特卡洛求解器，支持轨道杂化、非对角杂化阻尼等，是求解多轨道 Emery 杂质模型的利器。
  • 开源链接: https://github.com/TRIQS/cthyb
• TRIQS/maxent (Maximum Entropy Method): 用于对松原格林函数进行高精度的最大熵解析延拓。
  • 开源链接: https://github.com/TRIQS/maxent

3.2 CDMFT 核心复现逻辑及伪代码

以下是基于 TRIQS 逻辑复现 12-格点 Emery 模型自洽循环的核心 Python 伪代码骨架。读者可以参考其逻辑在本地集群或超算上配置并运行：

import numpy as np
from triqs.gf import *
from triqs.operators import *
from triqs.applications.impurity_solvers.cthyb import Solver

# 1. 物理参数设定
t_pp = 1.0
t_prime_pp = 1.0
t_pd = 1.5
U_d = 11.15         # 对应 3 号点
epsilon_d = 0.0
epsilon_p = 6.0 + 2*t_pp  # 重整化能级
beta = 50.0         # 逆温度 (T = 1/50)
mu = 14.8           # 调节化学势以控制掺杂度

# 定义 12-格点团簇结构 (4 Cu + 8 O)
# 这里简化定义杂质空间 GF 块
spin_names = ['up', 'down']
orb_names = [f'Cu_{i}' for i in range(4)] + [f'O_{j}' for j in range(8)]

# 2. 初始化自洽格林函数
g_loc = BlockGf(mesh=MeshImFreq(beta=beta, statistic='Fermion', n_max=1024), 
                gf_struct=[(spin, len(orb_names)) for spin in spin_names])

sigma_c = BlockGf(mesh=MeshImFreq(beta=beta, statistic='Fermion', n_max=1024), 
                  gf_struct=[(spin, len(orb_names)) for spin in spin_names])
sigma_c.zero()

# 3. 初始化 CT-HYB 求解器
# 实际运行中需要配置多核 MPI 并行
solver = Solver(beta=beta, gf_struct=[(spin, len(orb_names)) for spin in spin_names])

# 自洽循环主环路
max_iter = 40
tolerance = 1e-4
for iteration in range(max_iter):
    print(f"--- CDMFT Iteration {iteration} ---")
    
    # 步骤 A: 计算格子局域格林函数 (布里渊区数值积分)
    # 构造紧束缚哈密顿量矩阵 H_0(k)
    # g_loc_new = \int_BZ [ (iw + mu)*I - H_0(k) - \Sigma_c ]^-1 dk
    g_loc_new = calculate_lattice_green_function(H_0, sigma_c, mu, beta)
    
    # 步骤 B: 提取并设置杂质非相互作用格林函数 G_0^-1
    for spin, g_0 in solver.G0:
        g_0 << inv(inv(g_loc_new[spin]) + sigma_c[spin])
    
    # 步骤 C: 运行 CT-HYB 蒙特卡洛采样，求解相互作用格林函数
    # 配置在位库仑斥力 U_d（仅作用在 Cu 轨道上）
    h_int = calculate_local_interaction_operator(U_d, orb_names)
    solver.solve(h_int=h_int, n_cycles=1000000, length_cycle=200, warmup_cycles=10000)
    
    # 步骤 D: 计算新自能
    sigma_new = BlockGf(mesh=MeshImFreq(beta=beta, statistic='Fermion', n_max=1024), 
                        gf_struct=[(spin, len(orb_names)) for spin in spin_names])
    for spin, sig in sigma_new:
        sig << inv(solver.G0[spin]) - inv(solver.G[spin])
        
    # 步骤 E: 自能混合与收敛判断
    diff = calculate_difference(sigma_c, sigma_new)
    print(f"Self-energy difference: {diff}")
    if diff < tolerance:
        print("CDMFT Cycle Converged successfully!")
        break
    
    # 线性混合更新自能
    sigma_c << 0.3 * sigma_new + 0.7 * sigma_c

3.3 计算资源与收敛建议

• 计算开销：由于 Emery 模型的杂质格林函数矩阵维度为 $12 \times 12$，蒙特卡洛求解过程面临严重的维数灾难。复现该研究的数据点需要使用 MPI 并行。建议每个数据点配置 128 到 512 个 CPU 核心，单次 solve 循环至少需要迭代 $10^8$ 个蒙特卡洛步骤以获得统计误差极小的格林函数。
• 超导测量策略：当从正常态逐步降温至超导相变点附近时，由于存在热涨落，必须在杂质求解器中显式加入极微弱的超导配对场作为引线（seed），并在自洽循环收敛后撤去该场，观察超导序参量 $\Phi$ 是否能实现自维持自洽。

4. 关键引用文献，以及你对这项工作局限性的评论

4.1 关键引用文献综述

本项工作紧密建立在以下几项里程碑式研究的基础之上：

1. Emery, V. J. (PRL 1987) [33]：首次提出描述二氧化铜平面的三带 Hubbard 模型（即 Emery 模型），指明了氧空穴空腔在超导中的潜在主导地位。
2. Zaanen, J., Sawatzky, G. A., & Allen, J. W. (PRL 1985) [36]：确立了 ZSA 关联材料分类相图，为区分莫特绝缘体与电荷转移绝缘体提供了统一的标准。
3. Weber, C. et al. (EPL 2012) [21]：通过 LDA+DMFT 方法系统分析了高 $T_c$ 超导体的转变温度与电荷转移能的关系，指出减小电荷转移能会提升 $T_c$。
4. Kowalski, N. et al. (PNAS 2021) [24]：利用变分团簇近似深入探讨了氧空穴含量对多带超导体配对能隙的影响。
5. Fratino, L. et al. (PRB 2016) [22]：详细阐述了基于 CDMFT 方法在 CTI 区域计算超导转变温度与赝能隙边界的技术细节。

4.2 局限性深度评论（Critique）

尽管本工作在定量化 $T_c$ 的调控规律方面取得了极其优秀的突破，但作为面向前沿理论的学术工作，仍存在以下几点不容忽视的局限性：

1. 对 Kosterlitz-Thouless (BKT) 相位涨落的忽略

CDMFT 在本质上依然是一种“团簇级”的动力学平均场方法。它能够极好地处理团簇内部的短程空间关联与动力学自能，但由于其外部自洽浴依然是平均场性质的，计算中测得的超导温度 $T_c^{\text{CDMFT}}$ 实质上对应了配对幅度（Pairing amplitude）建立的平均场温度。 在真实的二维超导平面内，由于热起伏，超导的长程相干性受制于 BKT 物理（BKT phase fluctuations）。相位涨落会使实际超导转变温度 $T_{c,\text{BKT}}$ 显着低于 $T_c^{\text{CDMFT}}$。未来的工作必须结合大尺度路径积分蒙特卡洛，将相位涨落效应定量引入调控机制中。

2. 固定跃迁参数对真实物理复杂性的简化

为了解决三维多维空间优化的不可行性，研究将跃迁常数 $t'_{pp}$ 与 $t_{pd}$ 设为了固定常数（常数 1.0 和 1.5）。然而，在不同的高 $T_c$ 材料（如 LSCO、YBCO、Bi2212）中，晶格畸变、顶角氧原子的存在会显着改变这些跃迁强度的绝对值及相对比值。这使得该模型在定量上目前只能极好地拟合 LSCO 的物理行为（由于其氧空穴区间落在 $0.26 \sim 0.33$ 附近），而对 YBCO 等大空穴共价混合体系的定量预测偏保守。

3. 费米子符号问题对零温极限外推的限制

量子蒙特卡洛底层的符号问题导致研究无法直接计算低于 $T = 1/50$ 温度下的超导热力学性质。这迫使研究人员不得不依靠在极低掺杂区域的外推，或者直接使用最接近基态的虚轴数据作为替代品，一定程度上降低了临界区参数拟合的精度。

5. 补充物理剖析与未来展望

5.1 从多带 Emery 模型到单带 Hubbard 模型的映射边界

一个极其深刻的物理问题是：为什么我们不能仅仅使用更简单的单带 Hubbard 模型来优化 $T_c$？

在低能极限下，Zhang-Rice 单态（Zhang-Rice Singlet）理论允许我们将 Emery 模型的低能物理映射为等效的单带模型，其中氧 $2p$ 轨道上的空穴与铜 $3d$ 的局部自旋配对成非自旋的单态，并像单带空穴一样跃迁。然而，本项研究所揭示的最核心物理机制——氧空穴含量 $2p_p$ 是调控超导配对的最关键物理量，这一结论在单带 Hubbard 模型中是完全无法显式定义并观察到的。

当系统被掺杂时，空穴电荷并不会简单地呆在铜位上，而是极大地分布在氧位上。本项研究（图6 & 图8）的数据表明：

• 掺杂电荷的流向：空穴掺杂主要增加了氧 $2p$ 的空穴含量（曲线斜率极平缓，表明铜的空穴 $p_d$ 几乎不随掺杂改变，增加的掺杂电荷绝大部分注入了氧的 $2p$ 轨道）。
• 这说明超导配对的微观机制涉及极其强烈的、偏离 Zhang-Rice 极限的高能电荷涨落过程。氧位上的电荷涨落对于维持极高强度的 Cooper 配对势能具有不可磨灭的作用。

5.2 对超冷原子光学晶格模拟（Cold Atom Simulation）的指导

随着量子模拟技术的跨越式发展，利用超冷费米原子在光学晶格中人工合成 Emery 模型，已经在实验上成为现实（如近年来 Max Planck 研究所 Bloch 团队、Harvard 大学 Greiner 团队的系列突破）。

• 可调控性：在冷原子系统内，通过 Feshbach 共振可调节强关联斥力 $U_d$，通过空间激光调控可精确重构铜氧平面的晶格结构与杂化强度。
• 实验指南：本研究预测的 ZSA 相图以及 $T_c^{\text{max}}$ 的二维物理等高线（图12），为冷原子实验寻找最优高温超导区指明了明确的靶点。由于冷原子系统目前同样面临极低温度难以实现的困境，本项工作证明的“利用母体正常态物理量（$\Delta_I$、 $2p_p$）直接预测最佳 $T_c$”的普适规律，免去了冷原子实验必须先达到超低温并精确测量超导配对的烦恼，具有无与伦比的工程指导意义。