做超导的第一性原理计算最难受的时刻往往不是任务排不进队列而是几十万核时跑完算出来的晶格常数、声子频率、电声耦合常数跟实验差了一截最后套 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.1Tc 能从 23 K 涨到 27 K。这还只是结构误差一条线还没算电子结构本身在费米面附近的偏差。所以“泛函选型”在超导计算里不是一个保守与激进的风格问题而是一个直接决定结论量级的技术决策。你做结构预测、做相稳定性、做声子软模判断用的泛函必须能把晶格常数和力常数压到足够准的范围内后面所有环节才有意义。1.2 SCAN 在 Jacobs 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 几种常见泛函在超导相关量上的对照下面这张表是我按自己用过的体系整理的定性对照不是普适定律但方向性可以参考。数值区间来自我跑过的过渡金属间化合物和轻元素体系换体系会变。泛函晶格常数偏差相对实验声子频率系统性偏差费米面处态密度单步计算成本相对 PBE适合的场景LDA-1% 到 -2%低估偏高 3%~8%常高估1.0高压、致密体系历史数据对标PBE0.5% 到 2%高估偏低 2%~6%常略高估1.0大体系筛选、粗扫相图PBEsol0.2% 到 0.8%接近 PBE略好与 PBE 接近1.0固体晶格常数性价比高SCAN-0.3% 到 0.5%偏差通常 3%方向不固定与实验吻合较好2~3精算晶格、声子、电声耦合r2SCAN与 SCAN 相当个别体系更好与 SCAN 相当与 SCAN 相当1.5~2SCAN 收敛困难时的替代SCANrVV10层状体系显著改善层间模式更可靠影响较小3~4层状、含范德华间隙的超导体系这张表里最值得注意的一列是“声子频率系统性偏差”。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λ) / (λ - μ*(10.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。分母 λ - μ*(10.62λ) 1.0 - 0.1×1.62 0.838。指数项 -1.04×2/0.838 -2.482exp(-2.482) 0.0835。Λ₁ 2.46×1.38 3.395f₁ [1(1/3.395)^1.5]^{1/3} 1.051。Λ₂ 1.82×1.63×1.024 3.038f₂ 1 (1×0.6708)/(19.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.344exp 0.0959。f₁ 变成 1.058f₂ 变成 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 问题速查表下面这张表是我自己攒的按出现频率排序。左边是症状右边是优先排查的方向。症状最可能的原因优先动作SCF 不收敛能量振荡混合系数过大、τ 噪声混合系数降到 0.1~0.2FFT 网格加密应力波动大优化不收敛截断能不足、FFT 网格不足截断能提到 800 eV 以上网格 2 倍声子谱出现虚频但实验无相变超胞太小、位移幅度过大或过小超胞加倍位移改 0.03 Å正负双向费米能级附近能级排序随参数变k 网格和展宽不匹配先展宽 0.2 eV 稳定再窄展宽细算磁矩在迭代中跳变多个近简并磁态多初态扫描比较总能量λ 随 q 网格持续漂移q 网格不足至少 12×12×12用 PBE 先定网格Tc 与实验差一倍以上μ* 取值、ω_log 误差做 μ* 敏感性分析查低频模式同一结构两次算的总能量差几十 meV未完全收敛或初始 τ 不同用上一步的电荷和 τ 重启4.2 那个“看起来收敛了其实没收敛”的陷阱这是我在 SCAN 上栽得最狠的一次。当时算一个含轻重元素混合的金属间化合物SCF 在 60 步内收敛能量变化小于 1 meV力小于 0.02 eV/Å一切看起来完美。声子谱算出来有几支低频光学支出现软化我以为是真实的软模写了一堆分析。后来一位同行建议我把 FFT 网格加密一倍重算结果那几支软模消失了整个物理图像推翻重来。问题出在 τ 的网格表示上。电荷密度在实空间的分布相对平滑FFT 网格对它的采样余量很大τ 因为含波函数梯度的平方高频成分多网格不足时相当于给 τ 做了一次低通滤波交换关联势里的部分细节被抹掉。抹掉之后能量反而“更平滑”了SCF 收敛更快看上去更稳实际上物理已经错了。这个错误特别隐蔽因为所有常规的收敛判据能量、力、应力都可能通过。我的应对办法现在是固定的任何 SCAN 计算的最终结果都要做一次 FFT 网格双倍加密的复核只看关键量费米能、N(E_F)、最大力、声子频率的前几支如果变化超出阈值就说明原网格不够。这个复核的成本大概是主计算的 30%但能挡住大部分这类问题。4.3 成本控制什么时候该退到 r2SCAN 或者干脆用 PBESCAN 贵而且贵得不是线性的。当你的体系超过 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 在轻元素体系上的晶格常数改善是它最有说服力的地方之一。第三类是层状或者有范德华间隙的体系。不过这类要配 SCANrVV10光用 SCAN 不够因为半局域泛函天生描述不了长程关联。加了非局域项之后层间距和层间声子模式低频的褶皱模式会明显改善而这些模式在层状超导体的电声耦合里权重不小。第四类是高压下的致密体系。这一块历史上是 LDA 的地盘因为 LDA 低估体积的毛病在高压下反而“正好”补偿了。但 SCAN 在高压区的表现通常也稳而且对压力诱导的结构相变判断更可靠。如果你的研究涉及压力调控超导SCAN 值得一试。5.2 边界与失效场景SCAN 也不是万能药说完好处得说清楚它救不了什么否则容易误导。强关联体系是第一个边界。SCAN 毕竟还是半局域泛函对莫特绝缘体、近藤效应、重费米子这些本质上需要非局域自能描述的现象它无能为力。在重费米子超导体里f 电子的局域-巡游竞争是关键SCAN 给出的 f 态位置可能比 PBE 好一些但离定量描述还差得远。这类体系必须上 DFTDMFT 或者动力学平均场这类方法。自旋涨落驱动的非常规超导是第二个边界。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 都给出同样定性结论的部分你可以放心往前推凡是两者分歧的部分就是需要额外小心、需要实验或者更高等级方法仲裁的地方。这种分歧本身往往比任何一个单独的结果更有信息量。