16-pseudopotential 问答笔记¶
Q1：赝势这一步是不是只会影响 \(E_{xc}\) 的计算？¶
原文末尾的调试 checklist
不是。赝势影响的远不止 \(E_{xc}\)——它改变了整个待求解的问题（Hamiltonian、电子数、能量每一项都变了）。让我梳理清楚。
澄清可能的误解来源¶
你可能是这样联想的： - 化简 5 讨论 \(E_{xc}\) 泛函 - 化简 6 讨论赝势 - 上一题说"赝势和泛函强绑定" - → 产生"赝势只影响 XC"的印象
这是误解。赝势"携带了" XC 信息是一致性要求，但它的主要作用是重新定义整个电子问题。
赝势影响能量泛函的每一项¶
回忆 KS 总能量： $\(E = T_s + V_{ext}[\rho] + J[\rho] + E_{xc}[\rho]\)$
赝势对每一项都有影响：
| 能量项 | 有赝势 vs 无赝势的变化 | 影响程度 |
|---|---|---|
| \(V_{ext}\) 外势 | 真实库仑 \(-Z/r\) → 光滑赝势（含非局域算符） | 直接重定义——这是赝势的核心 |
| 电子数 \(N\) | 14 电子（Si 全电子）→ 4 电子（Si 价电子） | 大幅减少 |
| \(T_s\) 动能 | 含所有电子 → 只含价电子 | 少了一大块（内层贡献被冻结） |
| \(J[\rho]\) Hartree | 全电子密度 → 只价电子密度 | 内层库仑贡献打包进赝势 |
| \(E_{xc}[\rho]\) | 全电子密度 → 只价电子密度，配合 NLCC 处理内层-价电子交换 | 重定义计算方式 |
赝势的真正本质——换问题¶
没有赝势的问题： $\(\hat{H}_{\text{all-electron}} = -\tfrac{1}{2}\sum_{i=1}^{N_{\text{total}}} \nabla_i^2 + \sum_{i} \left(-\sum_I \frac{Z_I}{|\vec{r}_i - \vec{R}_I|}\right) + \sum_{i<j} \frac{1}{|\vec{r}_i - \vec{r}_j|}\)$
求解所有 \(N_{\text{total}}\) 个电子。
有赝势的问题： $\(\hat{H}_{\text{pseudo}} = -\tfrac{1}{2}\sum_{i=1}^{N_{\text{val}}} \nabla_i^2 + \sum_{i} V_{\text{pseudo}}(\vec{r}_i; \{\vec{R}_I\}) + \sum_{i<j} \frac{1}{|\vec{r}_i - \vec{r}_j|}\)$
只求解 \(N_{\text{val}}\) 个价电子，核+内层的效应塞进 \(V_{\text{pseudo}}\)。
注意关键变化： - 电子个数从 \(N_{\text{total}}\) 变成 \(N_{\text{val}}\)（Si 从 14 变 4） - \(V_{ext}\) 从 Coulomb 奇点变成光滑的 nonlocal 算符 - \(E_{xc}\) 的计算对象从"全电子密度"变成"价电子密度"——但计算\(E_{xc}\)本身的公式不变
赝势对其他化简步骤的间接影响¶
赝势不仅影响 \(E[\rho]\) 各项，还让其他化简步骤变可行：
| 化简步骤 | 没赝势 | 有赝势 | 关键影响 |
|---|---|---|---|
| 8 平面波基组 | ecutwfc > \(10^4\) Ry（不可行） | ecutwfc ~ 30–80 Ry | 决定性：这是赝势最大价值 |
| 4 KS 对角化 | 矩阵规模 \(\sim 10^8\) | 矩阵规模 \(\sim 10^3\)–\(10^4\) | 10^4 倍加速 |
| 11 SCF 收敛 | 内层波函数变化影响大 | 内层冻结，收敛更稳 | 收敛快 |
| 9 k 点采样 | BZ 结构复杂 | 不变 | 无关 |
化简 5 vs 化简 6：谁管 \(E_{xc}\)¶
| 化简 5（XC 泛函） | 化简 6（赝势） | |
|---|---|---|
| 主要目的 | 近似 \(E_{xc}[\rho]\) 的具体形式 | 替换全电子问题为价电子问题 |
| 决定 | LDA / PBE / SCAN / HSE06 | norm-conserving / USPP / PAW + 价电子数 |
| 对 \(E_{xc}\) 的影响 | 直接决定公式 | 只改变"输入 \(\rho\)"（价电子 vs 全电子） |
| 对 \(V_{ext}\) 的影响 | 无 | 完全重定义 |
| 对 \(N\) 的影响 | 无 | 大幅减少 |
结论：化简 5 专门管 \(E_{xc}\)；化简 6 影响整个 Hamiltonian 但主要改变 \(V_{ext}\) 和 \(N\)。
"赝势和 XC 泛函绑定"究竟意味着什么¶
（承接上一题的 Q2）
这个绑定关系是一致性要求——不是说"赝势决定 XC"：
- 生成赝势时用了某个 XC 泛函（如 PBE）去算真实原子的全电子结构
- 赝势的参数（截断半径、投影算符等）反映了这个泛函下的电子行为
- 跑 SCF 时如果 XC 换成 LDA，价电子感受的势就和它"原本应在的"世界不自洽
- 这会带来精度损失（~几十 meV/atom）
但赝势的主要功能不是"选 XC"——它是"把全电子问题变成价电子问题"。XC 信息只是附加标签，用来维持一致性。
类比（ML 视角）¶
把 DFT 整个求解流程比作训练 NN：
| 化简 | ML 对应 |
|---|---|
| 5 XC 泛函 | 选择损失函数形式（CE vs Focal vs MSE） |
| 6 赝势 | 选择预训练 backbone + 冻结策略（core 冻结 + valence 可训练） |
两件不同的事。赝势对损失函数形式只有"一致性"要求（backbone 预训练时用的归一化/损失函数最好和下游一致），但赝势本身不是损失函数。
一个具体例子：Si 的三种算法¶
| 方法 | 电子数 | ecutwfc | 单点能计算时间 | 能量精度 |
|---|---|---|---|---|
| 全电子 DFT（如 WIEN2k） | 14 | 等效 \(\sim 10^4\) Ry | 分钟–小时 | 参考 |
| PAW 赝势 + QE | 14 (frozen core + explicit valence) | 50 Ry | 秒–分钟 | ~1 meV/atom |
| norm-conserving 赝势 + QE | 4（仅价电子） | 40 Ry | 秒 | ~5–10 meV/atom |
看数字就知道：赝势影响的是能否算出来 + 多快算完——不只是 \(E_{xc}\) 这一项。
TritonDFT benchmark 中的影响¶
Benchmark 的每种材料评分都依赖：
1. 赝势选对（PBE 对应 PBE_standard 目录）
2. 截断能（ecutwfc）够大（关键！Ge 失败就在这里）
3. semi-core 处理合适（过渡金属的 3d 等）
赝势选错（如把 LDA 赝势用于 PBE 计算）可能导致： - 总能量偏差 ~几十 meV/atom - 晶格常数偏差 ~0.1–0.5% - 磁矩可能完全错
这些都不是 \(E_{xc}\) 单独的问题——是整个问题被"错误定义"的问题。
总结¶
- 不只影响 \(E_{xc}\)：赝势影响 \(V_{ext}\)、电子数 \(N\)、所有能量项
- 主要作用：把全电子问题换成价电子问题
- 附带要求：必须和 XC 泛函一致（但这不等于"只管 XC"）
- 对其他化简的影响：让 ecutwfc 可接受、KS 矩阵变小、SCF 收敛更易
- ML 类比：赝势 = 选 backbone + 冻结策略；XC 泛函 = 选损失函数。两件不同的事。
本轮 1 个问题速览¶
| # | 主题 | 核心 |
|---|---|---|
| Q1 | 赝势只影响 \(E_{xc}\) 吗？ | 不是——赝势影响整个 Hamiltonian（重定义 \(V_{ext}\)、减少电子数 \(N\)、让 ecutwfc 从 \(10^4\) Ry 降到 \(\sim 50\) Ry）；化简 5 专门管 \(E_{xc}\)，化简 6 管"换问题"；赝势与 XC 泛函的绑定是一致性要求而非"只管 XC" |