做超导的第一性原理计算，最难受的时刻往往不是任务排不进队列，而是几十万核时跑完，算出来的晶格常数、声子频率、电声耦合常数跟实验差了一截，最后套 McMillan 公式得到的 Tc 和实测值对不上，你还不知道是理论框架的问题还是泛函的问题。这几年我把密度泛函理论（DFT）用在超导相关体系上，越来越强烈的一个感受是：泛函选型对结果的支配程度，远超很多教程里“用 PBE 就行”那句话给人的印象。SCAN 泛函（Strongly Constrained and Appropriately Normed meta-GGA）之所以值得单独拎出来写一篇，是因为它在晶格常数、声子频率、金属体系费米面附近的电子结构这几件直接影响超导性判断的事情上，确实给出了与 PBE 明显不同的答案，而且这些差异不是小数点后第三位的抖动，是能把“这材料会不会超导、Tc 大概多少”的定性结论改掉的那种差异。

这篇东西写给谁看：手上正在用 DFT 算电子结构、想往电声耦合和超导方向延伸的硕博生；已经习惯 PBE 但被审稿人问过“为什么不用 meta-GGA”的同行；以及想搞清楚 SCAN 到底改善了哪部分物理、又在什么情况下会把计算成本翻三倍却换不来收益的实战派。我不打算把它写成泛函发展史的综述，而是按“为什么会偏、SCAN 改了什么、怎么算、坑在哪、结果怎么读”这条线走一遍，中间穿插可以直接抄的参数和收敛测试做法。

1. 超导计算里泛函选型这件事，为什么值得单独拎出来讲

1.1 从晶格常数到 Tc，误差是被逐级放大的

很多人对泛函误差的容忍度，来自能带计算的经验：带隙差个零点几电子伏，趋势还在。但超导这条链不一样。DFT 给出来的第一手量是总能量和力，由此得到平衡晶格常数和力常数；力常数进声子谱，声子谱进电声耦合矩阵元；耦合矩阵元进 Eliashberg 谱函数 α²F(ω)，最后进 Tc 公式。每一环都有误差，而且很多环节是乘性放大而不是加性叠加。

举个我自己算过的场景：PBE 把某类六角结构的 a 轴高估了约 0.8%，c 轴高估 1.2%。听起来不大，但声子频率对体积的敏感度在软模附近极高，尤其是低频声学支和部分光学支，频率偏移 5% 很常见。电声耦合常数 λ 的量纲关系大致是 λ ∝ α²F(ω)/ω 的积分，低频段权重被 1/ω 放大，声子频率系统性低估会同时把 α²F 的低频部分推高、把分母压小，两个效应叠加，λ 容易被高估 10% 以上。而 Tc 对 λ 是指数敏感的——我在第 3.5 节会用一个具体数字演示，λ 从 1.0 涨到 1.1，Tc 能从 23 K 涨到 27 K。这还只是结构误差一条线，还没算电子结构本身在费米面附近的偏差。

所以“泛函选型”在超导计算里不是一个保守与激进的风格问题，而是一个直接决定结论量级的技术决策。你做结构预测、做相稳定性、做声子软模判断，用的泛函必须能把晶格常数和力常数压到足够准的范围内，后面所有环节才有意义。

1.2 SCAN 在 Jacob's ladder 上的位置，以及它凭什么被认为“更准”

按 Perdew 那个著名的“雅各布天梯”分类，LDA 只看局域电子密度，GGA 加上密度梯度，meta-GGA 再加上动能密度 τ（或者等价地，占据轨道信息）。SCAN 属于第三级，2015 年提出，核心卖点是它是第一个（在构造时）满足全部已知精确约束的半局域泛函。

这里的“约束”不是数学装饰品，每一条都对应一类真实物理极限。比如均匀电子气极限、单电子体系的自相互作用消除、线性的密度响应、以及那个著名的 Lieb-Oxford 界。以前 GGA 时代，你满足这条就得牺牲那条，PBEsol 为了固体晶格常数优化了梯度展开系数，代价是对分子体系和非均匀区域的表现变差。SCAN 的思路是把这些极限情形作为插值端点，用动能密度构造一个判据量 α 来判断“当前这个点更接近单轨道还是更接近均匀电子气”，然后在这两个端点之间平滑插值。

