4.2 赝势与 PAW 方法

4.2 赝势与 PAW 方法

本节摘要：直接对全电子做平面波计算，原子核附近的波函数振荡会让基组需求爆炸。赝势与 PAW 用「冻结芯电子」的思路，把核与芯电子的复杂势替换成平滑的有效势，只显式处理价电子。本节讲清赝势的构造原理、范数守恒与超软赝势的区别、PAW 的线性变换思想，以及势库选择与验证方法。

本节你将学会

阅读完本节，你应当能够：

1. 解释为什么全电子平面波计算在核附近成本爆炸
2. 复述赝势的核心思想：冻结芯电子、用平滑势替代、保持散射性质
3. 对比范数守恒赝势与超软赝势的取舍（精度 vs 截断能）
4. 说出 PAW 方法与赝势的关系：线性变换重构全电子波函数
5. 知道如何选择与验证势库（VASP PAW、QE 赝势库）
6. 理解鬼态问题的来源与防范

问题与直觉

想象一个钨原子，74 个电子。其中内层电子紧紧贴着核，波函数在核附近剧烈振荡，能量很低。如果要对它做全电子平面波计算，就得用极其多的平面波去描述那些振荡——计算成本直接爆炸。

但化学直觉告诉我们：内层电子几乎不参与成键。钨的化学性质由价电子决定，芯电子只是「蹲在核附近」的背景。那为什么不让芯电子继续蹲着，只显式处理价电子？

这就是赝势的核心思想：把原子核 + 芯电子当成一个「离子实」，用一个精心构造的平滑有效势描述它对价电子的作用。价电子感受到的不再是又深又尖的库仑井，而是一个平滑的「赝势」。伪波函数在芯区不再剧烈振荡，平面波基组就能用较少的项描述它。

代价呢？赝势是近似，丢失了芯电子的信息。对大多数化学性质（成键、能带、磁性）无伤大雅，但想精确算核附近的性质（NMR、XAS 芯能级）就不够了——这时 PAW 或全电子方法才是正解。

值得强调：这个「只算价电子」的取舍不是偷懒，而是有明确的物理依据——化学键、能带、磁性都由价电子的重排驱动，芯电子在这些过程中几乎保持「冻结」。赝势只是把这个物理事实编码成了计算参数。理解了这条依据，你就能预判赝势在哪些场景会失效：一切需要「芯电子参与」的性质，比如核磁共振屏蔽、芯能级谱、极端压力下的芯态重叠。

核心原理

2.1 赝势的构造原理

一个好赝势必须满足几个条件：

1. 散射性质匹配：在给定能量下，赝势在芯区的散射相移必须与真实势一致——这是保证「价电子行为相同」的根本
2. 正交性：伪波函数要与被冻结的芯轨道正交，避免占据芯态
3. 可转移性：同一个元素的赝势，在不同化学环境（分子、固体、不同配位）都要准确——不要求每个环境重新生成
4. 平滑性：伪波函数在芯区平滑，截断半径 r_c 内不振荡

构造流程：先对自由原子做全电子计算，得到精确的芯轨道与价轨道；然后在截止半径 r_c 外要求伪波函数与全电子波函数一致，r_c 内平滑化；最后反推赝势。

伪波函数：r < r_c 平滑，r ≥ r_c 与全电子一致
赝势：    r ≥ r_c 与全电子势一致，r < r_c 平滑修正
​

这里「r_c」的选择本身是一门平衡艺术：r_c 越小，赝势越接近全电子、精度越高，但伪波函数越「硬」、需要的平面波越多；r_c 越大，赝势越软、计算越省，但丢失的真实信息越多。每个势文件背后都藏着这样一个「精度—效率」的抉择，理解它你就能看懂为什么同一个元素在不同软件里势文件长得不一样。

2.2 范数守恒 vs 超软赝势

范数守恒赝势（NCP）：要求 r_c 内伪波函数的电荷积分与全电子一致：

∫₀^r_c |ψ_pseudo|² r² dr = ∫₀^r_c |ψ_all|² r² dr
​

这个「范数守恒」保证了散射性质与可转移性都好，是赝势发展的里程碑。HSC、Troullier-Martins 是代表。缺点是：对 d/f 电子局域性强的重元素，范数守恒导致伪波函数仍偏「硬」，需要的截断能偏高。

超软赝势（USPP）：1990 年 Vanderbilt 提出，主动放弃范数守恒——允许 r_c 内电荷与真实不同，让伪波函数更平滑、截断能更低。损失的电荷用「补偿电荷密度」补回：

ρ_pseudo = Σ|ψ_pseudo|² + Σ Q_ij <ψ_pseudo|β_j><β_i|ψ_pseudo>
​

超软赝势把截断能要求大幅降低，尤其对过渡金属、稀土这类「硬」元素。代价是构造复杂、可能引入鬼态，使用要更谨慎。

2.3 PAW：全电子精度的「投影」

PAW（投影缀加波）由 Blöchl 1994 年提出，它不只是赝势，而是一个线性变换框架：

Ψ = T · Ψ_tilde
​

真实全电子波函数 Ψ 与平滑伪波函数 Ψ_tilde 之间通过变换算符 T 联系。在原子球外，两者相等（用平面波描述平滑部分）；原子球内，用预先算好的原子参考波函数重构出全电子行为：

Ψ(r) = Ψ_tilde(r) + Σ_i (φ_i - φ_tilde_i) <p_tilde_i | Ψ_tilde>
​

物理图像：把空间分成「原子球内」和「原子球外」。球外波函数平滑（平面波友好）；球内用原子轨道 + 投影函数精确重构。这样既保留了平面波的效率，又恢复了原子核附近的全部信息。

PAW 的优势：