实际效果层面，我个人的观察是：SCAN 对晶格常数的改善是稳健的，尤其是过渡金属和主族化合物的体积，通常能把 PBE 的 1%~2% 高估压到 0.5% 以内，个别体系甚至反过来略微低估。对金属体系，SCAN 给出的费米面附近态密度和能带色散与实验（ARPES、de Haas-van Alphen）的吻合度在多数体系里优于 PBE。这两点恰好是超导计算最吃紧的输入。

1.3 几种常见泛函在超导相关量上的对照

下面这张表是我按自己用过的体系整理的定性对照，不是普适定律，但方向性可以参考。数值区间来自我跑过的过渡金属间化合物和轻元素体系，换体系会变。

这张表里最值得注意的一列是“声子频率系统性偏差”。LDA 因为低估体积，声子频率系统性偏高，这在历史上反而让一些体系的 Tc 算得“碰巧”接近实验。PBE 反过来，体积大了，声子频率偏低，λ 偏高。这种“错得刚好抵消”的现象在文献里不罕见，也是为什么不能只看 Tc 的最终数值对不对，得看中间量。

提示：如果你接手一个老项目，发现前人用 PBE 算的 Tc 和实验吻合得很好，先别急着换 SCAN。把中间量（晶格常数、ω_log、λ）拉出来单独对比实验，确认那个“吻合”是不是误差抵消的结果。换成 SCAN 之后 Tc 反而偏离，这不一定是 SCAN 错了，很可能是之前抵消掉的两个误差被拆开了。

2. SCAN 到底改了什么：从约束条件讲到数值代价

2.1 从 LDA 到 GGA 再到 meta-GGA 的递进逻辑

要讲清楚 SCAN 的机制，得先接受一个前提：半局域泛函的表达能力受限于它“能看到多少信息”。LDA 只看到某一点的电子密度 n(r)，所以它对密度不均匀的区域天然无能为力——一块密度梯度剧烈的区域和一块均匀区域，只要密度值相同，LDA 给出的交换关联能一模一样，这显然不符合物理。

GGA 把 ∇n 加进来，于是能区分“密度高但平缓”和“密度高且陡变”。这一级解决了原子、分子的大部分问题，但引入了新的麻烦：梯度展开在高密度梯度区域会发散，需要截断，截断方式就是各种 GGA 的分歧所在（PBE、PBEsol、BLYP 的差别很大程度上就是这个截断函数不同）。更根本的问题是，GGA 无法区分“真正的单电子区域”和“多电子但密度恰好比较均匀的区域”，因为这两者的 n 和 ∇n 可以很像。

meta-GGA 加进来的 τ = (1/2)Σ|∇ψ_i|²，也就是动能密度。这个东西的价值在于它隐含了轨道占据信息：单个占据轨道的区域，τ 和 n 之间有严格的等号关系（von Weizsäcker 关系）；而在均匀电子气里，τ 是另一个确定值。于是你有了一个天然的判据来区分这两种极端情形。

2.2 α 插值：把两个极端端点连起来的那根绳子

SCAN 里最关键的一个量是

α = (τ - τ^W) / τ^unif

其中 τ^W 是 von Weizsäcker 动能密度，τ^unif 是同样密度下均匀电子气的托马斯-费米动能密度。α 的取值有明确边界：单轨道区域（如氢原子尾部、分子解离的渐近区）α = 0；均匀电子气 α = 1；一般的固体内部 α 在 0.5 到 1 之间浮动。

SCAN 的构造方式，就是在 α = 0 和 α = 1 这两个端点分别构造满足该极限下所有已知约束的交换增强因子 h_x^1(s) 和 h_x^0(s)，然后用一个 α 的函数 f_x(α) 做插值。这个插值函数不是随便选的，它必须保证插值结果在整个 α ∈ [0,1] 区间内都不违反约束，还得保证导数连续。这就是“appropriately normed”那半句话的含义——不只是满足约束，还得在约束允许的整个范围内都保持数值上的合理标定。

这个构造的直接后果是：SCAN 对自相互作用误差的修正比 PBE 好得多。在 PBE 里，一个被过度离域化的电子会把能量压得过低，导致分子解离极限错、过渡态能垒低估、金属氧化物带隙偏小。SCAN 因为 α 判据的存在，在单电子区域会自动切换到精确的 h_x^1 形式，把自相互作用大部分消化掉。对超导计算来说，这件事的间接价值在于：费米面附近的 d 电子或 f 电子局域化程度被描述得更接近真实，态密度形状更可信。

2.3 满足全部约束的代价：数值噪声和收敛难度

天下没有免费的午餐。α 依赖 τ，而 τ 是轨道梯度的平方和，它是一个比电荷密度“粗糙”得多的量。电荷密度是平滑的，在原子核附近发散但整体光滑；τ 在节点附近、在占据轨道切换处、在 FFT 网格上做傅里叶变换时，会带出大量高频分量。

实际表现就是：SCAN 对 FFT 网格密度和平面波截断能的要求比 PBE 高出一截，而且不是线性的。我做过一组测试，同一个结构，PBE 在截断能 500 eV、FFT 网格默认时应力收敛到 0.1 kBar 以内；SCAN 在同样参数下应力抖动能有几个 kBar，把截断能提到 800 eV、FFT 网格加密到 2 倍，抖动才降到可接受范围。这就直接吃掉了它“只比 PBE 贵 2 到 3 倍”的成本优势——算上必须加密的网格，端到端时间经常是 4 到 6 倍。

另一个代价是自洽场（SCF）收敛。因为交换关联势对 τ 的依赖是非线性的，SCF 迭代过程中电荷密度和动能密度会互相“追尾”，出现振荡。默认的混合方案在 SCAN 下经常不收敛，或者收敛到错误的局部极小。这个坑我在第 4 节会专门展开，先给个结论：混合系数得压到 0.1 到 0.2 之间，电子步上限提到 200 以上，最好开启对密度和 τ 的双重预条件。

3. 从结构到 Tc：一条完整计算链的实操拆解

3.1 收敛测试：截断能、FFT 网格、k 点该测到什么程度

这是最容易被跳过、又最影响后面所有结果的一步。我的做法是分三层测，不要一次全测。

第一层测截断能，固定一个中等 k 网格（比如 8×8×8 或按倒空间密度约 0.025 Å⁻¹ 折算），扫描截断能从 400 到 1000 eV，步长 100 eV，看总能量、最大力、以及应力张量的对角元。判据不是总能量变化小于 1 meV/atom（这个判据在 SCAN 下太松，能量收敛比力收敛快得多），而是最大残余力变化小于 5 meV/Å 且应力对角元变化小于 0.5 kBar。通常 700 到 900 eV 才能满足。

第二层测 FFT 网格。这一层在 SCAN 下不能省。做法是固定截断能，把网格密度从默认的 1 倍逐步加密到 1.5 倍、2 倍，观察费米能级附近本征值的变化。如果最高占据态和最低空态的能量差在网格加密后还能变动超过 10 meV，说明网格不够。我遇到过的极端情况是一个含重元素的金属间化合物，默认网格下费米能级附近的能级排序都是错的。

第三层测 k 点。金属体系必须用展宽（Methfessel-Paxton 或 Fermi smearing），展宽宽度取 0.1 到 0.2 eV 起步。k 网格密度的判据是：把网格加密一档，费米能级处态密度 N(E_F) 的变化小于 2%，同时总能量变化小于 1 meV/atom。N(E_F) 这个判据在超导计算里比总能量更重要，因为 λ 和 γ（Sommerfeld 系数）都直接正比于它。

注意：展宽宽度和 k 网格是耦合的。展宽太宽会把费米面附近的细节抹平，让你“看起来收敛得很快”；展宽太窄则 SCF 震荡。我的经验是先用一个较宽的展宽（0.2 eV）让 SCF 稳定收敛，再用较窄的（0.05 到 0.1 eV）跑最终静态，同时把 k 网格加密一档。两次的 N(E_F) 差在 3% 以内才算过关。

这里可以给一个参数计算的例子。假设你的倒空间想达到 0.02 Å⁻¹ 的密度，晶格常数 a = 5.0 Å，那么每个方向的 k 点数大约是 2π/(a × 0.02) ≈ 2π/(0.1) ≈ 62.8，取整并按对称性约化前大概是 63 个点，实际用 12×12×12 或者更密的网格。这个折算方法在六方、四方体系里要做倒格子度规修正，不能直接用实空间长度。

3.2 结构弛豫与初始磁态：SCAN 的磁性判断更“固执”