• 接近全电子精度：芯区信息完整，可算 NMR、XAS 等芯敏感性质
• 可转移性好：基于原子参考态构造
• 普适：对过渡金属、稀土等「硬」元素同样高效
• 统一框架：VASP 等软件里 PAW 已成默认，兼有赝势的效率与全电子的精度

从历史脉络看，PAW 可以理解为「把范数守恒赝势、超软赝势与全电子线性缀加平面波（LAPW）的思想统一起来」：它既像赝势一样只在价层显式求解、又像全电子方法一样保留芯区信息。VASP 从早期赝势切换到 PAW 后成为行业标准，很大程度上正是因为这套「两全其美」的框架。

2.4 鬼态与验证

鬼态是超软赝势/PAW 构造不当可能出现的假本征态——在价电子能量范围里冒出一个不该有的态，计算结果被污染。现代势库已大幅降低鬼态概率，但遇到异常结果仍要警惕。

验证一个势是否可靠的标准动作：

1. 原子验证：算该元素自由原子的能级、电离能，与全电子结果对比
2. 固体验证：算晶格常数、体模量、结合能，与实验或高精度结果对比
3. 可转移性验证：换不同配位环境（分子 vs 晶体）对比价电子性质
4. 收敛测试：截断能与 k 点的收敛曲线，确认进入平台区

「验证」这件事很容易被新手跳过，但它才是势文件质量的最终裁判。官方势库声明「已验证」只代表通过了开发者设定的测试集，你面对的体系未必在其中。多花一小时做原子/固体双验证，能避免整轮生产计算的返工——这是计算里最划算的时间投资之一。

工程实践要点

1. 默认用软件官方势库：VASP 的 PAW 势库、QE 的赝势库都是经过严格测试的，优先使用。别自己从头造势。
2. 注意势文件与泛函匹配：一个势文件通常绑定某个交换关联泛函（PBE 势、LDA 势、杂化势）。混用泛函与势会产生不一致。检查软件警告。
3. 重元素优先 PAW/超软：过渡金属、稀土用 PAW 或超软赝势，范数守恒势会要求过高截断能。
4. 第一次用某个势必须验证：即使官方库，也建议先算晶格常数/体模量核对，再放行到生产计算。
5. 芯电子敏感性质要换方法：NMR、XPS、XAS 这类依赖核附近电子密度的性质，普通赝势不够，用 PAW（VASP）或全电子（FHI-aims、WIEN2k）方法。

赝势与 PAW 的分层透视

FAQ：几个高频疑问

问：赝势和 PAW 到底哪个「更对」？

答：PAW 更「接近全电子」，因为它能重构核附近信息；但它也不是全电子，芯电子仍是近似处理（用参考原子态描述）。对绝大多数材料性质，两者结果几乎一致。差别主要体现在芯敏感性质（NMR、芯能级谱）和磁性细节上——这类问题优先 PAW 或全电子方法。

问：为什么同一个元素有不同版本的势文件？

答：势文件按精度/效率分档，常见维度有：价电子数（是全价还是把次外层也算进价）、泛函版本（PBE/LDA/HSE 各有对应势）、半芯状态处理（重元素常把半芯轨道显式化）。选「够用且匹配」的那档，VASP 势库里通常有 _sv、_pv、_s 等后缀说明价电子范围。

问：势文件选错了会怎样？

答：轻则结果偏差（晶格常数、带隙漂移），重则直接不收敛或给出明显错误。所以验证很重要。多数软件的势库里不同版本之间会有明显的标记约定，先读文档再选势，别凭文件名猜。

一个加深理解的问题

为什么「散射性质匹配」是赝势的根本要求？因为价电子在晶体里的行为，本质上是在离子实阵列中的散射与传播。赝势只要保证「价电子被这个离子实散射的方式」和真实势一致，就能让能带结构、波函数、化学键都正确。反过来，如果一个赝势只保证「总能量接近」而不保证散射匹配，它在不同环境下的表现就会漂移——这就是可转移性差的根源。理解「散射匹配」这个核心，你就理解了一大半赝势理论。

⚠️ 常见坑：把用了不同势库（甚至不同泛函）的两个计算结果直接对比绝对能量。绝对能量依赖势的选择，只有同势同泛函的结果才可比。对比前先确认势与泛函设置一致。

💡 关键直觉：赝势与 PAW 成功的前提是一条物理事实——化学性质主要由价电子决定。所有「只算价电子」的方法都是在为这条事实下注。这个判断也决定了方法的边界：当你关心的问题恰好落在「芯电子」上，就得换方法。

概念性伪代码：赝势验证

def validate_pseudopotential(pseudo_file, element):
    atom = run_all_electron(element)          # 全电子参考
    pseudo = run_pseudo(element, pseudo_file) # 赝势计算
    check_levels_close(atom.levels, pseudo.levels)   # 能级一致
    check_lattice_close(atom.lattice, pseudo.lattice) # 晶格一致
    return "pass"
​

收束：一张心智地图

• 动机：全电子平面波在核附近成本爆炸，芯电子又不参与成键
• 核心思想：冻结芯电子，用平滑有效势替代核+芯的作用
• 赝势条件：散射匹配、正交、可转移、平滑
• rc 平衡：rc 越小越准越贵，rc 越大越省越糙
• 范数守恒 NCP：电荷守恒，精度好，重元素截断能高
• 超软 USPP：放弃范数守恒，截断能低，重元素高效
• PAW：线性变换重构全电子波函数，精度接近全电子
• 势库与验证：用官方库，原子+固体双重验证
• 鬼态：构造不当的假本征态，异常结果要排查
• 边界：芯敏感性质（NMR/XAS）需要 PAW 或全电子方法