结构弛豫本身不复杂，SCAN 下要注意的是收敛阈值和磁态。

阈值方面，力收敛判据我一般设 0.01 eV/Å，应力判据 0.5 kBar。SCAN 的势能面比 PBE“尖锐”一些，同样的步长容易出现能量震荡，所以离子步的步长建议压小，或者用带信赖域的优化器而不是纯梯度下降。BFGS 类优化器在 SCAN 下表现比共轭梯度稳。

磁态是更大的坑。SCAN 对自旋极化的描述比 PBE 更“倾向明确”：如果一个体系在 PBE 下是弱磁性、磁矩 0.2 μ_B 左右的可疑状态，SCAN 经常要么把它推成明显的磁性（0.8 μ_B 以上），要么压成完全非磁。这个特性在判断超导母体的磁基态时非常关键——铁基、重费米子这类体系，磁性态和超导态是竞争关系，基态选错，后面的声子和费米面全错。

我的操作习惯是：对每个可疑体系，至少跑三个初始磁构型——非磁、铁磁（初始磁矩给 0.5 到 1.0 μ_B）、以及至少一种反铁磁排布。三个都收敛后比总能量。不要只跑一个然后看磁矩大小下结论。SCAN 下从非磁初态出发有时会“卡”在非磁解上，即使磁态能量更低，所以非磁解收敛后最好再人为施加一个小的初始磁矩重启一次，确认它是否真的不稳定。

3.3 电子结构：能带、态密度、费米面的计算细节

静态计算拿到电荷密度后，能带和态密度是常规操作，但 SCAN 下有两个细节值得说。

第一，能带计算不需要重新做 SCF，用收敛好的电荷密度和 τ 做非自洽即可，但要注意有些代码在非自洽步骤里对 τ 的处理和自洽步骤不完全一致，会出现能带整体平移的假象。判断方法是检查非自洽计算的费米能与自洽最后一步的费米能是否一致，差超过 20 meV 就要警惕。

第二，态密度计算用的 k 网格要比 SCF 密得多，尤其是金属。因为 N(E_F) 是要拿去算 λ 和 γ 的，误差直接传递。我一般把 SCF 用的网格加密 2 到 3 倍做 DOS，配合四面体方法而不是展宽。四面体方法在 SCAN 下也有个坑：τ 的网格噪声会让四面体插值在个别小四面体上给出非物理的负态密度，看到负值不要慌，先检查是不是 FFT 网格不够密。

费米面本身在超导机制分析里价值很大，尤其是判断 nesting 和可能的 nesting 驱动的不稳定性。SCAN 给出的费米面形貌和 PBE 在大拓扑上通常一致，但在小的 pocket 上会有差别——有的 pocket 在 PBE 下存在，SCAN 下消失。这种差别对判断“是不是某个特定的费米面片负责超导”很关键，值得花时间做。

3.4 声子与电声耦合：SCAN 下做 DFPT 的真实体验

这一段是整个链条里最贵的部分，也是 SCAN 的优势最能体现、同时坑最多的地方。

两条路线：密度泛函微扰理论（DFPT）和有限位移超胞法。DFPT 的优势是不需要超胞，q 网格可以取得比较密，对金属体系的电声耦合计算几乎是必需——因为 λ 需要对费米面上的大量 q 点做积分，超胞法在金属里收敛极慢。代价是 DFPT 对泛函的支持有限，meta-GGA 的 DFPT 实现比 GGA 复杂得多，因为要处理 τ 对波函数的一阶导数，非局域项的处理也更麻烦。目前主流平面波代码对 SCAN 的 DFPT 支持程度不一，有的只支持 Γ 点声子，有的支持全 q 网格但要求特定的赝势类型。上手之前务必先查清你手上代码的版本支持到什么程度，别写到一半发现算不了。

有限位移法在 SCAN 下反而是更稳妥的选择，尤其是做小体系的完整声子谱验证。做法是建超胞（2×2×2 或 3×3×3，取决于声学支的色散范围），逐个位移不等价原子，算力，提取力常数矩阵，再做傅里叶插值。SCAN 下超胞法的坑在于：位移幅度不能太小。PBE 下常用 0.01 Å 的位移，SCAN 因为力上的数值噪声更大，0.01 Å 的位移算出来的力差可能被噪声淹没。我一般用 0.03 Å，并且对正负位移都算，取中心差分，噪音能压下去一半。

电声耦合常数 λ 的 DFPT 计算，输出的是每个 q 点上的 α²F 贡献。这里必须强调 q 网格收敛：λ 对 q 网格的收敛速度比总能量慢得多，金属体系经常需要 12×12×12 甚至 16×16×16 的 q 网格才能收敛到 5% 以内。这个成本在 SCAN 下是实打实的几天到几周。有个实用的降本技巧：先用 PBE 扫一遍 q 网格收敛性，确定需要多密的网格，再把这个网格直接用到 SCAN 上，因为 q 空间的收敛行为和泛函关系不大，主要取决于费米面的形貌和声子色散的各向异性。

还有一个必须提的点：库仑赝势 μ*。它不出现在 DFT 里，是 Eliashberg 方程里的一个参数，物理上代表电子间的库仑排斥被声子延迟效应削弱后的残余。常见取值 0.1 到 0.15，很多计算直接取 0.13 了事。这个参数的选取对 Tc 的影响不比 λ 小，尤其是在 λ 较小的体系（λ < 0.6）里，μ* 从 0.1 变到 0.15 能让 Tc 掉一半。严谨一点的做法是做约束随机相位近似（cRPA）算 μ*，成本很高但值得在关键体系上做一次。省钱的做法是至少做敏感性分析：μ* 取 0.1、0.13、0.15 各算一次 Tc，把区间报出来，而不是只报一个数。

3.5 从 α²F(ω) 到 Tc：公式、参数和误差传递的量化演示

拿到 α²F(ω) 之后，走 Allen-Dynes 修正的 McMillan 公式：

Tc = (f₁ f₂ ω_log / 1.2) × exp[ -1.04(1+λ) / (λ - μ*(1+0.62λ)) ]

其中 ω_log 是对数平均声子频率，f₁ 和 f₂ 是两个修正因子，形式分别是

f₁ = [1 + (λ/Λ₁)^{3/2}]^{1/3}，Λ₁ = 2.46(1 + 3.8μ*) f₂ = 1 + λ²(1 - λ/Λ₂) / (λ² + Λ₂²)，Λ₂ = 1.82(1 + 6.3μ*)(1 + 0.024λ)

我用一组具体数字走一遍，方便你对照自己的结果。设 λ = 1.0，ω_log = 300 K，μ* = 0.1。分母 λ - μ*(1+0.62λ) = 1.0 - 0.1×1.62 = 0.838。指数项 = -1.04×2/0.838 = -2.482，exp(-2.482) = 0.0835。Λ₁ = 2.46×1.38 = 3.395，f₁ = [1+(1/3.395)^1.5]^{1/3} = 1.051。Λ₂ = 1.82×1.63×1.024 = 3.038，f₂ = 1 + (1×0.6708)/(1+9.229) = 1.066。最终 Tc = 1.051×1.066×(300/1.2)×0.0835 ≈ 23.4 K。

现在把 λ 抬高 10%，到 1.1，其他不动。分母变成 1.1 - 0.1×1.682 = 0.9318，指数项 -1.04×2.1/0.9318 = -2.344，exp = 0.0959。f₁ 变成 1.058，f₂ 变成 1.074。Tc = 1.058×1.074×(300/1.2)×0.0959 ≈ 27.2 K。λ 变 10%，Tc 变 16%。如果 λ 本身在 1.3 附近，同样的 10% 相对误差会把 Tc 从 35 K 推到 45 K 以上。

这就是为什么前面那三层收敛测试不是形式主义。声子频率低估 5%，通过 λ ∝ 1/ω 的关系（粗略地，α²F 的幅度也随频率变），λ 可能高估 10%，Tc 就偏 16%。SCAN 的价值在于它把结构误差从 PBE 的 1% 以上压到 0.5% 以内，这个改进直接对应声子频率的 2% 到 3% 改善，落到 Tc 上就是 5% 到 10% 的精度提升。听起来不多，但在区分“这材料是 15 K 还是 25 K”的时候，够了。

提示：Tc 计算里最容易被忽略的误差源不是 λ，而是 ω_log。因为 ω_log 是对数平均，低频声子模式对它的拖累极大，而这些低频模式恰恰是最难算准的。如果你的 α²F 在 50 cm⁻¹ 以下还有明显权重，务必检查那部分是不是数值噪声。

4. 踩坑实录：SCAN 跑超导体系的高频问题与排查路径

4.1 问题速查表

下面这张表是我自己攒的，按出现频率排序。左边是症状，右边是优先排查的方向。

4.2 那个“看起来收敛了其实没收敛”的陷阱

这是我在 SCAN 上栽得最狠的一次。当时算一个含轻重元素混合的金属间化合物，SCF 在 60 步内收敛，能量变化小于 1 meV，力小于 0.02 eV/Å，一切看起来完美。声子谱算出来有几支低频光学支出现软化，我以为是真实的软模，写了一堆分析。后来一位同行建议我把 FFT 网格加密一倍重算，结果那几支软模消失了，整个物理图像推翻重来。

问题出在 τ 的网格表示上。电荷密度在实空间的分布相对平滑，FFT 网格对它的采样余量很大；τ 因为含波函数梯度的平方，高频成分多，网格不足时相当于给 τ 做了一次低通滤波，交换关联势里的部分细节被抹掉。抹掉之后能量反而“更平滑”了，SCF 收敛更快，看上去更稳，实际上物理已经错了。这个错误特别隐蔽，因为所有常规的收敛判据（能量、力、应力）都可能通过。

我的应对办法现在是固定的：任何 SCAN 计算的最终结果，都要做一次 FFT 网格双倍加密的复核，只看关键量（费米能、N(E_F)、最大力、声子频率的前几支），如果变化超出阈值就说明原网格不够。这个复核的成本大概是主计算的 30%，但能挡住大部分这类问题。

4.3 成本控制：什么时候该退到 r2SCAN 或者干脆用 PBE

SCAN 贵，而且贵得不是线性的。当你的体系超过 100 个原子、又要做声子或电声耦合的时候，SCAN 的端到端成本可能是一个月以上的机时。这时候要有一个明确的降级策略，而不是硬扛。

r2SCAN 是我目前最常用的折中。它是 SCAN 的正则化版本，把 α 依赖里那些导致数值不稳定的项做了平滑处理，代价是牺牲掉一部分约束的严格满足。实测下来，r2SCAN 的数值稳定性明显更好，同样的截断能和 FFT 网格下应力和力抖动小一个量级，SCF 也更容易收敛。精度上，晶格常数和声子频率与 SCAN 的差异通常在 0.3% 以内，对 Tc 的影响在 2% 以内。如果你的体系规模大到 SCAN 跑不动，r2SCAN 是首选替代。

什么情况下该彻底退回 PBE？我的判断标准是：如果你的结论只关心趋势（比如一系列同构化合物的 Tc 排序），而且这个趋势在 PBE 和 SCAN 下一致，那么用 PBE 做大规模筛选、只在关键节点用 SCAN 复核，是合理的。但如果你的结论是“某材料在常压下是否超导”这种二元判断，而且 PBE 和 SCAN 给出了相反的结论，那么必须用 SCAN，而且要配合实验数据（如果存在）或者更高等级的方法（如扩散蒙特卡洛，虽然对固体很贵）来仲裁。

还有一个隐形成本容易被忽略：SCAN 的赝势兼容性。不是所有为 PBE 优化的赝势在 SCAN 下都好用，尤其是那些半芯态处理得比较激进的赝势。SCAN 对半芯态的敏感度更高，因为 τ 在芯区附近变化剧烈。建议优先选专门为 meta-GGA 测试过的赝势库，或者至少用全电子方法做一次小体系的交叉验证。

5. 结果怎么读：SCAN 能扭转什么，救不了什么

5.1 真正从 SCAN 受益的体系类型

按我的经验，SCAN 带来的改进不是均匀分布的，有几类体系受益特别明显。

第一类是含 3d 过渡金属的化合物和合金。PBE 对这类体系的晶格常数高估、对 d 带宽度的描述偏差，会直接影响费米面附近态密度的形状。SCAN 的自相互作用修正让 d 电子的局域化程度更接近真实，态密度的尖峰位置和高度都更可信。这对判断 Stoner 不稳定性、自旋涨落强度很关键。

第二类是轻元素体系，特别是含氢、锂、硼的化合物。这类体系的声子频率高，对结构误差的绝对敏感度大。PBE 高估晶格常数 1% 在这种轻质体系里能让最高光学支频率偏移 5% 以上，直接改变 ω_log 和 Tc。SCAN 在轻元素体系上的晶格常数改善是它最有说服力的地方之一。

第三类是层状或者有范德华间隙的体系。不过这类要配 SCAN+rVV10，光用 SCAN 不够，因为半局域泛函天生描述不了长程关联。加了非局域项之后，层间距和层间声子模式（低频的褶皱模式）会明显改善，而这些模式在层状超导体的电声耦合里权重不小。

第四类是高压下的致密体系。这一块历史上是 LDA 的地盘，因为 LDA 低估体积的毛病在高压下反而“正好”补偿了。但 SCAN 在高压区的表现通常也稳，而且对压力诱导的结构相变判断更可靠。如果你的研究涉及压力调控超导，SCAN 值得一试。

5.2 边界与失效场景：SCAN 也不是万能药

说完好处，得说清楚它救不了什么，否则容易误导。

强关联体系是第一个边界。SCAN 毕竟还是半局域泛函，对莫特绝缘体、近藤效应、重费米子这些本质上需要非局域自能描述的现象，它无能为力。在重费米子超导体里，f 电子的局域-巡游竞争是关键，SCAN 给出的 f 态位置可能比 PBE 好一些，但离定量描述还差得远。这类体系必须上 DFT+DMFT 或者动力学平均场这类方法。

自旋涨落驱动的非常规超导是第二个边界。SCAN 能给你更准的费米面和磁化率，但超导配对本身的对称性分析和能隙函数计算，超出了单粒子图像的能力范围。你可以用 SCAN 的输出作为输入，去做自旋涨落配对计算（比如基于随机相位近似或 fLEX），但不要指望 SCAN 本身告诉你 d 波还是 s± 波。

非绝热效应是第三个边界。当 λ 接近或超过 1，或者声子频率和电子能标接近时，Migdal-Eliashberg 理论的绝热近似开始失效。这类体系（氢化物是典型）需要非绝热修正，SCAN 提供的准粒子能带只是起点，还得额外处理 vertex 修正和电子-声子相互作用的非绝热项。这个方向上的定量计算目前还是活跃的研究领域，泛函选型只是其中一环。

第四个边界是 μ* 和掺杂效应。SCAN 改善的是基态电子结构和声子，不改变 μ* 需要额外计算这个事实。另外，如果你算的是掺杂或者无序体系，SCAN 对单个缺陷或者单个构型的精度提升，会被构型平均的不确定性淹没，这时候提高泛函精度的边际收益很低，不如把机时花在采样上。

6. 我个人在这条链上的一点使用习惯

最后分享几个我这两年固定下来的操作习惯，都是踩坑换来的。

第一，任何 SCAN 计算的元数据里，我一定记下三样东西：截断能、FFT 网格倍数、展宽宽度。这三样如果不一致，两次计算结果就没法直接比较，哪怕其他参数全一样。我见过不止一次因为换了计算节点、代码读到的默认网格变了，导致结果对不上，排查了两天才发现是这个原因。

第二，声子计算之前，先做一次“静水压测试”。做法是给结构施加正负 1% 的体积变化，看总能量曲线和力的响应是否平滑。如果 SCAN 在某个体积点上出现了明显的能量跳变或者力异常，说明那个点的 SCF 收敛到了不同的解，声子计算肯定出问题。这个测试成本很低，能提前挡掉一大类问题。

第三，Tc 永远报区间不报单值。λ 的误差、ω_log 的误差、μ* 的选取，三个一叠加，报出来的单值精度是假的。我现在写论文或者内部报告，都是给出 μ* 取 0.1 和 0.15 两个端点对应的 Tc，再加上网格收敛带来的区间。审稿人看到这种诚实的误差处理，反而更容易接受。

第四，如果条件允许，用两种泛函并行算关键体系。不是为了让它们互相验证对错（它们可能都错），而是为了看清楚“哪些结论是泛函无关的”。凡是 PBE 和 SCAN 都给出同样定性结论的部分，你可以放心往前推；凡是两者分歧的部分，就是需要额外小心、需要实验或者更高等级方法仲裁的地方。这种分歧本身，往往比任何一个单独的结果更有信息量。